岩体变形和破坏的各种应力途径

在真三轴应力状态下进行了各种应力途径的岩石变形和破坏实验。实验表明,三个主应力的变化(最大主应力的增加,最小主应力的减小,围压的减小以及中间主应力的增加或减小)都能引起岩石变形和破坏。而且岩石的强度、延性和脆性、体积变化和破坏前兆等都与应力途径有关。据根试验结果探讨了岩体工程中的应力途径问题。     

Numerical modeling of plastic deformation and failure around a wellbore in compaction and dilation modes

In many wellbore stability analyses, the ability to forecast both the occurrence and extent of plastic deformation and failure hinges upon a fundamental understanding of deformation mode and failure mechanism in the reservoir rock. This study focuses on analyzing plastic zones, localized deformations, and failures around a borehole drilled overbalanced or underbalanced through a highly porous rock formation. Based on several laboratory experiments, porous rocks are prone to deform under both shear-induced dilation and shear-enhanced compaction mechanisms depending on the stress state. The shapes of the deformation and failure patterns around the borehole are shown, depending on the initial stress state and the local stress paths. The inquiry of the local stress paths in the near-wellbore zone facilitates the understanding of the reasons for different types of failure mechanisms, including the mixed-mode and the plastic deformation structures. The modification of the 2D plane strain condition by imitating third stress in the numerical scheme helps us bring the stress paths closer to the real state of loading conditions. Our modeling reveals that the transition from isotropic to anisotropic stress state is accompanied by an increase in the deviatoric part of effective shear tensor that leads to the development of inelastic deformation, degradation, and subsequent rock failure. Particular interest is devoted to the modeling of strain localization especially in compaction mode around a wellbore and computing the amount of stress concentration at the tips of dog-eared breakouts. Stress concentration can result in a change in irreversible deformation mode from dilatancy to compaction, elucidating the formation of the shear-enhanced compaction phenomenon at the failure tips in the direction of the minimum horizontal stress.     

The influence of the intermediate principal stress on rock failure behaviour: A numerical study

The influence of the intermediate principal stress on rock strength has been studied comprehensively by previous researchers. However, the reason why rock strength firstly increases and subsequently decreases with the increase of intermediate principal stress is still unclear. In this paper, the mechanism of the intermediate principal stress effect on rock failure behaviour is revealed through a numerical method using the EPCA3D system (Elasto-Plastic Cellular Automaton). In this study, both homogeneous and heterogeneous rocks are considered. The heterogeneity of a rock specimen is modelled by introducing Weibull's statistical distribution. Two criteria, i.e. the Drucker–Prager and Mohr–Coulomb models, are used to determine whether a meso-scopic element in the rock specimen is in a failure state or not during the polyaxial stress loading process. The EPCA3D simulation reproduces the typical phenomenon of the intermediate principal stress effect that occurs in some rock experiments. By studying the EPCA3D simulated acoustic emission and complete stress–strain curves illustrating failure initiation, propagation and coalescence in the failure process of rocks, the essence of the intermediate principal stress effect is tracked. It is concluded that the heterogeneous stress distribution induced by the natural heterogeneity of rocks and the effect of the loading platen are two of the reasons producing the intermediate stress effect. Studies indicate that a moderate intermediate principal stress delays the onset of local failure, which in turn leads to an increase in the rock strength. However, once the intermediate principal stress reaches a certain value, local failure will be formed through the application of the intermediate principal stress. It is the number of failed elements in the pre-peak region that determines whether the rock strength decreases or not. The extent of rock strength reduction when the intermediate principal stress reaches a certain value is lessened with the increase in the minimum principal stress.     

三维应力作用下岩体层面效应分析

为得到层状岩体在真三轴压缩情况下的强度变形特征,利用数值分析方法建立计算模型,得到不同结构面倾角? 下试样的应力和变形响应,结果表明,(1) 加载初期,结构面倾角? = 0°和? = 90°试样对应的初始弹性模量最大;随着试样变形的增大,对应的弹性模量以? = 90°、60°、0°的顺序逐渐减小;? = 90°和? = 0°试样发生破坏的脆性特征比较明显;(2) 随着围压? 3 = σ2的增大,压缩强度? 1不断增大;围压对强度极值发生的位置影响不大;(3) 最小主应力与结构面走向平行情况下,随着中间主应力? 2的增大,试样压缩强度呈现先增大后减小的趋势;随着最小主应力? 3的增大,? 1逐渐增大;当? 3较小时,二者符合线性关系;随着? 3的增大,二者的关系逐渐呈现非线性特征;并且结构面滑移破坏型试样的非线性特征更加明显;(4) 最小主应力与结构面走向垂直的情况下层状岩体的压缩强度与中间主应力和最小主应力的关系均呈现显著的线性特征。     

