节理岩体蠕变特性研究

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

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

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

含两组交叉贯通节理岩体的强度及破坏特征分析

天然岩体中常含有多组相互交叉的贯通节理,它们的存在极大地削弱了岩体的力学性质。为了研究含两组交叉贯通节理岩体的强度及破坏特征,基于弹塑性数值流形方法,对不同应力状态下,节理倾角和节理间距不同的岩体压缩试验进行数值模拟。结果表明：岩体强度随节理倾角的变化曲线呈现出多波峰、多波谷特点,岩体强度随节理间距变化曲线符合负指数函数形式。根据节理状态,两组节理岩体的破坏模式可以分为3种：岩块破坏、沿一组节理滑移和沿两组节理滑移。两组节理均会影响岩体强度,一组起主要作用,另一组起次要作用,并且节理组之间存在相互影响。通过对数值计算结果进行回归分析,基于一组节理岩体强度预测模型,量化节理组之间的相互影响,提出适用于含两组交叉贯通节理岩体的强度预测模型。模型形式简单,使用方便,可为实际工程中正确评估岩体强度提供指导性意见。     

Modelling of intact and jointed mudstone samples under uniaxial and triaxial compression

Rock mass is a highly complex entity where the strength and deformation behaviour can be significantly affected by its secondary structures such as joints, fissures and bedding planes. Whilst many research works have been conducted to study the behaviour of a specific rock mass, a thorough understanding of its strength and deformation behaviour incorporating different joint sets has not been established. In this study, a comprehensive numerical modelling using a three-dimensional distinct element code, 3DEC, was undertaken to study the strength and deformation behaviour of a mudstone, locally found in Melbourne, in unconfined and confined states. The initial unconfined model established for intact mudstone was calibrated against the well-established laboratory-based empirical strength relationships and subsequently compared with some strength test data available for field samples. The intact unconfined model was then extended to study the strength behaviour in the confined state. The results obtained from this confined intact model were compared with existing strength criteria and were found in good agreement. The confined model was extended further to investigate the effects of joint sets and dip angles on the rock mass strength and deformation behaviour by incorporating two different joint configurations (one-joint and two-joint) with varying dip angles (0°–90°). This study found that the rock mass strength in a confined state varied significantly between the two joint configurations.     

基于DIC的3D打印交叉节理试件破裂机制研究

为研究交叉节理的存在对岩体破裂机制的影响,采用3D打印技术制备了可模拟岩体交叉结构面的节理模型,通过相似材料浇筑形成含预制交叉节理的试件,并基于数字图像相关技术(DIC)分析了单轴压缩下试件裂纹起裂、扩展及破坏模式。研究结果表明:采用3D打印技术制作的闭合节理模型,可以有效替代传统的切割、插缝等方法形成的张开型节理裂隙。基于此展开的试验结果表明,交叉节理存在会显著降低试件强度,且随着交叉角度的增大,试件强度先升高后降低,在45°和60°之间达到最大值,峰值应变与强度呈现相反的变化规律;裂纹扩展过程可分为微裂隙闭合阶段、微破裂发展阶段、主节理起裂扩展阶段和次节理迅速延伸阶段,这与应力-应变曲线各阶段一一对应。研究还发现,在交叉节理岩体中,次节理对破裂的影响主要体现在峰后阶段,结合最大畸变能理论表明,其对主节理尖端应力分布影响较小,主节理对岩体的破坏起绝对控制作用,这对岩体工程有一定的指导意义。     

Applicability of a Constitutive Model to Jointed Block Mass

Summary ?Deformational behaviour of jointed rocks is governed both by the properties of the joints and the intact rock. The applicability of a constitutive model to predict the stress-strain behaviour of a jointed block mass is discussed in this note. The stress-strain curves predicted by the model are compared with the experimental curves from physical model tests. It is observed that the model predicts the tangent modulus of the mass very closely. The failure strains are, however, underestimated, which may be due to the constant value of shear stiffness, resulting in linear stress-strain response by the model. The non-linear response of stress-strain curves can be obtained by using a variable value of the shear stiffness of the joints. A simple approach has been presented to assign variable stiffnesses, based on direct shear test results on joints. It is found that by assuming a variable stiffness the axial strains can also be predicted accurately. The transverse strain, however, remains highly underestimated by the model.     

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

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

多节理岩体裂纹扩展及变形破坏试验研究

为了研究不同倾角下多节理岩体力学行为,采用10 MN微机控制电液伺服大型多功能动静力三轴仪,开展包含较多预制非贯通节理类岩石试件的单轴压缩试验,研究了多节理岩体裂纹的特征、贯通模式、破坏模式、应力-应变特征等与节理倾角之间的关系。试验结果表明,（1）多节理岩体的裂纹类型主要有翼裂纹和次生共面裂纹,翼裂纹的扩展路径与单个节理情况下的扩展路径差异较大,翼裂纹起裂后沿起裂方向存在较长的扩展长度,直接与相邻节理或翼裂纹形成贯通,并且裂纹的贯通表现出四种不同的模式;（2）多节理岩体的破坏模式归纳为3种类型：平面破坏、块体转动式破坏和台阶式破坏;（3）根据多节理岩体的应力-应变曲线在应变软化阶段所表现出的不同非线性变形行为特征,可以将曲线归纳为4种类型;（4）多节理岩体的强度和变形各向异性特征非常显著,强度和弹性模量均在节理倾角30°时最小,90°时最大。     

