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

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

Discrete element modelling of an underground explosion in a jointed rock mass

The purpose of this paper is to investigate the application of the discrete element program UDEC for modelling underground explosions in a jointed rock mass. A field underground explosion test has been conducted and UDEC and AUTODYN were utilised in a coupled manner to predict the rock mass response due to the underground explosion. The dynamic equation for the state of the rock material obtained from impact tests was incorporated into the calibration modelling. The comparison of modelling results with test results and empirical formulas shows that UDEC is capable of modelling explosion wave propagation in jointed rock mass with high reliability.     

Theoretical Methods for Wave Propagation across Jointed Rock Masses

Different methods are presently available for the analysis of wave propagation across jointed rock masses with the consideration of multiple wave reflections between joints. These methods can be divided into two categories. One is based on the displacement discontinuity model for representing rock joints, where the displacements across a joint are discontinuous and the tractions are continuous, and the other is the equivalent medium method. For the first category, there are three methods, i.e., method of characteristics (MC), scattering matrix method (SMM) and virtual wave source method (VWS). MC solves the equation of motion by using the theory of characteristic curves. SMM is based on the definition of the scattering matrix in which the reflection and transmission coefficients of a set of joints are stored. VWS method replaces the joints in the rock mass with a virtual concept. For the second category, equivalent medium model treats the problem in the frame of continuum mechanics and simplifies it from an explicit wave propagation equation. The objective of this paper is to review and compare these theoretical methods. The comparison shows that the four solutions agree very well with each other. Some additional considerations about the advantages and disadvantages of these methods are also given in the paper.     

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.     

UDEC modeling of a field explosion test

A field explosion test is conducted to investigate the blast wave characteristics and propagation in jointed rock mass. A hybrid AUTODYN and UDEC scheme is used to predict the rock mass response due to the explosion. The original prediction is then modified by involving the actual test conditions and dynamic equation of state of the rock material. Correlating the prediction results with the field explosion test results shows that the hybrid AUTODYN and UDEC scheme is capable of modeling explosion problems with high reliability. It is observed that the modified prediction result agrees better with the field test result than the original prediction result, which implies that involving dynamic parameters in the modeling is very important to produce reliable prediction results.     

Parallelization of the distinct lattice spring model

The distinct lattice spring model (DLSM) is a newly developed numerical tool for modeling rock dynamics problems, i.e. dynamic failure and wave propagation. In this paper, parallelization of DLSM is presented. With the development of parallel computing technologies in both hardware and software, parallelization of a code is becoming easier than before. There are many available choices now. In this paper, Open Multi‐Processing (OpenMP) with multicore personal computer (PC) and message passing interface (MPI) with cluster are selected as the environments to parallelize DLSM. Performances of these parallel DLSM codes are tested on different computers. It is found that the parallel DLSM code with OpenMP can reach a maximum speed‐up of 4.68× on a quad‐core PC. The parallel DLSM code with MPI can achieve a speed‐up of 40.886× when 256 CPUs are used on a cluster. At the end of this paper, a high‐resolution model with four million particles, which is too big to handle by the serial code, is simulated by using the parallel DLSM code on a cluster. It is concluded that the parallelization of DLSM is successful. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

动载下节理岩体破坏过程的数值试验研究

节理岩体是工程中最常见的一类岩体,其在地震、爆炸等动载下的力学响应及破坏过程对相关工程安全性的影响至关重要。采用基于有限元应力分析和统计损伤理论开发的动态版RFPA2D数值模拟软件,对动载下节理岩体的动态破坏过程进行了模拟,重点讨论了节理条数、节理贯通度、节理倾角及应力波峰值对岩体动态破坏过程的影响规律。计算结果表明,断续节理岩体动态破坏过程及破坏强度与节理构造形态、应力波峰值密切相关。相同动载下,随着节理条数的增加,岩体破坏程度以及应力波能量损失增强,但当节理条数数超过一定值后,岩体破坏程度及应力波能量损失逐渐趋于稳定;节理贯通度较小时,岩体破坏程度较低且破坏单元自上而下均匀分布。随着节理贯通度的增加,岩体破坏增强,且破坏主要出现于节理上部岩体;节理倾角较小时,节理上部岩体破坏严重,易形成次生贯通裂纹。随着节理倾角增加,破坏范围逐渐变大,不易形成次生贯通裂纹;倾角为45°～60°时,岩体破坏效果最佳;动载荷的峰值越大,试样的破坏越严重。当峰值达到一定值时,节理附近发育出多条裂隙并向上下方不断发展而导致岩体完全破坏。在不同节理贯通度工况下与岩石霍布金森压杆（SHPB）试验结果进行比较,结论吻合,证明该数值模拟的合可行性和结论的可靠性。     

