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

CL-20基混合炸药的冲击起爆特征

皮铮迪 陈朗 刘丹阳 伍俊英

引用本文:
Citation:

CL-20基混合炸药的冲击起爆特征

    作者简介: 皮铮迪(1984—),男,博士;
    通讯作者: 陈朗, chenlang@bit.edu.cn
  • 中图分类号: O381

Shock initiation of CL-20 based explosives

    Corresponding author: Chen Lang, chenlang@bit.edu.cn ;
  • CLC number: O381

  • 摘要: 为了研究CL-20基混合炸药的冲击起爆特征,深入分析冲击波作用下CL-20基混合炸药的爆轰成长规律,采用炸药驱动飞片冲击起爆实验方法,对CL-20、CL-20/NTO和CL-20/FOX-7三种压装混合炸药进行了冲击起爆实验,通过嵌入在炸药内部不同位置处的锰铜压力传感器,获得了炸药内部压力的变化历程。依据实验结果标定了混合炸药的点火增长模型参数,其中,利用反应速率方程中的两个增长项,分别模拟CL-20/NTO和CL-20/FOX-7混合炸药中两种组分的反应增长过程,得到这两种混合炸药的反应速率方程参数。并通过数值模拟的方法得到了三种炸药的临界起爆阈值和POP关系。研究结果表明:三种CL-20基混合炸药中,CL-20/NTO混合炸药具有更高的临界起爆阈值,而CL-20/FOX-7混合炸药具有更长的爆轰成长距离;此外,利用此套拟合双组分混合炸药反应速率方程的方法,可以对新型配方炸药的冲击起爆过程进行预测性计算
  • 图 1  炸药驱动飞片冲击起爆实验装置示意图及照片

    Figure 1.  Diagram and photo of shock initiation experimental apparatus

    图 2  CL-20混合炸药在距冲击波入射面不同距离处的内部压力-时间曲线

    Figure 2.  Pressure curves of CL-20 explosive at different displacements from shock incidence plane

    图 3  CL-20/NTO混合炸药在距冲击波入射面不同距离处的内部压力-时间曲线

    Figure 3.  Pressure curves of CL-20/NTO explosive at different displacements from shock incidence plane

    图 4  CL-20/FOX-7混合炸药在距冲击波入射面不同距离处的内部压力-时间曲线

    Figure 4.  Pressure curves of CL-20/FOX-7 explosive at different displacements from shock incidence plane

    图 5  CL-20混合炸药计算与实验结果对比

    Figure 5.  Experimental and calculated pressure histories in shock initiation of CL-20 explosive

    图 6  CL-20/NTO混合炸药计算与实验结果对比

    Figure 6.  Experimental and calculated pressure histories in shock initiation of CL-20/NTO

    图 7  CL-20/FOX-7混合炸药计算与实验结果对比

    Figure 7.  Experimental and calculated pressure histories in shock initiation of CL-20/FOX-7

    图 8  三种炸药的POP关系图

    Figure 8.  Pop plots for the three explosives

    表 1  实验炸药配方及平均密度

    Table 1.  Formula and average density of the samples

    编号 炸药配方(质量比) 平均密度/(g·cm-3)
    1 CL-20 (95%CL-20, 5%Binder) 1.94
    2 CL-20/NTO (47%CL-20, 47%NTO, 6%FPM) 1.89
    3 CL-20/FOX-7 (47%CL-20, 47%FOX-7, 6%FPM) 1.88
    下载: 导出CSV

    表 2  8701和TNT的材料参数[7]

    Table 2.  Parameters of TNT and 8701 explosive

    炸药 ρ/(g·cm-3) D/(cm·μs-1) PCJ/GPa A/GPa B/GPa R1 R2 ω E0/GPa
    8701 1.70 0.831 5 29.50 854.5 20.49 4.60 1.35 0.25 8.5
    TNT 1.64 0.693 0 27.00 371.3 3.23 4.15 0.95 0.30 7.0
    下载: 导出CSV

    表 3  几种惰性材料的Grüneisen状态方程参数[7]

    Table 3.  Grüneisen parameters for inert materials

    惰性材料 ρ0/(g·cm-3) c/(mm·μs-1) S1 S2 S3 γ0 a
    PMMA 1.851 2.240 2.090 -1.12 0.0 0.850 0.00
    Steel 7.830 4.570 1.490 0.00 0.0 1.930 0.50
    Al 2.700 5.355 1.345 0.00 0.0 2.130 0.48
    Polyimide 1.414 2.470 1.629 0.00 0.0 1.395 0.00
    下载: 导出CSV

    表 4  CL-20、CL-20/NTO、CL-20/FOX-7混合炸药的点火增长模型参数

    Table 4.  Ignition and growth model parameters for CL-20, CL-20/NTO and CL-20/FOX-7

    参数 未反应JWL状态方程
    CL-20 CL-20/NTO CL-20/FOX-7
    ρ0/(g·cm-3) 1.94 1.89 1.88
    A/GPa 444 400 24 900 122 300
    B/GPa -5.13 -3.18 -6.91
    R1 13.5 11.3 11.3
    R2 1.35 1.13 1.13
    ω 0.869 5 0.869 5 0.869 5
    Cv/(GPa·K-1) 2.78×10-3 2.78×10-3 2.78×10-3
    T0/K 298 298 298
    剪切模量/GPa 4.54 4.00 4.00
    屈服强度/GPa 0.2 0.2 0.2
    参数 爆轰产物JWL状态方程
    CL-20 CL-20/NTO CL-20/FOX-7
    A/GPa 1 887.64 1 640.27 1 589.47
    B/GPa 162.40 22.26 18.34
    R1 6.50 5.93 5.78
    R2 2.75 1.76 1.54
    ω 0.547 0.700 0.700
    Cv/(GPa·K-1) 1.0×10-3 1.0×10-3 1.0×10-3
    E0/GPa 11.50 8.46 9.07
    参数 点火增长反应速率方程
    CL-20 CL-20/NTO CL-20/FOX-7
    I 7.43×1011 4.0×106 7.43×1011
    a 0.0 0.22 0.0
    b 0.667 0.667 0.667
    x 20.0 7.0 20.0
    G1 402 402 254
    c 0.667 0.667 0.667
    d 0.333 0.333 0.333
    y 2.0 2.0 2.0
    G2 400 0.6 140
    e 0.333 0.667 0.667
    g 1.0 0.111 0.333
    z 2.0 1.0 2.0
    Figmax 0.3 0.3 0.3
    FG1max 0.5 0.5 0.5
    FG2min 0.5 0.0 0.0
    下载: 导出CSV

    表 5  CL-20/NTO和CL-20/FOX-7两种混合炸药的配方和生成热数据

    Table 5.  Formula and heat of formation for the explosive of CL-20/NTO and CL-20/FOX-7

    炸药配方(质量比) 生成热/(J·g-1) 物质的量(分子量为100时)/mol
    C H N O F
    CL-20:NTO:FPM(47:47:6) -189.86 1.537 5 1.507 4 2.733 9 2.372 3 0.200 7
    CL-20:FOX-7:FPM(47:47:6) -141.53 1.449 5 2.054 6 2.558 0 2.558 0 0.200 7
    下载: 导出CSV

    表 6  三种混合炸药冲击起爆的临界阈值

    Table 6.  Shock initiation critical thresholds for three mixed explosives

    炸药(质量比) 密度/(g·cm-3) 临界压力/GPa p2τ/(Pa2·s) puτ/(J·m-2)
    CL-20(95/5) 1.94 1.39 253.67×1010 41.97×104
    CL-20/NTO/FPM(47/47/6) 1.89 1.46 279.86×1010 66.71×104
    CL-20/FOX-7/FPM(47/47/6) 1.88 1.40 257.33×1010 43.19×104
    下载: 导出CSV
  • [1] Nielsen A T, Nissan R A.Polynitropolyaza caged explosives[J].Naval Weapon Center Technical Publication, 1986(5):6692.
    [2] Simpson R L, Urtiew P A, Ornellas D L, et al.CL-20 performance exceeds that of HMX and its sensitivity is moderate[J].Propellants, Explosive, Pyrotechnics, 1997, 22(5):249-255. doi: 10.1002/(ISSN)1521-4087
    [3] Tarver C M, Simposon R L, Urtiew P A. Shock initiation of an ε-CL-20-estane formulation[C]//Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Metter. AIP Publishing LLC, 1996: 891-894.
    [4] Lee E L, Tarver C M. Phenomenological model of shock initiation in heterogeneous explosives[C]//Physics of Fluids. AIP Conference Proceedings. AIP Publishing LLC, 1980: 2262-2372.
    [5] Hallquist J O.LS-DYNA user's manual:nonlinear dynamic analysis of structures in three dimensions[M].California:University of California, 2001:823-826.
    [6] Lee E, Finger M, Collins W. JWL equation of state coefficients for high explosives: UCLD-16189[R], 1973.
    [7] 陈朗, 刘群, 伍俊英.受热炸药的冲击起爆特征[J].爆炸与冲击, 2013, 33(1):21-28. doi: 10.3969/j.issn.1001-1455.2013.01.003
    Chen Lang, Liu Qun, Wu Junying.On shock initiation of heated explosives[J].Explosion and Shock Waves, 2013, 33(1):21-28. doi: 10.3969/j.issn.1001-1455.2013.01.003
    [8] Daniel J S. An equation of state for ploy-methyl-meth acrylate: UCID-16982[R], 1975.
    [9] Mader C L.Numerical modeling of detonation[M].Berkeley:California Press, 1979.
    [10] Mader C L. Detonation properties of condensed explosives computed using the BKW equation of state: LA-2900[R]. Los Alamos Scientific Laboratory Report, 1963.
    [11] Tarver C M, Hallquist J O, Erickson L M. Modeling short pulse duration shock initiation of solid explosives[C]//Proceedings of the 8th International Symposium of Detonation. Albuquerque, USA, 1985: 951-961.
    [12] Urtiew P A, Vandersall K S, Tarver C M, et al.Shock initiation of composition B and C-4 explosives:experiments and modeling[J].Russian Journal of Physical Chemistry:B, 2008, 2:162-171. doi: 10.1134/S1990793108020036
    [13] Walker F E, Wasley R J.Critical energy for shock initiation of heterogeneous explosives[J].Explosives Stoffe, 1969, 17(1):9-13.
    [14] Ramsay J B, Popolato A. Analysis of shock wave and initiation data for solid explosives: LA-DC-6992;CONF-651003-3[R]. Los Alamos Scientific Lab., Univ. of California, N. Mex., 1965.
  • [1] 张力陈朗王晨伍俊英 . 不同晶型CL-20热分解反应机理计算分析. 爆炸与冲击, 2014, 34(2): 188-195. doi: 10.11883/1001-1455(2014)02-0188-07
    [2] 陈朗刘群伍俊英 . 受热炸药的冲击起爆特征. 爆炸与冲击, 2013, 33(1): 21-28. doi: 10.11883/1001-1455(2013)01-0021-08
    [3] 金朋刚常海陈智群 . FOX-7热分解动力学和机理研究. 爆炸与冲击, 2006, 26(6): 528-531. doi: 10.11883/1001-1455(2006)06-0528-04
    [4] 冯凇饶国宁彭金华汪斌 . CL-20基炸药水中爆炸气泡脉动实验研究. 爆炸与冲击, 2018, 38(4): 855-862. doi: 10.11883/bzycj-2017-0093
    [5] 梁增友黄风雷张震宇 . PBX-9404的化学反应速率方程及起爆特性. 爆炸与冲击, 2008, 28(1): 38-41. doi: 1.011883/1001-1455(2008)01-0038-06
    [6] 王桂吉赵同虎莫建军吴刚韩梅谭福利 . 一种以TATB/HMX为基的高聚物粘结炸药的短脉冲冲击起爆特性. 爆炸与冲击, 2007, 27(3): 230-235. doi: 10.11883/1001-1455(2007)03-0230-06
    [7] 陈少杰吴立志沈瑞琪叶迎华胡艳 . 激光驱动复合飞片冲击起爆HNS-Ⅳ. 爆炸与冲击, 2015, 35(2): 285-288. doi: 10.11883/1001-1455-(2015)02-0285-04
    [8] 张忠陈卫东杨文淼 . 非均质固体炸药冲击起爆的物质点法. 爆炸与冲击, 2011, 31(1): 25-30. doi: 10.11883/1001-1455(2011)01-0025-06
    [9] 陶为俊浣石黄风雷蒋国平 . 侧向稀疏波对非均质凝聚炸药冲击波起爆过程的影响. 爆炸与冲击, 2011, 31(4): 397-401. doi: 10.11883/1001-1455(2011)04-0397-05
    [10] 张涛刘雨生高志鹏杨佳刘艺谷岩 . 垫层对破片冲击起爆带壳炸药影响的数值模拟. 爆炸与冲击, 2018, 38(6): 1241-1246. doi: 10.11883/bzycj-2017-0154
    [11] 潘昊胡晓棉 . 钝感炸药的超压爆轰与冲击起爆过程数值模拟. 爆炸与冲击, 2006, 26(2): 174-178. doi: 10.11883/1001-1455(2006)02-0174-05
    [12] 李志鹏龙新平黄毅民何碧王蓉何松伟 . 用组合式电磁粒子速度计研究JOB-9003炸药的冲击起爆过程. 爆炸与冲击, 2006, 26(3): 269-272. doi: 10.11883/1001-1455(2006)03-0269-04
    [13] 王昕蒋建伟王树有门建兵 . 破片撞击起爆柱面带壳装药的临界速度修正判据. 爆炸与冲击, 2019, 39(1): 012302-1-012302-8. doi: 10.11883/bzycj-2017-0271
    [14] 胡宏伟严家佳陈朗郭炜宋浦 . 铝粉含量和粒度对CL-20含铝炸药水中爆炸反应特性的影响. 爆炸与冲击, 2017, 37(1): 157-161. doi: 10.11883/1001-1455(2017)01-0157-05
    [15] 王保国陈亚芳张景林高敏 . 亚微米级型CL-20的制备、表征与性能. 爆炸与冲击, 2009, 29(5): 550-554. doi: 10.11883/1001-1455(2009)05-0550-05
    [16] 董润安肖连阶姜德珍 . 高聚物粘结NTO的研究. 爆炸与冲击, 1995, 15(2): 116-121.
    [17] 张振宇卢芳云王志兵浣石 . PBX-9404炸药高压反应速率方程的研究. 爆炸与冲击, 1999, 19(4): 360-364.
    [18] 张琪敏张旭赵康舒俊翔张蓉钟斌 . TATB基钝感炸药JB-9014的冲击起爆反应增长规律. 爆炸与冲击, 2019, 39(4): 041405-1-041405-6. doi: 10.11883/bzycj-2018-0050
    [19] 宋江杰张振宇谭晓莉林华令成丽蓉 . 固体非均质炸药冲击点火与起爆模型研究进展. 爆炸与冲击, 2012, 32(2): 121-128. doi: 10.11883/1001-1455(2012)02-0121-08
    [20] . 爆炸与冲击1987年(第7卷)作者索引. 爆炸与冲击, 1987, 7(4): 382-384.
  • 加载中
