岩石损伤过程中的热-流-力耦合模型及其应用初探

岩石损伤过程热-流-力（THM）耦合问题的研究对于深部采矿等许多工程领域都具有重要的理论意义。以岩石的损伤为主线,在多场耦合分析方程中引入损伤变量,基于质量守恒和能量守恒原理,提出岩体损伤过程中的THM耦合模型。通过把均匀弹性介质THM耦合响应的模拟结果与理论分析结果进行对比,验证了程序及有限元实施的正确性。然后,用该耦合模型进行了不同地应力条件下流固耦合过程的数值模拟,探讨了水压力对于岩石损伤过程的作用机制。数值模拟表明,水压力导致了拉伸损伤范围的扩大和损伤程度的加剧,同时亦对剪切损伤具有抑制作用。     

基于数字图像技术的岩土材料有限元-离散元分析

将数字图像技术与Munjiza提出的有限元法-离散元法（FDEM）耦合分析方法结合,研制了可表征岩石真实非均质性的FDEM分析系统,为从细观角度对岩体进行建模以及研究岩体的破裂机制提供了新途径。该系统借助数字图像技术从岩石断面的图像获取岩石材料的真实细观结构,借助成熟的网格剖分技术,将其映射到FDEM计算网格中,从而克服了原有FDEM在考虑材料非均质性所存在的不足。利用该系统,进行了巴西圆盘劈裂试验的数值模拟,再现了花岗岩在荷载作用下的真实破裂过程。数值模拟结果表明,考虑非均质性的岩样应力分布呈现出非对称性,岩石的细观结构对岩石中裂纹的扩展及应力分布有重要影响。     

储层非均质性对水力压裂的影响

从岩石细观非均质性的特点出发,采用RFPA2D-Flow软件对单孔和双孔数值模型进行压裂计算,研究岩石非均质性对水力压裂的响应,重点探讨双孔模型孔间吸引效应对裂纹演化形态的影响。岩石细观单元的力学、水力学特性由统计分布生成以体现岩石的随机不均质性,水力压裂过程中流体压力传递通过单元渗流-损伤耦合迭代来实现。数值计算结果表明:(1)岩石非均质性影响裂缝的扩展形态,导致水力裂纹尖端微裂纹的分支。随着均质度的增加,水力裂纹的扩展形态变得更加平直光滑,单孔模型两侧裂纹更加对称,双孔间裂纹的连通性变差。(2)岩石的非均匀性对于岩石的起裂压力和地层破裂压力影响较大。随着均质度的增大,起裂应力和地层破裂应力增大,并且两者间的差值逐渐变小,在储层为均质的条件下,两者几乎相等。(3)相同的边界条件下,均质模型的应力分布曲线光滑连续,非均质模型的应力分布曲线呈现出明显波动,井眼对称剖面上的应力分布不尽相同,反映了细观单元强度非均匀性及裂缝扩展形态对应力分布的影响。(4)双孔模型孔间存在孔隙水压力增加带,孔间产生吸引效应,双孔方位影响临界压力。研究结果对水压裂试验设计和现场压裂施工具有一定的参考意义。     

含水状态下膏溶角砾岩破裂全程的细观力学试验研究

利用新型岩石细观力学试验系统,对石家庄-太原铁路专线太行山隧道工程6#斜井的膏溶角砾岩在单轴压缩荷载作用下的岩石细观损伤裂纹扩展直至破裂的全程进行了实时观测,得到了不同含水状态下岩石初始损伤微裂纹的萌生、扩展、连接、贯通直至宏观破裂的数字显微和全场实时图像。从岩石细观损伤机制出发,结合同时获得的应力-应变曲线,将损伤破裂过程分为4个阶段,得到了含水状态下岩石损伤发展演化的初步规律,并给出了应力损伤门槛值,发现水环境影响下造成的初始损伤微裂纹是水对膨胀软岩力学性质产生软化作用的重要因素。     

单轴压缩下岩石蠕变失稳破坏过程数值模拟

