A model of strain localization in porous sandstone as a function of tectonic setting,burial and material properties; new insight from Provence (southern France)

The analysis of three cataclastic band sets from Provence (France) reveals that the band density, their conjugate angles, their ratio of shear displacement to compaction, and the amount of cataclasis within the bands differ and can be expressed as functions of tectonic setting and petrophysical properties. We identify (1) a dense and closely spaced network of shear-enhanced (reverse) compaction bands; (2) a regularly spaced less dense network of reverse compactional shear bands; and (3) a localized network of normal shear bands. The field data show that strain localization is favored in an extensional regime and is characterized by shear bands with a large shear to compaction ratio and a small conjugate band angle. In contrast, distributed strain is favored in a contractional regime and is characterized by compactional bands with a low ratio of shear to compaction and a large conjugate band angle. To explain the mechanical origin of this strain localization, we quantify the yield strength and the stress evolution in extensional and contractional regimes in a frictional porous granular material. We propose a model of strain localization in porous sands as a function of tectonic stresses, burial depth, material properties, strain hardening and fluid pressure. Our model suggests that stress reduction, inherent to extensional regime, favors strain localization as shear bands, whereas stress increase during contraction favors development of compactional bands.     

基于梯度塑性本构理论的岩样侧向变形分析(Ⅰ):基本理论及本构参数对侧向变形的影响

研究了岩样在单轴压缩条件下轴向应力．侧向或环向变形的全程曲线特征。基于考虑峰值剪切强度后微小结构之间相互影响和作用的梯度塑性理论，得到了由于剪切局部化而引起的侧向塑性变形。利用虎克定律描述了试件的弹性变形，得到了轴向应力．侧向变形全程曲线的解析解。在软化阶段，试件中部侧向变形及对靠近试件上端或下端部位的侧向变形并不相同。与轴向应力．应变曲线可能出现的回跳现象类似。试件中部轴向应力．侧向应变曲线也可能出现回跳现象。在应变软化阶段，与应力．侧向应变曲线相比，应力．环向应变曲线不容易发生回跳现象。若在试件内部出现多条剪切带，则应该以等效剪切带宽度替代本文中的剪切带宽度。随着剪切带倾角、内部长度参数的降低、剪切模量的增加及弹性模量的降低，轴向应力．侧向应变曲线越陡；甚至能出现弹性回跳。     

Capturing strain localization in reinforced soils

Lade’s single hardening soil model with Cosserat rotation embodied in the finite element method is employed to investigate the behavior of geosynthetic reinforced soils with special attention to the development of shear banding. The ability of the finite element model to detect shear banding in a reinforced soil is examined against three high quality small-scale laboratory plane strain tests on Toyoura sand with and without reinforcement. These three tests were chosen because of the clear failure surfaces that developed in the soil during loading. The FEM analyses were able to reasonably simulate the plane strain laboratory tests including both unreinforced and reinforced cases. The FEM analyses gave reasonably good agreement with the experimental results in terms of global stress–strain relationships and shear band occurrences. Furthermore, and based on FE analyses of a hypothetical geosynthetic reinforced soil (GRS) retaining wall, it is shown that the geosynthetic reinforcements are very effective in hindering the formation of shear bands in GRS retaining walls when small spacing between the reinforcement layers was used. When used properly, the geosynthetic reinforcements made the soil behave as a truly reinforced mass of considerable stiffness and strength.     

不同应力路径下超固结黏土试样变形局部化分析

基于改进伏斯列夫面超固结黏土三维本构模型,利用有限元软件ABAQUS材料子程序接口,采用回映应力更新算法,实现了该模型在有限元分析中的应用。通过该模型与比奥固结理论的耦合,对超固结比为8的超固结黏土在三轴压缩、三轴伸长及平面应变应力条件下的变形局部化问题,进行了水-土耦合弹塑性有限元分析。分析结果表明：剪切带带内、带外点经历不同应力路径;剪切带带外单元经历了体缩、剪胀及被吸水体缩过程,而剪切带带内单元一直保持剪胀趋势;剪切带的形成伴随着剪胀,剪切带内、外出现了负的孔压,且孔压的分布也具有局部化特性。关于剪切带带内、外的孔隙水压及体变变化趋势与剪切速率有关,而平面应变介于三轴压缩与三轴伸长之间,但平面应变较早出现剪切带。孔隙水的迁移速度影响剪切带带内单元的剪胀,进而影响剪切带的形成及发展;而围压和弱单元位置也对剪切带的形成也有影响。     

