动载下3条断续裂隙岩样的裂缝贯通机制

采用含3条断续预制裂隙的类砂岩模型试样进行单轴动力加载试验,对不同裂隙空间位置条件下断续裂隙岩体中裂隙的贯通机制进行了研究。静、动荷载下的对比研究成果显示：不同空间位置的裂隙贯通方式存在较大差异,且对动载的响应不同;动载下分支裂纹扩展及贯通具有惯性效应,即动载下裂尖次生共面、次生倾斜裂纹起裂后易朝原起裂方向快速发展,动载下易在两预制裂隙内端部产生直接贯通。这与静载下岩桥处的贯通常通过分支裂纹拐折扩展、相连不同,这也是导致裂隙试样中低应变速率下强度增大（即速率效应）的主要原因。同时,试验结果也表明：含裂隙试样静、动荷载下裂隙间的多次贯通是导致其呈现出渐进破坏特征的主要原因。     

含裂隙岩体单轴压缩裂纹扩展机制离散元分析EI北大核心CSCD

将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。     

含裂隙岩体单轴压缩裂纹扩展机制离散元分析

将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。     

静动载下共面非贯通裂隙贯通机制分析

共面非贯通裂隙的贯通机制对于确定外载下岩质边坡的滑动面位置及滑动面综合抗剪强度至关重要。采用预制共面非贯通裂隙石膏模型试样单轴静动载对比试验,对不同裂隙倾角共面排列裂隙的扩展、贯通过程进行了观测,研究显示：共面非贯通裂隙不同倾角下的贯通模式存在较大差异,裂隙倾角为0o～35o时以裂隙面发生闭合变形为主;45o～65o时裂隙间较易出现剪切型破坏模式;75o～90o预制裂隙面较难产生滑动,裂隙试样主要产生劈裂形式的破坏。动载下预制裂隙试样裂尖翼裂纹及次生共面裂纹起裂后易朝原起裂方向快速发展;易在两预制裂隙内端部产生直接贯通,这与静载下岩桥处的贯通常通过分支裂纹拐折扩展、相连不同。含共面非贯通裂隙试样在裂隙倾角为35o左右时强度呈现最小值,这与贯通性裂隙试样裂隙倾角为60o左右呈现最小值相差较大,这是因为裂隙面摩擦强度没来得及发挥作用所致。因此,含非贯通节理裂隙岩体的综合抗剪强度公式应引入强度发挥系数,以充分考虑岩桥胶结强度与裂隙面摩擦强度不能同步发挥作用的破坏本质。     

基于CT-CA的天然裂隙岩石受拉破坏特性研究

岩石内部存在的裂隙、孔洞等天然损伤对岩石的力学性能和破坏过程有重要影响,依据细胞自动机理论结合CT无损识别技术实现了含天然裂隙岩石在劈裂条件下裂纹扩展和贯通全过程及其力学性能变化规律的研究。从裂隙砂岩的真实细观结构出发,构建了天然裂隙岩石的数值计算模型,运用CASRock数值计算软件完成了含不同裂隙倾角的砂岩劈裂破坏的数值试验,分析了裂隙倾角对砂岩的力学特性、裂纹扩展过程及能量演化的影响规律。研究表明:(1)天然裂隙砂岩的抗拉强度与裂隙倾角密切相关,随着裂隙倾角的增加,其抗拉强度呈现先减小后增加的趋势;(2)裂隙起裂于天然裂隙尖端,当裂隙倾角0°≤θ<48°时,岩样的破坏是由错开型裂纹引起,裂纹沿着与天然裂隙近垂直方向扩展;当裂隙倾角48°≤θ<94°时,岩样的破坏是由张开型裂纹引起,裂纹沿着与天然裂隙近平行方向扩展;(3)劈裂过程中裂纹尖端应力场存在拉应力区和压应力区,拉应力造成翼裂纹由天然裂隙尖端沿加载端方向萌生扩展,而压应力则引发次生裂纹沿天然裂隙方向扩展;(4)含天然裂隙砂岩劈裂破坏过程能量演化可划分为4个阶段,随裂隙倾角的增大,峰值点处的总能量密度、弹性能密度先缓慢减少再迅速增加,但对岩样耗散能影响不大。     