在岩石破裂过程分析（RFPA2D）系统的基础上,考虑岩石损伤过程的时间因素影响,引入岩石细观单元蠕变本构方程,建立了考虑流变效应的岩石破裂过程RFPA2D数值模型。应用该模型模拟了恒定荷载作用下岩石的蠕变破坏过程,得到了岩石蠕变破裂的3个典型阶段：初始蠕变阶段、稳定蠕变阶段和加速蠕变阶段,模拟结果同实验室试验所观察到的现象十分吻合。这表明考虑流变效应的RFPA2D数值模型适用于模拟岩石的蠕变破坏这一复杂的、非线性演化问题。此外,数值模拟还揭示了岩石的宏观蠕变破坏实质上是细观层次上单元损伤累计的结果,这些结论对岩石工程的长期稳定性研究具有重要的理论指导和实践意义。     

基于摩尔-库仑准则的断续节理岩体复合损伤本构模型

针对工程实际中断续节理裂隙岩体的损伤本构模型,假设岩石微元强度服从Weibull随机分布,以摩尔-库仑破坏准则作为描述微元强度的表示方法,推导出细观损伤变量。利用能量和断裂力学理论,综合考虑节理几何特征及力学特性,推导宏观损伤变量计算公式。基于Lemaitre应变等效假设,考虑宏细观缺陷耦合作用,推导出复合损伤变量,建立了基于摩尔-库仑准则的宏细观缺陷耦合作用的断续裂隙岩体损伤本构模型。研究结果表明：（1）采用摩尔-库仑准则作为描述微元强度的统计分布变量建立的损伤模型能够较好地反映岩石内部缺陷分布和变形特征,该模型真实地反映岩石微元强度受应力状态的影响。（2）该模型建立的理论曲线与断续节理岩体的试验曲线吻合较好。（3）节理裂隙岩体宏观损伤变量及峰值强度随节理倾角的变化规律与综合考虑宏细观耦合作用下的损伤变量及裂隙岩体峰值强度随节理倾角的变化规律基本一致。（4）宏细观耦合作用下的等效弹性模量与节理贯通率呈非线性负相关;在节理倾角一定的情况下,损伤变量与节理长度呈非线性正相关;在贯通率较小时,岩体的宏观损伤变量与内摩擦角的关系呈线性负相关变化,贯通率达到一定程度,线性关系变成非线性关系。     

基于大直径霍普金森压杆数值试验的非均匀介质动态破坏过程分析

基于细观损伤力学基础而开发的动态岩石破裂过程分析系统RFPA2D代替大直径SHPB试验技术对非均匀介质的动态破坏过程和动态性能进行数值模拟,分析了加载波形和介质的非均匀性对数值试样的动态性能,如应力-时间曲线、应变-时间曲线、应力-应变曲线和应变率-时间曲线的影响。分析表明,大直径SHPB弥散效应对试验结果的影响较大,选择合适的加载波形可以减小SHPB装置中应力波的弥散效应,得到较准确的试验结果,其中三角形波加载可以有效降低大直径SHPB动态测试中的应力波弥散,是岩石等非均匀材料SHPB动态测试的较理想加载波形;在相同加载条件下,岩石的非均匀性对波的弥散效应影响不大,但非均匀岩石的试验曲线比均匀岩石的曲线在波峰后都出现较大的振荡,这主要是由于不同均质度的岩石其单元具有不同的破坏分布造成的,不是波形弥散造成的。     

基于断裂及高温损伤的岩石蠕变模型研究

为了反映热-力耦合作用下岩石蠕变变形的全过程,依据断裂力学原理提出了岩石裂纹扩展的临界损伤应力和一个新的可描述岩石在稳态蠕变阶段与临界损伤应力相关的非线性黏性分量,在传统西原模型和Burgers模型的基础上,将指数形式的损伤变量、临界损伤应力以及与其有关的非线性黏性分量引入到流变微分方程,通过叠加原理推导了考虑温度效应的单轴和三轴压缩条件下岩石的流变本构关系,建立了岩石的热-力耦合损伤蠕变本构模型。利用不同温度、不同应力条件下花岗岩的三轴蠕变试验曲线和本文蠕变模型的计算曲线进行比较,结果表明本文蠕变模型能够较好地模拟岩石在初始瞬态、稳态和加速蠕变阶段全过程的变形规律,验证了所建模型的有效性和合理性。该模型为分析高温、高应力环境下岩石工程的长期变形和稳定情况提供了理论依据。     

