Distinct element simulation of ultimate bearing capacity in jointed rock foundations

In this paper, distinct element method numerical modeling is applied to evaluate bearing capacity of strip footing rested on anisotropic discontinuous rock mass. As yet, a little work has been carried out to investigate the effect of joint set orientation on the bearing capacity of rock mass. Generally, the overall behavior of rock mass under loading is very complicated and such analysis should include deformation determination, sliding along discontinuities and failure of rock material. Failure mechanism of rock mass depended on both geometrical parameters of joint sets and strength parameters of rock mass. In this research, it is assumed that rock mass contains one joint set, and therefore the anisotropy in bearing capacity and rock behavior is only due to the existence and orientation of the joint set. In this study, by assuming constant strength parameters and using Mohr–Coulomb failure criterion for the single joint set and nonlinear Hoek–Brown failure criterion for rock material, variation of the bearing capacity values and the type of failure mechanism of rock mass with different joint set dips is investigated. The obtained results indicate that the ultimate bearing capacity of rock mass containing one joint set varies between 27 and 86 % of intact rock.     

A thermomechanical anisotropic continuum model for geological materials with multiple joint sets

Joints in geological materials introduce elastic compliance and weak planes on which sliding can occur. Although these materials can have multiple joint sets, they often have preferred orientations that cause both elastic and inelastic anisotropic response even when the unjointed material is isotropic. Azimuthal variations in radial velocity and polarity of tangential motion have been observed in experimental data for wave propagation caused by an initially spherical source in a geological material with multiple joint sets. This observed tangential ground motion was found to be related to mechanical anisotropy caused by preferred orientations of joints in the rock. This paper describes thermomechanical continuum constitutive equations, which model the effects of multiple persistent joint sets. A number of quasi‐static examples are considered, which show that the proposed model predicts anisotropic effects of sliding on multiple joint sets similar to those exhibited by computationally expensive mesoscale calculations, which model joint sets explicitly.     

Slope failure analysis considering anisotropic characteristics of foliated rock masses

Weak planes affect the strength and deformational behaviors of rock slopes, and the anisotropic characteristics of rock mass should be considered in slope stability analysis. Effects of joint plane orientations on failure mechanism and strength response of inherently anisotropic rock samples were firstly investigated. The specimens with various orientations of joints were evaluated under uniaxial compression, Brazilian tensile, and direct shear tests. By treating the foliated rock as transversely isotropic materials, the relevant elastic constants and strength parameters were obtained from experimental results. The slope damage zone was then investigated using Comsol Multiphysics code based on Hoffman criterion. It is indicated that the failure mechanism and strength response depend highly on the inclination of specimens with respect to the loading direction. For disks with the same inclination angle, the value of tensile strength has an increasing trend with the total fracture length. Numerical results show that partial slope mass failed in single slope and no large-scale landslide occurred. The failure pattern in numerical results agrees well with the field observations. The cooperation between the experimental results and the numerical results allows an in-depth analysis of the experimental results and thus better understanding the dominant effect of joints on the deformation and failure of rock mass.     

非均质共面断续节理岩体拉伸剪切破裂机制研究

工程开挖面附近卸荷扰动区的岩体,受结构面和拉应力共同影响作用,其变形和破坏具有拉剪复合特征。为研究节理岩体的拉剪力学特性,基于颗粒离散元法针对共面断续节理岩体开展了系列数值模拟研究。通过假设粒间接触的力学参数服从Weibull分布表征岩体的非均质性,探讨了非均质性、均质度、法向拉应力和节理连通率对节理岩体拉剪强度和破坏模式的影响。研究表明:拉剪应力条件下非均质性节理岩体主要沿阶梯型破裂面破坏,剪应力-水平位移曲线可以分为线性变形阶段、非线性变形阶段、峰值及峰后阶段;随均质度提高,节理岩体的剪切强度逐渐增加且提升幅度逐渐减弱,趋于均质岩体,岩体中微裂纹由弥散型分布向破裂面集中;节理岩体峰值剪切强度和法向拉应力的大小呈非线性负相关关系;岩体剪切强度随节理连通率增加而显著降低。     

