Processing math: 100%
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

初始应力下岩体爆破损伤特性及破裂机理

马泗洲 刘科伟 杨家彩 李旭东 郭腾飞

胡宏伟, 宋浦, 王建灵, 郭炜, 徐洪涛, 金鹏刚, 任松涛. 炸药水中爆炸冲击因子的新型计算方法[J]. 爆炸与冲击, 2014, 34(1): 11-16.
引用本文: 马泗洲, 刘科伟, 杨家彩, 李旭东, 郭腾飞. 初始应力下岩体爆破损伤特性及破裂机理[J]. 爆炸与冲击, 2023, 43(10): 105201. doi: 10.11883/bzycj-2023-0151
Hu Hong-wei, Song Pu, Wang Jian-ling, Guo Wei, Xu Hong-tao, Jin Peng-gang, Ren Song-tao. A new calculation method for shock factor of underwater explosion[J]. Explosion And Shock Waves, 2014, 34(1): 11-16.
Citation: MA Sizhou, LIU Kewei, YANG Jiacai, LI Xudong, GUO Tengfei. Blast-induced damage characteristics and fracture mechanism of rock mass under initial stress[J]. Explosion And Shock Waves, 2023, 43(10): 105201. doi: 10.11883/bzycj-2023-0151

初始应力下岩体爆破损伤特性及破裂机理

doi: 10.11883/bzycj-2023-0151
基金项目: 国家自然科学基金(51974360)
详细信息
    作者简介:

    马泗洲(1995- ),男,硕士研究生,sizhou_ma@126.com

    通讯作者:

    刘科伟(1982- ),男,博士,教授, kewei_liu@csu.edu.cn

  • 中图分类号: O383