Anisotropic Deformation of a Circular Tunnel Excavated in a Rock Mass Containing Sets of Ubiquitous Joints: Theory Analysis and Numerical Modeling

Determining anisotropic deformation surrounding underground excavations for tunnels is an intuitional task that involves many difficulties due to the inherent anisotropies in the strength and deformability of natural rocks. This study investigates joint-induced anisotropic deformation surrounding a tunnel via a numerical simulation that accounts for the mechanical behavior of intact rock, the orientations of joint sets, and the mechanical behavior of joint planes; this numerical simulation can model the complete stress–strain relationship with anisotropic rock mass characteristics. Simulation results demonstrate that the well-known excavation-induced stress variation–decrease in the radial component and increase in the tangential component–decrease shear strength and increase shear stress for the joint plane tangential to the tunnel wall, resulting in joint sliding failure and considerable shear deformation. This joint sliding failure and significant shear deformation account for the joint-induced anisotropic deformation surrounding a tunnel. When a rock mass has two joint sets with unfavorable joint orientations, the area with joint sliding failure can deteriorate mutually, resulting in large anisotropic deformation. Additionally, for a rock mass containing three joint sets with well-distributed orientations, joint sliding in various joint sets and associated stress variations can counter balance each other, resulting in less anisotropic deformation than those of rock masses containing one or two joint sets.     

拉张盆地不同埋深复合地层中隧道围岩变形破坏数值模拟研究

不同埋深的软硬岩层叠置复合地层变形破坏形式复杂而使得软岩大变形和硬岩岩爆位置相关关系不够明确。在拉张盆地中自重应力为主,侧压系数一般较小且多数在1.0以内。本文以在水平和铅垂叠置复合地层中TBM开挖圆形断面隧道为例,采用有限差分程序FLAC^3D对拉张盆地中不同埋深、不同叠置型式复合地层中TBM开挖后的三维弹塑性位移变形、主应力和塑性破坏分布变化特征展开数值模拟研究。模拟结果表明,围岩变形主要发生在软岩层地层中,埋深超过800 m后沿隧道轴向软岩大变形藕节状分段开始显现;随埋深增加,硬岩稳定性变差,顶拱位移增加尤其明显;随埋深增加,软岩和硬岩地层之间主应力差异变小,硬岩中储能明显;埋深越大塑性区分布范围越大;埋深较小时岩层以拉张破坏为主,埋深较大时以剪切破坏为主,两种状态的转换埋深(临界深度)约为800 m。由此对拉张盆地中深埋界限值给予了理论验证。     

煤炭开采对煤层底板变形破坏及渗透性的影响

煤层底板变形破坏除受地质因素控制外,还受开采因素影响。通过试验和理论分析,系统研究了煤炭开采对回采工作面底板应力、应变和破坏及渗透性的影响。研究结果表明,不同岩性岩石的渗透性在全应力-应变过程中为应变的函数,在微裂隙闭合和弹性变形阶段,岩石的原生孔隙和裂隙容易被压密,岩石的渗透率随应力的增加由大变小明显,当应力增大至极限强度时岩石试件破坏形成贯穿裂隙,岩石的渗透率迅速增大至最大,不同岩性岩石存在一定差异性;随着回采工作面推进,煤层底板岩层在横向上划分为原岩应力区、超前压力压缩区、采动矿压直接破坏区和底板岩体应力恢复区4个区。煤层底板岩体的渗透性随着煤炭开采底板岩体变形破坏而呈规律性变化。      

砂质泥岩三轴卸荷流变力学特性试验研究