采动影响下巷道围岩变形与破坏的弱结构面效应

岩体中弱结构面的产状要素和力学参数是影响深部巷道围岩变形和破坏的重要因素。通过对金川矿区地下开采过程中发生的与弱结构面效应有关的巷道变形特征、破坏模式与机理分析得出,在断层附近或穿过断层开挖时,往往会引起断层活化,引发巷道变形和破坏。当在节理岩体中开挖时会使围岩中低强度的节理裂隙等弱结构面正应力减小,抗剪强度降低,从而发生相对滑移、变形。此外,与巷道开挖相比,矿体开采的影响不仅在持续时间上,还是在影响以及程度上都是造成巷道变形和破坏的弱结构面效应更加显著的主要原因。     

An anisotropic strength criterion for jointed rock masses and its application in wellbore stability analyses

In this paper, an anisotropic strength criterion is established for jointed rock masses. An orientation distribution function (ODF) of joint connectivity, is introduced to characterize the anisotropic strength of jointed rock masses related to directional distributed joint sets. Coulomb failure condition is formulated for each plane of jointed rock masses by joint connectivity, where the friction coefficient and cohesion of the jointed rock mass are related to those of the intact rock and joint and become orientation dependent. When approximating joint connectivity by its second‐order fabric tensor, an anisotropic strength criterion is derived through an approximate analytical solution to the critical plane problem. To demonstrate the effects of joint distribution on the anisotropic strength of jointed rock masses, the failure envelopes are worked out for different relative orientations of material anisotropy and principal stress axes. The anisotropic strength criterion is also applied to wellbore stability analyses. It is shown that a borehole drilled in the direction of the maximum principal in situ stress is not always the safest due to the anisotropic strength of the jointed rock mass. Copyright © 2007 John Wiley & Sons, Ltd.     

深埋顺层偏压隧道围岩破坏机理及规律研究——以郑万线某隧道为例

处于薄—中层倾斜层状岩体中的深埋隧道常会产生地质顺层偏压的问题,导致隧道局部塌方、偏压变形及支护结构破坏。本文以郑万线某隧道为例,采用理论分析、数值模拟方法对深埋顺层隧道的破坏机理及不同结构面参数下的破坏规律展开了研究。研究结果表明:（1）深埋顺层偏压隧道洞周围岩将根据其切向应力与结构面夹角的不同发生岩层拉裂破坏、结构面剪切破坏及岩体自身破坏,其中切向应力与结构面平行位置,即反倾侧拱腰及顺倾侧拱脚位置主要发生拉裂破坏,此处围岩塑性区范围最广,围岩位移最大,围岩处于极不稳定状态;（2）顺层偏压隧道的破坏规律与结构面强度参数有直接关系,围岩塑性区范围及围岩位移均随着结构面摩擦角的增大而降低,且降低趋势逐渐放缓,当结构面摩擦角达到岩体摩擦角后,结构面摩擦角继续增加对围岩稳定性影响较小;（3）围岩塑性区及围岩位移场偏压分布特征随结构面倾角的变化而整体旋转,且对于隧道底部而言,结构面最不利倾角为0°,此时隧底最大上鼓量大于其他倾角下的最大上鼓量;对于隧道拱部而言,最不利倾角为40°,此时洞周最大收敛值大于其他倾角下的最大收敛值,最不利位置位于反倾侧拱腰。     

节理岩体锚杆的综合变形分析

