On the role of joint roughness on the micromechanics of rock fracturing process: a numerical study

The deformation process and failure mechanism of rock mass with increased joint roughness subjected to unconfined compression are investigated in this study using discrete element method. A numerical model is developed using soft-bonded particle and validated to realistically replicate the mechanical response of the rock mass. The micro-parameters of the rock material are first determined, and the effects of the joint roughness on the macromechanical response and fracture growth mechanism are then investigated. Analyses are also performed to examine the tensile and shear crack distributions, acoustic emission (AE) characteristics, coordination number, and crack anisotropy to advance the current understanding of the role of joint roughness on the mechanical behavior and deformability of rock mass. The results show that strength and deformability of the jointed rocks are highly dependent on the joint orientation and roughness. Joint roughness is found to restrain the propagation and coalescence of microcracks and AE events from the interlocking of asperities. In addition, the spatial distribution of the contact forces allows for better understanding of the effect of joint inclination angle on the response of the investigated rock samples.

     

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.     

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

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

白鹤滩柱状节理岩体真三轴模型试验研究

在柱状节理玄武岩工程地质调查、现场精细结构描述等工作基础上,开展了柱状节理玄武岩体的脆性材料模型真三轴物理模拟试验。研究了复合柱状岩体在不同卸荷-加载应力途径下的破裂、破坏机制及不同方向三轴主应力作用下的应力-应变关系及其强度特性;开展的8种不同应力比和6种不同受力方向下各向异性试验结果表明,柱状节理玄武岩体强度、变形特征受三向不等围压和柱状节理特性影响明显,强度变形的各向异性特征显著。试验结果表明,应力比 对试样的强度、变形和破坏模式具有控制作用,随着应力比的变化,试样的极限强度呈现分段特性;较小应力比下试样的变形和破坏模式受应力比和原生节理的共同影响;较大的应力比对试样的变形和破坏模式的影响比原生节理的影响显著;随着应力比的增大,中主应力由约束侧向变形变为主动加载,这一作用的变化影响了极限荷载和破坏模式。     

On the stress-force-fabric equation in triaxial compressions: Some insights into the triaxial strength

The strength of granular materials during triaxial compression is investigated via a grain scale analysis in this paper. A 3D Discrete Element Method (DEM) program provides the triaxial strength data and helps to validate the micromechanical analysis. Some standard methods in statistics are employed first to quantitatively examine the assumptions made when deriving the stress-force-fabric (SFF) equation. After careful validation, a more concise format for the SFF equation is proposed for triaxial compressions. With this SFF equation, the strength is found to be jointly contributed by the magnitudes of the contact force anisotropy and fabric anisotropy. The influence of the initial void ratio, confining pressure and loading direction on the development of contact force anisotropy and fabric anisotropy is examined and presented. With similar techniques, the “force” term in the SFF equation is further decoupled, and an equation is obtained such that it explicitly links the contact force term with the friction coefficient between grains, a tensor defined as a statistic of the normal contact forces and a tensor defined as a statistic of the mobilisation status of contacts. Based on this equation, another equation regarding the stress ratio of granular assembly is obtained, and it clearly indicates two sources that contribute to the phenomenological friction nature of granular assembly. These two sources are caused by the contact force at the grain scale. The first is the anisotropy of the average normal contact forces, and the second is the mobilisation of contacts.     

基于剪切试验的预应力锚杆变形性能分析

为了研究锚杆对节理剪切性能的作用机制和模式,开展了锚固节理岩体的实验室剪切试验,模拟了不同强度的岩体在剪力作用下的变形和受力特征,对比了锚杆加固前、后岩体的剪切变形规律,分析了节理岩体强度、预应力及锚固方式对节理的抗剪能力的影响,试验过程中通过应变片测点监测锚杆轴力的变化规律,研究了节理剪切过程中锚杆的轴向受力机制和变形特性。剪切试验完成后,取出剪切变形后的锚杆,统计了变形段长度与各参数之间的关系。通过试验结果分析,剪力-位移曲线存在明显的3个阶段：弹性阶段、屈服阶段和塑性强化阶段,曲线近似呈现双直线特征,弹性段锚杆主要发挥销钉作用;屈服段锚杆的轴向作用开始调动,锚杆同时存在销钉和约束作用;塑性段锚杆不再发挥销钉作用,只是依靠其轴向约束作用限制岩体的变形;曲线表现出韧性增强的剪切性能,这就使得锚杆锚固节理岩体的破坏特性由脆性转变为塑性,从而提高了岩体的稳定性和安全度。节理面的滑动对锚杆产生剪切作用,由于切向变形而产生了附加的锚杆轴向变形,锚杆的轴向应力随着剪切位移的增大呈增长的趋势,变形剧烈区域集中于节理面附近,且距离节理面越近其轴向应力增大越多。     