基于Mohr-Coulomb准则的断续节理岩体复合损伤本构模型

针对工程实际中断续节理裂隙岩体的损伤本构模型,假设岩石微元强度服从Weibull随机分布,以摩尔–库仑破坏准则作为描述微元强度的表示方法,推导出细观损伤变量。利用能量和断裂力学理论,综合考虑节理几何特征及力学特性,推导宏观损伤变量计算公式。基于Lemaitre应变等效假设,考虑宏细观缺陷耦合作用,推导出复合损伤变量,建立了基于摩尔–库仑准则的宏细观缺陷耦合作用的断续裂隙岩体损伤本构模型。研究结果表明:(1)采用摩尔–库仑准则作为描述微元强度的统计分布变量建立的损伤模型能够较好地反映岩石内部缺陷分布和变形特征,该模型真实地反映岩石微元强度受应力状态的影响。(2)该模型建立的理论曲线与断续节理岩体的试验曲线吻合较好。(3)节理裂隙岩体宏观损伤变量及峰值强度随节理倾角的变化规律与综合考虑宏细观耦合作用下的损伤变量及裂隙岩体峰值强度随节理倾角的变化规律基本一致。(4)宏细观耦合作用下的等效弹性模量与节理贯通率呈非线性负相关;在节理倾角一定的情况下,损伤变量随节理长度呈非线性正相关;在贯通率较小时,岩体的宏观损伤变量与内摩擦角的关系呈线性负相关变化,贯通率达到一定程度,线性关系变成非线性关系。     

盐岩储气库温度-渗流-应力-损伤耦合模型研究

推导了含夹层盐岩储气库在地应力、注-采气压力和温度联合作用下的温度-渗流-应力（THM）-损伤耦合数学模型,该模型不仅包含了盐岩的热损伤效应,还考虑了气体渗透的Klinkenberg效应以及气体沿夹层层面流动对渗透规律的影响。在此基础上,结合金坛储气库建设,采用数值仿真技术,对储气库在注-采气过程中考虑THM耦合与不耦合情况下的气体渗透、温度分布以及应力分布规律进行了详细的对比分析和探讨,研究结果表明：含夹层盐岩储气库应力-渗透-温度耦合效应十分明显。该问题的研究,为全面深入分析盐岩储气库的安全性和实用性提供了思路及前提,具有重要的理论及实践意义。     

Micromechanical Model for Simulating the Fracture Process of Rock

Summary A micromechanical model is proposed to study the deformation and failure process of rock based on knowledge of heterogeneity of rock at the mesoscopic level. In this numerical model, the heterogeneity of rock at the mesoscopic level is considered by assuming the material properties in rock conform to the Weibull distribution. Elastic damage mechanics is used to describe the constitutive law of meso-level elements, the finite element method is employed as the basic stress analysis tool and the maximum tensile strain criterion as well as the Mohr-Coulomb criterion is utilized as the damage threshold. A simple method, similar to a smeared crack model, is used for tracing the crack propagation process and interaction of multiple cracks. Based on this model, a numerical simulation program named Rock Failure Process Analysis Code (RFPA) is developed. The influence of parameters that include the Weibull distribution parameters, constitutive parameters of meso-level elements and number of elements in the numerical model, are discussed in detail. It is shown that the homogeneity index is the most important factor to simulate material failure with this model. This model is able to capture the complete mechanical responses of rock, which includes the crack patterns associated with different loading stages and loading conditions, localization of deformation, stress redistribution and failure process. The numerical simulation of rock specimens under a variety of static loading conditions is presented, and the results compare well with experimental results.     

单轴压缩下横观各向同性岩石破裂过程的数值模拟

