Three‐phase medium model for filled rock joint and interaction with stress waves

A three‐phase medium model is proposed in describing the dynamic property of filled rock joints and an analytical study on longitudinal wave transmission normally across a three‐phase rock joint is presented. Parameters in the three‐phase medium model were determined by a series of modified split Hopkinson pressure bar (SHPB) tests, where a sand or clay layer was used to represent an artificially filled rock joint. The effect of the unloading path on the transmitted wave was discussed by comparing the analytical and SHPB test results. The derived wave transmission coefficients across the filled joint agreed very well with those from the test results. Both the analytical and the test results showed that the wave transmission coefficients were affected by the mechanical properties of the fillings. Parametric studies with respect to the volume ratios of water and air in the three‐phase medium and the type of filling material have also been performed. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Numerical study on S-wave transmission across a rough,filled discontinuity

This paper presents a numerical simulation of S-wave propagation across a rough, filled discontinuity using the universal distinct element code (UDEC). The ability of UDEC to simulate a stress wave across a smooth and planar discontinuity filled with an elastic material is validated through comparisons with analytical solutions. Next, the effect of the plastic deformation of the fill on the wave propagation is investigated. The model is extended to further study S-wave propagation across a filled discontinuity with rough interfaces, which is described using a sawtooth. The transmission coefficient defined by the energy is used to measure the wave attenuation. Finally, a parametric study is conducted to investigate the influences of the filled thickness, asperity angle, and incident amplitude on the transmission waves and transmission coefficients. The asperity angle and filled thickness together determine the transmitted waveform and transmission coefficient. The transmitted wave may be cut off when the incident wave amplitude exceeds a threshold value. The transmission coefficient decreases with a different trend with the incident wave amplitude increasing when the asperity angle varies. Compared with planar discontinuity, a filled discontinuity with rough interfaces is more sensitive to the amplitude of the incident wave. The causes of these phenomena are analyzed in detail. In addition, the deformation of the fill material is strongly related to the wave attenuation.     

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

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

Analysis of Stochastic Seismic Wave Interaction with a Slippery Rock Fault

Stochastic seismic wave interaction with a slippery rock fault is studied, based on the principle of conservation of momentum at the wave fronts along the fault. By using the displacement discontinuity method, the wave propagation equations are derived for incident longitudinal-(P-) and shear-(S-) waves, respectively. This is an extension of the study by Li and Ma (2010) for blast-induced wave propagation across a linear rock joint. Stochastic seismic waves are generated from a frequency spectrum and used to analyze the seismic wave interaction with a rock fault having a Coulomb-slip behavior. Parametric studies are carried out to investigate the effect of the intensity and impinging angle of the incident seismic waves on wave propagation across a slippery rock fault. Results show that the transmission of the incident P-wave is almost not affected by the fault, on the contrary, this is not the case for an incident S-wave, due to the occurrence of a relative slip which is related to the impinging angle of the incident S-wave. A quantitative study is presented which is of help in understanding the propagation and attenuation laws of seismic waves in discontinuous rock masses.     

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.     

位移不连续模型在黏弹性节理刚度计算的应用

考虑节理质量对应力波传播的影响,运用波的位移位函数推导了谐波在厚黏弹性节理的透、反射系数计算公式,采用波形相关系数描述子波穿过黏弹性节理的波形变化,讨论具有一定厚度的黏弹性节理简化位移不连续模型的适用条件。设厚黏弹性节理模型和位移不连续模型的透射波波形相关系数为0.9时对应的节理厚度为临界厚度,岩体与节理的阻抗比对临界厚度的影响很小;临界厚度随子波中心频率增大呈负指数减小,入射角越大,临界厚度随中心频率减小得越慢。试验数据分析表明：当节理厚度为0.03 m时,采用位移不连续模型和厚黏弹性模型计算得到的节理力学参数非常接近,随节理厚度和子波中心频率增加,运用位移不连续模型的计算结果偏差越大,试验结果与理论分析是一致的。     