含预制节理岩体卸荷条件下力学特性试验研究

节理对卸荷条件下岩体的力学性质有重要影响。通过含2条不同间距预制断续节理岩体的三轴卸荷破坏试验,研究了节理岩体在卸荷应力条件下的应力-应变特征、强度、变形特征、破坏规律及节理间距对岩体力学性质的影响。研究表明：相比于完整岩体,节理岩体卸荷破坏时从峰值强度跌落至残余强度过程中轴向应变较大,为完整岩体的3～4倍,岩体破坏时极限强度明显低于完整岩体,脆性特征不如完整岩体明显;节理岩体卸荷破坏时,变形模量有较大幅度的降低,其降低程度是同条件下完整岩体的6～7倍,节理间距越大,变形模量降低程度越大;与含预制节理岩样三轴加载试验结果相比,节理岩体卸荷条件下破坏程度更为强烈,除剪切破裂面外,沿最大主应力方向分布的不同级别的张性裂隙非常发育,预制节理的间距对岩体破坏形态影响不大。     

节理岩体超声测试及单轴压缩试验研究

为了获取节理岩体各向异性力学特性和声波传播规律的相互关系,通过制作圆柱形单一预制贯通节理试件,开展了7种不同倾角节理试件的超声波波速测试和单轴压缩试验。试验结果表明：（1）节理试件波速测试值离散性较大,但总体服从正态分布,波速平均值介于石膏试件和完整试件之间,波速随节理倾角增大呈线性递减;（2）不同倾角的节理面对节理试件力学特性的影响很大,在应力-应变关系、变形特征、强度特征和破坏模式上均表现出显著的各向异性特征;（3）节理试件的力学性质和超声波传播特性变化规律区别很大,采用超声波波速确定节理试件的力学参数会存在很大误差。     

类节理岩石直剪试验力学特性的数值模拟研究

节理岩体的剪切特性是主导岩体工程稳定性的关键因素。基于PFC2D离散元颗粒流程序,结合室内试验结果对比分析,选取合理的细观参数进行数值模拟,分别从细观角度研究了节理岩石的裂纹发展、能量转化及声发射现象等特性,从宏观角度研究了节理岩石的强度模型和破坏形态。结果表明：节理岩体主要呈现磨损和剪断两种破坏形态,不同的破坏形态对应不同的强度模型;随着剪切变形增加,岩体沿节理面发生破坏,弹性阶段以法向裂纹为主,而塑性阶段切向裂纹起主导作用,滑移区R、P裂纹贯通形成破碎带,节理面产生较大滑移;在应力达到峰值强度前,边界能主要转化为应变能,法向裂纹生成较多;越过峰值强度后,摩擦能快速增长,并伴随大量切向裂纹产生。与室内试验结果相比,PFC2D较好地模拟了节理岩体剪切力学特性,弥补了室内试验中无法进行细观特性研究的缺陷,对于节理岩体后期研究提供了一些参考。     

A micromechanics-based model for estimating localized failure with effects of fabric anisotropy

This paper presents a micromechanics-based approach to investigate the effects of fabric anisotropy on the behavior of localized failure in granular materials. Based on a micromechanical analysis, the origin of deviatoric stress is decomposed into two components: contact force anisotropy and fabric anisotropy. Using a micro–macro approach, the back stress is interpreted as an contribution to the change of the fabric’s principal direction. The evolution of the back stress is deduced from the stress–fabric relationship and determined with reference to the deviation of the principal directions between the rate of the reduced stress tensor and the actual reduced stress tensor. With this micro–macro framework, a mixed (isotropic–kinematic) hardening model is developed based on the classical isotropic hardening theory. A laboratory simple shear test is first analyzed to validate the proposed model and illustrate the kinematic-hardening mechanism which is usually displayed under non-proportional loading. The analysis further focuses on the anisotropic aspect of localized failure. It has been discovered that the fabric anisotropy can play an important role in the occurrence of shear banding. An increasing degree of fabric anisotropy tends to delay the initiation of the strain localization and result in higher failure strength. The effects of fabric anisotropy have also been illustrated by comparing the theoretical predictions and measured results on the shear band inclination angle, shear strain level and dilatancy at bifurcation.     

