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

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

基于图像测量的饱和细砂排水剪切试验研究

将计算机数字图像测量技术应用于土工静力三轴试验的土样变形测量,实现了土样变形的非接触测量。针对南宁饱和细砂,采用干装敲击法和控制不同初始成样含水率的湿装夯实法制备松砂和密砂试样,联合采用数字图像测量技术与常规三轴变形量测方法,在不同围压下进行了固结排水剪切试验。着重探讨了成样方法、初始成样含水率、密实程度、固结围压对两种量测方法所获得的应力-应变关系与邓肯-张模型中的切线模量5参数的影响。试验结果表明：成样方法、初始成样含水率、密实程度和围压决定土样的剪切破坏形式,影响应力-应变关系;剪切破坏的形式以及剪切带位置的不同,导致图像测量和常规三轴试验获得的轴向应变和偏应力数值均存在不同程度的差异;邓肯-张模型5参数不仅受到试样的成样方法、初始成样含水率、密实程度的影响,而且受到量测方法的影响。     

真三维状态下土体强度及变形分叉影响

通过引入应力Lode角建立合理的角隅函数描述偏平面形状,对三维Mohr-Coulmb强度准则进行修正,使之能合理反映中主应力比对土体峰值内摩擦角的影响。修正后的准则在偏平面上的强度线与Lade-Duncan准则相似,但更具灵活性,预测结果更符合松砂真三轴试验结果。由于土体在加载过程中常发生应变局部化现象,导致其强度得不到充分发挥。基于非共轴塑性模型的分叉分析研究了应变局部化发生对土体强度发挥的影响,通过与密砂真三轴试验结果对比表明,分叉强度能合理反映在不同中主应力比条件下应变局部化发生对土体强度的降低作用。     

基于三维数字图像相关技术的砂岩渗流-应力耦合变形局部化特性试验研究

为研究渗流-应力耦合作用下砂岩变形局部化破坏特征,将三维数字图像相关（three-dimensionaldigitalimage correlation,简称3D-DIC）技术运用于可视化三轴伺服控制试验系统中,对砂岩开展了不同渗透水压条件下的三轴压缩试验研究。研究结果表明：（1）在渗流-应力耦合作用下,砂岩表面变形场云图经历了从均匀到非均匀的演化过程。在峰前的变形较为均匀,在峰值处出现应变集中,在峰后较短时间内迅速形成变形局部化带;（2）随着渗透水压的增大,岩石的峰值强度、弹性模量呈现指数型非线性衰减,峰值强度弱化系数和弹性模量弱化系数均增大,且渗透率逐渐增大,最大渗透率出现的时间点越早;（3）随着渗透水压的增大,砂岩变形局部化带从剪切状过渡到张拉-剪切状,轴向和径向变形局部化带带内外应变在峰后短时间内出现较大差异,且带内应变远大于带外应变;（4）随着渗透水压的增大,砂岩变形局部化启动应力水平越高,启动的时刻点越早。渗透水压与砂岩变形局部化启动应力水平呈线性相关性,轴向变形局部化启动水平受渗透水压的影响比径向更为显著。     

颗粒物质在等比例应变加载下的分散性失稳模式

由地下水引起的静力液化可能是边坡失稳的隐含机制之一,松砂在不排水剪切条件下可能发生静力液化,密实的颗粒集合体在特定的应变路径下也会出现相似的现象,即试样整体发生急剧的失稳,应力状态尚处于峰值强度线以内。该种失稳模式称为分散性失稳,是为了强调失稳模式中没有出现应变局部化或者剪切带。采用连续-离散耦合分析方法,研究由不规则形状颗粒组成的密实集合体在等比例应变加载路径下的力学特性。根据Hill的材料失稳理论,当试样的应力增量 和应变增量 对应的2阶功 为负时,试样即发生不可逆的整体失稳破坏。以根据不同等比例应变路径得到 曲线为界,在 平面内将试样的应力状态分为剪缩区、剪胀-稳定区和剪胀-非稳定区,连接不同围压下试样发生分散性失稳时的应力状态形成失稳线发现,峰值强度线高于临界状态线,临界状态线高于失稳线。     

减围压平面应变压缩试验条件下丰浦砂中剪切带特性研究