Investigation into the effect of fault properties on wave transmission

The estimation of wave transmission across the fractured rock masses is of great importance for rock engineers to assess the stability of rock slopes in open pit mines. Presence of fault, as a major discontinuity, in the jointed rock mass can significantly impact on the peak particle velocity and transmission of blast waves, particularly where a fault contains a thick infilling with weak mechanical properties. This paper aims to study the effect of fault properties on transmission of blasting waves using the distinct element method. First, a validation study was carried out on the wave transmission across a single joint and different rock mediums through undertaking a comparative study against analytical models. Then, the transmission of blast wave across a fault with thick infilling in the Golgohar iron mine, Iran, was numerically studied, and the results were compared with the field measurements. The blast wave was numerically simulated using a hybrid finite element and finite difference code which then the outcome was used as the input for the distinct element method analysis. The measured uplift of hanging wall, as a result of wave transmission across the fault, in the numerical model agrees well with the recorded field measurement. Finally, the validated numerical model was used to study the effect of fault properties on wave transmission. It was found that the fault inclination angle is the most effective parameter on the peak particle velocity and uplift. Copyright © 2017 John Wiley & Sons, Ltd.     

Experimental Study of Ultrasonic Waves Propagating Through a Rock Mass with a Single Joint and Multiple Parallel Joints

Experiments were conducted to study the relationship between the transmission ratio (TR) and normal stress, joint roughness, joint number and frequency of incident waves, respectively, when ultrasonic waves pass across a rock mass with one joint and multiple parallel joints oriented normally. The ultrasonic waves were generated and received by pairs of piezoelectric transducers and recorded by an ultrasonic detector. The specimens were subjected to normal stress by a hydraulic jack and loading frame. The jointed rock mass was produced by superposing rock blocks in the study. Rough joints were produced by grooving notches on the planar joints formed by sawing directly. In the case of multiple parallel joints, the overall thickness of specimens was maintained while the joint number changed. Three pairs of P-wave transducers and one pair of S-wave transducers with different frequencies were, respectively, applied and all transducers emitted signals perpendicular to the joints in the experiment. The results indicate that TR increases with increasing normal stress while the increment rate decreases gradually. This is particularly so when the normal stress is high enough that TR will approximate 1 even if the rock mass has many joints. In addition, the experiments indicate that the higher the wave’s frequency, the lower its TR, and this phenomenon is gradually reduced as the normal stress increases. In response to S-waves, TR increases with increase in joint roughness; however, in response to P-waves, TR decreases gradually with increase in joint roughness. For multiple parallel joints in a fixed thickness rock mass with normally incident P-waves, TR does not always decrease with increase in the number of joints, and there is a threshold joint spacing for a certain incident wave: when the joint spacing is smaller than the threshold value, TR will increase with a decrease in joint spacing. The experimental results support similar conclusions based on analytical results drawn by Cai and Zhao (Int J Rock Mech Min Sci 37(4):661–682, 2000), Zhao et al. (Int J Rock Mech Min Sci 43(5):776–788, 2006b) and Zhu et al. (J Appl Geophys 73:283–288, 2011a).     

Analytical solutions for dynamic pressures of coupling fluid–porous medium–solid due to SV wave incidence

This paper presents the results of theoretical investigation on the dynamic coupling of an ideal fluid‐porous medium‐elastic half‐space system subjected to SV waves to study the effect of sediment on the seismic response of dams for reservoirs that are deposited with a significant amount of sediment after a long period of operation. The effects of the porous medium and the incident wave angle on dynamic pressures in the overlying ideal fluid are analyzed, and the reflection and transmission coefficients of the wave at the material interfaces are derived using an analytical solution in terms of displacement potentials. The numerical test of modeling shows that the dynamic pressures significantly depend on the properties of porous medium. The fully saturated porous medium reduces the response peaks slightly, while the partially saturated porous medium causes a considerable increase in the resonance peaks. Copyright © 2009 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.     