流变对工程岩体的长期变形稳定具有十分重要的影响,但目前关于卸荷流变的试验研究中通常只考虑恒轴压卸围压的应力路径,与隧洞开挖过程中围岩的应力调整过程存在一定的差距。以砂质泥岩为试验研究对象,设计进行了加轴压卸围压和恒轴压卸围压条件下的分级卸荷流变试验。试验结果表明：恒轴压卸围压和加轴压卸围压方案下岩样的流变变形趋势总体一致,但相同初始围压条件下,加轴压卸围压试样破坏的围压相对较高,偏应力相对较大,但长期强度与破坏应力的比值相对较小;在围压卸载至岩样临近破坏时,加轴压卸围压方案下岩样的流变应变增长速率明显较快,试样破坏的更加突然;恒轴压卸围压条件下岩样的破坏形态相对简单,一般只存在一条完整的剪切破坏面,而加轴压卸围压条件下岩样的破坏形态要复杂得多,除了控制性的剪切破坏面之外,还伴随有一定数量的次生剪裂纹和张拉裂纹,而且初始围压越大,试样的次生裂纹越多。因此,在隧洞围岩长期变形稳定分析过程中,单纯的恒轴压卸围压流变试验不能满足工程实际需求,应该尽量丰富岩石力学试验的应力路径,以便较好地模拟工程岩体应力的实际变化过程,研究成果为隧洞围岩的长期变形稳定性研究提供了较好的研究思路。     

泥质砂岩蠕变特性与锚固控制效应试验研究

通过泥质砂岩单轴压缩蠕变试验,揭示岩石蠕变变形与破坏的特点,反映了不同应力路径下的岩石蠕变速率与发展过程,建立了泥质砂岩的非线性蠕变模型,并通过对试件施加侧向刚性约束和锚固约束,以反映地下工程围岩在平面应变状态和支护状态下的蠕变特性。研究结果表明,锚固约束可以显著提高岩石发生稳定蠕变的应力阀值,减缓蠕变速率。根据改进的蠕变方程进行试验数据的拟合,确定了改进的西原模型蠕变方程中的各物理量。与单轴蠕变试验结果对比表明,该模型的计算结果与试验数据吻合较好,说明改进的西原模型可较好地反映岩石加锚后的蠕变特征。研究结果对深部地下工程围岩的变形与稳定控制研究具有重要的指导意义。     

水压影响岩石渐进破裂过程的试验研究

岩石的压缩过程伴随着裂纹扩展,在低、中围压条件下扩展的裂纹主要受到张拉作用而产生,张拉破坏是首要的作用机制。研究发现,岩石的渐进破裂过程受岩石的结构和构造影响,比如矿物成分、颗粒大小以及胶结情况等,另外,围压以及开挖扰动等外界因素对岩石的渐进破裂过程也有重要影响。基于试验方法探讨水压对岩石渐进破坏过程的影响,首先,阐述岩石渐进破裂过程各个阶段的特征;然后,对岩石渐进破裂指标--岩石的启裂强度 和损伤强度 的确定方法进行总结;最后,选取细粒的石英砂岩为研究对象,通过试验研究水压对岩石渐进破坏过程的影响。试验结果表明,相同围压条件下,随着岩样两端水压的增大,岩石的 有逐渐增大的趋势,而岩石的 和峰值强度 逐渐变小;随着围压逐渐增大,岩石的启裂强度 、损伤强度 及峰值强度 均逐渐增大。     

考虑软弱夹层厚度的岩体力学响应及破坏特征研究

作为一种典型的地质结构,软弱夹层与硬脆性岩体共同形成了围岩层状复合结构,进而显著影响着隧洞围岩的稳定。以往对含软弱夹层的复合岩石的研究多集中于单轴、双轴或常规三轴,对隧洞临空面处真三轴应力路径下的复合围岩力学性质和破坏特征缺乏分析讨论。通过制作的含不同厚度的软弱夹层复合岩样,探讨了软弱夹层厚度对隧洞临空面围岩力学响应和破坏特征的影响。研究表明：软弱夹层厚度显著影响着复合岩样峰值应力和应变,随着厚度的增大,软弱夹层上方岩块向临空面方向的滑移变形逐渐增大,软弱夹层压缩变形逐渐减小;复合岩样靠近临空面的岩石单元破坏模式随着软弱夹层厚度的增大逐渐由拉剪混合破坏转变为张拉破坏,且宏观裂隙数量和破坏范围均逐渐减小,而远离临空面的岩石单元则由剪切破坏逐渐转变为基本无损伤断裂;不同厚度的隧洞侧帮复合围岩的破坏区域均集中在软弱夹层及其上方围岩处,软弱夹层下方围岩则基本保持稳定;在应力分布方面,软弱夹层厚度越大,最大压应力越向深部软弱夹层处转移,而拉应力区分布范围越小,但拉应力区深度越大。     