不同于常规土体单元试验中恒定围压的应力路径,进行了减围压（卸载）应力路径下丰浦砂的平面应变压缩试验,研究了初始围压和初始相对密度对砂土的剪切强度和剪切带特性的影响。结合二维数字图像（DIC）相关方法,分析剪切带的形成过程,对剪切带的厚度和倾角及其影响因素进行定量地分析。研究结果表明,剪切带形成于峰值剪切强度之前;对于密实丰浦砂试样,随着围压的增加,其剪切强度(最大有效主应力比)减小且峰值强度出现越晚,剪切带厚度减小而剪切带倾角变化不大;丰浦砂试样的相对密度越高,剪切强度越大且峰值强度出现越早,剪切带厚度越小而剪切带倾角越高。此外,基于峰值内摩擦角的Coulomb公式可以较好地预测平面应变试验条件下丰浦砂的最大剪切带倾角。     

强度参数对初始随机缺陷岩样全部变形的影响

在单轴平面应变压缩条件下, 采用FLAC模拟了初始内聚力及内摩擦角对具有随机材料缺陷岩样轴向、侧向、体积变形及由侧向应变及轴向应变计算得到的泊松比的演变的影响。采用编写的FISH函数于试样内部规定随机缺陷并计算其全部变形特征。密实的岩石服从莫尔库仑剪破坏与拉破坏复合的破坏准则, 破坏之后呈现应变软化—理想塑性行为。随着强度参数(初始内聚力及内摩擦角)的增加, 应力峰值及对应的侧向应变的值提高, 体积应变的峰值增加, 应力峰值所对应的计算得到的泊松比稍有增加。由于被剪切带所分割的毗邻块体之间的相对滑动, 纵然扩容角为零, 试样在峰后的变形阶段, 仍然可以观察到体积膨胀现象, 这与作者的理论分析结果一致。当强度参数降低时, 侧向应变—轴向应变曲线、体积应变—轴向应变曲线及计算得到的泊松比—轴向应变曲线的峰前非线性部分变得短暂。剪切带倾角的数值解的上限低于Coulomb理论, 下限在Roscoe理论附近波动, Arthur理论的预测结果与本文剪切带倾角的数值解更接近.      

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

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

应变软化扩容对含随机缺陷岩石的渐进变形及破坏前兆特征的影响

在单轴平面应变压缩条件下, 采用FLAC模拟了剪切扩容对含随机缺陷岩石破坏前兆及变形特征的影响。密实的岩石服从莫尔库仑剪破坏与拉破坏复合的破坏准则, 破坏之后呈现应变软化-理想塑性行为。缺陷在破坏之后经历理想塑性行为。随着轴向应变的增加, 试样内部破坏的单元数目增加, 直到达到一个常数, 该常数随着扩容角的增加而增加。当扩容角较高时, 计算得到的泊松比在峰前就可以超过0.5;剪切扩容于峰前发生; 变形后试样的最终体积大于初始体积。剪切局部化(导致了毗邻块体之间的相对滑动)及剪切扩容(发生于剪切带内部)是非零扩容角试样峰后体积膨胀的原因。在峰前, 通过观察剪切应变增量、破坏的单元数目、侧向应变、计算得到的泊松比及体积应变可以发现, 扩容角越高, 试样破坏的前兆越明显。在低扩容角时, 由于弯曲的剪切带边界, 试样内部充分发展的剪切带的倾角比较分散, 剪切带的倾角更接近Arthur理论。      

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

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

Effects of Previous Deviatoric Strain Histories on the Undrained Behaviour of Hostun RF Loose Sand