Influence of relative density,particle shape,and stress path on the plane strain compression behavior of granular materials

Two types of stress path-controlled plane strain compression tests were performed on both loose and dense specimens of angular and sub-angular sands and two rounded glass beads with different particle sizes. Digital image correlation method was used to analyze local deformation developments, especially shear band patterns. The material behavior in response to shearing has been found to be dependent on the relative density, particle shape, and stress path. The results of analysis on local deformation developments showed that the onset of shear bands occurred prior to their peak strengths in both dense and loose specimens. The growth rates of local maximum shear strain along a shear band were approximately consistent with an increasing global axial strain after the onset of shear band. The shear band width was influenced by both the mean particle size and the particle shape. The measured shear band inclination angles were in between those estimated by Coulomb’s and Roscoe’s formulas.     

Strain analysis using quartz deformation bands

Quartz deformation bands are kink bands in quartz crystals. A deformation band develops as a region of localized crystal-plastic deformation with boundaries perpendicular to the slip plane and slip direction, which usually is along an -axis in the basal plane. Under cross-polarized light, the difference in crystallographic orientation between a deformation band and its host is indicated by a difference in extinction positions. The displacement between the c axis in a deformation band and the c axis in the host represents the angular shear of the deformation band in the direction of the c axis in the host grain. Assuming the deformation is homogeneous at the grain scale, the angular shear of the grain (the gauge) is calculated by multiplying the angular shear of the deformation band by the ratio of the sheared part to the whole grain. Using the strain-gauge method for three-dimensional infinitesimal strain analysis, a minimum number of five grains measured on universal stage is needed to solve for the deviatoric strain components of the aggregate if the strain is homogeneous in the aggregate. Data from more than five grains are used to find the best-fit strain components by a least-squares method. The principal strains and their orientations are found from these strain components by calculating the eigenvalues and eigenvectors. A 3-D strain ellipsoid also is obtained from strain ellipses in three perpendicular planes determined from the two-dimensional flat-stage measurements by the Wellman method. Both the strain-gauge method and the Wellman method are tested by using synthetic data sets and applied to a naturally deformed sample. Both methods give similar results; the established Wellman method thus confirms the strain-gauge calculation.     

Centrifuge model tests on sand specimen under extensional load

The appearance of shear banding in granular materials has been investigated intensively during the last decades and is still of ongoing importance in terms of understanding the stress–strain behaviour of the material, the localization phenomena and the interaction between soil and structure. Only less attention has been paid to the occurrence of systems of shear bands although such systems can be identified in geotechnical structures as well as in geological formations. In this paper we present results of experiments on sand specimens under extensional load in natural gravity as well as in increased gravity in the centrifuge where the influence of the stress level on the geometry of a shear band pattern, specified by the spacing of the shear bands and the angle between failure surfaces and minor stress direction, has been investigated. X‐ray technique has been used to visualize the failure zones inside the specimen, an optical measurement system called Digital Image Correlation has been applied to identify and observe the appearing deformation mechanism on the sides of the specimens in natural gravity as well as during the flight in the centrifuge. It can be shown that the geometry of the shear band pattern is sparsely influenced by the change of the stress level. Copyright © 2004 John Wiley & Sons, Ltd.     

单轴压缩黏性土剪切带相互作用及损伤试验研究

