不同加载速率下岩石强度特性的扰动状态概念理论分析

岩石的强度特性受很多因素的影响,加载速率是其中很重要的影响因素之一.利用扰动状态概念理论,建立了非线弹性岩石材料的本构模型.采用RMT-150B岩石力学试验机对砂岩试件进行了四种不同加载速率下的单轴压缩试验,得到了不同加载速率下砂岩的应力-应变全过程曲线,并分析了加载速率对岩石强度特性的影响.通过建立的扰动状态本构模型对岩石的应力-应变关系进行了模拟.相关研究及分析比较表明,该模型与试验结果较为一致.     

单裂隙砂岩单轴压缩的中等应变率效应颗粒流模拟

加载速率对裂隙岩体的力学性质及变形破坏均有重要影响。利用二维颗粒流程序PFC2D开展了不同倾角非贯通单裂隙砂岩试件的单轴压缩试验,研究了中等应变率对裂隙砂岩应力-应变曲线特征、裂隙尖端应力状态、特征应力状态、岩体损伤及裂隙扩展等力学响应的影响规律。裂隙岩体应力-应变曲线呈现明显的波动性,定义应力突变指标 对应力突变型波动剧烈程度进行了定量统计分析：随应变率的增加,曲线应力突变波动越剧烈,且峰后明显大于峰前;随裂隙倾角的增大,波动幅度峰前增大,而峰后减小。裂隙尖端破裂应力随应变率增大均有所提高,随裂隙倾角的增大,切向剪应力 总体上呈增加变化,而法向应力 明显减小。尖端破裂时岩样加载应力 、岩样临界扩容应力 及峰值应力 均随应变率增大而增大。裂隙尖端的破裂可立即引起岩体扩容,一般应变率越低,岩体裂隙尖端破裂点 和扩容点 越接近峰值强度 。随着应变率的提高,损伤裂纹及宏观裂隙类型越多,岩体试件损伤破裂程度越强,特别是试件端部效应愈显著。裂隙首先以I型翼裂纹在其尖端起裂,而I型翼裂纹的扩展长度与加载速率与裂隙倾角具有较强的相关性。     

重塑黄土的加载速率效应研究

为了揭示不同加载速率条件下重塑土的变形与强度特性对工程实践的影响,通过不同围压、含水率条件下重塑土样在4种加载速率条件下的三轴CU试验,对重塑土在不同状态下的加载速率效应进行分析,并从土的结构性入手对其与原状土加载速率效应的差异进行讨论。基于重塑土样的三轴试验成果,得到了加载速率对重塑土变形与力学特性的影响规律:重塑土的应力-应变曲线在不同加载速率条件下均表现出应变硬化特性;破坏强度随着加载速率的增大而表现出先增大后减小的特性,存在临界加载速率;随着加载速率的增大,黏聚力c先增大后减小,也存在临界加载速率,而内摩擦角φ则保持了小幅减小的趋势。     

基于颗粒离散元模型的不同加载速率下花岗岩数值试验研究

不同加载速率条件下岩石的力学特性,对于其动载下破裂内在机制的研究具有积极的意义。基于颗粒流理论,通过黏结颗粒模型（bonded particle model,简称BPM）虚拟实现不同加载速率0.001～0.500 m/s下花岗岩单轴压缩和巴西劈裂试验,定量分析加载速率对应力-应变、破裂形态、应变能率及声发射的影响。结果表明：单轴抗压强度和抗拉强度及其对应峰值应变随加载速率增加而非线性增长;单轴压缩作用下,随加载速率增加,试样由单一斜截面破坏向多斜截面破坏转变,且主控裂隙带宽度急剧增大,由裂纹数量及水平向高应变率区域变化规律可明显看出,试样破坏程度随着加载速率增加而逐渐加剧;巴西劈裂作用下试样从一条主控裂隙向多条主控裂隙转变,且裂纹向圆盘试样两侧边缘部分延伸,破坏程度加剧;单轴压缩和巴西劈裂作用下,声发射事件及应变能率均随加载速率增加而呈现出非线性增长趋势。     

岩石张拉裂缝与充填胶结裂隙的相互作用

