拉-压应力状态下脆性岩石强度及破坏机制颗粒流模拟

高地应力区地下岩体工程开挖常形成围岩拉-压应力状态,发生岩体张性破坏灾害。本文针对传统PFC平行黏结模型不能模拟脆性岩石高单轴压缩与拉伸强度比的问题,建立双抗拉强度参数的平行黏结强度准则,开展岩石拉-压数值模拟试验,得到了与物理试验接近的拉-压强度,实现了岩石高压拉强度比的模拟,并深入分析了破坏机制。研究结果表明随着围压的增加,破裂面倾角逐渐增大,由拉伸破裂转化为拉-剪破裂,发现了拉-压应力状态下破裂面处的雁行裂纹。根据细观颗粒位移场揭示了破裂面力学性质,随着围压的增加(破裂面倾角逐渐增大),破裂面张性逐渐减弱而剪性增强。可将拉-压应力状态下岩石损伤演化过程大致分为弹性变形阶段、稳定破裂发展阶段、不稳定破裂发展阶段和整体破裂阶段(峰后应力跌落及残余阶段)。围压较大时弹性变形和稳定破裂发展阶段相对较短,不稳定破裂发展阶段相对较长较剧烈,峰后残余阶段破裂面摩擦更强、应力波动较大。     

Cracking mechanism of coal under high-pressure water jet and its applications for enhanced coalbed methane drainage

This paper is concerned with the mechanism of coal breakage under high-pressure water jet (HPWJ) and its applications. A model of HPWJ impinging on coal target was established to study the cracking mechanism of coal under impact load. The characteristic and pressure distribution of HPWJ, the propagation characteristics of stress wave in coal, the mechanical properties of different coal particles, and the fracture characteristics of coal under HPWJ erosion were investigated theoretically and numerically. The results show that the shock wave and water wedge pressure are the main factors that cause coal breakage and crack propagation. The damage to the far-field coal particles affected by HPWJ is primarily caused by tensile stress, and the damage to the near-field coal particles affected by HPWJ is caused by the coupled effects of tensile stress and compressive stress. An erosion cavity is formed in the coal model with diameters of 1.25 to 2.5 times that of the jet at different depths. Meanwhile, the strong quasi-static pressure at the crack discontinuities further promotes the propagation of radial cracks around the erosion cavity to form a fracture zone, and the diameter of the fracture zone at different depths is 3.5 to 4.0 times that of the jet. In addition, the results of field application show that there is a significant difference between the methane parameters in the hydraulic flushing borehole and the conventional borehole; the average methane volume fraction and the average methane flow rate in hydraulic flushing boreholes are 3.85 and 3.67 times, respectively, that in conventional boreholes. Indicating hydraulic flushing can effectively promote the initiation and propagation of coal cracks. These results are of great significance to improve coalbed methane drainage technology and prevent gas disaster accidents in coal mines.     

Can tensile stress develop in fractured multilayers under compressive strain conditions?

