基于弹性应变能的岩石强度准则

为建立更符合岩石屈服与破坏机制的强度准则,基于能量转化是岩石屈服与破坏的本质属性,采用试验与理论分析相结合的方法,对岩石屈服与破坏准则进行了研究。以岩石强度与整体破坏准则为基础,通过引入弹性应变能释放分散系数,建立基于弹性应变能强度准则;分别采用M-C准则、Murrell准则、三剪能量准则、统一能量准则、三维H-B强度准则及基于弹性应变能岩石强度准则对盐岩和花岗岩的破坏强度进行了计算。结果表明,基于弹性应变能岩石强度准则的计算结果与试验值比较吻合（尤其是真三轴试验条件下）,并且分析了产生上述结果的内在机制。所建立的强度准则仅需测定常规岩石力学参数（单轴抗压强度与泊松比）,物理力学意义明确,对于定量描述岩石的屈服与破坏特性具有重要的意义。     

大型地下洞室高边墙的拉应变准则研究

高地应力等复杂环境下大跨度地下洞室群开挖容易诱发高边墙的脆性拉裂破坏,在Mohr-Coulomb准则基础之上,基于拉应变准则,提出大型地下厂房高边墙破坏的复合拉剪屈服准则。根据弹塑性理论,张拉和剪切屈服分别采用相关联和非相关联的塑性流动法则,详细推导了该准则在有限差分法中的计算格式,并对张拉和剪切屈服面的应力空间进行了分界。利用C语言程序编写的动态链接库文件,实现了该准则在通用软件的科学仿真计算。以高应力地区的拉西瓦水电站地下厂房为工程实例,与FLAC3D自带的拉剪屈服准则相比,基于拉应变准则的计算结果表明,大跨度高边墙拉破坏明显,塑性区以及应力松弛区的面积更小,能合理描述硬岩高边墙的脆性拉裂破坏特征。大型硬岩地下洞室群的计算分析验证了该准则的合理性和实用性     

平面应变条件下M-C材料屈服时的中主应力公式

基于Mohr-Coulomb准则,通过引入罗德角参数及Mohr-Coulomb准则平面应变等效广义Mises准则,推导了平面应变条件下材料屈服时的中主应力公式,并给出其与基于SMP准则（Spatial Mobiliaztion Piane）的平面应变中主应力公式之间的关系式。对平面应变条件下的条形地基进行了有限元数值模拟,计算结果验证了所推导的平面应变条件下M-C材料屈服时的中主应力公式的正确性,其公式为岩土体的强度和变形分析、中主应力的强度效应的进一步研究,以及Mohr-Coulomb准则的修正提供了有力的基础。     

考虑张剪破坏的边坡下限原理有限元法

在工程中发生的很多边坡破坏模式不仅表现为剪切破坏,还表现为在边坡后缘产生一定深度的张拉破坏。这种既考虑剪切破坏又考虑张拉破坏的屈服准则,是经典的Mohr-Coulomb（简称M-C）屈服准则难以模拟的。从空间方位离散的角度出发,对经典的M-C屈服面线性化的方式进行改造,建立了基于方位离散的线性化屈服准则,同时在屈服准则中考虑了张拉破坏准则,并引入伪黏聚力,保证在每一个离散方位平面上既不违背屈服准则,又不违背张拉破坏准则,建立了同时考虑张拉及剪切破坏的下限原理有限元法。算例表明了该方法的有效性,同时也证明了仅考虑剪切破坏时,会过高地估计边坡的安全系数。     

A yield criterion for isotropic and cross‐anisotropic cohesive‐frictional materials

This paper proposes a yield and failure criterion for cohesive and frictional materials. The function is given by the combination of a Lode dependence for the behaviour in the deviatoric plane and a meridian function for the pressure‐dependent behaviour. A variety of shapes can be achieved with the proposed criterion including Lode dependences which are able to reproduce the behaviour of isotropic and cross‐anisotropic materials in the deviatoric plane. The criterion is validated through the comparison with experimental data based on multiaxial experimental tests on clays, sands, rocks and concrete. Finally, the convexity of the criterion is analysed and discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

