含裂隙组岩体的稳定问题

借鉴单晶体塑性理论,构建了岩体的多屈服面模型,用以模拟岩体中初始存在的多组裂隙。通过重构屈服函数,及简单的两种情况表明,可以采用常用的屈服准则作为屈服函数。随后,使用该多屈服面模型,分析了包含初始裂隙组岩体的稳定问题。在算例中研究了硬化模量随局部化条带方向的变化规律,并分析了不同破坏模式对临界硬化模量和临界局部化条带方向的影响。分析表明,该模型用于分析裂隙岩体的稳定问题是可行的。     

Gradient plasticity modelling of strain localization in granular materials

The flow stress in the yield surface of plastic constitutive equation is modified with a higher order gradient term of the effective plastic strain to model the effect of inhomogeneous deformation in granular materials. The gradient constitutive model has been incorporated into the finite element code ABAQUS and used to simulate biaxial shear tests on dry sand. It is shown that the shape of the post-peak segment of the load displacement curve predicted by the numerical analysis is dependent on the mesh size when gradient term is not used. Use of an appropriate gradient coefficient is shown to correct this and predict a unique shape of the load displacement curve regardless of the mesh size. The gradient coefficient required turns out to be approximately inversely proportional to the mesh elemental area. Use of the strain gradient term is found to diffuse the concentration of plastic strains within shear band resulting in its consistent width. The coefficient of the higher gradient term appears as a function of the grain size, the mean confining stress, and the plastic softening modulus. Copyright © 2004 John Wiley & Sons, Ltd.     

单轴压缩下岩石应变局部化的应变梯度塑性解

通过在屈服函数中引入材料内部长度建立应变梯度塑性增量本构关系,对岩石材料在单轴压缩条件下应变局部化的二维情况进行了研究。首先,建立了在Mises屈服准则下平面应变情况的岩石材料应变局部化带的带宽与其倾角的关系式,讨论了其带宽的分布范围和最大值与最小值形成的条件,给出岩石材料破坏Ⅰ, Ⅱ类变形的条件,并推出破坏的临界角度;其次,对岩石材料Ⅰ类变形进行了详细地分析,讨论了在不同泊松比 和 条件下材料应变局部化带的倾角的变化范围;最后,利用所得的结论同材料破坏的库仑准则进行对比研究,对内摩擦系数给予了合理的解释。     

基于梯度塑性模型的多孔介质应变局部化分析的基本理论

讨论了引入梯度塑性模型情况下的含液多孔介质双重内尺度律的特征以及相互作用问题。在引入梯度塑性本构模型的情况下,一方面材料应变局部化分析数值结果的正则性得到保证;而另一方面,含液多孔介质会出现双重内尺度律参数的共同作用问题。本文给出了此时内尺度律预测的一个基本方法,并对稳定性问题进行了分析,讨论了不同情况下实波速存在的条件,给出了对于给定的渗透系数情况下实波速存在的波数区间,并对相关现象进行了解释。     

塑性梯度体元及应变局部化启动机制分析

利用代表性体元可较好地反映出并非理想连续均匀的岩土介质破坏前强度和变形行为。但当材料出现高的应变梯度(应变局部化)时,在体元上的应力、应变呈现高次非线性变化,代表性体元不能如实地反映出介质的力学性能。提出一种可包含局部化变形的研究基元--塑性梯度体元(Plasticity Gradient Volume Element) ,塑性应变在其内可用微分形式表达。该体元不但可反映介质的平均力学响应,而且可反映应变呈高次非线性变化时的影响。随着硬(软)化模量从正值变为负值,变形模式将由均匀变形模式变为局部化变形模式,这就是应变局部化的启动机制。     

Reaction localization and softening of texturally hardened mylonites in a reactivated fault zone, central Argentina