图(8)表(6)
计量
  • 文章访问数:  414
  • HTML全文浏览量:  45
  • PDF下载量:  346
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-04-26
  • 录用日期:  2016-07-26
  • 刊出日期:  2017-11-25

CL-20基混合炸药的冲击起爆特征

    作者简介:皮铮迪(1984—),男,博士
    通讯作者: 陈朗, chenlang@bit.edu.cn
  • 北京理工大学爆炸科学与技术国家重点实验室,北京 100081

摘要: 为了研究CL-20基混合炸药的冲击起爆特征,深入分析冲击波作用下CL-20基混合炸药的爆轰成长规律,采用炸药驱动飞片冲击起爆实验方法,对CL-20、CL-20/NTO和CL-20/FOX-7三种压装混合炸药进行了冲击起爆实验,通过嵌入在炸药内部不同位置处的锰铜压力传感器,获得了炸药内部压力的变化历程。依据实验结果标定了混合炸药的点火增长模型参数,其中,利用反应速率方程中的两个增长项,分别模拟CL-20/NTO和CL-20/FOX-7混合炸药中两种组分的反应增长过程,得到这两种混合炸药的反应速率方程参数。并通过数值模拟的方法得到了三种炸药的临界起爆阈值和POP关系。研究结果表明:三种CL-20基混合炸药中,CL-20/NTO混合炸药具有更高的临界起爆阈值,而CL-20/FOX-7混合炸药具有更长的爆轰成长距离;此外,利用此套拟合双组分混合炸药反应速率方程的方法,可以对新型配方炸药的冲击起爆过程进行预测性计算