Blast-induced damage characteristics and fracture mechanism of rock mass under initial stress

  • 摘要: 钻爆法在岩体开挖工程中应用广泛,其工作效率与岩体结构及地质环境密切相关。初始应力会显著影响岩体爆破裂纹的扩展行为及破坏特征,导致深部岩体超/欠挖等问题。本文采用理论分析和数值模拟相结合的方法研究了不同初始应力下岩体的爆破损伤特性及破裂机理。基于弹性力学建立单孔爆破动静组合理论模型,分析了静态应力分布及动态应力演化特征,揭示了初始应力作用下岩体爆破损伤机制。通过经验公式及动态力学试验对Riedel-Hiermaier-Thoma (RHT) 模型参数进行标定与修正,并结合室内爆破实验及理论结果对数值模型进行验证。此外,在单孔爆破数值模型中分析了爆炸压力演化规律、裂纹扩展行为及分形特征,结果表明数值模拟结果可以较好地论证理论模型的合理性。研究发现:环向拉应力是影响爆破裂纹扩展的重要因素,当初始应力较大时,合理调整环向应力分布可改善岩体爆破碎裂效果,实际工程中可结合预裂技术进行地压调控改变应力分布。
  • 水中兵器对舰船的破坏程度是由很多因素决定的,主要包括主装药的性质和质量、爆心与舰船的距离和位置、舰船的结构和性质等。在评价水中兵器的破坏威力时,引入了冲击因子Q这样一个衡量标准[1],冲击因子的物理意义是,对于同一舰船,若冲击因子相等,则认为其水下爆炸的冲击响应近似相等。以前,使用一种基于冲击波超压的冲击因子考核潜艇结构生命力,这也是冲击因子的最早形式。由大量的实验和数值计算可发现,该形式的冲击因子不能很好地反映结构的破坏程度。随着对水下爆炸现象的进一步研究,越来越倾向于使用另一种基于平面波假定的冲击因子,这种冲击因子是从作用在结构上的冲击波能的角度定义的,远场时具有较好的效果,但该冲击因子没有考虑近场冲击波形状和能量损耗。姚熊亮等[2]利用球面波理论对水中爆炸冲击因子进行了修正,考虑了冲击波形状对垂直投射到舰船结构的冲击波能的影响。但是,水中爆炸近场的能量损耗非常严重,如48%冲击波能损耗在25个装药半径的范围内[3],而且不同种类炸药的能量释放特性也存在一定的差异。然而,传统冲击因子中描述炸药性能的TNT当量只是总化学能的一个比,并不能完全表征炸药水中爆炸的能量输出与衰减特性,尤其是水中爆炸近场。本文中,通过冲击波峰值压力和冲量相似方程的乘积推导一种适用于水中爆炸冲击波因子的计算方法,从冲击波毁伤作用角度描述该计算方法的物理意义,通过水中爆炸实验进行验证,并与传统的水下爆炸冲击因子进行对比。

    水中爆炸冲击波的毁伤作用可用冲击因子表示,基于平面波的冲击因子的表达式为[4]:

    Q1=1+sinα2WnR

    式中:Q1为冲击因子;W 为炸药的质量(TNT当量),kg;R为目标距爆心的距离,m;α为冲击波的入射角;n为通过实验确定的装药指数,常取n=0.5。式(1)仅适用于炸药在水底以上爆炸,如果炸药被淤泥覆盖,则相当多的炸药能量消耗在淤泥中,冲击因子不再适用。

    对于潜艇目标,冲击波直接作用的结构表面即为垂直于冲击波传播方向的潜艇表面。此时,可认为α=90°,(1+sin90°)/2=1。当水面舰船遭受水下爆炸冲击时,冲击波是向各个方向传播的。其中,只有一部分冲击波能量对水面舰船起作用,这部分能量与水面舰船在垂直于冲击波方向的投影面积成正比,通常用(1+sinα)/2表示到达水面船只的冲击波方向的影响。在特定的攻击位置和方位,可把水底和水面的影响因素可看作一个常数,因此式(1)可简化为[5]:

    Q1=W0.5/R

    这种冲击因子实际上是基于平面波假设,并从结构遮挡冲击波能量的角度定义的,与冲击波能的关系可表述为[2]:

    Q21=4πEsρeηeSe

    式中:E为冲击波总能量,MJ/kg;Es为结构遮挡的冲击波能,MJ/kg;Se为结构在垂直于冲击波阵面上的投影面积,m2e为炸药的质量化学能,MJ/kg;ηe为炸药化学能转化为冲击波能的比例。

    由式(3)可知,当水下爆炸冲击波为平面波时,无论目标距爆心的距离大小如何,Se为常数。在爆炸远场,匀化冲击波可近似为平面波,而且冲击波能基本恒定,Q为Es的函数。在爆炸近场,冲击波能在不断衰减,冲击波的形状不能近似为平面波,此时,Q为ηe、Es的函数,在某一距离处,ηe与炸药的性能有关,Es与冲击波的形状和目标结构有关。因此,炸药的TNT当量并不能完全表征近场冲击波的特性,基于平面波的水中爆炸冲击因子具有一定局限性。而通过水中爆炸实验得到的冲击波相似方程包含了ηe的特性,反映了炸药水下爆炸的能量释放和衰减特性,下面以冲击波相似方程为基础对冲击因子进行推导。

    水中爆炸的相似方程为[3]:

    Y=K(W1/3/R)α

    式中:W 为炸药的质量,kg;R为测点距爆心的距离,m;K和α为与炸药相关的系数,K为相似常数,α为相似指数;冲击波参数Y 包括峰值压力pmax、比时间常数θ/W1/3、比冲量I/W1/3和比能流密度E/W1/3

    由水中爆炸实验得到的峰值压力的经验计算公式为:

    pmax=Kp(W1/3/R)αp

    冲量的计算公式为:

    I+=KIW1/3(W1/3/R)αI

    将峰值压力pmax和冲量I+相乘,得到:

    pmaxI+=KpKIWαp/3W1/3WαI/3RαpRαI

    整理得:

    (pmaxI+KpKI)1/αp+αI)=W(αp+αI+1)/3(αp+αI)R

    可简化为[6]:

    A(pmaxI+)m=WnR

    式中:m=1/(αpI),A=(KpKI)m,n=其中m和A为冲击波参数的影响因子,n为装药指数。

    通过水中爆炸实验得到的冲击波相似系数[7],计算了几种典型炸药的装药指数n,见表1

    表  1  几种典型炸药的相似系数和装药指数
    Table  1.  Exponentsα,K,nfor various high explosives
    下载: 导出CSV 
    | 显示表格

    表1可看出,在实验精度许可范围内,所有以TNT为基本组分的炸药,式(9)中装药指数n均可近似等于0.5。因此,冲击因子可表示为:

    Q2=A(pmaxI+)m=W0.5R

    式(10)将峰值压力pmax冲量I的毁伤作用双曲线与冲击因子的R-W关系曲线联系起来,从冲击波的毁伤作用角度表述了水中爆炸冲击因子的物理意义。

    TNT、RDX基和HMX基含铝炸药,全部为圆柱形压装炸药,配方组成和爆热见表2。药柱的长径比为(1.0~1.2)∶1,一端带雷管孔。Al粉为球形,直径为4.5~5.5μm。黏结剂中,w(F2603)=1.5%,w(蜡)=3%,w(石墨)=0.5%。试样质量均为25g,采用8号铜电雷管端面起爆。

    表  2  炸药配方及爆热
    Table  2.  Explosive formulation and explosion heat
    下载: 导出CSV 
    | 显示表格

    水池∅3.2m×2.6m,水深2.4m,池底和池壁均由8mm钢板焊接而成。入水深度1.6m,为总水深的2/3,满足冲击波和气泡的测试要求,可以消除边界效应的影响[8]。爆心与传感器的距离分别为23、32、44、58和90cm,测量的是炸药柱中心轴向的径向冲击波,实验布局如图1所示。

    测试系统包括138系列ICP型压电式电气石水下激波传感器和482A型信号适配器、高低频数据记录仪。冲击波信号的采样频率为10MHz。

    图  1  实验布局图
    Figure  1.  Experimental layout

    通过实验测量的冲击波参数,拟合得到了不同组分炸药的相似系数和装药指数,见表3

    表3可知,由pmaxI+推导的冲击因子计算的理想炸药、RDX基和HMX基含铝炸药的装药指数n都非常接近于0.5,平行性也非常好。可见,基于平面波假设的冲击因子的装药指数同样适合于RDX、HMX基的含铝炸药。

    表  3  不同组分炸药的相似系数和装药指数
    Table  3.  Exponentsα,K,nfor various explosives compositions
    下载: 导出CSV 
    | 显示表格

    为了对比基于平面波的水中爆炸冲击因子与由峰值压力与冲量的乘积pmaxI+推导的冲击因子计算公式的差异,依据实验数据计算了测距23、58和90cm的冲击因子,3个距离的水中爆炸冲击波参数见表4,不同距离处的冲击因子见表5图2

    表  4  水中爆炸冲击波峰值压力和冲量
    Table  4.  Shock wave peak pressure and impluse of underwater explosion
    下载: 导出CSV 
    | 显示表格

    表5图2可知,在爆炸远场,Q1与Q2一致性很好,但随着测距的减小,Q1与Q2的误差逐渐增大。例如在90cm处,Q1与Q2非常接近,但在23cm处,Q1基本都大于Q2,并且随着铝含量的增加误差越来越大。这是由于,冲击因子Q1仅包含炸药质量、爆距等参数,没有包含冲击波形状、炸药的能量输出结构和冲击波衰减特性对冲击因子的影响。在近场,冲击波形状不能完全使用球面波理论近似[9],应根据炸药的形状和爆距确定。另外,铝粉含量也会影响冲击波能、气泡能的分配比例和冲击波能量的衰减,例如,对于铝含量为0%~30%的RDX基含铝炸药,冲击波能占总化学能的比例在40%~60%,药柱18倍半径处,初始冲击波能损失了约50%~60%[10]。因此,炸药的TNT当量难以全面反映水下爆炸冲击波的特性,结合冲击形状的影响,水中爆炸冲击因子可用下式计算:

    表  5  不同距离处的冲击因子
    Table  5.  The shock factors at different distances
    下载: 导出CSV 
    | 显示表格

    Q2 =AB(pmaxI+)m (11)式中:B为装药形状影响因子[11]。在近场,球形装药和L/D=1的圆柱形装药,B=1.00;L/D=6的圆柱形装药,B=1.08~1.10;锥形装药,B=1.02~1.03。在远场时,B≈1.00,式(11)可简化成式(10)。

    图  2  不同距离处的冲击因子
    Figure  2.  The shock factors at different distances

    (1)水中爆炸冲击因子装药指数n=0.5,不仅适合所有以TNT为基本组分的炸药,也适合于RDX、HMX基的含铝炸药。

    (2)由峰值压力与冲量的乘积pmaxI+推导的冲击因子计算公式,从冲击波的毁伤作用的角度表述了水中爆炸冲击因子的物理意义,结合冲击形状对冲击波能的影响,计算近场冲击因子时具有更高的准确性。

  • 图  1  理论分析几何模型

    Figure  1.  Geometrical model for theoretical analysis

    图  2  炮孔周边静态应力分布

    Figure  2.  Stress distribution around the borehole under static load

    图  3  不同位置处岩体静态应力的变化

    Figure  3.  Stress variation of rock mass at different distances under static load

    图  4  不同距离处岩体的动态应力波演化

    Figure  4.  Dynamic stress waves evolution in rock mass at various distances

    图  5  RHT模型破坏面的相关参数拟合

    Figure  5.  Fitting curve of failure surface parameters for RHT model

    图  6  数值模拟与试验中岩体冲击破坏过程比较

    Figure  6.  Comparison of failure process between numerical simulation and test

    图  7  数值模拟与试验中岩体应力时程曲线的对比

    Figure  7.  Comparison of stress with time between numerical simulation and test

    图  8  数值模型及局部网格

    Figure  8.  Configuration of numerical model and local mesh

    图  9  岩体爆破裂纹形态理论、试验与模拟结果比较

    Figure  9.  Comparison of rock blasting crack patterns among theoretical, experimental and simulated results

    图  10  压力衰减曲线比较

    Figure  10.  Comparison of attenuation curve of peak pressure

    图  11  数值模型及压力测点布置

    Figure  11.  Numerical model and arrangement of the measurement points

    图  12  初始应力下岩体爆破计算流程

    Figure  12.  Flow chart for the rock blast under initial stress

    图  13  应力初始化验证

    Figure  13.  Verification of stress initialization

    图  14  应力初始化模拟结果 (K=0)

    Figure  14.  Numerical results of the stress initialization (K=0)

    图  15  应力初始化模拟结果 (K=0.5)

    Figure  15.  Numerical results of the stress initialization (K=0.5)

    图  16  应力初始化模拟结果 (K=1.0)

    Figure  16.  Numerical results of the stress initialization (K=1.0)

    图  17  应力初始化模拟结果 (K=1.5)

    Figure  17.  Numerical results of the stress initialization (K=1.5)

    图  18  应力初始化模拟结果 (K=2.0)

    Figure  18.  Numerical results of the stress initialization (K=2.0)

    图  19  无初始压力下环向应力的演化过程 (E-1)

    Figure  19.  Evolution of hoop stress without initial pressure (E-1)

    图  20  各向同性压力下环向应力的演化过程 (E-4)

    Figure  20.  Evolution of hoop stress under equibiaxial pressure (E-4)

    图  21  各向异性压力下环向应力的演化过程 (A-4)

    Figure  21.  Evolution of hoop stress under anisotropic pressure (A-4)

    图  22  不同压力条件下环向应力时程曲线

    Figure  22.  Time histories of hoop stress under various pressure conditions

    图  23  不同距离峰值环向应力变化

    Figure  23.  Peak hoop stress variation at various distances

    图  24  各向同性压力下爆破裂纹扩展特征

    Figure  24.  Characteristics of blasting crack propagation under different equibiaxial pressures

    图  25  各向异性压力下爆破裂纹扩展特征

    Figure  25.  Characteristics of blasting crack propagation under different anisotropic pressures

    图  26  爆破裂纹盒子划分

    Figure  26.  Box division for blasting cracks

    图  27  各向同性压力下爆破裂纹分形维数拟合线

    Figure  27.  Fractal dimension of blasting cracks and its fitting lines under equibiaxial pressure

    图  28  各向异性压力下爆破裂纹分形维数拟合线

    Figure  28.  Fractal dimension of blasting cracks and its fitting lines under anisotropic pressure

    表  1  不同围压下花岗岩力学参数

    Table  1.   Mechanical parameters of granite sample under various confining pressure

    σ2/MPa σ3/MPa σ1/MPa p0 σf σ2/MPa σ3/MPa σ1/MPa p0 σf
    0 0 162 0.33 1.00 60 60 569 1.42 3.14
    10 10 265 0.59 1.58 70 70 616 1.56 3.37
    20 20 342 0.79 1.99 80 80 660 1.69 3.58
    30 30 408 0.96 2.33 90 90 703 1.82 3.79
    40 40 466 1.12 2.63 100 100 466 1.94 3.98
    50 50 520 1.28 2.90
    下载: 导出CSV

    表  2  岩石RHT模型材料参数

    Table  2.   RHT model parameters for rock mass

    参数名称 符号 来源 参数名称 符号 来源
    密度 ρ0 2620 kg/m3 试验测定 损伤因子 D1 0.04 参考文献[14]
    初始孔隙度 α 1.00 状态方程参数 B0 1.22
    抗压强度 fc 162 MPa 状态方程参数 B1 1.22
    参考压缩应变率 ˙εc0 3.0×10−5 s−1 模型给定 拉伸体积塑性应变分数 Pft 0.001
    参考拉伸应变率 ˙εt0 3.0×10−6 s−1 孔隙压实压力 Pcomp 6.00 GPa 理论计算
    破坏压缩应变率 ˙εc 3.0×1025 s−1 压缩应变率指数 βc 0.008
    破坏拉伸应变率 ˙εt 3.0×1025 s−1 拉伸应变率指数 βt 0.011
    损伤因子 D2 1.00 状态方程参数 T1 33.95 GPa
    侵蚀塑性应变 εfp 2.00 默认取值 状态方程参数 T2 0.00 GPa
    孔隙度指数 NP 3.0 Hugoniot多项式系数 A1 33.95 GPa
    最小损伤残余应变 εmp 0.012 Hugoniot多项式系数 A2 41.42 GPa
    剪切模量减小因子 ξ 0.50 Hugoniot多项式系数 A3 8.71 GPa
    相对抗剪强度 Fs 0.18 试验优化 孔隙坍塌压力 pcrush 108
    相对抗拉强度 Ft 0.06 洛德角相关因子 Q0 0.68
    压缩屈服面参数 Gc 0.50 洛德角相关因子 B 0.05
    拉伸屈服面参数 Gt 0.70 破坏面参数 A 2.48
    残余面参数 Af 1.62 破坏面参数 N 0.79
    残余面参数 Nf 0.62 弹性剪切模量 G 21.9 GPa
    下载: 导出CSV

    表  3  炸药模型材料参数

    Table  3.   Parameters for the explosive material

    ρe/(kg·m−3) vd /(m·s−1) pCJ/GPa Ee0/(J·m−3) Ae/GPa Be/GPa R1 R2 ω
    1320 6690 16 7.38×109 586 21.6 5.81 1.77 0.282
    下载: 导出CSV

    表  4  空气模型材料参数

    Table  4.   Mateiral parameters for the air

    ρa/(kg·m−3) C0 C1 C2 C3 C4 C5 C6 Ea0/(kJ·m−3) V0
    1.29 0 0 0 0 0.4 0.4 0 250 1.0
    下载: 导出CSV

    表  5  初始应力加载条件

    Table  5.   Initial stress conditions in numerical simulation

    应力状态 工况 埋深/m σx/MPa σy/MPa 应力状态 工况 埋深/m σx/MPa σy/MPa
    各向同性初始压力 E-1 0 0 0 各向异性初始压力 A-1 750 5 20
    E-2 375 10 10 A-2 750 10 20
    E-3 750 20 20 A-3 1125 30 20
    E-4 1125 30 30 A-4 1500 40 20
    下载: 导出CSV
  • [1] LI X B, GONG F Q, TAO M, et al. Failure mechanism and coupled static-dynamic loading theory in deep hard rock mining: a review [J]. Journal of Rock Mechanics and Geotechnical Engineering, 2017, 9(4): 767–782. DOI: 10.1016/j.jrmge.2017.04.004.
    [2] 杨建华, 孙文彬, 姚池, 等. 高地应力岩体多孔爆破破岩机制 [J]. 爆炸与冲击, 2020, 40(7): 075202. DOI: 10.11883/bzycj-2019-0427.

    YANG J H, SUN W B, YAO C, et al. Mechanism of rock fragmentation by multi-hole blasting in highly-stressed rock masses [J]. Explosion and Shock Waves, 2020, 40(7): 075202. DOI:10.11883/bayzj-2019-0427. DOI: 10.11883/bzycj-2019-0427.
    [3] 朱万成, 唐春安, 左宇军. 深部岩体动态损伤与破裂过程 [M]. 北京: 科学出版社, 2014: 112–121.
    [4] NICHOLLS H R, DUVALL W I. Presplitting rock in the presence of a static stress field [R]. Washington DC: U. S. Department of the Interior, Bureau of Mines, 1966.
    [5] KUTTER H K, FAIRHURST C. On the fracture process in blasting [J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1971, 8(3): 181–202. DOI: 10.1016/0148-9062(71)90018-0.
    [6] 肖正学, 张志呈, 李端明. 初始应力场对爆破效果的影响 [J]. 煤炭学报, 1996, 21(5): 497–501. DOI: 10.13225/j.cnki.jccs.1996.05.011.

    XIAO Z X, ZHANG Z C, LI D M. The influence of initial stress field on blasting [J]. Journal of China Coal Society, 1996, 21(5): 497–501. DOI: 10.13225/j.cnki.jccs.1996.05.011.
    [7] 刘殿书, 王万富, 杨吕俊, 等. 初始应力条件下爆破机理的动光弹实验研究 [J]. 煤炭学报, 1999, 24(6): 612–614. DOI: 10.13225/j.cnki.jccs.1999.06.012.

    LIU D S, WANG W F, YANG L J, et al. Holophotoelasticity study on mechanism of blasting under initiative stress field [J]. Journal of China Coal Society, 1999, 24(6): 612–614. DOI: 10.13225/j.cnki.jccs.1999.06.012.
    [8] YANG R S, DING C X, LI Y L, et al. Crack propagation behavior in slit charge blasting under high static stress conditions [J]. International Journal of Rock Mechanics and Mining Sciences, 2019, 119: 117–123. DOI: 10.1016/J.ijrmms.2019.05.002.
    [9] 杨仁树, 肖成龙, 丁晨曦, 等. 空孔与运动裂纹相互作用的动焦散线实验研究 [J]. 爆炸与冲击, 2020, 40(5): 052202. DOI: 10.11883/bzycj-2019-0091.

    YANG R S, XIAO C L, DING C X, et al. Experimental study on dynamic caustics of interaction between void and running crack [J]. Explosion and Shock Waves, 2020, 40(5): 052202. DOI: 10.11883/bzycj-2019-0091.
    [10] 岳中文, 田世颖, 张士春, 等. 单向围压作用下切缝药包爆破爆生裂纹扩展规律的研究 [J]. 振动与冲击, 2019, 38(23): 186–195. DOI: 10.13465/j.cnki.jvs.2019.23.027.

    YUE Z W, TIAN S Y, ZHANG S C, et al. Expanding law of cracks formed by slotted cartridge blast under unidirectional confining pressure [J]. Journal of Vibration and Shock, 2019, 38(23): 186–195. DOI: 10.13465/j.cnki.jvs.2019.23.027.
    [11] YANG L Y, HUANG C, BAO S J, et al. Model experimental study on controlled blasting of slit charge in deep rock mass [J]. Soil Dynamics and Earthquake Engineering, 2020, 138: 106318. DOI: 10.1016/j.soildyn.2020.106318.
    [12] MA G W, AN X M. Numerical simulation of blasting-induced rock fractures [J]. International Journal of Rock Mechanics and Mining Sciences, 2008, 45(6): 966–975. DOI: 10.1016/j.ijrmms.2007.12.002.
    [13] XIE L X, LU W B, ZHANG Q B, et al. Damage evolution mechanisms of rock in deep tunnels induced by cut blasting [J]. Tunnelling and Underground Space Technology, 2016, 58: 257–270. DOI: 10.1016/j.tust.2016.06.004.
    [14] XIE L X, LU W B, ZHANG Q B, et al. Analysis of damage mechanisms and optimization of cut blasting design under high in-situ stresses [J]. Tunnelling and Underground Space Technology, 2017, 66: 19–33. DOI: 10.1016/j.tust.2017.03.009.
    [15] LI X D, LIU K W, YANG J C, et al. Numerical study on blast-induced fragmentation in deep rock mass [J]. International Journal of Impact Engineering, 2022, 170: 104367. DOI: 10.1016/j.ijimpeng.2022.104367.
    [16] LI X D, LIU K W, SHA Y Y, et al. Numerical investigation on rock fragmentation under decoupled charge blasting [J]. Computers and Geotechnics, 2023, 157: 105312. DOI: 10.1016/j.compgeo.2023.105312.
    [17] KIRSCH G. Die theorie der elastizitat und die bedurfnisse der festigkeitslehre [J]. Zantralblatt Verlin Deutscher Ingenieure, 1898, 42: 797–807.
    [18] MIKLOWITZ J, KAUL R K. The theory of elastic waves and waveguides [J]. Journal of Applied Mechanics, 1979, 46(4): 969. DOI: 10.1115/1.3424709.
    [19] YI C P, JOHANSSON D, GREBERG J. Effects of in-situ stresses on the fracturing of rock by blasting [J]. Computers and Geotechnics, 2018, 104: 321–330. DOI: 10.1016/j.compgeo.2017.12.004.
    [20] TAO J, YANG X G, LI H T, et al. Effects of in-situ stresses on dynamic rock responses under blast loading [J]. Mechanics of Materials, 2020, 145: 103374. DOI: 10.1016/j.mechmat.2020.103374.
    [21] LI X H, ZHU Z M, WANG M, et al. Numerical study on the behavior of blasting in deep rock masses [J]. Tunnelling and Underground Space Technology, 2021, 113: 103968. DOI: 10.1016/j.tust.2021.103968.
    [22] WANG Z L, WANG H C, WANG J G, et al. Finite element analyses of constitutive models performance in the simulation of blast-induced rock cracks [J]. Computers and Geotechnics, 2021, 135: 104172. DOI: 10.1016/j.compgeo.2021.104172.
    [23] 俞茂宏. 双剪理论及其应用 [M]. 北京: 科学出版社, 1998: 849–860.
    [24] BORRVALL T, SWEDEN L, RIEDEL W. The RHT concrete model in LS-DYNA [C]// Proceedings of the 8th European LS-DYNA Conference. Strasbourg, France, 2011.
    [25] 李洪超, 陈勇, 刘殿书, 等. 岩石RHT模型主要参数敏感性及确定方法研究 [J]. 北京理工大学学报, 2018, 38(8): 779–785. DOI: 10.15918/j.tbit1001-0645.2018.08.002.

    LI H C, CHEN Y, LIU D S, et al. Sensitivity analysis determination and optimization of rock RHT parameters [J]. Transactions of Beijing Institute of Technology, 2018, 38(8): 779–785. DOI: 10.15918/j.tbit1001-0645.2018.08.002.
    [26] MA S Z, LIU K W, GUO T F, et al. Experimental and numerical investigation on the mechanical characteristics and failure mechanism of cracked coal & rock-like combined sample under uniaxial compression [J]. Theoretical and Applied Fracture Mechanics, 2022, 122: 103583. DOI: 10.1016/j.tafmec.2022.103583.
    [27] BANADAKI M M D, MOHANTY B. Numerical simulation of stress wave induced fractures in rock [J]. International Journal of Impact Engineering, 2012, 40/41: 16–25. DOI: 10.1016/j.ijimpeng.2011.08.010.
    [28] YANG J C, LIU K W, LI X D, et al. Stress initialization methods for dynamic numerical simulation of rock mass with high in-situ stress [J]. Journal of Central South University, 2020, 27(10): 3149–3162. DOI: 10.1007/s11771-020-4535-3.
    [29] 谢和平, 高峰, 周宏伟, 等. 岩石断裂和破碎的分形研究 [J]. 防灾减灾工程学报, 2003, 23(4): 1–9. DOI: 10.3969/j.issn.1672-2132.2003.04.001.

    XIE H P, GAO F, ZHOU H W, et al. Fractal fracture and fragmentation in rocks [J]. Journal of Disaster Prevention and Mitigation Engineering, 2003, 23(4): 1–9. DOI: 10.3969/j.issn.1672-2132.2003.04.001.
    [30] DING C X, YANG R S, YANG L Y. Experimental results of blast-induced cracking fractal characteristics and propagation behavior in deep rock mass [J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 142: 104772. DOI: 10.1016/J.ijrmms.2021.104772.
  • 加载中
图(28) / 表(5)
计量
  • 文章访问数:  264
  • HTML全文浏览量:  59
  • PDF下载量:  106
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-04-26
  • 修回日期:  2023-06-26
  • 刊出日期:  2023-10-27

目录

/

返回文章
返回