综合考虑宏微观复合损伤的节理岩体本构模型

针对目前节理岩体损伤本构模型中仅考虑节理等宏观缺陷造成的损伤,而没有考虑岩块中微裂纹等微观损伤的不足,提出了综合考虑宏微观损伤的节理岩体本构模型。其中微观损伤模型采用基于应变强度理论和岩石微元强度服从Weibull分布假定的统计损伤模型,把其应用于被节理切割而成的岩块。而宏观损伤模型主要基于连续介质力学原理,把节理对岩体性质的劣化作为损伤来考虑。在考虑节理传压及传剪系数的基础上,提出了综合考虑宏微观复合损伤的节理岩体本构模型。算例表明宏观节理的存在大大削弱了岩体的强度,降低了其刚度,增大了其柔性。所提出的模型能够较好地反映宏微观两类损伤对岩体应力-应变关系及强度的影响,较为合理。      

Fully coupled hydro-mechanical numerical manifold modeling of porous rock with dominant fractures

Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. In this study, we developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. Finally, we tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses.     

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.     

Strength and dilatancy of jointed rocks with granular fill

It is well recognised that the strength of rock masses depends upon the strain history, extent of discontinuities, orientation of plane of weakness, condition of joints, fill material in closely packed joints and extent of confinement. Several solutions are available for strength of jointed rock mass with a set of discontinuities. There is a great multiplicity in the proposed relationships for the strength of jointed rocks. In the present study, the author conceives the effect of increasing stresses to induce permanent strains. This permanent strain appears as micro crack, macro crack and fracture. A fully developed network of permanent deformations forms joint. The joint may contain deposits of hydraulic and hydrothermal origin commonly known as gouge. The joint factor numerically captures varied engineering possibilities of joints in a rock mass. The joints grow as an effect of loading. The growth of the joints is progressive in nature. It increases the joint factor, which modifies the failure stresses. The dilatancy explains the progressive failure of granular media. Hence, a mutual relationship conjoins effectively the strength of jointed rock and a dilatancy-dependent parameter known as relative dilatancy. This study provides a simple and integral solution for strength of jointed rocks, interpreted in relation to the commonly used soil, and rock parameters, used for a realistic design of structure on rock masses. It has scope for prediction of an equivalent strength for tri-axial and plane strain conditions for unconfined and confined rock masses using a simple technique.     

节理岩体蠕变特性研究

在已知岩石和节理蠕变规律的前提下,推导了节理岩体蠕变模型的一般表达式。假定岩石体积变形和节理法向压缩变形为弹性变形,忽略节理的剪胀现象,认为只有岩石畸变和节理剪切滑移与时间有关,从而推导了含三组相交节理的岩体蠕变模型及其参数。根据反演出的岩石和节理蠕变模型,计算了含三组相交节理的岩体在单轴应力作用下一些节理参数对岩体单轴蠕变的影响。分析表明,节理间距、剪胀系数以及节理夹角都对岩体的单轴蠕变变形有明显影响。节理间距越大,剪胀系数越大,节理夹角越小,节理岩体的单轴蠕变柔量也就越小,岩体的蠕变变形也越小。     

Modeling strategy for jointed rock masses reinforced by rock bolts in tunneling practice

This paper deals with the modeling of jointed rock masses reinforced by rock bolts. It is well known that rock bolts are extremely effective in reinforcing jointed rocks. However, if a continuum approach is adopted for modeling jointed rock masses, it is often misleading to evaluate the effectiveness of the rock bolts by numerical analyses such as the finite element method. This may be due to the fact that since no more joints exist in the continuum, the effectiveness of the rock bolts in constraining the relative displacements along the joints cannot be evaluated properly. In order to investigate the reinforcement effect of rock bolts, physical model tests were performed in the laboratory. The test results revealed that jointed rock masses should be modeled as an equivalent continuum after the installation of rock bolts and that the mechanical parameters of the equivalent continuum should be evaluated by considering the reinforcement effect of the rock bolts. Therefore, the values of the mechanical parameters differ from place to place in accordance with the relation between joint orientation and rock bolt direction, even though joint systems are homogeneous. In conclusion, in the continuum approach for modeling jointed rock masses, it should be emphasized that rock bolts and jointed rock masses should not be modeled separately but should be modeled simultaneously by considering the reinforcement effect of the rock bolts in constraining joint movement. The modeling of shotcrete reinforced by steel ribs is also discussed in comparison to the modeling by rock bolts.     

基于摩尔-库仑准则的断续节理岩体复合损伤本构模型

