THE EFFECTS OF THREE-DIMENSIONAL PENNY-SHAPED CRACKS ON ZONAL DISINTEGRATION OF THE SURROUNDING ROCK MASSES AROUND A DEEP CIRCULAR TUNNEL

In this study, it was assumed that three-dimensional penny-shaped cracks existed in deep rock masses. A new non-Euclidean model was established, in which the effects of pennyshaped cracks and axial in-situ stress on zonal disintegration of deep rock masses were taken into account. Based on the non-Euclidean model, the stress intensity factors at tips of the pennyshaped cracks were determined. The strain energy density factor was applied to investigate the occurrence of fractured zones. It was observed from the numerical results that the magnitude and location of fractured zones were sensitive to micro- and macro-mechanical parameters, as well as the value of in-situ stress. The numerical results were in good agreement with the experimental data.     

THE ZONAL DISINTEGRATION MECHANISM OF SURROUNDING ROCK AROUND DEEP SPHERICAL TUNNELS UNDER HYDROSTATIC PRESSURE CONDITION: A NON-EUCLIDEAN CONTINUUM DAMAGE MODEL

A new non-Euclidean continuum damage model is proposed to investigate the zonal disintegration phenomenon of the surrounding rocks around deep spherical tunnels under hydrostatic pressure condition as well as the total elastic stress field distributions.The elastic stress fields of the surrounding rocks around deep spherical tunnels under hydrostatic pressure condition are obtained.If the elastic stresses of the surrounding rocks satisfy the strength criterion of the deep rock masses,the number,size and location of fractured and nonfractured zones are determined.The efect of physico-mechanical parameters of the surrounding rocks on the zonal disintegration phenomenon is studied and numerical computation is carried out.It is found from numerical results that the number,size and location of fractured and non-fractured zones are sensitive to the physico-mechanical parameters of the surrounding rocks.     

深部岩体损伤对分区破裂化效应的影响

深部岩体是一种具有初始损伤的非连续介质,其内部存在着大量节理和裂隙。深部岩体分区破裂化效应与节理、裂隙的扩展、连接和汇合密切相关。由于节理、裂隙的扩展、连接和汇合可能会导致深部岩体的内部空间由欧氏几何空间向非欧几何空间转化。本文利用损伤变量表示含节理和裂隙岩体的损伤,根据非欧几何模型,获得了静水压力和非静水压力情况下深部圆形洞室损伤围岩的应力场,确定了损伤变量对圆形洞室围岩的应力场和分区破裂化效应的影响。采用Hoek-Brown准则,获得了静水压力和非静水压力情况下损伤围岩的破裂区与非破裂区的分布规律。本文的研究结果为理解深部岩体分区破裂模式提供了有价值的参考。     

ZONAL DISINTEGRATION MECHANISM OF DEEP CRACK-WEAKENED ROCK MASSES UNDER DYNAMIC UNLOADING

Size and quantity of fractured zone and non-fractured zone are controlled by cracks contained in deep rock masses. Zonal disintegration mechanism is strongly dependent on the interaction among cracks. The strong interaction among cracks is investigated using stress superposition principle and the Chebyshev polynomials expansion of the pseudo-traction. It is found from numerical results that crack nucleation, growth and coalescence lead to failure of deep crack- weakened rock masses. The stress redistribution around the surrounding rock mass induced by unloading excavation is studied. The effect of the excavation time on nucleation, growth, interaction and coalescence of cracks was analyzed. Moreover, the influence of the excavation time on the size and quantity of fractured zone and non-fractured zone was given. When the excavation time is short, zonal disintegration phenomenon may occur in deep rock masses. It is shown from numerical results that the size and quantity of fractured zone increase with decreasing excavation time, and the size and quantity of fractured zone increase with the increasing value of in-situ geostress.     

开挖卸荷条件下深部圆形洞室各向同性围岩分区破裂化机理

论文提出了一种静水压力条件下动态开挖卸荷对深埋圆形洞室各向同性围岩分区破裂化影响的力学模型,获得了动态开挖卸荷条件下深埋圆形洞室围岩的弹性应力场,确定了卸荷速率和岩体动态力学参数对深埋洞室各向同性围岩分区破裂化现象的影响.利用动态Hoek-Brown准则,获得了围岩破裂区和非破裂区的数量和宽度.数值计算结果表明卸荷速率和岩体动态力学参数对破裂区与非破裂区的数量和宽度有较大影响.     

Rock burst of deep circular tunnels surrounded by weakened rock mass with cracks

