A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Reef exploration in the East Sengkang Basin,Sulawesi, Indonesia

Reconnaissance seismic shot in 1971/72 showed a number of well defined seismic anomalies within the East Sengkang Basin which were interpreted as buried reefs. Subsequent fieldwork revealed that Upper Miocene reefs outcropped along the southern margin of the basin. A drilling programme in 1975 and 1976 proved the presence of shallow, gas-bearing, Upper Miocene reefs in the northern part of the basin. Seismic acquisition and drilling during 1981 confirmed the economic significance of these discoveries, with four separate accumulations containing about 750 × 10⁹ cubic feet of dry gas in place at an average depth of 700 m. Kampung Baru is the largest field and contains over half the total, both reservoir quality and gas deliverability are excellent. Deposition in the East Sengkang Basin probably started during the Early Miocene. A sequence of Lower Miocene mudstones and limestones unconformably overlies acoustic basement which consists of Eocene volcanics. During the tectonically active Middle Miocene, deposition was interrupted by two periods of deformation and erosion. Carbonate deposition became established in the Late Miocene with widespread development of platform limestones throughout the East Sengkang Basin. Thick pinnacle reef complexes developed in the areas where reef growth could keep pace with the relative rise in sea level. Most reef growth ceased at the end of the Miocene and subsequent renewed clastic sedimentation covered the irregular limestone surface. Late Pliocene regression culminated in the Holocene with erosion. The Walanae fault zone, part of a major regional sinistral strike-slip system, separates the East and West Sengkang Basins. Both normal and reverse faulting are inferred from seismic data and post Late Pliocene reverse faulting is seen in outcrop.     

Multiple longshore sand bars in the upper Chesapeake Bay

In the upper Chesapeake Bay (Maryland, U.S.A.) field surveys were conducted at 18 multiple longshore sand bar sites. The multiple bar systems were found in water depths less than approximately 2 m (mean sea level), and exhibited mild bottom slopes of 0·0052 or less. The number of bars composing each system ranged from four to 17 and the spacing between the crests typically increased in the offshore direction, ranging from 12 to 70 m. Bar height also typically increased with distance offshore and ranged from 0·03 to 0·61 m. A grain size analysis of crest and trough sediment did not reveal any significant differences and the sediment was categorized as ‘fine sand’. A review of the literature data indicated that the Chesapeake Bay multiple bars possessed similar characteristics to those found in Gelding Bay (Baltic Sea); similarities in fetch, wave height and tidal range between the two bays may account for this finding. The surf-scaling parameter indicated that the multiple bar systems were extremely dissipative with regard to wave energy, and wave height appeared to be an important factor in controlling bar spacing and bar height. A multiple wave break point hypothesis was discussed as a possible mechanism for the formation of Chesapeake Bay multiple longshore bars, and limited observational evidence appeared to support such a mechanism.     

Microstructures of co-axially folded vein segments and crenulation cleavage： Evidence for dissolution phenomena in the Chamba Thrust Sheet, western Himalayas

The metasediments in the Chamba region experienced three phases of deformation: DF1, DF2 and DF3.Folded quartz veins are co-folded with the F2 crenulation folds. Their geometric and tectonic significance is studied by microstructures and shortening adjacent to the discrete crenulation cleavage, S2. Microstructures of folded vein segments, their geometric changes and truncation to cleavage (S2) are mainly due to pressure-solution phenomena and the estimated volume loss from reconstructed vein segments range from 16 to 25.5%,which is closely related to volume decrease (26%) calculated from the polydeformed slates of North Wales areas.     

Organotemplate silica deposition in Neoproterozoic deep-marine environments: evidence from the Penganga Group, Adilabad, India

Silica in bedded cherts interstratified with manganese carbonates in a deep‐water carbonate ramp succession of the Neoproterozoic Penganga Group, India, displays film‐like, spherical to rod‐shaped, and tubular branching microstructures. The microstructures resemble mineralized extracellular polysaccharides, biofilms and bacterial morphologies. The microstructures suggest silicification by nucleation of silica on organic‐templates or indirectly as sorbed species accumulating on organic templates. Given that similar microstructures have also been documented in Archean cherts it is proposed that organotemplates might have been an important sink for the deposition of silica in Precambrian deep‐water marine environments.     

A critical state model for sands dependent on stress and density

An elastoplastic model for sands is presented in this paper, which can describe stress–strain behaviour dependent on mean effective stress level and void ratio. The main features of the proposed model are: (a) a new state parameter, which is dependent on the initial void ratio and initial mean stress, is proposed and applied to the yield function in order to predict the plastic deformation for very loose sands; and (b) another new state parameter, which is used to determine the peak strength and describe the critical state behaviour of sands during shearing, is proposed in order to predict simply negative/positive dilatancy and the hardening/softening behaviour of medium or dense sands. In addition, the proposed model can also predict the stress–strain behaviour of sands under three-dimensional stress conditions by using a transformed stress tensor instead of ordinary stress tensor. Copyright © 2004 John Wiley & Sons, Ltd.