The effects of previous deviatoric strain histories on the undrained behaviour of very loose and saturated Hostun RF clean sand are investigated in this paper. From an initial isotropic stress state in the triaxial plane, strain histories are generated by isotropic consolidation followed by standard drained triaxial preshear or presheared cycle, either in compression or in extension, up to a desired value of axial strain or mobilized stress ratio. Deviatoric strain histories are achieved by having nearly the same void ratio at the beginning of the undrained shearing for all tested samples. Subsequent undrained behaviour in triaxial compression and extension is analyzed in detail. Previous deviatoric strain histories can progressively transform the compressive and unstable behaviour of loose sand into a dilative and stable behaviour of dense-like sand, while being loose. Experiments show a common response induced by recent strain histories in terms of effective stress paths, independently of the axial strain attained during the drained presheared cycle, a unique initial gradient of the effective stress paths, a progressive appearance of dilatancy, an evolution the undrained behaviour and a systematic partial static liquefaction associated with softer behaviour when sheared in the opposite direction of the initial presheared direction. This paper offers a comprehensive understanding of the mechanisms of a specialized induced anisotropy created by simple linear stress paths in the classical triaxial plane.     

Undrained shear band in water saturated granular media: A critical revisiting with numerical examples

By taking into account the continuity of stresses and pore pressure across the boundaries of a shear band, it is theoretically shown that incipient localization may take place in saturated loose contractive sand in the regime of decreasing deviator stress under undrained conditions. The undrained shear band orientation primarily depends on the Poisson's ratio and the dilatancy characteristics of the material. However, under strict isochoric constraint, localization is precluded in dense dilative sand with deviator stress increasing only. If any local volume change takes place with a certain mechanism, undrained shear band would be possible in dense sand, most likely in the regime of decreasing friction angle after the peak friction angle is mobilized. Numerical examples are given to demonstrate the variation of shear band inclinations with void ratios and initial consolidation pressures. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

常剪应力路径下含气砂土的三轴试验

天然气水合物完全分解时,产生的气体使得能源土孔隙压力急速增加,有效应力减小,进而引起土体液化破坏。此时深海能源土斜坡的应力状态与静力液化失稳过程可简化为含气土在常剪应力排水(或不排水)应力路径下的破坏问题。以此为背景,提出了制备含气砂土试样的改进充气管法,并开展了含气砂土的常剪应力路径三轴试验。22组试验结果表明:同一孔隙比的含气密砂在不同围压与常剪应力下具有相同的失稳线;含气砂土试样失稳时的应力比和体变均随初始相对密实度的增大而增大;含气密砂在常剪应力路径下饱和度对失稳特征影响的规律性在排水与不排水条件下均不明显,但在不排水条件下含气砂土的孔压(或体变)对变形的敏感性降低;含气密砂在常剪应力路径到达失稳点之后,排水条件下是瞬变的液化鼓胀破坏,不排水条件下是渐变的剪切破坏。     

Exploring the undrained induced anisotropy of Hostun RF loose sand

The effects of recent history on the undrained behaviour of very loose and saturated Hostun RF clean sand are investigated in this paper. From an initial isotropic stress state in the triaxial plane, recent histories are generated by isotropic consolidation followed by standard drained triaxial preshear in compression, up to a desired value of axial strain or mobilized stress ratio, and unloading to an initial stress ratio. Subsequent undrained behaviour in triaxial compression is analysed in detail. This paper contributes to the traits explaining the progressive transformation of a compressive and unstable behaviour of loose sand into a dilative and stable behaviour of dense-like sand by previous history, while remaining in the same state of loose density. Experiments show a large pseudo-elastic domain induced by recent history in terms of effective stress paths, function of the initially mobilized stress ratio level, a unique initial gradient of the effective stress paths depending on the stress ratio at the beginning of the undrained shearing, a progressive appearance of dilatancy and a surprising evolution the undrained behaviour of loose sand. Experimental results evidence the important role of the recent deviatoric strain history, from any initial isotropic or anisotropic stress state. This paper also offers a comprehensive understanding of the history mechanisms created by simple linear stress paths with fixed direction in the classical triaxial plane.     

Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Non‐linear analysis of shear band formation in sand

The paper investigates the incipience of shear band with an incrementally non‐linear constitutive equation. Necessary conditions for the emergence of shear band are derived. The lower bound solution is obtained by taking the strain rate inside and outside the shear band into consideration. Numerical results of localized bifurcation for general stress and strain are presented and compared with experiments. In the principal stress space, the stresses at the onset of shear band form a surface, which is partially enclosed by the failure surface for homogeneous straining. The significance of the analysis for identification of the material parameters and verification of the constitutive model against experiments is discussed. Copyright © 2000 John Wiley & Sons, Ltd.