工程开挖扰动、卸荷和水力破岩产生的工程裂缝与先存裂隙的相互作用在岩石工程中十分普遍,建立适用于工程裂缝遇先存裂隙的扩展准则是解释工程裂缝扩展机制和研究岩体弱面剪切滑移的关键。采用两种不同强度的混凝土充填材料模拟强、弱胶结的先存裂隙,在只考虑准静态张拉应力作用下,依据半圆三点弯法测定了充填胶结型岩石的Ⅰ型静态断裂韧度;通过分析工程裂缝扩展轨迹,获得了临界启裂角与应力逼近角的关系;进一步探讨了工程裂缝与先存裂隙的相互作用扩展模式与判别准则。结果表明：张拉应力作用下工程裂缝与先存裂隙的相互作用受应力逼近角、启裂逼近距离和充填物胶结强度的共同影响;充填胶结砂岩的静态断裂韧度随应力逼近角增大先增大后减小;启裂逼近距离较大时,断裂韧度随应力逼近角变化不大,而启裂逼近距离较小时,断裂韧度随应力逼近角增大呈现先增大后减小的规律;应力逼近角影响充填胶结岩石的断裂韧度,但影响程度不如启裂逼近距离和充填物胶结强度,且存在一个极限影响距离。即使在张拉应力作用下,岩石材料的裂隙前端同样存在剪切局部化,因此,工程裂缝是否穿过先存裂隙取决于应力逼近角、启裂逼近距离和充填物胶结强度共同作用下裂隙面的抗剪强度和先存裂隙...     

含交叉裂隙节理岩体单轴压缩破坏机制研究

对含交叉裂隙相似材料试件进行单轴压缩试验,用Abaqus对含交叉裂隙节理岩体进行应力分析,研究了含交叉裂隙节理岩体单轴压缩下的破坏机制。结果表明：交叉裂隙的主裂隙与加载方向呈30°、45°时,主裂隙是裂隙扩展及破坏的控制裂隙,此时含交叉裂隙岩体强度高于含单向裂隙岩体,是由于主裂隙裂尖最大环向应力σθmax低于单向裂隙;主裂隙与加载方向呈0°或90°时,次裂隙是控制裂隙;主裂隙与加载方向呈0°时,大部分次裂隙裂尖σθmax及裂尖最大应力强度因子KⅡmax高于单向裂隙,因此,此时大部分含交叉裂隙岩体强度低于含单向裂隙岩体;在主裂隙与加载方向呈90°工况组中,主、次裂隙夹角为45°时,次裂隙裂尖KⅡmax最大,所以此工况组中此时强度最低。     

岩石破裂行为的实验研究

岩样的破裂行为、破坏过程和参数测试是裂隙断裂构造研究的基础和依据。实验岩石样品采自四川东部和新疆北部地区,为测试准确起见,对岩样进行了应力等值线的有限元法计算。通过单轴和三轴实验的岩样破坏观察和应力应变曲线对比,将岩石的破裂行为、应力应变划分为四个阶段,即裂隙压密阶段、弹性变形阶段、微观劈裂阶段和宏观破裂阶段。基于单轴抗压实验岩石劈裂—破裂—碎裂发展过程的微观分析,可以看出宏观破裂主要是沿岩样原有的隐裂隙、临界裂隙发育的,许多新裂隙则主要是在宏观破裂阶段产生的。     

低应力下岩石裂隙水力特性

针对地表附近30～50 m的岩体裂隙处于一种低围压状态（约1 MPa左右）这一情况,提出了低应力状态下裂隙变形的两点假设：（1）裂隙面处于线弹性状态;（2）隙宽减小量相对于初始力学隙宽是一个微小量。结合Barton公式,建立了低法向应力下裂隙渗流的半经验理论公式,并从实验的角度验证了该公式的合理性。实验结果和理论分析表明：低法向应力作用下,裂隙渗透系数随法向应力的增加而线性降低。此外,给出了一种间接确定裂隙初始力学隙宽的方法。     

不同加载速率Ⅰ型预制裂纹花岗岩断裂特征研究

