• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

多破片对柱壳装药冲击起爆速度阈值的数值模拟研究

郭淳 郭尚生 钱建平 顾文彬

引用本文:
Citation:

多破片对柱壳装药冲击起爆速度阈值的数值模拟研究

    作者简介: 郭 淳(1988- ),男,博士研究生,53chun@163.com;
    通讯作者: 钱建平, 13951837475@139.com
  • 中图分类号: O389; TJ410.3

Numerical simulation on shock critical initiation velocity of cylindrical covered charge by multiple fragment impacts

    Corresponding author: QIAN Jianping, 13951837475@139.com ;
  • CLC number: O389; TJ410.3

  • 摘要: 为了研究实战环境中多个钨球破片对导弹战斗部(柱壳装药)的冲击起爆问题,采用AUTODYN-3D数值模拟软件,基于单破片撞击柱壳装药模型,建立多破片撞击柱壳装药的模型,开展了不同钨球个数、空间碰撞位置间隔(撞击角θ、轴向球心距l)及时间间隔对冲击起爆特性影响的数值模拟,获得了带壳B炸药的起爆速度阈值。结果表明:相同条件下,随着钨球个数的增加、空间碰撞位置间隔的减小,起爆速度阈值逐渐减小,6个钨球同时撞击的起爆速度阈值约为单个作用下的50%;双钨球作用下,柱壳装药相较于平板装药更难以起爆;双钨球间隔撞击柱壳装药时,起爆速度阈值均随着撞击时间间隔的增大而先减小后增大,最小起爆速度阈值约为同时撞击时的95%,且|θ2|<|θ1| (θ1为第1个钨球的撞击角,θ2为第2个钨球的撞击角)时更容易起爆柱壳装药。
  • 图 1  单钨球冲击柱壳装药的物理模型和有限元模型

    Figure 1.  Physical model and finite element model of single tungsten fragment impacting on cylindrical charge with steel casing

    图 2  双钨球撞击柱壳装药的简化模型

    Figure 2.  Simplified model of double tungsten fragments impacting on cylindrical charge with casing

    图 3  截面位置示意图

    Figure 3.  Schematic diagram of cross section

    图 4  起爆速度阈值随撞击角的变化曲线

    Figure 4.  Curve of critical velocity vs. impact angle θ

    图 5  多钨球撞击位置方案

    Figure 5.  Multiple tungsten sphere fragment impact options

    图 6  起爆速度阈值随钨球个数的变化

    Figure 6.  Critical velocity vs. tungsten sphere number

    图 7  起爆速度阈值随装药半径倒数的变化

    Figure 7.  Critical velocity vs. l/R

    图 8  双钨球起爆速度阈值随撞击角的变化

    Figure 8.  Critical velocity vs. impact angle θ (double tungsten sphere fragments)

    图 9  柱壳装药和平板装药起爆速度阈值随球心距的变化

    Figure 9.  Critical velocities of cylindrical and plate charge

    图 10  双钨球起爆速度阈值随轴向球心距的变化

    Figure 10.  Critical velocity vs. axial ball center distance l/r(double tungsten sphere fragments)

    图 11  4个钨球打击位置方案

    Figure 11.  Four tungsten fragments combat options

    图 12  “矩形”结构起爆速度阈值随轴向球心距的变化

    Figure 12.  Critical velocity vs. axial ball center distance l/r (rectangle)

    图 13  “菱形”结构起爆速度阈值随轴向球心距的变化

    Figure 13.  Critical velocity vs. axial ball center distance l/r (rhombus)

    图 14  l12=0时起爆速度阈值随时间间隔的变化

    Figure 14.  Critical velocity vs. time interval h2/r (l12=0)

    图 15  观测点θ2-5-2的压力时间曲线

    Figure 15.  Pressure-time curves at observation point θ2-5-2 (v=2 800 m/s)

    图 16  观测点θ2-0-2的压力时间曲线

    Figure 16.  Pressure-time curves of θ2-0-2 (v=2 800 m/s)

    图 17  起爆速度阈值随时间间隔的变化

    Figure 17.  Critical velocities vs. time interval h2/r (l12/r=2)

    表 1  B炸药材料参数

    Table 1.  Material parameters of composition B

    I/μs−1baxG1cdyG2egz
    440.2220.0144140.2220.6672.00000
    下载: 导出CSV

    表 2  破片、壳体材料模型

    Table 2.  Material model of fragment and casing

    部件材料状态方程强度模型失效模型
    破片钨合金ShockJohnson-CookGeometric strain
    壳体4340合金钢LinearJohnson-CookGeometric strain
    下载: 导出CSV
  • [1] 居仙春. 炮射布阵式悬浮弹幕特性及其协同反导效能研究[D]. 南京: 南京理工大学, 2014: 7−8; 34.
    [2] JAMES H R. Critical energy criterion for the shock initiation of explosives by projectile impact [J]. Propellants, Explosives, Pyrotechnics, 1988, 13(2): 35–41. DOI: 10.1002/prep.19880130202.
    [3] HELD M. Initiation phenomena with shaped charge jets [C] // 9th International Symposium on Detonation. Portland, Oregon, US: OCNR, 1989: 1416−1426.
    [4] 陈卫东, 张忠, 刘家良. 破片对屏蔽炸药冲击起爆的数值模拟和分析 [J]. 兵工学报, 2009, 30(9): 1187–1191. DOI: 10.3321/j.issn:1000-1093.2009.09.007.
    CHEN W D, ZHANG Z, LIU J L. Numerical simulation and analysis of shock initiation of shielded explosive impacted by fragments [J]. Acta Armamentarii, 2009, 30(9): 1187–1191. DOI: 10.3321/j.issn:1000-1093.2009.09.007.
    [5] 方青, 卫玉章, 张克明, 等. 射弹倾斜撞击带盖板炸药引发爆轰的条件 [J]. 爆炸与冲击, 1997, 17(2): 153–158.
    FANG Q, WEI Y Z, ZHANG K M, et al. On the projectile oblique-impact initiation conditions for explosive covered with a plate [J]. Explosion and Shock Waves, 1997, 17(2): 153–158.
    [6] 张先锋, 赵有守, 陈惠武. 射弹冲击引爆带壳炸药临界条件 [J]. 弹道学报, 2006, 18(4): 57–59. DOI: 10.3969/j.issn.1004-499X.2006.04.016.
    ZHANG X F, ZHAO Y S, CHEN H W. The critical condition of shelled explosive initiated by projectile [J]. Journal of Ballistics, 2006, 18(4): 57–59. DOI: 10.3969/j.issn.1004-499X.2006.04.016.
    [7] 王树山, 李朝君, 马晓飞, 等. 钨合金破片对屏蔽装药撞击起爆的实验研究 [J]. 兵工学报, 2001, 22(2): 189–191. DOI: 10.3321/j.issn:1000-1093.2001.02.012.
    WANG S S, LI C J, MA X F, et al. An experimental study on the initiation of covered charge impacted by tungsten alloy fragments [J]. Acta Armamentarii, 2001, 22(2): 189–191. DOI: 10.3321/j.issn:1000-1093.2001.02.012.
    [8] 梁争峰, 袁宝慧. 破片撞击起爆屏蔽B炸药的数值模拟和实验 [J]. 火炸药学报, 2006, 29(1): 5–9. DOI: 10.3969/j.issn.1007-7812.2006.01.002.
    LIANG Z F, YUAN B H. Numerical simulation and experimental study of the initiation of shielded composition B impacted by fragment [J]. Chinese Journal of Explosives & Propellants, 2006, 29(1): 5–9. DOI: 10.3969/j.issn.1007-7812.2006.01.002.
    [9] 李文彬, 王晓鸣, 赵国志, 等. 多破片命中时炸药的冲击起爆研究 [J]. 南京理工大学学报(自然科学版), 2004, 28(1): 5–8. DOI: 10.3969/j.issn.1005-9830.2004.01.002.
    LI W B, WANG XM, ZHAO G Z, et al. Study on shock initiation of explosive by the impact of multi-fragment [J]. Journal of Nanjing University of Science and Technology, 2004, 28(1): 5–8. DOI: 10.3969/j.issn.1005-9830.2004.01.002.
    [10] LUECK M, HEINE A, WICKERT M. Numerical analysis of the initiation of high explosives by interacting shock waves due to multiple fragment impact [C] // 26th International Symposium Ballistics. Miami, FL, USA: NDIA, 2011: 73−81.
    [11] 贾宪振, 陈松, 杨建, 等. 双破片同时撞击对B炸药冲击起爆的数值模拟研究 [J]. 高压物理学报, 2011, 25(5): 469–474. DOI: 10.11858/gywlxb.2011.05.014.
    JIA X Z, CHEN S, YANG J, et al. Numerical study of explosives initiation by simultaneous impact from two fragments [J]. Chinese Journal of High Pressure Physics, 2011, 25(5): 469–474. DOI: 10.11858/gywlxb.2011.05.014.
    [12] 梁斌, 冯高鹏, 魏雪婷. 多枚破片冲击引爆带盖板炸药数值模拟分析 [J]. 弹箭与制导学报, 2013, 33(6): 62–66; 69. DOI: 10.3969/j.issn.1673-9728.2013.06.018.
    LIANG B, FENG G P, WEI X T. Numerical simulation on shock initiation of composition explosive of cover board subjected to multi-fragment [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2013, 33(6): 62–66; 69. DOI: 10.3969/j.issn.1673-9728.2013.06.018.
    [13] 宋浦, 梁安定. 破片对柱壳装药的撞击毁伤试验研究 [J]. 弹箭与制导学报, 2006, 26(1): 87–88, 92. DOI: 10.3969/j.issn.1673-9728.2006.01.029.
    SONG P, LIANG A D. Experimental investigation on charges covered shell damaged by fragments impact [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2006, 26(1): 87–88, 92. DOI: 10.3969/j.issn.1673-9728.2006.01.029.
    [14] 江增荣, 李向荣, 李世才, 等. 预制破片对战斗部冲击起爆数值模拟 [J]. 弹道学报, 2009, 21(1): 9–13.
    JIANG Z R, LI X R, LI S C, et al. Numerical simulation on shock initiation of performed fragment to warhead [J]. Journal of Ballistics, 2009, 21(1): 9–13.
    [15] 辛建国, 徐豫新, 李超, 等. 破片冲击柱面薄壳装药实验 [J]. 兵工学报, 2014, 35(S2): 222–227.
    XIN J G, XU Y X, LI C, et al. Experiment of fragment impact on cylinder charge covered with thin shell [J]. Acta Armamentarri, 2014, 35(S2): 222–227.
    [16] 王昕, 蒋建伟, 王树有, 等. 钨球对柱面带壳装药的冲击起爆数值模拟研究 [J]. 兵工学报, 2017, 38(8): 1498–1505. DOI: 10.3969/j.issn.1000-1093.2017.08.006.
    WANG X, JIANG J W, WANG S Y, et al. Numerical simulation on the initiation of cylindrical covered charge impacted by tungsten sphere fragment [J]. Acta Armamentarri, 2017, 38(8): 1498–1505. DOI: 10.3969/j.issn.1000-1093.2017.08.006.
    [17] AUTODYN, Users manual [Z]. California: Century Dynamics Corporation, 2005.
  • [1] 王昕蒋建伟王树有门建兵 . 破片撞击起爆柱面带壳装药的临界速度修正判据. 爆炸与冲击, 2019, 39(1): 012302-1-012302-8. doi: 10.11883/bzycj-2017-0271
    [2] 张涛刘雨生高志鹏杨佳刘艺谷岩 . 垫层对破片冲击起爆带壳炸药影响的数值模拟. 爆炸与冲击, 2018, 38(6): 1241-1246. doi: 10.11883/bzycj-2017-0154
    [3] 章猛华王鹏新余永刚阮文俊王健宁惠君 . 冲击起爆弹丸内装药延迟起爆数值模拟. 爆炸与冲击, 2016, 36(5): 728-733. doi: 10.11883/1001-1455(2016)05-0728-06
    [4] 韩璐韩庆杨爽 . 多破片高速冲击下飞机油箱水锤效应数值模拟. 爆炸与冲击, 2018, 38(3): 473-484. doi: 10.11883/bzycj-2017-0230
    [5] 王桂吉赵同虎莫建军吴刚韩梅谭福利 . 一种以TATB/HMX为基的高聚物粘结炸药的短脉冲冲击起爆特性. 爆炸与冲击, 2007, 27(3): 230-235. doi: 10.11883/1001-1455(2007)03-0230-06
    [6] 陈朗刘群伍俊英 . 受热炸药的冲击起爆特征. 爆炸与冲击, 2013, 33(1): 21-28. doi: 10.11883/1001-1455(2013)01-0021-08
    [7] 皮铮迪陈朗刘丹阳伍俊英 . CL-20基混合炸药的冲击起爆特征. 爆炸与冲击, 2017, 37(6): 915-923. doi: 10.11883/1001-1455(2017)06-0915-09
    [8] 陈少杰吴立志沈瑞琪叶迎华胡艳 . 激光驱动复合飞片冲击起爆HNS-Ⅳ. 爆炸与冲击, 2015, 35(2): 285-288. doi: 10.11883/1001-1455-(2015)02-0285-04
    [9] 李志鹏龙新平黄毅民何碧王蓉何松伟 . 用组合式电磁粒子速度计研究JOB-9003炸药的冲击起爆过程. 爆炸与冲击, 2006, 26(3): 269-272. doi: 10.11883/1001-1455(2006)03-0269-04
    [10] 白志玲段卓平温丽晶张震宇欧卓成黄风雷 . 多元混合PBX炸药冲击起爆的多元Duan-Zhang-Kim反应速率模型研究. 爆炸与冲击, 2019, 39(11): 112101-1-112101-11. doi: 10.11883/bzycj-2018-0410
    [11] 刘海庆段卓平白志玲温丽晶欧卓成黄风雷 . 孔隙度对PBX炸药冲击起爆影响的实验研究. 爆炸与冲击, 2019, 39(7): 072302-1-072302-6. doi: 10.11883/bzycj-2018-0226
    [12] 张忠陈卫东杨文淼 . 非均质固体炸药冲击起爆的物质点法. 爆炸与冲击, 2011, 31(1): 25-30. doi: 10.11883/1001-1455(2011)01-0025-06
    [13] 陶为俊浣石黄风雷蒋国平 . 侧向稀疏波对非均质凝聚炸药冲击波起爆过程的影响. 爆炸与冲击, 2011, 31(4): 397-401. doi: 10.11883/1001-1455(2011)04-0397-05
    [14] 潘昊胡晓棉 . 钝感炸药的超压爆轰与冲击起爆过程数值模拟. 爆炸与冲击, 2006, 26(2): 174-178. doi: 10.11883/1001-1455(2006)02-0174-05
    [15] 庞嵩林陈雄许进升王永平 . 聚能射流对固体火箭发动机的冲击起爆. 爆炸与冲击, 2020, 40(7): 1-10. doi: 10.11883/bzycj-2019-0469
    [16] 梁增友黄风雷张震宇 . PBX-9404的化学反应速率方程及起爆特性. 爆炸与冲击, 2008, 28(1): 38-41. doi: 1.011883/1001-1455(2008)01-0038-06
    [17] 沈慧铭李伟兵王晓鸣李文彬董晓亮 . 圆柱壳体装药偏心多点起爆下破片速度的分布. 爆炸与冲击, 2017, 37(6): 1039-1045. doi: 10.11883/1001-1455(2017)06-1039-07
    [18] 吕军军曾庆轩李明愉周利存 . 起爆高密度TATB炸药的飞片速度阈值. 爆炸与冲击, 2014, 34(1): 125-128. doi: 10.11883/1001-1455(2014)01-0125-04
    [19] 刘明涛汤铁钢胡海波2李庆忠胡秀章李永池 . 不同起爆方式下炸药驱动柱壳膨胀断裂的数值模拟. 爆炸与冲击, 2014, 34(4): 415-420. doi: 10.11883/1001-1455(2014)04-0415-06
    [20] 谯志强聂福德杨光成张娟 . 纳米TATB微结构与复合物冲击起爆阈值间的关系. 爆炸与冲击, 2010, 30(1): 75-79. doi: 10.11883/1001-1455(2010)01-0075-05
  • 加载中
图(17)表(2)
计量
  • 文章访问数:  1450
  • HTML全文浏览量:  963
  • PDF下载量:  19
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-12
  • 录用日期:  2020-03-05
  • 网络出版日期:  2020-05-25
  • 刊出日期:  2020-06-01

多破片对柱壳装药冲击起爆速度阈值的数值模拟研究

    作者简介:郭 淳(1988- ),男,博士研究生,53chun@163.com
    通讯作者: 钱建平, 13951837475@139.com
  • 1. 南京理工大学机械工程学院,江苏 南京 210094
  • 2. 辽沈工业集团有限公司产品研发中心,辽宁 沈阳 110045
  • 3. 中国人民解放军陆军工程大学,江苏 南京 210007

摘要: 为了研究实战环境中多个钨球破片对导弹战斗部(柱壳装药)的冲击起爆问题,采用AUTODYN-3D数值模拟软件,基于单破片撞击柱壳装药模型,建立多破片撞击柱壳装药的模型,开展了不同钨球个数、空间碰撞位置间隔(撞击角θ、轴向球心距l)及时间间隔对冲击起爆特性影响的数值模拟,获得了带壳B炸药的起爆速度阈值。结果表明:相同条件下,随着钨球个数的增加、空间碰撞位置间隔的减小,起爆速度阈值逐渐减小,6个钨球同时撞击的起爆速度阈值约为单个作用下的50%;双钨球作用下,柱壳装药相较于平板装药更难以起爆;双钨球间隔撞击柱壳装药时,起爆速度阈值均随着撞击时间间隔的增大而先减小后增大,最小起爆速度阈值约为同时撞击时的95%,且|θ2|<|θ1| (θ1为第1个钨球的撞击角,θ2为第2个钨球的撞击角)时更容易起爆柱壳装药。

English Abstract

  • 如何有效拦截来袭导弹,一直是反导弹药研究的热点问题。悬浮弹幕协同反导弹药是利用其多个子弹药战斗部协同作用产生的高速破片场毁伤来袭导弹[1],此破片场由成百上千个飞散方向任意的破片(预制破片或爆炸成型弹丸)组成。应用具有一定分布密度的破片撞击来袭导弹并引发导弹战斗部起爆是毁伤导弹最有效的方法之一。

    针对破片冲击起爆战斗部的问题,大多数学者是将其等效为单破片或多破片冲击起爆平板装药的问题进行研究。单破片冲击起爆平板装药的问题,学者们已通过理论、数值模拟和试验等方法进行了大量研究[2-8]。在单破片冲击起爆平板装药的基础上,一些学者也曾开展了多破片冲击起爆平板装药的研究。其中李文彬等[9]研究了双破片冲击起爆裸装药的机理;Lueck等[10]利用4组不同的参数,论证了Lee-Tarver点火增长模型模拟多破片冲击起爆炸药的可行性;贾宪振等[11]、梁斌等[12]通过数值模拟研究了多破片撞击平板装药的情况,得到了起爆速度阈值的变化规律。

    然而,现实中战斗部多为圆柱结构,一些学者曾在试验和数值模拟过程中考虑了战斗部为圆柱结构的情况,但都是采用单破片进行研究。其中宋浦等[13]针对某型破片式战斗部对带壳柱装药模拟靶的撞击毁伤试验研究,得到柱壳装药的毁伤阈值范围;江增荣等[14]通过数值模拟得到重金属材料破片可以提高毁伤战斗部的能力;辛建国等[15]研究得到破片及破碎的柱面薄壳对壳内装药的撞击、摩擦等机械作用是引燃/引爆薄壳装药的主要原因;王昕等[16]通过数值模拟分析了单个破片冲击起爆柱壳装药的规律,得到命中位置和装药曲率半径对柱壳装药冲击起爆的影响。

    在实战环境中,战斗部很有可能受到多个破片的累积毁伤作用,从现有文献可以看出,关于多破片冲击起爆柱壳装药问题的研究尚不深入全面,因此本文中将采用数值模拟的方法研究多钨球破片对柱壳装药的冲击起爆问题。采用AUTODYN-3D数值模拟软件,在单钨球破片撞击柱壳装药模型的基础上,建立多钨球破片撞击柱壳装药的模型,分析当存在一定空间间隔(撞击角、轴向球心距)和时间间隔的条件下,多钨球破片冲击柱壳装药的起爆过程,得到相关参数(钨球个数、装药半径、撞击角、轴向球心距和时间间隔)对起爆柱壳装药速度阈值的影响规律。研究结果可为反导弹药战斗部的设计和多破片累积毁伤评估提供参考。

    • 为了建立多钨球破片撞击柱壳装药的模型,先建立单个钨合金球破片冲击起爆带钢壳的圆柱B炸药的模型,如图1(a) 所示,进行数值模拟。其中钨球质量为m,半径为r;炸药装药半径为R,长度为L;壳体厚度为H。假设钨球垂直于带壳装药轴线入射。

      图  1  单钨球冲击柱壳装药的物理模型和有限元模型

      Figure 1.  Physical model and finite element model of single tungsten fragment impacting on cylindrical charge with steel casing

      采用AUTODYN-3D程序,简化计算,建立1/2模型,通过对称性实现全模型,如图1(b) 所示。其中,钨球、钢壳和B炸药均采用Lagrange网格,网格数量分别为2 048、28 800、192 000。

      钨球、钢壳和B炸药的材料参数均从AUTODYN材料库[17]选取。B炸药的状态方程采用Lee-Tarver三项反应过程炸药起爆模型,该模型可以很好地模拟非均质炸药的冲击起爆特性:

      $ \frac{{\partial F}}{{\partial t}} = I{\left( {1 - F} \right)^b}{\left( {\frac{\rho }{{{\rho _0}}} - 1 - a} \right)^x} + {G_1}{\left( {1 - F} \right)^c}{F^d}{p^y} + {G_2}{\left( {1 - F} \right)^e}{F^g}{p^z} $

      式中:F为反应比(炸药气体质量与炸药总质量比),ρρ0为炸药当前密度和初始密度,p为炸药爆轰压力,其余参数为常数。式(1)中右边第1项为点火项,描述了热点燃烧阶段;第2项为成长项,描述了燃烧从热点向内部和外部的颗粒传递阶段;第3项为快反应项,描述了燃烧快速向爆轰转变的过程。炸药反应物和未反应物均采用JWL状态方程。B炸药材料参数见表1。破片和柱壳的材料状态方程、强度模型和失效方程见表2

      I/μs−1baxG1cdyG2egz
      440.2220.0144140.2220.6672.00000

      表 1  B炸药材料参数

      Table 1.  Material parameters of composition B

      部件材料状态方程强度模型失效模型
      破片钨合金ShockJohnson-CookGeometric strain
      壳体4340合金钢LinearJohnson-CookGeometric strain

      表 2  破片、壳体材料模型

      Table 2.  Material model of fragment and casing

    • 在单破片冲击起爆柱壳装药模拟模型的基础上,增加钨球破片数量,建立如图2所示的双钨球撞击柱壳装药的简化模型,其中双钨球完全相同,质量均为m=3 g;两球间平行于装药轴线方向的球心距即轴向球心距为l。假设钨球均垂直于带壳装药轴线入射。图3为任意一个钨球碰撞点的截面示意图。其中,AA′、BB′为圆柱装药相互垂直的两个对称面,双破片撞击时,假设碰撞点在AA′面的同一侧;O0Oi0i=1,2,···,最大值为钨球总个数)分别为装药的轴心和初始钨球球心,Oi为碰撞瞬间的钨球球心;钨球的初始速度为vi,钨球从初始位置运动到碰撞点的距离为hi,忽略钨球在飞行过程中的速度衰减,即碰撞速度亦为viO0Oi0的连线与vi的夹角为撞击角θi;设θi在∠AO0B范围内取正值,以步长为5°,从0°到60°设置一系列观测序列,列号设为θi-j(j=0~12),在每列上设置间距为5 mm的观测点θi-j-kk=1~9);设θi在∠AO0B范围内取负值,以步长为−5°,从0°到−60°设置一系列观测序列,列号设为θi-j(j=0~−12),在每列上设置间距为5 mm的观测点θi-j-kk=1~9),其中点θi-j-1在装药表面,点θi-j-9O0重合。

      图  2  双钨球撞击柱壳装药的简化模型

      Figure 2.  Simplified model of double tungsten fragments impacting on cylindrical charge with casing

      图  3  截面位置示意图

      Figure 3.  Schematic diagram of cross section

      综上,通过l/rθi表征两钨球的空间碰撞位置间隔,hi/r表征两钨球的时间间隔。

    • 为了研究多钨球累积冲击起爆柱壳装药的问题,先对单钨球在不同位置冲击起爆柱壳装药的情况进行数值模拟,以10 m/s为步长,得到装药半径R=40 mm时,单钨球起爆柱壳装药的速度阈值随撞击角θ的变化如图4所示。由图4可以看出,当0°≤θ<25°时,起爆装药的速度阈值随θ的增大而缓慢增长;当25°≤θ<55°时,起爆装药的速度阈值随θ的增大而迅速增长;当θ≥55°时,钨球发生跳飞,本文中不考虑跳飞的情况。简单来说,起爆速度阈值随撞击角增大呈指数增长是因为钨球破片速度方向上的壳体厚度的不同引起的。撞击角度越大,钨球速度方向上的壳体越厚,因此需要更大的起爆速度。

      图  4  起爆速度阈值随撞击角的变化曲线

      Figure 4.  Curve of critical velocity vs. impact angle θ

      本文中采用的单钨球模型的撞击角度与起爆柱壳装药的速度阈值间的变化关系,采用最小二乘法得到拟合公式为:

      $ {v_{{i}}} = 2\;809 + 60.93{{\rm e}^{{\theta }/{{18.1}}}}\;\;\;\;\;\;\;\;0^\circ {\text{≤}} \theta {\text{<}} 55^\circ $

      式(2)可作为研究多钨球起爆柱壳装药的基础,方便获得其中任意一个钨球的碰撞位置对应此位置的单钨球冲击起爆柱壳装药的速度阈值。

      θ=0°时,钨球沿柱壳装药的径向碰撞时的起爆速度阈值为2 840 m/s,与文献[16]中数值的误差为0.5%,表明本文中采用的模拟模型具有较好的吻合性,采用冲击起爆机制开展引爆机理研究是合理的。当d=7.6 mm、H=6 mm时,数值模拟得到起爆速度阈值为2 700 m/s,对比文献[14]的试验值2 676 m/s,误差为0.9%,表明在改变钨球尺寸时,该模型仍有较好的适用性。因此可在单破片冲击起爆柱壳装药模拟模型的基础上,增加钨球破片数量,分析多钨球破片对柱壳装药冲击起爆的影响。

    • 在实战环境中,导弹战斗部可能受到多个高速破片的同时打击,图5展示了1、2、3、4、6个钨球的撞击位置方案,钨球关于柱壳装药的法线成轴对称,两两钨球间的球心距为2r图6为5种方案下多钨球同时冲击柱壳装药的起爆速度阈值的变化曲线,可以看出当钨球个数从1个增加到4个时,起爆柱壳装药的速度阈值变化明显,4个钨球作用时的速度阈值约是单个的55%;当从4个增加到6个时速度阈值的变化曲线趋于平缓,6个钨球作用时的速度阈值约是单个的50%。随着钨球个数的增加,钨球撞击柱壳装药产生的冲击波数量增加,使得多道冲击波作用下更容易叠加成更强的冲击波,所以随着钨球个数的增加,起爆速度阈值逐渐减小。而随着钨球个数的增加,钨球球心距离撞击位置的轴对称中心更远,在到达叠加位置时每道冲击波衰减得更多,所以钨球个数从4个增加到6个时,个数对速度阈值的影响已经不明显。

      图  5  多钨球撞击位置方案

      Figure 5.  Multiple tungsten sphere fragment impact options

      图  6  起爆速度阈值随钨球个数的变化

      Figure 6.  Critical velocity vs. tungsten sphere number

    • 由于柱壳装药柱面的存在,在研究双钨球紧靠在一起的情况时分为两种情况:一是关于BB′平面对称,并沿着平行此平面的方向撞击柱壳装药的情况;二是双钨球同时碰撞柱壳装药同一母线的情况。得到这两种情况下的柱壳装药起爆速度阈值随装药半径倒数的变化,如图7所示。由于选取的计算步长为10 m/s,可以看出针对不同的装药半径,起爆柱壳装药的速度阈值没有明显差异,这是由于在双钨球紧靠在一起时,不同装药半径的带壳装药在两种对称情况下的钨球速度方向的壳体厚度与H近似相等,因此可以认为此时装药半径对速度阈值没影响。

      图  7  起爆速度阈值随装药半径倒数的变化

      Figure 7.  Critical velocity vs. l/R

    • 双钨球冲击起爆柱壳装药时,双钨球的空间碰撞位置可分为两种情形:一是两个碰撞位置关于柱壳装药的任意对称面而对称;二是两个碰撞位置不关于柱壳装药的任意对称面对称。

      当双钨球同时撞击柱壳装药,碰撞位置关于平面BB′对称(双钨球的撞击角θ1=−θ2、双钨球轴向球心距l =l12=0)时,装药起爆速度阈值随撞击角的增大而增大,起爆速度阈值的最小值约为对应位置单个钨球作用的76%,如图8所示。当|θ1|=|θ2|≥20°时,随着撞击角的增大起爆速度阈值已近似等于单个钨球的速度阈值。

      图  8  双钨球起爆速度阈值随撞击角的变化

      Figure 8.  Critical velocity vs. impact angle θ (double tungsten sphere fragments)

      双钨球同时碰撞柱壳装药同一母线(θ1=θ2=0°)时,装药起爆速度阈值随钨球球心距离的变化如图9所示。由前文知,当双钨球紧靠在一起时(即l=2r),柱壳装药与平板装药具有相同的阈值起爆速度。随着间隔距离的增大,速度阈值增大,并最终趋近于单个钨球的速度阈值,起爆速度阈值的最小值约为单个钨球作用的77%。当球心距在7倍钨球半径范围内,即速度小于单个钨球撞击柱壳装药的起爆速度阈值时,相同的球心距条件下,起爆柱壳装药的速度阈值不小于平板装药,即柱壳装药相对于同厚度的平板装药更难以起爆,这是由于此时冲击波在柱壳装药中传播的损耗比在平板装药中多。对于柱壳装药,球心距为4倍钨球半径时,双钨球的起爆速度阈值约为单个的92%,即l<4r时双钨球起爆柱壳装药的叠加作用明显。

      图  9  柱壳装药和平板装药起爆速度阈值随球心距的变化

      Figure 9.  Critical velocities of cylindrical and plate charge

      而当双钨球碰撞位置不关于柱壳装药的任意对称面而对称时,选取双钨球分别以撞击角0°和20°撞击柱壳装药时,装药起爆速度阈值随球心距的变化如图10所示,双钨球的起爆速度阈值的最小值约为单钨球以0°撞击时的95%。可以看出当|θ1|≠|θ2|、l≠0时,即当两个钨球撞击角度一定时,圆柱装药的起爆速度阈值随球心距的增加而增大,最终趋于撞击角度绝对值更小的钨球起爆装药的速度阈值。

      图  10  双钨球起爆速度阈值随轴向球心距的变化

      Figure 10.  Critical velocity vs. axial ball center distance l/r(double tungsten sphere fragments)

    • 4个钨球同时撞击柱壳装药位置的方案如图11所示,假设4个钨球同时撞击,撞击形状近似为“矩形”和“菱形”。

      图  11  4个钨球打击位置方案

      Figure 11.  Four tungsten fragments combat options

      图12是4个钨球呈“矩形”撞击柱壳装药在不同的空间间隔作用下,装药的起爆速度阈值随轴向球心距(钨球1、2的轴向球心距与钨球3、4的轴向球心距相等时,即l=l12=l34)的变化。当撞击角度相同时,随着轴向球心距l的增大,起爆速度阈值增大。而当钨球轴向球心距l相同时,随着撞击角的增大,起爆速度阈值增大。由图6图12可以得出,4个钨球以“矩形”撞击时的起爆速度阈值的最小值约为撞击角度绝对值最小的单钨球作用时的54%。

      图  12  “矩形”结构起爆速度阈值随轴向球心距的变化

      Figure 12.  Critical velocity vs. axial ball center distance l/r (rectangle)

      图13是4个钨球呈“菱形”撞击柱壳装药在不同的空间间隔作用下,装药的起爆速度阈值随钨球2、4轴向球心距(l=l24)的变化。当撞击角度相同时,随着轴向球心距l的增大,起爆速度阈值增大;而当钨球间轴向球心距l相同时,随着撞击角的增大,起爆速度阈值增大。可以看出,4个钨球间排得越紧密,即空间间隔越小,装药的起爆速度阈值越小,4个钨球以“菱形”撞击时的起爆速度阈值的最小值约为撞击角度绝对值最小的单钨球作用时的55%。

      图  13  “菱形”结构起爆速度阈值随轴向球心距的变化

      Figure 13.  Critical velocity vs. axial ball center distance l/r (rhombus)

      因此,当多钨球破片同时撞击柱壳装药时,起爆速度阈值随空间碰撞位置间隔的增大而增大,这是由于间隔的增大造成每个破片产生的冲击波到达叠加位置时衰减得更多。

    • 为了研究时间间隔对起爆速度阈值的影响,也可分为两种情况分析:一是双钨球破片的球心在同一个圆柱截面;二是双钨球破片的球心不在同一个圆柱截面。

      假设双钨球的球心在同一个圆柱截面,即l12=0。分别研究不同撞击角下,速度阈值随时间间隔(即h2/r)的变化;当2个钨球碰撞位置相同、而碰撞顺序不同时,速度阈值随h2/r的变化。

      图14可以看出,3种情况均是钨球1先碰撞柱壳装药,钨球2后碰撞,当双钨球撞击角分别为θ1=0°、θ2=10°时,装药起爆速度阈值随撞击时间间隔的增大而增大,最后达到单个钨球的速度阈值;而在θ1=0°、θ2=20°时,装药起爆速度阈值随碰撞时间间隔的增大先减小后增大;在θ1=20°、θ2=0°时,装药起爆速度阈值随碰撞时间间隔的增大先减小后增大。θ1=0°、θ2=20°与θ1=20°、θ2=0°的2条变化曲线的趋势基本保持一致,都是先减小后增大,最小的起爆速度阈值约为同时撞击的97%。

      图  14  l12=0时起爆速度阈值随时间间隔的变化

      Figure 14.  Critical velocity vs. time interval h2/r (l12=0)

      θ1=0°、θ2=10°时,由于两个钨球距离较近,钨球2撞击瞬间,钨球1已经侵彻一段距离,而钨球1对于钨球2的侵彻过程造成干扰,所以随着撞击时间间隔的增大需要更高的起爆速度阈值。本文不侧重研究这种互相干扰的情况。当θ1=0°、θ2=20°、v=2 800 m/s时,选取h2=7,9,14 mm,提取不同计算文件中同一观测点θ2-5-2的压力-时间曲线,如图15所示。可以看出,当h2=9 mm时,峰值压力最大,其次是h2=7 mm,而h2=14 mm时压力最小,即随着h2的增大,峰值压力先增大后减小。当θ1=20°、θ2=0°、v=2 800 m/s时,同样选取h2=7,9,14 mm,提取不同计算文件中同一观测点θ2-0-2的压力-时间曲线,如图16所示。从图16可以看出,当h2=9 mm时峰值压力也是最大,同样可以得到随着h2的增大,峰值压力先增大后减小。对比图15图16,在相同时间间隔下,图16中观测点的峰值压力均高于图15中对应点的峰值压力。因此,采用相同的撞击速度v=2 800 m/s时,无论钨球撞击的先后顺序怎样,随着时间间隔的增大,炸药内同一位置的压力值均先增大后减小,因此,起爆装药的速度阈值随时间间隔的增大先减小后增大。考虑钨球撞击先后顺序时,在|θ2|<|θ1|的情况下更容易起爆柱壳装药,产生这种变化趋势的原因,可以借鉴一维平面冲击波的相互作用来解释。先减小是由于第2个钨球破片撞击时,第1个钨球破片的冲击波没扫过第2个钨球的撞击区域,所以两道冲击波迎面相遇发生叠加。后增大是由于第2个钨球破片撞击时,第1个钨球破片的冲击波扫过第2个破片的撞击区域,而冲击波的传播相对于波前介质是超音速的,而相对于波后介质是亚音速的,因此第2个钨球破片形成的冲击波必定追赶上第1道冲击波,叠加成更强的一道冲击波。而追赶过程相较于迎面相遇过程,两破片的冲击波衰减得更多,所以就需要更高的起爆速度阈值。在考虑撞击先后顺序时,以相同速度撞击时,后撞击的撞击角小于先撞击的角度时,后一个钨球速度方向的壳体厚度更薄,冲击波传播过程衰减相对少,更容易与前一道冲击波叠加成更强的一道冲击波,所以|θ2|<|θ1|时更容易起爆柱壳装药。

      图  15  观测点θ2-5-2的压力时间曲线

      Figure 15.  Pressure-time curves at observation point θ2-5-2 (v=2 800 m/s)

      图  16  观测点θ2-0-2的压力时间曲线

      Figure 16.  Pressure-time curves of θ2-0-2 (v=2 800 m/s)

      当两个钨球球心不在同一圆柱截面时,选取l12/r=2进行研究,钨球1先碰撞柱壳装药,钨球2后碰撞。当θ1=0°、θ2=20°、l12/r=2与θ1=20°、θ2=0°、l12/r=2时,装药起爆速度阈值均随碰撞时间间隔的增大先减小后增大,最小的起爆速度阈值约为同时撞击的95%,如图17所示。当2个钨球的撞击位置相同、撞击顺序不同时,曲线的整体趋势是相同的,都是先减小后增大。在h/r<7范围内,|θ2|<|θ1|时的起爆速度阈值小于|θ2|>|θ1|时的速度阈值,即当2个钨球撞击角度一定时,后撞击的撞击角小于先撞击的角度时,更容易起爆柱壳装药。

      图  17  起爆速度阈值随时间间隔的变化

      Figure 17.  Critical velocities vs. time interval h2/r (l12/r=2)

    • 通过AUTODYN-3D软件对多钨球破片冲击起爆柱壳装药特性的研究,得到以下结论。

      (1)相同球心距情况下,以柱壳装药的法线为对称轴,随着钨球个数的增加,起爆柱壳装药的速度阈值明显减小,4个钨球作用时的速度阈值约是单个钨球以0°撞击的55%,6个钨球作用时的速度阈值约是单个钨球以0°撞击的50%。

      (2)双钨球紧靠在一起撞击柱壳装药时,装药半径对装药起爆速度阈值影响不明显。当以θ1=θ2=0°撞击时,柱壳装药和平板装药有明显的差别,当l≤7r时,柱壳装药相较于平板装药更难以起爆,且l<4r时双钨球起爆柱壳装药的叠加作用明显。

      (3)多钨球同时撞击柱壳装药的起爆速度阈值均随着撞击角、轴向球心距的增大而增大,最终趋近于多钨球破片中撞击角绝对值最小的那枚破片碰撞位置对应的单钨球起爆装药的速度阈值。双钨球作用时起爆速度阈值的最小值约为撞击角绝对值最小的单钨球作用的76%,4个钨球作用时起爆速度阈值的最小值约为撞击角绝对值最小的单钨球作用时的54%。

      (4)双钨球间隔撞击柱壳装药的起爆速度阈值均随撞击时间间隔的增大而先减小后增大,最小起爆速度阈值约为同时撞击的95%,且|θ2|<|θ1|的方案更容易起爆柱壳装药。

参考文献 (17)

目录

    /

    返回文章
    返回