不同饱和度土层分界面上剪切波的反射与透射

基于多相孔隙介质弹性理论,给出了非饱和土中不同弹性波的传播方程。根据分界面上的边界条件,建立了各势函数波幅值之间的关系式,讨论了入射剪切波在不同饱和度土层分界面上的反射与透射问题。在无限空间非饱和土体中存在3种压缩波和1种剪切波,因此,当剪切波传播到不同饱和度的非饱和土层分界面上将分别在上、下土层激发产生4种反射波和4种透射波。推导出不同反射波和透射波的振幅比例系数和能量比例系数的理论表达式,并且在此基础上进行数值分析。在数值算例中分别研究了各反射波与透射波的能量比例系数（即能量反射率和能量透射率）受入射频率、入射角度以及上、下土层土体饱和度变化的影响情况。计算结果表明：各能量反射率和能量透射率不仅与入射角和入射频率有关,而且其受上、下土层饱和度变化的影响也同样不能忽视。     

酸性干湿循环下充填节理岩石劣化性能试验研究

为探究酸性干湿循环作用对充填节理岩石的损伤劣化规律,首先对预制充填节理岩石进行酸性和中性环境下不同次数的干湿循环前处理,再采用万能试验机和分离式Hopkinson杆（split Hopkinson pressure bar,简称SHPB）装置对其进行单轴压缩试验和动态冲击试验,对充填节理岩石纵波波速、静态抗压强度以及动态抗压强度的劣化效应进行分析,并进一步分析酸性和中性环境干湿循环作用下充填节理岩石动态冲击的破坏形态。结果表明,随着干湿循环次数的增加,无论是酸性环境还是中性环境,充填节理岩石的波速值、静态抗压强度以及动态抗压强度不断减小,劣化程度不断加剧,其中酸性干湿循环作用对静态抗压强度的劣化最显著;此外发现干湿循环作用下充填节理岩石静态、动态抗压强度与波速值之间均存在一致性较强的线性关系;酸性环境干湿循环作用下充填节理层粉化程度和两侧岩石破碎程度更严重。研究成果可为酸性环境干湿循环作用下工程岩体稳定性分析提供科学依据。     

应力波斜入射对岩体节理端部的影响

岩体节理在受到应力波扰动时端部受力会发生不同程度的连续性变化,为具体分析这种动态变化与应力波入射角度之间的关系,引用岩石非线性法向本构关系与线性切向本构关系的组合模型以及相应的在P波斜入射节理的应力波传播方程,结合岩石断裂力学中Ⅰ、Ⅱ型裂纹尖端应力和位移场的相关理论,得出组合形式下的节理端部应力场和位移场随质点振速变化的计算公式。通过不同入射角的模拟脉冲信号作用对端部应力位移的变化分析：应力波斜入射节理导致节理端部上、下两侧应力与位移场非对称分布,随着入射波质点振速的增减变化,应力集中位置会出现变化;计算并整理模型中节理端部上、下两侧0.005 m位置的数据,节理法向刚度由入射波质点振速带来的变化直接影响到应力波的透射与反射,进而会导致端部的应力与位移出现滞后效应;节理端部横向位移值与入射角度并非单调变化,而竖向位移会随着入射角的增大呈现下降趋势。     

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

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

饱和土半空间中圆柱形孔洞对平面P波的散射

在 Biot饱和多孔介质动力学理论的基础上,首次建立了求解饱和土半空间中圆柱形孔洞对平面 P 波散射问题的波函数展开法。首先。分析了具有圆柱形孔洞的饱和土半空间场地在平面 P 波入射下产生散射波系,并将入射波和散射波的波函数在圆柱坐标下展开。然后,引入边界条件,求出散射波函数的待定系数,从而,得到饱和土半空间中圆柱形孔洞对平面P波的散射问题的解析解。根据所得的波函数的解,可求解区域内的位移、应力的值,同时,分析了入射波频率、入射角对柱面上的应力集中因子的影响。     

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