Graphics processing unit based parallelization of the distinct lattice spring model

This study reports a Graphics Processing Unit (GPU)-based parallelization of the Distinct Lattice Spring Model (DLSM) for geomechanics simulation. The DLSM is a newly developed numerical model for rock dynamics problems, i.e., dynamic failure and wave propagation. Despite its applicability, one of the drawbacks of this model is the high computational load for practical simulations. To tackle this problem, a GPU with a Compute Unified Device Architecture (CUDA) is adopted to parallelize the DLSM code. The performance of the GPU DLSM code is tested on two computers equipped with modern GPU cards. The results show that significant performance improvements are gained from GPU parallelization of the DLSM code (the maximum speed up achieved was 23×).     

Considerations of the discontinuous deformation analysis on wave propagation problems

In rock engineering, the damage criteria of the rock mass under dynamic loads are generally governed by the threshold values of wave amplitudes, such as the peak particle velocity and the peak particle acceleration. Therefore, the prediction of wave attenuation across fractured rock mass is important on assessing the stability and damage of rock mass under dynamic loads. This paper aims to investigate the applications of the discontinuous deformation analysis (DDA) for modeling wave propagation problems in rock mass. Parametric studies are carried out to obtain an insight into the influencing factors on the accuracy of wave propagations, in terms of the block size, the boundary condition and the incident wave frequency. The reflected and transmitted waves from the interface between two materials are also numerically simulated. To study the tensile failure induced by the reflected wave, the spalling phenomena are modeled under various loading frequencies. The numerical results show that the DDA is capable of modeling the wave propagation in jointed rock mass with a good accuracy. Copyright © 2009 John Wiley & Sons, Ltd.     

预制节理岩体试件强度及破坏模式的试验研究

采用相似材料模型试验对不同节理倾角、节理贯通度、节理组数、载荷应变率、试件长径比、节理充填物厚度及类型等7种工况下的预制节理岩体在单轴压缩下的峰值强度及破坏模式进行了研究。结果表明：节理岩体的破坏模式及峰值强度与节理构造形态密切相关。贯通节理岩体将产生沿节理面的剪切破坏或穿切节理面破坏,且与第1种破坏模式对应的岩体峰值强度更低。非贯通节理岩体的强度介于完整岩体和贯通节理岩体之间。随着平行节理组数的增加,岩体峰值强度逐渐下降。随着载荷应变率的增加,岩体峰值强度逐渐增大,相应地试件的破坏模式也变得更加复杂。试件长径比基本没有改变其破坏模式,完整试件仍主要是以张拉破坏为主,而节理试件仍以剪切破坏为主。随着长径比增加,试件峰值强度逐渐增大。随着节理充填物厚度增加,试件峰值强度降低。不同节理填充物对试件峰值强度也有一定影响。     

Propagation Characteristics of Blasting Stress Waves in Layered and Jointed Rock Caverns

The existence of joint fissures makes explosive actions between rock masses more complex. Therefore, it is of great significance to carry out experiments studying blasting stress waves propagating in jointed rock masses. Based on the Froude Similarity principle, the geological mechanical models of intact rock masses and jointed rocks have been proposed. A blasting vibration experiment was carried out and demonstrated that the propagation of the blasting stress waves and changing structures have an important relationship. A numerical simulation of the blasting stress wave propagation law in a layered jointed rock mass was carried out. This study found that with an increase in the joint angle, the peak velocity of blasting stress wave, transmission coefficient and reflection coefficient all gradually increased, while the attenuation coefficient gradually decreased. With an increase in joint spacing, the attenuation rate of the blasting stress waves increased.     

节理对爆炸波传播影响的数值研究

采用加入无反射边界条件的DDA程序,研究了节理面对应力波传播的影响。结果表明,节理面能阻碍波的传播,有利于波的衰减,节理面越多,波的反射越强,而波的透射越弱。模拟了一个现场爆炸试验,研究爆炸产生的应力波在节理岩体中传播、衰减的规律,模拟结果与现场试验结果比较吻合。研究表明,DDA方法可以模拟节理面对应力波传播的阻碍作用,用它来模拟爆炸波在节理岩体中的传播是适用的。     

斜入射线弹性节理应力波传播特征的数值模拟