采用基于细观损伤力学基础上开发的RFPA2D数值模拟软件,用2种不同的岩石材料来组成7个不同岩层倾角的横观各向同性的岩石试件,通过单轴加载数值模拟试验,模拟横观各向同性岩石渐进破裂的整个过程,分析了岩层与最大主应力之间的倾角和强度之间的关系,讨论了不同岩层倾角的横观各向同性岩体的不同破裂模式及其破坏准则。     

Numerical analysis of the failure process around a circular opening in rock

Damage and fracture propagation around underground excavations are important issues in rock engineering. The analysis of quasi-brittle materials can be performed using constitutive laws based upon damage mechanics. The finite element code RFPA2D (Rock Failure Process Analysis) based on damage mechanics was used to simulate a loading-type failure process around an underground excavation (model tunnel) in brittle rock. One of the features of RFPA2D is the capability of modeling heterogeneous materials. In the current model, the effect of the homogeneous index (m) of rock on the failure modes of a model tunnel in rock was studied. In addition, by recording the number of damaged elements and the associated amount of energy released, RFPA2D is able to simulate acoustic activities around circular openings in rock. The results of a numerical simulation of a model tunnel were in very good agreement with the experimental test using the acoustic emission technique. Finally, the influence of the lateral confining pressure on the failure mechanism of the rock around the model tunnel was also investigated by numerical simulations.     

含预制裂纹巴西盘试样破裂模式的数值模拟

利用岩石破裂过程分析系统（RFPA）软件, 进行了岩石含中心裂纹的巴西盘在单轴压缩荷载作用下破裂过程的数值模拟。数值模拟再现了含不同预制裂纹角度时巴西盘试样所表现出不同的破裂模式以及主裂纹的产生、扩展、贯通、次生裂纹的产生等过程,得到的破裂模式结果与有关文献中的实验结果十分吻合。除此以外,数值试验还能够给出在物理实验中不能观察到的应力场信息及岩石的破裂机制。     

动载荷作用下岩石破坏过程的数值试验研究

采用基于细观损伤力学基础上开发的动态版RFPA2D数值模拟软件,对动载荷作用下应力波延续时间、应力波峰值和围压对岩石试样破坏的影响进行了数值研究,结果表明,应力波延续时间较短,则尾随应力波波前的高应力区范围较窄,应力波衰减较快;相反,应力波延续时间较长,则紧跟应力波波前的高应力区范围较大,岩石处于破坏状态的时间延长,岩石的破碎程度加大。此外,存在一个合适的应力波延续时间,过分地加大应力波延续时间,反而不利于岩石裂隙的发育。动载荷的峰值越大,试样的破坏程度越大,当峰值达到一定值时,试样顶部呈现粉碎状,试样从上到下破坏程度逐渐减弱。在冲击载荷作用下的岩石随着围压的增加更难破碎,但当围压增大到一定程度时,岩石会突然失稳破坏。     

Numerical analysis of slope stability based on the gravity increase method

A micromechanical model is proposed for studying the stability and failure process of slopes based on the gravity increase method (GIM). In this numerical model the heterogeneity of rock at a mesoscopic level is considered by assuming that the material properties conform to the Weibull distribution. Elastic damage mechanics is a method used for describing the constitutive law of the meso-level element, the finite element method (FEM) is employed as the basic stress analysis tool, and the maximum tensile strain criterion and the Mohr–Coulomb criterion are utilised as the damage threshold. The numerical model is implemented into the Realistic Failure Process Analysis (RFPA) code using finite element programming, and an extended version of RFPA, i.e., RFPA-GIM, is developed to analyse the failure process and stability of slopes. In the numerical modelling with RFPA-GIM, the critical failure surface of slopes is obtained by increasing the gravity gradually but keeping material properties constant. The acoustic emission (AE) event rate is employed as the criterion for slope failure. The salient feature of the RFPA-GIM in stability analysis of slopes is that the critical failure surface as well as the safety factor can be obtained without any presumption for the shape and location of the failure surface. Several numerical tests have been conducted to demonstrate the feasibility of RFPA-GIM. Numerical results agree well with experimental results and those predicted using the FEM strength reduction method and conventional limit equilibrium analysis. Furthermore it is shown that selection of the AE rate as the criterion for slope failure is reasonable and effective. Finally, the RFPA-GIM is applied to several more complex cases, including slopes in jointed rock masses and layered rock formations. The results indicate that the RFPA-GIM is capable of capturing the mechanism of slope failure and has the potential for application in a larger range of geo-engineering.     