在总结国内外对节理岩体中锚杆加固机制的试验研究和理论探讨基础上,综合考虑锚杆的切向和轴向变形能力,建立节理锚固锚杆在剪切荷载作用下的变形模型,将节理锚固锚杆的变形区划分为弹性变形段和挤压破坏段,引入表征挤压破坏段长度的变量,对锚杆与岩体的相互作用机制进行理论分析,推导了剪切荷载与剪切位移和轴向荷载与轴向位移的关系。通过分析锚杆的屈服破坏形式,得到了确定挤压破坏段长度的方法。最后,通过算例分析了挤压破坏段长度与锚杆直径、岩体强度、锚固角度等参数的关系,得到了以下结论：（1）节理锚固锚杆抗剪作用的实质是锚杆调动岩体的抗压强度抵抗节理切向荷载。在抗压强度较高的硬岩中,挤压破坏段局限于节理面附近,锚杆影响范围小;而在抗压强度较低的软岩中,挤压破坏段较大,而且会产生较大的剪切变形,锚杆影响范围较大。（2）锚杆屈服破坏形式与岩质和锚杆直径有关。硬质岩体发生剪切屈服,而较软岩体中容易发生弯曲屈服;小直径锚杆一般直接剪切屈服,而大直径锚杆可能发生弯曲屈服。锚杆屈服破坏后出现塑性铰,挤压破坏段范围在节理一侧约为直径的1～2倍,继续增加剪切荷载,挤压破坏段长度不再增大。（3）随岩质的不同,锚杆锚固节理的最优锚固角变化较大。岩质较硬时,最优锚固角度较小,反之则较大。     

深埋隧道层状围岩变形特征分析

层状岩体在地下工程中经常遇到的,它具有明显的各向异性力学性质。结合共和隧道现场监测和数值模拟相结合的方法对层状岩体的破坏特征进行了分析,研究结果表明：围岩的变形位移、破坏区都主要集中隧道拱顶右侧,即靠河侧大于靠山侧。隧道围岩变形破坏区不在最大主应力方向上,而是在岩体层理垂直方向。层状岩体中洞室变形破坏特征除了因地形产生的偏压影响外,更重要的是受地层结构特征的影响,即与层状岩体的力学性质极大相关,其结果可为指导隧道的施工和设计提供有效依据。     

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

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

基于3D雕刻技术的岩体结构面剪切各向异性研究

相同形貌结构面重复性剪切试验和各向异性剪切试验一直是岩体结构面剪切力学特性试验研究中的难点,而该问题对于工程岩体的开挖力学响应和稳定性分析、评价与控制至关重要,问题的关键在于同一种岩石的结构面形貌的再现。为此,基于3D扫描和3D雕刻技术,重复制备了3种不同粗糙度的大理岩结构面试样,开展了不同法向应力下相同形貌结构面的各向异性剪切试验。试验结果发现：（1）同一法向应力、不同剪切方向下,相同形貌结构面的剪切强度、剪胀和剪切破坏特征均呈现明显各向异性;（2）结构面粗糙度各向异性很大程度上决定了剪切强度的各向异性,两者具有较好的正相关性;（3）随着法向应力的升高,结构面剪切特征的各向异性有逐渐弱化的趋势。同时,研究还充分表明,3D雕刻技术是系统开展结构面剪切力学特征各向异性研究的可靠手段,可在将来研究中发挥更为重要的作用。     

高应力区隧道围岩变形破坏的数值模拟及物理模拟研究

隧道失稳和维护困难是高地应力隧道的普遍问题,对隧道的支护设计提出了更高的要求。研究从地下工程岩体应力环境变化和岩体强度变化的角度探讨了高应力隧道围岩的变形破坏机制。根据重庆某深埋隧道围岩实际情况,运用FLAC3D三维显式有限差分法分析软件,建立了摩尔-库仑剪破坏与拉破坏复合的应变软化模型。通过隧道的三维数值计算,分析了高应力环境下隧道周边塑性区分布、应力场、位移场等的分布特点,得到了高应力隧道围岩在高地应力环境下的破坏规律。通过物理模型验证了高应力隧道围岩的破坏特点,并进行了超载试验,将其与数值模拟进行对比,进一步验证了所建数值模型的科学性。     

岩溶区隧道围岩--支护体系稳定性研究

岩溶是地下水地质作用的过程和结果，是隧道建设面临的不良地质现象之一。正在建设中的某高速公路隧道出口段岩溶地质作用强烈，溶隙、孔洞发育，且受构造变形及后期淋滤溶蚀作用强烈，力学强度低，岩性复杂。其角砾状粘土岩和水云母粘土岩多已泥化成粘土，遇水易膨胀。岩溶区隧道围岩——支护体系的稳定性是该隧道建设面临的主要问题。文章结合隧道施工过程中的围岩变形监控量测和三维弹塑性有限元数值仿真模拟对隧道围岩——支护体系的稳定性进行了研究。围岩水平收敛和拱顶下沉位移及速率变化特性表明，在监测时段内隧道围岩无趋于稳定的迹象。施工动态力学数值仿真模拟表明，在初期支护条件下，围岩出现了较大范围的塑性破坏区，两隧道之间(间距45．0m)围岩破坏接近度达到0．70，最大竖向位移达到3．0mm，表明隧道全断面施工时，两隧道相互作用影响程度较大。为确保隧道围岩——支护体系的稳定，选择合理的施工方法，制定切实可行的支护方案提供了信息。     