依据弹性波理论,应力波斜入射线弹性节理时会发生波场分解。根据入射、透射及反射各波形的不同到时,运用离散元软件UDEC模拟应力波在含倾斜节理岩体中的传播并计算其透射、反射系数,并分析其波型转换规律。应力波斜入射单节理时,模拟得到的透射、反射系数随节理刚度、入射角度的变化规律,与已有的理论解是吻合的。应力波斜入射一组平行节理时,随着节理间距的增大,其同类波的透射系数Tpp、Tss先增大后减小,最后趋于稳定值;节理条数越多,Tpp、Tss越小。此外,不同条数的节理,透射系数达到最大值的临界节理间距值基本一致,但趋于稳定时的节理间距值随节理条数的增加而逐渐增大。     

Long-Wavelength Elastic Wave Propagation Across Naturally Fractured Rock Masses

Geophysical site investigation techniques based on elastic waves have been widely used to characterize rock masses. However, characterizing jointed rock masses by using such techniques remains challenging because of a lack of knowledge about elastic wave propagation in multi-jointed rock masses. In this paper, the roughness of naturally fractured rock joint surfaces is estimated by using a three-dimensional (3D) image-processing technique. The classification of the joint roughness coefficient (JRC) is enhanced by introducing the scan line technique. The peak-to-valley height is selected as a key indicator for JRC classification. Long-wavelength P-wave and torsional S-wave propagation across rock masses containing naturally fractured joints are simulated through the quasi-static resonant column (QSRC) test. In general, as the JRC increases, the S-wave velocity increases within the range of stress levels considered in this paper, whereas the P-wave velocity and the damping ratio of the shear wave decrease. In particular, the two-dimensional joint specimen underestimates the S-wave velocity while overestimating the P-wave velocity. This suggests that 3D joint surfaces should be implicated to obtain the reliable elastic wave velocity in jointed rock masses. The contact characteristic and degree of roughness and waviness of the joint surface are identified as a factor influencing P-wave and S-wave propagation in multi-jointed rock masses. The results indicate a need for a better understanding of the sensitivity of contact area alterations to the elastic wave velocity induced by changes in normal stress. This paper’s framework can be a reference for future research on elastic wave propagation in naturally multi-jointed rock masses.     

岩体裂隙模拟的扩展有限元法应用研究

岩体中普遍存在着断层﹑节理和裂隙等结构面,这些结构面的存在和发展对岩体的整体强度﹑变形及稳定性有极大的影响。因此,研究岩体中原生结构面的萌生﹑发展以及贯通演化过程对评估岩体工程安全性和可靠性具有非常重要的理论与现实意义。扩展有限元法（XFEM）作为一种求解不连续问题的有效数值方法,模拟裂隙时独立于网格,因此,在模拟岩体裂隙扩展﹑水力劈裂等方面具有独特优势。针对扩展有限元法的基本理论及其在岩体裂隙扩展模拟中的应用展开了研究,建立了扩展有限元法求解岩体裂隙摩擦接触、岩体裂隙破坏等问题的数值模型,并将计算模型应用于岩质边坡稳定性分析和重力坝坝基断裂破坏等工程问题。     

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

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

Evaluation of equivalent medium methods for stress wave propagation in jointed rock mass

This paper evaluates the existing equivalent medium methods for jointed rock mass and further develops the equivalent viscoelastic medium method proposed by the authors. The advantages and limitations of different equivalences to the discontinuous rock mass are discussed. Theoretical derivation of stress wave propagation through the equivalent viscoelastic medium is carried out by adopting the Fourier transformation method, and the parameters of the equivalent viscoelastic medium method are determined analytically. The frequency dependence and the wave attenuation phenomenon can be properly described when the imaginary terms of the complex moduli of the rock mass are included. The results show that the equivalent viscoelastic medium method is able to predict the effective velocity and the stress wave transmission coefficient in a rock mass more accurately than the conventional effective elastic moduli methods. An example of the stress wave propagation through rock mass with parallel joints shows that the equivalent viscoelastic medium method is promising and worthy to be further explored for application in practical rock engineering problems. Copyright © 2011 John Wiley & Sons, Ltd.     

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

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

黏弹性节理岩体的位移不连续模型研究

选用三参数标准线性固体作为岩石本构,提出了一种考虑岩体黏弹性的位移不连续模型;根据一维黏弹性波的特征线法,推导了节理处质点速度、应力和应变递推公式。首先,基于分离式霍普金森压杆（SHPB）对砂层进行试验,得到其应力-应变关系,并换算出砂层节理的法向刚度;接着,通过一维强间断黏弹性波的波速公式、高频波衰减系数以及任一频率下的衰减系数,确定出数值算法中的三参数。最后,基于自制的摆锤装置,探讨了一维应力波在节理岩体中的传播规律,试验中以两根长1 000 mm、直径为68.50 m的岩杆作为入射和透射杆,以3 mm砂层模拟节理。试验和数值结果吻合度良好,进一步验证了该方法的可靠性。