The Tres Arboles ductile fault zone in the Eastern Sierras Pampeanas, central Argentina, experienced multiple ductile deformation and faulting events that involved a variety of textural and reaction hardening and softening processes. Much of the fault zone is characterized by a (D2) ultramylonite, composed of fine‐grained biotite + plagioclase, that lacks a well‐defined preferred orientation. The D2 fabric consists of a strong network of intergrown and interlocking grains that show little textural evidence for dislocation or dissolution creep. These ultramylonites contain gneissic rock fragments and porphyroclasts of plagioclase, sillimanite and garnet inherited from the gneissic and migmatitic protolith (D1) of the hangingwall. The assemblage of garnet + sillimanite + biotite suggests that D1‐related fabrics developed under upper amphibolite facies conditions, and the persistence of biotite + garnet + sillimanite + plagioclase suggests that the ultramylonite of D2 developed under middle amphibolite facies conditions. Greenschist facies, mylonitic shear bands (D3) locally overprint D2 ultramylonites. Fine‐grained folia of muscovite + chlorite ± biotite truncate earlier biotite + plagioclase textures, and coarser‐grained muscovite partially replaces relic sillimanite grains. Anorthite content of shear band (D3) plagioclase is c. An30, distinct from D1 and D2 plagioclase (c. An35). The anorthite content of D3 plagioclase is consistent with a pervasive grain boundary fluid that facilitated partial replacement of plagioclase by muscovite. Biotite is partially replaced by muscovite and/or chlorite, particularly in areas of inferred high strain. Quartz precipitated in porphyroclast pressure shadows and ribbons that help define the mylonitic fabric. All D3 reactions require the introduction of H⁺ and/or H₂O, indicating an open system, and typically result in a volume decrease. Syntectonic D3 muscovite + quartz + chlorite preferentially grew in an orientation favourable for strain localization, which produced a strong textural softening. Strain localization occurred only where reactions progressed with the infiltration of aqueous fluids, on a scale of hundreds of micrometre. Local fracturing and microseismicity may have induced reactivation of the fault zone and the initial introduction of fluids. However, the predominant greenschist facies deformation (D3) along discrete shear bands was primarily a consequence of the localization of replacement reactions in a partially open system.     

The geotechnical properties and behaviour of a pro-glacial lake clay, and its cementitious stabilization

The Tees Laminated Clay forms a continuous deposit around the Tees estuary in north east England and was deposited in a proglacial lake which occupied the Teesside area during late Devensian times, disappearing about 13 000 years ago. It overlies till, with which it has a sharp contact and reaches a maximum thickness in parts of central Middlesbrough of just over 9 m. The laminations consist of alternating thin layers of clay and silt, and represent seasonal increments of sediment. Kaolinite and illite, with lesser amounts of chlorite, dominate the mineral composition. Quartz tends to account for less than 35%. The notable carbonate content may be diagenetic in origin. The birefringence ratio of the clay minerals indicates a medium degree of particle orientation which suggests a partially flocculated fabric.This laminated clay has a high plasticity. However, the liquid and plastic limits of the silty layers differ significantly from those of the clay layers or the bulk samples. The clay fraction has a greater influence on these limits than the silty fraction. In addition, the linear shrinkage, liquidity indices and consistency indices of the silty layers differ appreciably from those of the clay layers and the bulk samples. The laminated clay is inactive. The orientation of the laminations in relation to the direction of maximum principal stress influences the shear strength, the lowest strengths being developed when these orientations are between 30° and 60°. The stress paths, anisotropic ratios and values of Af indicate that the Tees Laminated Clay is normally consolidated or lightly overconsolidated. The McLamore–Gray concept possibly can be applied to the anisotropic strength behaviour.Tees Laminated Clay can give rise to problems during construction operations since it can lose strength significantly on remoulding. Accordingly, these soils were treated with cementitious materials to determine the extent to which their engineering behaviour was enhanced by the addition of such materials. © Rapid Science Ltd. 1998     

边坡稳定分析中应变局部化的简化计算

采用分岔理论作为局部化判断条件,提出通过弹粘塑性软化模型来捕捉水工结构中的应变局部化剪切带的有限元简化算法。方法的特点是在每个时步,可以得到作为分岔理论所必须的弹塑性矩阵,仅用其作局部化判断,而不影响程序本身的粘塑性迭代。与其他方法相比,该方法理论简洁,易于程序的实现。通过对边坡的稳定分析,结果表明：该方法是合理可行的。     

梯度塑性理论及其在地质灾害预测研究中的应用

岩石结构稳定性问题是岩土工程实践中迫切需要解决的重要课题之一。文章以地震为例，讨论了地质灾害的局部化特征，介绍了经典弹、塑性理论的某些缺陷及梯度塑性理论的几点优越性：(1)控制方程总是适定的；(2)病态网格依赖性消失；(3)局部化带宽度由材料的内部长度完全确定；(4)可对岩石结构的尺寸效应及失稳回跳进行合理解释和预测；(5)对预测宏观及微观问题均比较有效。介绍了基于梯度塑性理论的岩石变形、破坏及稳定性研究进展。梯度塑性理论可用于研究单轴压缩条件下，岩石试件发生剪切破坏时全程应力一应变曲线、尺寸效应、剪切带倾角尺寸效应及失稳回跳等问题。它们对土木工程及岩土工程均十分重要。若将单轴压缩岩样比拟为矿柱，则该失稳判据即为矿柱岩爆准则。梯度塑性理论可用于研究韧性断层带内部应变、应变率分布规律、断层带错动位移及带内孔隙度分布规律，为韧性断层带定量分析提供了新的手段。此外，该理论还可对直剪试验机——岩样系统不稳定性进行分析，系统失稳可比拟为断层岩爆。发展基于梯度塑性理论尺度律及失稳判据等解析解一方面可加深对岩石变形、破坏的理解；另一方面。还可用来检验数值结果的正确性。     