Numerical study of excavation induced fractures using an extended rigid block spring method

This paper presents a numerical study of fracturing process induced by excavation around a gallery using an extended rigid block spring method (RBSM). The surrounding rock mass is characterized by an assembly of rigid blocks based on a degraded Voronoi diagram. The macroscopic mechanical behavior of rock is related to that of interfaces between blocks. The mechanical behavior of each interface is described by its elastic stiffness and failure criterion. The failure process of interface is controlled by both normal stress and shear stress. Both tensile and shear failures are considered. The macroscopic fracturing process is described by the coalescence of cracked interfaces. The rock structural anisotropy is taken into account through a spatial variation of elastic stiffness and failure strength of interfaces. A series of sensitivity studies are performed to investigate effects of gallery orientation, failure strength of interfaces and rock structural anisotropy on gallery deformation and fracturing. Numerical results are compared with in-situ observations in terms of fracture patterns.     

泥化弱胶结软岩地层中矩形巷道的变形破坏过程分析

在西部地区,一定数量的矿区处于泥化弱胶结软岩地层,此类软岩胶结性差、强度低、遇水泥化。矩形是采区巷道的常用型式,但其断面受力不均、稳定性差。在上述软岩地层中的矩形巷道承载力低、变形量大、变形持续时间长,给煤矿的安全生产带来极大困难。以内蒙古新上海一号煤矿皮带顺槽矩形巷道为背景,运用FLAC3D软件中的Cvisc黏弹塑性模型,对矩形巷道的变形破坏进行了数值模拟,并将模拟结果与现场监测结果对比分析。结果表明:巷道开挖支护后,受断面形状影响,矩形巷道四角出现压应力集中和顶板受拉区,巷道顶板下沉量大,底板底臌严重,两帮向巷道挤出;受围岩岩性影响,围岩进入塑性的时间快短、范围大,塑性区超出了支护体的作用范围,造成锚杆(索)的锚固效果难以完全发挥,围岩出现整体滑动的现象;巷道变形呈现出流变变形的特性,变形量随时间持续增加,持续的蠕变变形超出了支护体的可控范围,最终引起巷道的失稳破坏。     

昆明白泥井3号隧道围岩稳定性FLAC3D模拟

应用FLAC^3D数值模拟软件对昆明白泥井3号隧道围岩稳定性进行了分析和评价。结果表明隧道支护后围岩的最大变形出现在隧道右侧拱顶，为1．00cm；隧道支护后左侧拱和右侧拱脚处剪应变明显增大，为应力集中区；围岩可以通过自身应力、应变和能量的调整，使整个隧道围岩处于稳定状态；同时也表明FLAC^3D用于隧道围岩变形分析是可行的。     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.     

含二阶起伏体的模拟岩体节理试样剪切特性试验研究

利用规则的小尺寸锯齿构造二阶起伏体,对不同二阶起伏体高度的人工节理进行常法向应力下的直剪试验,研究了节理剪切力学特性。对节理的剪切强度、变形特征进行了分析,并对破坏特征做出解释。试验结果表明：二阶起伏体对节理剪切力学特性有重要影响,含二阶起伏体节理的剪切力学特性与只含有一阶起伏体节理的剪切力学特性不同,只含一阶起伏体的节理剪应力只有一个峰值,含二阶起伏体的节理会出现波浪状的剪胀曲线和多峰值剪应力,且峰值剪应力随剪切位移增大依次出现并逐个减小;随着二阶起伏体高度增大,节理峰值剪切强度增大,节理依次出现磨损破坏、多次性剪断破坏、一次性剪断破坏。     

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

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