裂隙充填影响巴西圆盘抗拉力学特性试验研究

裂隙是影响岩石力学特性和破裂模式最重要的因素之一,通过试验、数字照相技术和声发射监测对完整和含不同裂隙倾角 单裂隙、不同岩桥与裂隙夹角 双裂隙充填与非充填类岩石材料巴西圆盘的抗拉强度和破裂模式进行研究,探讨了随着 、 的不断变化和裂隙充填与否对试样最终破坏机制的影响。试验结果表明：（1）随着裂隙倾角 不断增大,单裂隙试样的抗拉强度先减小后增大,然后又减小;预制裂隙尖端萌生的翼裂纹贯通造成试样破坏;（2）双裂隙试样的抗拉强度随着 的增大先减小后增大;预制裂隙尖端萌生的翼裂纹和远场裂纹的扩展导致试样破坏;（3）在 、 相同的情况下,充填试样抗拉强度明显要高于非充填试样;预制裂隙充填与否对试样裂纹扩展的时间和裂纹数目影响较大;（4）加载初期,声发射较为平稳;宏观裂纹出现或者抗拉强度跌落时声发射将会变的异常活跃;在 、 相同的情况下充填试样声发射的起伏变化更为剧烈。     

高渗压条件下压剪岩石裂纹断裂损伤演化机制研究

探讨了高渗透压作用下压剪岩石裂纹的起裂规律及分支裂纹尖端应力强度因子的演变规律,建立了高渗压下裂隙岩体发生拉剪破坏的临界水压力值和初裂强度判据,分析了不同渗压作用下裂纹的扩展情况表明,渗透压的存在加剧了分支裂纹的扩展,高渗透压作用下分支裂纹扩展由稳定扩展变成不稳定扩展,并导致分支裂纹尖端岩桥剪切破坏,同时考虑分支裂纹的相互作用,建立了高渗透压作用下压剪岩石裂纹体岩桥剪切贯通的断裂破坏力学模型,最后依据裂隙岩体的损伤力学效应研究了岩体的初始损伤及损伤演化柔度张量,提出了高渗压下压剪岩石裂纹渐进破坏的损伤演化方程。该理论为定量研究高渗压下裂隙岩体的失稳破坏提供了理论依据。     

基于DIC技术的3D打印节理试件破裂机制研究

为了准确表征不同角度预制节理岩石在单轴压缩下的变形破坏模式,基于3D打印技术制作了节理模型用于模拟岩体中的结构面,通过水泥砂浆的浇筑得到含不同角度预制节理的岩石试件并进行单轴压缩试验,同时采用数字图像相关技术（DIC）观测、分析试验过程中试件裂纹产生、扩展以及贯通过程。结果表明：随着预制节理从0°增加到90°,试件强度与峰值应变均呈现先降低后升高的变化趋势,0°和45°试件弹性模量相对于完整试件有所降低。基于DIC检测结果,0°、30°、45°及60°试件裂纹皆从预制节理尖端部位起裂,各角度试件的起裂应力与试件强度变化规律一致。各角度试样起裂时在剪应力控制下以剪切翼型裂纹形式起裂,0°与45°试件裂纹在扩展过程由剪切发展为张拉型裂纹,30°和60°试件以剪切裂纹形式贯穿始终,90°试件从底部起裂并最终表现为张拉破坏。研究还发现,下翼起裂角θ2和上翼起裂角θ1之间存在明显的线性正相关关系,关系式为θ2 =0.828 6θ1 +12.185,且起裂应力大小变化与峰值应力变化一致,皆随节理角度的增加先减小后增大。     

裂隙岩石单轴压缩损伤扩展细观机理CT分析初探

完成了单一裂纹的裂隙砂岩单轴压缩条件下细观损伤破坏机理CT实时试验,得到了裂纹萌生、发展、宏观裂纹形成、破坏等各阶段的CT图像、CT数和CT数方差。结果表明,与无预制裂隙的岩石试样相比,已有预制裂纹对新的裂纹的起裂位置及贯通性宏观破坏裂纹的形成具有重要影响,预制裂纹的存在导致裂隙砂岩试样的扩容量大于完整试样破坏时的扩容量。     

单轴压缩条件下含三叉裂隙类岩石试样力学特性的细观研究