Analytical and numerical investigation of uniqueness and localization in saturated porous media

This paper presents first the applications of uniqueness and strain localization analysis of saturated porous media, where localization of deformation into well defined narrow zones in a saturated porous medium is studied in terms of discontinuous bifurcation theory. A generalized plasticity constitutive model and a Mohr–Coulomb model are used in both the theoretical and numerical analyses of shear band formations. The critical hardening moduli and shear band angle for localization are computed, and quantitative results are given for both constitutive models. Numerical results previously obtained and new ones are confirmed by this analytical and numerical investigation. Copyright © 2002 John Wiley & Sons, Ltd.     

Comparison of continuous and discontinuous constitutive models to simulate concrete behaviour under mixed‐mode failure conditions

The paper presents detailed FE simulation results of concrete elements under mixed‐mode failure conditions according to the so‐called shear‐tension test by Nooru‐Mohamed, characterized by curved cracks. A continuous and discontinuous numerical two‐dimensional approach was used. In order to describe the concrete's behaviour within continuum mechanics, two different constitutive models were used. First, an elasto‐plastic model with isotropic hardening and softening was assumed. In a compression regime, a Drucker–Prager criterion with a non‐associated flow rule was used. In turn, in a tensile regime, a Rankine criterion with an associated flow rule was adopted. Second, an isotropic damage constitutive model was applied with a single scalar damage parameter and different definitions of the equivalent strain. Both constitutive laws were enriched by a characteristic length of micro‐structure to capture properly strain localization. As an alternative approach, the extended finite element method was used. Our results were compared with the experimental ones and with results of other FE simulations reported in the literature. Copyright © 2015 John Wiley & Sons, Ltd.     

Interaction between different internal length scales in strain localization analysis of fully and partially saturated porous media—the 1‐D case

The paper analyses the interaction between two internal length scales during dynamic strain localization in multiphase porous materials. The first internal length is introduced in the mathematical model by the gradient‐dependent plasticity for the solid skeleton, while the second one is naturally contained in the multiphase model and is due to the seepage process of the water via Darcy's law, which induces a rate‐dependent behaviour of the solid skeleton. Numerical results of a one‐dimensional example of water saturated porous medium demonstrate the competing effect between these two length scales. The porous medium is here treated as a multiphase continuum, with the pores filled by water and air, the last one at constant atmospheric pressure. Copyright © 2005 John Wiley & Sons, Ltd.     

基于梯度塑性本构理论的岩样侧向变形分析(II): 尺寸效应及弹性回跳

提出了利用不同尺寸试件的轴向应力-应变曲线得到轴向应力-侧向应变曲线的一种方法，并研究了结构尺寸对轴向应力-侧向应变曲线的影响。根据尺寸不同试件的轴向应力-轴向应变曲线的实验结果，并基于梯度塑性理论的解析解，在应变软化阶段，确定出了各种试件的剪切带条数。由此求出了不同宽度试样的轴向应力-侧向应变曲线。剪切带条数与试件宽度之比是决定轴向应力-侧向应变曲线特征的关键指标。若该比率为常量，则上述曲线不具有尺寸效应。当宽试件及窄试件在局部化启动后，在试件内部都出现一条剪切带时，随着试件宽度的增加，该曲线越陡，甚至出现回跳。上述方法也可用于分析轴向应力-环向应变曲线的尺寸效应问题。尺寸效应的原因是局部化，但局部化并非总引起尺寸效应。     

Compaction localization in the Earth and the laboratory: state of the research and research directions

Localized compaction in porous rocks is a recently recognized phenomenon that has been shown to reduce permeability dramatically. Consequently, the phenomenon is relevant to a variety of technologies involving fluid injection or withdrawal. This article summarizes current understanding of localized compaction and impediments to further progress. The article is based on discussions at a small workshop on localized compaction sponsored by the Office of Science, U. S. Department of Energy.     

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

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

考虑塑性应变率梯度的单轴压缩岩样轴向响应

基于剪切应变率梯度格式,采用解析方式研究了岩石材料在单轴压缩条件下的应变软化的结构响应。根据非局部连续介质模型,提出了一维二阶剪切应变率梯度格式。非局部剪切应变率与局部剪切应变率及其二阶梯度有关。将经典塑性理论中的局部剪切应变率替换为非局部剪切应变率,可以直接得到局部剪切应变率的封闭解析解,而不必通过将局部剪切应变对时间求导获得。通过对局部剪切应变率积分,得到了沿剪切带方向的相对剪切速度。试件峰值强度后的端部速度由弹性及塑性两部分构成。前一部分由虎克定律描述;后一部分与相对剪切速度有关。对弹性及塑性两部分速度求和,得到了单轴压缩岩样剪切破坏问题轴向响应的解析式。研究表明：试样高度越大、内部长度越小、剪切软化模量越大及泊松比越小,则岩样的轴向响应倾向于脆性。根据岩样与矿柱的相似性,岩样响应倾向于脆性,意味着矿柱将失去稳定性,发生矿柱岩爆。目前的基于剪切应变率梯度格式的主要优点是简洁。     