Slope stability analysis by SRM-based rock failure process analysis (RFPA)

The fundamental principles of the strength reduction method (SRM) are incorporated into the rock failure process analysis (RFPA) code to produce an RFPA–SRM method for analysing the failure process and stability of rock and soil slopes. The RFPA–SRM method not only satisfies the global equilibrium, strain-consistent, and non-linear constitutive relationship of rock and soil materials but also takes into account the heterogeneous characteristics of materials on the micro- and macro-scales. When the proposed method is used for slope stability analysis, both the critical failure surface and the safety factor can be obtained directly without any assumptions regarding the shape and location of the failure surface. The numerical results agree well with those obtained using conventional limit equilibrium and other FEM strength reduction methods. The proposed technique is applied to a number of more complex cases, including slopes in mixed rock–soil formations, rock layer formations, and highly jointed rock masses. It is shown that the RFPA–SRM method can describe the mechanism of failure of slopes and has potential applications in a large range of geoengineering problems.     

Numerical studies on the failure process and associated microseismicity in rock under triaxial compression

In this paper, firstly the mesoscopic elemental mechanical model for elastic damage is developed and implemented into the rock and tool interaction code (R-T^2D). Then the failure processes of a heterogeneous rock specimen subjected to a wide variety of confining pressures (0–80 MPa) are numerically investigated using the R-T^2D code. According to the simulated results, on the one hand, the numerical simulation reproduced some of the well-known phenomena observed by previous researchers in triaxial tests. Under uniaxial compression, rock failure is caused by a combination of axial splitting and shearing. Dilatancy and a post-failure stage with a descending load bearing capacity are the prominent characteristics of the failure. As the confining pressure increases, the extension of the failed sites is suppressed, but the individual failure sites become dense and link with each other to form a shear fracture plane. Correspondingly, the peak strength, the residual strength and the shear fracture plane angle increase, but the brittleness decreases. When the confining pressure is high enough, the specimen behaves in a plastic manner and a narrow shear fracture plane leads to its failure. The prominent characteristics are volume condensation, ductile cataclastic failure, and a constant load bearing capacity with increasing strain. On the other hand, the numerical simulation revealed some new phenomena. The highest microseismicity events occur in the post-failure stage instead of the maximal stress, and most of the microseismicity energies are released in the failure localization process. As the confining pressure increases, the microseismicity events in the non-linear deformation stage increase dramatically and the ratio between the energies dissipated at the non-linear deformation stage and those dissipated in the whole loading process increases correspondingly. Therefore, it is concluded that the developed mesoscopic elemental mechanical model for elastic damage is able to reproduce accurately the failure characteristics in loading rock specimens under triaxial conditions, and the numerical modelling can furthermore obtain some new clarifications of the rock fracture process.     

深部断续节理岩体中渗流对巷道稳定性影响的数值分析

应用自主开发的岩石破坏过程渗流与应力耦合分析软件F-RFPA2D,通过对高应力节理岩体中巷道的失稳破坏过程所进行的数值模拟分析,研究了深部断续节理岩体中地应力及渗流对巷道稳定性的影响。模拟结果不但直观形象地给出了巷道的渗流场、应力场、破坏区分布,而且得到了巷道周边关键部位的位移变化情况,并与无地下水的稳定性进行了比较。结果表明,岩体中断续节理与地应力及渗流相互作用,对巷道失稳的演化起了重要作用,如巷道与节理的相对方位决定了巷道坍塌的方向;岩桥对巷道的稳定起了重要的作用;虽然渗流应力加速了巷道的失稳破坏,但没有改变高地应力对巷道的主要破坏模式。该研究对于进一步理解深部岩体中节理、地应力及渗流对巷道稳定性的影响进而采取有效的加固措施具有重要的理论和实践意义