考虑岩体应力的结构面中应力波传播特性分析

应力波在含结构面岩体中的传播一直是岩石动力学的一个重要课题。深部岩体应力会影响岩体结构面的力学特性,进而影响应力波的传播规律。基于时域递归分析方法,首先对考虑岩体应力情况下应力波与非线性岩体结构面的相互作用进行了分析,推导出了波的传播方程。研究中假定岩体结构面法向满足非线性双曲变形模型（Barton-Bandis模型,简称B-B模型）,切向满足库仑滑移模型。其次,对一些特殊情况,如结构面水平或垂直分布、波阵面水平分布等,进行了计算模型和波的传播方程的简化分析。通过对比,发现岩体应力对应力波的传播特性有明显的影响。最后,参数研究分析了岩体应力特性、节理特性及入射波频率等对透、反射系数的影响规律。     

基于离散颗粒模型的岩体结构面动-静态刚度系数对比的数值模拟研究

岩体工程计算分析中结构面刚度系数是至关重要的力学参数,计算分析的精度和可靠程度与这个参数密不可分,结构面刚度系数取值仍然是一个难点.岩体中应力波传播至结构面处将会发生反射和透射现象,利用应力波透射系数可反演结构面动态刚度系数.本文从细观力学角度运用颗粒离散元方法,开发分段线性接触模型及应力波吸收边界模型,开展宏观岩体中...     

水平层状岩体边坡动力响应中的结构面效应研究

根据应力波传播原理分析了水平层状岩体边坡中应力波传播特征,建立了应力波在该类边坡中传播的模型.利用离散元软件UDEC分析了不同频率垂向压缩应力波作用下边坡动力响应规律中的结构面效应.结果表明:边坡中的水平层面对坡顶的动力响应有明显影响.低频应力波作用下,水平层状岩体边坡坡顶的垂向峰值速度较均质坡体相同部位的峰值速度的增...     

岩体中应力波幅值随时空衰减的关系

把岩石抽象为Kelvin-Voigt黏弹性体,应用复数理论建立了应力波在岩石中传播时的振动频率和波动频率的关系,在此基础上,研究了岩石和含有一组平行结构面的岩体中应力波幅值随时间和随空间衰减的关系。研究表明,应力波在岩石的传播过程中,频率是影响应力波幅值衰减的的主要因素;应力波在含有1组平行结构面岩体的传播过程中,由于结构面的存在,影响了应力波能量的传播,其幅值随时空衰减系数主要由结构面间距及其反射系数决定,结构面间距减小和反射系数增加,幅值随传播距离的衰减系数增加和随时间的衰减系数减小,而结构面的黏性系数增加,幅值随传播距离的衰减系数增加,但随时间的衰减系数变化很小。上述结论对岩体结构探测具有重要的指导意义。     

块系岩体非协调动力响应特征试验研究

岩体在块系构造状态下的动力传播存在块体和块体间软弱介质的非协调变形现象，对岩体结构的动力稳定产生重要影响。目前块系岩体的非协调动力传播响应研究尚不充分，导致不能进行有效的动力响应特征识别。通过试验对块系岩体非协调动力传播过程岩块动力响应特征进行研究，分别从动力传播速度、岩块加速度振幅衰减、岩块振动位移响应、岩块振动持续时间、块体动能显现、块体振动频域响应特性等方面对比纵波传播，给出块系岩体非协调动力传播的特征分析。通过分析可知，块系介质非协调变形运动产生的动力传播现象相比于纵波传播是一种低频低速波，并伴随着岩块的大位移运动，同时岩块振动持续时间较长，蕴含着动能与势能的相互转换过程，这些现象构成了块系岩体非协调动力传播的重要特征。该研究将为识别岩体非协调动力传播特征及块系岩体的冲击预警波动特征识别提供参考。     

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

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

应力腐蚀和压力溶解引起的裂隙开度-刚度变化对THMM耦合过程影响的数值模拟

在使用Yasuhara等建立的裂隙开度的应力腐蚀和压力溶解模型的基础上,考虑裂隙闭合量对裂隙刚度的影响,针对一个假设的位于非饱和双重孔隙-裂隙岩体中且有核素泄漏的高放废物地质处置模型,拟定两种计算工况：（1）裂隙的刚度系数是裂隙闭合量的线性函数;（2）裂隙刚度为常数,进行热-水-应力-迁移耦合的二维有限元数值模拟,考察了岩体中的温度、裂隙开度的闭合速率、闭合量、孔（裂）隙水压力、地下水流速、核素浓度、裂隙刚度和正应力的变化、分布等情况。结果表明,两种工况岩体中的温度场、孔（裂）隙水压力、地下水流速、核素浓度无明显差别;裂隙闭合基本由应力腐蚀产生;在相同计算时间内两种工况的裂隙闭合量较为接近,工况1略大;工况1中离玻璃固化体越近,裂隙刚度值越高,并且在玻璃固化体附近的应力值较大,且集中程度较高。     

岩体夹层应力波能量演化及应力响应特征分析

