Blast-induced damage characteristics and fracture mechanism of rock mass under initial stress
-
摘要: 钻爆法在岩体开挖工程中应用广泛,其工作效率与岩体结构及地质环境密切相关。初始应力会显著影响岩体爆破裂纹的扩展行为及破坏特征,导致深部岩体超/欠挖等问题。本文采用理论分析和数值模拟相结合的方法研究了不同初始应力下岩体的爆破损伤特性及破裂机理。基于弹性力学建立单孔爆破动静组合理论模型,分析了静态应力分布及动态应力演化特征,揭示了初始应力作用下岩体爆破损伤机制。通过经验公式及动态力学试验对Riedel-Hiermaier-Thoma (RHT) 模型参数进行标定与修正,并结合室内爆破实验及理论结果对数值模型进行验证。此外,在单孔爆破数值模型中分析了爆炸压力演化规律、裂纹扩展行为及分形特征,结果表明数值模拟结果可以较好地论证理论模型的合理性。研究发现:环向拉应力是影响爆破裂纹扩展的重要因素,当初始应力较大时,合理调整环向应力分布可改善岩体爆破碎裂效果,实际工程中可结合预裂技术进行地压调控改变应力分布。Abstract: The technique of drilling and blasting is widely applied in rock excavation, and its working efficiency mainly relies on the rock mass structure and geological conditions. The initial stress plays a crucial role in the blast-induced cracking behavior and failure characteristics, and troubles such as over/underbreak and insufficient fragmentation may be arisen in deep rock masses, leading to other issues related to the waste of mineral resources and production costs. In this paper, the damage features and fracture mechanism of rock blasting under various initial stresses were studied using combined theoretical analysis and numerical simulation. Considering the effect of initial stress on the mechanical response of rock blasting, a theoretical model for single-hole blasting was developed based on elastic mechanics, and then the features of static stress distribution and dynamic pressure evolution were analyzed separately, revealing the fracture mechanism of rock blast-induced damage under initial stress. In addition, the Riedel-Hiermaier-Thoma (RHT) model parameters of rock were determined and adjusted based on a series of empirical formulas as well as dynamic mechanical tests. After calibrating the numerical model against the blasting crack pattern and attenuation of peak pressure by combining test and theoretical results, the single-hole geometric model was created and meshed in the commercial software ANSYS to study the dynamic tangential stress evolution as well as cracking behavior, and the fractal features of rock blasting cracks under different initial stresses were also discussed. The results show that the rationality of the theoretical model can be proved by the numerical simulation: the tangential tensile stress is a critical factor affecting the blasting crack propagation and the rock fragmentation can be improved by adjusting the distribution of hoop stress reasonably when the initial stress is large. In practical engineering, the stress distribution can be changed by controlling the field of geo-stress with presplitting technology.
-
Key words:
- initial stress /
- rock blasting /
- theoretical model /
- numerical simulation /
- crack propagation /
- fractal dimension
-
水中兵器对舰船的破坏程度是由很多因素决定的,主要包括主装药的性质和质量、爆心与舰船的距离和位置、舰船的结构和性质等。在评价水中兵器的破坏威力时,引入了冲击因子Q这样一个衡量标准[1],冲击因子的物理意义是,对于同一舰船,若冲击因子相等,则认为其水下爆炸的冲击响应近似相等。以前,使用一种基于冲击波超压的冲击因子考核潜艇结构生命力,这也是冲击因子的最早形式。由大量的实验和数值计算可发现,该形式的冲击因子不能很好地反映结构的破坏程度。随着对水下爆炸现象的进一步研究,越来越倾向于使用另一种基于平面波假定的冲击因子,这种冲击因子是从作用在结构上的冲击波能的角度定义的,远场时具有较好的效果,但该冲击因子没有考虑近场冲击波形状和能量损耗。姚熊亮等[2]利用球面波理论对水中爆炸冲击因子进行了修正,考虑了冲击波形状对垂直投射到舰船结构的冲击波能的影响。但是,水中爆炸近场的能量损耗非常严重,如48%冲击波能损耗在25个装药半径的范围内[3],而且不同种类炸药的能量释放特性也存在一定的差异。然而,传统冲击因子中描述炸药性能的TNT当量只是总化学能的一个比,并不能完全表征炸药水中爆炸的能量输出与衰减特性,尤其是水中爆炸近场。本文中,通过冲击波峰值压力和冲量相似方程的乘积推导一种适用于水中爆炸冲击波因子的计算方法,从冲击波毁伤作用角度描述该计算方法的物理意义,通过水中爆炸实验进行验证,并与传统的水下爆炸冲击因子进行对比。
1. 基于平面波的水中爆炸冲击因子
水中爆炸冲击波的毁伤作用可用冲击因子表示,基于平面波的冲击因子的表达式为[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为结构在垂直于冲击波阵面上的投影面积,m2;ρe为炸药的质量化学能,MJ/kg;ηe为炸药化学能转化为冲击波能的比例。
由式(3)可知,当水下爆炸冲击波为平面波时,无论目标距爆心的距离大小如何,Se为常数。在爆炸远场,匀化冲击波可近似为平面波,而且冲击波能基本恒定,Q为Es的函数。在爆炸近场,冲击波能在不断衰减,冲击波的形状不能近似为平面波,此时,Q为ηe、Es的函数,在某一距离处,ηe与炸药的性能有关,Es与冲击波的形状和目标结构有关。因此,炸药的TNT当量并不能完全表征近场冲击波的特性,基于平面波的水中爆炸冲击因子具有一定局限性。而通过水中爆炸实验得到的冲击波相似方程包含了ηe的特性,反映了炸药水下爆炸的能量释放和衰减特性,下面以冲击波相似方程为基础对冲击因子进行推导。
2. 冲击因子的推导
2.1 水中爆炸相似方程
水中爆炸的相似方程为[3]:
Y=K(W1/3/R)α 式中:W 为炸药的质量,kg;R为测点距爆心的距离,m;K和α为与炸药相关的系数,K为相似常数,α为相似指数;冲击波参数Y 包括峰值压力pmax、比时间常数θ/W1/3、比冲量I/W1/3和比能流密度E/W1/3。
2.2 压力和冲量乘积推导
由水中爆炸实验得到的峰值压力的经验计算公式为:
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/(αp+αI),A=(KpKI)m,n=其中m和A为冲击波参数的影响因子,n为装药指数。
通过水中爆炸实验得到的冲击波相似系数[7],计算了几种典型炸药的装药指数n,见表1。
表 1 几种典型炸药的相似系数和装药指数Table 1. Exponentsα,K,nfor various high explosives由表1可看出,在实验精度许可范围内,所有以TNT为基本组分的炸药,式(9)中装药指数n均可近似等于0.5。因此,冲击因子可表示为:
Q2=A(pmaxI+)m=W0.5R 式(10)将峰值压力pmax冲量I的毁伤作用双曲线与冲击因子的R-W关系曲线联系起来,从冲击波的毁伤作用角度表述了水中爆炸冲击因子的物理意义。
3. 实 验
3.1 样 品
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 heat3.2 装 置
水池∅3.2m×2.6m,水深2.4m,池底和池壁均由8mm钢板焊接而成。入水深度1.6m,为总水深的2/3,满足冲击波和气泡的测试要求,可以消除边界效应的影响[8]。爆心与传感器的距离分别为23、32、44、58和90cm,测量的是炸药柱中心轴向的径向冲击波,实验布局如图1所示。
测试系统包括138系列ICP型压电式电气石水下激波传感器和482A型信号适配器、高低频数据记录仪。冲击波信号的采样频率为10MHz。
4. 结果与分析
4.1 装药指数
通过实验测量的冲击波参数,拟合得到了不同组分炸药的相似系数和装药指数,见表3。
由表3可知,由pmaxI+推导的冲击因子计算的理想炸药、RDX基和HMX基含铝炸药的装药指数n都非常接近于0.5,平行性也非常好。可见,基于平面波假设的冲击因子的装药指数同样适合于RDX、HMX基的含铝炸药。
表 3 不同组分炸药的相似系数和装药指数Table 3. Exponentsα,K,nfor various explosives compositions4.2 冲击因子
为了对比基于平面波的水中爆炸冲击因子与由峰值压力与冲量的乘积pmaxI+推导的冲击因子计算公式的差异,依据实验数据计算了测距23、58和90cm的冲击因子,3个距离的水中爆炸冲击波参数见表4,不同距离处的冲击因子见表5和图2。
表 4 水中爆炸冲击波峰值压力和冲量Table 4. Shock wave peak pressure and impluse of underwater explosion由表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 distancesQ2 =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)。
5. 结 论
(1)水中爆炸冲击因子装药指数n=0.5,不仅适合所有以TNT为基本组分的炸药,也适合于RDX、HMX基的含铝炸药。
(2)由峰值压力与冲量的乘积pmaxI+推导的冲击因子计算公式,从冲击波的毁伤作用的角度表述了水中爆炸冲击因子的物理意义,结合冲击形状对冲击波能的影响,计算近场冲击因子时具有更高的准确性。
-
表 1 不同围压下花岗岩力学参数
Table 1. Mechanical parameters of granite sample under various confining pressure
σ2/MPa σ3/MPa σ1/MPa p∗0 σ∗f σ2/MPa σ3/MPa σ1/MPa p∗0 σ∗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 表 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 相对抗剪强度 F∗s 0.18 试验优化 孔隙坍塌压力 pcrush 108 相对抗拉强度 F∗t 0.06 洛德角相关因子 Q0 0.68 压缩屈服面参数 G∗c 0.50 洛德角相关因子 B 0.05 拉伸屈服面参数 G∗t 0.70 破坏面参数 A 2.48 残余面参数 Af 1.62 破坏面参数 N 0.79 残余面参数 Nf 0.62 弹性剪切模量 G 21.9 GPa 表 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 表 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 表 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 -
[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. -