月壤双轴试验的剪切带离散元数值分析

土体的破坏问题一般都是从其剪切带入手进行研究的。针对真实月壤所处的环境（无水、低重力场、低气压等）和内摩擦角较大的特点,采用最近提出的1种考虑粒间抗转动作用与粒间范德华力2个因素的月壤散粒体力学接触模型,用离散单元法模拟了柔性边界条件下月球环境（含范德华力）与地球环境（不含范德华力）2种试样的双轴压缩试验,通过试样局部变形、速度场、孔隙比及转动场的变化情况研究了2种试样剪切带形成与发展。研究结果表明,剪切带的形成是试样内部应变局部化的结果,同时也伴随着试样内部颗粒相对转动的局部化,月面环境对试样的破坏形式与性状（剪切带的倾角及厚度）有显著的影响,探月工程中必须考虑到月面环境对月壤力学性质的影响。     

轴对称单轴压缩岩石材料应变局部化剪切带数值模拟研究

利用FLAC3D数值软件中的应变硬化—软化本构模型,模拟轴对称单轴压缩条件下岩石试件应变局部化剪切带的形成过程,同时对不同端面效应下长径比效应对应变局部化剪切带的影响进行分析。结果表明,岩石试件的应变局部化剪切带的启动、发展到最终形成是一个渐进的演化过程。端面效应和长径比效应均影响应变局部化剪切带的形状和岩石强度。光滑边界应变局部化剪切带形成于试件上下端面,摩擦边界形成于试件中部区域,且摩擦边界时试件的峰值强度高于光滑边界,峰值强度后的软化阶段陡峭于光滑边界。随着岩石试件长径比的增加,岩石的峰值强度在逐步降低,达到峰值强度所需的时步却呈增加的趋势。这一研究成果有助于加深理解岩石破坏过程。     

平面应变状态下岩石剪切带网络数值模拟研究

研究了围压、端面约束、试件高度及界面特性对剪切带网络形成的影响。首先对影响剪切带网络形成的因素（包括试件高度、围压及温度、试件端部轴向应变大小的计算方法的差异）进行了分析。之后采用FLAC 3D对平面剪切带网络进行了数值模拟,其中摩擦角及内聚力为应变软化。探讨了倾向于出现剪切带网络的若干条件。若端面缚束较强且存在一定的侧压力,试件中部出现多重剪切带。若端面缚束较弱,试件端部易于形成剪切带网络,且随试件高度的增加剪切带条数增加。界面法向和切向刚度越大,剪切带穿越断层并按其固有方向延伸能力越强,剪切带网络格局越明显。获得的数值结果可在实验中找到佐证,并且可以用来解释地震中的一些剪切应变局部化现象。     

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

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

岩土材料塑性损伤模型及拱坝的变形局部化分析

常规的弹塑性模型难以模拟材料破坏时由于内部损伤累积引起的变形局部化带形成过程。将岩土材料抽象为含有孔洞的球形体胞单元的集合体,采用一个标量--孔洞体积百分比来刻划单元的损伤程度,基于细观损伤力学提出了一个适用于岩土材料的塑性损伤模型。将该模型通过用户子程序嵌入到大型有限元商业软件MRAC中,研究岩土材料的局部化剪切带形成过程。采用该模型,成功地对平面应力条件下岩土材料单轴受压和单轴受拉试件的局部化破坏分别进行了数值模拟。采用该模型,对拱坝在超水载作用下的变形局部化进行了研究,数值计算结果表明,上游坝踵处最先发生变形局部化现象。     

孔隙流体对岩体变形局部化的影响及数值模拟研究

本文首先叙述了孔隙流体对固体变形局部化影响的若干研究进展,并利用三维连续体快速拉格朗日分析,模拟了孔隙压力对岩体变形局部化应变场(剪切带)的影响。数值模拟结果表明,孔隙压力对变形局部化应变场有较大影响,随着孔隙压力的增加剪切带趋于不明显,岩体试件倾向于发生拉伸破坏。      

结构性黄土二元介质本构模型在局部化剪切带中的应用