剪切带损伤研究对于理解材料破坏机制和建立剪切带本构模型具有重要意义。为了研究单轴压缩黏性土样的剪切带损伤演化规律,根据纵向应变较高时清晰剪切带位置布置切向测线,对数字图像相关方法获得的结果进行双三次样条插值,从而获得光滑性较好的各种应变场。将土样整体的损伤变量-纵向应变曲线与各条剪切带的损伤演化曲线进行了对比。研究发现：（1）总体上,土样整体的损伤变量演化曲线呈线性,而各条剪切带的损伤变量演化曲线均上凹,这表明随着纵向应变的增加,各条剪切带的损伤发展越来越快;（2）各条剪切带的损伤变量演化曲线的轮廓线呈马尾形,这说明随着剪切带的逐渐发育,各条剪切带的相互影响和作用规律变强;（3）对于含水率较低的土样,平行或共线剪切带的损伤变量变化基本同步,特别是在剪切带充分发育之后,但两条剪切带的距离应较小;在一定时期,共轭或交叉剪切带中剪切带的损伤占优,这与剪切带的相互竞争有关,但若两条剪切带达到独立发展的程度,则二者的损伤变量变化可以同步;对于含水率较高的土样,各条平行或共轭剪切带的损伤变量变化基本同步,剪切带的相互作用不明显。     

Finite strain analysis of nonuniform deformation inside shear bands in sands

A methodology has been developed to extend the incremental (Eulerian) Digital Image Correlation (DIC) technique to enable a Lagrangian‐based large‐strain analysis framework to examine the nature of strain and kinematic nonuniformity within shear bands in sands. Plane strain compression tests are performed on dense sands in an apparatus that promotes unconstrained persistent shear band formation. DIC is used to capture incremental, grain‐scale displacements in and around shear bands. The performance of the developed accumulation algorithm is validated by comparing accumulated displacements with two sources of reference measurements. A comparison between large and infinitesimal rotation is performed, demonstrating the nature of straining within shear bands in sands and the necessity of using a finite strain formulation to characterize ensuing behavior. Volumetric strain variation along the shear band is analyzed throughout macroscopic postpeak deformation. During softening, volumetric activity within the shear band is purely dilative. During the global critical state, the shear band material is seen on the average to deform at zero volumetric strain; however, locally, the sand is seen to exhibit significant nonzero volumetric strain, putting into question the current definition of critical state. At the softening‐critical state transition, a spatially periodic pattern of alternating contraction and dilation along the shear band is evidenced, and a preliminary evaluation indicates that the periodicity appears to be a physical phenomenon dictated only in part by median grain size. Copyright © 2011 John Wiley & Sons, Ltd.     

不同应力路径下剪切带的数值模拟

采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。     

Effect of consolidation on strain localization of soft clays

A finite element formulation developed by Belytschko and his coworkers [1] is used to examine localized deformation as it exists for laboratory specimens of compressible Chicago clays. The element is based on an assumed strain formulation wherein localized and non-localized zones are embedded into an element when material based bifurcation is detected. A plane strain compression test of a natural clay specimen in which deformations localized into a single shear band during undrained shear is simulated using this element. Localization is initiated by imposing a non-uniform applied displacement consistent with that measured experimentally. Results indicate that to insure localization to a single band, considerations must be made for non-uniformities developed during the consolidation phase.     

超固结粘土的剪切带数值模拟

基于中井的子负荷面本构关系模型,采用高精度的隐式应力积分算法,研究了平面应变试验中超固结粘土试样的变形局部化问题。模拟了不同的加载速度和边界约束条件,得到了超固结比为8的藤森粘土土体内部明显可见的X型剪切带和单一型剪切带。结果表明：土体的剪胀特性、孔隙水在土体内的移动规律以及边界约束条件对剪切带的形成和发展起着控制作用。     

Stress-driven Melt Segregation and Strain Partitioning in Partially Molten Rocks: Effects of Stress and Strain

The evolution of melt segregation in deforming partially moltenolivine-rich rocks has been studied in a series of laboratoryexperiments. During deformation, melt segregates into networksof anastamosing channels (or ‘bands’) surroundinglenses of melt-depleted material. We quantify the nature ofthe melt distribution in the samples, including thickness, angle,spacing, volume fraction, and melt fraction of melt-rich bands,to understand the dynamics of melt-network organization. Twoseries of experiments were designed to isolate the effects of(1) increasing shear strain (at similar stress levels), and(2) varying stress levels (deformed to similar shear strains).Melt-rich bands develop by a shear strain of unity. In samplesdeformed at varying stress levels, higher stress produces smallercharacteristic band spacings. We relate these variations tothe compaction length, _c, which varies only as a result of thereduction of matrix viscosity with increasing stress. Simpleapproaches to scaling from experimental to mantle conditionssuggest that stress-driven melt segregation can occur in theasthenosphere; if so, it will significantly affect rheological,transport and seismic properties, with enticing consequencesfor our understanding of plate–mantle interactions. KEY WORDS: melt segregation; rock rheology; magma transport; self-organization; mid-ocean ridges     