中主应力对硬岩破裂机制的影响

采用数值方法模拟硬岩的三轴压缩试验,应用4个本构模型,即Mohr-Coulomb模型、Drucker-Prager模型、应变软化模型及考虑变形模量劣化的应变软化模型,研究中主应力对均质及非均质硬岩破裂机制的影响。结果表明：不论是均质还是非均质岩体,中主应力对采用Drucker-Prager模型的岩体强度影响较大,而对采用其他本构模型的岩体强度影响不大;当非均质岩体采用Mohr-Coulomb模型或Drucker-Prager模型时,中主应力对岩体破坏过程影响不大,但对其破坏模式有较大影响;当非均质岩体采用应变软化模型时,中主应力对岩体破坏过程及其模式均有较大影响。针对工程算例,采用不同的本构模型获得的均质或非均质岩体强度相差很大,除Drucker-Prager模型外,同一本构模型的均质岩体强度远大于非均质岩体。实际工程岩体对中主应力的响应是不同的,故在地下结构设计过程中应根据岩体特性选择合理的本构模型以保证工程安全。     

盐岩三轴卸荷力学特性试验研究

结合金坛地下盐穴储气库工程,以腔体围岩实际受力状态为研究依据,进行了盐岩卸围压力学特性试验,得到了盐岩卸围压过程的应力-应变关系、变形特征及其规律。试验结果表明,在卸围压初始阶段,试样的轴向和径向应变增加相对缓慢,且应变值和围压基本呈线性,随着围压继续降低,轴向与径向应变值急剧增加;卸载曲线与加载曲线相比,在最大主应力 相等的条件下,当卸载围压达到与常规三轴压缩围压相对应值时,对应的轴向和径向应变值较三轴压缩时应变值要大,表明卸荷试验能引起试样更大的变形,卸荷产生的扩容量比加载时扩容量更大,更容易导致试样变形破坏;盐岩卸围压表现为塑性变形特征,与其他硬岩卸围压属脆性破坏有较大区别。研究结果对盐穴储气库工程的稳定性评价和注采气过程具有一定的指导意义和参考价值。     

Characteristics of the unloading process of rocks under high initial stress

The unloading process of rocks under high initial stress is complex, and verifying the mechanism of the unloading process in the field or in a laboratory is not straightforward. In this study, the unloading process of rocks under high initial stress was characterised by a mathematical physics model, which was then implemented in the finite element program LS-DYNA for analysis. In particular, the implicit and explicit methods were performed in sequence in the finite element simulation of rocks with initial stresses. In the numerical simulation, the characteristics of the dynamic unloading process of rocks were investigated for various peak initial stresses, initial stress release paths and initial stress release rates (ISRRs). The numerical results indicated that the rock failure could be induced by the release of the initial stress; furthermore, there is a relationship between the magnitude of the unloading failure and the peak initial stresses, the initial stress release paths and the ISRRs. When the initial stresses were at the same level, the equivalent initial stress release rate (EISRR) was introduced to quantitatively describe the characteristics of the unloading process. Using the numerical results, the unloading failure process was characterised, and a method for the static stress initialisation-dynamic unloading of rock was developed.     

岩腔后退近水平砂-泥岩互层危岩稳定性研究

岩腔后退是三峡库区重庆段近水平砂-泥岩危岩稳定性改变的关键。基于岩腔平行后退理论,提出了岩腔泥岩后退的剪切破坏和拉破坏力学机制,考虑岩腔后退参数提出了3类危岩体稳定性计算方法,结合重庆万州太白岩危岩稳定性计算得到如下结论:采用本文方法计算得出滑塌式危岩W12、倾倒式危岩W22、坠落式危岩W4处于稳定状态,倾倒式危岩W59处于欠稳定状态,而坠落式危岩W16处于基本稳定状态,这与2008~2013年现场监测情况相吻合。针对滑塌式危岩W12讨论得出,最小主应力随时间呈线性降低趋势,随着最小主应力持续降低,泥岩先产生剪切破坏后发生拉破坏,泥岩破裂角随着最小主应力降低不断降低;随着岩腔深度加大,稳定系数呈凹曲线降低趋势,岩腔后退是危岩失稳的根本原因;存在一临界岩腔深度,此岩腔深度下危岩将发生失稳破坏,并可据此预测危岩失稳时间。     

采矿巷道围岩变形机制数值模拟研究