基于梯度塑性理论的岩样拉压剪破坏统一失稳判据

提出了单轴拉伸,直剪及单轴压缩剪切破坏失稳判据的统一形式。首先对de Borst R及Muhlhaus H B关于单轴拉伸试件梯度塑性理论解析解方面的开创性工作进行了回顾、评论,并指出了他们工作的缺憾。受de Borst R等工作的启发,对直接剪切试件进行了分析,得到了快速回跳条件。并发现:岩样的快速回跳条件就是系统(由剪切带及带外弹性岩石构成)的失稳判据。又对单轴压缩岩样剪切破坏进行了分析,得到了其失稳判据。通过对单轴拉伸失稳判据、直接剪切失稳判据及单轴压缩剪切破坏失稳判据的类比,得到了失稳判据的统一形式。统一失稳判据表明:岩石材料的局部化带宽度及压缩或剪切弹性模量越小,软化模量及等效试件高度越大,试件越容易发生失稳破坏。     

Coupled damage and plasticity modelling in transient dynamic analysis of concrete

In a concrete structure subjected to an explosion, for example a concrete slab, the material is subjected to various states of stress which lead to many modes of rupture. Closer to the explosive, a state of strong hydrostatic compression is observed. This state of stress produces an irreversible compaction of the material. Away from the zone of explosion, confinement decreases and the material undergoes compression with a state of stress, which is slightly triaxial. Finally, the compression wave can be reflected on a free surface and becomes a tensile wave, which by interaction with the compression wave, produces scabbing. We present, in this paper, a model aimed at describing these three failure modes. It is based on visco‐plasticity and rate dependent damage in which a homogenization method is used in order to include the variation of the material porosity due to compaction. The model predictions are compared with several experiments performed on the same concrete. Computations of split Hopkinson tests on confined concrete, a tensile test with scabbing, and an explosion on a concrete slab are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Biological soil crust development and its topsoil properties in the process of dune stabilization, Inner Mongolia, China

This study was undertaken at Horqin Sand Land, Inner Mongolia, Northern China. Field samples of biological soil crusts (BSCs) and underlying topsoil (0–5 cm under BSC) were taken in areas of different dune stabilization stages, and their physicochemical properties were analyzed, including particle size distribution, bulk density, organic matter, nitrogen, phosphorus, electrical conductivity (EC), pH, and CaCO₃ content. The results revealed that semi-mobile dunes, semi-fixed dunes and fixed dunes had developed a physical crust, algae crust and moss crust, respectively. The thickness, hardness, water content, fine fraction and nutrient contents of BSCs were gradually increasing along the dune stabilization gradient. Meanwhile, BSC establishment and development enhanced the bulk density, silt and clay content, and nutrients of the topsoil under it, in an increasing tend from semi-mobile dune to fixed dune. Organic matter concentrations and other nutrients in the 0–5 cm topsoil layer under BSCs were significantly higher compared to unconsolidated soil (control). Moreover, there were strong significant positive correlations between topsoil and BSCs’ organic matter, total nitrogen, available nitrogen, available phosphorus, CaCO₃, and <0.05 mm particle content, suggesting that BSCs have an influence on some of the properties of the underlying topsoil.     

Shear localization,velocity weakening behavior,and development of cataclastic foliation in experimental granite gouge

Microstructural aspects of room-temperature deformation in experimental Westerly granite gouge were studied by a set of velocity stepping rotary-shear experiments at 25 MPa normal stress. The experiments were terminated at: (a) 44 mm, (b) 79 mm, and (c) 387 mm of sliding, all involving variable-amplitude fluctuations in friction. Microstructural attributes of the gouge were studied using scanning (SEM) and scanning transmission electron microscopy (STEM), image processing, and energy dispersive X-ray (EDX) analyses. The gouge was velocity weakening at sliding distances >10 mm as a core of cataclasites along a through-going shear zone developed within a mantle of less deformed gouge in all experiments. Unlike in experiment (a), the cataclasites in experiments (b) and (c) progressively developed a foliation defined by stacks of shear bands. The individual bands showed an asymmetric particle-size grading normal to shearing direction. These microstructures were subsequently disrupted and reworked by high-angle Riedel shears. While the microstructural evolution affected the effective thickness and frictional strength of the gouge, it did not affect its overall velocity dependence behavior. We suggest that the foliation resulted from competing shear localization and frictional slip hardening and that the velocity dependence of natural fault gouge depends upon compositional as well as microstructural evolution of the gouge.