针对工程实际中断续节理裂隙岩体的损伤本构模型,假设岩石微元强度服从Weibull随机分布,以摩尔-库仑破坏准则作为描述微元强度的表示方法,推导出细观损伤变量。利用能量和断裂力学理论,综合考虑节理几何特征及力学特性,推导宏观损伤变量计算公式。基于Lemaitre应变等效假设,考虑宏细观缺陷耦合作用,推导出复合损伤变量,建立了基于摩尔-库仑准则的宏细观缺陷耦合作用的断续裂隙岩体损伤本构模型。研究结果表明：（1）采用摩尔-库仑准则作为描述微元强度的统计分布变量建立的损伤模型能够较好地反映岩石内部缺陷分布和变形特征,该模型真实地反映岩石微元强度受应力状态的影响。（2）该模型建立的理论曲线与断续节理岩体的试验曲线吻合较好。（3）节理裂隙岩体宏观损伤变量及峰值强度随节理倾角的变化规律与综合考虑宏细观耦合作用下的损伤变量及裂隙岩体峰值强度随节理倾角的变化规律基本一致。（4）宏细观耦合作用下的等效弹性模量与节理贯通率呈非线性负相关;在节理倾角一定的情况下,损伤变量与节理长度呈非线性正相关;在贯通率较小时,岩体的宏观损伤变量与内摩擦角的关系呈线性负相关变化,贯通率达到一定程度,线性关系变成非线性关系。     

柱状节理岩体原位变形试验力学浅析与模拟

在白鹤滩坝址区柱状节理岩体原位变形试验成果的基础上,针对节理刚度的取值进行讨论,利用柔性中心孔试验值反算得到柱状节理岩体节理面刚度值。岩体变形是由岩块和节理两部分组成,分析含单节理岩体概念模型的单轴加、卸载规律性和其变形特性,发现其滞回特性是节理力学行为造成的;推导得到了岩体变形计算公式,可用于多组相交节理岩体等效弹性模量的计算,能够反映出节理岩体的各向异性特征;建立了柱状节理随机模型,采用离散元方法和反算的刚度值进行了仿真分析,得到等效变形参数在试验结果范围内。     

基于宏细观损伤耦合的非贯通裂隙岩体本构模型

针对非贯通裂隙岩体工程结构中的受荷岩体,提出受荷细观损伤与裂隙宏观损伤的概念。以完整岩石的初始损伤状态作为基准损伤状态,综合考虑裂隙宏观缺陷的存在、微裂纹细观缺陷在受荷下的损伤扩展以及宏细观缺陷在受荷过程中的耦合,基于Lemaitre应变等效假设,推导了考虑宏细观缺陷耦合的复合损伤变量,并给出同时考虑试件尺寸、裂隙几何与力学特性的宏观损伤变量的计算公式,从而建立了基于宏细观缺陷耦合的非贯通裂隙岩体在荷载作用下的损伤本构模型。用宏细观损伤耦合的本构模型来描述非贯通裂隙岩体在受荷过程中的细观损伤演化与宏观损伤行为,与非贯通裂隙岩体实际受荷情况符合较好。研究结果表明：（1）完整岩样和裂隙岩样的应力-应变行为在峰值强度之前差异较大,峰值强度以后差异逐渐减小,最后趋于一致,二者具有相近的残余强度;（2）裂隙岩体强度随裂隙贯通率的增加而增大,随裂隙倾角的变化具有明显的各向异性,同时还与裂隙面的内摩擦角有关;（3）裂隙倾角为90°时,裂隙岩样的峰值强度最高;张开型裂隙岩样的裂隙倾角为45°时,峰值强度最低;（4）非贯通裂隙岩体工程结构中的受荷岩体,其力学性能由受荷细观损伤与裂隙宏观损伤及其耦合效应所决定,基于宏细观损伤耦合的复合损伤变量可以较好地反映非贯通裂隙岩样的力学特性。     

Stability analysis of tunnels driven in jointed rocks by means of a homogenized limit analysis approach

Conceived as a potential alternative to the classical design methods employed for analyzing the stability of underground structures driven in jointed rocks, the homogenization approach stems from the heuristic idea that, from a macroscopic point of view, a rock mass cut by a network of joints may be perceived as a homogenized continuum. The strength properties of the latter can be theoretically obtained from the failure conditions of its individual constituents: rock matrix and joint interfaces. At the material level, the limit analysis reasoning is used in the context of homogenization to formulate the homogenized strength criterion of a jointed rock mass in the particular situation of a single set of parallel joints. As it could be expected, the obtained closed‐form expressions show the strength anisotropy induced by joint preferential orientation. The support functions (π functions) associated with the homogenized strength criterion are also determined in both plane strain and three‐dimensional cases. This criterion is then applied to the investigation of stability analysis of a tunnel excavated in a jointed rock mass. Upper bounds estimated of the stability factor are derived from the implementation of the kinematic approach directly on the homogenized underground structure. Finally, the approach is applied to analyze and discuss the collapse of the Pinheiros subway station (São Paulo, Brazil). Copyright © 2014 John Wiley & Sons, Ltd.     

