裂隙岩体局部化破坏分析的多重势面分岔模型

采用多重势面弹塑性分岔理论对裂隙岩体的局部化破坏进行计算模拟．基于多重势面弹塑性理论分析局部化问题，构造了适用于裂隙岩体破坏的多重势面不连续分岔模型，并使用数值方法求解局部化方阶在有限元方法的基础上，使用该模型计算裂隙岩体的局部化破坏条带．算例分析表明这一模型用于分析裂隙岩体的局部化破坏是有效的．     

深部岩体损伤对分区破裂化效应的影响

深部岩体是一种具有初始损伤的非连续介质，其内部存在着大量节理和裂隙。深部岩体分区破裂化效应与节理、裂隙的扩展、连接和汇合密切相关。由于节理、裂隙的扩展、连接和汇合可能会导致深部岩体的内部空间由欧氏几何空间向非欧几何空间转化。本文利用损伤变量表示含节理和裂隙岩体的损伤，根据非欧几何模型，获得了静水压力和非静水压力情况下深部圆形洞室损伤围岩的应力场，确定了损伤变量对圆形洞室围岩的应力场和分区破裂化效应的影响。采用Hoek-Brown准则，获得了静水压力和非静水压力情况下损伤围岩的破裂区与非破裂区的分布规律。本文的研究结果为理解深部岩体分区破裂模式提供了有价值的参考。     

不连续岩体变形破坏过程中AE的空间演化与意义

通过一些典型节理化岩体模型变形破坏过程中的AE事件的定位分析和破裂面形态的研究,探讨了不连续岩体变形破坏过程的空间演化特征。结果表明：单节理岩体变形破坏的空间动态特征受节理角控制,随着节理角从小变大,岩体的变形破坏方式从AE在局部集中的突发式失稳向AE弥漫于整个节理层面的类似渐进式的破坏形式发展；失稳破裂面从简单的整体状岩-膏界面变成形态复杂的层内复合破裂面；复杂节理岩体的变形破坏受一条（组）主导构造控制,强AE事件集中发生在不连续构造交汇部位,在包含多条主构造的岩体模型中最终失稳破裂面只与最后阶段AE的位置一致。     

非贯通裂隙岩体损伤演化率相关性及变形特征

含非贯通裂隙岩体是自然界中岩体的主要赋存形式，其裂隙几何特征对岩体的强度及变形均产生显著影响。应变率对岩体的损伤演化及黏滞效应也具有显著的率相关性。首先，运用模型元件的方法，将非贯通裂隙岩体动态破坏过程视为具复合损伤、静态弹性特性、动态黏滞特性的非均质点组成，对黏弹性响应的Maxwell体进行改进，将细观损伤体与裂隙损伤演化的宏观损伤体根据等效应变假设并联组成宏细观复合损伤体，构建综合考虑岩体宏细观缺陷的动态损伤模型；其次，基于断裂力学及应变能理论，对岩体宏观裂隙动态扩展的能量机制进行分析，综合考虑初始裂隙应变能、裂隙动态损伤演化过程应变能、裂隙闭合应变能，得到裂隙岩体宏观动态损伤变量计算公式；最后，将模型计算结果与实验结果进行比较，模型计算结果与实验结果吻合较好，证明了模型的合理性，同时利用模型讨论了裂隙倾角、应变率、岩石性质对岩体变形特征的影响规律。     

基于TCK模型的非贯通节理岩体动态损伤本构模型

提出在岩体动态损伤本构模型中应同时考虑宏、细观缺陷；基于能量原理和断裂力学理论推导得出了同时考虑节理几何及力学特征的宏观损伤变量(张量)的计算公式；基于综合考虑宏、细观缺陷的复合损伤变量(张量)及完整岩石动态损伤Taylor-Chen-Kuszmaul(TCK)模型，建立了相应的单轴压缩下节理岩体动态损伤本构模型；利用该模型讨论了节理内摩擦角及节理长度对岩体动态力学特性的影响规律。研究表明，试件动态峰值强度随着节理内摩擦角的增大而增大，随着节理长度的增加而减小。     

隧道围岩强度不均地段塌方成因及其处理方法

