Microplane damage model for jointed rock masses

The paper presents a new microplane constitutive model for the inelastic behavior of jointed rock masses that takes into account the mechanical behavior and geometric characteristics of cracks and joints. The basic idea is that the microplane modeling of rock masses under general triaxial loading, including compression, requires the isotropic rock matrix and the joints to be considered as two distinct phases coupled in parallel. A joint continuity factor is defined as a microplane damage variable to represent the stress‐carrying area fraction of the joint phase. Based on the assumption of parallel coupling between the rock joint and the rock matrix, the overall mechanical behavior of the rock is characterized by microplane constitutive laws for the rock matrix and for the rock joints, along with an evolution law for the microplane joint continuity factor. The inelastic response of the rock matrix and the rock joints is controlled on the microplane level by the stress–strain boundaries. Based on the arguments enunciated in developing the new microplane model M7 for concrete, the previously used volumetric–deviatoric splits of the elastic strains and of the tensile boundary are avoided. The boundaries are tensile normal, compressive normal, and shear. The numerical simulations demonstrate satisfactory fits of published triaxial test data on sandstone and on jointed plaster mortar, including quintessential features such as the strain softening and dilatancy under low confining pressure, as well as the brittle–ductile transition under higher confining pressure, and the decrease of jointed rock strength and Young's modulus with an increasing dip angle of the joint. Copyright © 2014 John Wiley & Sons, Ltd.     

The Influence of Shearing Velocity on Shear Behavior of Artificial Joints

In this paper, the effects of shear velocity on the shearing behavior of artificial joints have been studied at different normal stress levels. Here, artificial joints with planar and rough surfaces were prepared with the plaster (simulating soft rock joints) and concrete (medium-hard rock joints) materials. The rough joints had triangular shaped asperities with 10° and 20° inclination angles. Direct shear tests were performed on these joints under various shear velocities in the range of 0.3–30 mm/min. The planar plaster–plaster and planer concrete–concrete joints were sheared at three levels of normal stress under constant normal load boundary condition. Also, the rough plaster–plaster and concrete–concrete joints were sheared at one level of normal stress under constant normal stiffness boundary condition. The results of the shear tests show that the shearing parameters of joints, such as shear strength, shear stiffness and friction angle, are related to the shear velocity. Shear strength of planar and rough plaster–plaster joints were decreased when the shear velocity was increased. Shear strength of concrete joints, except for rough joints with 10° inclination, increased with increasing shear velocity. Regardless of the normal stress level, shear stiffness of both planar plaster–plaster and concrete–concrete joints were decreased when the shear velocity was increased.     

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

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

NON-LINEAR FAULT BEHAVIOUR NEAR UNDERGROUND EXCAVATIONS — A BOUNDARY ELEMENT APPROACH

A boundary element model for stress/stability analysis of underground excavations in the vicinity of faults is presented. The boundary element formulation adopts the fictitious stress method for the simulation of excavation boundaries and the displacement discontinuity method for the representation of faults. The numerical model employs the Barton–Bandis non-linear joint model for the modelling of the fault behaviour and linear elastic behaviour for the rock. An incremental-iterative in situ stress relaxation algorithm is implemented for the non-linear analysis of the faults. Both deformation and peak strength models of Barton–Bandis are incorporated for modelling the mechanical behaviour of the fault. The non-linear deformation of fault considers the effects of coupling between shear and normal stresses and displacement, joint closure, joint separation, hardening followed by post-peak or residual behaviour. The peak strength model employs a mobilized non-linear shear strength envelope. The differences between linear and non-linear simulation of the fault models are discussed. A comparison of model predictions with the classical Mohr–Coulomb peak strength model with constant joint stiffness is presented. The numerical model is used for a case study of Canadian hard rock underground mine. The shear and normal displacements along the fault during four mining sequences with backfill simulation are presented and discussed.     

A Constitutive Model for Shear Behavior of Rock Joints Based on Three-Dimensional Quantification of Joint Roughness