Strength anisotropy of granular material consisting of perfectly round particles

The strength anisotropy of granular materials deposited under gravity has mostly been attributed to elongated particles' tendency to align long axes along the bedding plane direction. However, recent experiments on near‐spherical glass beads, for which preferred particle alignment is inapplicable, have exhibited surprisingly strong strength anisotropy. This study tests the hypothesis that certain amount of fabric anisotropy caused by the anisotropic stress during deposition under gravity can be locked in a circular‐particle deposit. Such locked‐in fabric anisotropy can withstand isotropic consolidation and leads to significant strength anisotropy. 2D discrete element method simulations of direct shear tests on circular‐particle deposits are conducted in this study, allowing for the monitoring of both stress and fabric. Simulations on both monodispersed and polydispersed circular‐particle samples generated under downward gravitational acceleration exhibit clear anisotropy in shear strength, thereby proving the hypothesis. When using contact normal‐based and void‐based fabric tensors to quantify fabric anisotropy in the material, we find that the intensity of anisotropy is discernible but low prior to shearing and is dependent on the consolidation process and the dispersity of the sample. The fact that samples with very low anisotropy intensity measurements still exhibit fairly strong strength anisotropy suggests that current typical contact normal‐based and void‐based second‐order fabric tensor formulations may not be very effective in reflecting the anisotropic peak shear strength of granular materials. Copyright © 2017 John Wiley & Sons, Ltd.     

含横向节理的柱状节理岩体力学特性试验研究

针对柱状节理岩体的构造特性,采用水泥砂浆等材料制作出具有不同柱体倾角且含有横向节理分布的柱状节理岩体模型,并通过单轴压缩试验,研究柱体倾角和横向节理对岩体各向异性力学特性及破坏机制的影响。结果表明：试样峰值强度和变形模量随柱体倾角变化曲线均近似为U形,并体现出典型的各向异性特征;随着柱体倾角的变化,试样单轴压缩破坏模式可分为垂直于柱体轴向的劈裂破坏、沿纵向节理面的剪切滑移破坏、滑移破坏与劈裂破坏同时发生的复合破坏和平行于柱体轴向的劈裂破坏4类;横向节理的存在降低了柱状节理岩体沿柱轴方向的完整性,并影响了劈裂破坏裂纹的分布与扩展,因此,横向节理切割柱体是影响柱状节理岩体承载能力的重要因素。     

An Anisotropic Geomechanical Model for Jointed Rock Masses at Taishan Nuclear Power Station, Southern China

An anisotropic geomechanical model for jointed rock mass is presented. Simultaneously with deriving the orthotropic anisotropy elastic parameters along the positive axis, the equivalent compliance matrix for the deflection axis orthotropic anisotropy was derived through a three-dimensional coordinate transformation. In addition, Singh’s analysis of the stress concentration effects of intermittent joints was adopted, based on two groups of intermittent joints and a set of cross-cutting joints in the jointed rock mass. The stress concentration effects caused by intermittent joints and the coupling effect of cross-cutting joints along the deflection-axis are also considered. The proposed anisotropic mechanics parameters method is applied to determine the deformation parameters of jointed granite at the Taishan Nuclear Power Station. Combined with the deterministic mechanical parameters of rock blocks and joints, the deformation parameters and their variability in jointed rock masses are estimated quantitatively. The computed results show that jointed granite at the Taishan Nuclear Power Station exhibits typical anisotropic mechanical characteristics; the elastic moduli in the two horizontal directions were similar, but the elastic modulus in the vertical direction was much greater. Jointed rock elastic moduli in the two horizontal and vertical directions were respectively about 24% and 37% of the core of rock, showing weakly orthotropic anisotropy; the ratio of elastic moduli in the vertical and horizontal directions was 1.53, clearly indicating the transversely isotropic rock mass mechanical characteristics. The method can be popularized to solve other rock mechanics problems in nuclear power engineering.     