采用微面模型理论和损伤力学方法，建立了节理岩体的弹塑性损伤耦合微面模型. 在节理岩体的微 面上，将岩体视为由节理面与岩石组成的二元介质，以节理连通率作为岩体沿该方向的面积 损伤变量，考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制，基于塑性 理论建立了节理岩体的微面塑性损伤增量本构关系. 采用微面物理量与宏观物理量的几何约 束模型，根据微面方向积分导出了节理岩体的宏观弹塑性增量本构关系. 编制了节理岩体微 面模型的MARC有限元子程序，对节理岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳 定问题进行了数值模拟研究. 数值计算结果表明，该模型能很好地揭示载荷作用下节理岩体 的各向异性非弹性变形和次生节理演化过程.     

深部节理岩体塑性损伤耦合微面模型

采用微面模型理论和损伤力学方法,建立了节理岩体的弹甥性损伤耦合微面模型.在节理岩体的微面上,将岩体视为由节理面与岩石组成的二元介质,以节理连通率作为岩体沿该方向的面积损伤变量,考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制,基于塑性理论建立了节理岩体的微面塑性损伤增量本构关系.采用微面物理量与宏观物理量的几何约束模型,根据微面方向积分导出了节理岩体的宏观弹颦性增量本构关系.编制了节理岩体微面模型的MARC有限元子程序,对节珲岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳定问题进行了数值模拟研究.数值计算结果表明,该模型能很好地揭示载荷作用下节理岩体的各向异性非弹性变形和次生节理演化过程.     

节理裂隙岩体损伤力学研究中的若干问题

节理裂隙岩体损伤力学的研究与应用有着很好的发展前景。但作为一种经历并隐含了复杂的应力、变形及损伤历史的地质体，节理裂隙岩体在物理、力学特征上较一般工程材料具有更显著的随机性、特殊性和复杂性，因而，在其损伤力学研究中有不少基本的、重要的问题还没有解决。本文根据笔者近年来的研究体会，列举了该研究领域的若干重要问题，并就其可能的解决途径与发展趋势提出了一些初步的看法。损伤力学，节理岩体，岩体力学，损伤，断裂     

单轴压缩下断续节理岩体动态损伤本构模型

断续节理将对工程岩体的强度及变形等力学特性产生显著影响，损伤力学中视节理为岩体的一种宏观损伤，因而采用损伤张量来刻画其对岩体的影响。目前学术界提出了用节理的几何、强度及变形等3类参数来描述节理的物理力学性质，而目前的岩体损伤张量计算方法都只涉及前2类参数，均没有涉及其变形参数即法向及切向刚度。为此，在前人研究的基础上，基于断裂及损伤理论提出了考虑节理法向及切向刚度的单轴压缩下单条断续节理引起的损伤张量计算公式，进而通过考虑节理间相互作用给出了单组单排或多排节理岩体损伤张量计算公式。其次，以岩石细观动态损伤模型为基础，结合宏细观损伤耦合观点提出了一个能够同时考虑节理几何、强度及变形参数的断续节理岩体动态损伤本构模型。最后，利用该模型讨论了节理参数及载荷应变率等对岩体动态力学特性的影响，认为节理长度减小及摩擦角增大将导致岩体动态峰值强度及弹性模量增大；岩体动态峰值强度及弹性模量则随着节理法向及切向刚度的增大分别减小或增大；而当节理法向及切向刚度按照同一比例增大时，岩体动态峰值强度及弹性模量则是增大的。岩体动态峰值强度与载荷应变率呈正相关。     

国家自然科学基金重大项目“材料损伤、断裂机理和宏微观力学理论”子课题变形与损伤的局部化理论（1994—1996年工作总结）

国家自然科学基金重大项目“材料损伤、断裂机理和宏微观力学理论”子课题变形与损伤的局部化理论（１９９４—１９９６年工作总结）１损伤演化的物理方程现有的大多数损伤演化方程均不是基于演化的物理机制，而是一种经验形式，因此讨论损伤局部化时只能就事论事．为了使...     

Relationship of crack fabric tensors of different orders

