Optimization and Experimental Study of Pre-Splitting Blasting Parameters in a Certain Open-Pit Mine
-
摘要: 为了得到确保露天矿山边坡安全稳定的最佳预裂爆破参数组合,采用数值模拟结合相似模型试验的方法,探讨了岩石双孔爆破下裂纹扩展规律,并进行了现场预裂爆破试验。结果表明:裂纹先由炮孔壁产生,且炮孔附近的原岩破坏程度较大;随着不耦合系数由1.33增大至3.00,孔壁所承受的冲击压力显著减小,粉碎区半径逐渐减小。采用极差分析法,研究了不同因素对试验指标的影响程度,根据评价指标确定试验因素的最优组合。模型试验结果表明:孔距与孔径之比为9、不耦合系数为3.00、起爆延期时间为12 ms、最大单响药量为5.4 g(三孔一响)的组合条件下,模型试件预裂成缝效果最佳。根据数值模拟及相似模型试验结果,结合理论计算,选取孔距为0.8 m(孔距与孔径之比约为9)、不耦合系数为3.00、延期时间为12 ms、最大单响药量为11.25 kg(三孔一响)的预裂爆破参数进行现场试验,取得了良好的预裂效果。研究结果可以为露天矿山现场预裂爆破参数的设计提供参考,对于实现精确控制爆破、确保露天矿山边坡安全稳定具有重要意义。Abstract: In order to obtain the optimal combination of pre-splitting blasting parameters to ensure the safety and stability of open-pit mine slope, the law of crack propagation under rock double-hole blasting was discussed by numerical simulation combined with similar model experiment, and the field pre-splitting blasting experiment was carried out. The results show that the cracks are first generated by the blasthole wall and the original rock near the blasthole is more damaged. As the decoupling coefficient increases from 1.33 to 3.00, the impact pressure on the hole wall decreases significantly, and the radius of the crushing zone decreases gradually. The range analysis method was used to study the influence of different factors on the experiment index, and the optimal combination of experiment factors was determined according to the evaluation index. The results of the model experiment indicate that the pre-splitting effect of the model specimen is the best under the combined conditions of the hole spacing-to-aperture ratio of 9, the decoupling coefficient of 3.00, the initiation delay time of 12 ms, and the maximum charge weight per delayed interval of 5.4 g (three holes and one shot). According to the results of numerical simulation and similar model experiment, combined with theoretical calculation, the pre-splitting blasting parameters with hole spacing of 0.8 m (the spacing-to-aperture ratio of about 9), decoupling coefficient of 3.00, delay time of 12 ms and maximum charge weight per delayed interval of 11.25 kg (three holes and one shot) were selected for field experiment, and good pre-splitting effect was achieved. The research results can provide reference for the design of pre-splitting blasting parameters in open-pit mines, and it is of great significance to realize accurate control blasting and ensure the safety and stability of open-pit mine slopes.
-
1. 引言
自1996年German-French Research Institute Saint Louis (ISL)实验室第一次提出横向效应增强型侵彻体(Penetrator with Enhanced Lateral Effect, PELE)概念以来, 国内外学者对PELE弹丸的作用过程进行了广泛研究。在PELE弹丸壳体断裂的机理方面, Kesberg等人[1]认为壳体侵彻目标的内核因密度较低而被压缩并沿径向膨胀, 导致壳体内壁压力剧增至吉帕级, 弹丸穿透目标后能量瞬间释放, 导致壳体破裂成大量破片。Paulus等人[2]结合Recht-Ipson模型[3]和Mott模型[4]计算出弹体的剩余速度、装填物的压力以及壳体破片的径向飞散速度, 并且用弱激波理论描述了PELE弹丸壳体的断裂过程。对于PELE弹丸壳体断裂的研究, 除了可从力学和能量的角度定性地解释PELE弹丸壳体断裂过程外, 还可以从损伤断裂的角度, 具体分析其在侵彻过程中的断裂机理。
在研究材料动态拉伸断裂方面有许多模型, Seaman、Curran和Shockey[5]从实验出发, 基于微损伤的统计性分析结果, 提出了成核与长大模型(Nucleation and Growth, NAG)。但是NAG模型计算需要大量的实验参数, 实用性比较差。封加波等人[6-7]在NAG模型的基础上利用单孔洞近似, 从能量平衡的角度出发, 提出了简单实用的损伤度函数模型。本研究基于损伤度模型, 结合PELE弹丸侵彻过程, 将PELE弹丸壳体的膨胀断裂过程分为加速膨胀和匀速膨胀两个阶段, 推导出这两个阶段PELE弹丸壳体的损伤断裂方程, 并根据损伤断裂方程分析影响PELE弹丸横向效应的内在机理。
2. 金属壳体动态损伤模型
2.1 金属壳体的动态本构关系
PELE弹丸壳体一般为高密度金属或合金。金属或合金材料在高应变率下的变形过程比较复杂, 理论上通常把动态变形过程划分为几个阶段, 如弹性阶段、塑性阶段和失效阶段, 对不同变形阶段需要用不同的本构关系表达式进行描述。动态塑性本构方程可以写为
σ=φ(ep,˙ep) (1) 式中:σ为材料的动态屈服应力, ep为径向塑性自然应变,
为径向塑性自然应变率。
根据文献[7], 金属壳体的动态本构关系可以表示为
σ=σ0+Kep+η˙ep (2) 式中:σ0为静态屈服强度, K为材料强化系数, η为材料黏度系数。
对于柱壳结构有
ep=ln(R(t)R0) (3) ˙ep=u(t)/R(t) (4) εp=R(t)−R0R0 (5) 式中:R0为壳体初始半径, R(t)为壳体实时半径, R0、R(t)均是壳体内径和外径的平均值, εp为径向塑性应变, u(t)为壳体径向膨胀速度。结合(3)式和(5)式可以得到ep与εp的关系
ep=ln(1+εp) (6) 将(6)式代入(2)式, 可以将动态本构关系改写成以下形式
σ=σ0+Kln(1+εP)+η˙ep (7) 2.2 金属壳体的动态损伤模型
本研究所采用的断裂准则建立在封加波提出的损伤度模型的基础上。金属材料的损伤度包括在应力波作用下材料初始损伤缺陷及其发展, 以及新产生的损伤缺陷及其发展, 而且材料的损伤度具有累积效应。在PELE弹丸壳体动态断裂的过程中, 为了考虑材料的损伤累积效应, 引入损伤变量D(t), D(t)定义为材料某一截面上微孔洞所占的面积比。参照文献[6], 有如下形式的方程
dDdt=[(1−D)(σ−σ0)]2ηK (8) 求解方程(8)式需要知道σ(或εp和
随时间变化的规律。从(4)式和(5)式看出, 需要知道壳体膨胀的实时半径R(t)和壳体径向膨胀速度u(t)的函数形式。在处理实际问题时, 很难求出解析表达式, 一般需要结合实验结果求解。
3. PELE弹丸壳体膨胀规律
3.1 PELE作用机理
PELE弹丸是一种基于物理效应、不需要装填炸药和引信的新概念弹药, 主要用于防空、反导、城市作战等领域。PELE弹丸由壳体和内核组成, 壳体一般为密度较大的金属材料, 内核是密度及强度比较低的金属或聚合物。高密度金属壳体对靶板具有良好的侵彻性能, 而内核材料在侵彻过程中被封闭压缩, 巨大的压力在内核中迅速积聚, 弹丸穿透靶板后能量瞬间释放导致外壳破裂解体。
PELE的横向增强效应作用过程分为3个主要时段和一个后效作用时段[1]。如图 1所示, 第1时段, 当弹丸击中目标后, 外层高密度壳体侵彻目标, 内核因材料强度低而无法侵彻, 将会在靶前停滞而被压缩; 第2时段, 封闭压缩过程中产生的高压将使内核材料沿径向膨胀, 导致壳体内壁压力急剧增加, 可达到数吉帕; 第3时段, 当应力应变达到壳体的破坏极限时, 出靶后的卸载应力波致使外层弹体破裂成大小不同的破片。后效作用时段则是壳体破片的径向飞散以及内核和塞块间继续相互作用的过程。
3.2 PELE侵彻过程数值计算
图 2给出了PELE弹丸计算模型, 其中壳体材料为35CrMnSi, 长度为105 mm, 外径为30 mm, 内径为18 mm, 弹底部厚为15 mm; 内核材料为聚乙烯, 长度为90 mm, 直径为18 mm; 靶板材料为装甲钢, 厚度为8 mm。网格平均尺寸为壳体1 mm, 内核1.2 mm, 靶板采用变步长网格(0.5~10 mm)以提高计算效率。为了获得侵彻过程中的时间历程参数, 在弹丸模型的壳体中心线沿轴向均匀分布4个观测点(Gauge点), 每个观测点间隔20 mm, 如图 2所示。
表 1、表 2分别为材料模型和材料参数。数值计算采用拉格朗日方法, 材料均采用侵蚀算法, 靶板采用材料应变失效模型, 壳体使用有效应力失效模型。采用Stochastic模型[8]模拟壳体随机断裂, 其中材料弱化点的破坏概率在[0, 1]范围为
表 1 材料模型Table 1. Material modelComponent Material species Equation of state Strength model Invalidation model Shell 35CrMnSi Linear Johnson-Cook Principal stress Core Polyethylene Shock von-Mises - Target plate RHA Shock von-Mises Material strain 表 2 材料参数Table 2. Material parametersMaterial ρ/(g/cm3) E/(GPa) σ/(MPa) G/(GPa) ν 35CrMnSi 7.8 210 1 270 81 0.33 Polyethylene 0.94 1.5 0.26 0.55 0.45 RHA 7.8 171 1400 64.1 0.33 P=1−exp(−ceγε/γ) (9) 式中:P为应变为ε时的网格破坏概率, c、γ取决于材料性质, 根据文献[8], 取c=0.467, γ=10。
PELE弹丸以700 m/s的着靶速度垂直侵彻8 mm厚的装甲钢靶板。壳体中Gauge点沿壳体径向速度变化曲线如图 3所示。从图 3中看出, PELE侵彻过程第2时段, 壳体径向先经历短暂加速膨胀阶段, 膨胀速度达到一个较大值; PELE侵彻过程第3时段, 壳体膨胀速度小幅震荡并趋于稳定。
3.3 实验过程与结果
PELE实验弹丸参数:外径为30 mm, 内径为18 mm, 长度为105 mm; 壳体材料采用35CrMnSi合金钢, 内核为聚乙烯, 长度为90 mm, 靶板为装甲钢, 厚度为8 mm; 装配全弹质量为350 g。
实验弹采用轻气炮发射。弹丸着靶和出靶分别设置了测速点探针, 以测量着靶速度和弹丸穿透靶板的时间。为了获得弹丸头部出靶时的膨胀速度, 设计了一套环形电探针测量装置。实验装置示意图如图 4所示, 实验装置各部分如图 5所示。
表 3分别统计了3发有效实验的弹丸着靶速度、弹丸头部从着靶到出靶的穿靶时间, 以及出靶时弹丸头部的径向速度。应该指出的是, 由于弹丸穿靶时间极短, 弹丸出靶时头部的径向速度比较大, 探针测试存在一定的误差。
表 3 着靶速度与穿靶时间统计Table 3. Statistics of the impact velocity and the penetrating timeExperimental No. Velocity/(m/s) Time/(μs) Radical velocity/(m/s) 1 643 8.4 127 2 694 7.2 142 3 762 6.5 153 3.4 PELE侵彻过程壳体的膨胀规律
壳体膨胀断裂主要处于PELE侵彻过程的第2时段和第3时段。壳体动态断裂过程可以描述为:应力波从内核传入壳体内壁, 壳壁开始加速膨胀, 在应力波多次反射效应下膨胀速度振荡衰减, 直至最终趋于平稳。其衰减规律基本上按余弦函数变化[9-10], 所以径向速度公式可以写为
u(t)={Ctt0⩽t<t1Ae−λtcos(ψt)+u0t1⩽t<tc (10) 式中:A、ψ、C、λ为与总体模型有关的参数, u0为稳态时的径向速度值。t0为应力波传到壳体外壁的时刻, t1为壳体膨胀加速结束的时刻, tc为壳体发生贯穿性断裂的时刻。
根据侵彻过程壳体的膨胀规律, 可以将PELE壳体断裂过程分为壳体加速膨胀阶段和壳体匀速膨胀阶段。调整数值模拟中弹丸的着靶速度, 使其工况与3发实验相吻合, 对比数值模拟结果与实验结果发现, PELE壳体加速膨胀的时间与壳体从着靶到出靶的时间相吻合。可以看出当壳体处于靶板内时, 壳体处于加速状态; 当壳体出靶后, 壳体匀速膨胀。
4. PELE弹丸壳体断裂过程求解
根据PELE侵彻过程壳体的膨胀规律, 将PELE壳体断裂过程分为两个阶段:壳体加速膨胀阶段和壳体匀速膨胀阶段。
第1阶段:壳体加速膨胀阶段。假定自然应变率可以近似看作常数
˙ep=u(t)/R(t)=c (11) 另外
dt=dep/˙ep (12) 对(12)式积分得
t1−t0=ep1/˙ep (13) 将(7)代入(8), 在[t0, t1]内积分得
D1=1−3ω1(ω1+ep1)3+3ω1−ω31 (14) 式中:D1为t1时刻的损伤度。(14)式中
ω1=(K/η)˙ep (15) 通过实验结果可以得到
代入(13)式可以求出
将
代入(14)式求出D1。
第2阶段:壳体匀速膨胀阶段。第1阶段结束时刻即为第2阶段初始时刻。依据PELE壳体膨胀规律, 给出描述壳体膨胀速度的变化规律为
u(t)=Ae−λtcos(ψt)+u0 (16) (16) 式反映应力波多次反射引起的衰减振荡, 代表了壳体膨胀速度的真实情况, 但直接代入损伤度方程难以求解。为了方便积分, 这里按照能量等效原理取均值, 假设柱壳膨胀速度u(t)等于常数
u(t)=u0 (17) 由于
˙ep=u(t)/R(t) (18) ˙ep2=u(t)/R2 (19) u(t)=dR(t)/dt (20) 将(5)式、(18)式、(19)式代入(20)式, 得
dt=dR(t)/u(t)=dεp/˙ep2 (21) 式中:
为第2阶段的初始径向塑性自然应变率, R2为第2阶段壳体的初始半径。
将(7)式、(17)式、(21)式代入(8)式, 在[t1, tc]积分得
∫DcD1dD(1−D)2=∫εpcεp11ω2[ln(1+εp)+ω21+εp]dεp (22) 式中:
为壳体发生贯穿性破坏的损伤度。经过整理, 最终得
εpc−εp1(1+εpc)(1+εp1)ω22+[(epc)2−(ep1)2−Dc−D1(1−Dc)(1−D1)]ω2+(1+εpc)[(epc)2−2epc+2]−(1+εp1)[(ep1)2−2ep1+2]=0 (23) 式中:
为第一阶段壳体材料径向塑形应变。由(6)式知
与
满足以下关系式
ep1=ln(1+εp1) (24) (23) 式为PELE弹丸壳体的损伤断裂方程。第1阶段的自然应变率
和第2阶段的初始半径R2、径向速度u0可以通过实验得到, 通过R2和u0可以求出
而壳体材料参数σ0、K、η和Dc均可以通过查找手册和文献获得, 最终可以得到一个确定的临界应变
作为损伤断裂判据。
将(17)式、(19)式代入(21)式, 在[t1, tc]内积分可得壳体发生贯穿断裂的时刻tc
tc−t1=R2(εpc−εp1)/u0 (25) 5. PELE壳体断裂实例
对于35CrMnSi, 参照文献[10], 壳体特征参数取η=-12 kPa·s, K=1.66 GPa, ωλ=3.5 kJ/m2, D0=10-5, Dc=0.2, σ0=0.72 GPa。当初始半径R0=15 mm的PELE弹丸以1 km/s速度着靶时, 结合实验数据, 得到第1阶段的自然应变率
和第2阶段的初始半径R2及径向膨胀速度u0, 见表 4。
表 4 3个时刻膨胀断裂过程的参数Table 4. Parameters of expansion and fracture at three time pointst/(μs) ˙ep/(104 s−1) R2/(mm) u0/(m/s) 8.4 2.62 17.7 127 7.2 1.91 18.6 142 6.5 1.71 19.9 153 将第2阶段的初始半径R2及径向膨胀速度u0代入(19)式, 得到第2阶段的初始应变率
将第一阶段
代入(13)式可以求出
将
代入(14)式求出D1; 同时将
代入(24)式求出
将
代入(23)式, 得到3个典型时刻临界应变
如图 6所示。从图 6中可以看出, 在一定速度范围内, 随着弹丸着靶速度的增加, 壳体断裂的临界应变明显降低。
6. 结论
(1) PELE弹丸在侵彻过程中壳体膨胀断裂可以分为加速膨胀和匀速膨胀两个阶段, 在封加波提出的损伤度模型的基础上, 得到了PELE弹丸壳体的损伤断裂方程和断裂判据。
(2) 根据PELE弹丸壳体损伤断裂方程, 计算出PELE弹丸壳体膨胀断裂过程的参数, 分析了弹丸着靶速度与壳体断裂之间的关系, 可为PELE弹丸设计提供理论指导。
致谢: 感谢胡玉涛博士和曹雷博士在公式推导和实验方面给予的帮助。 -
表 1 大理岩参数
Table 1. Marble parameters
ρ/(kg·m−3) E/GPa ν σbc/MPa Rm/MPa G/GPa K/GPa Cρρr/(kg·m−2·s−1) 2700 55.8 0.31 74.7 8.5 22.3 37.2 3.6×106 表 2 岩石RHT模型参数
Table 2. Rock RHT model parameters
fc/GPa N α0 pel/GPa βc βt A1/GPa A2/GPa A3/GPa B0 B1 0.1678 4.00 1.0 0.016 0.0244 0.0294 45.39 40.85 41.98 0.9 0.9 T1/GPa T2/GPa ˙εc0/s−1 ˙εt0/s−1 ˙εc/s−1 ˙εt/s−1 D1 D2 B g*t g*c 45.39 0 3.0×10–5 3.0×10–6 3.0×1025 3.0×1025 0.037 1.0 0.0105 0.7 0.78 A n f*s f*t Q0 ζ εmp Af nf pcomp/MPa 2.51 0.72 0.21 0.04 0.68 0.5 0.015 0.25 0.62 8.00 表 3 炸药材料及JWL状态方程参数
Table 3. Explosive material and the JWL equation of state parameters
ρb/(kg·m−3) Cd/(m·s−1) Ae/GPa Be/GPa R1 R2 ω 1100 4500 625.3 23.29 5.25 1.6 0.28 表 4 模型材料指标
Table 4. Model material indicators
ρm/(kg·m−3) σm/MPa Em/GPa Cρρr/(kg·m−2·s−1) 2680 30.02 28.1 5.06×106 表 5 量纲分析
Table 5. Dimensional analysis
ρ σbc E Cρρr H Cd ρb FL−4T2 FL−2 FL−2 FL−3T L LT−1 FL−4T2 表 6 正交试验方案设计
Table 6. Orthogonal experimental scheme design
Test No. k Δt/ms Q/g Error term Test No. k Δt/ms Q/g Error term 1 1.33 8 1.8 1 9 2.00 8 5.4 2 2 1.33 10 3.6 2 10 2.00 10 1.8* 1 3 1.33 12 5.4 3 11 2.00 12 1.8 4 4 1.33 15 1.8* 4 12 2.00 15 3.6 3 5 1.67 8 3.6 4 13 3.00 8 1.8* 3 6 1.67 10 1.8 3 14 3.00 10 5.4 4 7 1.67 12 1.8* 2 15 3.00 12 3.6 1 8 1.67 15 5.4 1 16 3.00 15 1.8 2 表 7 爆破效果统计
Table 7. Statistics of blasting effects
Test No. Evaluating indicator Test No. Evaluating indicator δ/% b/mm δ/% b/mm 1 45 51 9 65 16 2 48 63 10 76 15 3 51 80 11 70 10 4 41 54 12 65 12 5 56 41 13 76 12 6 30 45 14 90 21 7 78 62 15 87 17 8 62 51 16 86 15 表 8 正交试验极差分析(半孔率)
Table 8. Extreme analysis for orthogonal experiments (half porosity)
Factor ¯K1 ¯K2 ¯K3 ¯K4 R k 46.25 56.50 71.50 84.75 38.50 Δt 60.50 61.00 71.50 66.00 11.00 Q 57.75 66.50 67.00 67.75 10.00 表 9 正交试验极差分析(预裂缝宽度)
Table 9. Extreme analysis for orthogonal experiments (pre-split crack width)
Factor ¯K1 ¯K2 ¯K3 ¯K4 R k 62.00 49.75 13.25 16.25 48.75 Δt 30.00 36.00 42.25 33.00 12.25 Q 30.25 33.25 42.00 35.75 11.75 -
[1] 刘昭桃. 岩石爆破技术的现状与发展浅谈[C]//中国铁道学会工程分会. 爆破工程技术交流论文集. 北京: 中国铁道出版社, 2018: 85−87.LIU Z T. Status and development of rock blasting technology [C]//Engineering Branch of China Railway Society. Proceedings of Blasting Engineering Technology Exchange, Beijing: China Railway Publishing House, 2018: 85−87. [2] 许守信, 黄绍威, 李二宝, 等. 复杂破碎岩体矩形聚能药包预裂爆破试验研究 [J]. 金属矿山, 2021(11): 55–63.XU S X, HUANG S W, LI E B, et al. Experimental study on pre-splitting blasting of rectangular shaped charge in complex fractured rock mass [J]. Metal Mine, 2021(11): 55–63. [3] BENDEZU M, ROMANEL C, ROEHL D. Finite element analysis of blast-induced fracture propagation in hard rocks [J]. Computers and Structures, 2017, 182: 1–13. [4] 王和平, 郭连军, 张大宁, 等. 大孤山铁矿预裂爆破研究与应用 [J]. 金属矿山, 2015(10): 18–23.WANG H P, GUO L J, ZHANG D N, et al. Application and research of pre-splitting blasting in dagushan iron mine [J]. Metal Mine, 2015(10): 18–23. [5] 杨仁树, 苏洪. 爆炸荷载下含预裂缝的裂纹扩展实验研究 [J]. 煤炭学报, 2019, 44(2): 482–489.YANG R S, SU H. Experimental study on crack propagation with pre-crack under explosion load [J]. Journal of China Coal Society, 2019, 44(2): 482–489. [6] 朱必勇, 焦文宇, 寇向宇, 等. 基于数值模拟的预裂爆破参数优化研究 [J]. 有色金属(矿山部分), 2019, 71(4): 32–36.ZHU B Y, JIAO W Y, KOU X Y, et al. Parameters optimization of pre-split blasting based on numerical simulation [J]. Non-Ferrous Metals (Mine Part), 2019, 71(4): 32–36. [7] MA J, LI X L, WANG J G, et al. Experimental study on vibration reduction technology of hole-by-hole presplitting blasting [J]. Geofluids, 2021, 2021: 1–10. [8] 叶海旺, 唐可, 万涛, 等. 时序控制预裂爆破参数优化及应用 [J]. 爆炸与冲击, 2017, 37(3): 502–509. doi: 10.11883/1001-1455(2017)03-0502-08YE H W, TANG K, WAN T, et al. Optimization of time sequence controlled pre-splitting blasting parameters and its application [J]. Explosion and Shock Waves, 2017, 37(3): 502–509. doi: 10.11883/1001-1455(2017)03-0502-08 [9] ZHANG B S, YANG Z P, YANG X M, et al. Presplit blasting technique in treating hard overlying strata: from numerical simulation to field practice [J]. Advances in Civil Engineering, 2021: 1–20. [10] 宫嘉辰, 陈士海. 隧道爆破力学模型相似材料配比的正交试验 [J]. 华侨大学学报(自然科学版), 2020, 41(2): 164–170.GONG J C, CHEN S H. Orthogonal test of similar materials ratio in tunnel blasting mechanical model [J]. Journal of Huaqiao University (Natural Science Edition), 2020, 41(2): 164–170. [11] 单仁亮, 黄宝龙, 蔚振廷, 等. 岩巷掘进准直眼掏槽爆破模型试验研究 [J]. 岩石力学与工程学报, 2012, 31(2): 256–264. doi: 10.3969/j.issn.1000-6915.2012.02.004SHAN R L, HUANG B L, WEI Z T, et al. Model test of quasi-parallel cut blasting in rock drivage [J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(2): 256–264. doi: 10.3969/j.issn.1000-6915.2012.02.004 [12] 王正煜. 凤凰山石灰岩矿区边坡控制爆破参数优化研究[D]. 太原: 太原理工大学, 2021.WANG Z Y. Study on optimization of controlled blasting parameters of slope in Fenghuangshan limestone mining area [D]. Taiyuan: Taiyuan University of Technology, 2021. [13] 胡启文. 精确延时数(单)孔一响预裂爆破试验研究[D]. 昆明: 昆明理工大学, 2018.HU Q W. Study on preburst test of exact delay number (single) hole [D]. Kunming: Kunming University of Science and Technology, 2018. [14] 蒲传金, 杨鑫, 肖定军, 等. 爆炸载荷下双孔裂纹扩展的数值模拟研究 [J]. 振动与冲击, 2022, 41(15): 300–311.PU C J, YANG X, XIAO D J, et al. Numerical simulation study of double-hole crack propagation under explosive load [J]. Journal of Vibration and Shock, 2022, 41(15): 300–311. [15] 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: 1–12. [16] 李洪超. 岩石RHT模型理论及主要参数确定方法研究[D]. 北京: 中国矿业大学(北京), 2016.LI H C. Study on rock RHT model theory and determination method of main parameters [D]. Beijing: China University of Mining and Technology (Beijing), 2016. [17] 赵星宇, 白春华, 姚箭,等. 燃料空气炸药爆轰产物JWL状态方程参数计算 [J]. 兵工学报, 2020, 41(10): 1921–1929.ZHAO X Y, BAI C H, YAO J, et al. Parameters calculation of JWL EOS of FAE detonation products [J]. Acta Armamentarii, 2020, 41(10): 1921–1929. [18] 卢文波, 耿祥, 陈明, 等. 深埋地下厂房开挖程序及轮廓爆破方式比选研究 [J]. 岩石力学与工程学报, 2011, 30(8): 1531–1539.LU W B, GENG X, CHEN M, et al. Study of selection of excavation procedure and contour blasting method for deep underground powerhouse [J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(8): 1531–1539. [19] 徐颖, 孟益平, 程玉生. 装药不耦合系数对爆破裂纹控制的试验研究 [J]. 岩石力学与工程学报, 2002(12): 1843–1847.XU Y, MENG Y P, CHENG Y S. Study on control of blast crack by decoupling charge index [J]. Chinese Journal of Rock Mechanics and Engineering, 2002(12): 1843–1847. [20] 何理, 钟冬望. 微差爆破地震波沿高程传播特性的试验研究 [J]. 化工矿物与加工, 2015, 44(5): 36–40. doi: 10.16283/j.cnki.hgkwyjg.2015.05.011HE L, ZHONG D W. Experimental study on propagation characteristics of millisecond blasting seismic waves along elevation [J]. Chemical Minerals and Processing, 2015, 44(5): 36–40. doi: 10.16283/j.cnki.hgkwyjg.2015.05.011 [21] 司剑峰, 钟冬望, 黄小武. 钻孔爆破孔间最佳延时时间模型试验研究 [J]. 金属矿山, 2015(6): 19–23.SI J F, ZHONG D W, HUANG X W. Experimental model of the optimal delay time in drilling blasting [J]. Metal Mine, 2015(6): 19–23. [22] 夏祥. 爆炸荷载作用下岩体损伤特征及安全阀值研究[D]. 武汉: 中国科学院研究生院, 2006.XIA X. Study on rock mass damage characteristics and safety threshold value under the action of explosion load [D]. Wuhan: Graduate School of Chinese Academy of Sciences, 2006. [23] YANG L J, YANG A Y, CHEN S Y, et al. Model experimental study on the effects of in situ stresses on pre-splitting blasting damage and strain development [J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 138: 1–9. [24] 李洪伟, 吴延梦, 吴立辉, 等. 电子雷管起爆条件下隧道掏槽孔与辅助孔的延时优化试验研究 [J]. 高压物理学报, 2023, 37(1): 015301. doi: 10.11858/gywlxb.20220638LI H W, WU Y M, WU L H, et al. Experimental study on delay time optimization of tunnel cutting holes and caving holes under electronic detonator initiation condition [J]. Chinese Journal of High Pressure Physics, 2023, 37(1): 015301. doi: 10.11858/gywlxb.20220638 -