考虑应变局部化的粗粒料剪切损伤力学机制

针对粗粒料的应变软化、剪胀等力学特性,通过考虑以剪切带为标志的应变局部化现象,建立了具有广泛适用性的剪切损伤力学模型。损伤模型采用了包体理论中的剪切带数学简化,基于应变等价原理、Weibull分布,推导了粗粒料的应力-应变关系方程。从剪胀作用的机制出发,提出可以描述剪胀弱化的轴向塑性应变和体积塑性应变的非线性函数关系。结合粗粒料三轴压缩试验中的伺服过程,提出了基于遗传算法的损伤模型参数确定方法。通过开展不同围压下的粗粒料三轴压缩试验,对剪切损伤力学模型进行验证,进一步分析了参数演化对粗粒料强度和变形特征的影响。研究结果表明,考虑应变局部化特征的剪切损伤力学模型可以高精度的模拟粗粒料的应变软化和剪胀等特征,有效揭示剪切带内部变形对试样整体宏观变形的影响机制,模型中剪切带参数和围压的关系与粗粒料细观机制一致,计算得到强度组成与颗粒破碎、重组特征较为吻合。     

Development of quartz c-axis crossed/single girdles under simple-pure shear deformation: Results of visco-plastic self-consistent modeling

Quartz c-axis fabrics are widely used to determine the shear plane in ductile shear zones, based upon an assumption that the shear plane is perpendicular to both the central segment of quartz c-axis crossed girdle and single girdle. In this paper the development of quartz c-axis fabric under simple-pure shear deformation is simulated using the visco-plastic self-consistent (VPSC) model so as to re-examine this assumption. In the case of no or weak dynamic recrystallization, the simulated crossed girdles have a central segment perpendicular or nearly perpendicular to the maximum principal finite strain direction (X) and the XY finite strain plane, and at a variable angle relative to the imposed kinematic framework that is dependent on the modeled flow vorticity and finite strain. These crossed girdles have a symmetrical skeleton with respect to the finite strain axes, regardless of the bulk strain and the kinematic vorticity, and rotate in a way similar to the shear sense with increasing bulk strain ratio. The larger the vorticity number the more asymmetrical their legs tend to be. In the case of strong dynamic recrystallization and large bulk strain, under simple shear the crossed girdle switches into single girdles, sub-perpendicular to the shear plane, by losing the weak legs. The numerical results in our models do not confirm the above-mentioned assumption.     

材料缺陷对岩样变形局部化影响的数值模拟

在平面应变压缩条件下,采用拉格朗日元法研究了材料缺陷对岩样应变局部化及宏观力学行为的影响。在数值计算中,采用了莫尔-库仑与拉破坏复合的破坏准则,峰后岩石的本构关系为线性应变软化。对于理想岩样（不含任何缺陷）,从剪切应变率的等值线图及变形网格图发现,变形场相对试样垂直对称轴对称。在试样的边界上设置材料缺陷后,变形场的对称性被打破。材料缺陷附近是局部化现象启动的主要位置。与不含材料缺陷试样相比,含缺陷试样的局部化提前启动。当缺陷位于试样左边界中部附近时,试样内部出现多条剪切带,发生韧性剪切破坏,峰后应力-轴向应变曲线和应力-侧向应变曲线均倾向于韧性,岩样的稳定性增强;当缺陷位于试样上、下端面附近时,仅出现一条贯穿试样左右边界的剪切带,发生脆性剪切破坏,峰后两种曲线倾向于脆性,易发生失稳破坏。从数值结果中还观测到了剪切带的跳跃或迁移现象及剪切带的相互竞争、此消彼涨的现象。当岩样中不包含任何缺陷时,试样应力-变形曲线的峰值强度最高。随着缺陷的上移,峰值强度下降,直到峰值强度基本保持不变。缺陷越接近于岩样的下端,对岩样应力-变形曲线的峰值强度的影响越小。     