基于幂强化本构模型的轴对称圆巷弹塑性解

轴对称圆巷的弹塑性求解的关键是选择合适的屈服准则。已经有诸多学者选择Mohr-Coulomb准则、Drucker-Prager准则和Hoek-Brown准则等,进行了相应的求解。为了探讨更符合工程实际需要的准则和求解,在考虑岩石材料的应变强化效应的条件下,建立了轴对称圆巷的幂强化本构模型和基于Drucker-Prager屈服准则的幂强化-理想塑性模型,并进行了弹塑性求解。以工程实例为计算条件,将幂强化-理想塑性模型的计算结果与基于Mohr-Coulomb准则、Drucker-Prager准则的理想塑性模型和幂强化模型的计算结果分别进行了对比,分析幂强化参数对围岩弹塑性解的影响。研究表明,应变强化效应对围岩稳定性有较大影响,对于应变强化效应较强的岩石材料,采用幂强化模型分析更接近工程实际。     

考虑渗透力和原始Hoek-Brown屈服准则时圆形洞室解析解

采用原始Hoek-Brown非线性屈服准则,推导了渗透力作用下圆形洞室弹塑性解析表达式。根据解析公式,绘制了渗透力作用下基于原始Hoek-Brown屈服准则的围岩特性曲线、塑性区半径与洞壁支护力关系曲线、塑性区半径与围岩自重应力关系曲线,并与基于Mohr-Coulomb屈服准则以及不考虑渗透力的图形进行了对比研究,研究表明：应用原始Hoek-Brown屈服准则求出的洞室径向位移和塑性区半径都大于应用Mohr-Coulomb屈服准则求出的结果;在围岩情况较好的条件下基于原始Hoek-Brown屈服准则计算的塑性区半径远大于基于Mohr-Coulomb屈服准则算出的值;渗流效应对地下洞室的塑性区半径和洞壁径向位移的影响十分显著。     

Mises屈服准则及其在钻探工具强度设计中的应用

在地质钻探工具设计和使用过程中,需要用到强度理论或屈服准则对遇到的钻采设备、工具失效等问题进行分析。为帮助人们对Mises屈服准则（应力）的理解,从理论上阐述了Mises屈服准则与畸变能密度理论的关系,其实质就是第四强度理论。在弹塑性范围内,Mises屈服准则统一了单向拉压状态的第一强度理论、复杂应力状态的第三强度理论和第四强度理论。Mises屈服准则在有限元数值分析中获得了更为广泛的应用。     

极限分析有限元法讲座--Ⅱ有限元强度折减法中边坡失稳的判据探讨

边坡失稳,滑体滑出,滑体由稳定静止状态变为运动状态,同时产生很大的且无限发展的位移,这就是边坡破坏的特征。有限元中通过强度折减使边坡达到极限破坏状态,滑动面上的位移和塑性应变将产生突变,且此位移和塑性应变的大小不再是一个定值,有限元程序无法从有限元方程组中找到一个既能满足静力平衡又能满足应力-应变关系和强度准则的解,此时,不管是从力的收敛标准,还是从位移的收敛标准来判断有限元计算都不收敛。塑性区从坡脚到坡顶贯通并不一定意味着边坡破坏,塑性区贯通是破坏的必要条件,但不是充分条件,还要看是否产生很大的且无限发展的塑性变形和位移,有限元计算中表现为塑性应变和位移产生突变。在突变前计算收敛,突变之后计算不收敛,表征滑面上土体无限流动,因此可把有限元静力平衡方程组是否有解,有限元计算是否收敛作为边坡破坏的依据。-     

基于非线性破坏准则的边坡稳定性极限分析

