节理岩体中应力波传播规律研究的进展

讨论了位移不连续模型与特征值法结合的方法,对垂直于单一线性变形节理、多个平行线性变形节理、单一非线性变形节理,分别采用线性变形不连续模型、双曲线变形不连续模型(BB模型)进行的理论与相应的实验研究。用离散元程序UDEC,对上述问题进行数值模拟;用三维离散元程序3DEC,对节理岩体中三维波的传播进行的数值模拟。并用UDEC与有限差分程序AUTODYN-2D耦合,模拟爆炸过程及爆炸波在节理岩体中的传播。同时提出了几个有待研究的问题。     

节理网络模拟的程序实现及其在地下洞室中的动态响应分析

随着岩体工程的日益发展,人们愈发认识到岩体裂隙发育情况对工程影响的重要性,而岩体节理裂隙网络模拟是研究该类问题的重要手段。基于节理网络模拟结果,建立岩体力学分析模型,然后开展进一步的数值计算分析,在岩体工程中具有良好的发展前景。在岩体结构统计理论的基础上,自行编制了节理岩体裂隙网络模拟分析程序,直观再现了岩体内部的节理分布情况;将程序生成的可执行数据文件与离散元软件相融合;以含节理岩体的地下洞室为例,分别计算了静力和地震荷载作用下地下洞室围岩在不同时段的应力及位移变化情况。计算方法与结果对含节理岩体地下洞室围岩变形和破坏机制以及抗震工程设计具有指导意义     

三峡永久船闸高边坡开挖三维离散元数值模拟

采用面－面接触的三维离散元刚性块体模型，从实测节现面中取出其中的三组，按照其倾向，倾角和节理间距将三峡永久船闸未开挖的区域划分为10^5个离散单元，通过施加力边界条件，给出了与实测初始地应力场接近的数值模拟结果；然后，分4步模拟了永久船闸的开挖过程。计算结果表明：开挖过程会引起节理面出现张开趋势，个别岩体还会沿着节理面滑移。岩体位移的不对称现象较为自然地说明了由节理引起的岩体各向异性特征。     

爆破荷载作用下岩体振动特征的数值模拟

根据福建牛头山水电站地基岩体爆破开挖监测,运用离散元方法模拟了节理岩体距爆源不同距离处质点的振动速度和频率的变化特征,由此确定岩体质点最大振动速度和振动主频随爆源距离的衰减规律,并得到了距爆源一定距离处质点最大振动速度和振动主频与爆破药量的关系。数值模拟与现场实测的结果表明,用离散元软件UDEC计算得到的岩体振动特征和衰减规律与现场监测结果是基本符合的,误差在工程应用的允许范围之内,因此UDEC用于对岩体动态响应的数值模拟是适合的。     

基于MatDEM的砂土侧限压缩试验离散元模拟研究

离散元法基于非连续介质力学理论,尤其适用于砂土等离散介质体的数值模拟研究。利用岩土体离散元模拟软件MatDEM的二次开发功能,研发了砂土侧限压缩试验三维离散元模拟器。对三个不同级配砂土试样进行了侧限压缩试验,并且进行了与之相对应的数值模拟,通过分析对比试验结果与数值模拟结果验证了所开发模拟器的有效性。模拟结果表明:离散元法可以很好地反应砂土压缩过程中的配位数变化;每个模拟样品中,粒径较小的单元受到较大的平均压力,导致平均位移较大;数值计算结果的主要误差是由离散元颗粒自身的泊松比引起的。研究突破了常规土力学研究方法的局限性,为今后岩土工程离散元模拟研究提供了参考。      

连续屈服节理模型对剪切地震波传播的影响

地震波在跨越岩体中节理时的透反射特性是岩石地下工程抗震分析问题中的重要基础,目前开展的研究工作中尚未有针对连续屈服（CY）非线性节理模型透反射特性的研究报导。在前人提出的非线性节理透反射系数时域递归法解答的基础上,将这一方法拓展至CY模型。并采用3DEC离散元软件,对比了CY模型的时域递归解与数值模拟结果的差异。在此基础上以时域递归解为手段,辅以离散元数值模拟,开展了CY模型的参数影响研究。最终比较了库仑（MC）模型与CY模型在正弦脉冲激励和真实地震动激励下表现的差异。结果表明：CY模型的时域递归解与3DEC数值模拟结果显示了良好的一致性,证明了时域递归解用以进行后续参数研究的合理性。法向应力、入射波幅值、节理初始刚度、节理间距等参数对CY模型具有显著影响。相比库仑模型,连续屈服模型可以更好地反映地震波穿越节理时发生的复杂力学现象,如切向刚度退化、抗剪强度劣化、法向应力依赖性、滞回现象等。研究成果可为岩石地下工程的抗震设计与分析工作提供一定参考。     

麻柳嘴滑坡成因机制离散元数值模拟

在坡体中考虑地下水的影响下 ,用离散单元法对麻柳嘴滑坡的成因机制进行数值模拟。模拟结果与实际位移量值基本一致 ,表明离散单元法适用于节理边坡岩体变形和破坏机制的研究     

节理玄武岩体变形模量的尺寸效应和各向异性