柱状节理玄武岩隧洞破坏模式及其力学机制模拟

柱状节理岩体由于其内部赋存大量的隐节理面,开挖卸荷后极易出现隐节理面开裂松弛等特征,导致其破坏模式异于一般岩体。其破坏模式主要受到异常发育的节理面和较高地应力的共同作用影响。由于柱状节理岩体节理面发育,岩体结构控制型破坏占主要部分,包括单临空面节理面滑移（塌方）、多临空面节理面滑移（塌方）、与错动带、断层等弱面相组合的坍塌等破坏模式;应力控制型破坏主要为河谷侧顶拱喷层开裂;应力-结构面型破坏主要为岩性交界处的节理岩体塌落等。柱状节理岩体表层主要发生柱内竖直隐节理面和柱间节理面的拉破坏,而围岩内部的柱状节理主要发生柱间节理面的剪切破坏。因此,现场柱状节理的支护也应主要包括两个方面：以喷射混凝土钢纤维来阻止柱状节理岩体表层的张性开裂,以系统锚杆来控制柱状节理岩体内部的剪切破坏。     

基于离散元法的节理岩体边坡稳定性分析

节理岩体边坡的稳定性在很大程度上取决于节理的强度及其分布形式。由于节理岩体边坡的失稳破坏具有大变形和非连续的特点,因此,离散单元法成为研究节理岩体边坡破坏机制的最有效方法之一。通过采用离散元软件PFC2D进行数值模拟,对完整岩石及节理的力学性能进行研究,并建立含密集节理的岩体边坡模型,讨论了节理连通率对边坡破坏形式的影响。结果表明,节理岩体边坡的失稳破坏是一个渐进的过程;在多组节理密集分布的岩体边坡中,连通率越大,其稳定性越差;随着连通率的减小,边坡的破坏形式由大范围的滑坡转变为局部崩塌的形式     

Composite damage constitutive model of jointed rock mass considering crack propagation length and joint friction effect

In order to study the damage constitutive model of rock mass with non-persistent joints in engineering practice, the assumption that the rock mesoscopic elements strength obeys the Weibull distribution function of random, the Drucker–Prager criterion is used as a representation method to describe the strength of the mesoscopic elements, and mesoscopic damage variable is deduced. Combined with the energy principle and the fracture damage theory, the macroscopic damage variable formula is deduced considering the crack propagation length and the friction effect of joint closure in the rock mass. Finally, based on the strain equivalent hypothesis of Lemaitre and considering the coupling of macro and micro defects, the composite damage variables are derived. A macro–meso composite damage constitutive model of rock mass with non-persistent joints is established based on the Drucker–Prager criterion. The theoretical constitutive curves of the model are in good agreement with the experimental constitutive curves of the non-persistent jointed rock masses. The considering effects of crack propagation length and joint closure friction effect are compared with those without considering the crack propagation length and joint closure friction effect, which finds that the former is superior to the latter, and the rationality and validity of the model is verified.     

A finite element slideline approach for calculating tunnel response in jointed rock

The general purpose finite element code DYNA2D is used to calculate the response of an unlined tunnel in jointed rock, with sliding interfaces used to represent the joints. The model problem selected for analysis closely follows a calculation using the discrete element code UDEC. Results for two cases [(1) jointed rock with lubricated joints (no joint friction) and (2) jointed rock with frictional resistance to joint motion] are presented and compared with each other and the UDEC results.     

基于离散裂隙网络模型的节理岩体渗透张量及特性分析

节理岩体几何结构非常复杂,研究其渗流特性对于指导含水岩层稳定性分析具有重要价值。应用离散裂隙网络模型DFN方法,基于VC++6.0软件平台,建立了平面渗流分析方法,分析了节理岩体不同几何分布情况下的渗透率张量特征,通过定义渗流定向性系数对岩体渗流的定向性特征进行了定量分析。结果表明：单组节理岩体渗流具有明显的各向异性特征,渗流定向性随着节理角度变化显著;节理随着节理贯通性增加,节理渗透率呈现对数增加趋势;两组节理情况下,各向异性特征随着节理组间夹角变化;两组节理岩体渗流特征研究中,正交分布下,岩体仍存在各向异性,但渗流定向性系数较低;当节理倾角服从正态分布时,随着节理倾角标准差增大,渗透率增加;两组节理夹角不同时,节理渗透主方向倾角随着夹角增大而相应增大,基本沿两组节理夹角方向的角平分线方向。     