上限定理是求解岩土工程问题的有效工具。以上限定理为理论基础,分析边坡的稳定性问题,并考虑了岩土材料破坏准则的非线性特性。在非线性Mohr-Coulomb破坏准则下,采用条分法与极限分析上限法相结合的方法,对边坡稳定性进行分析。通过切线法引入非线性强度参数 、 ,推导了岩土材料在非线性破坏准则下的相关联流动法则,建立功能方程,推导了边坡直线滑裂面、折线滑裂面和光滑曲线滑裂面安全系数F的计算方程。采用数学规划法计算后与工程实际中常用的边坡稳定性分析方法进行对比,并获得安全系数F =1.0时的稳定性系数Ns。与已有的研究成果进行比较分析,结果表明了该方法的正确性及优越性     

Further Development of the Three-Dimensional Hoek–Brown Yield Criterion

The original Hoek–Brown failure criterion was first introduced in 1980. Since then, it has been further developed and used in different rock mechanics applications. The main advantage of the Hoek–Brown criterion is that its input data are determined from the uniaxial compressive strength, mineralogy and structural properties of the rock. However, its application in the petroleum industry has been limited, partly because it does not take account of the intermediate principal stress in a three-dimensional stress state. In 2005, Priest introduced comprehensive and simplified versions of a new three-dimensional Hoek–Brown yield criterion, which was developed by combining the two-dimensional Hoek–Brown and three-dimensional Drucker–Prager criteria. Until now, direct analytical solutions have not been developed to solve the governing equations and numerical solution strategies have been applied. This paper presents and validates an explicit solution for the comprehensive three-dimensional Hoek–Brown yield criterion, which can be combined with numerical stress modelling codes.     

平面应变状态下边坡安全系数的近似求解

Mohr—Coulomb屈服准则不考虑中间主应力对材料强度的贡献,是一个保守的强度屈服准则,在近似平面应变状态下,匹配于Mohr-Coulomb屈服准则的Drucker-Prager屈服准则能给出边坡安全系数偏保守一些的估计;相对而言,Matsuoka—Nakai和Lade-Duncan屈服准则均考虑了中间主应力对土体强度的影响,在近似平面应变状态下,匹配于Matsuoka-Nakai和Lade-Duncan屈服准则的Drucker-Prager屈服准则能给出岩土类材料较为合理的破坏强度预测,土体的抗剪强度能被更有效地发挥出来,因而借助这些转化的Drucker-Prager屈服准则能给出平面应变状态下土坡安全系数较为合理和真实的预测。文中的边坡稳定算例证实了这些观点。     

Single hardening constitutive model for frictional materials III. Comparisons with experimental data

A unified constitutive model for the behavior of frictional materials is described. The model is based on concepts from elasticity and plasticity theories. In addition to Hooke's law for the elastic behavior, the framework for the plastic behavior consists of a failure criterion, a nonassociated flow rule, a yield criterion that describes contours of equal plastic work, and a work-hardening/softening law. The functions that describe these components are all expressed in terms of stress invariants. The model incorporates twelve parameters which can all be determined from simple experiments such as isotropic compression and conventional triaxial compression tests. Validation of the model is achieved by comparison of predicted and measured stress-strain curves for various two- and three-dimensional stress-paths obtained for different types of frictional materials.     

Kinematical analysis of 3D passive earth pressure with nonlinear yield criterion

Conventional methods for calculation of passive earth pressure were mainly based on the assumptions of the linear Mohr‐Coulomb yield condition and plane strain failure mechanism. However, both theoretical and experimental studies have shown that such assumptions are not satisfied in some geotechnical projects. Herein, a novel method incorporating a kinematically admissible 3‐dimensional (3D) rotational failure mechanism and the nonlinear power‐law yield criterion is proposed to compute the passive earth pressure acting on the inclined retaining walls. Instead of using the nonlinear yield criterion directly, a straight line tangential to the nonlinear yield curve is employed to represent the strength of soils, and therefore, the nonlinear problem is transformed into the traditional linear problem. The 3D failure mechanism is generated through rotating a circle defined by 2 log‐spirals, and a plane strain block is inserted into the mechanism to consider the retaining walls with different widths. Earthquake effects are taken into account by using quasi‐static representation, and the horizontal seismic coefficient concept is adopted for the estimation of passive earth pressure under seismic conditions. An analytical expression about the 3D passive earth pressure is educed by means of the upper bound theorem of limit analysis. Numerical results for different practical parameters are obtained from an optimization scheme where the minimum of passive earth pressure is sought. Compared with available 2‐dimensional and 3D solutions, the proposed method is validated. A parametric study is conducted to investigate the effects of different parameters on the 3D static and seismic passive earth pressure.     