Continuum modelling of jointed porous rock

Constitutive equations for the mechanical and hydraulic behaviour of saturated porous rock with joint sets of specified orientations are developed by superposing continuum representations for the mechanical and hydraulic properties of the intact rock and each of the joint sets. The resulting continuum theory allows for fluid diffusion through and between interconnected rock pores and joint sets of specified orientation, and also accounts for the anisotropy of the mechanical properties due to joint stiffnesses. The accuracy and reliability of this model are verified by finite element simulation of example problems. The first example considers joint orientation-dependent rock deformation in a hypothetical porous medium with one joint set of different dip angles. More realistic examples related to rock slope stability and reservoir-induced seismicity are also considered in which the constitutive law's utility for modelling time-dependent fluid pressures is illustrated.     

节理化岩体复合式滑移破坏的模型试验研究

针对节理化岩体的自身特点,提出沿岩体层面产生的应力集中以及层面与节理之间的残余完整岩桥破坏,将可能导致岩体产生复合式滑移破坏。并在室内利用岩石三轴刚性伺服机,通过模型试验对岩体的层面和一组与层面正交的节理的情况进行模拟研究,分析了结构面性状和围压对复合式滑移破坏的影响,给出了在贯通性结构面(层面)倾角、非贯通性结构面以及围压等因素共同作用下,节理化岩体发生复合式滑移破坏的条件以及破坏面特征。在与前人的试验资料进行对比分析的基础上得到可靠的结论,期望为在将来的实践工作中正确评估节理化岩体的破坏模式提供指导性意见。     

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

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

A three-dimensional equivalent continuum constitutive model for jointed rock masses containing up to three random joint sets

A three dimensional constitutive model is formulated for deformation analysis of jointed rock masses containing up to three joint sets with arbitrary spatial configurations. A representative elementary volume (REV) that represents the deformational response of the rock mass is defined and the constitutive relationships are developed based on superposition of deformations of the REV components. By representing the constitutive relationships in a tensorial form, the model is able to implement deformation anisotropy of jointed rock masses. The Mohr-Coulomb failure criterion with tension cut-off is used for the intact rock and the joint sets. The model is implemented in FLAC^3D and the deformations and strength values calculated by the model are compared with the results from a 3DEC model and analytical solutions. The model results are in good agreement with those obtained from 3DEC.     

Elastic stress and strain in jointed rock masses by means of crack tensor analysis

Summary An elastic stress-strain relation is formulated in terms of crack tensors which makes it possible to take into account explicitly the effect of joints on elastic behavior of rock masses. The present study is to discuss some related topics which may be encountered in its practical application. Two problems are solved by incorporating the elastic stress-strain relation into a program for three-dimensional finite element analyses; i. e., stress concentration by surface loading and displacement by excavation of an intersecting tunnel. Validity of the results is checked by comparing them with a laboratory model test and a field test, with the following conclusions: The overall distribution of stress definitely depends on a joint stiffness ratio (i. e., normal stiffness to shear stiffness). If the ratio is chosen as unity, the stress concentration occurs mainly in the direction parallel to major joints. If the ratio is high, say 10, then the stress concentrates along the perpendicular as well as the parallel directions to major joints. It can be said, on the basis of the fairly good agreement of the calculations using the high stiffness ratio with the field and laboratory measurements, that the elastic solution by crack tensors provides a practical tool for estimating the stress and strain in strongly jointed rock masses.     

节理化岩体复合式滑移破坏的模型试验研究

针对节理化岩体的自身特点,提出沿岩体层面产生的应力集中以及层面与节理之间的残余完整岩桥破坏,将可能导致岩体产生复合式滑移破坏。在室内利用岩石三轴刚性伺服机,通过模型试验对岩体的层面和一组与层面正交的节理的情况进行模拟研究,分析了结构面性状和围压对复合式滑移破坏的影响,给出了在贯通性结构面（层面）倾角、非贯通性结构面的以及围压等因素共同作用下,节理化岩体发生复合式滑移破坏的条件以及破坏面特征。最后在与前人的试验资料进行对比分析的基础上得到可靠的结论,期望为在将来的实践工作中正确评估节理化岩体的破坏模式提供指导性意见。     

Analysis of Wave Propagation Through a Filled Rock Joint

An analytical and experimental study on a longitudinal wave (P-wave) transmission normally across a filled rock joint is presented in this paper. The dynamic property of the filling material for the artificial rock joints is derived from a series of modified split Hopkinson pressure bar (SHPB) tests. The incident and transmitted waves in granitic pressure bars are calculated by wave separations of the strain gauge readings. The incident wave is approximated by a series of half-sinusoidal waves, and an analytical model on wave propagation across a filled rock joint is then deduced. The derived wave transmission coefficients across the filled joint agree very well with those from the test results. Both the analytical and test results show that the wave transmission coefficients are influenced by the mechanical properties and the input energy of the incident waves. Analytical parametric studies with respect to pre-compaction of the filling material, the frequency and amplitude of the incident wave have also been conducted.     