三叉裂隙是自然界普遍存在的一种岩体缺陷形式,其对岩体的力学特性有重要影响。对含预制三叉裂隙的水泥砂浆试样进行室内单轴压缩试验,配合使用摄像机拍摄裂纹的起裂、扩展、贯通过程,通过数字图像技术处理获取试样的应变场云图,并结合PFC2D程序研究不同?、? 条件下试样的强度特征、裂纹模式和裂纹演化扩展规律。研究表明：三叉裂隙对试样单轴抗压强度有明显的削弱作用。当? 恒定为120°时,试样在? = 30°时达到最大抗压强度;当? 恒定为90°时,随?增大,试样抗压强度呈先减小后增大的趋势,且当? = 45°时达到最大抗压强度。试样产生的裂纹可分为3类,分别是张拉型裂纹(Ⅰ型裂纹)、剪切型裂纹(Ⅱ型裂纹)、混合型裂纹(Ⅲ型裂纹)。这3类裂纹通常从裂隙尖端开始产生,并且Ⅰ型裂纹沿加载方向扩展,通常未扩展至试样边界;Ⅱ型和Ⅲ型裂纹通常与加载方向呈一定角度扩展至试样边界。通过对裂纹的几何形态和组成宏观裂纹的微裂纹成分的分析,得知导致含三叉裂隙试样在单轴压缩条件下失效的是张拉破坏。数字图像技术得到的应变云图表明,当载荷达到一定阶段,裂隙尖端出现应力集中,微破裂开始发育并聚集成微破裂区,微破裂区扩大产生宏观裂纹。通过对主应变和剪应变云图分析,发现导致试样失效的是张拉破坏,剪应变在裂纹扩展过程中的影响较小。     

Numerical Study on Coalescence of Pre-Existing Flaw Pairs in Rock-Like Material

The present numerical study, which is an extension of our previous numerical analysis on cracking processes of a single pre-existing flaw, focuses on the coalescence of two pre-existing parallel open flaws in rock subjected to a uniaxial compressive loading. To facilitate a systematic investigation, the arrangements of the flaw pair are classified into 11 categories. Simulations engaging AUTODYN are conducted on each category. The numerical results are compared with some published physical experimental test results. Eleven typical coalescence patterns are obtained, which are in good agreement with the experimental results, which include two coalescence patterns obtained in flaw pair arrangements (II) and (VIII″) not being reported in previous studies. The information gathered in the simulations helps identify the type (tensile/shear) of each crack segment involved in the coalescence. Most of the coalescence cracks initiate at or around the flaw tips, except those in flaw pair arrangements (II) and (IX′) with a very short ligament length, in which the coalescence cracks initiate on the flaw surfaces away from the flaw tip regions. Based on the numerical simulation results, the properties of the 11 coalescence patterns are obtained. Except those in flaw pair arrangements (II) and (IX′), the other coalescence patterns can be interpreted with respect to the basic crack types—tensile wing crack, horsetail crack and anti-wing crack. In addition, based on the type of crack segments involved in coalescence, namely tensile and shear, the coalescence can be classified into T mode (tensile mode), S mode (shear mode) and TS mode (mixed tensile–shear mode).     

Cracking Processes in Rock-Like Material Containing a Single Flaw Under Uniaxial Compression: A Numerical Study Based on Parallel Bonded-Particle Model Approach

Cracking processes have been extensively studied in brittle rock and rock-like materials. Due to the experimental limitations and the complexity of rock texture, details of the cracking processes could not always be observed and assessed comprehensively. To contribute to this field of research, a numerical approach based on the particle element model was used in present study. It would give us insights into what is happening to crack initiation, propagation and coalescence. Parallel bond model, a type of bonded-particle model, was used to numerically simulate the cracking process in rock-like material containing a single flaw under uniaxial vertical compression. The single flaw’s inclinations varied from 0° to 75° measured from the horizontal. As the uniaxial compression load was increased, multiple new microcracks initiated in the rock, which later propagated and eventually coalesced into longer macrocracks. The inclination of the pre-existing flaw was found to have a strong influence on the crack initiation and propagation patterns. The simulations replicated most of the phenomena observed in the physical experiments, such as the type, the initiation location and the initiate angle of the first cracks, as well as the development of hair-line cracks, which later evolved to macrocracks. Analyses of the parallel bond forces and displacement fields revealed some important mechanisms of the cracking processes. The first cracks typically initiated from the tensile stress concentration regions, in which the tensile stress was partially released after their initiation. The tensile stress concentration regions subsequently shifted outwards close to the propagating tips of the first cracks. The initiation and propagation of the first cracks would not significantly influence the compressive stress singularity at the flaw tips, which was the driving force of the initiation of secondary cracks. The initiation of microcracking zone consisting almost exclusively of micro-tensile cracks, and that of microcracking zone consisting of micro-tensile cracks and mixed micro-tensile and shear cracks, were found to be correlated with two distinct types of displacement fields, namely type I (DF_I) and type II (DF_II), respectively.     