以CCD相机配合高速相机搭建试验变形观测系统,开展不同加载速率下含Ⅰ型预制裂纹花岗岩试件的三点弯试验。利用数字散斑相关方法对试验过程采集的数字散斑图像进行分析,对不同加载速率作用下含Ⅰ型预制裂纹花岗岩试件裂纹扩展过程中的变形演化特征、裂纹扩展速率、裂纹尖端张开位移、裂纹张开角以及初始断裂应力强度因子进行研究,研究结果表明：（1）花岗岩试件初始裂纹扩展速率随着加载速率的增加总体呈线性增加趋势;最大扩展速率随着加载速率增加呈现先较快增加后缓慢增加的变化趋势。（2）不同加载速率作用下,花岗岩试件预制裂纹尖端张开位移均呈现非线性缓慢增长-迅速增长-线性增长的变形趋势,且线性增长阶段曲线斜率随着加载速率的增加而增大。（3）随着加载速率的增加,花岗岩试件初始裂纹张开角先逐渐减小后稳定在0.1o～0.13o。（4）随着加载速率的增加,花岗岩试件初始断裂应力强度因子呈指数增长。     

含折线型裂隙砂岩试件翼型裂纹起裂与扩展机制研究

岩体内部赋存的裂隙很多表现为折线型,为探究这类岩体的断裂机制,制备含折线型裂隙砂岩试件并对其进行单轴压缩试验。采用数字图像相关(DIC)方法计算加载过程中的变形场演化,根据新生裂纹两侧的位移差异识别裂纹类型;运用扩展有限元法(XFEM)模拟断裂过程,根据应力分布特征解释翼型裂纹起裂与扩展机制。DIC计算结果表明,新生裂纹处出现应变局部化带,裂纹两侧发生相对分离;含直线型和折线型裂隙砂岩试件的翼型裂纹分别萌生于预制裂隙端部以及折角处,这是因为裂隙几何形态会改变拉应力集中位置;含折线型裂隙砂岩试件的起裂应力小于含直线型裂隙砂岩试件,这是因为相同加载条件下前者的最大拉应力值更大;这2类试件的裂纹扩展均是由于裂纹尖端集中的拉应力引起的,裂纹依然呈张开状态;裂隙几何形态未改变试件的最终破坏模式,均表现为对角剪切破坏。     

Initiation Stress of Crack

In this paper, through analyzing the subcritical crack growth, a new formula is presented which can be used to calculate the crack fracture initiation stress at low and middle stress rates. By way of integration of the stress intensity factortime relation of crack under impact load, is attained a formnla from which the crack fracture initiation stress at middle and high stress rates cen be calculated. Combination of the two formulae leads to a general expression for caloulating the crack fracture initiation stress at all stress rate, low or high. The critical unstable fracture of material which contains many cracks is also analyzed in this paper, and a new viewpoint is put forward that the σfi-σ curve of material at high stress rate is an envelop of all σfi-σ curves of crack in the material, which leads to a conclusion that crack fracture initiation stress is proportional to the cube root of loading rate.     

煤层水力压裂应力与裂隙演化的细观规律

水力压裂作为煤层强化增透技术的一种,其应力演化特征及裂隙形态与扩展范围的判断尤为重要。采用离散元数值方法,以导向压裂为背景,建立水力压裂流固耦合模型;通过应力路径、裂纹热点图等手段,探究水力压裂过程中压裂排量、泊松比、天然裂隙密度对应力演化和裂隙演化的影响及其细观规律。结果表明:不同压裂排量下的应力演化方向及最终应力路径曲线形状有着明显的不同,低排量下裂隙附近的应力比值逐渐增大,而在高排量下先增大后减小;煤层泊松比越大,平均压裂半径越低,但对起裂时间及裂隙的扩展形态影响不明显;天然裂隙的发育情况对水力裂隙的扩展起着关键性作用,高裂隙发育煤层水力裂隙扩展的方向性无法预测,应力演化方向会出现反转现象;压裂过程中不同区域的应力演化特征能够反映出裂隙的扩展状态,现场可通过监测压裂区域附近应力变化,判断水力压裂缝网的扩展范围。      

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.     

Experimental investigation of the strain rate effect on crack initiation and crack damage thresholds of hard rock under quasi-static compression

Acta Geotechnica - The crack closure stress σcc, crack initiation stress σci and crack damage stress σcd are important mechanical parameters to describe the brittle failure mechanism...     

华北东部深部岩溶及煤矿岩溶水害特征