基于破损力学理论,将结构性黄土抽象成具有一定结构强度的结构块和摩擦带组成的二元介质模型。对结构性黄土体来说,局部化剪切带问题也是土体的破损问题,剪切带萌生发展的实质就是结构块向摩擦带转化的动态过程。应用结构性土的双参数破损率二元介质本构模型,采用数值分析方法模拟了平面应变压缩条件下结构性土中局部化剪切带萌生、扩展的过程,研究了不同缺陷方案下局部化剪切带的形态、特性与规律,发现结构性土中局部化剪切带的发展起初是由一段段不连续的微小局部破坏区域在外荷载逐步作用下渐进扩展连接贯通而形成整体剪切带的破坏形式。二元介质本构模型和常规有限元的结合,形象生动地再现了局部化剪切带萌生、发展的过程。     

基于Cosserat连续体的颗粒破碎影响研究

颗粒破碎对颗粒材料宏观力学行为有重要影响。 结合Hardin的破碎经验公式,将表征破碎程度的破碎参量与Cosserat连续体的内部长度参数相关联,形成一个基于Cosserat连续体且能考虑颗粒破碎的弹塑性模型。数值算例主要考察了颗粒破碎对颗粒材料承载能力、塑性应变及局部化行为的影响,数值结果表明,颗粒破碎主要发生在剪切带内,颗粒破碎使得剪切带明显变窄且剪切带内外等效塑性应变梯度明显增大。     

基于数字图像测量系统的饱和细砂渐进破坏特性研究

基于亚像素角点检测的数字图像测量系统能够记录三轴试样表面方块角点的位移,从而获得试样表面每一时刻的局部应变及应变场。分析了试样不同位置的局部径向应变在剪切过程中的变化,获得了以不同初始成样含水率制备的松砂和密砂试样在不同特征时刻的应力与应变的特征值。通过轴向应变场的变化分析了从应变局部化出现到剪切带发育、形成的这一完整的渐进破坏过程,总结了剪切带形成时的局部最大轴向应变特征值,并定性地分析了剪切带内、外土体在渐进破坏过程中不同的轴向应变增长率。试验结果表明：在应变场中试样应变局部化明显,并可以依此确定应变局部化的出现、剪切带的形成;对于密砂及初始含水率为0%制备的松砂试样,应力在应变局部化出现之初即达到峰值,剪切带形成时应力已经开始下降,进入应变软化阶段;以初始含水率为6%和12%制备的松砂试样在达到应力峰值时剪切带已经形成;剪切带内土体的局部轴向应变增长幅度比剪切带外的土体大得多。试样整体轴向应变的增大主要是由剪切带内土体剪切破坏产生的较大轴向应变所致。     

基于应变梯度理论的韧性剪切带理论研究

地质领域中存在的大量事实都可以归为应变局部化范畴,韧性剪切带是剪切应变局部化带。本文首次应用应变梯度塑性理论对韧性剪切带进行了理论分析。获得了韧性剪切带内部局部剪应变(率)分布规律,剪切带错动总位移及剪切带法向位移的理论表达式,以及韧性剪切带内部体积应变(率)和孔隙率分布规律。剪应变局部化导致了体积应变和孔隙率局部化。随着应变软化的加剧,剪应变、体积应变和孔隙率局部化加剧;随着应力率绝对值的增加,剪应变率、体积应变率集中程度增强。      

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.     

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

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

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.     

Numerical modeling of combined effects of upward and downward propagation of shear bands on stability of slopes with sensitive clay

Modeling of progressive development of zones of large inelastic shear deformation (shear band) that results from strain‐softening behavior of sensitive clays could explain the failure mechanisms of large landslides. Because of toe erosion, a shear band can be initiated and propagated upward (inward) from the river bank. On the other hand, upslope surcharge loading could generate shear bands that might propagate down towards the river bank. In the present study, upward and downward propagation of shear bands and failure of sensitive clay slopes are modeled using the Coupled Eulerian Lagrangian approach in Abaqus finite element (FE) software. It is shown that the formation and propagation of shear bands are significantly influenced by kinematic constraints that change with displacements of the soil masses, and therefore the propagation of an existing shear band might be stopped and new shear bands could be formed. The main advantages of the present FE modeling are: (i) extremely large strains in the shear bands can be successfully simulated without numerical issues, (ii) a priori definition of shearing zones is not required to tackle severe strains; instead, the FE program automatically identifies the critical locations for shear band formation and propagation. Toe erosion could significantly increase the slope failure potential because of upslope surcharge loading. FE analyses with a thick and thin sensitive clay layers show that the global failure could occur at lower surcharge loads in the former as compared to the latter cases. Copyright © 2016 John Wiley & Sons, Ltd.     