An anisotropically elastoplastic solution to excavation responses of a circular opening considering three-dimensional strength

Circular opening is commonly encountered in wellbore drilling of petroleum engineering, boring for cast-in situ pile installation, and tunneling excavation. This paper presents a rigorous solution for the elastoplastic responses of the anisotropic soft soil mass around a circular opening excavated under undrained and drained conditions. Both the anisotropic elastoplastic behavior and the 3D strength of the soft clay are incorporated in the present solutions. The well-established anisotropic critical state elastoplastic model S-CLAY1, which can represent the initial fabric anisotropy and stress-induced anisotropy of soft soil, is further modified by the Spatially Mobilized Plane criterion to consider the 3D strength of geomaterials. Then, the investigated problems, excavation of a circular opening under both short-term (undrained) and long-term (drained) conditions, are formulated as a system of first-order differential equations and are solved as initial value problems. The distributions of stress components and anisotropy parameters around the opening, the stress trajectory of a soil particle at the opening wall, as well as the stress–displacement curve at the opening wall are presented to investigate the elastoplastic responses of the opening. Extensive parameters show that the overconsolidated ratio and coefficient of earth pressure at rest (K₀) have remarkable effects on the elastoplastic responses around a circular opening.

     

考虑统计损伤模型的非均质节理岩体拉剪破坏机制研究

拉剪应力状态极易导致岩体破坏乃至失稳,为研究节理岩体拉剪破坏规律,开展了拉剪荷载下共面非贯通节理岩体变形破坏的理论与数值计算研究。通过自定义考虑岩石统计损伤演化的Mohr-Coulomb和最大拉应力准则模型,编写力学参数服从Weibull分布的fish函数,研究了拉剪条件下非均质节理岩体的破坏模式及破坏规律,讨论了岩石均质度、法向拉应力及剪切速率对岩体破坏模式及其力学性质的影响。结果表明,(1) 拉剪应力状态下节理岩体的破坏模式以张拉破坏为主,加载初期破坏位置分布散乱,随着加载和损伤演化逐渐形成带状破裂面,岩体宏观力学性质明显降低;(2) 非均质性对岩体破坏影响显著,主要表现为均质度的增加,岩体由弥散型破坏向集中型破坏转变,破裂面起伏度增大,同时岩体的宏观力学性质增强并最终趋向于均质岩体;(3) 低应力水平下拉应力增大不改变节理岩体以拉张破坏为主的破裂模式,但剪切破坏比例明显减少,同时岩体抗剪强度降低,破裂面的粗糙度增大;(4) 剪切速率对岩体力学性质的影响显著,静态加载范围内岩体抗剪强度随剪切速率的增大而增大,且增幅越来越小。     

Anisotropic mechanical properties and brittleness evaluation of layered phylliteEI北大核心

为了研究层状千枚岩的力学特性与各向异性特征,开展了不同层理倾角与不同围压下的千枚岩力学试验。对比分析了千枚岩试样强度、变形、脆性与破坏模式等各向异性特征。结果表明:(1)随着层理倾角β增加,岩样的强度、变形特征曲线形状呈U型;随围压增加,岩样强度及塑性增强,各向异性度逐渐减弱稳定。(2)采用多种强度准则描述岩样强度各向异性,其中Saeidi准则和改进Ramamurthy准则能很好地预测岩样在不同层理倾角下的抗压强度。(3)基于岩样峰前应力-应变曲线与能量特征提出了综合脆性评价指标,在层理倾角β=45°左右时,岩样脆性指标较低,更易发生剪切滑移破坏,得出脆性下降顺序为:沿层理面拉伸劈裂>穿层理面拉伸劈裂>沿层理面剪切>穿层理面剪切。(4)千枚岩的破坏模式与层理倾角和围压有关,单轴条件下,岩样劈裂破坏后易形成复杂裂纹网络;高围压下,岩样破裂后多形成单一的沿层理面或贯穿多层理面的剪切破坏。     

Investigation on strength and stability of jointed rock mass using three‐dimensional numerical manifold method