A Numerical Study on Wave Transmission Across Multiple Intersecting Joint Sets in Rock Masses with UDEC

This paper presents a numerical study on wave transmission across jointed rock masses with UDEC, where multiple intersecting joint sets exist. The capability of UDEC of studying wave transmission across rock joints is validated through comparison with analytical solutions and experimental data. Through parametric studies on wave transmission across jointed rock masses, it is found that joint mechanical and spatial parameters including joint normal and shear stiffnesses, nondimensional joint spacing, joint spacing ratio, joint intersecting angle, incident angle, and number of joint sets together determine the wave transmission. And for P wave incidence, compared with other parameters, joint normal stiffness, nondimensional joint spacing, and joint intersecting angle have more significant effects on wave transmission. The physical reasons lying behind those phenomena are explained in detail. Engineering applications and indications of the modeling results are also mentioned.     

单轴压缩荷载下非贯通闭合节理岩体损伤本构模型

目前损伤力学已被认为是研究节理岩体力学行为的有效工具,但是在目前的节理岩体损伤变量定义中大多仅考虑节理几何特征而未考虑节理内摩擦角等力学参数,这显然不能很好地反映节理岩体的力学特征。为此,拟推导出一个能够综合考虑节理几何及力学参数的损伤变量（张量）,并由此建立单轴压缩荷载下岩体损伤本构模型。首先,基于断裂力学的由于单个节理存在引起的附加应变能增量与损伤力学的损伤应变能释放量相关联的观点,推导出了含非贯通节理岩体的损伤变量计算公式;其次,根据断裂力学理论对单轴压缩荷载下的单个节理尖端应力强度因子计算方法进行了研究,得出了应力强度因子K_Ⅰ、K_Ⅱ的计算公式;同时考虑多节理间的相互作用给出了单组单排及多排非贯通节理应力强度因子计算公式。最后,利用该模型对含单条非贯通节理的岩体在单轴压缩荷载作用下的峰值强度及损伤变量进行了分析计算。结果表明,当节理倾角小于其内摩擦角时,岩体强度与完整岩石相同,岩体损伤为零,而后随着节理倾角增加,岩体强度、损伤随节理倾角的变化分别呈开口向上及向下的抛物线,当节理倾角约为60°时,岩体损伤最大,强度最低。随着节理长度增加,岩体损伤增加,而随着节理内摩擦角的增加,岩体损伤则减小。      

Rock joint modeling using a visco-plastic multilaminate model at constant normal load condition

Rock joints play an important role in the behavior of rock masses under normal and shear loading conditions. Numerical simulation of the behavior of jointed rock masses is not an easy task due to complexities involved in the problem such as joint roughness, joint shear strength, hardening and softening phenomenon and mesh dependency. In this study for modeling purposes, a visco-plastic multilaminate model considering hardening and softening effects has been employed. For providing the necessary data for numerical simulation, a series of laboratory experiments have been carried out on regular tooth-shape asperities made by gypsum, under constant normal load conditions. Shear stress–shear displacement and normal displacement–shear displacement of artificial joint specimens are simulated using the proposed numerical model at constant normal load condition (CNL). The results indicate the capability of the model for simulating rock joints behavior in both strength and deformation field. Although the numerical model has been developed for simulating the behavior of artificial joints, the concept of the method can also be used for natural rock joints.     

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

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

High lateral strain ratio in jointed rock masses

Rock masses encountered in the field are jointed and highly anisotropic. An important manifestation of this characteristic of rock masses is that the blocks, when slipping, separate and create voids in the mass resulting in high lateral strain even under low axial stress. It has been shown in this study that the ratio of lateral to axial strain may be very high, especially, if the joints are critically oriented. The assumption of isotropic linearly elastic material is not applicable in such situations. This observation is based on the outcome of an extensive laboratory testing programme, in which a large number of specimens of a jointed rock mass with various joint configurations were tested under uniaxial loading conditions. The trends of experimental results for both lateral strain ratio and rock mass strength have also been verified through distinct element modelling. The reason for high lateral strains has been attributed to the creation of voids and also to the fact that permanent deformations due to slip commence along rock joints right from the start of loading process. A simple mechanistic model has also been suggested to explain the high values of lateral strain for rough and dilatant rock joints. The lateral strains in such situations are important in the design of rock bolts. An example from a recently completed hydroelectric project has also been discussed where high lateral strain was found to be responsible for excessive wall closure and failure of rock bolts.