近场爆炸冲击波对屏蔽压装TNT的冲击引爆试验和仿真

陈兴旺 王金相 唐奎 陈日明 周莲 郝春杰

引用本文:
Citation:

近场爆炸冲击波对屏蔽压装TNT的冲击引爆试验和仿真

    作者简介: 陈兴旺(1994-),男,博士研究生,主要从事屏蔽装药结构安全性研究. E-mail: m18851199809@163.com;
    通讯作者: 王金相, wjx@njust.edu.cn
  • 中图分类号: O383.1

Experimental and Numerical Study of Shock Initiation of Covered TNT by Near-Field Shock Wave

    Corresponding author: WANG Jinxiang, wjx@njust.edu.cn ;
  • CLC number: O383.1

  • 摘要: 为研究近场强冲击波对屏蔽装药的冲击引爆效应,通过试验和LS-DYNA3D非线性有限元计算程序对屏蔽压装TNT的冲击引爆进行研究,得到了接触爆炸时压装TNT被引爆的临界屏蔽板厚度和非接触爆炸时屏蔽压装TNT的殉爆距离,分析了非接触爆炸时屏蔽板厚度对殉爆距离的影响,并通过非线性最小二乘法拟合得到屏蔽板厚度与殉爆距离的函数关系。结果表明:数值计算结果与试验结果基本一致;接触爆炸时临界起爆的屏蔽板厚度在20~23 mm之间;屏蔽板采用3 mm厚的45钢时,非接触爆炸作用下压装TNT的殉爆距离在12~15 mm之间;非接触爆炸时殉爆距离随着屏蔽板厚度的增加而减小,当无屏蔽板时,压装TNT的殉爆距离为79 mm,当屏蔽板厚度从1 mm加大到9 mm时,殉爆距离从51 mm减为1.5 mm,被发装药的屏蔽板对冲击波有显著的防护作用。
  • 图 1  屏蔽装药冲击引爆数值计算模型图

    Figure 1.  Numerical calculation model of covered charge explosive detonation

    图 2  屏蔽装药冲击引爆实物照片

    Figure 2.  Photo of shock initiation of covered charge

    图 3  套筒高度为10 mm时不同时刻起爆装置的高速摄影照片

    Figure 3.  High speed camera photos of experimental setup with 10 mm PPR sleeve at different time

    图 4  套筒高度为15 mm时不同时刻起爆装置的高速摄影照片

    Figure 4.  High speed camera photos of experimental setup with 15 mm PPR sleeve at different time

    图 5  接触爆炸试验后被发装药的见证板变形

    Figure 5.  Deformation of the witness plates under the acceptor with contact explosion

    图 6  非接触爆炸试验后被发装药的见证板变形

    Figure 6.  Deformation of the witness plates under the acceptor with non-contact explosion

    图 7  被发装药爆炸后见证板变形的数值模拟结果

    Figure 7.  Numerical simulation results of witness plate deformation

    图 8  13.0 mm距离下炸药不同时刻的变形

    Figure 8.  Calculated pressure contours at various time for the configuration with 13 mm separation

    图 9  13.5 mm距离下炸药不同时刻的变形

    Figure 9.  Calculated pressure contours at various times for the configuration with 13.5 mm separation

    图 10  13.0 mm间隔高度下被发装药的爆轰成长过程

    Figure 10.  Detonation development process with spacing distance of 13.0 mm

    图 11  13.5 mm间隔高度下被发装药的熄爆过程

    Figure 11.  Extinguished process with spacing distance of 13.5 mm

    图 13  殉爆距离随屏蔽板厚度变化曲线

    Figure 13.  Variation of sympathetic detonation distance with shield thickness

    图 12  屏蔽装药冲击起爆示意

    Figure 12.  Shock initiation of covered charge

    表 1  主发炸药和屏蔽板的材料模型参数[11-13]

    Table 1.  Parameters for donor charge and shell model[11-13]

    Material ρ0/(g·cm-3 D/(m·s-1 pCJ/GPa A/GPa B/GPa R1 R2 ω E/(109 J·m-3)
    TNT 1.60 6 812 18.56 370 3.23 4.15 0.95 0.3 7
    Material ρ0/(g·cm-3 G/GPa A1/GPa B1/GPa C1 n m Tm/K Tr/K
    45 steel 7.85 82.3 0.507 0.320 0.280 0.064 1.06 1 765 298
    Q235 steel 7.83 80.0 0.235 0.250 0.022 1.000 0.36 1 793 294
    下载: 导出CSV

    表 2  TNT点火增长反应模型参数[13]

    Table 2.  Ignition and growth reactive flow model parameters for TNT[13]

    I/μs–1 G1/( Pa–3·s–1) G2/( Pa–3·s–1) a b c d e g x y z
    8×108 4.20×10–25 2.60×10–25 0.111 0.667 0.667 0.667 0.333 1 6 3 3
    下载: 导出CSV

    表 3  屏蔽装药冲击引爆试验结果

    Table 3.  Results of shock initiation of covered charge

    Explosion type Test No. d/mm h/mm Explosion situation of acceptor charge Pit depth/mm
    Contact explosion 1 20 Exploded 4.22
    2 26 Unexploded 0
    3 23 Unexploded 0
    Non-contact explosion 4 3 10 Exploded 4.20
    5 3 19 Unexploded 0
    6 3 15 Unexploded 0
    7 3 12 Part of the explosion 2.20
    下载: 导出CSV

    表 4  屏蔽装药冲击引爆试验和仿真计算结果

    Table 4.  Test and simulation results of shock initiation of covered charge

    Explosion type Test results Simulation results
    d/mm Explosion situation d/mm Explosion situation
    Contact explosion 20 Exploded 24.0 Exploded
    23 Unexploded 24.5 Unexploded
    Explosion type Test results Simulation results
    h/mm Explosion situation h/mm Explosion situation
    Non-contact explosion 12 Partial explosion 13.0 Exploded
    15 Unexploded 13.5 Unexploded
    下载: 导出CSV
  • [1] 徐松林, 阳世清, 李玉龙. 弹药与爆炸品的安全标准体系 [J]. 安全与环境工程, 2012, 19(5): 85–89 doi: 10.3969/j.issn.1671-1556.2012.05.021
    XU S L, YANG S Q, LI Y L. Safety standard system of ammunitions and explosives [J]. Safety and Environmental Engineering, 2012, 19(5): 85–89 doi: 10.3969/j.issn.1671-1556.2012.05.021
    [2] ROSENBERG Z, DEKEL E. Terminal ballistics [M]. Berlin: Springer, 2012.
    [3] 陈卫东, 张忠, 刘家良. 破片对屏蔽炸药冲击起爆的数值模拟和分析 [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
    [4] 童宗保, 王金相, 彭楚才, 等. 预制破片对屏蔽炸药冲击引爆研究 [J]. 科学技术与工程, 2014, 14(7): 173–177 doi: 10.3969/j.issn.1671-1815.2014.07.038
    TONG Z B, WANG J X, PENG C C, et al. Study on the initiation of shielded explosive impacted by prefabricated fragment [J]. Science Technology and Engineering, 2014, 14(7): 173–177 doi: 10.3969/j.issn.1671-1815.2014.07.038
    [5] 杨正才, 廖昕, 李晓刚, 等. 某隔板起爆器冲击起爆过程的数值模拟与实验研究 [J]. 含能材料, 2011, 19(2): 221–225 doi: 10.3969/j.issn.1006-9941.2011.02.022
    YANG Z C, LIAO X, LI X G, et al. Numerical simulation and experimental investigation for shock initiation of bulkhead initiator [J]. Chinese Journal of Energetic Materials, 2011, 19(2): 221–225 doi: 10.3969/j.issn.1006-9941.2011.02.022
    [6] 周杰, 何勇, 何源, 等. 含能毁伤元冲击引爆模拟战斗部试验研究 [J]. 含能材料, 2016, 24(11): 1048–1056 doi: 10.11943/j.issn.1006-9941.2016.11.003
    ZHOU J, HE Y, HE Y, et al. Experimental study on shock initiation of simulative warhead by energetic kill element [J]. Chinese Journal of Energetic Materials, 2016, 24(11): 1048–1056 doi: 10.11943/j.issn.1006-9941.2016.11.003
    [7] WALKER F E, WASLEY R J. Critical energy shock initiation of heteropterous explosive [C]//Explosive Pyrotechnics, 1968.
    [8] LU J P, LI J, KENNEDY D L. Simulation of sympathetic reaction rests for PBXN-109 [C]//Proceedings of 13th Symposium (International) on Detonation. New York, USA, 2006: 1338–1349.
    [9] 王晨, 伍俊英, 陈朗, 等. 壳装炸药殉爆实验和数值模拟 [J]. 爆炸与冲击, 2010, 30(2): 152–158 doi: 10.11883/1001-1455(2010)02-0152-07
    WANG C, WU J Y, CHEN L, et al. Experiments and numerical simulations of sympathetic detonation of explosives in shell [J]. Explosion and Shock Waves, 2010, 30(2): 152–158 doi: 10.11883/1001-1455(2010)02-0152-07
    [10] 时党勇, 李裕春, 张胜民. 基于ANSYS/LS-DYNA8.1进行显示动力分析 [M]. 北京: 清华大学出版社, 2005.
    [11] 陈刚, 陈忠富, 陶俊林, 等. 45钢的Johnson-Cook失效模型 [C]//中国力学学会学术大会, 2005.
    [12] 张磊. Q235钢拉伸性能的有限元仿真与试验 [D]. 长春: 吉林大学, 2012.
    [13] 蒋国平. 非均质凝聚炸药二维冲击起爆实验与Lagrange分析研究 [D]. 长沙: 湖南大学, 2005.
    [14] 北京工业学院八系. 爆炸及其作用(上) [M]. 北京: 国防工业出版社, 1982.
  • [1] 李金河赵继波谭多望王彦平张远平 . 不同起爆方式对含铝炸药水中爆炸近场冲击波压力的影响. 高压物理学报, 2012, 26(3): 289-293. doi: 10.11858/gywlxb.2012.03.007
    [2] 王长利周刚马坤陈春林赵南冯娜 . 聚能装药水下爆炸冲击波载荷规律. 高压物理学报, 2017, 31(4): 453-461. doi: 10.11858/gywlxb.2017.04.014
    [3] 池家春刘雨生龚晏清李强 . JB9014炸药在常温和-54 ℃低温下冲击引爆压力场发展的实验研究. 高压物理学报, 2001, 15(1): 39-47 . doi: 10.11858/gywlxb.2001.01.006
    [4] 蒋海燕李芝绒张玉磊苏健军 . 运动装药空中爆炸冲击波特性研究. 高压物理学报, 2017, 31(3): 286-294. doi: 10.11858/gywlxb.2017.03.010
    [5] 章冠人 . 冲击波作用下铁电体的击穿. 高压物理学报, 1995, 9(1): 1-4 . doi: 10.11858/gywlxb.1995.01.001
    [6] 尤文斌马铁华丁永红崔敏张晋业 . 冲击波测试系统的动态建模及应用. 高压物理学报, 2014, 28(4): 429-434. doi: 10.11858/gywlxb.2014.04.007
    [7] 章冠人 . 冲击波后温度测量的理论问题. 高压物理学报, 1992, 6(4): 241-246 . doi: 10.11858/gywlxb.1992.04.001
    [8] 金孝刚王弘 . 冲击波处理45钢的微结构研究. 高压物理学报, 1993, 7(4): 254-259 . doi: 10.11858/gywlxb.1993.04.003
    [9] 章冠人 . 冲击波作用下材料的极化率. 高压物理学报, 1990, 4(3): 161-166 . doi: 10.11858/gywlxb.1990.03.001
    [10] 姚熊亮杨文山初文华张阿漫 . 水中空气隔层衰减冲击波性能研究. 高压物理学报, 2011, 25(2): 165-172 . doi: 10.11858/gywlxb.2011.02.013
    [11] 胡宏伟宋浦郭炜冯海云张立建 . 地面爆炸冲击波的相互作用. 高压物理学报, 2014, 28(3): 353-357. doi: 10.11858/gywlxb.2014.03.014
    [12] 贺红亮 . 立方氮化硅及其冲击波合成. 高压物理学报, 2007, 21(1): 1-7 . doi: 10.11858/gywlxb.2007.01.001
    [13] 孙悦袁长迎张修路吴国栋 . 铝阳极化冲击波探针的工艺探索. 高压物理学报, 2004, 18(2): 157-162 . doi: 10.11858/gywlxb.2004.02.011
    [14] 卢芳云蔡清裕周新贵陈刚周安郴 . 氮化硅粉末的冲击波活化研究. 高压物理学报, 1999, 13(2): 115-119 . doi: 10.11858/gywlxb.1999.02.006
    [15] 王同权张若棋汤文辉肖亚斌 . 强流质子束辐照形成的冲击波. 高压物理学报, 2003, 17(2): 117-121 . doi: 10.11858/gywlxb.2003.02.007
    [16] 温殿英林其文 . 冲击波压缩PVDF膜的电响应研究. 高压物理学报, 2000, 14(4): 291-297 . doi: 10.11858/gywlxb.2000.04.010
    [17] 章冠人 . 冲击波和分形集团的凝聚. 高压物理学报, 1997, 11(4): 241-244 . doi: 10.11858/gywlxb.1997.04.001
    [18] 谭华 . 金属的冲击波温度测量(Ⅱ)──界面卸载近似. 高压物理学报, 1996, 10(3): 161-169 . doi: 10.11858/gywlxb.1996.03.001
    [19] 谭华 . 金属的冲击波温度测量(Ⅲ)基板/样品界面间隙对辐射法测量冲击波温度的影响. 高压物理学报, 1999, 13(3): 161-168 . doi: 10.11858/gywlxb.1999.03.001
    [20] 肖强强黄正祥顾晓辉 . 冲击波影响下的聚能射流侵彻扩孔方程. 高压物理学报, 2011, 25(4): 333-338 . doi: 10.11858/gywlxb.2011.04.008
  • 加载中
图(13)表(4)
计量
  • 文章访问数:  493
  • 阅读全文浏览量:  500
  • PDF下载量:  8
出版历程
  • 收稿日期:  2018-07-23
  • 录用日期:  2018-08-23
  • 网络出版日期:  2019-01-17
  • 刊出日期:  2019-02-01

近场爆炸冲击波对屏蔽压装TNT的冲击引爆试验和仿真

    作者简介:陈兴旺(1994-),男,博士研究生,主要从事屏蔽装药结构安全性研究. E-mail: m18851199809@163.com
    通讯作者: 王金相, wjx@njust.edu.cn
  • 1. 南京理工大学瞬态物理国家重点实验室,江苏 南京 210094
  • 2. 中国北方车辆研究所,北京 100072
  • 3. 晋西工业集团,山西 太原 030027

摘要: 为研究近场强冲击波对屏蔽装药的冲击引爆效应,通过试验和LS-DYNA3D非线性有限元计算程序对屏蔽压装TNT的冲击引爆进行研究,得到了接触爆炸时压装TNT被引爆的临界屏蔽板厚度和非接触爆炸时屏蔽压装TNT的殉爆距离,分析了非接触爆炸时屏蔽板厚度对殉爆距离的影响,并通过非线性最小二乘法拟合得到屏蔽板厚度与殉爆距离的函数关系。结果表明:数值计算结果与试验结果基本一致;接触爆炸时临界起爆的屏蔽板厚度在20~23 mm之间;屏蔽板采用3 mm厚的45钢时,非接触爆炸作用下压装TNT的殉爆距离在12~15 mm之间;非接触爆炸时殉爆距离随着屏蔽板厚度的增加而减小,当无屏蔽板时,压装TNT的殉爆距离为79 mm,当屏蔽板厚度从1 mm加大到9 mm时,殉爆距离从51 mm减为1.5 mm,被发装药的屏蔽板对冲击波有显著的防护作用。

English Abstract

  • 屏蔽装药结构在外界强载荷作用下的安全性对于武器弹药研制具有重要意义。冲击波和高速破片是战场环境中两种最常见的威胁,研究含装药结构在上述外界能量作用下的冲击引爆机理和毁伤效应可为含装药弹药的安全性设计提供依据[1-2]。目前很多学者对破片或破片群对屏蔽装药结构的冲击引爆开展了研究[3-6]。对于炸药的冲击引爆,Walker等[7]以热起爆理论为基础,推导出著名的非均相炸药的起爆判据p2t=C(其中p为炸药界面的冲击压力,t为压力持续作用时间,C是与炸药有关的参数),此判据为不敏感弹药的冲击引爆研究提供了理论基础。Lu等[8]对裸装PBXN-109炸药殉爆试验进行了数值仿真,研究了爆炸冲击波对裸露炸药的冲击引爆作用。王晨等[9]对壳装固黑铝炸药殉爆试验进行了数值仿真,分析了主发装药壳体厚度对殉爆距离的影响。然而,作为弹药不敏感研究的重要方面,屏蔽装药结构在超近场强冲击波作用下的响应研究却鲜见报道,因此有必要开展爆炸冲击波对屏蔽装药的冲击引爆研究,从而为弹药安全性设计提供指导。

    本研究以压装TNT为例,屏蔽板采用45钢,通过试验和数值模拟方法开展接触和非接触爆炸两种情况下压装TNT在近场爆炸冲击波作用下的冲击引爆效应研究,重点分析屏蔽板厚度和起爆距离对冲击引爆的影响,以期确定临界起爆条件。

    • 采用非线性有限元计算程序LS-DYNA对屏蔽装药的冲击引爆进行数值模拟。如图1所示,有限元计算模型包括:主发装药、屏蔽板、被发装药、见证板和空气域,具体模型尺寸与试验对应。为简化计算模型,不考虑雷管、PPR套筒和传爆药。由于试验装置为轴对称结构,为减少计算量,采用1/4模型进行计算。数值计算采用ALE(Arbitrary Lagrange-Euler)算法,主发装药和空气域采用Euler网格,屏蔽板和被发装药采用Lagrange网格,屏蔽板与被发装药以及被发装药与见证板之间的接触均采用滑移接触,以便观察被发装药的起爆过程和见证板的变形。起爆点设置在主发装药表面中心点处。

      图  1  屏蔽装药冲击引爆数值计算模型图

      Figure 1.  Numerical calculation model of covered charge explosive detonation

    • 空气域采用无偏应力动力模型(Null)和Linear Polynomial状态方程[10]。主发装药采用HIGH_EXPLOSIVE_BURN本构模型和JWL状态方程,JWL状态方程如下

      式中:p为反应物的压力,V为反应物的比容,E为单位体积炸药内能,ωABR1R2为表征炸药特性的常数。屏蔽板和见证板材料模型均采用Johnson-Cook模型和Grüneisen状态方程,屏蔽板和主发装药的具体材料参数如表1所示,其中:ρ0为密度;D为爆速;pCJ为CJ压力;G为剪切模量;A1B1C1mn为材料常数,分别为材料的初始屈服应力、应变强化系数、应变速率强化系数、温度软化系数、应变强化指数;Tr为常温,Tm为材料的熔化温度。

      Material ρ0/(g·cm-3 D/(m·s-1 pCJ/GPa A/GPa B/GPa R1 R2 ω E/(109 J·m-3)
      TNT 1.60 6 812 18.56 370 3.23 4.15 0.95 0.3 7
      Material ρ0/(g·cm-3 G/GPa A1/GPa B1/GPa C1 n m Tm/K Tr/K
      45 steel 7.85 82.3 0.507 0.320 0.280 0.064 1.06 1 765 298
      Q235 steel 7.83 80.0 0.235 0.250 0.022 1.000 0.36 1 793 294

      表 1  主发炸药和屏蔽板的材料模型参数[11-13]

      Table 1.  Parameters for donor charge and shell model[11-13]

      被发装药采用点火增长反应模型,该模型可以反映受冲击炸药内部的爆轰发展过程,具体形式为

      式中:λ为反应度,即已反应炸药量与炸药总量的比值,它在模拟爆轰过程中控制着炸药化学能的释放;IG1G2abcdegxyz为材料参数,具体取值如表2所示。(2)式等号右侧第1项代表部分炸药在冲击压缩下被点火,第2项代表热点的增长,第3项代表主要反应后相对缓慢的扩散控制反应。

      I/μs–1 G1/( Pa–3·s–1) G2/( Pa–3·s–1) a b c d e g x y z
      8×108 4.20×10–25 2.60×10–25 0.111 0.667 0.667 0.667 0.333 1 6 3 3

      表 2  TNT点火增长反应模型参数[13]

      Table 2.  Ignition and growth reactive flow model parameters for TNT[13]

    • 本工作主要研究近场强冲击波对屏蔽装药的冲击引爆效应,并获得接触爆炸时压装TNT的临界屏蔽板厚度和非接触爆炸时屏蔽压装TNT的殉爆距离。为避免大量试验,依据仿真结果和初步理论计算[14]设计了3发接触爆炸试验和4发非接触爆炸试验。

    • 图2为接触爆炸屏蔽装药冲击引爆实物照片。共进行7发屏蔽装药冲击引爆试验。其中,接触爆炸试验3发:试验装置由雷管、传爆药、主发装药、45钢屏蔽板、被发装药和见证板组成;主发装药和被发装药均为柱形压装TNT,且尺寸相同,直径为35 mm,高度为30 mm;传爆药为柱形聚黑-14,其密度为1.136 g/cm3,直径为20 mm,厚度为5 mm,放置在主发装药的正上方;主发装药与45钢屏蔽板直接接触,屏蔽板厚度分别为20、23、26 mm;见证板为Q235钢,厚度为20 mm。非接触爆炸试验4发:试验装置比接触爆炸试验装置多一个PPR套筒,用来隔离主发装药与45钢屏蔽板;PPR套筒内径33.7 mm,外径40.0 mm,高度分别为10、12、15、19 mm;屏蔽板厚度均为3 mm。

      图  2  屏蔽装药冲击引爆实物照片

      Figure 2.  Photo of shock initiation of covered charge

      接触爆炸试验中,首先由雷管起爆传爆药,进而起爆主发装药,主发装药爆炸产生的冲击波直接经45钢屏蔽板衰减后作用于被发装药;在非接触屏蔽装药冲击引爆试验中,主发装药爆炸产生的冲击波分别经过空气和45钢屏蔽板衰减后作用于压装TNT。

      通过改变屏蔽板厚度d和PPR套筒高度h,使冲击波以不同幅值冲击引爆被发装药;根据见证板的变形情况和高速摄像的记录过程,判断被发装药的引爆情况。

    • 以非接触爆炸为例,通过高速摄影记录被发装药的起爆和熄爆过程。图3图4分别是套筒高度h为10和15 mm时不同时刻起爆装置的高速摄影照片。对比图3图4发现,主发装药TNT爆炸后,图3(a)图4(a)中装置上端均产生亮光。图3中,被发装药爆炸后,屏蔽板下面的被发装药有新的亮光产生,且爆炸后产生的冲击波和爆轰产物将地面黑色尘土颗粒掀起。图4中,被发装药未爆炸,屏蔽板下面的被发炸药没有产生新的亮光,爆轰剧烈程度相对较小。

      图  3  套筒高度为10 mm时不同时刻起爆装置的高速摄影照片

      Figure 3.  High speed camera photos of experimental setup with 10 mm PPR sleeve at different time

      图  4  套筒高度为15 mm时不同时刻起爆装置的高速摄影照片

      Figure 4.  High speed camera photos of experimental setup with 15 mm PPR sleeve at different time

      图5图6分别显示了接触和非接触爆炸后被发装药的见证板变形情况。对比图5图6发现:当被发装药爆炸后,试验后无剩余的TNT炸药粉末,见证板被炸出凹坑,坑深为4.2 mm,凹坑直径为40 mm,与炸药直径基本相同;当被发装药未起爆,试验后见证板上遗留大量的白色粉末,没有发生任何变形。

      图  5  接触爆炸试验后被发装药的见证板变形

      Figure 5.  Deformation of the witness plates under the acceptor with contact explosion

      图  6  非接触爆炸试验后被发装药的见证板变形

      Figure 6.  Deformation of the witness plates under the acceptor with non-contact explosion

      表3列出了屏蔽装药冲击引爆试验结果。接触爆炸试验中:当屏蔽板厚度d=20 mm时,被发炸药完全爆轰;当屏蔽板厚度d为23和26 mm时,炸药未爆炸。因此,压装TNT的临界屏蔽板厚度在20~23 mm之间。非接触爆炸试验中,PPR套筒高度为主发装药和45钢屏蔽板之间的垂直距离。当套筒高度h=10 mm时,被发炸药完全爆炸;h=12 mm时,见证板的坑深较浅,且残留炸药粉末,被发炸药仅部分爆炸;h=15 mm和h=19 mm时,炸药未爆炸。由此得到,屏蔽板厚度d=3 mm的压装TNT装药的临界起爆距离在12~15 mm之间。

      Explosion type Test No. d/mm h/mm Explosion situation of acceptor charge Pit depth/mm
      Contact explosion 1 20 Exploded 4.22
      2 26 Unexploded 0
      3 23 Unexploded 0
      Non-contact explosion 4 3 10 Exploded 4.20
      5 3 19 Unexploded 0
      6 3 15 Unexploded 0
      7 3 12 Part of the explosion 2.20

      表 3  屏蔽装药冲击引爆试验结果

      Table 3.  Results of shock initiation of covered charge

    • 为了验证数值模拟的有效性,对临界屏蔽板厚度、临界起爆距离、见证板变形情况的试验结果和数值模拟结果进行对比,如表4所示。可见,在接触爆炸试验中,压装TNT临界起爆的屏蔽板厚度在20~23 mm之间,仿真计算结果为24.0 mm,数值模拟结果与试验结果的相对误差不超过16.7%。当屏蔽板采用3 mm厚的45钢时,在非接触爆炸作用下压装TNT殉爆距离的试验结果在12~15 mm之间,数值计算结果为13.0 mm,相对误差不超过15.4%。数值计算结果与试验结果基本一致。

      Explosion type Test results Simulation results
      d/mm Explosion situation d/mm Explosion situation
      Contact explosion 20 Exploded 24.0 Exploded
      23 Unexploded 24.5 Unexploded
      Explosion type Test results Simulation results
      h/mm Explosion situation h/mm Explosion situation
      Non-contact explosion 12 Partial explosion 13.0 Exploded
      15 Unexploded 13.5 Unexploded

      表 4  屏蔽装药冲击引爆试验和仿真计算结果

      Table 4.  Test and simulation results of shock initiation of covered charge

      图7显示了h=10 mm时被发装药爆炸后数值模拟的见证板变形情况。被发装药爆炸后,数值计算得到的Q235钢见证板的坑深为3.80 mm,凹坑直径为42.0 mm,而试验所得坑深为4.20 mm,直径为40.0 mm。数值模拟结果与试验结果的相对误差不超过9.52%,二者基本一致。

      图  7  被发装药爆炸后见证板变形的数值模拟结果

      Figure 7.  Numerical simulation results of witness plate deformation

    • 为研究被发装药的冲击起爆过程,以非接触爆炸为例分析被发装药的爆轰成长过程和熄爆过程。图8为数值计算得到的h=13.0 mm时炸药不同时刻的变形图。由图8可见:t=8 μs时,主发炸药爆炸产生的冲击波刚好传到被发炸药上表面;t=11 μs时,被发装药爆炸,爆轰波由上向下传播;16 μs以后,炸药完全起爆,见证板有明显的凹坑。图9显示了h=13.5 mm时炸药不同时刻的变形情况。由图9可见:t=12 μs时,冲击波已传入被发装药一段时间,但是由于冲击波压力较小,并未引爆被发装药,冲击波继续由上向下传播;当t=16 μs以后,冲击波从非反射界面流出,被发装药仍未爆炸,只是由于爆轰产物膨胀和屏蔽板压缩发生了一定的变形,且见证板也未出现任何变化。

      图  8  13.0 mm距离下炸药不同时刻的变形

      Figure 8.  Calculated pressure contours at various time for the configuration with 13 mm separation

      图  9  13.5 mm距离下炸药不同时刻的变形

      Figure 9.  Calculated pressure contours at various times for the configuration with 13.5 mm separation

      图10图11显示h分别为13.0和13.5 mm时被发装药的爆轰成长过程和熄爆过程,图中AK为被发装药的不同位置。对比图10图11发现:h=13.0 mm时,被发装药表面所受冲击波压力达到2.1 GPa后,炸药开始反应,压力瞬间升高,最高压力达到9.8 GPa左右;随着爆轰的成长,在距离炸药表面22 mm处,最大压力达到18 GPa,并在以后的距离下基本保持稳定,达到稳定爆轰状态。当h=13.5 mm时,被发装药在受到峰值压力为1.8 GPa的冲击波作用时,未发生任何反应,在侧向和轴向稀疏波的共同作用下,随着传播深度的增加,压力逐渐衰减,被发装药未爆炸。

      图  10  13.0 mm间隔高度下被发装药的爆轰成长过程

      Figure 10.  Detonation development process with spacing distance of 13.0 mm

      图  11  13.5 mm间隔高度下被发装药的熄爆过程

      Figure 11.  Extinguished process with spacing distance of 13.5 mm

      对于不敏感炸药TNT,主要是由于爆轰产物的冲击使被发装药产生热点,从而导致爆炸或爆轰。因此,采用炸药临界起爆压力判据,得到压装TNT炸药的临界起爆压力约为2.1 GPa,爆轰成长距离约为22 mm。

    • 屏蔽板厚度对被发装药的冲击起爆起着关键作用。图12为屏蔽装药冲击引爆示意图。主发炸药爆炸后在空气中产生冲击波,冲击波经过一段距离的传播到达屏蔽板,在空气与屏蔽板的接触面上发生反射和透射,由于空气的波阻抗小于屏蔽板的波阻抗,因此在空气中反射冲击波,并在屏蔽板内透射冲击波。当屏蔽板中的冲击波传播到屏蔽板与被发炸药界面处时,由于被发炸药的波阻抗远小于屏蔽板的波阻抗,在屏蔽板与被发炸药界面处冲击波产生卸载反射,向屏蔽板反射稀疏波,而向被发炸药透射冲击波。对于远场爆炸,通过空气冲击波峰值超压的经验公式,求得冲击波到达屏蔽板界面的冲击波参数;再利用冲击波前后质量守恒、动量守恒和连续条件,可求得进入被发装药的冲击波压力和质点速度;最后,根据炸药临界起爆压力判据,判断被发装药是否起爆[14]。对于近场爆炸,由于主发装药与屏蔽板的距离太小,爆轰产物的影响不可忽略,目前尚无成熟理论和经验公式对冲击波衰减后到达屏蔽板界面的冲击波参数进行求解,因此无法通过理论计算得到进入炸药的冲击波压力以及屏蔽板厚度与殉爆距离之间的关系。为了分析近场非接触爆炸时屏蔽板厚度对殉爆距离的影响,本研究结合试验数据,通过非线性最小二乘法拟合得到二者的关系,具体形式如下

      式中:d为屏蔽板厚度(mm),h为殉爆距离(mm),abmnC为表征冲击引爆的特性参数,通过拟合图12中的数据得到各参数值:a=4.45,b=1.83,m=3,n=1,C=8 033。由(3)式可知,随着屏蔽板厚度的增大,殉爆距离呈幂指数衰减,拟合得到的屏蔽板厚度与殉爆距离的函数关系可以很好地描述屏蔽板的防护特性以及殉爆距离对被发装药冲击引爆的影响。如图13所示,当无屏蔽板时,压装TNT的殉爆距离为79 mm;当屏蔽板厚度从1 mm逐渐加大到9 mm时,屏蔽压装TNT的殉爆距离从51 mm降为1.5 m,屏蔽板对冲击波有显著的防护作用;屏蔽板厚度在3 mm以内时,厚度加大对冲击引爆效应的影响较大,大于3 mm时,进一步增大屏蔽板厚度对提高炸药抗冲击引爆能力的影响显著变小。

      图  13  殉爆距离随屏蔽板厚度变化曲线

      Figure 13.  Variation of sympathetic detonation distance with shield thickness

      图  12  屏蔽装药冲击起爆示意

      Figure 12.  Shock initiation of covered charge

    • (1)对于装药尺寸为ø35 mm×30 mm的压装TNT,当屏蔽板为45钢时:在接触爆炸作用下,试验和数值模拟得到的临界起爆的屏蔽板厚度为20~23 mm和24 mm,相对误差小于16.67%;在非接触爆炸作用下,实验得到的殉爆距离在12~15 mm之间,数值模拟得到的临界殉爆距离为13 mm,与试验结果的相对误差小于15.4%。数值计算结果与试验结果基本一致,表明采用ALE算法的计算模型能够有效描述屏蔽装药冲击引爆试验。

      (2)本试验中,压装TNT炸药的临界起爆压力约为2.1 GPa,爆轰成长距离约为22 mm。

      (3)非接触爆炸条件下:殉爆距离随着屏蔽板厚度的增加而减小;当无屏蔽板时,压装TNT的殉爆距离为79 mm;当屏蔽板厚度从1 mm逐渐加大到9 mm时,屏蔽压装TNT的殉爆距离从51 mm降为1.5 mm,屏蔽板对冲击波有显著的防护作用。

      (4)通过非线性最小二乘法拟合得到屏蔽板厚度与殉爆距离的函数关系,公式可以在一定范围内很好地描述屏蔽板的防护特性以及殉爆距离对被发装药冲击引爆的影响,为屏蔽装药的冲击引爆试验研究和后续的可靠性研究提供一定的理论参考。

参考文献 (14)

目录

    /

    返回文章
    返回