研究了大冶铁矿龙洞-62 m、-74 m水平采矿巷道开挖后的二次应力分布及巷道变形机制。首先根据现场工程地质勘查和室内岩石力学试验对巷道围岩进行了工程地质分组和岩石力学参数确定;在此基础上运用FLAC3D数值模拟软件研究了巷道开挖后的应力应变状态,分析了围岩变形机制;并根据-74 m水平采矿巷道的收敛监测数据对比验证了数值模拟结果。研究结果表明,围岩条件不同的采矿巷道其二次应力分布影响范围有所差异,但围岩主应力总体上表现为由巷道边墙中下部位的压应力集中带逐步过渡到拱顶、底板一定范围内的拉应力集中带。     

法向卸荷下贯通裂隙岩样压剪破坏特征研究

开展贯通裂隙岩样在法向卸荷下的压剪破坏研究,有助于更好地理解岩体裂隙在卸荷条件下的失稳特征。对单裂隙贯通岩样开展剪切应力峰前一系列法向卸荷破坏试验,分析了法向卸荷下试样裂隙的变形及能量演化规律。结果表明,法向卸荷下裂隙压剪失稳时的法向应力大于常规直剪试验的对应值,裂隙的抗剪能力减弱;剪切位移随法向应力减小而增大;卸荷变形比K随裂隙轮廓面积比Rs的增加而增大;卸荷过程中总变形能U0由减小转变为增加的现象可用于预测岩石裂隙卸荷失稳。研究结果对理解贯通裂隙岩体压剪卸荷失稳破坏具有一定参考意义。     

不同应力路径下土石混填体变形力学特性大型三轴试验研究

为更加全面地掌握土石混填体的变形、力学特性,采用大型三轴剪切仪分别在常规路径、等应力比路径和等p路径下对不同含石量的土石混填体进行了一系列试验研究。试验结果表明,在低围压下土石混填体的应变软化特征并不明显,其黏聚力普遍较低,而内摩擦角则较高,随试样的含石量从25%增加到70%,其内摩擦角和变形模量亦会随之增加。土石混填体在常规路径和等p路径下的应力-应变关系曲线比较相似,其所求得的抗剪强度指标亦比较接近;而等应力比路径下的应力-应变关系曲线的形态和发展过程均与之存在明显差别。土石混填体在常规路径下具有低压剪胀和高压剪缩的特性,含石量越高,低压剪胀性越明显,反之,则高压剪缩性更为显著;等应力比路径下试样的体积应变与轴向应变之间近似呈线性关系,其斜率随应力比增加而减小;等p路径下试样则是以剪胀性为主,尤其在平均主应力较低时更为显著。     

A nonlinear deformation modulus based on rock mass classification

Summary This paper describes the development of an empirical nonlinear stress dependent expression for the deformation modulus of rock masses based on rock mass classification. The expression defines the deformation modulus as the ratio of the deviator stress at failure to the major principal strain at failure. The Hoek and Brown failure criterion is used to predict the deviator stress at failure. Research was directed toward developing a failure criterion defining the major principal strain at failure. The expression for the deformation modulus was extended to rock mass conditions through correlations with observed deformations from case history studies and predicted deformations from finite element analyses.     

冻融循环下花岗岩剪切蠕变试验与模型研究

为探讨冻融对寒区工程岩石剪切蠕变特性的影响,以吉林省辉白隧道花岗岩为研究对象,对经历不同冻融循环次数的试样开展细观特征分析和剪切蠕变试验。试验结果表明：（1）随着冻融次数的增加,试样裂隙、孔隙不断扩展,岩石表面损伤现象愈发明显;（2）试样主要以中小孔隙为主,孔隙度随着冻融次数的增加呈非线性增长趋势;（3）随着冻融循环次数增加,蠕变变形量和蠕变速率逐渐增大,而蠕变时长、破坏应力和长期强度均呈现明显降低趋势。根据试验结果,进行冻融岩石非定常蠕变参数的表达,提出了冻融岩石损伤黏性元件,构建了花岗岩冻融剪切蠕变本构模型。将蠕变试验曲线和理论模型拟合曲线进行对比,验证了模型的正确性和适用性。通过对蠕变参数进行敏感性分析,研究了其对花岗岩蠕变变形的影响,并给出了蠕变参数随冻融循环次数的变化规律。该研究结果对于寒区岩体工程长期稳定性评价具有指导意义。