An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: Case of hydrostatic pore‐water pressure

The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd.     

A geochronological and petrological study of anatectic paragneiss and associated granite dykes from the Day Nui Con Voi metamorphic core complex,North Vietnam: constraints on the timing of metamorphism within the Red River shear zone

The Red River shear zone (RRSZ) is a major left‐lateral strike‐slip shear zone, containing a ductilely deformed metamorphic core bounded by brittle strike‐slip and normal faults, which stretches for >1000 km from Tibet through Yunnan and North Vietnam to the South China Sea. The RRSZ exposes four high‐grade metamorphic core complexes along its length. Various lithologies from the southernmost core complex, the Day Nui Con Voi (DNCV), North Vietnam, provide new constraints on the tectonic and metamorphic evolution of this region prior to and following the initial India–Asia collision. Analysis of a weakly deformed anatectic paragneiss using P–T pseudosections constructed in the MnO–Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (MnNCKFMASHTO) system provides prograde, peak and retrograde metamorphic conditions, and in situ U–Th–Pb geochronology of metamorphic monazite yields texturally controlled age constraints. Tertiary metamorphism and deformation, overprinting earlier Triassic metamorphism associated with the Indosinian orogeny and possible Cretaceous metamorphism, are characterized by peak metamorphic conditions of ~805 °C and ~8.5 kbar between c. 38 and 34 Ma. Exhumation occurred along a steep retrograde P–T path with final melt crystallizing at the solidus at ≥~5.5 kbar at ~790 °C. Further exhumation at ~640–700 °C and ~4–5 kbar at c. 31 Ma occurred at subsolidus conditions. U–Pb geochronological analysis of monazite from a strongly deformed pre‐kinematic granite dyke from the flank of the DNCV provides further evidence for exhumation at this time. Magmatic grains suggest initial emplacement at 66.0 ± 1.0 Ma prior to the India–Asia collision, whereas grains with metamorphic characteristics indicate later growth at 30.6 ± 0.4 Ma. Monazite grains from a cross‐cutting post‐kinematic dyke within the core of the DNCV antiform provide a minimum age constraint of 25.2 ± 1.4 Ma for the termination of fabric development. A separate and significant episode of monazite growth at c. 83–69 Ma is suggested to be the result of fluid‐assisted recrystallization following the emplacement of magmatic units.     

On the stability of stratified plane couette flow

Abstract

Unbounded stratified plane Couette flow is shown to be stable against small amplitude disturbances. The Brunt-Väisälä frequency is assumed to be constant. Both viscosity and thermal diffusion are included, and shown to be stabilizing.     

Coupled finite‐element simulation of injection well testing in unconsolidated oil sands reservoir

This paper presents a finite‐element (FE) model for simulating injection well testing in unconsolidated oil sands reservoir. In injection well testing, the bottom‐hole pressure (BHP) is monitored during the injection and shut‐in period. The flow characteristics of a reservoir can be determined from transient BHP data using conventional reservoir or well‐testing analysis. However, conventional reservoir or well‐testing analysis does not consider geomechanics coupling effects. This simplified assumption has limitations when applied to unconsolidated (uncemented) oil sands reservoirs because oil sands deform and dilate subjected to pressure variation. In addition, hydraulic fracturing may occur in unconsolidated oil sands when high water injection rate is used. This research is motivated in numerical modeling of injection well testing in unconsolidated oil sands reservoir considering the geomechanics coupling effects including hydraulic fracturing. To simulate the strong anisotropy in mechanical and hydraulic behaviour of unconsolidated oil sands induced by fluid injection in injection well testing, a nonlinear stress‐dependent poro‐elasto‐plastic constitutive model together with a strain‐induced anisotropic permeability model are formulated and implemented into a 3D FE simulator. The 3D FE model is used to history match the BHP response measured from an injection well in an oil sands reservoir. Copyright © 2013 John Wiley & Sons, Ltd.     

Nonlinear analysis of laterally loaded rigid piles in cohesive soil