Evolution of shear zones using numerical analysis at the May 12th Wenchuan Earthquake Site,China

Landslides and collapses occurred during the May 12th earthquake in Wenchuan County. Most of these landslides and collapses were caused by shear bands. Shear band triggers instability on mountain slopes, resulting in debris flow, landslides, or collapse. According to experimental results, there was only one shear band forming in the soil layer prior to the initiation of debris flow under shear load, although several fine shear bands appeared. The development of shear bands in saturated soils is numerically investigated in this paper using the in situ soil from the Weijia Gully, Beichuan County. The evolution of shear banding from several finite amplitude disturbances (FADs) in pore pressure has been studied. The numerical analysis revealed that the FADs evolved into a fully developed shear band. It is shown that the shear banding process consists of two stages: inhomogeneous shearing and true shear banding.     

Influence of Lode angle‐dependent failure criteria on shear localization analysis in sand

Geotechnical experiments show that Lode angle‐dependent constitutive formulations are appropriate to describe the failure of geomaterials. In the present study, we have adopted one such class of failure criteria along with a versatile constitutive relationship to theoretically analyze the effects of Lode angle on localized shear deformation or shear band formation in loose sand for both drained and undrained conditions. We determine the variation in the possible stress states for shear localization due to the introduction of Lode angle by considering the localized deformation as a bifurcation problem. Further, similar bifurcation analysis is performed for the stress states along a specific loading path, namely, plane strain compression at the constitutive level. In addition, the plane strain compression tests have been simulated as a boundary value finite element problem to see how Lode angle affects the post‐localization response. Results show that the inclusion of a Lode angle parameter within the failure criterion has considerable effects on the onset, plastic strain, and propagation of shear localization in loose sand specimens. For drained condition, we notice early inception of shear localization and multiple band formation when the Lode angle‐dependent failure criterion is used. Undrained localization characteristics, however, found to be independent of Lode angle consideration.     

Use of acoustic emissions to investigate localization in high-porosity sandstone subjected to true triaxial stresses

Acoustic emissions were used to investigate the evolution of damage and strain localization in Castlegate sandstone specimens subjected to true triaxial states of stress, where the intermediate principal stress ranged from equal to minimum compression to equal to maximum compression. Three failure modes were observed: shear band formation at low mean stresses, compaction localization at intermediate stresses, and no localization at high mean stresses. For shear bands, the onset of localization consistently occurred prior to peak stress, while compaction localization initiated at the beginning of a stress plateau. The band angle (defined as the angle between the band normal and the maximum compression direction) determined by fitting a plane through the localized acoustic emission events corresponded well with the shear band angle expressed on the specimen jacket. As expected, the band angle decreased with increasing mean stress. The theoretically predicted dependence of band angle on deviatoric stress state was not confirmed; however, data scatter due to natural variations in material could obscure such a trend. Each failure mode displayed a unique acoustic emission rate response; therefore, this rate response alone can be used to determine the failure mode and the onset of localization in Castlegate sandstone.     

Elasto-viscoplastic finite element study of the effect of degradation on bearing capacity of footing on clay ground

The bearing capacity of footing has been studied by both conventional and numerical methods by many researchers. However, degradation of the microstructure of material, that is, a change in the microstructure of the soil, has not been adequately taken into account. Degradation of microstructure causes strain softening of materials and it leads to strain localization such as shear bands and slip bands. From an engineering point of view the strain localization is crucial because it is a precursor of failure. In the present study, finite element analyses of the bearing capacity of a shallow foundation on homogeneous and inhomogeneous saturated clay strata have been conducted using an elasto-viscoplastic soil constitutive model of microstructure change. A series of analyses of footing on clay deposit with different microstructure parameters have been carried out. Numerical results show that strain localization can be predicted during the loading of rigid footing on highly structured soil and strain localization affects the footing–soil interaction. The effects of footing roughness on the failure mechanism are also discussed in the study.