A new constitutive model to describe the shear behavior of rock joints under constant normal stiffness (CNS) and constant normal load (CNL) conditions is proposed. The model was developed using an empirical approach based on the results of a total of 362 direct shear tests on tensile fractured rock joints and replicas of tensile joints and on a new quantitative roughness parameter. This parameter, the active roughness coefficient C _r, is derived from the features of the effective roughness mobilized at the contact areas during shearing. The model involves a shear strength criterion and the relations between stresses and displacements in the normal and shear directions, where the effects of the boundary conditions and joint properties are considered by the shape indices C _d and C _f. The model can be used to predict the shear behavior under CNS as well as CNL conditions. The shear behavior obtained from the experimental results is generally in good agreement with that estimated by the proposed model, and the effects of joint roughness, initial normal stress, and normal stiffness are reasonably reflected in the model.     

Critical review on shear strength models for soil-infilled joints

An infilled rock joint is likely to be the weakest plane in a rock mass. The presence of infill material within the joint significantly reduces the friction of the discontinuity boundaries (i.e. rock to rock contact of the joint walls). The thicker the infill, the smaller the shear strength of the rock joint. Once the infill reaches a critical thickness, the infill material governs the overall shear strength, and the joint walls (rock) play no significant role. Several models have been proposed to predict the peak shear strength of soil-infilled joints under both constant normal load (CNL) and constant normal stiffness (CNS) boundary conditions, taking into account the ratio of infill thickness (t) to the height of the joint wall asperity (a). CNS models provide a more realistic picture of the soil-infilled joint behaviour in the field. This paper presents a critical review on the existing mathematical models for predicting the shear strength of soil-infilled rock joint and verifies the normalised peak shear stress model with further laboratory investigations carried out on idealised saw-tooth rock joints at the University of Wollongong. Based on the prediction of the experimental data, the normalised peak shear stress model is slightly modified by the authors. A simplified approach for using this model in practice is presented and a new expression for prediction of dilatation at peak shear stress is suggested.     

不同边界条件下剪切速率对类岩石节理剪切力学特性的影响

为研究不同边界条件下剪切速率对岩石节理剪切力学特性的影响,采用RDS-200型岩石节理剪切试验系统对人工浇筑的具有相同节理形貌的不规则水泥节理试样进行了常法向应力和常法向刚度2种边界条件下5种剪切速率的直剪试验。结果表明：（1）常法向应力边界条件下,随剪切速率增大,相同法向应力下的类岩石节理峰前剪切刚度减速增大,峰值剪切强度和残余剪切强度呈对数降低;随剪切速率增大,类岩石节理黏聚力减速增大,内摩擦角呈对数降低。（2）常法向刚度边界条件下,随剪切速率增大,相同法向应力的类岩石节理峰前剪切刚度减速增大,峰值剪切强度呈对数降低,较高法向应力下的残余剪切强度先增大后减小;随剪切速率增大,类岩石节理黏聚力呈对数降低,内摩擦角减速增大。（3）与常法向应力边界条件相比,常法向刚度条件下,节理黏聚力平均增加了115.85%,内摩擦角平均降低了8.44%;相同初始法向应力和剪切速率下,峰前剪切刚度、峰值剪切强度和残余剪切强度分别平均增加了11.96%、19.47%和32.32%,峰值法向位移平均降低了40.12%。该研究结论可为不同剪切速率下地表和地下工程岩体节理的剪切失稳评价提供一定参考。     

不同粗糙度的节理直剪颗粒流分析