This article presents a method for the nonlinear analysis of laterally loaded rigid piles in cohesive soil. The method considers the force and the moment equilibrium to derive the system equations for a rigid pile under a lateral eccentric load. The system equations are then solved using an iteration scheme to obtain the response of the pile. The method considers the nonlinear variation of the ultimate lateral soil resistance with depth and uses a new closed‐form expression proposed in this article to determine the lateral bearing factor. The method also considers the horizontal shear resistance at the pile base, and a bilinear relationship between the shear resistance and the displacement is used. For simplicity, the modulus of horizontal subgrade reaction is assumed to be constant with depth, which is applicable to piles in overconsolidated clay. The nonlinearity of the modulus of horizontal subgrade reaction with pile displacement at ground surface is also considered. The validity of the developed method is demonstrated by comparing its results with those of 3D finite element analysis. The applications of the developed method to analyze five field test piles also show good agreement between the predictions and the experimental results. The developed method offers an alternative approach for simple and effective analysis of laterally loaded rigid piles in cohesive soil. Copyright © 2011 John Wiley & Sons, Ltd.     

Garnet crystal plasticity in the continental crust,new example from south Madagascar

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.     

A mixing-length model for predicting vertical velocity distribution in flows through emergent vegetation

Abstract

The vertical profiles of streamwise velocities are computed on flood plains vegetated with trees. The calculations were made based on a newly developed one-dimensional model, taking into account the relevant forces acting on the volumetric element surrounding the considered vegetation elements. A modified mixing length concept was used in the model. An important by-product of the model is the method for evaluating the friction velocities, and consequently bed shear stresses, in a vegetated channel. The model results were compared with the relevant experimental results obtained in a laboratory flume in which flood plains were covered by simulated vegetation.     

宽窄相间河道水流紊动特性试验研究

为了探索宽窄相间河道的水流紊动特性，以西南地区宝兴河上游宽窄相间河段为研究对象，基于室内概化模型试验，采用多普勒声学流速仪（ADV）测量了室内模型典型断面上的三维瞬时流速，分析典型断面上的纵向时均流速、紊动强度、雷诺切应力和紊动能的分布规律。试验结果显示：宽窄相间水槽中，扩散段边壁的紊动强度大于中心区域的紊动强度，最大值位于0.2倍水深处；扩散段两侧坡脚处紊动能最大；侧壁区的平面和立面雷诺切应力最大值出现在扩散段内，中心区域最大雷诺切应力位于两槽间的中间断面处；扩散段内水流紊乱，两侧出现旋涡和涡脱，易造成侧壁侵蚀加强，引起河道拓宽。深入分析了宽窄相间河道水流的紊动特性，可为山区河流治理和自然灾害防治提供参考。     

大别变质地体的构造样式及变形序列

大别山区晚太古代大别群变质岩系,经历了复杂的变形过程。在大别变质地块内部识别出两套大型剪切带,它们在该区构造格架形成过程中起着重要的作用。本文还分析了变质的上壳岩、片麻岩和剪切带内的变形特点。实际资料说明,剪切带内的褶皱是由于剪切的不稳定性及平行面理局部缩短形成的,这些大型地壳剪切作用多伴随有挤压流动或伸展流动。根据褶皱叠加关系以及剪切带、新成体交切关系建立了大别群总体构造序列。该区提供了一个研究中下地壳作用的窗口。     

Shear resistance of jointed rock masses and stability calculation of rock slopes

Summary Based on an extension to a model originally proposed by Zhang (1990), a method for estimating the shear strength of a jointed rock mass containing two joint sets is developed. Compared to the original method in which only a single en-echelon joint set was considered, the new method has the following advantages: (1) it allows a much more realistic joint pattern to be incorporated; (2) the effects of joint connectivity and intersection angle between joint sets on shear strength are included. Use of the new method is demonstrated by application to the calculation of the stability of slopes in jointed rock. The importance of joint connectivity, joint persistence and joint set intersection angle on slope stability, as predicted by the model, is revealed.