目前系统研究柱状节理玄武岩工程力学特性的成果仅见于美国汉佛德玄武岩核废料埋藏工程和中国西南白鹤滩水电工程。列举了国内外相似工程柱状节理玄武岩的原位试验成果,分析了我国不规则柱状节理玄武岩体的柱体偏转特性和多级结构面发育特征,采用可变形离散元法建立了柱状节理玄武岩体的三维离散元数值模型,通过数值模拟不同尺寸的承压板试验,探讨了尺寸效应和各向异性对试验成果的影响。数值分析表明:承压板的尺寸效应是产生柱状节理岩体水平和铅直原位试验各向异性的主要原因,随着承压板直径的增大,各向异性效应减小,当承压板直径大于6m后,柱状节理玄武岩体可视为各向同性体,研究成果对利用和改造柱状节理岩体具有较好的工程意义。     

基于离散单元法的裂隙岩体渗流应力耦合分析

文章从离散单元法的基本原理和岩体裂隙网络渗流的特点出发,以节理开度变化为纽带,对裂隙岩体渗流场与应力场耦合的离散网络模型进行分析,并用离散元数值方法结合算例分别进行了渗流场与应力场的耦合、非耦合计算。结果表明:耦合计算所得到的应力、位移和渗流量与非耦合计算得到的结果相比,都有不同程度的增大。因此,应从渗流场应力场两场耦合的角度考虑,才能提供更加符合实际的满足安全性的解答。     

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

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

节理岩体冲击破坏的数值流形方法模拟

节理的存在是岩体明显区别于岩石的一个重要特征,也是导致岩体具有非均质性及各向异性的一个重要因素。节理岩体在爆炸冲击荷载作用下的破坏特征不仅是岩体力学研究的一个热点问题,也是岩石爆破领域经常遇到的一个工程问题。目前,常用的有限元和离散元等数值方法都不能很客观地反映岩体中节理分布的实际规律,而最新出现的数值流形方法则完全能够按照岩体中节理的实际分布规律进行建模,并能够很好地模拟岩体在外力作用下新裂纹的产生及已有节理的开裂。对3种不同节理岩体,即无节理岩体、均布水平节理岩体、均布垂直节理岩体等在均布垂直于圆周方向上的冲击荷载作用下的破坏过程进行了模拟。由模拟结果可以明显看出,节理的存在对岩体的破坏型式起到了关键性的控制作用。对于完整岩体,其破坏型式呈现出了很好的对称性;而对于水平及垂直节理岩体,其破坏型式在很大程度上受到节理分布的控制。研究结果对岩石爆破工程设计具有一定的指导意义。     

基于复杂结构精细描述的岩质高边坡稳定性分析

用数值方法对岩质高边坡进行稳定性分析时,描述岩体结构对的精细程度会影响分析结果,但常见的有限单元法程序仍难以对复杂节理岩体进行精细建模。为解决这一问题,将结构面网络模拟与离散单元法相结合,在UDEC软件中,利用FISH语言编写网络模拟程序,依据结构面统计资料和结构面分级,实现对复杂岩体结构的精细描述。以某大型水电工程边坡为例,在对岩体结构进行精细描述的基础上,采用离散元强度折减法对边坡进行稳定性分析。通过与极限平衡法和一般离散元结果的对比,表明基于复杂岩体结构精细描述确定的边坡潜在滑动面和安全系数是合理的,为复杂岩质边坡破坏模式和稳定性的分析提供了新的途径。     

Analysis of Critical Excavation Depth for a Jointed Rock Slope Using a Face-to-Face Discrete Element Method

Summary The critical excavation depth of a jointed rock slope is an important problem in rock engineering. This paper studies the critical excavation depth for two idealized jointed rock slopes by employing a face-to-face discrete element method (DEM). The DEM is based on the discontinuity analysis which can consider anisotropic and discontinuous deformations due to joints and their orientations. It uses four lump-points at each surface of rock blocks to describe their interactions. The relationship between the critical excavation depth D _s and the natural slope angle α, the joint inclination angle θ as well as the strength parameters of the joints c _r ,φ _r is analyzed, and the critical excavation depth obtained with this DEM and the limit equilibrium method (LEM) is compared. Furthermore, effects of joints on the failure modes are compared between DEM simulations and experimental observations. It is found that the DEM predicts a lower critical excavation depth than the LEM if the joint structures in the rock mass are not ignored.     

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

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

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.     

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.     

离散元法及其在岩土工程中的应用综述

离散元法是基于不连续性假设的数值方法,它特别适合于求解节理岩体中的非连续性问题。在介绍离散元法基本原理的基础上,着重对离散元法在岩土工程领域的应用现状作了叙述和分析,并对其发展趋势进行了探讨。     

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.     

On the role of pre-existing discontinuities on the micromechanical behavior of confined rock samples: a numerical study

The deformation process and failure mechanism of rock mass with increased density of initial joints subjected to confined stress state are investigated in this study using discrete element method (DEM). A numerical model of standard size granite samples is developed and validated using experimental data for both intact and jointed rocks. The micro-parameters of the rock material are first determined, and the effects of the rock discontinuity on strength, deformability, stress–strain relationship, and failure modes are then investigated at the macro-scale level. Analyses are also performed to examine the tensile and shear crack distributions, fragmentation characteristics, particle kinematics, and energy dissipation to advance the current understanding of the deformation processes and failure mechanisms of jointed rock masses. The microscopic evolutions in the fabric and force anisotropy during loading and distributions of contact forces provide insights into the influence of increasing initial jointing on the macroscopic deformational behavior of the rock. The results show how the deceleration in the growth of fabric and contact force anisotropies develops and confirms that the increase in initial jointing and the associated changes in microstructure can restrain the development of anisotropy, thereby reducing significantly the strength of the rock samples.

     

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.