夹层是常见的地质结构,在地震或爆破荷载作用下,夹层对应力波的传播以及对岩体的响应具有重要的影响。以往对应力波在夹层介质中传播的研究多集中于夹层对应力波的隔振或透射性能,而对应力波在夹层中的多次折、反射过程中能量的演化规律缺乏讨论,对夹层介质的应力响应与破坏没有开展较好地分析。通过理论方法对应力波在夹层内部传播过程中的能量系数的变化规律进行了研究,分析了岩体介质波阻抗和应力波入射角对夹层内外介质中累积波动能量系数的影响规律,以及平面型边坡中软弱夹层的应力响应特征和动态稳定性。研究发现,应力波在夹层内部往复反射过程中,夹层内剩余应力波能量随折、反射发生次数呈指数曲线下降,第4次折、反射后产生的应力波能量可以忽略;夹层内外介质中应力波的累积能量系数的差异随着夹层内外介质波阻抗的相对差异的增大而增大。在平面谐波入射下,边坡内部的夹层中的剪应力和抗剪强度呈波动变化;相对P波,SV波入射会产生较高水平的剪应力,对边坡稳定性影响最大。且SV波入射时,边坡的安全系数对夹层的倾角变化更为敏感,随倾角增大而迅速降低。     

应力波与可滑移岩石界面间的相互作用研究

研究了应力波与用节理刚度系数描述的线性滑移岩石界面间的相互作用，给出了其透、反射波幅的一般表达式。通过应力波垂直人射界面时的具体分析，反映出岩石界面的节理刚度对应力波传播的重要影响，并揭示了岩石界面对应力波的高频滤波作用。     

结构面对程潮铁矿西区地表变形的影响分析

以程潮铁矿西区为例,通过对现场结构面和裂缝分布特征的调查,以及结合地表变形监测数据的分析,揭示了矿区结构面对地表变形的影响。研究结果表明,在矿区较大的水平应力条件下,岩体结构面改变了地表拉伸变形分布和破坏形态,加剧了岩体变形,在ⅰ区(剖面Ⅲ以东区域),目前最外侧裂缝以内的岩体沿着NNW结构面发生倾倒滑移破坏,地表变形以快速变形为主,而最外侧裂缝以外的岩体沿着NNW结构面发生弯曲变形,地表变形以线性稳定增长为主,进入了倾倒破坏阶段,该区域的岩体主要发生水平位移;在ⅱ区(剖面Ⅲ以西区域),岩体在开采沉陷引起的南北向应力作用下,分离成平行的块体,地表变形以快速变形为主,在局部区域,产生的一部分平行块体在近东西向的应力作用下,沿着NNW结构面产生倾倒破坏。所得成果对类似金属矿山的地表征地及安全高效生产具有实际指导意义。     

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

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

Stress‐dependent elastic wave velocity of microfractured sandstone

A model for the stress‐dependent elastic wave velocity response of fractured rock mass is proposed based on experimental evidence of stress‐dependent fracture normal and shear stiffness. Previously proposed models and previous experimental studies on stress‐dependent fracture stiffness have been reviewed to provide a basis for the new model. Most of the existing stress‐dependent elastic wave velocity models are empirical, with model parameters that do not have clear physical meanings. To propose the new model, the rock mass is assumed to have randomly oriented microscopic fractures. In addition, the characteristic length of microfractures is assumed to be sufficiently short compared to the rock mass dimensions. The macroscopic stress‐dependent elastic wave velocity response is assumed to be attributed to the stress dependency of fracture stiffness. The stress‐dependent fracture normal stiffness is defined as a generalized power law function of effective normal stress, which is a modification of the Goodman's model. On the other hand, the stress dependency of fracture shear stiffness is modeled as a linear function of normal stress based on experimental data. Ultrasonic wave velocity responses of a dry core sample of Berea sandstone were tested at effective stresses ranging from 2 to 55 MPa. Visual observation of thin sections obtained from the Berea sandstone confirms that the assumptions made for microstructure of rock mass model are appropriate. It is shown that the model can describe the stress‐dependent ultrasonic wave velocity responses of dry Berea sandstone with a set of reasonable material parameter values. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

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.     

Validation study of the distinct lattice spring model (DLSM) on P-wave propagation across multiple parallel joints

A validation study of the distinct lattice spring model (DLSM) for wave propagation problems is performed. DLSM is a microstructure-based numerical model, which is meshless and has advantages in modelling dynamic problems where stress wave propagation is important. To verify the applicability of DLSM to modelling wave propagation through a discontinuous medium, the virtual wave source (VWS) method is used to obtain analytical solutions for wave propagation across a jointed rock mass. Numerical modelling results of the commercial code UDEC are selected as the reference. The effects of particle size and lattice rotation angle on wave propagation are first studied. Then, the results of wave transmission across a single joint with a different joint stiffness and across multiple parallel joints with different joint spacings are derived with DLSM, UDEC and VWS. These results are in good agreement with each other. Therefore, the capability of DLSM to model P-wave propagation across jointed rock mass is verified, which provides confidence for the further application of DLSM to modelling more complex problems.     

Analysis of Blast Wave Interaction with a Rock Joint

The interaction between rock joints and blast waves is crucial in rock engineering when rock mass is suffered from artificial or accidental explosions, bursts or weapon attacks. Based on the conservation of momentum at the wave fronts and the displacement discontinuity method, quantitative analysis for the interaction between obliquely incident P- or S-blast wave and a linear elastic rock joint is carried out in the present study, so as to deduce a wave propagation equation. For some special cases, such as normal or tangential incidence, rigid or weak joint, the analytical solution of the stress wave interaction with a rock joint is obtained by simplifying the wave propagation equation. By verification, it is found that the transmission and reflection coefficients from the wave propagation equation agree very well with the existing results. Parametric studies are then conducted to evaluate the effects of the joint stiffness and incident waves on wave transmission and reflection. The wave propagation equation derived in the present study can be straightforwardly extended for different incident waveforms and nonlinear rock joints to calculate the transmitted and reflected waves without mathematical methods such as the Fourier and inverse Fourier transforms.