动载确定方法对岩石动态断裂韧度测试的影响

为了考察不同方法确定的动态载荷对测试岩石动态断裂韧度的影响,在SHPB压杆系统上动态冲击直径80 mm的大理岩圆孔裂缝平台巴西圆盘,获得了弹性压杆上的应力波形,间接计算得到3种不同的作用在圆盘端部的动态载荷。将载荷输入ANSYS动态有限元模型中,求得了相应的动态应力强度因子,并根据试验-数值分析方法确定了岩石的动态断裂韧度测试值。结果表明,在加载速率约为4.0×104 MPa•m1/2/s的条件下,采用三波法确定的大理岩的平均动态断裂韧度为 3.92 MPa•m1/2,采用一波法比三波法计算的结果偏低11.22%,采用二波法比三波法计算的结果偏高20.15%,3种方法得到的结果差异较大。应力波在传播过程中,通过圆盘表面和预制裂缝面发生散射,部分能量不断发生释放是造成圆盘试件两端加载载荷不相等的主要原因。三波法是3种方法中比较理想的动态载荷确定方法,但需要考察试件的动态应力平衡性。     

动高压冲击下岩石破裂机理及相变研究

用高压气体炮作为原动力,在不同的冲击应力加载条件下,加速平板飞片撞击角闪岩样品。对回收样品应用SEM、XRD及RAM谱仪进行了宏观和微观研究,确定了样品出现的破裂起始阶段应力强度以及其裂纹的扩展形态,从而对岩石样的破坏过程进行了量化研究,并阐释了裂纹形成的机理和相变的存在。实验表明,动高压冲击下岩石破裂机理及相变的研究,可为矿山开采、建筑工程等爆破及应用提供相关的理论依据。     

爆炸气体驱动下岩体破裂的有限元-离散元模拟

在原有有限元/离散元（FEM/DEM）耦合分析方法中,实现了一种新的爆破计算模型。该模型考虑了在爆生气体的作用下,随着裂隙的扩展,气体占据的体积不断增大,气体压力逐渐减小这一问题。同时考虑了气体嵌入与爆腔联通的裂隙对裂隙的作用力。克服了原有FEM/DEM方法中的爆破模型仅仅将压力施加于爆腔四周的岩壁上,无法考虑爆生气体嵌入生成的裂隙对裂隙的作用。提出了一种新颖的贯通裂隙网络形成的递归搜索算法,只需通过编写一个简单的递归函数,即可实现复杂裂隙网络的搜索,采用一种简洁的方法完成了对复杂问题的处理。最后通过一个爆破算例,结果表明FEM/DEM方法可以对爆炸过程中应力波的传播及岩体中裂纹的萌生、扩展进行全程捕捉,展现了该方法用于爆破模拟的潜力。     

Thermal shock weakening of granite rock under dynamic loading: 3D numerical modeling based on embedded discontinuity finite elements

This paper presents a numerical study of thermal shock weakening of granite rock under dynamic loading. A fully 3D numerical scheme based on a combined continuum viscodamage-embedded discontinuity model and an explicit scheme to solve the underlying thermomechanical problem was developed and validated through numerical examples. First, the dynamic Brazilian disc test is simulated on intact numerical rock. Then, thermal shock-induced cracking due to a moving external heat flux boundary condition, mimicking experiments based on plasma jet treatment, is numerically predicted. Finally, numerical Brazilian disc test is conducted on the thermal shocked numerical samples. The predicted and experimental weakening effects are in good agreement demonstrating that the present modeling approach has good predictive capabilities. The practical significance of the results is that heat shock pretreatment can substantially enhance rock gravel and rubble crushing.     

Numerical modeling of dynamic rock fracture with a combined 3D continuum viscodamage‐embedded discontinuity model

This paper deals with numerical modeling of dynamic failure phenomena in rate‐sensitive quasi‐brittle materials, such as rocks, with initial microcrack populations. To this end, a continuum viscodamage‐embedded discontinuity model is developed and tested in full 3D setting. The model describes the pre‐peak nonlinear and rate‐sensitive hardening response of the material behavior, representing the fracture‐process zone creation, by a rate‐dependent continuum damage model. The post‐peak response, involving the macrocrack creation accompanied by exponential softening, is formulated by using an embedded displacement discontinuity model. The finite element implementation of this model relies upon the linear tetrahedral element, which seems appropriate for explicit dynamic analyses involving stress wave propagation. The problems of crack locking and spreading typical of embedded discontinuity models are addressed in this paper. A combination of two remedies, the inclusion of viscosity in the spirit of Wang's viscoplastic consistency approach and introduction of isotropic damaging into the embedded discontinuity model, is shown to be effective in the present explicit dynamics setting. The model performance is illustrated by several numerical simulations. In particular, the dynamic Brazilian disc test and the Kalthoff–Winkler experiment show that the present model provides realistic predictions with the correct failure modes and rate‐dependent tensile strengths of rock at different loading rates. The ability of initial embedded discontinuity populations to model the initial microcrack populations in rocks is also successfully tested. Copyright © 2016 John Wiley & Sons, Ltd.     

