An investigation of energy consumption variation in rock blasting breaking with the resistance line
-
摘要: 针对爆炸荷载下岩体破碎块度和有用功能耗及能耗利用率问题,运用断裂力学、分形基础理论分析和模型实验等方法,对爆炸荷载下岩体破碎块度和能耗利用率随最小抵抗线的变化规律开展了系统的分析研究。研究结果表明:在模型实验条件下,破碎块度分形维数在1.2~1.7之间,随最小抵抗线增大呈现较好的线性衰减趋势;破碎能耗随最小抵抗线呈现先增加后降低的趋势,爆炸能量利用率在4.57%~12.51%之间,趋势与能耗值一致,模型实验中能耗利用率存在最大值;破碎块度与能耗利用率变化趋势相反,说明在最小抵抗线变化过程中存在一个最佳值,使得破碎块度和能耗利用率均处于最优状态,模型实验中这一最佳值为160 mm,是装药直径的26.7倍。该研究结果可为提高爆炸能利用率理论分析及工程中降低大块率的设计和施工提供理论依据。Abstract: The topics such as fragmentation degree, active energy consumption and energy consumption efficiency of rock mass under explosive load have attracted increasing attention in recent years. However, it is very difficult to conduct such research due to its instantaneity, high temperature and high pressure characteristics. Systematic analysis and research on broken blocks of rock mass and the variation of energy utilization under explosion load with the different minimum resistance lines have been carried out. Plain concrete material was used to construct the model and carry out the model experiment. Theory of energy consumption in fracture mechanics was used to calculate the crushing energy. The basic fractal theory was used to calculate and analyze the fragmentation distribution law. Research results indicate that: the fractal dimension of broken blocks is between 1.2 and 1.7, exhibiting a good linear attenuation trend with the increase of the minimum resistance line from 120 mm to 200 mm; the crushing energy consumption increases first and then decreases. More specifically, the crushing energy consumption is 440.0 J at 120 mm, and increases to the maximum of 1 106.5 J at 180 mm, and then decreases to 1 084.8 J at 200 mm. The explosive energy utilization rate is between 4.57% and 12.51% and the maximum value corresponds to the minimum resistance line of 180 mm, the variation trend is consistent with that of the energy consumption value. The trend of broken fragmentation and utilization rate of energy consumption is opposite. There is an optimum minimum resistance line, corresponding to the optimized fragmentation degree and energy consumption utilization rate, in the model experiment is 160 mm which is 26.7 times the diameter of the charge. The research results can provide a theoretical basis for improving the explosive energy utilization rate and guide the design and construction processes in future engineering applications.
-
Key words:
- fragmentation /
- energy explosion /
- energy efficiency /
- 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. Physical and mechanical parameters of materials
密度/(kg∙m−3) 纵波波速/(m∙s−1) 泊松比 抗压强度/MPa 弹性模量/GPa 断裂韧性/(MPa∙m1/2) 1 850 2 326 0.235 8.38 10.02 1.50 表 2 炸药性能
Table 2. Explosive performance
炸药类别 线密度/(kg∙m−1) 爆热/(kJ∙kg−1) 爆力/mL 爆速/(m∙s−1) 黑索金 0.025 5600 480 8300 表 3 爆破参数
Table 3. Parameters of blasting
编号 孔深/mm 药卷直径/mm 药量/g 装药长度/mm 单耗/(kg·m−3) 最小抵抗线/mm 1 225 6 1.58 40 0.49 120 2 245 0.33 140 3 265 0.24 160 4 285 0.17 180 5 305 0.13 200 表 4 岩块累计质量占比
Table 4. Cumulative mass ratio of rock blocks
抵抗线/mm 爆堆质量/kg 不同尺寸破碎块体长度分布占比 <50 mm <70 mm <90 mm <110 mm <130 mm <150 mm 120 20.60 0.17 0.39 0.51 0.62 0.79 1 140 28.90 0.13 0.30 0.48 0.59 0.78 1 160 42.04 0.14 0.29 0.46 0.58 0.77 1 180 55.39 0.13 0.31 0.45 0.57 0.77 1 200 57.66 0.10 0.25 0.35 0.55 0.75 1 表 5 拟合函数
Table 5. Fitting function
抵抗线/mm 拟合函数 D R2 120 ln yi=1.38 ln(xi/12) 1.62 0.97 140 ln yi=1.49 ln(xi/12) 1.51 0.96 160 ln yi=1.55 ln(xi/12) 1.45 0.97 180 ln yi=1.59 ln(xi/12) 1.41 0.98 200 ln yi=1.72 ln(xi/12) 1.28 0.96 表 6 破碎能分布参数
Table 6. Parameter of fragmentation energy
抵抗线/
mm破碎质量/
kg原表面积/
m2岩块表面积/
m2新增表面积/
m2破碎能/
J破碎体积/
(10−3 m3)破碎能耗密度/
(kJ·m−3)破碎能利用率/
%120 20.60 0.15 1.13 0.98 440.0 11.13 39.51 4.97 140 28.90 0.19 1.50 1.31 589.5 15.62 37.73 6.66 160 42.04 0.24 2.13 1.89 848.7 22.72 37.35 9.59 180 55.39 0.29 2.75 2.46 1106.5 29.94 36.96 12.51 200 57.66 0.30 2.71 2.42 1084.8 31.16 34.80 12.26 -
[1] WHITTLES D N, KINGMAN S, LOWNDES I, et al. Laboratory and numerical investigation into the characteristics of rock fragmentation [J]. Minerals Engineering, 2006, 19(14): 1418–1429. DOI: 10.1016/j.mineng.2006.02.004. [2] HAMDI E, DU MOUZA J, FLEURISSON J A. Evaluation of the part of blasting energy used for rock mass fragmentation [J]. Fragblast, 2001, 5(3): 180–193. DOI: 10.1076/frag.5.3.180.7386. [3] GRADY D E. Length scales and size distributions in dynamic fragmentation [J]. International Journal of Fracture, 2010, 163(1−2): 85–99. DOI: 10.1007/s10704-009-9418-4. [4] DARYADEL S S, MANTENA P R, KIM K, et al. Dynamic response of glass under low-velocity impact and high strain-rate SHPB compression loading [J]. Journal of Non-Crystalline Solids, 2016, 432: 432–439. DOI: 10.1016/j.jnoncrysol.2015.10.043. [5] 杨仁树, 许鹏. 爆炸作用下介质损伤破坏的分形研究 [J]. 煤炭学报, 2017, 42(12): 3065–3071. DOI: 10.13225/j.cnki.jccs.2017.0107.YANG R S, XU P. Fractal study of media damage under blasting loading [J]. Journal of China Coal Society, 2017, 42(12): 3065–3071. DOI: 10.13225/j.cnki.jccs.2017.0107. [6] 杨仁树, 李炜煜, 杨国梁, 等. 炸药类型对富铁矿爆破效果影响的试验研究 [J]. 爆炸与冲击, 2020, 40(6): 065201. DOI: 10.11883/bzycj-2019-0396.YANG R S, LI W Y, YANG G L, et al. Experimental study on the blasting effects of rich-iron ore with different explosives [J]. Explosion and Shock Waves, 2020, 40(6): 065201. DOI: 10.11883/bzycj-2019-0396. [7] 李清, 王平虎, 杨仁树, 等. 切槽孔爆破动态力学特征的动焦散线实验 [J]. 爆炸与冲击, 2009, 29(4): 413–418. DOI: 10.11883/1001-1455(2009)04-0413-06.LI Q, WANG P H, YANG R S, et al. Experimental investigation on dynamic mechanical behaviors of cracks induced by V-notch borehole blasting with dynamic caustics [J]. Explosion and Shock Waves, 2009, 29(4): 413–418. DOI: 10.11883/1001-1455(2009)04-0413-06. [8] 吴亮, 卢文波, 宗琦. 岩石中柱状装药爆炸能量分布 [J]. 岩土力学, 2006, 27(5): 735–739. DOI: 10.3969/j.issn.1000-7598.2006.05.010.WU L, LU W B, ZONG Q. Distribution of explosive energy consumed by column charge in rock [J]. Rock and Soil Mechanics, 2006, 27(5): 735–739. DOI: 10.3969/j.issn.1000-7598.2006.05.010. [9] 冷振东, 卢文波, 范勇, 等. 侧向起爆条件下的爆炸能量分布及其对破岩效果的影响 [J]. 爆炸与冲击, 2017, 37(4): 661–669. DOI: 10.11883/1001-1455(2017)04-0661-09.LENG Z D, LU W B, FAN Y, et al. Explosion energy distribution by side initiation and its effects on rock fragmentation [J]. Explosion and Shock Waves, 2017, 37(4): 661–669. DOI: 10.11883/1001-1455(2017)04-0661-09. [10] 胡振中, 庄亚明, 蔡天意, 等. 单颗粒煤岩冲击破碎能耗与粒度分布特性试验研究 [J]. 煤炭学报, 2015, 40(S1): 230–234. DOI: 10.13225/j.cnki.jccs.2014.1179.HU Z Z, ZHUANG Y M, CAI T Y, et al. Experimental study on energy consumption and particle size distribution of single particle coal under impact crushing [J]. Journal of China Coal Society, 2015, 40(S1): 230–234. DOI: 10.13225/j.cnki.jccs.2014.1179. [11] 李祥龙, 王建国, 张智宇, 等. 应变率及节理倾角对岩石模拟材料动力特性的影响 [J]. 爆炸与冲击, 2016, 36(4): 483–490. DOI: 10.11883/1001-1455(2016)04-0483-08.LI X L, WANG J G, ZHANG Z Y, et al. Experimental study for effects of strain rates and joint angles on dynamic responses of simulated rock materials [J]. Explosion and Shock Waves, 2016, 36(4): 483–490. DOI: 10.11883/1001-1455(2016)04-0483-08. [12] 甘德清, 刘志义, 李占金, 等. 冲击载荷作用下磁铁矿石破碎能耗特征 [J]. 岩石力学与工程学报, 2018, 37(S1): 3500–3506. DOI: 10.13722/j.cnki.jrme.2016.1150.GAN D Q, LIU Z Y, LI Z J, et al. Broken energy dissipation characteristics of magnetite under impact loads [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(S1): 3500–3506. DOI: 10.13722/j.cnki.jrme.2016.1150. [13] 武仁杰, 李海波, 李晓锋, 等. 冲击载荷作用下层状岩石破碎能耗及块度特征 [J]. 煤炭学报, 2020, 45(3): 1053–1060. DOI: 10.13225/j.cnki.jccs.2019.0266.WU R J, LI H B, LI X F, et al. Broken energy dissipation and fragmentation characteristics of layered rock under impact loading [J]. Journal of China Coal Society, 2020, 45(3): 1053–1060. DOI: 10.13225/j.cnki.jccs.2019.0266. [14] 叶洲元, 李夕兵, 万国香, 等. 受三维静载压缩岩石对冲击能的吸收效应 [J]. 爆炸与冲击, 2009, 29(4): 419–424. DOI: 10.11883/1001-1455(2009)04-0419-06.YE Z Y, LI X B, WAN G X, et al. Impact energy-absorption property of rock under tri-axial compression [J]. Explosion and Shock Waves, 2009, 29(4): 419–424. DOI: 10.11883/1001-1455(2009)04-0419-06. [15] 祝文化, 明锋, 宋成梓. 爆破荷载作用下岩体损伤破坏的分形研究 [J]. 岩土力学, 2011, 32(10): 3131–3135. DOI: 10.3969/j.issn.1000-7598.2011.10.040.ZHU W H, MING F, SONG C Z. Fractal study of rock damage under blasting loading [J]. Rock and Soil Mechanics, 2011, 32(10): 3131–3135. DOI: 10.3969/j.issn.1000-7598.2011.10.040. [16] 于永江, 王来贵, 何峰. 煤体爆堆块度分布的测试 [J]. 煤炭学报, 2005, 30(3): 337–339. DOI: 10.3321/j.issn:0253-9993.2005.03.015.YU Y J, WANG L G, HE F. The fragmentation distribution testing of rock blasting [J]. Journal of China Coal Society, 2005, 30(3): 337–339. DOI: 10.3321/j.issn:0253-9993.2005.03.015. [17] BAO R H, ZHANG L C, YAO Q Y et al. Estimating the peak indentation force of the edge chipping of rocks using single point-attack pick [J]. Rock Mechanics and Rock Engineering, 2011, 44(3): 339–347. DOI: 10.1007/s00603-010-0133-2. [18] 许金余, 刘石. 大理岩冲击加载试验碎块的分形特征分析 [J]. 岩土力学, 2012, 33(11): 3225–3229.XU J Y, LIU S. Research on fractal characteristics of marble fragments subjected to impact loading [J]. Rock and Soil Mechanics, 2012, 33(11): 3225–3229. -