Influence of Detonation Modes on Energy Release Characteristics of a Charge with a Non-Circular Cross-Sectional Structure
-
摘要: 为了研究不同起爆方式下非圆截面装药结构的释能规律,采用AUTODYN软件开展了非圆截面装药结构在不同起爆方式下的释能特性数值模拟,分析了起爆方式对爆轰波形演变、破片质量、破片初速的影响。结果表明:由于装药结构的特殊性,采用端部单点起爆时装药能量分布不均匀,部分区域产生大量的无效小质量破片,且不同位置处的破片初速波动较大;采用端部两点和端部三点起爆时,能够对爆轰能量起到匀化效果,减少无效破片数量,提升破片初速的一致性。由此证明通过调整起爆方式可以对非圆截面装药结构的能量输出结构进行有效调控,对其周向能量场起到匀化效果。Abstract: In order to study the energy release law of a charge with a non-circular cross-sectional structure under the influence of initiation modes, the energy release characteristics of the charge with a non-circular cross-sectional structure under different initiation modes were calculated by AUTODYN simulation software, then effects of initiation modes on the evolution of detonation waveforms, and masses and initial velocities of fragments were analyzed. The results show that, due to particularity of the charge’s structure, energy distribution of the charge is uneven when using single-endpoint initiation, which results in a large number of invalid fragments with small masses in some areas, and a large fluctuation of initial velocities of fragments at different positions. When the two-endpoint and three-endpoint initiation modes are adopted, the detonation energy is homogenized, the number of invalid fragments is reduced, and the consistency of fragments’ initial velocities is improved. It proves that the energy output structure of a charge with a non-circular cross-sectional structure could be effectively regulated by adjusting the initiation mode, and the circumferential energy field could be homogenized.
-
近年来,随着新型战斗部技术的不断发展,出现了一类非圆截面结构战斗部,该战斗部的外形及内部装药结构并非传统的轴对称形状,因而其能量释放特性也不同于传统轴对称战斗部[1-2]。传统战斗部的毁伤元在环向上可近似为均匀分布,而非圆截面战斗部由于装药结构的特殊性,能量输出结构与传统战斗部存在明显差异。因此,对于非圆截面结构战斗部,如何有效控制其爆炸过程,进而优化其能量释放结构,减少装药能量的损失,实现装药爆炸能量的充分释放及高效利用,是目前急需解决的问题。
改变起爆方式是一种优化战斗部能量输出结构的常见方法,在聚能射流[3–5]、爆炸成形弹丸(EFP)[6–7]、聚能杆式侵彻体[8]、破片杀伤[9–12]、复合战斗部[13]、水下战斗部[14]、FAE战斗部[15]等领域具有广泛应用,但有关起爆方式对非圆截面装药结构释能特性的影响鲜有报道。因此,本研究针对一种典型的非圆截面结构战斗部,采用AUTODYN数值模拟方法,分析不同起爆方式对装药结构能量输出特性的影响,以期为提升非圆截面装药的能量利用率,进而提升该类战斗部的综合毁伤威力提供参考。
1. 计算模型
1.1 几何模型
非圆截面装药结构战斗部截面如图1所示。战斗部由壳体和炸药组成,其中壳体由两段圆弧拼接而成:上段圆弧内径R1为52 mm,圆心角为180°;下段圆弧内径R2为67 mm,圆心角为102°;壳体壁厚3 mm,战斗部高度60 mm。
起爆方式分别为端部单点、端部两点、端部三点起爆。起爆点位置的选择基于以下原则:端部单点起爆时,起爆点位于装药截面的几何中心;端部两点起爆时,起爆点关于装药截面对称轴对称分布,以左侧起爆点为例,起爆点的选择基于在截面内到左侧壳体及截面对称轴所有点的距离的方差最小原则确定;端部三点起爆时,1个起爆点位于装药截面的几何中心,另外2个起爆点关于装药截面对称轴对称分布,起爆位置在两点起爆位置的基础上向对称轴两侧各偏移5 mm。具体起爆点位置如图2所示。
1.2 有限元模型
采用非线性动力分析软件AUTODYN和基于Mott破片分布理论的Stochastic随机模型进行计算。为提高计算效率,建立三维1/2模型。壳体采用Lagrange算法,炸药及空气采用Euler算法,空气模型设置流出边界,整体采用流固耦合。有限元模型如图3所示,网格划分如表1所示。
表 1 不同部件的网格划分Table 1. Meshing of different partsPart Grid size/mm Elements number Nodes number Grid quality Shell 3 9 702 12 320 ≥0.95 Charge 3 66 943 87 695 ≥0.90 Air 1 1 450 000 ≥0.95 1.3 材料模型
空气材料采用IDEAL AIR 状态方程,其表达式为
p=(γ−1)ρe−pshift (1) 式中:p为压力,
γ 为绝热指数,ρ 为密度,e为内能,pshift为一个较小的初始压力。具体参数如表2所示。表 2 空气参数Table 2. Parameters of the airγ ρ/(g·cm−3) e/kJ pshift/GPa 1.4 0.001 225 206.8 0 炸药材料为Comp. B,采用JWL状态方程,表达式为
p=A(1−ωR1V)e−R1V+B(1−ωR2V)e−R2V+ωEV (2) 式中:A、B、R1、R2、
ω 为JWL状态方程参数,V为相对体积,E为单位体积内能。具体参数如表3所示。其中,DCJ为CJ爆轰速度,pCJ为CJ爆轰压力,E0为单位体积内能初始值。表 3 炸药参数Table 3. Parameters of the explosiveDCJ/(km·s−1) pCJ/GPa A/GPa B/GPa R1 R2 ω E0/(kJ·cm−3) 7.98 29.5 524.23 7.678 4.2 1.1 0.34 8.5 壳体材料为Steel4340,采用Johnson-Cook强度模型和Linear状态方程描述,失效模型采用基于Mott破片分布理论的Stochastic随机破坏模型[16],失效准则为Principal Stress。强度模型表达式为
Y=(A0+B0εnp)(1+Cln˙εp˙ε0)[1−(T−T0Tm−T0)m] (3) 式中:Y为剪切强度,A0、B0、C、n和m为材料常数,T0为参考温度,Tm为熔点,
εp 和˙εp 分别为等效塑性应变和等效塑性应变率,˙ε0 为参考应变率。壳体材料参数如表4所示。表 4 壳体参数Table 4. Parameters of the shellA0/MPa B0/MPa C n m T0/K Tm/K ˙ε0/s−1 792 510 0.014 0.26 1.03 294 1 793 1 上述所有材料模型参数均取自AUTODYN软件的材料数据库。
2. 结果与分析
2.1 模型验证
为了对模型进行验证,开展了端部单点起爆非圆截面装药结构的爆轰波形观测试验,试验所用样品的技术状态与端部单点起爆模型的技术状态一致。试验仪器采用SJZ-15型转镜式高速分幅相机,摄影频率设置为106 s−1,试验布局见图4。
数值模拟与试验结果的对比如图5所示,可以看出,数值模拟得到的爆轰波形(图5中虚线所示)演变过程与试验结果基本吻合,证明该模型可用于开展非圆截面装药的释能特性研究。
2.2 波形演变
不同起爆方式下的爆轰波形演变如图6所示。当采用端部单点起爆时,炸药起爆以后,爆轰波呈球形传播,并于5 μs左右传播至壳体表面,驱动该处壳体运动、破裂,直至产生破片,随后爆轰波继续向左右两侧传播,并逐渐传播至上下两段圆弧形成的夹角处。采用端部两点起爆时,起爆点处的爆轰波首先以球形传播,约2 μs时,两列波阵面相遇,产生汇聚叠加;在随后的传播过程中,波形不断匀化,并与壳体外形更加贴合,爆轰能量的空间分布也更加均匀,从而减少了爆轰波作用于壳体同截面上不同位置处的时间差,保证壳体破裂的同步性。采用端部三点起爆时的爆轰波形演变与端部两点起爆类似,也包含爆轰波的传播-汇聚-匀化过程,由于增加了起爆点,爆轰波的匀化效果相比端部两点起爆更加充分,波形与壳体也更加贴合,但随着爆轰波的不断传播、演变,其波形逐渐与端部两点起爆趋于一致,至5 μs左右,即爆轰波作用于壳体表面,开始对其做功加速时,二者的波形已基本一致。
2.3 壳体破裂过程
不同起爆方式下的壳体膨胀、破裂、形成破片过程如图7所示。20 μs左右时,壳体开始破裂,随着时间延长,至80 μs时破裂基本完成,形成自然破片。
对80 μs时破片的空间分布状态进行分析,可以看出,轴线方向上,由于受到端部稀疏波的影响,壳体中心部位的破片速度高于两端,反映到图7中,表现为不同速度的破片形成中心部位向外凸出的弧形面结构。径向上,对于端部单点起爆,破片的空间分布形状与截面形状类似,仍然呈椭圆形;端部两点及三点起爆时的破片空间分布形状近乎呈圆形,这主要是由于采用多点起爆时,爆轰波的对心碰撞会在起爆点的对称平面区域形成马赫发射[17],进而在该区域产生比CJ爆压更强的局部超压,使得该方向的抛射压力增大,破片速度增加,从而导致破片空间分布近乎呈圆形。
此外,从图7还可以看出,无论采用哪种起爆方式,在两段圆弧相交位置均存在15°~20°的破片稀疏区,该区域的破片密度很小,这是由本研究中壳体结构的截面特性导致的。
2.4 破片质量及速度
2.4.1 破片质量
不同起爆方式下的破片质量分布如表5所示。采用单点起爆时,破片总数为518,其中大于0.4 g的破片数为264,占破片总数的51.0%;两点起爆时,破片总数为528,大于0.4 g的破片数为330,占破片总数的62.5%;三点起爆时,破片总数为532,大于0.4 g的破片数为346,占破片总数的65.0%。由此可知,起爆方式对破片质量分布具有一定的影响,随着起爆点的增加,形成的有效破片数也逐渐增加。分析认为,采用不同起爆方式时,由于改变了装药的空间能量分布,致使壳体不同位置处的有效装药量产生变化,进而导致有效破片总数存在差别。
表 5 破片质量分布Table 5. Mass distribution of fragmentsMass range/g Fragments number Single-endpoint Two-endpoint Three-endpoint 4.0–8.0 12 10 8 1.0–4.0 86 112 120 0.8–1.0 40 52 62 0.6–0.8 32 52 54 0.4–0.6 94 104 102 <0.4 254 198 186 Total 518 528 532 根据装药的结构特征,将其截面划分为A、B、C、D共4个区域,如图8所示,其中,B、D区域可近似等效为圆截面装药结构的一部分,该区域破片的驱动过程可等效为圆截面装药对破片的驱动,而A、C区域为非圆截面装药结构。计算得出A、B、C、D区域内单位长度内炸药质量m0与壳体质量M的比值(m0/M)分别为1.56、1.36、1.56、1.35,由此可知,作用于A、C区域单位质量破片上的装药量高于B、D区域,并且由于壳体厚度仅有3 mm,因此,A、C区域的破片在炸药作用下更容易碎化,即更容易产生质量较小的破片。而采用多点起爆时,通过2.2节分析可知,会对爆轰波形产生匀化效应,虽然爆轰波对心碰撞会产生马赫反射,造成局部高压,但是由于受到装药结构尺寸限制,马赫反射区域小于图8所示的A、C区域。以两点起爆为例,图9虚线部位所示即为该起爆方式下的马赫反射区。由此可见,采用多点起爆会使作用于壳体不同位置上的炸药能量分布更加均匀,从而避免某区域壳体由于能量过高导致的破片碎化,减少无效破片的数量。统计计算结果:单点起爆时,A、C区域中质量小于0.4 g的破片数为202,占该类破片总数的79.5%,而A、C区域壳体质量占壳体总质量的68.1%;两点和三点起爆时,该类破片数分别为128和130,占该类破片总数的64.6%和69.9%。这也从侧面证明了相对于端部单点起爆方式,通过增加起爆点数量,可以对作用于壳体上的炸药能量进行匀化。
2.4.2 破片速度
以装药结构中部截面上的破片为研究对象,提取不同位置的破片初速沿径向的分布,如图10所示。可以看出,采用端部单点起爆时,两圆弧交界处(图8中的A、C区域)的破片初速波动较大,且相对其他位置处的破片初速较低。根据2.4.1节的分析可知,该区域单位质量破片的装药量高于B、D区域,因此,依据格尼公式,该区域的破片初速应高于B、D区域,但实际结果却正好相反。经分析,这主要是由于依据格尼公式计算破片初速时,忽略了反应区后稀疏波的影响,而本装药结构由于截面的非对称性,爆炸时B、D区域的破片被提前驱动加速,继而稀疏波介入,而此时爆轰波还未传播至A、C区域,因此,稀疏波的介入会影响爆轰波对A、C区域破片的驱动,导致该区域破片初速相比B、D区域较低。由此可知,对于本研究中的非圆截面装药结构,稀疏波对破片初速,特别是A、C区域破片初速的影响不可忽略,若采用格尼公式对该装药结构的破片初速进行预估,会与实际结果存在较大误差。采用端部两点和端部三点起爆时,由于对爆炸能量进行了匀化,缩短了爆轰波传播至截面不同位置处的时间差,使得破片几乎被同时驱动,因此,沿径向不同位置的破片初速较平均。同时,在截面的对称轴附近,由于爆轰波的对心碰撞会产生马赫反射,在该区域形成局部高压,因而该区域的破片初速高于其他区域。通过对图10的数据进行统计,得到不同起爆方式下的破片初速分布,如表6所示。
从表6可以看出,不同起爆方式下的破片速度均值相差不大,但是端部两点和三点起爆时的速度标准差明显低于端部单点起爆,仅为端部单点起爆速度标准差的57.8%和60.9%,与图10反映的内容一致。这也证明采用端部单点起爆时不同位置的破片初速存在较大波动,而采用多点起爆方式能够使破片速度匀化,减小不同位置处的破片速度差。通过表6还可以看出,端部两点与端部三点起爆时的破片速度标准差几乎一致,证明这两种起爆方式下的破片初速波动相差不大。
表 6 破片速度分布Table 6. Velocity distribution of fragmentsDetonation mode Velocity range/(m·s−1) Average velocity/(m·s−1) Velocity standard deviation/(m·s−1) Single-endpoint 1 551–1 758 1 660 64 Two-endpoint 1 572–1 735 1 649 37 Three-endpoint 1 563–1 760 1 657 39 3. 结 论
(1) 随着起爆点数的增加,非圆截面装药结构的爆轰波形逐渐匀化,爆轰波与装药壳体更加贴合,爆轰能量的空间分布更加均匀。
(2) 端部单点起爆时,由于部分区域内装药与壳体质量的比值偏高,导致该区域壳体更容易碎化成质量较小的无效破片;增加起爆点后,通过对装药能量进行匀化,能够有效减少该区域及整个壳体的无效破片数。
(3) 不同起爆方式下的破片初速均值相差不大,但是端部单点起爆时不同位置的破片初速波动较大,其标准差为64 m/s,端部两点和端部三点起爆时的破片初速标准差为端部单点起爆时的57.8%和60.9%,证明增加起爆点后,能够提高不同位置的破片初速的一致性,减小其初速波动。
-
表 1 不同部件的网格划分
Table 1. Meshing of different parts
Part Grid size/mm Elements number Nodes number Grid quality Shell 3 9 702 12 320 ≥0.95 Charge 3 66 943 87 695 ≥0.90 Air 1 1 450 000 ≥0.95 表 2 空气参数
Table 2. Parameters of the air
γ ρ/(g·cm−3) e/kJ pshift/GPa 1.4 0.001 225 206.8 0 表 3 炸药参数
Table 3. Parameters of the explosive
DCJ/(km·s−1) pCJ/GPa A/GPa B/GPa R1 R2 ω E0/(kJ·cm−3) 7.98 29.5 524.23 7.678 4.2 1.1 0.34 8.5 表 4 壳体参数
Table 4. Parameters of the shell
A0/MPa B0/MPa C n m T0/K Tm/K ˙ε0/s−1 792 510 0.014 0.26 1.03 294 1 793 1 表 5 破片质量分布
Table 5. Mass distribution of fragments
Mass range/g Fragments number Single-endpoint Two-endpoint Three-endpoint 4.0–8.0 12 10 8 1.0–4.0 86 112 120 0.8–1.0 40 52 62 0.6–0.8 32 52 54 0.4–0.6 94 104 102 <0.4 254 198 186 Total 518 528 532 表 6 破片速度分布
Table 6. Velocity distribution of fragments
Detonation mode Velocity range/(m·s−1) Average velocity/(m·s−1) Velocity standard deviation/(m·s−1) Single-endpoint 1 551–1 758 1 660 64 Two-endpoint 1 572–1 735 1 649 37 Three-endpoint 1 563–1 760 1 657 39 -
[1] 李元, 赵倩, 熊诗辉, 等. 一种异面棱柱战斗部威力特性的数值模拟 [J]. 含能材料, 2019, 27(2): 97–103. doi: 10.11943/CJEM2018143LI Y, ZHAO Q, XIONG S H, et al. Numerical modeling on lethality of a faceted prismatic warhead [J]. Chinese Journal of Energetic Materials, 2019, 27(2): 97–103. doi: 10.11943/CJEM2018143 [2] 诸德放, 冯长根, 李友, 等. 基于Gurney假设的一种非对称型战斗部破片初速计算 [J]. 弹箭与制导学报, 2006, 26(1): 74–76. doi: 10.3969/j.issn.1673-9728.2006.01.024ZHU D F, FENG C G, LI Y, et al. Calculation of first velocity of the unsymmetrical warheadʼ fragment based on Gurney hypothesis [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2006, 26(1): 74–76. doi: 10.3969/j.issn.1673-9728.2006.01.024 [3] 王成, 付晓磊, 宁建国. 起爆方式对聚能射流性能影响的数值分析 [J]. 北京理工大学学报, 2006, 26(5): 401–404. doi: 10.3969/j.issn.1001-0645.2006.05.006WANG C, FU X L, NING J G. Numerical simulation of shaped charge jet formation under different ways of initiation [J]. Transactions of Beijing Institute of Technology, 2006, 26(5): 401–404. doi: 10.3969/j.issn.1001-0645.2006.05.006 [4] 武双章, 顾文彬, 李裕春, 等. 环形向内切割聚能装药射流成型性能研究 [J]. 含能材料, 2017, 25(11): 896–902. doi: 10.11943/j.issn.1006-9941.2017.11.003WU S Z, GU W B, LI Y C, et al. Jet formation performance of circular shaped charge of cutting inward [J]. Chinese Journal of Energetic Material, 2017, 25(11): 896–902. doi: 10.11943/j.issn.1006-9941.2017.11.003 [5] 陈伟, 马宏昊, 沈兆武, 等. 起爆方式对环向聚能射流成型影响的数值模拟 [J]. 高压物理学报, 2015, 29(6): 419–424. doi: 10.11858/gywlxb.2015.06.003CHEN W, MA H H, SHEN Z W, et al. Numerical simulation of influence of different modes of initiation on the forming of radial shaped jet [J]. Chinese Journal of High Pressure Physics, 2015, 29(6): 419–424. doi: 10.11858/gywlxb.2015.06.003 [6] 刘健峰, 龙源, 纪冲, 等. 含偏心起爆对EFP战斗部飞行特性的影响 [J]. 爆炸与冲击, 2015, 35(3): 335–342. doi: 10.11883/1001-1455-(2015)03-0335-08LIU J F, LONG Y, JI C, et al. Effect of eccentric initiation on the flight characteristics and ballistic dispersion of EFP [J]. Explosion and Shock Waves, 2015, 35(3): 335–342. doi: 10.11883/1001-1455-(2015)03-0335-08 [7] RONDOT F. Tracking the “poorness” of a linear shaped charge [C]//23rd International Symposium on Ballistic. Tarragona, Spain: International Ballistic Committee, 2007: 169–176. [8] 黄正祥, 张先锋, 陈惠武. 起爆方式对聚能杆式侵彻体成型的影响 [J]. 兵工学报, 2004, 25(3): 289–291. doi: 10.3321/j.issn:1000-1093.2004.03.008HUANG Z X, ZHANG X F, CHEN H W. Influence of modes of detonation on the mechanism of jetting projectile charge [J]. Acta Armamentarii, 2004, 25(3): 289–291. doi: 10.3321/j.issn:1000-1093.2004.03.008 [9] 刘明涛, 汤铁钢, 胡海波, 等. 不同起爆方式下炸药驱动柱壳膨胀断裂的数值模拟 [J]. 爆炸与冲击, 2014, 34(4): 415–420. doi: 10.11883/1001-1455(2014)04-0415-06LIU M T, TANG T G, HU H B, et al. Numerical studies of explosion induced cylindrical shell fracture under different detonating modes [J]. Explosion and Shock Waves, 2014, 34(4): 415–420. doi: 10.11883/1001-1455(2014)04-0415-06 [10] 范中波, 周淑荣, 杭义洪, 等. 爆炸变形式定向战斗部的数值仿真研究 [J]. 兵工学报, 2001, 22(3): 334–337. doi: 10.3321/j.issn:1000-1093.2001.03.013FAN Z B, ZHOU S R, HANG Y H, et al. Numerical investigation of detonatively deformable aimable warhead [J]. Acta Armamentarii, 2001, 22(3): 334–337. doi: 10.3321/j.issn:1000-1093.2001.03.013 [11] 刘琛, 李元, 李燕华, 等. 偏心起爆方式对棱柱形定向战斗部破片飞散规律的影响 [J]. 含能材料, 2017, 25(1): 63–68. doi: 10.11943/j.issn.1006-9941.2017.01.011LIU C, LI Y, LI Y H, et al. Influence of eccentric initiation ways on fragment dispersion rule of prismatic aimable warhead [J]. Chinese Journal of Energetic Materials, 2017, 25(1): 63–68. doi: 10.11943/j.issn.1006-9941.2017.01.011 [12] HELD M. Velocity enhanced warheads [J]. Journal of Explosives and Propellants, 2001, 17(2): 1–12. [13] 周唯潇, 王雅君, 于佳鑫, 等. 起爆方式对复合战斗部毁伤输出的影响 [J]. 高压物理学报, 2021, 35(1): 015101. doi: 10.11858/gywlxb.20200593ZHOU W X, WANG Y J, YU J X, et al. Effect of initial detonation method on damage power of composite warhead [J]. Chinese Journal of High Pressure Physics, 2021, 35(1): 015101. doi: 10.11858/gywlxb.20200593 [14] 鲁忠宝, 黎勤, 哈海荣. 不同能量输出结构战斗部水下爆炸毁伤威力试验研究 [J]. 水下无人系统学报, 2019, 27(1): 71–77.LU Z B, LI Q, HA H R. Experimental study on underwater explosion damage power of warhead with different energy output configuration [J]. Journal of Unmanned Undersea Systems, 2019, 27(1): 71–77. [15] 李贝, 余文力, 王文欣, 等. 中心装药起爆方式对FAE影响的仿真研究 [J]. 兵工自动化, 2014, 33(11): 30–33. doi: 10.7690/bgzdh.2014.11.009LI B, YU W L, WANG W X, et al. Simulation and research on influence of center charge ignition method on FAE [J]. Ordnance Industry Automation, 2014, 33(11): 30–33. doi: 10.7690/bgzdh.2014.11.009 [16] JONATHAN P G, GREG F, COLIN H. Numerical simulation of fragmentation using AUTODYN-2D & 3D in explosive ordnance safety assessment [C]//Proceedings of 6th PARARI International Explosive Ordnance Symposium. Canberra, Australia, 2003. [17] 王力, 韩峰, 陈放, 等. 偏心对称起爆战斗部破片初速的增益 [J]. 爆炸与冲击, 2016, 36(1): 69–74. doi: 10.11883/1001-1455(2016)01-0069-06WANG L, HAN F, CHEN F, et al. Fragmentsʼ velocity of eccentric warhead with double symmetric detonators [J]. Explosion and Shock Waves, 2016, 36(1): 69–74. doi: 10.11883/1001-1455(2016)01-0069-06 -