上海粘性土剪切带倾角的试验研究

采用可量测局部侧向应变的平面应变仪，对上海地区第（2）层褐黄色粉质粘土和第（6）层暗绿色粉质粘土的原状土样，进行固结不排水平面应变压缩试验，观察土样的局部化变形过程，量测剪切带倾角，研究剪切带倾角随固结压力的变化情况，计算剪切带倾角的传统土力学理论解并与实测值进行比较，根据剪切带内局部剪胀角的近似估算，初步研究剪切带内的变形机理，认为剪切带的形成与孔隙水的运动密切相关。     

东天山大型韧性剪切带基本特征与金矿预测

东天山大型韧性剪切带即秋格明塔什-黄山韧性剪切带,规模宏大,东西长逾600km,宽5～20km。变形标志明显,分四期变形。应变测量、岩组分析、剪切位移量计算、同位素测年等定量分析表明：①以平面单剪应变为主,稍晚向伸长应变转化;②存在逆冲剪切和脆-韧性变形转换;③剪切位移量在75km以上;④变形时代在300～250Ma之间,高峰期在285～265Ma,属海西中晚期。它是康-黄断裂带的重要组成体,控制着康古尔塔格金矿带的空间展布。据现有大中型金矿研究建立了区域脆-韧性变形转换带结构模式和金矿预测模式,可有效指导找矿勘探。     

剪切带图案的复杂性及应力-应变曲线的离散性

在平面应变压缩条件下,通过在岩样左边界靠近下端面的位置及右边界的不同位置预制材料缺陷,采用拉格朗日元法研究了材料缺陷位置不同时的剪切带图案启动、演变及试样的宏观力学行为。在数值计算中,采用了莫尔-库仑与拉破坏复合的破坏准则,峰后岩石的本构关系为线性应变软化。当试样左、右边界上的缺陷距离较远时,剪切带仅出现1条,以通过两缺陷形式贯通试样两边界;当两缺陷距离较近时,启动于缺陷附近的剪切带按各自其固有方向发展,剪切带图案十分复杂。当右缺陷距离试样上端面较近时,由于上端面的约束,仅位于缺陷下方的剪切带能得到充分的发展。若剪切带数量较少且贯通试样两边界时,试样发生脆性剪切破坏,应力-变形曲线软化段比较陡峭;反之,发生韧性剪切破坏,软化段比较平缓。当剪切带的前端迎面遇到试样上、下端面附近的上、下三角形区域时,由于端面的约束,剪切带的方向将改变。当剪切带的前端与试样的纵向对称线距离较近时,剪切带发生折射现象;当距离较远时,剪切带将发生反射现象。     

Kinematic modelling of shear band localization using discrete finite elements

Modelling shear band is an important problem in analysing failure of earth structures in soil mechanics. Shear banding is the result of localization of deformation in soil masses. Most finite element schemes are unable to model discrete shear band formation and propagation due to the difficulties in modelling strain and displacement discontinuities. In this paper, a framework to generate shear band elements automatically and continuously is developed. The propagating shear band is modelled using discrete shear band elements by splitting the original finite element mesh. The location or orientation of the shear band is not predetermined in the original finite element mesh. Based on the elasto‐perfect plasticity with an associated flow rule, empirical bifurcation and location criteria are proposed which make band propagation as realistic as possible. Using the Mohr–Coulomb material model, various results from numerical simulations of biaxial tests and passive earth pressure problems have shown that the proposed framework is able to display actual patterns of shear banding in geomaterials. In the numerical examples, the occurrence of multiple shear bands in biaxial test and in the passive earth pressure problem is confirmed by field and laboratory observations. The effects of mesh density and mesh alignment on the shear band patterns and limit loads are also investigated. Copyright © 2003 John Wiley & Sons, Ltd.