Coupling effect of fuel property parameters on gas/coal dust composite explosion
-
摘要: 为研究瓦斯煤尘复合爆炸影响因素的耦合规律，对煤粉质量浓度、甲烷体积分数、煤粉粒径、煤粉种类等4种影响因素进行了多因素与单因素实验分析。通过正交实验，将各因素对爆炸的影响进行了定量分析，结果表明，4个因素对最大爆炸压力pmax的影响由强到弱依次为：甲烷体积分数、煤粉质量浓度、煤粉种类、煤粉粒径；对最大爆炸压力上升速率(dp/dt)max的影响程度由强到弱依次为：甲烷体积分数、煤粉质量浓度、煤粉粒径、煤粉种类。对于体积分数为9%、11%的甲烷，复合体系的pmax随煤粉的质量增加而减小。当煤粉质量浓度增加到100、200 g/m3时，在与体积分数为6%的甲烷耦合作用下，会产生更强的“激励”作用，且煤粉浓度较大时，挥发分低的煤种最佳瓦斯浓度会降低。甲烷体积分数存在临界值，该临界值会改变挥发分因素的影响方式：低于此临界值时，高挥发分煤尘体系的(dp/dt)max更高，(dp/dt)max来临时间更短；高于此临界值时，低挥发分体系具有更高的爆炸强度。粒径影响挥发分的作用，粒径越大，挥发分的影响差异越明显。当甲烷体积分数为11%时，挥发分高的煤尘更容易受粒径的影响，直径小的煤尘体系，爆炸系数Kst更小；而低挥发分煤粉在甲烷体积分数接近当量时受粒径影响更明显。Abstract: In order to reveal the coupling law of the influencing factors in gas/coal dust composite explosion, experiments were conducted in a 20 L spherical explosion device. Multi-factor and single-factor experimental analyses were conducted on four influencing factors according to explosion parameters, including coal concentration, methane volume fraction, coal particle size, and coal type. By conducting orthogonal experiments on the influencing factors, and using the explosion parameters as indicators to carry out range analysis and variance analysis, the influence of each factor on the explosion had been quantified. The experimental results show that the degree of influence of the four factor on pmax are (from strong to weak): methane volume fraction, coal dust mass concentration, coal dust type, and coal dust particle size; the effects on (dp/dt)max are (from strong to weak): methane volume fraction, coal dust mass concentration, coal dust particle size, and coal dust type. For methane concentrations of 9% and 11%, the value of pmax of the composite system decreases with increasing mass of coal dust. When the mass concentration of coal dust increases to 100 g/m3 and 200 g/m3, coupled with 6% volume fraction of methane, it will produce a stronger “incentive” effect, and when the concentration of coal dust is larger, the low volatility will reduce the optimal gas concentration. The existence of a certain critical value of methane concentration will change the influence mode of volatile factors. Below this value, the value of (dp/dt)max of the high volatile coal dust system is higher, and its arrival time is shorter, while above this value, the low volatile fraction system has a higher explosion intensity. Particle size implicates the influence of volatile fraction, the larger the particle diameter, the more obvious the difference produced by the volatile fraction factor. At 11% methane concentration, coal dust with high volatile fraction is more susceptible to the effect of particle size, and the smaller diameter coal dust system has smaller Kst value; at methane concentration close to equivalent, low volatile fraction coal dust is more significantly affected by the particle size factor.
-
Key words:
- gas /
- coal dust /
- gas/dust explosion /
- influencing factor /
- coupling
-
据统计，现代战争中有70%以上的士兵伤亡是由破片和枪弹所致[1]。对第二次世界大战后的历次现代战争的战伤统计数据表明，80%以上的战伤是由手榴弹、迫击弹和其他爆炸物的破片所致[2]。大面积的杀伤破片已成为威胁士兵安全的主要因素。而自防弹衣问世以后，士兵伤亡率大幅降低。这对单兵武器的毁伤元提出了更高的要求。如何设计合适的破片使其贯穿防弹衣后对人体目标造成有效或致命伤害成为单兵破片战斗部设计的重点。
单兵破片战斗部要求体积小，威力大，机动性强，而钨合金材料密度大、强度高，球形小破片体积小、数量多、存速高，可有效提高杀伤威力。因此，在单兵战斗部上采用小质量钨合金球形破片是未来发展的方向之一，既能提高单位面积的破片数量又能保持良好的机动性能。目前，国内外关于防弹衣加人体等效靶的侵彻研究大多以标准枪弹为主，如Lidén等[3]研究了9 mm子弹侵彻带有软防护活体麻醉猪的非贯穿性损伤，结果表明非贯穿损伤能对软防护后的胸膛造成严重伤害；Roberts等[4-5]模拟计算了9 mm手枪弹侵彻带有软质防弹衣的拟人体上躯干，并与假人靶标试验结果进行了对比，分析了非贯穿性弹道冲击对人体躯干内部器官的影响；Merkle等[6]采用试验与数值模拟方法，研究了9 mm手枪弹侵彻带有软质防弹衣的假人体躯干，分析了人体内部器官在非贯穿损伤下的压力分布情况；Gilson等[7]采用试验和数值模拟方法，研究了9、10.9 mm枪弹非贯穿侵彻带软防护的明胶靶的弹道冲击响应；董萍等[8]利用CT扫描和MIMICS软件重建了人体躯干三维模型，并对9 mm手枪弹侵彻带防弹衣的拟人体躯干靶进行了数值模拟，研究了拟人体躯干主要脏器在弹头非贯穿冲击下的瞬态响应特性；韩瑞国等[9]通过试验和数值模拟研究了某步枪弹对带软、硬复合防护明胶靶标的侵彻机制；刘坤等[10]研究了9 mm全铜弹和5.8 mm手枪弹对带软防护明胶靶标的侵彻机理；除此之外，针对92a铅芯弹[11]、7.62 mm[12]和9 mm[13]枪弹分别侵彻带软防护明胶靶的动态力学特性也作了相应的研究。而关于小尺寸破片对带有防弹衣的人体等效靶的侵彻研究却鲜有报道。
本文采用国际惯用标准25 mm厚红松木靶作为人体等效靶，采用实验与数值模拟方法研究小钨球对三级软体防弹衣加25 mm红松靶的侵彻过程及破坏机理，探讨钨球质量变化对弹道极限及靶板能量吸收的影响。在此基础上，利用量纲分析的方法研究钨球的穿靶能量与初速及直径的关系，并建立钨球的穿靶能量公式及弹道极限公式，以期为单兵破片战斗部的设计提供参考。
1. 实验研究
1.1 实验布置
实验采用12.7 mm弹道枪发射置于尼龙弹托中的小钨球，通过调整药筒中装药量来控制钨球的速度。钨球质量为（0.210+0.003） g，直径为（2.80+0.02） mm，密度为18.1 g/cm3，其化学成分及力学性能如表1所示。靶板为FDY3R-01型三级软体防弹衣加25 mm厚红松木，两者紧密贴合组成防弹衣加红松木复合靶。防弹衣由衣套和防弹层构成，衣套为涤纶，防弹层为凯夫拉材料，共45层，总厚度9 mm。复合靶用专用夹具固定在钢靶架上。为测量着靶前钨球速度及穿透靶后的剩余速度，靶前及靶后分别设置梳状通断靶，测速装置采用南京理工大学机械工程学院设计的NLG202-Z型六路测速仪，精度为0.1 μs。图1为实验所用的钨球、弹托及药筒，图2为弹道实验示意图。本次实验的侵彻均为垂直侵彻。
表 1 钨球成分及其力学性能Table 1. Chemical composition and mechanical properties of tungsten sphereρ/（g·cm−3） 质量分数/% 变形量/% HRC硬度 Ra/μm W Ni Fe Co 18.1 95.30 3.15 1.35 0.20 ≤40 ≥26 ≤1.6 注：ρ为密度，Ra为表面粗糙度，变形量为7 kN下的变形量。 1.2 钨球侵彻防弹衣加松木靶的实验结果及分析
1.2.1 实验结果
0.21 g、直径2.8 mm的钨球侵彻防弹衣加红松木复合靶的实验结果如表2所示。图3为实验后的防弹衣加红松木复合靶，其中图3(a)为单独的防弹衣，图3(b)为防弹衣加红松木复合靶。
表 2 钨球侵彻防弹衣加红松木复合靶的实验结果Table 2. Experimental results of tungsten spheres penetrating into pine target covered with body armor序号 着靶速度/（m·s−1） 剩余速度/（m·s−1） 结果 1 681.2 − 穿透防弹衣，嵌入松木靶 2 700.4 74.8 穿透防弹衣和松木靶 3 711.7 111.6 穿透防弹衣和松木靶 4 725.5 160.4 穿透防弹衣和松木靶 5 744.5 194.2 穿透防弹衣和松木靶 6 748.0 204.9 穿透防弹衣和松木靶 7 753.4 220.9 穿透防弹衣和松木靶 8 775.7 250.6 穿透防弹衣和松木靶 1.2.2 实验结果处理与分析
本文利用Recht等[14]提出的R-I公式以获得钨球侵彻防弹衣加红松木复合靶的弹道极限，R-I公式为：
vr=a(vpi−vpbl)(1/p) (1) 式中：vi为破片着靶速度，m/s；vr为破片剩余速度，m/s；vbl为弹道极限，m/s；a、p为模型参数。
依据文献[15]，p值设定为2。a和vbl可根据实验数据通过最小二乘法拟合得到。图4给出了钨球侵彻防弹衣加红松木复合靶的剩余速度与着靶速度关系曲线。由图4可看出，通过式(1)拟合出的关系曲线拟合效果较好。表3给出了通过式(1)计算所得的弹道极限及模型参数。
图 4 钨球剩余速度-着靶速度曲线Figure 4. Curve of residual velocity-impact velocity of tungsten spheres表 3 弹道极限及模型参数Table 3. Ballistic limit and model parameters钨球规格 a p vbl/（m·s−1） 质量0.21 g/直径2.8 mm 0.73 2 692.9 2. 数值模拟研究
2.1 数值模型及其参数
利用TrueGrid软件建模和划分网格，模型均选用8节点6面体单元。为节约计算时间，考虑到模型的对称性，模型简化为1/4模型。钨球尺寸及靶板厚度与实验状态一致，靶板采用圆形靶，其中防弹衣采用分层建模，每层0.2 mm厚，共45层。靶板半径设为30 mm（大于钨球直径10倍），以减小边界效应对侵彻过程的影响。为兼顾计算的时长与精度，网格采用渐进式，弹着点中心8倍半径区域加密，向外逐渐稀疏。密集区网格尺寸控制在0.10～0.15 mm。稀疏区网格尺寸控制在0.15～1.27 mm，钨球最小网格尺寸为0.01 mm，有限元模型如图5所示。
数值模拟选用LSDYNA-3D软件，单位制设为cm-g-μs，算法采用Lagrange算法。根据模型结构的对称性，在模型对称面添加对称边界条件，在靶板边缘添加无反射边界条件。弹靶之间的接触定义为面面侵蚀接触，防弹衣纤维层与层之间的接触定义为固连失效接触，防弹衣与松木靶之间的接触定义为自动面面接触。
钨球选用考虑应变率效应的弹塑性材料模型（MAT_PLASTIC_KINEMATIC），其材料模型参数[16]见表4，其中ρ为密度，E为弹性模量，μ为泊松比，σy为屈服强度，Et为切线模量，β为硬化参数，R1与R2为应变率参数，εf为失效应变。
表 4 钨球材料模型参数Table 4. Material model parameters of tungsten sphereρ/（g·cm−3） E/GPa μ σy/MPa Et/MPa β R1 R2 εf 18.1 367 0.303 1506 792 1 3.9 6 1.2 防弹衣属于层合结构，选用基于经典层合理论和Chang-Chang失效准则[17-18]的复合材料损伤模型（MAT_COMPOSITE_DAMAGE），具体材料模型参数[19]见表5，其中E1、E2、E3为各向弹性模量，μ21、μ32、μ31为各向泊松比，G12、G23、G31为各向剪切模量，κf为损坏材料体积模量，Gs为面内剪切强度，Tx、Ty分别为纵向和横向拉伸强度，Cy为横向压缩强度，α为非线性剪切应力修正系数，Tn为法向拉伸强度，Gyz和Gzx为横向剪切强度。
表 5 凯夫拉材料模型参数Table 5. Material model parameters of kevlarρ/（g·cm−3） E1/GPa E2/GPa E3/GPa μ21 μ31 μ32 1.35 21 21 4.6 0.31 0.14 0.14 G12/GPa G23/GPa G31/GPa κf/GPa Gs/GPa Tx/GPa Ty/GPa 1.2 1.2 1.2 2 0.35 1.0 1.0 Cy/GPa α Tn/GPa Gyz/GPa Gzx/GPa 0.8 0.5 0.55 0.55 0.55 文献[20]指出，红松木在比例极限下可近似看成弹性，当应力超过比例极限后，红松木发生变形和破坏，故红松木可选用MAT_ELASTIC材料模型并添加单元控制失效算法（ADD_EROSION）予以描述，其材料模型参数[21]见表6，其中σf为失效应力。
表 6 红松木材料模型参数Table 6. Material model parameters of pineρ/（g·cm−3） E/GPa μ σf/MPa 0.46 11.68 0.31 294 2.2 数值模型的验证
利用实验数据对仿真模型进行验证，结果见表7和图6。由表7可以看出，剩余速度的仿真值与实验值最大相对误差不超过10%，满足工程误差要求。同时，根据式(1)对计算得到的剩余速度与着靶速度的拟合可得：钨球侵彻防弹衣加红松木复合靶的弹道极限为690.5 m/s，与实验所得的弹道极限相对误差仅为−0.25%。可见该数值模型及其参数可信。
表 7 模拟值与实验值的对比Table 7. Comparison between simulated results and experimental resultsvi/（m·s−1） vr/（m·s−1） 相对误差/% 实验值 计算值 681.2 0 0 0 700.4 74.8 80.9 8.16 711.7 111.6 119.2 6.81 725.5 160.4 168.3 4.93 744.5 194.2 202.8 4.43 748.0 204.9 209.2 2.10 753.4 220.9 218.4 −1.13 775.7 250.6 253.9 1.32 
图 6 计算值与实验值的对比Figure 6. Comparison between simulated results and experimental results2.3 侵彻过程及破坏机理分析
图7给出了不同着靶速度下钨球的速度v（图7(a)）和加速度a变化曲线（图7(b)）。从图7可看出，不同侵彻速度下钨球的速度和加速度变化趋势相近。以vi=748.0 m/s为例分析侵彻过程，其速度和加速度随时间的变化分别如图8(a)和图8(b)所示。在图8(a)和图8(b)中分别提取5个特征点并用相同字母不同的角标表示：其中A0和A1点表示钨球与防弹衣初始接触时刻；B0和B1点表示钨球侵彻防弹衣阻力达到最大时刻；C0和C1点表示钨球刚穿透防弹衣并开始侵彻松木靶时刻；D0和D1点表示钨球开始稳定侵彻松木靶时刻；E0和E1点表示钨球在松木靶内部稳定侵彻结束时刻；F0和F1点表示钨球完全穿透松木靶时刻。根据钨球速度的衰减过程，可大致将钨球侵彻过程分为侵彻防弹衣（A0C0段）和侵彻松木靶（C0F0段）两阶段。
图 7 不同着靶速度下钨球速度与加速度变化曲线Figure 7. Variation curves of velocity and acceleration of tungsten sphere at different impact velocities
图 8 钨球速度与加速度变化曲线（vi=748.0 m/s）Figure 8. Variation curves of velocity and acceleration of tungsten sphere (vi=748.0 m/s)（1）侵彻防弹衣阶段。
A0B0段，在侵彻防弹衣初期，钨球速度大，动能高，当钨球与防弹衣初始接触时，由于防弹衣的阻碍作用，钨球速度急剧下降，并对防弹衣产生压力作用，此时纤维受到压缩和剪切作用；当剪应力超过极限时，纤维发生剪切断裂破坏。同时，钨球冲击防弹衣产生的应力波沿纤维轴向和靶体纵向两个方向传播；在纤维轴向上，应力波通过基体连接传递至周围其他纤维，二维面积上产生基体开裂、纤维拉伸变形和断裂等现象；在靶体纵向上生成压缩波，压缩波沿靶体纵向传递至基体与纤维的分界面时发生反复透射和反射并形成拉伸波，当拉伸波强度超过基体与纤维间或纤维层间结合强度时，基体与纤维脱粘，纤维层产生分层破坏[22]。侵彻10 μs左右，钨球阻力达到最大，此时穿透17层纤维。B0C0段，随着侵彻的进行，由于断裂破坏的纤维不断增多，且未断裂的纤维已受到一定的拉伸作用使其强度降低，钨球侵彻阻力逐渐减小；虽然钨球速度持续下降，但降低幅度有所减缓，直至贯穿防弹衣到达C0点。
（2）侵彻红松靶阶段。
C0D0段，穿透防弹衣后，钨球主要受松木靶阻力作用；由于松木材质较软，屈服强度较低，因此，钨球侵彻松木靶的阻力较防弹衣的弱，表现在速度衰减较A0C0段缓慢，其侵彻阻力继续减小直至侵彻1.1 mm厚到达D0点。D0E0段，钨球侵入松木靶并进入稳定侵彻阶段，此时松木靶在钨球压力作用下发生剪切破坏。由于木靶材质较均匀，钨球速度几乎呈线性降低，即速度衰减率基本不变，在85 μs左右，稳定侵彻到达E0点，此时钨球侵彻松木靶23.3 mm厚，还剩1.7 mm厚未穿透，松木靶背面产生微小裂纹并发生部分剥落现象。E0F0段：钨球进入穿透松木靶阶段，其头部已经露出松木靶；受靶板背面自由边界的影响，钨球侵彻阻力迅速减小，直到完全穿透松木靶达到F0点，松木靶在钨球的作用下，主要发生剪切和冲塞剥落破坏。
结合图4和图7(a)可看出，随着钨球着靶速度的提高，剩余速度增加但增加斜率逐渐减小。这主要与复合靶的吸能有关。在弹道极限附近的速度范围内，随着着靶速度的提高，钨球侵彻能力增强，弹靶作用时间相对缩短，这样拉伸波不能及时沿靶体纵向传播，导致纤维的分层损伤降低，纤维分层吸收钨球动能的能力下降，故纤维拉伸分层吸能占总吸能比例逐渐降低。图9给出了不同着靶速度下松木靶背面的损伤形貌。从图9可知，随着速度的提高，松木靶背面剥落程度降低，冲塞剥落吸能占总吸能比例降低，复合靶能量吸收不稳定，因而剩余速度曲线随着靶速度的提高变得陡峭。当着靶速度进一步提高时，基体开裂、纤维剪切断裂和松木靶剪切破坏成为复合靶主要吸能方式，复合靶吸收能量比较稳定，剩余速度曲线更加平缓。从图7(b)上看，由于松木靶破坏模式的转变，在钨球穿透松木靶阶段，随着着靶速度的提高，钨球负加速度下降更快。
图 9 不同侵彻速度下松木靶背面的损伤形貌Figure 9. Damage morphologies of the back of pine target at different impact velocities图10给出了防弹衣正面与背面的典型von Mises应力变化过程。
图 10 防弹衣正面与背面的典型von Mises应力变化Figure 10. Typical von Mises stress variation of front and back of body armor由图10可知，破片冲击产生的应力波在纤维层面内通过基体连接传向周围其他纤维，纤维轴向上的波阵面形状类似双纽线[10]，着靶处纤维在应力波作用下产生类似“十”字型的损伤，宏观表现为基体开裂，图11的实验结果也验证了这一点。
2.4 质量变化对侵彻的影响
2.4.1 质量变化对弹道极限的影响
表8给出了不同质量（m）钨球侵彻防弹衣加红松木复合靶弹道极限（vbl）的计算结果。图12给出了弹道极限随钨球质量变化的关系曲线。
表 8 不同质量钨球侵彻防弹衣加红松木复合靶弹道极限的模拟结果Table 8. Simulated results of ballistic limit of tungsten spheres with different mass penetrating into pine target covered with body armorm/g vbl/（m·s−1） 0.21 690.5 0.26 647.0 0.31 619.0 0.36 595.0 0.41 575.0 0.46 556.0 从图12可看出，随着钨球质量的增加，钨球侵彻防弹衣加红松木复合靶的弹道极限呈幂函数递减趋势。从减轻破片战斗部质量提高射程或增加破片数量提高杀伤威力等角度考虑，应存在合适的破片质量范围，既能增加破片数量又能提高破片对人体的杀伤。
2.4.2 质量变化对靶板能量吸收的影响
在侵彻靶板过程中，破片部分动能被靶板吸收，靶板的能量吸收可用能量吸收率来表征[23]，定义靶板的能量吸收率η为：
η=EabsEk=0.5m(v2i−v2r)0.5mv2i=1−v2rv2i (2) 式中：Eabs为靶板吸收的能量，即破片损失的动能；Ek为破片初始动能。
表9给出了不同质量钨球侵彻防弹衣加红松木复合靶的仿真结果。图13给出了靶板能量吸收率η在不同着靶速度下的变化曲线。由图13可知，当钨球嵌入靶板时，靶板的能量吸收率为1。若钨球穿透靶板，在同一着靶速度下，靶板的能量吸收率随钨球质量增加而降低。并且随着着靶速度的增加，靶板能量吸收率急剧下降。
表 9 不同质量钨球侵彻防弹衣加红松木复合靶的仿真结果Table 9. Simulated results of tungsten spheres with different mass penetrating into pine target covered with body armor着靶速度/（m·s−1） m=0.21 g m=0.26 g m=0.31 g 剩余速度/（m·s−1） 靶板能量吸收率 剩余速度/（m·s−1） 靶板能量吸收率 剩余速度/（m·s−1） 靶板能量吸收率 550 0 1 0 1 0 1 600 0 1 0 1 0 1 650 0 1 26.2 0.998 156.7 0.942 700 79.2 0.987 188.3 0.928 254.0 0.868 750 212.6 0.920 281.6 0.859 322.5 0.815 800 294.3 0.865 350.2 0.808 390.3 0.762 着靶速度/（m·s−1） m=0.36 g m=0.41 g m=0.46 g 剩余速度/（m·s−1） 靶板能量吸收率 剩余速度/（m·s−1） 靶板能量吸收率 剩余速度/（m·s−1） 靶板能量吸收率 550 0 1 0 1 0 1 600 65.1 0.988 121.0 0.959 172.8 0.917 650 197.6 0.908 232.4 0.872 267.1 0.831 700 287.4 0.831 315.0 0.798 345.0 0.757 750 354.9 0.776 376.7 0.748 400.3 0.715 800 416.8 0.729 437.1 0.701 457.0 0.674 
图 13 靶板能量吸收率-着靶速度曲线Figure 13. Curve of energy absorption efficiency of target-impact velocity3. 量纲分析
穿靶能量Ec，即消耗破片的动能，是反映破片侵彻靶板特性的一个重要参数。为进一步定量研究钨球穿靶能量与着靶速度以及直径之间的关系，本文依据量纲分析建立小钨球侵彻防弹衣加红松木复合靶的穿靶能量计算公式。
通过理论分析，影响钨球侵彻防弹衣加红松木复合靶穿靶能量的主要物理量如表10所示。
表 10 影响穿靶能量的主要的物理量Table 10. Main physical quantities affecting the energy of penetrating into target材料名称 物理量 量符号 量纲式 钨球 着靶速度 vi LT−1 密度 ρp ML−3 直径 Dp L 弹性模量 Ep L−1MT−2 屈服强度 σsp L−1MT−2 特征应变 εp 1 声速 cp LT−1 防弹衣 密度 ρf ML−3 厚度 hf L 弹性模量 Ef L−1MT−2 抗压强度 σsf L−1MT−2 抗剪强度 στf L−1MT−2 抗拉强度 σff L−1MT−2 特征应变 εf 1 声速 cf LT−1 红松木 密度 ρs ML−3 厚度 hs L 弹性模量 Es L−1MT−2 失效应力 σss L−1MT−2 特征应变 εs 1 声速 cs LT−1 可见，穿靶能量Ec是关于上述物理量的函数：
Ec=f(vi,ρp,Dp,Ep,σsp,εp,cp,ρf,hf,Ef,σsf,στf,σff,εf,cf,ρs,hs,Es,σss,εs,cs) (3) 选取ρp、Dp和σsp为量纲独立变量，根据量纲齐次原则，其他导出量可写成以下无量纲形式：
Π=EcD3pσsp,Π1=vi√σsp/ρp,Π2=Epσsp,Π3=εp,Π4=cp√σsp/ρp,Π5=ρfρp, Π6=hfDp,Π7=Efσsp,Π8=σsfσsp,Π9=στfσsp,Π10=σffσsp,Π11=εf,Π12=cf√σsp/ρp, Π13=ρsρp,Π14=hsDp,Π15=Esσsp,Π16=σssσsp,Π17=εs,Π18=cs√σsp/ρp 根据Π定理，式(3)可写为：
Π=EcD3pσsp=f(Π1,Π2,Π3,⋯,Π18) (4) 在弹靶材料不变的情况下，除了Π、Π1、Π6和Π14外，其余导出量均为常数，式(4)可简化为：
Π=EcD3pσsp=f(vi√σsp/ρp,hfDp,hsDp)=C0(vi√σsp/ρp)α(hfDp)β(hsDp)γ (5) 式中：C0、α、β、γ均为待定常数。
将式（5）展开：
Ec=C0σsp1−0.5αρp0.5αhfβhsγviα(1Dp)β+γ−3 (6) 由于本实验中的防弹衣和红松木材料及厚度不变，钨球材料不变，令：
δ=β+γ−3 C1=C0σsp1−0.5αρp0.5αhfβhsγ 则式（6）可写为：
Ec=C1viα(1Dp)δ (7) 式中：C1、α、δ为待定常数。
由式（7）可知，在防弹衣和红松木材料及其厚度不变以及钨球材料不变的情况下，穿靶能量Ec只与着靶速度vi和钨球直径Dp有关。
为方便求解，将式（7）两端同时取对数：
lnEc=αlnvi+δln1Dp+lnC1 (8) 令
y=lnEc ，x1=lnvi ，x2=ln1Dp ，k=lnC1 ，则式（8）为：y=αx1+δx2+k (9) 利用表9中的数据对式（9）进行二元一次线性回归拟合，可得：
α=0.94,δ=−2.04,k=9.77 将α、δ和k回代至式（7），有：
Ec=e9.77v0.94iDp2.04 (10) 为验证式（10）的有效性，对0.17 g、直径2.6 mm和0.44 g、直径3.6 mm两种小钨球侵彻防弹衣加红松木复合靶进行了侵彻实验。表11对比了通过实验和式（10）计算得到的穿靶能量。由表11可知，利用式（10）计算的穿靶能量与实验所得的穿靶能量最大相对误差不超过10%，满足工程误差要求。无需剩余速度，式（10）仅用着靶速度和破片直径就能预测不同质量钨球在不同着靶速度下的穿靶能量（适用范围为m＜0.46 g）。
表 11 不同方法计算的穿靶能量的对比Table 11. Comparison of energy of penetrating into target calculated by different methods钨球规格 vi/（m·s−1） vr/（m·s−1） Ec/J 相对误差/% 实验结果 量纲分析结果 质量0.17 g
直径2.6 mm753.6 124.6 46.95 47.22 0.58 771.3 187.7 47.57 48.26 1.45 786.0 230.1 48.01 49.12 2.31 809.1 296.8 48.16 50.48 4.82 834.9 357.9 48.36 51.99 7.51 869.6 410.8 49.93 54.02 8.19 质量0.44 g
直径3.6 mm578.4 78.2 72.25 71.52 −1.01 597.7 141.8 74.17 73.76 −0.55 610.9 187.2 74.39 75.29 1.21 639.5 248.0 76.44 78.60 2.83 660.3 260.9 80.94 81.00 0.07 698.3 335.1 82.57 85.37 3.39 值得注意的是，当着靶速度
vi 为弹道极限vbl 时，穿靶能量为极限穿靶能量：Ec=0.5mv2bl (11) 将式（11）代入式（7）可推导出计算弹道极限的另一个公式：
vbl=[2C1m(1Dp)δ]12−α=[2e9.77mDp2.04]11.06 (12) 根据试验数据，利用式(12)计算弹道极限，并与利用式(1) R-I公式计算得到的弹道极限进行对比，结果如表12所示。由表12可知，与式(1)公式相比，式(12)仅利用破片的直径和质量就能计算弹道极限，且相对误差不超过5%，满足工程应用要求。因此，式(12)可作为计算不同质量钨球侵彻防弹衣加红松木复合靶弹道极限的一个经验公式（适用范围为m＜0.46 g）。
表 12 式（1）和式（12）计算的弹道极限的对比Table 12. Comparison of ballistic limit calculated by formula (1) and formula (12)钨球规格 vbl/（m·s−1） 相对误差/% 式(1) 式(12) 质量0.17 g/直径2.6 mm 742.3 738.1 −0.57 质量0.21 g/直径2.8 mm 692.9 697.4 0.65 质量0.44 g/直径3.6 mm 570.1 562.9 −1.26 4. 结 论
采用实验和数值模拟方法研究了小钨球对三级软体防弹衣加25 mm红松靶的侵彻过程及破坏机理；探讨了钨球质量变化对弹道极限及靶板能量吸收的影响；利用量纲分析研究了钨球穿靶能量与其初速及直径的关系，建立了钨球的穿靶能量公式及弹道极限公式。得出研究结论如下：
(1) 0.17 g，直径2.6 mm，0.21 g、直径2.8 mm以及0.44 g、直径3.6 mm的小钨球侵彻防弹衣加红松木复合靶的弹道极限分别为742.3、692.9和570.1 m/s；
(2) 在钨球侵彻防弹衣加红松木复合靶的过程中，防弹衣的破坏形式主要表现为基体开裂、纤维断裂和拉伸分层，纤维层面内出现类似“十”字型的损伤，而松木靶则以剪切和冲塞剥落破坏为主；
(3) 随着钨球质量的增加，钨球侵彻防弹衣加红松木复合靶的弹道极限呈幂函数减小趋势，靶板的能量吸收率逐渐降低；
(4) 利用量纲分析法建立了钨球侵彻防弹衣加红松木复合靶的穿靶能量计算公式及弹道极限计算公式，可分别用于计算不同侵彻速度下的穿靶能量和不同质量钨球的弹道极限。
-
图 3 不同煤粉质量浓度下pmax随甲烷体积分数的变化
Figure 3. Variations of pmax with methane volume fraction at different coal dust concentrations
图 4 不同甲烷体积分数下，pmax和(dp/dt)max随煤粉质量浓度的变化
Figure 4. Variations of pmax and (dp/dt)max with the mass concentration of coal dust at different methane volume fractions
图 5 瓦斯煤尘爆炸(dp/dt)max及其来临时刻（煤粉质量浓度60 g/m3）
Figure 5. (dp/dt)max of methane coal dust explosion and its moment of arrival (coal dust mass concentration: 60 g/m3)
图 6 不同粒径煤粉的Kst随甲烷体积分数的变化（煤粉浓度100 g/m3）
Figure 6. Variation of Kst with methane volume fraction for different particle sizes of coal dust (coal dust concentration: 100 g/m3)
表 1 煤粉的工业分析
Table 1. Industrial analysis of coal dust
煤样 ω（挥发分）/% ω（水分）/% ω（灰分）/% ω（固定碳）/% 褐煤 38.02 6.68 4.98 54.75 无烟煤 7.97 2.58 10.59 79.91 
下载: 导出CSV
表 2 正交实验各因素及水平参数
Table 2. Factors and level parameters of orthogonal experiment
水平 煤粉质量浓度/(g·m−3) 甲烷体积分数/% 煤粉粒径/μm 煤粉种类 1 30 6 39 褐煤 2 60 9 123 无烟煤 3 100 11 283 − 4 200 − − −
下载: 导出CSV
表 3 正交实验的极差分析结果
Table 3. Results of extreme difference analysis of orthogonal experiments
指标 因变量 煤粉质量浓度/(g·m−3) 甲烷体积分数/% 煤粉粒径/μm 煤粉种类 K1 pmax 6.119 7.752 8.048 12.256 (dp/dt)max 238.406 181.905 283.525 446.191 K2 pmax 6.132 8.489 8.153 12.041 (dp/dt)max 228.572 386.572 301.977 437.548 K3 pmax 6.137 8.056 8.096 — (dp/dt)max 220.595 315.262 298.237 — K4 pmax 5.909 — — — (dp/dt)max 196.166 — — — k1 pmax 0.680 0.646 0.671 0.681 (dp/dt)max 26.489 15.159 23.627 24.788 k2 pmax 0.681 0.707 0.679 0.669 (dp/dt)max 25.397 32.214 25.165 24.308 k3 pmax 0.682 0.671 0.675 — (dp/dt)max 24.511 26.272 24.853 — k4 pmax 0.657 — — — (dp/dt)max 21.796 — — — 极差R pmax 0.025 0.061 0.009 0.012 (dp/dt)max 4.693 17.055 1.538 0.480 最佳水平 pmax 100 9 138 褐煤 (dp/dt)max 30 9 138 褐煤 因素主次 pmax 甲烷体积分数、煤粉质量浓度、煤粉种类、煤粉粒径 (dp/dt)max 甲烷体积分数、煤粉质量浓度、煤粉种类、煤粉粒径
下载: 导出CSV
表 4 正交实验的方差分析结果
Table 4. Results of the variance analysis of orthogonal experiments
试验因素 因变量 平方和 均方 F 影响程度 煤粉质量浓度 pmax 0.005 0.002 1.307 — (dp/dt)max 108.575 36.192 3.597 显著 甲烷体积分数 pmax 0.023 0.012 10.037 高度显著 (dp/dt)max 1798.827 899.414 89.378 高度显著 煤粉粒径 pmax 0.001 0.001 0.201 — (dp/dt)max 15.859 7.930 0.788 — 煤粉种类 pmax 0.001 0.001 1.094 — (dp/dt)max 2.075 2.075 0.206 — 随机误差E pmax 0.031 0.001 — — (dp/dt)max 271.694 10.063 — —
下载: 导出CSV
-
[1] 张巍, 张帆, 张军, 等. 与新能源耦合发展 推动现代煤化工绿色低碳转型的思考与建议 [J]. 中国煤炭, 2021, 47(11): 56–60. DOI: 10.19880/j.cnki.ccm.2021.11.009.ZHANG W, ZHANG F, ZHANG J, et al. Coupling development with new energy and thinking and suggestions on promoting the green and low-carbon transformation of modern coal chemical industry [J]. China Coal, 2021, 47(11): 56–60. DOI: 10.19880/j.cnki.ccm.2021.11.009. [2] 毛健雄. 燃煤耦合生物质发电 [J]. 分布式能源, 2017, 2(5): 47–54. DOI: 10.16513/j.cnki.10-1427/tk.2017.05.008.MAO J X. Co-firing biomass with coal for power generation [J]. Distributed Energy, 2017, 2(5): 47–54. DOI: 10.16513/j.cnki.10-1427/tk.2017.05.008. [3] 马仑, 方庆艳, 张成, 等. 深度空气分级下煤粉耦合氨燃烧及NO生成特性 [J]. 洁净煤技术, 2022, 28(3): 201–213. DOI: 10.13226/j.issn.1006-6772.CC21101401.MA L, FANG Q Y, ZHANG C, et al. Combustion and NO formation characteristics of pulverized coal co-firing with ammonia in a deep-air staging condition [J]. Clean Coal Technology, 2022, 28(3): 201–213. DOI: 10.13226/j.issn.1006-6772.CC21101401. [4] 王德明, 邵振鲁, 朱云飞. 煤矿热动力重大灾害中的几个科学问题 [J]. 煤炭学报, 2021, 46(1): 57–64. DOI: 10.13225/j.cnki.jccs.YG20.1798.WANG D M, SHAO Z L, ZHU Y F. Several scientific issues on major thermodynamic disasters in coal mines [J]. Journal of China Coal Society, 2021, 46(1): 57–64. DOI: 10.13225/j.cnki.jccs.YG20.1798. [5] CASHDOLLAR K L. Overview of dust explosibility characteristics [J]. Journal of Loss Prevention in the Process Industries, 2000, 13(3): 183–199. DOI: 10.1016/S0950-4230(99)00039-X. [6] GOING J E, CHATRATHI K, CASHDOLLAR K L. Flammability limit measurements for dusts in 20-L and 1-m3 vessels [J]. Journal of Loss Prevention in the Process Industries, 2000, 13(3): 209–219. DOI: 10.1016/S0950-4230(99)00043-1. [7] KUNDU S K, ZANGANEH J, ESCHEBACH D, et al. Explosion severity of methane–coal dust hybrid mixtures in a ducted spherical vessel [J]. Powder Technology, 2018, 323: 95–102. DOI: 10.1016/j.powtec.2017.09.041. [8] MISHRA D P, AZAM S. Experimental investigation on effects of particle size, dust concentration and dust-dispersion-air pressure on minimum ignition temperature and combustion process of coal dust clouds in a G-G furnace [J]. Fuel, 2018, 227: 424–433. DOI: 10.1016/j.fuel.2018.04.122. [9] CLONEY C T, RIPLEY R C, PEGG M J, et al. Lower flammability limits of hybrid mixtures containing 10 micron coal dust particles and methane gas [J]. Process Safety and Environmental Protection, 2018, 120: 215–226. DOI: 10.1016/j.psep.2018.09.004. [10] CLONEY C T, RIPLEY R C, PEGG M J, et al. Role of particle diameter in the lower flammability limits of hybrid mixtures containing coal dust and methane gas [J]. Journal of Loss Prevention in the Process Industries, 2019, 61: 206–212. DOI: 10.1016/j.jlp.2019.06.010. [11] 李庆钊, 翟成, 吴海进, 等. 基于20 L球形爆炸装置的煤尘爆炸特性研究 [J]. 煤炭学报, 2011, 36(S1): 119–124. DOI: 10.13225/j.cnki.jccs.2011.s1.031.LI Q Z, ZHAI C, WU H J, et al. Investigation on coal dust explosion characteristics using 20 L explosion sphere vessels [J]. Journal of China Coal Society, 2011, 36(S1): 119–124. DOI: 10.13225/j.cnki.jccs.2011.s1.031. [12] LI Y, XU H L, WANG X S. Experimental study on the influence of initial pressure on explosion of methane-coal dust mixtures [J]. Procedia Engineering, 2013, 62: 980–984. DOI: 10.1016/j.proeng.2013.08.151. [13] SONG S X, CHENG Y F, MENG X R, et al. Hybrid CH4/coal dust explosions in a 20-L spherical vessel [J]. Process Safety and Environmental Protection, 2019, 122: 281–287. DOI: 10.1016/j.psep.2018.12.023. [14] LIU Q M, BAI C H, LI X D, et al. Coal dust/air explosions in a large-scale tube [J]. Fuel, 2010, 89(2): 329–335. DOI: 10.1016/j.fuel.2009.07.010. [15] 司荣军, 李润之, 苏岱峰. 煤尘云质量浓度对瓦斯爆炸压力影响的试验研究 [J]. 安全与环境学报, 2018, 18(5): 1796–1798. DOI: 10.13637/j.issn.1009-6094.2018.05.025.SI R J, LI R Z, SU D F. Investigation of the influence of the coal dust cloud on the gas explosion pressure [J]. Journal of Safety and Environment, 2018, 18(5): 1796–1798. DOI: 10.13637/j.issn.1009-6094.2018.05.025. [16] 李润之. 瓦斯煤尘共存条件下的煤尘云爆炸下限 [J]. 爆炸与冲击, 2018, 38(4): 913–917. DOI: 10.11883/bzycj-2016-0331.LI R Z. Minimum explosive concentration of coal dust cloud in the coexistence of gas and coal dust [J]. Explosion and Shock Waves, 2018, 38(4): 913–917. DOI: 10.11883/bzycj-2016-0331. [17] 景国勋, 张胜旗, 段新伟, 等. 竖直管道内煤尘浓度对瓦斯爆炸特性影响研究 [J]. 中国安全科学学报, 2020, 30(3): 15–20. DOI: 10.16265/j.cnki.issn1003-3033.2020.03.003.JING G X, ZHANG S Q, DUAN X W, et al. Impact of coal dust concentration on gas explosion in vertical pipelines [J]. China Safety Science Journal, 2020, 30(3): 15–20. DOI: 10.16265/j.cnki.issn1003-3033.2020.03.003. [18] 景国勋, 彭乐, 班涛, 等. 甲烷煤尘耦合爆炸传播特性及伤害研究 [J]. 中国安全科学学报, 2022, 32(1): 72–78. DOI: 10.16265/j.cnki.issn1003-3033.2022.01.010.JING G X, PENG L, BAN T, et al. Research on pressure propagation and injury of methane and coal dust coupled explosion [J]. China Safety Science Journal, 2022, 32(1): 72–78. DOI: 10.16265/j.cnki.issn1003-3033.2022.01.010. [19] CAO W G, QIN Q F, CAO W, et al. Experimental and numerical studies on the explosion severities of coal dust/air mixtures in a 20-L spherical vessel [J]. Powder Technology, 2017, 310: 17–23. DOI: 10.1016/j.powtec.2017.01.019. [20] 裴蓓, 朱知印, 余明高, 等. 瓦斯/煤尘爆炸初期复合火焰加速及灾害强化机制分析 [J]. 工程热物理学报, 2021, 42(7): 1879–1886.PEI B, ZHU Z Y, YU M G, et al. Analysis on the acceleration of composite flame and the strengthening mechanism of disaster in the initial stage of gas/coal dust explosion [J]. Journal of Engineering Thermophysics, 2021, 42(7): 1879–1886. [21] 陈晓坤, 张自军, 王秋红, 等. 20 L近球形容器中微米级铝粉的爆炸特性 [J]. 爆炸与冲击, 2018, 38(5): 1130–1136. DOI: 10.11883/bzycj-2017-0101.CHEN X K, ZHANG Z J, WANG Q H, et al. Explosion characteristics of micro-sized aluminum dust in 20L spherical vessel [J]. Explosion and Shock Waves, 2018, 38(5): 1130–1136. DOI: 10.11883/bzycj-2017-0101. [22] 覃小玲, 李晓泉. NH4H2PO4对蔗糖粉尘爆炸的抑制作用试验研究 [J]. 中国安全科学学报, 2020, 30(4): 41–46. DOI: 10.16265/j.cnki.issn1003-3033.2020.04.007.QIN X L, LI X Q. Experimental research on suppression of NH4H2PO4 on sucrose dust explosion [J]. China Safety Science Journal, 2020, 30(4): 41–46. DOI: 10.16265/j.cnki.issn1003-3033.2020.04.007. [23] WU X L, XU S, PANG A M, et al. Hazard evaluation of ignition sensitivity and explosion severity for three typical MH2 (M= Mg, Ti, Zr) of energetic materials [J]. Defence Technology, 2021, 17(4): 1262–1268. DOI: 10.1016/j.dt.2020.06.011. [24] PICO P, RATKOVICH N, MUÑOZ F, et al. Analysis of the explosion behaviour of wheat starch/pyrolysis gases hybrid mixtures through experimentation and CFD-DPM simulations [J]. Powder Technology, 2020, 374: 330–347. DOI: 10.1016/j.powtec.2020.07.016. [25] WANG X W, WANG Z R, NI L, et al. Explosion characteristics of aluminum powder in different mixed gas environments [J]. Powder Technology, 2020, 369: 53–71. DOI: 10.1016/j.powtec.2020.04.056. [26] 徐伟巍, 王家祎, 覃欣欣, 等. 酒精蒸气-烟草粉尘耦合体系燃爆猛度特性研究 [J]. 工业安全与环保, 2021, 47(5): 9–12. DOI: 10.3969/j.issn.1001-425X.2021.05.003.XU W W, WANG J Y, QIN X X, et al. Study on the explosion fierce characteristics of alcohol vapor-tobacco dust coupling [J]. Industrial Safety and Environmental Protection, 2021, 47(5): 9–12. DOI: 10.3969/j.issn.1001-425X.2021.05.003. [27] 孙超伦, 张一民, 裴蓓, 等. 惰气/赤泥两相抑爆剂抑制瓦斯爆炸试验研究 [J]. 中国安全科学学报, 2020, 30(10): 112–118. DOI: 10.16265/j.cnki.issn1003-3033.2020.10.016.SUN C L, ZHANG Y M, PEI B, et al. Experimental study on suppression effects of inert gas/red mud two-phase inhibitors on gas explosion [J]. China Safety Science Journal, 2020, 30(10): 112–118. DOI: 10.16265/j.cnki.issn1003-3033.2020.10.016. [28] ABBAS Z, ZINKE R, GABEL D, et al. Theoretical evaluation of lower explosion limit of hybrid mixtures [J]. Journal of Loss Prevention in the Process Industries, 2019, 60: 296–302. DOI: 10.1016/j.jlp.2019.05.014. [29] ABBAS Z, GABEL D, KRIETSCH A, et al. Quasi-static dispersion of dusts for the determination of lower explosion limits of hybrid mixtures [J]. Journal of Loss Prevention in the Process Industries, 2022, 74: 104640. DOI: 10.1016/j.jlp.2021.104640. [30] JIANG J J, LIU Y, MANNAN M S. A correlation of the lower flammability limit for hybrid mixtures [J]. Journal of Loss Prevention in the Process Industries, 2014, 32: 120–126. DOI: 10.1016/j.jlp.2014.07.014. [31] ADDAI E K, CLOUTHIER M, AMYOTTE P, et al. Experimental investigation of limiting oxygen concentration of hybrid mixtures [J]. Journal of Loss Prevention in the Process Industries, 2019, 57: 120–130. DOI: 10.1016/j.jlp.2018.11.016. [32] ADDAI E K, GABEL D, KAMAL M, et al. Minimum ignition energy of hybrid mixtures of combustible dusts and gases [J]. Process Safety and Environmental Protection, 2016, 102: 503–512. DOI: 10.1016/j.psep.2016.05.005. [33] 杨虎, 钟波, 刘琼荪. 应用数理统计 [M]. 北京: 清华大学出版社, 2006: 123-158.YANG H, ZHONG B, LIU Q S. Application mathematical statistic study [M]. Beijing: Tsinghua University Press, 2006: 123-158. [34] 王亚军, 徐秀艳, 秦宪礼. 煤粉对泡沫金属抑制爆炸火焰波性能的影响规律 [J]. 煤炭学报, 2017, 42(11): 2885–2891. DOI: 10.13225/j.cnki.jccs.2017.0476.WANG Y J, XU X Y, QIN X L. Study on the influence regulation after adding coal dust for inhibition of flame wave of gas explosion by foam metal [J]. Journal of China Coal Society, 2017, 42(11): 2885–2891. DOI: 10.13225/j.cnki.jccs.2017.0476. [35] GARCIA-AGREDA A, DI BENEDETTO A, RUSSO P, et al. Dust/gas mixtures explosion regimes [J]. Powder Technology, 2011, 205(1): 81–86. DOI: 10.1016/j.powtec.2010.08.069. [36] 周永浩, 甘波, 姜海鹏, 等. 甲烷/煤尘复合爆炸火焰的传播特性 [J]. 爆炸与冲击, 2022, 42(1): 015402. DOI: 10.11883/bzycj-2021-0064.ZHOU Y H, GAN B, JIANG H P, et al. Investigations on the flame propagation characteristics in methane and coal dust hybrid explosions [J]. Explosion and Shock Waves, 2022, 42(1): 015402. DOI: 10.11883/bzycj-2021-0064. [37] NISHIMURA I, MOGI T, DOBASHI R. Simple method for predicting pressure behavior during gas explosions in confined spaces considering flame instabilities [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(2): 351–354. DOI: 10.1016/j.jlp.2011.08.009. [38] 汤其建, 秦汝祥, 戴广龙. 索特平均直径对煤粉及其在瓦斯气氛下爆炸特性的影响 [J]. 煤炭学报, 2021, 46(2): 489–497. DOI: 10.13225/j.cnki.jccs.2019.1539.TANG Q J, QIN R X, DAI G L. Effect of Sauter mean diameter of coal dust on its explosibility with and without methane gas [J]. Journal of China Coal Society, 2021, 46(2): 489–497. DOI: 10.13225/j.cnki.jccs.2019.1539. 期刊类型引用(4)
1. 付世杰, 邵伟平, 郝永平. 基于自注意与图神经网络的弹体侵彻靶板材质检测方法研究. 弹箭与制导学报. 2025(03) 
百度学术2. 张磊，石少卿，王昊康，田镇华，王文康，廖瑜. 玻璃纤维掩蔽部顶盖构件抗爆炸性能分析. 合成纤维. 2024(07): 57-61+78 . 百度学术3. 唱忠良，张志俊，吴晨晨，邹广平，吴松阳，闫安石，李寅辰. 平纹机织可展收层合壳弹道冲击特性研究. 机械强度. 2024(04): 853-858 . 百度学术4. 董方栋，杨耿，王舒，王志军，冯志威. 刚性破片侵彻层合板运动模型. 兵工学报. 2022(09): 2350-2359 . 百度学术其他类型引用(0)
-
下载:
百度学术
计量
- 文章访问数: 565
- HTML全文浏览量: 153
- PDF下载量: 57
- 被引次数: 4