利用二维颗粒流程序生成5种不同粗糙程度的节理模型,并对5种节理模型进行了5种不同法向恒定荷载作用下的直剪试验,从细观角度分析了不同粗糙程度的节理模型在法向荷载下的形貌损伤情况和裂纹演化机制。与此同时,分析了节理JRC值和节理面颗粒摩擦系数对节理抗剪强度影响,并反推出了节理面抗剪强度参数Cj与?j与JRC值的关系。结果为：法向恒定荷载越大时,节理峰值抗剪应力越大,剪胀现象越小,节理形貌损伤范围越大。随着剪切的进行,上下节理面接触范围减小,微裂纹开始主要沿节理面产生,随着剪切位移的继续增加微裂纹数量显著增加,并且不局限于节理面附近而深入到模型内部。随着节理粗糙程度（JRC值）和节理面颗粒摩擦系数的增加节理峰值抗剪应力也增大。节理抗剪强度参数Cj与?j随着JRC值的增大而增大。所得结果可以为室内试验和工程应用提供参考和依据。     

峰前循环剪切作用下岩石节理剪切力学特性

微震、工程爆破等低应力循环剪切荷载作用对节理岩体工程失稳破坏具有重要影响。为研究峰前循环剪切加卸载作用下岩石节理剪切力学特性,采用RDS-200型岩石节理剪切试验系统对人工劈裂黄砂岩节理进行了峰前循环剪切下的直剪试验。通过与未进行峰前循环剪切加卸载时岩石节理力学参数预测值对比,得到峰前循环剪切加卸载作用对峰前剪切刚度、峰值剪切强度、峰值剪切位移与残余剪切强度的影响。结果表明：（1）峰前循环剪切加卸载后,当法向应力为2 MPa时,岩石节理峰前剪切刚度增大,当法向应力为4～10 MPa时,岩石节理峰前剪切刚度在循环剪切应力幅值范围内增大,在超出循环剪切应力幅值时减小;（2）峰前循环剪切加卸载后,峰值剪切强度降低了10%～20%,降低百分比随法向应力增大整体呈对数函数增大;峰值剪切位移增加了2%～40%,增加百分比随法向应力增大整体呈对数函数减小;（3）峰前循环剪切加卸载后,岩石节理残余剪切强度无明显变化,峰值剪切强度与残余剪切强度差值减小,峰后剪切应力做功损失百分比降低。     

考虑峰后力学特性的岩石节理渗流广义立方定理

通过建立考虑峰后剪胀及剪切软化效应的岩石节理法向变形模型,基于渗流扩散能等效原理,提出了一个岩石节理渗流广义立方定理,较为合理地描述了节理在较大法向闭合变形条件下的渐进渗流特性以及节理出现剪胀及剪切软化条件下的峰后渗流特性。采用已有压剪渗流耦合试验对广义立方定理进行了验证,表明采用常系数修正的立方定理能较好表征节理的峰后渗透特性。     

非贯通类节理岩体剪切变形规律分析

通过直剪条件下的模型试验研究具有相同连通率、含不同起伏角的共面闭合节理在不同法向应力作用下的剪切变形特征。试验结果表明：在达到峰值剪切强度前,剪切应力-位移曲线具有明显的线性段与非线性段;岩桥在法向与剪切荷载的共同作用下会产生一定程度的弱化,积累到一定程度时宏观表现为剪切应力软化。详细阐述了试验过程中裂纹发展的4个阶段并描述了剪切应力-剪切变形曲线变化的特征,提出一个能反映剪切应力-剪切变形曲线全过程的经验本构模型。新模型采用分段函数描述曲线变化特征。直剪试验数据分析表明,该模型能够拟合剪切应力-剪切变形变化发展趋势且具有较高的拟合精度,验证了模型的正确性。最后,对其存在的不足亦进行了简要分析     

岩体结构面剪切流变试验及模型参数反演分析

针对锦屏二级水电站地下洞室群围岩富含节理的实际情况,采用弱面直剪流变仪对组成隧洞围岩的主要岩样（板岩、大理岩）进行了含软弱夹层岩石的剪切流变试验,得出岩石夹层的剪应力-剪切位移时程曲线。通过对岩石夹层试验结果的分析,得出了岩石的长期强度参数,与快剪试验获得抗剪强度参数相比较发现长期抗剪强度参数有所降低,且凝聚力对剪切流变特性的影响略高于内摩擦系数。最后通过对6种模型的拟合分析,得出了适合岩石夹层的标准线性体流变模型,在此基础上通过反演分析得到了岩石夹层标准线性体黏弹-塑性剪切流变模型,并反演了该模型的参数。     