A lower-order crack fabric tensor can be determined accurately by a higher-order one, but in general the reverse does not hold. In this paper, the approximate dependence of a higher-order fabric tensor on a lower-order one is established based on the properties of orientation distribution functions (ODF). As a demonstration of its application, the approximate relationship is used to simplify the fabric-tensor dependent compliance increment due to the presence of the cracks.     

基于Oda裂隙结构张量的流纹岩各向异性试验研究

岩石内天然存在长度、倾角和形态不同的裂隙,造成岩石的各向异性特征。为揭示岩石内天然随机裂隙发育特征对岩石物理力学特性的影响规律,以泥巴山隧址区采集裂隙性流纹岩为研究对象,首先对试样裂隙进行素描统计分析;然后基于Oda裂隙结构张量,获得天然随机分布裂隙的几何统计参数;最后对裂隙性流纹岩试样分别进行单轴和常规三轴压缩试验,得到不同应力路径下流纹岩的应力-应变曲线及物理力学参数。分析Oda裂隙结构张量定义的各向异性参数与试验获得的力学参数之间的规律,研究结果表明:(1)Oda裂隙结构张量适用于天然随机分布裂隙的几何统计分析,各向异性参数A(F)越大,裂隙优势方向越明显;(2)单轴压缩下,随着各向异性参数I1和A(F)的增大,流纹岩各向异性程度增大,弹性模量减小,泊松比增大;(3)常规三轴压缩下,流纹岩弹性模量和泊松比随各向异性参数改变的规律较不明显,Oda裂隙结构张量不再适用。     

A fracture damage model for Jointed rock masses and its application to the stability analysis of dam abutments

Based on continuum damage mechanics, for jointed rock masses, a fracture damage model is presented in this paper. First, the damage tensors are defined through the elastic-flexibility of intact rock and the equivalent elastic-damage flexibility for rock mass. Then, by the self-consistent principle of solid mechanics, the equivalent elastic-damage flexibility tensors involving the interaction between multicracks are deduced. The damage evolution law is proposed involving the mechanism of crack propagation process: frictional sliding, crack kinking, growing of branched tension cracks, interlinking of the microcracks near branched crack tips leading to the breakthrough of macro-cracks and finally the failure of rock mass. Thus the evolution of damage variables reasonably unified with the process of crack propagation is given. Finally, a plastic-brittle damage constitutive relation including brittle coupled strain rate, developed and applied to the stability analysis of complicated rock foundation of a dam in China, is described in this paper.     

岩石材料动态破坏的近场动力学方法数值模拟

裂纹的聚集、扩展、分叉是岩石等脆性材料破坏失效的关键因素，本文在验证了近场动力学方法在研究岩石类材料裂纹动态扩展方面的有效性之后，采用近场动力学方法分别对冲击载荷作用下含有双裂纹岩石材料和单轴压缩作用下含单斜裂纹的岩石材料进行数值模拟．结果表明，对于双垂直裂纹，其裂纹扩展路径大致与预制裂纹成70°夹角；对于单裂纹，裂纹的扩展路径随裂纹倾角的变化而变化，最终导致构件的整体破坏．数值模拟结果表明近场动力学方法可以很好地模拟岩石等脆性材料的裂纹扩展直至破坏的过程，反映裂纹扩展的物理机理；其作为一种新的基于非局部理论的数值方法，在地下岩体工程方面及页岩气的开采方面会有很好的发展前景．     

Coupled modeling of damage growth and permeability variation in brittle rocks

This paper presents a coupled model for anisotropic damage and permeability evolution by using a micro–macro approach. The damage state is represented by a second order tensor. The evolution of damage is determined from a crack propagation criterion. The free enthalpy function of cracked material is obtained by using micromechanical considerations. It is assumed that cracks exhibit normal aperture which is coupled with the crack growth due to asperities of crack faces. By using Darcy’s law for macroscopic fluid flow and assuming laminar flow in microcracks, the overall permeability of the RVE is obtained by a volume averaging procedure taking into account crack aperture in each orientation.     

基于正则结构的各向异性损伤演化律