To contribute to the understanding of how opening-mode fractures (joints) form and open or close at depth in layered rocks, we present a 2D numerical study aiming to determine whether tensile stress can develop in pre-fractured elastic multilayers submitted to biaxial compressive strain conditions.First, we investigate the role of the elastic and geometrical properties of the layers on the development of tensile stress in models with five bonded layers and containing one open fracture in the central layer. Our results indicate that, in absence of elastic contrast (in Young's modulus) between the layers, no tensile stress develops in the models. However, when the fractured layer is stiffer than the two adjacent layers directly above and below, a lobe of horizontal tensile stress develops centered on the pre-existing fracture. The creation of this tensile stress is contingent upon the partial closing of the fracture. The levels of tensile stress and the thickness of the lobe of tensile stress increase logarithmically with an increase in the elastic contrast and are systematically larger for a larger Soft/Stiff ratio (ratio of the total thickness of the soft layers with the total thickness of the stiff layers).Second, we investigate the role of fracture interaction in the development of tensile stress in models containing a pair of open fractures. We observe that the levels of tensile stress in the region between the fractures are systematically higher than those observed in identical models containing a single fracture. This increase in tensile stress is very large for small elastic contrasts between the layers but diminishes when the elastic contrast increases. Furthermore, the spacing between the pre-existing fractures plays an important role in the stress distribution in the region between them. When the fracture spacing is equal to or lower than 1.15 times the height of the fractured layer for the experimental conditions chosen, the lobes of tensile stress centered on the fractures coalesce. This results in the formation of vast areas of tensile stress in models under remote compressive loading conditions. Such tensile areas are likely to allow the initiation and propagation of subsequent opening-mode fractures.The results obtained provide new insights into the formation of joints in layered rocks in compressive environments, with important consequences on fluid flow.     

卸荷条件下岩石变形特征及本构模型研究

基于岩石试件的卸荷试验,研究了卸荷条件下岩石的应力—应变全过程曲线和破裂特征,研究表明：卸荷过程中岩石向卸荷方向回弹变形强烈、扩容显著、脆性破坏特征明显;卸荷条件下岩石破坏具有较强的张性破裂特征,各种级别的张裂隙发育,其剪性破裂面追随张拉裂隙发展。基于试验认识,将卸荷条件下岩石的本构模型分为4段：弹性段、卸荷屈服段、峰后脆性段及残余理想塑性段;通过体积应变εv将Griffith和Mohr Coulomb屈服准则结合起来,建立卸荷条件下岩石的屈服准则,并在应变空间建立相应的卸荷非线性屈服段本构方程;运用应变软化理论,建立卸荷岩石峰后脆性段本构方程。     

宁夏羊场湾煤矿浅埋煤层开采地面塌陷发育规律及形成机理

以宁夏羊场湾煤矿Y110207工作面为研究对象,采用无人机遥感技术、野外调查与有限差分软件模拟方法研究浅埋煤层开采的地面塌陷类型、发育规律及其形成机理。(1)浅埋煤层开采地面塌陷以地表裂缝发育为主,地表破坏严重。(2)平行切眼裂缝间隔性出现,展布于整个工作面内,间隔距离为10～120m,局部裂缝形成错台高度约为15cm。平行顺槽裂缝为拉张型裂缝,发育在顺槽至外围一定范围。(3)采煤活动导致地下形成采空区,上覆岩层发生移动破坏,破坏区分为剪切破坏区、拉张破坏区及剪-拉破坏区,分别对压应力区、拉应力区和压-拉转化区。(4)当应力扰动传递至地表,应力值超过覆盖层抗拉强度时地表产生裂缝。随着工作面推进,覆岩内部裂缝带上行裂缝与地表下行裂缝贯通,形成错台。研究成果丰富了该区浅埋煤层的地面塌陷理论知识,为地面塌陷防治提供了理论依据。     

略论破裂构造的形成机制

“破裂”一词比之“断裂”有着更为广泛的涵义,其论述必然牵涉到破裂构造的许多问题。本文试图在前人工作基础上,通过实验观察,着重讨论剪切破裂的形成机制和模式。 岩石破裂是由于本身存在着弹性应变能的缘故。本世纪六十年代初,罗贝尔思(J.C. Roberts)等人,已将其原理应用到节理裂隙力学机制的分析上。     

基于连续-非连续方法的裂隙破坏与相互作用研究

岩体内裂隙等非连续结构面对岩体的强度及变形等力学特性有着显著的影响,研究岩体裂隙起裂、扩展、相互作用和贯通机制,对工程岩体力学行为的表征和工程性能的评价十分重要。本文基于连续介质力学模型的离散元方法,通过考虑裂隙分布、模型加载条件及其与裂隙产状的关系,建立了一系列裂隙力学计算模型,研究了不同模型裂隙扩展演化特征和岩体破裂机制,分析了岩体裂隙扩展规律及其对岩体破坏路径和强度的影响,研究结果表明:（1）裂隙岩体模型加载条件下的破坏起裂点、最终贯通破坏特征及损伤分布受控于裂隙的产状及其与最大主压应力取向角度大小及围压大小。（2）裂隙弱面走向与最大主压应力取向斜交时,裂隙弱面在加载条件下其端部裂隙扩展、贯通破坏表现比较明显,反之,当裂隙弱面走向与最大主压应力取向一致时,裂隙弱面被动影响裂隙模型内新生裂隙的萌生、扩展和贯通模式,自身未出现新的扩展破坏。（3）裂隙数目的增多和围压的增大会显著增加模型内部剪切裂缝的数量和模型破坏后的破碎程度,模型内部的损伤区域主要围绕破裂面呈滑移线型交叉分布,非破裂面区域损伤呈条带状X型分布。（4）裂隙弱面走向与最大主压应力取向斜交时,裂隙对岩体模型强度的弱化程度高于裂隙弱面走向与最大主压应力取向一致的情况,而裂隙模型破坏后的残余强度则正好相反。     

Transitional–tensile fracture propagation: a status report

One model for the development of hybrid shear fractures is transitional–tensile fracture propagation, a process described as the in-plane propagation of a crack subject to a shear traction while held open by a tensile normal stress. Presumably, such propagation leads to a brittle structure that is the hybrid of a joint and a shear fracture. Crack–seal veins with oblique fibers are possible candidates. While these veins clearly show shear offset, this is not conclusive evidence that a shear traction was present at the time of initial crack propagation. Many recent structural geology textbooks use a parabolic Coulomb–Mohr failure envelope to explain the mechanics of transitional–tensile fracturing. However, the laboratory experiments cited as demonstrating transitional–tensile behavior fail to produce the fracture orientation predicted by a parabolic failure envelope. Additional attempts at verification include field examples of conjugate joint sets with small acute angles, but these conjugate joints form neither simultaneously nor in the stress field required by the transitional–tensile model. Finally, linear elastic fracture mechanics provides strong theoretical grounds for rejecting the notion that individual cracks propagate in their own plane when subject to a shear traction. These observations suggest that transitional–tensile fracture propagation is unlikely to occur in homogeneous, isotropic rock, and that it is not explained by a parabolic Coulomb–Mohr failure envelope as several recent structural geology textbooks have suggested.     

干湿交替作用后砂岩破裂过程实时观测与分析

采用岩石破裂全过程的细观力学试验系统进行了天然状态及干湿交替作用后的完整和含预制圆孔砂岩试件的单轴压缩破裂试验。利用显微镜观察了砂岩裂纹萌生、扩展以及试件的失稳破坏过程。经过干湿交替作用后岩样更容易在较低的应力状态下发生裂纹的扩展,并且裂纹扩展的初始阶段一般表现为拉破裂。完整试件的裂纹是随机性的三维扩展,而有预制圆孔试件一般是从预制圆孔周边开始扩展。天然状态下的预制圆孔试件在起裂时从圆孔的压应力集中区产生剪切裂纹;而干湿交替作用后,起裂时可能先从拉应力集中区产生拉裂纹,并且裂纹演化过程更为复杂一些,破裂形式多样化。与天然状态相比,经过干湿循环作用后的砂岩破坏脆性降低,强度等力学参数也随干湿交替次数增加而逐渐降低。     

Effects of pressure on the micromechanics of faulting in San Marcos gabbro

Samples of San Marcos gabbro were deformed to both pre- and post-failure stages at confining pressures of 100, 250 and 350 MPa. Scanning electron microscope observations were performed on the samples after ion-milling. The unstressed gabbro has insignificant crack porosity but a high volume fraction of healed and sealed cracks. Extensive stress-induced cracking at high angles to the maximum compressive stress is observed at the two higher pressures. Such high-angle cracks are mostly associated with the pre-existing healed cracks or cleavage planes, and they control the geometry and length scale of crack networks in plagioclase. Intensive kinking in biotite is widely observed at 350 MPa pressure, which together with the crack networks control the brittle fracture process. Stress-induced cracks in samples deformed at 100 MPa pressure are predominantly subparallel to the maximum compressive stress, and their coalescence leads to shear localization.Quantitative measurements of crack density as a function of differential stress indicate that crack interaction is appreciable beyond about 90% peak stress for all three pressures. The data also show significant decreases in stress-induced anisotropy as pressure increases, suggesting that there may be a threshold pressure beyond which stress-induced anisotropy vanishes with the onset of cataclasis.     

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

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

Life-time prediction for rocks under static compressive and tensile loads: a new simulation approach

The life time or time to failure of rocks under load is governed by microstructural defects, like microcracks, voids etc. The life time can be predicted either by empirical exponential laws or physical laws based on damage and fracture mechanics. The proposed numerical model is based on subcritical crack growth using the linear elastic fracture mechanical approach and is implemented as a numerical cellular automate. The algorithm considers both tensile and shear fracturing. Each cell contains a microcrack of random length according to a given probability function. Fracture growth is controlled by the Charles equation. Macroscopic cracks are the results of the coalescence of growing microcracks. Within the numerical approach elasto-plastic stress redistributions take place. If the stress intensity factors have reached the critical values or the microcrack has reached the zone dimension, the zone is considered as fractured and residual strength values are assigned. The proposed approach was applied to rock samples under uniaxial compressive and tensile loads (creep tests). Successful results were obtained in respect to the predicted life time, damage evolution and the fracture pattern. Conclusions for further improvements and extensions of this methodology were drawn.     

节理倾角对黄土力学特性影响试验研究

黄土节理种类多样、角度多变,是影响黄土力学性质的重要因素.为研究节理对黄土力学性质影响,开展含预制节理黄土试样单轴与三轴压缩强度试验,分析了非贯通节理对试样剪切破坏模式影响,探讨了节理倾角对黄土强度与变形特性影响规律.结果表明:节理试样剪切破坏均为压剪破坏,破坏模式可分为节理面与破裂面贯通型和节理面与破裂面斜交型,压力...     

锥形聚晶金刚石复合片齿破岩特征与机制研究

锥形聚晶金刚石复合片(polycrystalline diamond compact,PDC)齿是一种具有较强抗冲击性和耐磨性的新型PDC齿,在坚硬、强研磨性和软硬交错地层中取得了非常好的钻井提速效果。为了揭示锥形齿破碎硬岩机制,开展锥形齿破碎花岗岩试验与数值模拟研究,分析了切削深度和前倾角对锥形齿切削力和破岩比能的影响规律,采用高速摄像机和透明K9玻璃观测了锥形齿作用下岩屑形成过程及微裂纹萌生与扩展过程,通过数值模拟分析了破岩过程中岩石应力响应与损伤演化特性,结合对切削槽和大尺寸块体岩屑表面形貌及断口微观特征分析,建立了锥形齿破碎花岗岩机制模型。结果表明,锥形齿破碎花岗岩的过程可以分为挤压成核和块体崩裂两个阶段,前倾角对岩石破碎过程影响较小,切削深度的影响显著;锥形齿周围的裂纹主要由压实核、纵向裂纹和横向裂纹组成,纵向裂纹和横向裂纹扩展的最大深度分别为切削深度的6.69倍和4.53倍;齿尖周围压应力集中,岩石发生压剪破坏,压应力区外围形成弧形条带状拉应力区,并在齿尖及压应力区边界处诱导出拉伸微裂纹;微裂纹向齿前扩展形成弧形拉伸主裂纹,发生块体岩屑崩裂,提高破岩效率,向岩石内部扩展劣化...     

Numerical models of hydraulic fracturing and the interpretation of syntectonic veins

The hydraulic fracture is modelled as an ellipse in an infinite elastic medium with an internal fluid pressure and loaded under biaxial stresses at infinity. The available stress function for this model has been evaluated numerically, and the magnitudes of the stresses generated around the crack calculated for a variety of loading conditions and crack orientations. Fracture initiation is predicted from the Griffith maximum tensile stress criterion. The location of the maximum tensile stress around the crack is recorded and it is found that for many conditions of applied stresses and crack fluid pressure, the hydraulic shear fracture has a symmetrically developed maximum tensile stress and fracture initiation will occur by growth along the direction of the crack. It is also predicted that fracture initiation will occur when the ratio of fluid pressure to applied least principal stress is considerably less than one. The elastic strain energy fields around elliptical hydraulic flaws have been calculated, and in particular, the change in strain energy upon introduction of a small flaw, and the change in strain energy upon growth of this flaw, have been investigated. The results allow an evaluation of the second part of the Griffith criterion-that fracture growth is accompanied by a decrease in strain energy-for hydraulic fractures. Changes in strain energy with small increases in fluid pressure provide a physical basis for dilatancy hardening and fracture instability. Quasi-static growth from a flaw is modelled by calculating changes in strain energy for unit increases in half length. The distinction between fractures which show an increasing and a decreasing rate of change in strain energy with increasing length, and between fractures which may only extend spontaneously for short distances and those which may show extensive spontaneous growth on the basis of the rate of change of strain energy with length, is made. A gradual drop in crack fluid pressure once the threshold for fracture initiation has been passed may promote the extent to which spontaneous crack growth occurs.The formation of syntectonic veins, particularly in rocks being deformed under low grade metamorphic conditions, is often the most abundant evidence of natural hydraulic fracturing in rocks. Commonly observed geometric features of syntectonic veins-length, simple tapering, symmetric and asymmetric forking, branching, irregular zig-zag traces, en échelon patterns—are discussed primarily with reference to the strain energy model for growth established, and the geometric variation is interpreted in terms of variation in applied stress and fluid pressure conditions and the rate of change of stored strain energy with crack growth. In particular, terminal branching arises when the minimum stress changes from a symmetric to an asymmetric location at the tip of a growing shear fracture, and terminal forking results when there is an increase in the energy release rate during crack growth, and may be symmetric or asymmetric depending on the location of the minimum stress at the crack tip at the time of forking.     

Primary and remote fracture around underground cavities

The development of fractures around cavities in rock is examined using physical and numerical models. In uniaxial compression, fracture proceeds from primary through remote to compressional (slabbing) fracture. There is a wide stress interval where the only event is the stable propagation of the primary fracture. When the crack length becomes equal to about the size of the diameter of the cavity, the critical stress condition shifts to either side of the primary fracture, leading to the nucleation of remote fractures. Almost simultaneously, failure at the compressive stress concentration follows. With increasing stress, fracture in the remote region is seen to extend away from and towards the compressive region, eventually merging with the slabbing fractures. Increasing the confining pressure shortens the stress interval during which primary fracture is the only process and shifts failure to the remote and the compressional fracture zone.     

热应力影响下干热岩水压致裂数值模拟

干热岩水压致裂过程中低温诱导热应力与注入水压共同影响裂缝的萌生与扩展。首先通过THM耦合分析了低温压裂液注入过程中注入水压与热应力的相互作用及其对裂缝萌生的影响,随后建立描述岩石细观结构的THMD耦合模型对热应力影响下高温岩石水压致裂过程进行初探。结果表明：低温压裂液注入高温岩石产生的热应力包括岩石自身温度梯度形成的热应力与岩石颗粒非均匀膨胀导致的热应力,并在井筒周围呈现为拉应力。高注入压力将抑制热应力导致的多裂缝萌生,井筒附近热应力的存在对注入压力也具有削弱作用。基岩温度升高,裂缝萌生阶段更多裂缝在井筒附近起裂,缝网沿最大地应力方向的扩展速度减慢,但改造规模增加,同时多裂缝的存在也使得裂缝延伸压力增加。     

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

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

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

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

Microplane damage model for jointed rock masses

The paper presents a new microplane constitutive model for the inelastic behavior of jointed rock masses that takes into account the mechanical behavior and geometric characteristics of cracks and joints. The basic idea is that the microplane modeling of rock masses under general triaxial loading, including compression, requires the isotropic rock matrix and the joints to be considered as two distinct phases coupled in parallel. A joint continuity factor is defined as a microplane damage variable to represent the stress‐carrying area fraction of the joint phase. Based on the assumption of parallel coupling between the rock joint and the rock matrix, the overall mechanical behavior of the rock is characterized by microplane constitutive laws for the rock matrix and for the rock joints, along with an evolution law for the microplane joint continuity factor. The inelastic response of the rock matrix and the rock joints is controlled on the microplane level by the stress–strain boundaries. Based on the arguments enunciated in developing the new microplane model M7 for concrete, the previously used volumetric–deviatoric splits of the elastic strains and of the tensile boundary are avoided. The boundaries are tensile normal, compressive normal, and shear. The numerical simulations demonstrate satisfactory fits of published triaxial test data on sandstone and on jointed plaster mortar, including quintessential features such as the strain softening and dilatancy under low confining pressure, as well as the brittle–ductile transition under higher confining pressure, and the decrease of jointed rock strength and Young's modulus with an increasing dip angle of the joint. Copyright © 2014 John Wiley & Sons, Ltd.