基于Hoek-Brown准则的节理岩体能量参数估算

岩体变形和破坏可以看作能量耗散和释放的过程。由于节理岩体结构复杂且难以开展室内试验,因此无法通过试验直接求取其能量学参数。基于Hoek-Brown准则和岩石能量理论,提出了节理岩体在临界状态能量学参数的估算方法。针对含贯通节理(或层面)岩体,通过修正岩块单轴抗压强度以体现贯通节理的方向效应。采用PFC~(3D)分别模拟小尺寸岩样(Φ50 mm×100 mm)和大尺寸岩体(Φ2 m×2 m)的三轴压缩试验,通过岩石三轴试验结果拟合岩石数值模拟的细观参数并应用于节理岩体的模拟。根据节理岩体模拟得到应力应变曲线和能量流,验证了Hoek-Brown准则对节理岩体能量参数估算的合理性。     

非贯通节理岩体损伤变量计算方法研究

作为岩体组成部分的非贯通节理对岩体力学特性有着重要影响,然而几乎目前所有的岩体损伤变量计算方法都仅考虑节理几何参数对岩体力学特性的影响。对含单组非贯通节理的岩体力学特性进行研究,提出一个能够同时考虑节理几何及强度参数对岩体力学特性影响的新的岩体损伤变量计算方法。首先,采用弹性余能等效假设代替Lemaitre应变等效假设研究由节理引起的岩体各向异性损伤,并基于断裂力学中单个节理引起的附加应变能增量与损伤力学应变能释放量相关联的观点,推导出由单条节理引起的损伤变量计算公式。其次,根据断裂力学理论获得了单轴压缩下单条节理尖端应力强度因子（SIF）KⅠ、KⅡ的计算公式。最后,通过考虑节理间的相互作用给出了单组单排或多排节理尖端应力强度因子KⅠ、KⅡ的计算公式,得到了单组节理引起的岩体损伤变量计算公式,并与已有试验结果的对比分析证明了该公式的合理性。     

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

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

采用Laplace变换分析桩基应力波传播问题

将土体中的混凝土桩基简化为黏滞介质中的一维黏弹性杆,桩顶部受锤头冲击产生内部应力波。根据杆内微元应力平衡条件建立杆中一维黏弹性应力波传播的控制方程,结合桩顶锤头冲击条件和桩底弹-黏性约束条件给出桩基两端的耦合边界条件。对控制方程和定解条件作Laplace变换并求解变换后的常微分方程,得到变换域的应力像函数解析解。采用数值反变换技术将像函数转变为时间域的应力波形。应用此方法可以较方便地分析桩基中应力波的产生、传播、反射及相互作用过程。     

基于岩体三轴压缩试验的节理力学参数确定方法

岩体节理面的刚度系数、黏结力和内摩擦角是反映节理面力学特性的重要参数。基于含不同倾角节理岩样的三轴压缩试验,通过对标准岩样的应力-位移分析,由互等功原理,建立节理岩样应变与位移新的几何关系,推导出弹性状态下节理刚度系数解析解,根据摩尔-库仑强度理论得出节理面的黏结力和内摩擦角。为了检验方法的正确性,对含节理面的同类岩体进行直接剪切试验,测出节理面的法向刚度、切向刚度、黏结力和内摩擦角。结果表明,两种方法所得结果基本吻合,说明所提出的方法可以通过节理岩样三轴压缩试验及理论分析间接得到岩体节理面的刚度系数、黏结力和内摩擦角,实现了节理岩体力学参数的简捷、快速测定。     

基于复合损伤的节理岩体动态本构模型研究

基于运用霍普金森压杆（SHPB）装置对节理岩体动载试验得出的数据,从节理面倾角、贯通度、厚度、组数、填充物及应变率等不同方面分析各因素对节理岩体力学特性的影响,通过对节理岩体在高应变率下的损伤机制和破坏形式进行分析的基础上,基于复合损伤理论,对广义Bingham模型进行改进,构造了节理岩体材料在不同应变率下动态响应的本构模型,对模型计算和试验结果的比较表明,该模型能够很好地描述节理岩体动荷载下初始弹性变形阶段、稳态塑性变形阶段和加速变形破坏阶段的应力-应变关系,并且理论与试验结果吻合较好,从而证明了该模型的正确性和合理性。     

一种模拟岩石蠕变的数值流形方法

为了保障岩土工程结构能长期正常使用,需要对其蠕变变形进行分析。“时步-初应变”法是一种常用的计算岩石蠕变的方法。数值流形方法是一种新兴的数值计算方法,常用于计算节理岩体的变形,但尚未被试用于计算蠕变变形。在原数值流形方法的程序中增加了基于“时步-初应变”法的计算模块,通过对广义开尔文模型进行的模拟,显示新程序可以正确反映岩石的黏弹性蠕变趋势,并能够计算包含节理的岩体的蠕变变形,改进后的数值流形方法不但能够模拟岩石的线弹性变形,而且可以模拟岩石的黏弹性蠕变,比原流形方法更能全面地模拟岩石的变形,扩展了数值流形方法在岩土工程中的使用范围。     