Fractal and geostatistical methods for modeling of a fracture network

The modeling of fracture networks is useful for fluid flow and rock mechanics studies. About 6600 fracture traces were recorded on drifts of a uranium mine in a granite massif. The traces have an extension of 0.20–20 m. The network was studied by fractal and by geostatistical methods but can be considered neither as a fractal with a constant dimension nor a set of purely randomly located fractures. Two kinds of generalization of conventional models can still provide more flexibility for the characterization of the network: (a) a nonscaling fractal model with variable similarity dimension (for a 2-D network of traces, the dimension varying from 2 for the 10-m scale to 1 for the centimeter scale, (b) a parent-daughter model with a regionalized density; the geostatistical study allows a 3-D model to be established where: fractures are assumed to be discs; fractures are grouped in clusters or swarms; and fracturation density is regionalized (with two ranges at about 30 and 300 m). The fractal model is easy to fit and to simulate along a line, but 2-D and 3-D simulations are more difficult. The geostatistical model is more complex, but easy to simulate, even in 3-D.This paper was presented at Emerging Concepts, MGUS-87 Conference, Redwood City, California, 13–15 April 1987.     

基于P-CCNBD试样的岩石动态断裂韧度测试方法

预裂的人字形切槽巴西圆盘(Pre-cracked chevron notched Brazilian disc,简称P-CCNBD)是将人字形切槽巴西圆盘(cracked chevron notched Brazilian disc,简称CCNBD)的切槽尖端再稍加切削制成直裂纹前沿的试样。利用霍普金森压杆对P-CCNBD砂岩试样进行径向冲击,完成I型动态断裂试验后再做数值分析得到岩石的动态断裂韧度。为了验证数值模拟的可靠性,先进行了无限平面中一条有限尺寸裂纹表面受冲击拉伸作用的动态有限元分析,结果表明,数值模拟的结果与Shi得到的结果非常吻合。将试验-数值法和他人的准静态法分别确定的砂岩的动态起裂韧度进行对比,两种方法得到的结果有一定的差异。采用试验-数值法,将比较成熟的直裂纹巴西圆盘(cracked straight-through Brazilian disc,简称CSTBD)和P-CCNBD两种试样测得的结果进行对比,两者吻合较好。得到的动态起裂韧度都有随着加载率的增加而增大的加载率效应。分析了准静态法的缺陷,认为试验-数值法得到的结果更为合理。     

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

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

基于离散裂隙网络模型的节理岩体渗透张量及特性分析

节理岩体几何结构非常复杂,研究其渗流特性对于指导含水岩层稳定性分析具有重要价值。应用离散裂隙网络模型DFN方法,基于VC++6.0软件平台,建立了平面渗流分析方法,分析了节理岩体不同几何分布情况下的渗透率张量特征,通过定义渗流定向性系数对岩体渗流的定向性特征进行了定量分析。结果表明：单组节理岩体渗流具有明显的各向异性特征,渗流定向性随着节理角度变化显著;节理随着节理贯通性增加,节理渗透率呈现对数增加趋势;两组节理情况下,各向异性特征随着节理组间夹角变化;两组节理岩体渗流特征研究中,正交分布下,岩体仍存在各向异性,但渗流定向性系数较低;当节理倾角服从正态分布时,随着节理倾角标准差增大,渗透率增加;两组节理夹角不同时,节理渗透主方向倾角随着夹角增大而相应增大,基本沿两组节理夹角方向的角平分线方向。     

动静组合作用下刀具破岩机制数值分析

基于细观损伤有限元方法,模拟分析了刀具在单一动载、动静联合荷载、静态围压条件下动静联合荷载3种情况下岩体破碎的全过程。模型采用黏弹性人工边界剔除了边界应力波反射对模拟结果准确性的影响。数值模拟结果表明：在弹性情况下,静压的存在对岩体内部最小主应力值影响不大,却显著提高了材料内部最大主应力水平,增大了剪应力的大小,导致剪切破坏可能性增加;当有围压存在时,岩体内部受拉区域减少,岩体强度有所提高。单一动载和动静联合荷载破岩时,岩体内部除刀头附近呈现少量压破坏外,破坏均以拉破坏为主;而围压条件下,岩体破碎面积相对减小,裂纹在围压的作用下向两侧自由面延伸,岩体内部破坏形式则趋于多样化,压破坏比重明显增大,整体表现为拉压复合作用。模拟结果还表明,刀头侵入量主要受动载力大小影响,在相同幅值增量下,动载力增加导致的刀头侵入量远大于静压增加所导致的侵入量。相对单一动载和静压作用下的岩石破碎机制来说,动静组合加载破岩的研究还需更为深入的探讨。研究结果可望对岩体破坏机制以及地下工程作业等实际应用提供一定的参考。     

杂交边界点方法求解动态断裂力学问题

结合杂交边界点法和双互易法则,推导出求解动力问题的纯边界类型无网格方法--双互易杂交边界点方法,并将该方法用于求解含中心裂纹的方板受瞬态载荷作用的问题。该方法将问题的解分为通解和特解两部分,通解使用杂交边界点法求解,特解则利用局部径向基函数近似,域内布点仅仅为了径向基插值,因此仍然是一种纯边界类型的无网格方法。同时,将移动最小二乘近似中的基函数扩充,使该方法能更好地模拟裂纹尖端应力场的奇异性,具有后处理简单、精度高的优点。数值算例表明了该方法的稳定性和有效性。     

含顺层断续节理岩质边坡地震作用下的破坏模式与动力响应研究

基于二维颗粒流软件PFC2D的人工合成岩体技术（SRM）,研究了岩桥倾角和节理间距不同组合形式的含顺层断续节理岩质边坡在地震作用下的破坏模式与动力响应规律。研究结果显示：在地震动力作用下,含单潜在滑动面的顺层断续节理岩质边坡呈现出滑移-倾倒的混合破坏特征,含多潜在滑动面的顺层断续节理岩质边坡则主要发生倾倒破坏;由顺层断续节理以及岩桥交替连接所组成的潜在滑动面是控制边坡动力稳定性的关键因素。在地震动力作用下,最靠近坡脚的岩桥段首先萌生翼裂纹,使得拉应力得到释放,随后各节理相继萌生裂纹并扩展、贯通,最终导致坡体发生阶梯状整体失稳。裂纹扩展受顺层断续节理控制,萌生裂纹中以张拉裂纹为主,且裂纹数量与输入地震波的加速度曲线具有同步性。另一方面,节理面的存在对边坡动力响应产生明显影响,沿坡表以及沿水平方向上的峰值速度、峰值位移随着岩桥倾角的增大、节理间距的减小而增大,同时节理间距和岩桥倾角对于峰值加速度（PGA）放大系数的影响范围主要集中在坡表、坡肩;沿竖直方向上,峰值位移随着岩桥倾角、节理间距的增大而减小,PGA放大系数曲线随高程变化总体呈现U型分布特征。     

不连续岩体-拱坝系统的动力过程数值模拟方法

提出了一种分析三维不连续变形块体系统动力学过程和建坝的数值方法,该方法不但能够模拟位移不连续面在地震载荷作用下的张开和大尺度的摩擦滑移,而且还能给出不连续面上安全系数的分布和随时间的变化.通过有2个块体组成的系统验证了所提出的方法的正确性,同时用它分析了一个具有断层和伸缩缝的拱坝系统模型在地震载荷下的动力反应.     

单轴压缩非贯通节理岩石尖端塑性区及弹塑性断裂模型研究

基于Drucker-Prager（下简称D-P）准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。     

岩石非理想裂纹圆盘试件动态断裂韧性测试的有限元分析及试验研究

采用中心裂纹圆盘-霍普金森压杆(CCCD-SHPB)试验系统对大理岩中心裂纹圆盘试件进行不同加载速率下的纯Ⅰ型加载试验,研究岩石材料动态断裂韧性的加载速率相关性。考虑到试验中所使用的圆盘试件是带有中心切口的非理想裂纹,结合试件的实际加工情况,提出采用切口尖端形状为小圆弧形的预制裂纹,介绍了试件的制作方法,并利用有限元计算论证了采用这种形状非理想裂纹试件的试验可行性。结果表明,采用圆弧形裂尖、裂纹宽度为1 mm的中心切口非理想裂纹圆盘试件来代替理想裂纹圆盘试件是可行的,误差不超过2.13%。分析试验数据得出：非理想裂纹圆盘试件在纯Ⅰ型加载条件下,其断裂韧度表现出明显的速率相关性,且随着加载速率的增大而增大。     

Development and validation of computer models for structures in jointed rock

An extensive program was conducted to develop and validate computer models for the deformation of structures in jointed rock. Both discrete element and finite element models were investigated. The program comprised a sequence of seven calculational and experimental studies. In the final step, model predictions were compared with experimental data from a well-instrumented jointed-rock precision test. The results show that the sophisticated models developed were able to predict some elements of structural deformation in jointed rock, but that it is difficult to make high-confidence predictions of all aspects of structural deformation with these models. Copyright © 1999 John Wiley & Sons, Ltd.     

Improvements for the smooth joint contact model of the particle flow code and its applications

This paper deals with two shortcomings of the smooth-joint contact model (SJCM) used in the particle flow code (PFC). The first shortcoming is the increase of the shear strength of the joint when the shear displacement of the joint exceeds a specific value that is related to the particle size. This problem is named as the interlocking problem, which is caused by the interlocking particles. It occurs due to a shortcoming of the updating procedure in the PFC software related to the contact conditions of the particles that lie around the intended joint plane during high shear displacements. This problem also increases the dilation angle and creates unwanted fractures around the intended joint plane. To solve this problem two new approaches are proposed in this paper: (1) joint plane checking (JPC) approach and (2) joint sides checking (JSC) approach. These approaches and the regular approach are used to model: (a) the direct shear test using the PFC^2D and PFC^3D, (b) the biaxial test on a sample having a persistent joint with a dip angle varying from 0° to 90° at an interval of 15° using the PFC^2D and (c) the polyaxial test on two samples, one of them having a joint with a dip direction of 0° and the dip angle varying from 0° to 90° at an interval of 15°, and the other sample having a joint with a dip angle of 60° and the dip direction varying from 0° to 90° at an interval of 15° using the PFC^3D. All numerical results show that the JPC and JSC approaches can solve the interlocking problem. Also, they proved to be more consistent with the theory compared to the regular approach. However, the JPC approach leads to a slightly softer joint. Therefore, the JSC approach is suggested for jointed rock modeling using the PFC. The other shortcoming of the SJCM dealt within this paper is its inability to capture the non-linear behavior of the joint closure varying with the joint normal stress. This problem is solved in this paper by proposing a new modified smooth-joint contact model (MSJCM). MSJCM uses a linear relation between the joint normal stiffness and the normal contact stress to model the non-linear relation between the joint normal deformation and the joint normal stress observed in the compression joint normal stiffness test. A good agreement obtained between the results from the experimental test and the numerical modeling of the compression joint normal test shows the accuracy of this new model.     

隧道岩爆的动力学机理及其控制的实验研究

隧道爆破开挖过程对隧道围岩的扰动,可以用弹性动力学的Lamb问题来描述。利用动态光弹性方法,对Lamb模型进行了验证,研究了冲击荷载下模型内应力波的传播、相互作用以及刻槽的拦截效应。结果表明：（1）冲击荷载作用下,模型内传播的有P波、S波和Rayleigh波,与Lamb模型的力学解答是一致的;（2）等差条纹图的分布规律,可以合理解释岩爆高发区段与开挖面的位置关系,认为P波和Rayleigh波的致裂作用和传播规律是岩爆发生的重要因素;（3）在模型表面的刻槽,可以实现对岩爆有突出贡献的Rayleigh波的大量拦截,说明在隧道表面切槽的岩爆控制方法是可行的。     

Damage‐viscoplastic consistency model for rock fracture in heterogeneous rocks under dynamic loading

This paper presents a damage‐viscoplastic consistency model for numerical simulation of brittle fracture in heterogeneous rocks. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept with separate damage variables in tension and compression. This approach does not suffer from ill‐posedness, caused by strain softening, of the underlying boundary/initial value problem since viscoplasticity provides the regularization by introducing a length scale effect under dynamic loading conditions. The model uses the Mohr–Coulomb yield criterion with the Rankine criterion as a tensile cut‐off. The damage law in compression is calibrated via the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition occurring in confined compression at a certain level of confinement. The heterogeneity of rock is accounted for by the statistical approach based on the Weibull distribution. Numerical simulations of confined compression test in plane strain conditions demonstrate a good agreement with the experiments at both the material point and structural levels as the fracture modes are realistically predicted. Copyright © 2009 John Wiley & Sons, Ltd.