Rock masses containing pre-existing cracks are considered as non-homogeneous geomaterials. During excavation of tunnels, pre-existing cracks may nucleate, grow and propagate through rock matrix, then secondary cracks may appear. The stress concentration at the tips of secondary cracks is comparatively large, which may lead to the unstable growth and coalescence of secondary cracks, and consequently the occurrence of fractured zones. For brittle rocks, the dissipative energy of slip and growth of pre-existing cracks and secondary crack growth is small, but the elastic strain energy storing in rock masses may be larger than the dissipative energy of slip and growth of pre-existing macrocracks and secondary crack growth. The sudden release of the residual elastic strain energy may lead to rock burst in crack-weakened rock masses. Based on this understanding, the criteria of rock burst in crack-weakened rock masses are established. The influences of the in situ stresses, micromechanical parameters and physico-mechanical parameters on the distribution of rock burst zones and area of rock burst zones are investigated in detail.     

Hoek–Brown criterion applied to circular tunnel using elastoplasticity and in situ axial stress

Based on the nonlinear Hoek–Brown failure criterion, elastoplastic analytical solutions are developed for the elastoplastic stresses, strains and plastic zones around a circular tunnel subjected to different value of the axial in situ stress. Effects of the transverse in situ stress, the axial in situ stress and the strength parameters of rock masses on the elastoplastic stresses, strains and plastic zones in the surrounding rock masses are investigated. It is found from the numerical results that the stresses, strains, and plastic zones in the surrounding rock depend not only on the transverse in situ stress but also on the axial in situ stress as well as the mechanical parameters of rock masses.     

Zonal fracturing mechanism in deep crack-weakened rock masses

The mechanical behaviors of deep crack-weakened rock masses are different from those of shallow crack-weakened rock masses. The surrounding rock in shallow crack-weakened rock mass engineering is classified into loose zone, plastic zone and elastic zone, while the surrounding rock in deep crack-weakened rock mass engineering is classified into fractured zone and non-fractured zone, which occur alternatively. It is assumed that the deep rock masses contain one joint set, in which the probability density function describing the distribution of sizes is assumed to follow the Rayleigh distribution, and the probability density function describing the distribution of spacing is assumed to follow the Weibull distribution. On the basis of strength criterion of deep rock mass, the near-field stress redistribution around circular opening induced by excavation is determined. The strong interaction among cracks is investigated by using the dislocation model. The nucleation, growth, interaction and coalescence of cracks were analyzed based on the strain energy density factor theory. When cracks coalesce, failure of deep crack-weakened rock masses occurs, fractured zone is formed. Then, size and quantity of fractured zone and non-fractured zone are given out. The size and quantity of fractured zone increase with decreasing strength of rock mass. The size and quantity of fractured zone increase with increasing in situ stress. Zonal fracturing phenomenon occurs once value of in situ stress is larger than the unaxial compressive strength of rock masses. The size and quantity of fractured zone decrease with increasing λ when p₂ > p₁. The size and quantity of fractured zone increase with increasing λ when p₂ < p₁.     

Zonal disintegration of surrounding rock mass around the diversion tunnels in Jinping II Hydropower Station, Southwestern China

By means of numerical simulation, the special phenomenon of zonal disintegration of surrounding rock mass around the diversion tunnels of Jinping II Hydropower Station is analyzed in this paper. In order to model the growth and coalescence of cracks within rock mass in Jinping II Hydropower Station, the weak-element is adopted. When cracks coalesce, failure of deep crack-weakened rock masses occurs and fractured zone is formed. The present result is different from the one obtained by the traditional elasto-plastic theory. The numerical results show that the slip-line zonal fracture is created within rock mass around the diversion tunnels in Jinping II Hydropower Station. Meanwhile, the magnitude and distributions of fractured zones are determined by numerical simulation. It is shown that the present results are in good agreement with the one observed by model tests. Through sensitivity analysis, the effects of stress condition, cohesion and the angle of internal friction on the phenomenon of zonal disintegration is determined.     

2021-07期目录

在陈宗基院士关于地下硐室长期稳定性力学问题研究的基础上,采用Sadovsky院士关于复杂地质岩体的等级构造学说,围绕深部岩体非均匀构造与封闭应力固有的统计力学属性,研究了岩体非均匀变形与封闭应力特性,以及深部硐室围岩的长期稳定性等两个力学问题。给出了岩体非均匀构造与封闭应力的数学表征;根据质量守恒定律得到了计算深部硐室围岩长期变形的一般公式;得到了围岩变形中劈裂扩容变形占主要部分的结论,并且阐明了深部围岩卸荷时更易出现劈裂破坏的原因。给出了劈裂破坏形态的演进序列与扩容位移的计算方法。将围岩松动圈范围、破裂区位置和边壁位移的计算结果与锦屏一级电站厂房现有监测数据进行了对比,两者相当吻合。     