岩桥直剪细观破坏特征与剪胀效应研究

节理岩体的剪切贯通机制影响着边坡的稳定性。为揭示锁固段型非贯通节理岩体在不同连通率和法向应力下的破坏特征,在室内直剪试验中结合高速摄影与AE特征参数分析其剪切全过程及剪胀效应。结果表明:节理岩体直剪试验中,法向应力的增大及节理连通率的下降会致使峰值剪切应力及峰值剪切位移增大;节理连通率与法向应力对其破坏特征具显著影响,表现为节理连通率较高且法向应力较小时呈直接剪断的特性,节理连通率降低后呈拉剪复合破坏,出现剪胀现象,而法向应力的增大使得剪胀效应呈波动现象;AE特征与岩桥贯通过程一致,事件数峰值随节理连通率的降低及法向应力的增大而增大且位于峰后。试验得到的岩桥细观破坏特征与剪胀效应对研究锁固段型岩质边坡的贯通破坏机制具指导意义。     

Importance of Tensile Strength on the Shear Behavior of Discontinuities

In this study, the shear behavior of discontinuities possessing two different rock wall types with distinct separate compressive strengths was investigated. The designed profiles consisted of regular artificial joints molded by five types of plaster mortars, each representing a distinct uniaxial compressive strength. The compressive strengths of plaster specimens ranged from 5.9 to 19.5 MPa. These specimens were molded considering a regular triangular asperity profile and were designed so as to achieve joint walls with different strength material combinations. The results showed that the shear behavior of discontinuities possessing different joint wall compressive strengths (DDJCS) tested under constant normal load (CNL) conditions is the same as those possessing identical joint wall strengths, but the shear strength of DDJCS is governed by minor joint wall compressive strength. In addition, it was measured that the predicted values obtained by Barton’s empirical criterion are greater than the experimental results. The finding indicates that there is a correlation between the joint roughness coefficient (JRC), normal stress, and mechanical strength. It was observed that the mode of failure of asperities is either pure tensile, pure shear, or a combination of both. Therefore, Barton’s strength criterion, which considers the compressive strength of joint walls, was modified by substituting the compressive strength with the tensile strength. The validity of the modified criterion was examined by the comparison of the predicted shear values with the laboratory shear test results reported by Grasselli (Ph.D. thesis n.2404, Civil Engineering Department, EPFL, Lausanne, Switzerland, 2001). These comparisons infer that the modified criterion can predict the shear strength of joints more precisely.     

Applicability of existing models to predict the behavior of replicas of natural fractures of welded tuff under different boundary conditions

Assessing the shear behavior of intact rock and rock fractures is an important issue in the design of a potential nuclear waste repository at Yucca Mountain, Nevada. Cyclic direct shear experiments were conducted on replicas of three natural fractures and a laboratory-developed tensile fracture of welded tuff. The tests were carried out under constant normal loads or constant normal stiffnesses with different initial normal load levels. Each test consisted of five cycles of forward and reverse shear motion. In this paper, the results of the constant normal load shear experiments are analyzed using several constitutive models proposed in the rock mechanics literature for joint shear strength, dilatancy, and joint surface damage. It is shown that some of the existing models have limitations. New constitutive models are proposed and are included in a mathematical analysis tool that can be used to predict joint behavior under various boundary conditions. © Rapid Science Ltd. 1998     

Effects of Cyclic Loading on the Shear Behaviour of Infilled Rock Joints Under Constant Normal Stiffness Conditions