Coupled flow network and discrete element modeling of injection-induced crack propagation and coalescence in brittle rock

We present a numerical analysis on injection-induced crack propagation and coalescence in brittle rock. The DEM network coupling model in PFC is modified to capture the evolution of fracture geometry. An improved fluid flow model for fractured porous media is proposed and coupled with a bond-based DEM model to simulate the interactions among cracks induced by injecting fluid in two nearby flaws at identical injection rates. The material parameters are calibrated based on the macro-properties of Lac du Bonnet granite and KGD solution. A grain-based model, which generates larger grains from assembles of particles bonded together, is calibrated to identify the microscopic mechanical and hydraulic parameters of Lac du Bonnet granite such that the DEM model yields a ratio between the compressive and tensile strength consistent with experiments. The simulations of fluid injection reveal that the initial flaw direction plays a crucial role in crack interaction and coalescence pattern. When two initial flaws are aligned, cracks generally propagate faster. Some geometrical measures from graph theory are used to analyze the geometry and connectivity of the crack network. The results reveal that initial flaws in the same direction may lead to a well-connected crack network with higher global efficiency.

     

断续三裂隙砂岩强度破坏和裂纹扩展特征研究

利用岩石力学伺服试验机,对尺寸为80 mm×160 mm×30 mm的断续三裂隙砂岩试样进行了单轴压缩试验,研究了岩桥倾角? 2对断续三裂隙砂岩（? = 30°和? 1 = 60°）强度破坏和裂纹扩展特征的影响规律。与完整砂岩试样相比,断续三裂隙砂岩试样应力-应变曲线显示了较多的应力跌落,其峰值强度也呈显著降低趋势,但降低程度与? 2密切相关,随着? 2从75°增加到90°,峰值强度从82.04 MPa 降低到77.82 MPa,而当? 2从90°增加到120°,其峰值强度无明显变化。完整砂岩试样呈现轴向劈裂脆性破坏,而断续三裂隙砂岩试样是由许多从裂隙尖端产生的裂纹扩展与汇合,导致了其失稳破坏。通过照相量测技术,探讨了? 2对断续三裂隙砂岩试样裂纹扩展特征的影响：? 2为75°、90°和105°的断续三裂隙试样中裂隙①、③和②、③之间均出现了贯通,而裂隙①、②之间无任何贯通;? 2为120°的断续三裂隙试样中裂隙①、③和①、②之间均出现了贯通,而②、③之间无任何贯通。最后给出了断续三裂隙砂岩试样宏观变形特性与裂纹扩展过程之间的关系。     

A three-dimensional numerical investigation of the fracture of rock specimens containing a pre-existing surface flaw

Three-dimensional surface crack initiation and propagation in two kinds of heterogeneous rocks were numerically investigated via parallel finite element analysis using a supercomputer. Numerically simulated rock specimens containing a pre-existing flaw were subjected to uniaxial compression until failure. The initiation and propagation of wing cracks, anti-wing cracks, and shell-like cracks were reproduced by numerical simulations. The numerically simulated results demonstrate that the further propagation of wing cracks and shell-like cracks stop due to their wrapping (curving) behavior in three-dimensional spaces, even if the applied loads continue to increase. Furthermore, rock heterogeneity could significantly influence crack propagation patterns and the peak uniaxial compressive strengths of rock specimens. Moreover, anti-wing cracks only appeared in relatively heterogeneous rocks, and the peak uniaxial compressive strengths of the specimens were observed to depend on the inclination of the pre-existing flaw. Finally, the mechanism of surface crack propagation is discussed in the context of numerically simulated anti-plane loading tests, wherein it was identified that Mode III loading (anti-plane loading) does not lead to Mode III fracture in rocks due to their high ratio of uniaxial compressive strength to tensile strength. This finding could explain the lateral growth of an existing flaw in its own plane, which is a phenomenon that has not been observed in laboratory experiments.     