深部岩体变形破坏的特征能量因子与应用

深部岩体在高地应力作用下储存了大量的弹性应变能。在开挖或爆破扰动作用下,原有的平衡状态被打破,围岩中形成了有势场和不平衡应力场。在不平衡力场和扰动场的共同作用下,岩体的变形与破坏表现出了诸如分区破裂化、大变形、岩爆以及人工地震等非线性行为。传统的连续介质理论并不能考虑岩体的构造特性与含能特性,因此无法很好地解释深部岩体的特殊非线性力学现象。特征能量因子从能量的角度出发,结合统计物理学观点,为分析深部岩体在动静荷载组合作用下的变形和破坏过程提供了有效的理论支撑。本文主要对特征能量因子进行了简要介绍,并回顾了其在深部岩体分区破裂以及动力诱发围岩不可逆变形等非线性工程灾害现象中的应用。     

Localization of deformation and stress–strain relation for mesoscopic heterogeneous brittle rock materials under unloading

Stress redistribution induced by excavation of underground engineering and slope engineering results in the unloading zone in parts of surrounding rock masses. The mechanical behaviors of crack-weakened rock masses under unloading are different from those of crack-weakened rock masses under loading. A micromechanics-based model has been proposed for brittle rock material undergoing irreversible changes of their microscopic structures due to microcrack growth when axial stress is held constant while lateral confinement is reduced. The basic idea of the present model is to classify the constitution relation of rock material into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The closed-form explicit expression for the complete stress–strain relation of rock materials containing cracks under unloading is obtained. The results show that the complete stress–strain relation and the strength of rock materials under unloading depend on the crack spacing, the fracture toughness of rock materials, orientation of the cracks, the crack half-length and the crack density parameter.     

岩体变形参数渐变取值模型及工程应用

分析了当前岩体现场试验、变形参数取值和常规计算模式存在的问题。针对河谷、斜坡岩体从表部到深部 ,岩体风化、岩体波速和岩体地应力等彼此均呈现渐变性特征 ,提出了岩体变形参数渐变取值模型 ,并以黄河上游李家峡水电站初期蓄水监测资料 ,论证了岩体变形参数空间上连续性渐变取值的合理性和科学性 ,对今后其它岩体工程的力学参数取值和计算模型的概化提供了指导意义.     

基于岩石蠕变及D-P准则的深部巷道围岩弹塑性分析

为深入分析岩石蠕变及中间主应力对深部高应力软岩巷道围岩稳定性及分区变形特性的影响机理, 通过剖析长期荷载作用下深部巷道围岩的变形特征, 揭示了基于围岩蠕变特性的深部巷道围岩四分区变形机理, 并依据Druck-Prager准则及关联性流动法则考虑了中间主应力、岩体扩容及应变软化的影响, 推导出巷道围岩各分区应力、变形和半径解析解. 结合工程算例, 通过对现场监测数据以及不同力学模型的计算结果进行对比分析, 论证了本文理论的科学性, 并揭示了不同强度准则和围岩参数对深部巷道围岩各分区形态的影响规律. 结果表明, 忽视深部巷道围岩的蠕变特性将导致初始黏聚力的取值大于实际值, 从而造成围岩的理论承载能力大于实际承载能力; 在[0, 0.7]区间内提高中间主应力系数可有效控制围岩塑性区及破碎区范围的扩展; 当初始黏聚力及内摩擦角减小时, 中间主应力系数对围岩塑性区、破碎区范围及巷道变形的影响力显著提高. 研究成果可为地下工程支护设计及围岩稳定性评估提供理论参考.      

裂隙岩体水力特征数值模拟试验的初步结果

本文利用文献中的数据进行数值模拟试验,结果表明：流量等效和渗透作用力变形效果等效互不包容;不同级别结构面体系难以用连续介质方法处理,同级别结构面体系在工程作用力变化范围内,可由其静态水力特征推演动态特征;工程岩体结构面间充分连通的假设过于粗糙,常规试验手段获取的量值倾向于局部特征的反映。     

Upper and lower bounds for constitutive relation of crack-weakened rock masses under dynamic compressive loads