在Rice的正则结构框架下，推导出基于共轭力的各向异性损伤演化律。其中损伤变量采用二阶裂隙张量，它是固体内微裂纹的一个宏观测度。推导过程不涉及自由能的具体形式，主要结果包括损伤势函数及演化方程的解析表达式。在唯象的损伤力学模型里，损伤演化方程经常以唯象方程的形式出现。研究了唯象方程成立的条件及损伤特征张量的解析表达式。引入了广义裂隙张量及脆性指数的概念，并介绍了它们的作用和意义。     

含微裂纹材料的损伤理论

本文从含微裂纹材料的变形能出发引出了裂纹的方位张量。在考虑裂纹受压闭合与滑动摩擦的基础上,给出了损伤张量、损伤应变及有效弹性常数。文中给出了损伤机构离散化的方法,并对方位密度给出了演化方程。最后给出一个单向拉压的应力应变关系例子,并揭示了裂纹扩展时的应力突跌现象。     

Discontinuous bifurcation analysis of thermodynamically consistent gradient poroplastic materials

In this work, analytical and numerical solutions of the condition for discontinuous bifurcation of thermodynamically consistent gradient-based poroplastic materials are obtained and evaluated. The main aim is the analysis of the potentials for localized failure modes in the form of discontinuous bifurcation in partially saturated gradient-based poroplastic materials as well as the dependence of these potentials on the current hydraulic and stress conditions. Also the main differences with the localization conditions of the related local theory for poroplastic materials are evaluated to perfectly understand the regularization capabilities of the non-local gradient-based one. Firstly, the condition for discontinuous bifurcation is formulated from wave propagation analyses in poroplastic media. The material formulation employed in this work for the spectral properties evaluation of the discontinuous bifurcation condition is the thermodynamically consistent, gradient-based modified Cam Clay model for partially saturated porous media previously proposed by the authors. The main and novel feature of this constitutive theory is the inclusion of a gradient internal length of the porous phase which, together with the characteristic length of the solid skeleton, comprehensively defined the non-local characteristics of the represented porous material. After presenting the fundamental equations of the thermodynamically consistent gradient based poroplastic constitutive model, the analytical expressions of the critical hardening/softening modulus for discontinuous bifurcation under both drained and undrained conditions are obtained. As a particular case, the related local constitutive model is also evaluated from the discontinuous bifurcation condition stand point. Then, the localization analysis of the thermodynamically consistent non-local and local poroplastic Cam Clay theories is performed. The results demonstrate, on the one hand and related to the local poroplastic material, the decisive role of the pore pressure and of the volumetric non-associativity degree on the location of the transition point between ductile and brittle failure regimes in the stress space. On the other hand, the results demonstrate as well the regularization capabilities of the non-local gradient-based poroplastic theory, with exception of a particular stress condition which is also evaluated in this work. Finally, it is also shown that, due to dependence of the characteristic lengths for the pore and skeleton phases on the hydraulic and stress conditions, the non-local theory is able to reproduce the strong reduction of failure diffusion that takes place under both, low confinement and low pore pressure of partially saturated porous materials, without loosing, however, the ellipticity of the related differential equations.     

DEM analysis of the effect of joint geometry on the shear behavior of rocks

In order to comprehensively investigate the effect of different joint geometries on the shear behavior of rocks, the Distinct Element Method (DEM) was utilized with a new bond contact model. A series of direct shear tests on coplanar and non-coplanar jointed rocks was simulated using the PFC2D software, which incorporates our bond contact model. Both coplanar jointed rocks with different joint persistence and non-coplanar ones with different joint inclinations were simulated and investigated numerically. The numerical results were compared and discussed with relevant laboratory tests as well as some reported numerical works. The results show that for coplanar jointed rocks, the peak shear stress decreases nonlinearly with the joint persistence, and the failure process can be divided into four stages: elastic shearing phase, crack propagation, failure of rock bridges, and residual phase. For non-coplanar jointed rocks, as the absolute value of the inclination angle of the rock joints increases, its shear strength increases, changing the failure patterns and the length of new fractures between existing cracks. When the absolute value increases from 15° to 30°, the average shear capacity increases the most as 39%, while the shear capacity increases the least as 2.9% when the absolute value changes from 45° to 60°. There is a good consistency of the failure patterns obtained from experiments and numerical tests. All these demonstrate that the DEM can be further applied to rock mechanics and practical rock engineering with confidence in the future.