English Abstract

  • 1987年, A.T.Nielsen等[1]首次合成六硝基六氮杂异伍兹烷(CL-20),由于其具有较高的密度和能量[2],得到了广泛的关注和研究。随着现代武器系统的发展,对炸药使用的安全性要求逐渐提高,其中,在冲击加载作用下炸药起爆并成长为爆轰的过程是了解炸药安全性能的核心问题,因此,研究CL-20基混合炸药的冲击起爆规律具有重要意义。实验中,通常利用Lagrange量计测量冲击波后反应流场的变化,获得冲击作用下炸药的爆轰成长过程,C.M.Tarver等[3]对LX-19炸药(95.2%wt.ε-CL-20, 4.8% Estane)进行了冲击起爆实验研究,通过嵌入在炸药内部的锰铜压力传感器获得了距离待测炸药起爆面不同位置处的压力变化历程,根据实验结果,分析了LX-19炸药的冲击波感度。由于炸药冲击起爆过程的瞬时性以及炸药爆炸的破坏性,在实验过程中得到的数据仍然比较有限,因此,可以借助数值模拟的方法深入了解炸药的起爆和爆轰成长过程,其中,炸药在受到冲击加载后的反应速率是数值模拟研究的核心问题。基于经验的反应速率方程,已经可以较好地描述这一过程,其中,应用比较广泛的是Lee和Tarver设计的点火增长反应速率方程[4],此方程基于爆轰波的ZND模型和热点理论假设,描述炸药受冲击后的点火-增长过程,通过匹配平面一维反应冲击波阵面后反应流场的实验结果,拟合反应速率方程中的待定参数。可是,该反应速率方程的应用范围比较有限,一旦炸药的配方发生变化,就需要基于新的实验结果重新调整反应速率方程中的参数。

    本文中采用炸药驱动飞片的实验方法,对CL-20、CL-20/FOX-7和CL-20/NTO三种CL-20基混合炸药进行冲击起爆实验研究,利用嵌入在炸药内部的锰铜压力传感器,测量距离起爆面不同位置处的压力变化历程,获得三种CL-20基混合炸药的爆轰成长过程。采用点火增长反应速率模型,对炸药的冲击起爆过程进行数值模拟,根据实验结果拟合三种炸药的点火增长反应速率方程参数,在拟合CL-20/FOX-7和CL-20/NTO两种混合炸药的模型参数时,利用反应速率方程中两个独立的增长项分别描述混合炸药中两种组分炸药的反应过程。建立飞片撞击炸药的计算模型,得到三种CL-20基混合炸药的冲击起爆阈值和入射冲击压力与爆轰成长距离关系(Pop plot),并对三种炸药的冲击感度进行定量分析。

    • 炸药驱动飞片冲击起爆实验装置如图 1所示,装置由雷管、平面波炸药透镜、加载炸药、有机玻璃隔板、钢飞片、有机玻璃支架、铝隔板、炸药试件和见证板组成。实验时,通过雷管起爆炸药平面波透镜和加载炸药,产生一束平面波,经有机玻璃隔板衰减后作用于钢飞片并驱动其向下运动,飞片撞击铝隔板,产生的平面冲击波经铝隔板衰减后最终起爆待测炸药。通过改变有机玻璃隔板和铝隔板的厚度来调节起爆待测炸药的冲击波强度。平面波炸药透镜和加载炸药的直径均为50 mm,有机玻璃隔板直径为60 mm,钢飞片的直径为50 mm、厚3 mm,铝隔板的直径为150 mm。

      图  1  炸药驱动飞片冲击起爆实验装置示意图及照片

      Figure 1.  Diagram and photo of shock initiation experimental apparatus

      实验炸药为CL-20、CL-20/NTO和CL-20/FOX-7压装混合炸药,详细配方如表 1所示,按照预定的压制工艺将其制成直径40 mm,高度分别为25 mm和1~4 mm的药柱,药柱间密度极差不超过0.02 g/cm3,三种炸药的平均密度列于表 1中。在装配炸药试件时,将高25mm的药柱作为底座置于最底端,在其上方依次放置三个厚度1~4 mm的药柱,在药柱间和试件的最上方嵌入锰铜压力传感器,其中,传感器被封装在两层0.1 mm厚的聚四氟乙烯薄膜之中,按此结构,就可获得从冲击波入射面至炸药内部0~8 mm之间4个位置处的压力变化历程。

      编号 炸药配方(质量比) 平均密度/(g·cm-3)
      1 CL-20 (95%CL-20, 5%Binder) 1.94
      2 CL-20/NTO (47%CL-20, 47%NTO, 6%FPM) 1.89
      3 CL-20/FOX-7 (47%CL-20, 47%FOX-7, 6%FPM) 1.88

      表 1  实验炸药配方及平均密度

      Table 1.  Formula and average density of the samples

      利用上述实验装置,对CL-20、CL-20/NTO和CL-20/FOX-7混合炸药进行了冲击起爆实验,每次实验获得冲击波入射面及炸药内部共4个位置处压力随时间的变化曲线,通过改变有机玻璃和铝隔板的厚度,调节入射冲击波压力。

    • 图 2是CL-20混合炸药在p0=4.18 GPa的入射压力下,距冲击波入射面不同距离x处的压力-时间变化曲线,从图中可以看出,当冲击波传至2 mm位置时,炸药反应释放的能量已能够支持到爆轰波阵面,但支持度很小,只是使波阵面压力有小幅度的提高,达到4.29 GPa,随着时间的推移,波后炸药反应释能的压力峰逐渐追赶前导冲击波,在4 mm处爆轰波阵面压力已达到25 GPa左右。

      图  2  CL-20混合炸药在距冲击波入射面不同距离处的内部压力-时间曲线

      Figure 2.  Pressure curves of CL-20 explosive at different displacements from shock incidence plane

      图 3是3.49 GPa的入射压力下,CL-20/NTO混合炸药内部不同位置的压力-时间曲线,可以看到,3.49 GPa压力的冲击波引发了炸药的反应,但反应相对缓慢,压力成长较慢,在5 mm之内炸药反应释放的能量没有支持到爆轰波阵面。另外,几乎每个前导波阵面后都出现了压力的卸载,即压力到达峰值后又迅速下降。这可能是多种机制引起的,例如混合压制炸药中存在的孔穴,在压力波的作用下闭合,出现压力卸载的现象;也可能是由于NTO的屈服强度相对较低,在冲击波的作用下发生了较大的形变,从而导致了压力的卸载。

      图  3  CL-20/NTO混合炸药在距冲击波入射面不同距离处的内部压力-时间曲线

      Figure 3.  Pressure curves of CL-20/NTO explosive at different displacements from shock incidence plane

      在入射压力为3.40 GPa时,CL-20/FOX-7混合炸药内部不同位置的压力变化历程如图 4所示,可以看出,在有效信号范围内,炸药释放的能量未能支持到爆轰波阵面,表明在此范围的入射压力下,波阵面只激发了少量的炸药反应,并且反应增长速率也相对缓慢。

      图  4  CL-20/FOX-7混合炸药在距冲击波入射面不同距离处的内部压力-时间曲线

      Figure 4.  Pressure curves of CL-20/FOX-7 explosive at different displacements from shock incidence plane

      CL-20、CL-20/NTO和CL-20/FOX-7混合炸药在相同冲击加载条件下的实验结果如图 2~4中所示,这三次实验中有机玻璃隔板厚度均为14 mm,铝隔板厚度均为10 mm,由于CL-20混合炸药的密度高于后两者,因此,其入射压力相对较高,为4.18 GPa,而后两种混合炸药的入射冲击波压力分别为3.49、3.40 GPa。通过比较三种混合炸药的压力成长趋势,可以明显看出,向CL-20中分别加入NTO和FOX-7,均可有效降低其冲击波感度。

    • 根据炸药驱动飞片冲击起爆实验装置,建立计算模型。由于装置具有对称性,因此采用二维轴对称模型,模型中各组件的尺寸与实验装置一致。在炸药内部嵌入了锰铜压力传感器,且传感器的上下表面均贴敷聚四氟乙烯薄膜,这些结构对冲击波在炸药内部的传播会产生一定的影响,因此在计算中也建立了传感器模型。采用非线性有限元计算方法[5]对炸药的起爆过程进行模拟。

      平面波炸药透镜由TNT和8701炸药组成,加载炸药为TNT,均采用高能炸药材料模型和JWL状态方程[6]进行描述。两种炸药的材料参数如表 2所示,表中ρ为炸药的初始密度,D为爆速,PCJ为C-J点爆轰压力。

      炸药 ρ/(g·cm-3) D/(cm·μs-1) PCJ/GPa A/GPa B/GPa R1 R2 ω E0/GPa
      8701 1.70 0.831 5 29.50 854.5 20.49 4.60 1.35 0.25 8.5
      TNT 1.64 0.693 0 27.00 371.3 3.23 4.15 0.95 0.30 7.0

      表 2  8701和TNT的材料参数[7]

      Table 2.  Parameters of TNT and 8701 explosive

      计算中有机玻璃隔板、支架、铝隔板、钢飞片、包裹传感器的聚四氟乙烯薄膜均为惰性材料,采用弹塑性流体力学材料模型和Grüneisen状态方程描述,Grüneisen状态方程[8]的具体形式为:

      $ P = \frac{{{\rho _0}{C^2}\mu \left[ {1 + \left( {1 - \frac{{{\gamma _0}}}{2}} \right)\mu - \frac{a}{2}{\mu ^2}} \right]}}{{{{\left[ {1 - \left( {{S_1} - 1} \right)\mu - {S_2}\frac{{{\mu ^2}}}{{\mu + 1}} - {S_3}\frac{{{\mu ^3}}}{{{{\left( {\mu + 1} \right)}^2}}}} \right]}^2}}} + \left( {{\lambda _0} + \mathit{a}\mu } \right)E $

      式中:μρ/ρ0-1,ρ是密度,ρ0是材料初始密度,Cus-up曲线截距。S1S2S3us-up曲线斜率的系数,γ0是Grüneisen系数,a是对γ0的一阶体积修正。表 3是几种惰性材料的Grüneisen状态方程参数。

      惰性材料 ρ0/(g·cm-3) c/(mm·μs-1) S1 S2 S3 γ0 a
      PMMA 1.851 2.240 2.090 -1.12 0.0 0.850 0.00
      Steel 7.830 4.570 1.490 0.00 0.0 1.930 0.50
      Al 2.700 5.355 1.345 0.00 0.0 2.130 0.48
      Polyimide 1.414 2.470 1.629 0.00 0.0 1.395 0.00

      表 3  几种惰性材料的Grüneisen状态方程参数[7]

      Table 3.  Grüneisen parameters for inert materials

      利用点火增长模型模拟凝聚炸药的冲击起爆过程,模型中使用JWL(Jones-Wikins-Lee)状态方程[9]来描述未反应炸药和爆轰产物,具体形式为:

      $ {P_{\rm{E}}} = A{{\rm{e}}^{ - {R_{{1^v}{\rm{E}}}}}} + B{{\rm{e}}^{ - {R_{{2^v}{\rm{E}}}}}} + \frac{{\omega {c_V}{T_0}}}{{{v_{\rm{E}}}}} $

      $ {P_{\rm{P}}} = A{{\rm{e}}^{ - {R_{{1^v}{\rm{P}}}}}} + B{{\rm{e}}^{ - {R_{{2^v}{\rm{P}}}}}} + \frac{{\omega {c_V}{T_{\rm{p}}}}}{{{v_{\rm{p}}}}} $

      式中:PEPp为炸药初始压力和产物压力,vEvp为炸药初始相对比容和产物相对比容,cV为比定压热容,T0Tp为炸药初始温度和产物温度,ABR1R2、和ω为待定参数。

      通过匹配冲击波传入炸药中最初几个位置的压力变化趋势,拟合出CL-20/NTO和CL-20/FOX-7两种混合炸药的未反应JWL状态方参数,由于CL-20混合炸药的压制密度与LX-19[3]接近,利用已有的LX-19未反应JWL状态方程参数进行试算,计算结果可以较好地模拟冲击波传入炸药中最初几个位置的压力波形,因此,CL-20混合炸药的未反应JWL状态方程选用已有的LX-19参数。三种混合炸药的未反应JWL状态方程参数如表 4中所示。利用CL-20混合炸药圆筒实验,得到密度为1.94g/cm3的CL-20混合炸药爆轰产物JWL状态方程,方程参数如表 4中所示。采用BKW状态方程[9]和计算程序[10]来计算CL-20/NTO和CL-20/FOX-7两种混合炸药爆轰产物的等熵膨胀,两种混合炸药的配方和生成热数据见表 5。利用BKW状态方程计算得到的炸药爆轰产物等熵膨胀数据,拟合出CL-20/NTO和CL-20/FOX-7两种混合炸药爆轰产物的JWL状态方程参数,详见表 4中所示。

      参数 未反应JWL状态方程
      CL-20 CL-20/NTO CL-20/FOX-7
      ρ0/(g·cm-3) 1.94 1.89 1.88
      A/GPa 444 400 24 900 122 300
      B/GPa -5.13 -3.18 -6.91
      R1 13.5 11.3 11.3
      R2 1.35 1.13 1.13
      ω 0.869 5 0.869 5 0.869 5
      Cv/(GPa·K-1) 2.78×10-3 2.78×10-3 2.78×10-3
      T0/K 298 298 298
      剪切模量/GPa 4.54 4.00 4.00
      屈服强度/GPa 0.2 0.2 0.2
      参数 爆轰产物JWL状态方程
      CL-20 CL-20/NTO CL-20/FOX-7
      A/GPa 1 887.64 1 640.27 1 589.47
      B/GPa 162.40 22.26 18.34
      R1 6.50 5.93 5.78
      R2 2.75 1.76 1.54
      ω 0.547 0.700 0.700
      Cv/(GPa·K-1) 1.0×10-3 1.0×10-3 1.0×10-3
      E0/GPa 11.50 8.46 9.07
      参数 点火增长反应速率方程
      CL-20 CL-20/NTO CL-20/FOX-7
      I 7.43×1011 4.0×106 7.43×1011
      a 0.0 0.22 0.0
      b 0.667 0.667 0.667
      x 20.0 7.0 20.0
      G1 402 402 254
      c 0.667 0.667 0.667
      d 0.333 0.333 0.333
      y 2.0 2.0 2.0
      G2 400 0.6 140
      e 0.333 0.667 0.667
      g 1.0 0.111 0.333
      z 2.0 1.0 2.0
      Figmax 0.3 0.3 0.3
      FG1max 0.5 0.5 0.5
      FG2min 0.5 0.0 0.0

      表 4  CL-20、CL-20/NTO、CL-20/FOX-7混合炸药的点火增长模型参数

      Table 4.  Ignition and growth model parameters for CL-20, CL-20/NTO and CL-20/FOX-7

      点火增长反应速率方程最初由两项组成[4],点火项和燃烧项,后来,为了更好地模拟高压、短脉冲下的冲击起爆过程,在两项式的基础上又增加了快速反应项[11],三项式点火增长反应速率方程为:

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

      式中:λ是炸药反应度,t是时间,ρ是密度,ρ0是初始密度,P是压力,IG1G2abxcdyegz是可调系数。

      采用点火增长反应速率方程模拟CL-20混合炸药的起爆过程,根据冲击起爆实验结果,标定出它的反应速率方程参数,如表 5中所示,模拟计算与实验结果的对比如图 5所示。

      炸药配方(质量比) 生成热/(J·g-1) 物质的量(分子量为100时)/mol
      C H N O F
      CL-20:NTO:FPM(47:47:6) -189.86 1.537 5 1.507 4 2.733 9 2.372 3 0.200 7
      CL-20:FOX-7:FPM(47:47:6) -141.53 1.449 5 2.054 6 2.558 0 2.558 0 0.200 7

      表 5  CL-20/NTO和CL-20/FOX-7两种混合炸药的配方和生成热数据

      Table 5.  Formula and heat of formation for the explosive of CL-20/NTO and CL-20/FOX-7

      图  5  CL-20混合炸药计算与实验结果对比

      Figure 5.  Experimental and calculated pressure histories in shock initiation of CL-20 explosive

      使用三项式点火增长反应速率方程模拟CL-20/NTO和CL-20/FOX-7两种混合炸药的冲击起爆过程,利用方程中第二项来描述CL-20的反应增长过程,而第三项来模拟两种混合炸药中反应相对较慢的NTO和FOX-7的反应释能过程。由于CL-20的反应较快,且两种混合炸药中CL-20的质量分数均为47%,因此,假设在反应度达到0.5时混合炸药中的CL-20就已经完成了全部反应,而NTO和FOX-7的反应速率相对较慢,它们的反应进程应贯穿始终。本次实验已得到CL-20单组分炸药的反应速率方程参数,而NTO和FOX-7则参考与其冲击感度相近的RDX和TATB炸药的模型参数[12],在得到每种组分炸药的两项式点火增长反应速率方程参数后,即可“拼”出双组分混合炸药的反应速率方程,其中,还需要根据实际情况对方程参数进行一定的调整。尝试采用两种不同的方法获得CL-20/NTO和CL-20/FOX-7混合炸药的反应速率方程参数,第1种方法选用反应较慢组分的点火项作为反应速率方程的第一项,选用两种组分的增长项作为反应速率方程的第2~3项,利用此方法,得到了CL-20/NTO混合炸药的参数,如表 5中所示。另一方法中,在CL-20混合炸药反应速率方程的基础上添加第3项,描述FOX-7的反应过程,然而,添加此项后,会使初期的反应速率增大,这明显与实际情况不符,因此,尝试对CL-20增长项系数G进行修正来解决此问题,参数的修正公式如下:

      $ \begin{array}{l} {G_{{\rm{new}}}} = {G_{{\rm{CL}} - 20}}{\eta _{{\rm{CL}} - 20}} + {G_{{\rm{FOX}} - 7}}{\eta _{{\rm{FOX}} - 7}}\\ \end{array} $

      式中:GCL-20GFOX-7分别为CL-20和FOX-7的增长项系数,ηCL-20和ηFOX-7为两种组分的质量分数,Gnew为修正后的系数。当GCL-20为402、GFOX-7为140、ηCL-20和ηFOX-7均为47%时,按公式(5)计算出修正后的增长项系数为254,CL-20/FOX-7混合炸药的反应速率方程参数如表 5中所示。

      图 6~7是CL-20/NTO和CL-20/FOX-7混合炸药的计算与实验结果对比,可以看出,计算结果较好地模拟出了两种混合炸药中不同位置处压力的成长趋势,表明利用上述方法得到的反应速率方程参数,可以较好地描述双组分混合炸药CL-20/NTO和CL-20/FOX-7的冲击起爆过程。因此,利用上述拟合双组分混合炸药点火增长反应速率方程的方法,可以对不同配比混合炸药的冲击起爆过程进行预测性计算,从而降低实验量及成本。

      图  6  CL-20/NTO混合炸药计算与实验结果对比

      Figure 6.  Experimental and calculated pressure histories in shock initiation of CL-20/NTO

      图  7  CL-20/FOX-7混合炸药计算与实验结果对比

      Figure 7.  Experimental and calculated pressure histories in shock initiation of CL-20/FOX-7

      非均质炸药冲击起爆的p2τ判据和临界能量流判据[13],是研究炸药冲击感度的重要参数,两种判据的具体形式为:

      $ {P^2}\tau = {\rm{con}}\;{\rm{st}} $

      $ {E_{\rm{c}}} = pu\tau = {\rm{con}}\;{\rm{st}} $

      式中:p为入射冲击波压力,u为冲击波后粒子速度,τ为入射冲击波脉冲宽度,Ec为临界起爆能量。此外,炸药在不同入射压力下的爆轰成长距离,即POP关系[14],也是反映炸药冲击感度的重要参数。

      这几种重要判据可以利用飞片撞击有隔板炸药的实验获得,因此,建立了飞片冲击起爆计算模型,尝试采用数值模拟的方法,计算CL-20、CL-20/NTO和CL-20/FOX-7三种混合炸药的冲击感度相关参数。计算模型由飞片、隔板、待测炸药组成,赋予飞片一定的速度,使其撞击隔板产生冲击波并起爆待测炸药,通过改变飞片的速度以获得不同的入射冲击波压力,计算分析待测炸药在不同入射压力下的起爆规律。钢飞片和铝隔板采用弹塑性流体力学材料模型和Grüneisen状态方程描述,使用点火增长模型描述待测炸药。钢飞片的厚度为3 mm,铝隔板厚6 mm,待测炸药高40 mm,三者直径均为50 mm。

      利用此计算模型对CL-20、CL-20/NTO和CL-20/FOX-7三种混合炸药进行了飞片冲击起爆数值模拟,计算得到了使炸药起爆的最低入射压力(临界压力)p,入射冲击波后粒子速度u以及飞片撞击隔板产生冲击波的脉宽τ,最终得到三种混合炸药冲击起爆的临界阈值,如表 6所示。通过改变飞片撞击隔板的速度,得到待测炸药在不同入射冲击压力下的爆轰成长距离,并在对数坐标下对样本点进行线性拟合,得到炸药的冲击起爆POP关系,如图 8所示。由以上结果可以看出,与CL-20混合炸药相比,CL-20/FOX-7混合炸药冲击起爆的p2τ判据和临界能量流判据只有少许提高,而添加了NTO的CL-20/NTO混合炸药的临界起爆阈值却有较大提升;从POP关系可以看出,向CL-20中分别添加FOX-7和NTO均可增加其爆轰成长距离,其中,在相同入射压力下,CL-20/FOX-7混合炸药的爆轰成长距离要大于CL-20/NTO混合炸药。综合这两种结果,认为在相同加载条件下,CL-20/NTO混合炸药内部形成的“热点”数量要少于CL-20、CL-20/FOX-7两种混合炸药,因此,其临界起爆阈值较高。然而,在“热点”的增长阶段,FOX-7炸药的反应速率要低于NTO炸药,从而导致CL-20/FOX-7混合炸药比CL-20/NTO混合炸药的爆轰成长距离更大。

      炸药(质量比) 密度/(g·cm-3) 临界压力/GPa p2τ/(Pa2·s) puτ/(J·m-2)
      CL-20(95/5) 1.94 1.39 253.67×1010 41.97×104
      CL-20/NTO/FPM(47/47/6) 1.89 1.46 279.86×1010 66.71×104
      CL-20/FOX-7/FPM(47/47/6) 1.88 1.40 257.33×1010 43.19×104

      表 6  三种混合炸药冲击起爆的临界阈值

      Table 6.  Shock initiation critical thresholds for three mixed explosives

      图  8  三种炸药的POP关系图

      Figure 8.  Pop plots for the three explosives

    • 通过在炸药内部嵌入锰铜压力传感器,并采用炸药驱动飞片冲击起爆的实验方法,获得了CL-20、CL-20/NTO和CL-20/FOX-7三种混合炸药的内部压力变化历程,利用点火增长模型,对炸药的冲击起爆过程进行数值模拟。研究发现,三种CL-20基混合炸药相比,CL-20/NTO混合炸药具有更高的临界起爆阈值,而在相同冲击加载条件下CL-20/FOX-7混合炸药具有更长的爆轰成长距离。在拟合双组分混合炸药CL-20/NTO和CL-20/FOX-7的反应速率方程过程中,使用两个增长项对两种炸药的反应分开描述,通过对比计算与实验结果,发现此方法得到的反应速率方程可以有效地模拟混合炸药的起爆过程。因此,利用此套拟合双组分混合炸药反应速率方程的方法,可以对新型配方炸药的冲击起爆过程进行预测性计算。

参考文献 (14)

目录

    /

    返回文章
    返回