Crack Initiation, Propagation and Coalescence in Rock-Like Material Containing Two Flaws: a Numerical Study Based on Bonded-Particle Model Approach

Cracking and coalescence behavior in a rectangular rock-like specimen containing two parallel (stepped and coplanar) pre-existing open flaws under uniaxial compression load has been numerically studied by a parallel bonded-particle model, which is a type of bonded-particle model. Crack initiation and propagation from two flaws replicate most of the phenomena observed in prior physical experiments, such as the type (tensile/shear) and the initiation stress of the first crack, as well as the coalescence pattern. Eight crack coalescence categories representing different crack types and trajectories are identified. New coalescence categories namely “New 1” and “New 2”, which are first observed in the present simulation, are incorporated into categories 3 and 4, and category 5 previously proposed by the MIT Rock Mechanics Research Group, respectively. The flaw inclination angle (β), the ligament length (L) (spacing between two flaws) and the bridging angle (α) (inclination of a line linking up the inner flaw tips, between two flaws) have different effects on the coalescence patterns, coalescence stresses (before, at or post the peak stress) as well as peak strength of specimens. Some insights on the coalescence processes, such as the initiation of cracks in the intact part of specimens at a distance away from the flaw tips, and coalescence due to the development and linkage of a number of steeply inclined to vertical macro-tensile cracks are revealed by the present numerical study.     

模拟岩石中裂纹扩展连接的近场动力学方法

脆性岩石材料在压应力作用下常出现两类裂纹：翼型张拉裂纹和次生剪切裂纹。近场动力学是一种新型的无网格数值计算方法。在近场动力学理论中,采用积分形式的控制方程代替微分形式的控制方程使得该数值算法在断裂问题上具有独特的优势。将Mohr-Coulomb准则和最大主应力准则引入非普通“态”基近场动力学理论中,分别用于模拟材料常见的压剪和张拉破坏。这种扩展的非普通“态”基近场动力学可以有效地模拟脆性岩石材料在多种受力状态下的裂纹起裂、扩展和连接问题。通过5个不同的数值算例说明该数值算法在处理脆性岩石材料断裂问题的有效性和准确性。首先,通过模拟含圆孔的弹性板拉伸数值试验说明该数值算法的有效性和准确性。其次,数值模拟了简单三点弯曲试验以及不使用其他外部准则条件下动荷载作用下裂纹的分叉试验,所得结果与其他试验结果或数值结果相吻合,从而验证了该理论的有效性。然后,模拟了包含斜裂纹的巴西圆盘试验,裂纹扩展路径和计算所得的断裂韧度同样吻合于试验结果。最后,模拟了单轴压缩状态下,预制裂纹试样的裂纹扩展和连接问题。将该数值算法与试验结果对比表明,所提出的数值方法可以模拟和预测岩石类材料的张拉和压剪裂纹的起裂、扩展和连接行为。     

含预制裂纹的脆性岩石单轴压缩下渐进性破坏过程的试验研究

对含预制裂纹的花岗岩进行单轴压缩试验研究预制裂纹倾角α对脆性岩石渐进性破坏过程的影响。首先埘破坏过程的轴向应力-横向应变曲线进行总结和讨论,然后分析预制裂纹与加载方向夹角α埘岩石的应应力门槛值：裂纹起始应力σci、裂纹扩展应力σed、峰值强度矿σf,由应变片记录的应力-应变曲线和试样的表嘶裂纹扩展情况的影响机制。结果表明,含有预制裂纹的岩石试样进行加载试验过程中,顶制裂纹倾角α的变化成了决定脆性岩石破裂办式的主要因素。故在对含节理、裂隙的脆性岩石的工程应用上,通过对岩休的轴向应力-横向应变向应变典线进行分析,可以对地下开挖工程起到指导设计开挖方式及支护形式的作用。