A micro-mechanics-based model is proposed to investigate the rate-dependent constitutive relation for crack-weakened rock masses subjected to dynamic compressive loads. The present micro-mechanical model reveals that the nucleation, growth and coalescence of sliding cracks dominate the failure and macroscopic properties of crack-weakened rock masses subjected to dynamic compressive loads. The interactions among multiple parallel sliding cracks in crack-weakened rock masses subjected to dynamic compressive loads are examined asymptotically in an explicit and quantitative manner in order to reveal fully their so-called shielding and magnification effects on the stress–strain relation. Based on the micro-mechanical framework and the asymptotic analysis, analytical upper and lower bounds are proposed for the rate-relation for rock masses containing multiple rows of echelon cracks subjected to dynamic compressive loads. The factors that affect the rate-dependent properties of crack-weakened rock masses have been analyzed. The strain energy density factor approach, which is related to crack growth velocity and dynamic fracture toughness of rock material, is employed in the analysis. The rate-dependent constitutive relation of crack-weakened rock masses is derived from micro-mechanical framework and the asymptotic analysis. The closed-form explicit expression for the rate-dependent constitutive relation of rock masses containing echelon cracks subjected to dynamic compressive loads is obtained. Finally, the present model is used to analyze the complete stress–strain relation and strength for jointed rock masses at shiplock slope of the Three Gorges Dam.     

深部节理岩体塑性损伤耦合微面模型

采用微面模型理论和损伤力学方法,建立了节理岩体的弹塑性损伤耦合微面模型. 在节理岩体的微面上,将岩体视为由节理面与岩石组成的二元介质,以节理连通率作为岩体沿该方向的面积损伤变量,考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制,基于塑性理论建立了节理岩体的微面塑性损伤增量本构关系. 采用微面物理量与宏观物理量的几何约束模型,根据微面方向积分导出了节理岩体的宏观弹塑性增量本构关系. 编制了节理岩体微面模型的MARC有限元子程序,对节理岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳定问题进行了数值模拟研究. 数值计算结果表明,该模型能很好地揭示载荷作用下节理岩体的各向异性非弹性变形和次生节理演化过程.      

深部岩体圆形巷道围岩松动圈形成机理数值试验研究

采用损伤力学和统计理论的单元本构模型,探讨深部岩体巷道围岩破坏规律, 利用岩石破裂过程分析软件RFPA$^{\rm 2D}$, 对深部 岩体中的圆形巷道的变形及非线性渐进破坏特征、巷道周边关键部位的位移和应力变化进行 了分析,研究了圆形巷道围岩松动圈的形成机理. 研究表明：巷道开挖后,其周边形成应力 集中带,随着围岩压力的持续作用,巷道周边产生塑性变形区,并沿径向形成裂纹;随着裂 纹的不断扩展,巷道周边出现松动破裂区,即松动圈. 产生破裂区后,应力集中程度减弱, 应力高峰点向远处转移.     

高地应力岩体多孔爆破破岩机制

深部岩体爆破破岩是爆炸荷载与高地应力共同作用的结果。基于一些简化假设,建立了一个高地应力岩体双孔爆破计算模型,采用光滑粒子流体力学-有限元方法耦合数值模拟方法,研究了高地应力作用下炮孔间裂纹的传播及贯通过程,分析了炮孔周围应力场动态演化过程与分布特征。研究结果表明:爆破引起的岩体开裂主要是环向动拉应力所致,地应力对岩体的压缩降低了炮孔周围环向动拉应力、缩短了环向动拉应力的作用时间,因而对爆炸致裂起抑制作用;静水地应力条件下多孔爆破时,垂直于炮孔连线方向传播的爆生裂纹更易受到地应力的抑制;对于高地应力岩体爆破,炮孔间的裂纹扩展长度随地应力水平的提高而减小,裂纹主要沿最大主应力方向扩展,因此沿最大主应力方向布置炮孔、缩短炮孔间距有利于炮孔间裂纹的连接贯通,形成良好的爆破开挖面。     

深埋地下洞室断裂型岩爆机理的数值模拟

针对深部岩体中由断层、节理等不连续性结构面引发的岩爆地质灾害, 根据深埋地下隧洞中潜在发震断裂的分布特征和几何形态建立数值分析模型, 采用离散元单元法模拟存在刚性平直断裂的深部围岩的开挖响应, 并分别考察开挖接近并通过断裂附近时围岩应力状态的变化特征。通过探讨断裂的存在对围岩应力状态改变的作用机理, 揭示出断裂型岩爆是开挖面附近一定范围内存在的断裂构造在高应力作用下发生错动, 导致能量突然释放, 对围岩造成强烈冲击作用的结果, 基本与地震的断层“粘滑”机制相类似。