岩体工程计算分析中结构面刚度系数是至关重要的力学参数,计算分析的精度和可靠程度与这个参数密不可分,结构面刚度系数取值仍然是一个难点。岩体中应力波传播至结构面处将会发生反射和透射现象,利用应力波透射系数可反演结构面动态刚度系数。本文从细观力学角度运用颗粒离散元方法,开发分段线性接触模型及应力波吸收边界模型,开展宏观岩体中应力波传播的模拟,结合准静态压缩试验模拟,研究了较为平直的岩体结构面分别在不同正应力条件下的动、静态刚度系数的变化特征。模拟结果表明:(1)利用C++语言开发的分段线性接触模型很好地实现了结构面非线性变形特征的模拟;(2)基于颗粒离散元方法能够准确地反映岩体中应力波传播规律,应力波通过不同刚度结构面的透射系数与理论解一致;(3)在离散颗粒模型中加入黏滞吸收边界条件很好地实现了在有限尺寸模型中远场应力波传播模拟;(4)在岩体模型中结构面接触部位运用分段线性接触模型,通过模拟应力波传播与单轴压缩试验分别获得了一致性较好的结构面动、静态刚度系数,结构面动/静态刚度系数之比值约为1.0。本文对岩体中结构面刚度的测试和取值具有重要的指导意义。     

充填物性质影响结构面剪切特性试验研究

通过自主研发的煤岩剪切-渗流耦合试验装置,开展了无充填结构面和充填石膏、岩屑、黄泥结构面的剪切试验。利用三维光学扫描技术,研究了结构面的三维形貌特征和裂隙开度演化。结果表明：充填结构面的峰值剪应力和法向位移由充填石膏、岩屑到黄泥依次递减,受充填物强度影响较大;无充填结构面由于直接接触,剪切后表面发生明显磨损,结构面粗糙度系数（JRC）整体减小且变化最为显著,充填石膏和岩屑结构面的JRC变化次之,充填黄泥结构面几乎没有变化;充填物破坏和结构面磨损主要受裂隙开度演化和局部应力集中的影响。随着剪切位移的增大,在裂隙开度较小的区域,应力分布较大,充填物易沿此处破碎,无充填结构面沿此处严重磨损。     

不同充填度岩石分形节理抗剪强度的数值模拟

在对岩石节理进行分形模拟的基础上,建立了岩石分形充填节理直剪试验数值分析模型,采用数值模拟方法,研究了充填度与岩石分形充填节理抗剪强度之间的关系。结果表明,随着充填度的增大,分形节理抗剪强度总体上呈非线性减小。当充填度小于1时,随着充填度的增大,分形节理抗剪强度迅速减小,节理抗剪强度受节理表面形态、充填物力学性质以及壁岩力学性质的共同影响;当充填度大于1时,节理抗剪强度主要受充填物力学性质的控制,强度参数接近于节理充填物的强度参数。      

节理面透射模型及其隔振性能研究

振动波入射到节理界面上,与软弱断层一样会产生波场的分解,能够起到隔振作用。根据线性变形节理的透反射规律研究发现,隔振后各点的波场可由透射P波与透射S波叠加而成,其能量分配服从Zoeppritz方程,不连续位移造成了透射波的衰减;法向刚度增加透射纵波系数会增加,而切向刚度与透射纵波无关,法向与切向刚度对透射S波均有影响,刚度增加透射S波振幅逐渐降低,逐步逼近于0;由于透射后波场相互干涉,不同波场延时不同,透射合成波与入射波相比,振幅降低,持续时间增加,波形尾段有畸变;总能量系数与节理面的刚度无关,而随入射角增加不断降低。