三向应力作用下节理岩体渗流-应力耦合特性

节理岩体渗流-应力耦合特性研究是岩石力学领域的重点和难点。水利水电工程与地下工程中天然节理岩体往往处于三向不等压应力状态,针对三向应力状态下节理岩体渗流-应力耦合特性研究较少。其中,单节理岩体渗流-应力耦合特性研究是基础。首先分析与探究了单节理岩体侧向应力作用效果及渗流作用机制,在单节理岩体渗透性-应力耦合分析模型基础上推导了单节理岩体渗流-应力耦合的理论公式。通过室内真三轴渗流-应力耦合试验,基于全局优化算法的非线性拟合及块体离散元数值仿真,验证了所提出理论公式的合理性。基于电路原理,给出了模拟电路与渗流问题的相似物理量,通过引入等效渗流阻方法,探究了典型的含正六边形主干节理网络的渗流-应力耦合特性。该研究结果对实际工程岩体的渗流分析有一定的参考价值。     

共面闭合断续节理岩体直剪强度特性研究

介绍了节理面和岩桥各自的抗剪强度机制,引入法向变形协调条件,基于Mohr-Coulomb理论,推导了共面闭合断续节理岩体的直剪强度公式。模型试验发现,剪切破坏面以拉剪复合破坏为主,同时岩块中伴随大量的拉张微裂隙。试样的强度和变形具有明显的阶段性,全应力应变曲线主要经历了线弹性增长、节理面错动、次生裂纹起裂稳态扩展、节理面贯通破坏和残余强度5个阶段。对比发现,理论计算结果与试验值吻合较好。     

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

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

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.     

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

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

风化花岗岩声谱特征分析

岩石材料可视为声波的天然低通滤波器 ,声波在不同力学性质岩石材料中传播时 ,由于岩石对声波吸收程度不等因而高频声波滤除状况也不同。接收波中频率丰富程度与岩 石 的力学性质密切相关,岩石力学强度愈高,声谱中高频成分就愈丰富。反之,声谱中主要为 低频,高频分量因衰减而缺失。实验采用全数字化DB4型多波参数分析仪和计算机,对波形 数据进行FFT分析,通过接收谱和频率响应的特征来研究岩石的力学性质。声衰减研究可为 岩石、岩体强度和综合评价提供更多、更可靠的信息,是解决工程地质风化岩问题研究的一 种可行方法。     

Numerical Homogenization of Jointed Rock Masses Using Wave Propagation Simulation

Homogenization in fractured rock analyses is essentially based on the calculation of equivalent elastic parameters. In this paper, a new numerical homogenization method that was programmed by means of a MATLAB code, called HLA-Dissim, is presented. The developed approach simulates a discontinuity network of real rock masses based on the International Society of Rock Mechanics (ISRM) scanline field mapping methodology. Then, it evaluates a series of classic joint parameters to characterize density (RQD, specific length of discontinuities). A pulse wave, characterized by its amplitude, central frequency, and duration, is propagated from a source point to a receiver point of the simulated jointed rock mass using a complex recursive method for evaluating the transmission and reflection coefficient for each simulated discontinuity. The seismic parameters, such as delay, velocity, and attenuation, are then calculated. Finally, the equivalent medium model parameters of the rock mass are computed numerically while taking into account the natural discontinuity distribution. This methodology was applied to 17 bench fronts from six aggregate quarries located in Tunisia, Spain, Austria, and Sweden. It allowed characterizing the rock mass discontinuity network, the resulting seismic performance, and the equivalent medium stiffness. The relationship between the equivalent Young’s modulus and rock discontinuity parameters was also analyzed. For these different bench fronts, the proposed numerical approach was also compared to several empirical formulas, based on RQD and fracture density values, published in previous research studies, showing its usefulness and efficiency in estimating rapidly the Young’s modulus of equivalent medium for wave propagation analysis.     

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.     

Kelvin—Voigt黏弹性介质地震波场数值模拟与衰减特征

利用高阶交错网格有限差分模拟Kelvin-Voigt黏弹性介质中传播的地震波,同时将完全匹配层吸收边界条件引入到其边界处理中。数值模拟结果表明,完全匹配层吸收边界效果好,高阶有限差分能模拟得到的黏弹性介质波场精度较高。对模拟的黏弹性波场进行分析,表明介质的粘滞性使地震反射波的能量变弱,高频衰减明显,并比低频衰减得快,主频向低频方向移动,有效频带变窄,即降低了地震波的分辨率;并且反射转换波比反射纵波要衰减得快;而且还随着传播距离的增加,其峰值频率也逐渐降低。通过数值模拟分析具有不同的粘滞系数介质对地震波的吸收和衰减,结果表明随着粘滞系数的增大,地下介质对地震波的吸收衰减更明显。     

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

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

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.