基于离散元的重力坝多滑面深层抗滑稳定分析

重力坝的深层抗滑稳定分析多采用刚体极限平衡法和非线性有限元法,但刚体极限平衡法不能反映坝基岩体渐进失稳过程和破坏的力学机制,非线性有限元法模拟岩体不连续效率低,且还未有统一的失稳判据。基于离散块体边界应力计算结果,提出重力坝坝基多滑面抗滑稳定安全系数计算公式,并探讨了以坝体-坝基系统能量突变作为坝基失稳判据的物理意义。通过与非线性有限单元法、刚体极限平衡法计算结果的比较,验证了所提方法及判据的可靠性及合理性。结合向家坝水电站重力坝泄12#坝段进行计算分析,结果表明：基于离散元强度储备系数法搜索的坝基失稳通道由坝基岩体屈服区和结构面滑移带组成,所提计算方法与判据适用于实际工程的坝基深层抗滑稳定计算分析,且计算结果偏安全。     

Multi-surface Cyclic Plasticity Sand Model with Lode Angle Effect

A Drucker-Prager J ₂ multi-surface-plasticity sand model is modified to employ the Lade-Duncan failure criterion as the yield function. This function includes the first and third stress invariants to account for the dependence of cyclic shear stress–strain behavior on confining pressure and the Lode angle. Related modifications to the flow rule and hardening rule are described. Dependence of dilatancy on confinement is also included. Salient features of the model performance are presented under general three-dimensional (3D) loading conditions, where the yield function provides a more accurate representation of nonlinear shear response. Dynamic response analyses of a mildly inclined infinite slope are performed to illustrate the influence of excitation direction on the accumulation of liquefaction-induced lateral ground deformation.     

Application of a unified linear yield criterion in plane strain to study the strength theory effect in geotechnical engineering

The plane strain condition is a common, but polyaxial stress state for geotechnical structure designs, in which the selection of an appropriate yield or failure criterion is crucial to reasonably account for the intermediate principal stress. Under plane strain condition, a unified linear yield criterion for seven commonly used geotechnical yield criteria is presented in conjunction with the inductive method. These seven yield criteria considered in this study are the Mohr–Coulomb, Tresca, Drucker–Prager, Mogi–Coulomb, Extended Matsuoka–Nakai, Extended Lade–Duncan criteria, and the Unified Strength Theory. The generalized analytical solutions for earth pressure of retaining walls, critical load of strip foundations as well as stress and displacement of circular tunnels are derived on the basis of the proposed unified yield criterion, and their respective theoretical significance is analyzed. Thereafter, the critical load of strip foundations obtained herein is compared with two numerical results from the literature. Furthermore, the effect of strength theory on result differences of the three typical geotechnical problems by simply selecting constants, which conform to different yield criteria, is explored through a parametric study. It is found that the proposed unified yield criterion is convenient for investigating analytical solutions of the aforementioned geotechnical structures. The strength theory effect due to adopting different yield criteria is considerably significant, which cannot be ignored. Additionally, recommendations are provided on how to make use of these seven yield criteria for an optimum design.     

A new rock mass failure criterion for biaxial loading conditions