Anisotropic Deformation of a Circular Tunnel Excavated in a Rock Mass Containing Sets of Ubiquitous Joints: Theory Analysis and Numerical Modeling

Determining anisotropic deformation surrounding underground excavations for tunnels is an intuitional task that involves many difficulties due to the inherent anisotropies in the strength and deformability of natural rocks. This study investigates joint-induced anisotropic deformation surrounding a tunnel via a numerical simulation that accounts for the mechanical behavior of intact rock, the orientations of joint sets, and the mechanical behavior of joint planes; this numerical simulation can model the complete stress–strain relationship with anisotropic rock mass characteristics. Simulation results demonstrate that the well-known excavation-induced stress variation–decrease in the radial component and increase in the tangential component–decrease shear strength and increase shear stress for the joint plane tangential to the tunnel wall, resulting in joint sliding failure and considerable shear deformation. This joint sliding failure and significant shear deformation account for the joint-induced anisotropic deformation surrounding a tunnel. When a rock mass has two joint sets with unfavorable joint orientations, the area with joint sliding failure can deteriorate mutually, resulting in large anisotropic deformation. Additionally, for a rock mass containing three joint sets with well-distributed orientations, joint sliding in various joint sets and associated stress variations can counter balance each other, resulting in less anisotropic deformation than those of rock masses containing one or two joint sets.     

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.     

Numerical investigation of the effect of joint geometrical parameters on the mechanical properties of a non-persistent jointed rock mass under uniaxial compression

The strength and deformability of rock masses transected by non-persistent joints are controlled by complex interactions of joints and intact rock bridges. The emergence of synthetic rock mass (SRM) numerical modelling offers a promising approach to the analysis of rock masses, but has not been rigorously compared with actual physical experiments. In this work, SRM modelling by the discrete element software PFC3D is used to investigate the effect of geometric parameters of joints on the rock mass failure mechanism, unconfined compressive strength and deformation modulus. Firstly, a validation study is undertaken to investigate the ability of SRM modelling to reproduce rock mass failure modes and strength as determined by uniaxial and biaxial compression testing in the laboratory. The numerical analyses agree well with physical experimentation at low confining pressure. A sensitivity study is then undertaken of the effect of joint configuration parameters on the failure mode, unconfined compressive strength and deformation modulus of the rock mass. Five failure modes are predicted to occur: intact rock, planar, block rotation, step-path and semi-block generation. It is found that the failure mode is determined principally by joint orientation and step angle and the joint orientation with respect to principal stress direction is the parameter with the greatest influence on rock mass properties.     

双轴应力下非贯通节理岩体压缩损伤本构模型

针对目前节理岩体损伤变量定义中大多仅考虑节理长度、倾角等几何性质,而未考虑节理抗剪强度等力学性质的不足,基于断裂力学中的由于单个节理存在引起的附加应变能增量与损伤力学中的损伤应变能释放量相关联的观点,推导出了在双轴应力下含单条非贯通闭合节理岩体的损伤变量计算公式,并根据断裂力学理论对双轴压缩荷载下的单个节理尖端应力强度因子计算方法进行了研究,得出了应力强度因子KⅠ、KⅡ的计算公式;同时考虑多节理间的相互作用,给出了单组单排及单组多排非贯通节理尖端应力强度因子计算公式,由此建立了相应的节理岩体双轴压缩损伤本构模型,并利用该模型进行了相应的算例分析。结果表明：对含单条非贯通闭合节理的岩体而言,当节理倾角小于其内摩擦角时,岩体强度与完整岩石相同,岩体损伤为0,而后随着节理倾角增加,岩体强度、损伤随节理倾角的变化分别呈开口向上及向下的抛物线,当节理倾角约为60°时,岩体损伤最大,强度最低。随着节理长度增加,岩体损伤增加,而随着节理内摩擦角的增加,岩体损伤则减小;对含单组单排非贯通闭合节理的岩体而言,当节理总长度一定时,随着单条节理长度的减小及节理条数的增加,岩体损伤则逐渐减小,但其减小幅度与节理条数并不呈线性关系。     

裂隙网络模拟与REV尺度研究

研究了利用平面四边形模拟节理岩体三维裂隙网络的方法。在生成裂隙网络时,同时考虑结构面几何参数和力学参数的随机性。利用裂隙网络研究了确定岩体REV尺度的指标。在此基础上,编制了岩体裂隙网络模拟程序、裂隙网络图形输出程序和岩体REV尺度指标分析程序。通过算例分析,验证了研究成果及程序的合理性,并得出了岩体REV尺度约为各组裂隙中最大迹长期望值3~4倍的结论。研究成果为后续计算岩体等效力学参数奠定了基础。