This paper proposes a numerical model for jointed rock masses within the 3‐D numerical manifold method (NMM) framework equipped with a customized contact algorithm. The strength of rock sample containing a few sets of discontinuities is first investigated. The results of models with simple geometries are compared with the available analytical solutions to verify the developed computer code, whereas models with complex geometries are simulated to better understand the fundamental behavior and failure mechanism of jointed rock mass. Furthermore, the stability of jointed rock mass in an underground excavation is studied, where rock failure process is determined by the 3‐D NMM simulation. The simulation results provide valuable guidance on excavation process design and stabilization design in rock engineering practice. Copyright © 2012 John Wiley & Sons, Ltd.     

峰前卸荷对节理岩体力学特性影响试验研究

实际工程中,岩体的应力变化状态非常复杂,并不只是简单的加载或者卸荷,如水利工程中的大坝,当坝基开挖时其力学条件主要表现为开挖方向的卸荷(其他方向有可能卸荷,也有可能加载),而坝体修筑以后,力学条件又主要表现为加载。因此,对卸荷损伤后的岩体再次加载时的力学性质进行研究具有十分重要的工程意义。通过三轴卸荷及卸荷损伤后加载破坏试验,研究了不同围压下,含断续预制节理岩体经过不同程度卸荷损伤后再次加载时的力学特性。研究结果表明:围压对于含断续节理岩体变形模量影响不明显,岩体弹性模量与围压没有明显关系,相同围压下节理岩体弹性模量和变形模量均显著低于完整岩体:节理岩样峰值应变随着围压的增大线性增大,二者之间有较好的相关性,相比于完整岩样,节理岩体峰值应变随围压增大增长更快:围压越低,卸荷对岩体强度的影响越大,随着围压增加,卸荷量不同造成岩样的损伤差异逐渐减小。     

川西北茂县群千枚岩各向异性力学特性

利用MTS815Teststar岩石力学试验系统,对川西北茂县群千枚岩开展单轴和三轴试验,研究千枚岩各向异性力学特性。研究结果表明：千枚岩在达到峰值强度50%以前就普遍出现很大的变形和局部破坏,表明围岩产生大变形或出现早期扩容并不意味着其丧失承载能力,如及时进行针对性支护或加强围岩保护仍具有较强的强度。损伤扩容是千枚岩变形破坏的典型特征,并对变形产生不同的影响,横向应变对损伤扩容比轴向应变更敏感,能更灵敏地反映岩石材料内部的屈服、弱化。千枚岩各向异性随含水率与围压水平表现出不同的变化规律,强度各向异性随含水率和围压水平提高而降低,弹性模量和泊松比则与加载方向和转化围压有关。不同含水状态下干燥千枚岩的单轴抗压强度各向异性最为显著,大部分千枚岩的强度各向异性率在0.80以下,其中绢云母千枚岩最强,为0.48,其强弱对比依次为绢云母千枚岩、绿泥石千枚岩、石英千枚岩、炭质千枚岩。千枚岩表现出不同的破坏模式,并随围压或含水率的增加而变化;剪切机制单独控制着各向异性力学特性,千枚岩统一沿软弱面剪切破坏。     

A particle mechanics approach for the dynamic strength model of the jointed rock mass considering the joint orientation

In nature, there exist several forms of anisotropy in rock masses due to the presence of bedding planes, joints, and weak layers. It is well understood that the anisotropic properties of jointed rock masses significantly affect the stability of surface and underground excavations. However, these critical anisotropic characteristics are often ignored in existing uniaxial dynamic failure criteria. This study investigates the effect of a pre-existing persistent joint on the rate-dependent mechanical behaviours of a rock mass using a particle mechanics approach, namely, bonded particle model (BPM), to realistically replicate the mechanical response of the rock mass. Firstly, in order to capture the rate-dependent response of the jointed rock mass, the BPM model is validated using published experimental data. Then, a dynamic strength model is proposed based on the Jaeger criterion and simulation results. To further investigate the dynamic behaviours, the dynamic uniaxial compressive strength (UCS) for anisotropic rock masses with various joint orientations is investigated by subjecting the BPM models to uniaxial compression numerical tests with various strain rate. The proposed dynamic strength model is validated based on numerical simulation results. Finally, the fragmentation characteristics of the jointed rock masses are analysed, which demonstrate that the failure mode affects the dynamic UCS. This is further confirmed by the analysis of the orientations of microscopic cracks generated by the compression loading.