超前区域注浆是煤层底板灰岩水害防治的有效手段,而注浆工程中的劈裂注浆过程是决定注浆效果的关键环节,但由于对受注灰岩地层裂隙起裂机制认识不清,致使劈裂注浆过程中注浆压力、浆液水灰比等参数缺少有效控制,难以保证注浆效果。针对上述问题,利用颗粒元(Particle Flow Code,PFC)数值模拟软件,考虑浆液水灰比、地应力水平、弱面和裂隙的倾角和宽度等因素,开展奥陶系灰岩顶部劈裂注浆数值模拟计算。结果表明：在弱面和裂隙条件下起裂压力均随浆液水灰比(1∶1、2∶1、3∶1)的增大而减小,随最大主应力与最小主应力差值(9、12、15 MPa)的减小而增大,随弱面或裂隙宽度(3、8、15 mm)及其与最大主应力夹角(30°、60°、90°)的增大而减小;起裂裂隙沿平行于最大主应力方向延展;起裂压力值随弱面与基质强度比(0.30、0.03)的减小而减小,当弱面强度较高时,起裂压力大于裂隙条件下的起裂压力,而且沿着基质起裂;当注浆孔置于弱面两端或裂隙中间时,起裂裂隙沿弱面或裂隙的两端起裂,当注浆孔置于弱面中间位置时,起裂裂隙沿中间位置的基质起裂。研究结果有助于底板水害超前区域改造劈裂注浆的有效控制,指导注浆参数的选取,为解放深部煤炭资源提供技术支撑。

     

Modeling soil desiccation cracking by analytical and numerical approaches

An energy approach is proposed as a complement to the stress approach commonly considered for investigating soil desiccation cracking. The elastic strain energies before and after crack initiation are estimated by both numerical and analytical solutions. The energy released by cracking is then compared with the fracture energy to discuss crack initiation conditions. This leads to combined energy and stress conditions for crack initiation following Leguillon's theory. An approximate analytical solution is derived from a variational formulation of the porous elastic body equations. A cohesive zone model and finite element code are used to simulate crack propagation in an unsaturated porous body. This analysis shows that the energy criterion is reached before the stress criterion, and this can explain unstable crack propagation at the beginning. The approximate analytical solution allows predicting correctly the crack depth and opening in its initiation stage.     

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

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

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

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

Mapping the state of fracture around cavities

The state of stress around cavities is heterogeneous. Consequently, the state of fracture may also vary from point to point. Under compressive loading, cavity may be in one of four possible states: pre-fracture (pre-microfracture initiation), microfracture propagation (initiation to the onset of dilata dilatancy to failure) and post-failure. These four states are separated by the crack initiation stress, the crack damage stress and the failure stress. fracture events with confining pressure is examined with reference to three intact rocks: a brittle granite, a semi-brittle limestone and a ductile sal

The maximum principal stress (σ₁) at crack initiation, crack damage (onset of dilatancy), yielding and failure are established as a function stress (σ₃). For a single intact rock, all four fracture events can be represented using one function (the Rocker function) with a single fac various fracture states.

The proposed fracture criteria, based on the experimental data, are combined with the existing state of stress to prepare a fracture map around an elli intact Lac du Bonnet granite. The state of stability is expressed through a newly defined stability factor, the unconfined strength ratio (USR), wh to the traditional safety factor. In contrast to the conventional safety factor in rock mechanics (SFstrength/σ₁), which is σ₃ 3–σ₁, space of the stability diagram.     

Similarity solutions for wedge-shaped hydraulic fractures driven into a permeable medium by a constant inlet pressure

Similarity solutions are derived for wedge-shaped hydraulic fractures driven by a constant inlet pressure P₀ into a permeable medium under a uniform confining stress σ. These results describe the seepage-dominated regime in which most of the injected fluid is lost into the permeable walls of the fracture; they complement previous results for the capacitance-dominated regime in which seepage is negligible. Fracture propagation velocity is obtained as an analytical function of fracture length, driving pressure, confining stress, material properties and a single separation constant or eigenvalue which is determined numerically. Self-similar profiles of pressure, opening displacement and fluid velocity along the fracture are presented, together with the self-similar isobars of the two-dimensional pressure field within the permeable medium. Comprehensive results are reported for laminar or turbulent flow of a constant-compressibility liquid or an ideal gas driven by overpressures (P₀?σ)/σ ranging from 10^?2 to 10².