The variation of the shear strength of infilled rock joints under cyclic loading and constant normal stiffness conditions is studied. To simulate the joints, triangular asperities inclined at angles of 9.5° and 18.5° to the shear movement were cast using high-strength gypsum plaster and infilled with clayey sand. These joints were sheared cyclically under constant normal stiffness conditions. It was found that, for a particular normal stiffness, the shear strength is a function of the initial normal stress, initial asperity angle, joint surface friction angle, infill thickness, infill friction angle, loading direction and number of loading cycles. Based on the experimental results, a mathematical model is proposed to evaluate the shear strength of infilled rock joints in cyclic loading conditions. The proposed model takes into consideration different initial asperity angles, initial normal stresses and ratios of infill thickness to asperity height.     

非规则砂岩节理剪切变形本构关系试验研究

岩石节理剪切变形对岩体工程的安全性和稳定性具有重要影响。为研究常法向应力下岩石节理剪切变形本构关系,采用RDS-200型岩石直剪仪对非规则砂岩节理进行了不同法向应力下的直剪试验。根据岩石节理剪切应力?位移全程曲线形状特征,将其依次划分为峰前压密阶段、线性阶段、屈服阶段和峰后软化阶段;根据剪切应力在峰后软化阶段降低幅度和速率大小,将岩石节理剪切应力?位移曲线划分为3种类型：峰后平台型、峰后缓降型和峰后跌落型。基于岩石节理剪切应力?位移曲线各阶段剪切变形特征,采用分段函数建立了岩石节理剪切变形本构模型。与其他模型相比,新提出的岩石节理剪切变形本构模型对试验数据拟合精度更高,更好地描述了岩石节理剪切应力?位移全程曲线。另外,在通过岩石节理直剪试验由经验公式确定模型参数之后,所提出本构模型可在不同法向应力下实现对不同粗糙度岩石节理剪切应力?位移曲线的预测。研究结果对岩石节理剪切变形的数值模拟和工程估计具有一定的实用价值。     

New Methodology to Characterize Shear Behavior of Joints by Combination of Direct Shear Box Testing and Numerical Simulations

A new procedure is presented, which combines big shear box tests on rocks and corresponding numerical simulations with explicit consideration of joint roughness to get deeper insight into the shear behavior of rock joints. The procedure consists of three parts: (1) constant normal load- or CNS-shear box tests with registration of shear- and normal-components of stress and displacements and deduction of basis rock mechanical parameters; (2) high resolution 3D-scanning of joint surface to deduce joint topography; and (3) set-up, run and evaluation of 3-dimensional numerical model with explicit duplication of joint roughness as back-analysis of shear box tests. The numerical back-analysis provides deeper insight into the joint behavior at the micro-scale. Several parameters can be deduced, like micro-slope angle distribution, aperture size distribution, local normal stress distribution and detailed analysis of dilation in relation to shear direction. The potential of the new procedure is illustrated exemplary by shear box tests on slate.     

Failure Characteristics of Complicated Random Jointed Rock Mass Under Compressive-Shear Loading

Most natural rock masses contain a large number of random joints and fissures, and most of the rock masses at the rock engineering are commonly in both compression and shear stress environment. However, the research on the failure characteristics of complex random jointed rock mass under compressive-shear loading is still limited. To address this gap, this paper uses the particle flow code 2D to establish a discrete fractured rock mass model and carry out a series of numerical tests with different compressive-shear angles (α) and different joint geometric parameters. The effects of compressive-shear angle and joint geometric parameters on the strength and failure characteristics of fractured rock masses are studied. The results indicate that with the increase of α, the peak strength of the specimen decreases gradually, and the failure mode changes from the composite shear failure mode (Mode-I) to a plane shear failure mode (Mode-II) and then to intact shear failure mode (Mode-III). Specifically, the three failure modes occur in the specimens with α?=?15°, 30° or 45°, 60°, respectively. The existence of joints affects stress distribution on rock mass during the loading process. Furthermore, the stress at the joint tip is relatively concentrated, while on both sides of the joint is smaller. Three kinds of crack coalescence patterns are observed: tensile, shear, and tensile-shear mixed coalescence. The inclination angle of the rock bridge between adjacent joints affects the specific type of coalescence.