To simulate brittle rocks, a mixture of glastone, sand and water was used as a model material. Thin galvanized sheets of thickness 0.254 mm were used to create joints in blocks made out of the model material. To investigate the failure modes and strength, both the intact material blocks as well as jointed model material blocks of size 35.6 × 17.8 × 2.5 cm having different joint geometry configurations were subjected to uniaxial and biaxial compressive loadings. A new intact rock failure criterion is proposed at the 3-D level. This criterion is validated for biaxial loading through laboratory experimental results obtained on intact model material blocks. Results obtained from both the intact and jointed model material blocks are used to develop a strongly non-linear new rock mass failure criterion for biaxial loading. In this failure criterion, the fracture tensor component is used to incorporate the directional effect of fracture geometry system on jointed block strength. The failure criterion shows the important role, the intermediate principal stress plays on rock mass strength.     

基于系列Drucker-Prager破坏准则评述土坡的稳定性

各种在? 平面上等效面积的Drucker-Prager屈服准则,如DP4(Drucker-Prager与Mohr-Coulomb两屈服准则在? 平面上的投影面积相等)、DP5(Drucker-Prager与Matsuoka-Nakai两屈服准则在? 平面上的投影面积相等)和DP6(Drucker-Prager与Lade-Duncan两屈服准则在? 平面上的投影面积相等)等模型都不同程度地模拟了平面应变条件下土体的破坏强度,这就是为什么等效面积的Drucker-Prager屈服准则在二维土坡稳定分析中被推广应用的原因。结合砂土的真三轴强度破坏试验结果,从量化角度论证了DP4模型求出的土坡安全系数可能会保守一些,而DP5和DP6两模型求出的土坡安全系数可能会更真实一些,能更好地发挥材料的承载潜力,因而进一步改善了土坡稳定安全系数的计算精度。     

Failure Loci of Some Igneous and Metamorphic Rocks

Summary Experimental evidence from true triaxial tests on dense rocks are analysed with emphasis on the failure modes of these materials under multiaxial loading, ambient temperature and external pressure. The strong dependence of the modes of fracture on the secondary components of applied stresses, and especially on the intermediate principal stress, indicated that the failure surface of these brittle materials may be appropriately described by a failure tensor polynomial criterion. As such, the elliptic paraboloid failure criterion was found to conveniently describe their mode of failure, by considering also the severe influence of anisotropy of the material.  For this purpose, a method developed recently (Theocaris and Panagiotopoulos, 1995a, 1995b) was applied, defining anisotropic hardening plasticity through an appropriate sequence of anisotropic elasticity problems. Assuming a particular path of loading or unloading, we measured the instantaneous tension and compression yield stresses along the transient principal-stress directions. These parameters completely define the instantaneous state of anisotropy of the body for the corresponding loading step, by applying the theory of the elliptic paraboloid failure locus (EPFS) (Theocaris, 1989a). A parameter identification problem was formulated on the constitutive expressions for this most general failure criterion. Then, by applying convenient constraints derived from the EPFS theory, which serve as filters throughout the whole procedure, the characteristic values of terms defining the variable components of the failure tensor polynomial were calculated, as the material was continuously loaded from the elastic into the plastic region and up to the ultimate failure load. Accurate simple tests in uniaxial tension and compression provided sufficient data for the definition of the yield loci of the material, at the considered loading step. These tests may be complemented with biaxial and triaxial modes of loading of the specimens. The results improve the accuracy and sensitivity of the method. All such data were used as input values, for establishing the mode of plastic deformation of the body during particular loading paths.  Moreover, the method employed allows the complete definition of the components of the failure, H, and the strength differential effect, h, tensors at each loading step. These quantities define completely the failure tensor polynomial for each material. Therefore, it presents the important advantage over other experimental methods by clearly indicating the parts contributed to the failure mode (either by plasticity, or by the strength differential effect) and their evolution during plastic deformation.  As convenient prototype materials for testing the method, specimens of metamorphic rocks such as Westerly granite (G), or quartzite (Q) were selected. Interesting results concerning the mechanical and especially the failure modes of such materials were obtained. Furthermore, the mechanical tests indicated clearly some basic properties of these materials as concerns the mode of their structure.