利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

杨舒棋 张旭 彭文杨 舒俊翔 覃双 钟斌

引用本文:
Citation:

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

    作者简介: 杨舒棋(1993–),女,硕士研究生,主要从事冲击起爆研究. E-mail:yangshuqi77@126.com;
    通讯作者: 张旭, caepzx@sohu.com
  • 中图分类号: O384

Impact Initiation Characteristics of TATB Based Insensitive Explosives Mixed with HMX by Electromagnetic Velocity Gauges

    Corresponding author: ZHANG Xu, caepzx@sohu.com ;
  • CLC number: O384

  • 摘要: 为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。
  • 图 1  组合式电磁粒子速度计实验测试系统

    Figure 1.  Measurement system of aluminum-based multiple electromagnetic particle velocity gauge

    图 2  铝基组合式电磁粒子速度计的安装

    Figure 2.  Installation of aluminum-based multiple electromagnetic particle velocity gauge

    图 4  两条拟合直线交点细节

    Figure 4.  Details of two fitted straight lines’ intersection

    图 5  未反应 PBX-3炸药的Hugoniot关系

    Figure 5.  Hugoniot relationship of uncreated PBX-3 explosive

    图 6  PBX-3炸药的Pop关系

    Figure 6.  Pop-plot of PBX-3 explosive

    图 7  PBX-3和TATB-1的Pop关系

    Figure 7.  Pop-plots of PBX-3 and TATB-1 explosive

    图 8  爆轰反应区结构

    Figure 8.  Structure of the detonation reaction zone

    图 9  6个速度曲线修正

    Figure 9.  Correction of 6 speed curves

    表 1  PBX-3炸药平面冲击实验参数

    Table 1.  Parameters of plane impact experiments on PBX-3 explosive

    Shot No.m/gρ0/(g·cm–3) DS/(km·s–1)up/(km·s–1)p0/GPa
    01 5001.9003.9660.724 5.456
    02 6001.9004.2830.830 6.754
    03 7001.9004.5161.019 8.743
    04 8001.9004.4921.096 9.354
    05 9001.9004.8081.13210.341
    061 1001.9005.2811.29212.964
    下载: 导出CSV

    表 2  PBX-3炸药的 DS-up实验数据

    Table 2.  Experimental results of DS-up for PBX-3 explosive

    Depth/mmShot 01Shot 02Shot 03
    up/(km·s–1)DS/(km·s–1)up/(km·s–1)DS/(km·s–1)up/(km·s–1)DS/(km·s–1)
    00.7244.1890.8304.2241.0193.489
    11.1043.489
    21.0494.259
    30.7764.1890.8944.2241.0334.686
    40.7373.9560.9114.5211.1184.726
    50.7323.7960.9384.2291.2824.546
    60.7663.8641.0344.3181.4744.799
    70.8333.9501.1874.5791.6825.252
    80.8734.0071.2794.6041.9575.165
    90.9374.1901.4505.470
    100.9834.2801.5984.960
    Depth/mmShot 04Shot 05Shot 06
    up/(km·s–1)DS/(km·s–1)up/(km·s–1)DS /(km·s–1)up/(km·s–1)DS/(km·s–1)
    01.0964.1641.1314.9551.2925.152
    11.1134.164
    21.1224.119
    31.1374.6801.3474.9551.6345.152
    41.2324.6821.5454.6991.9985.774
    51.4074.6641.6544.9472.2306.187
    61.6374.9901.8675.1442.0667.353
    71.9625.2132.1045.6562.2827.082
    82.5105.8682.1377.0832.1207.043
    92.3086.746
    102.3186.6992.2837.374
    下载: 导出CSV

    表 3  PBX-3炸药Pop关系相关参数

    Table 3.  Related parameters of PBX-3 explosive’s Pop-plot

    Shot No.ρ0/(g·cm–3)p0/GPaxD/mmtD/μs
    011.900 5.456
    021.900 6.7549.4392.224
    031.900 8.7437.3291.674
    041.900 9.3546.2571.421
    051.90010.3416.2341.292
    061.90012.9644.1430.785
    下载: 导出CSV

    表 4  PBX-3炸药化学反应区实验参数

    Table 4.  Parameters of chemical reaction zone of PBX-3 explosive

    Depth/mmup/(km·s–1)t/μsx/mmDepth/mmup/(km·s–1)t/μsx/mm
    51.5390.2261.429 81.4250.2261.454
    61.5400.2261.428 91.3700.2261.467
    71.4360.2261.452101.3860.2261.463
    下载: 导出CSV

    表 5  TATB基炸药反应区相关参数

    Table 5.  Related parameters of reaction zone of TATB-based explosive

    Shot No.x/mmt/μsSource
    PBX-31.449 ± 0.2000.226 ± 0.030This work
    PBX95022.10.28Ref.[8]
    JB-90141.5 ± 0.20.26 ± 0.02Ref.[9]
    JB-90141.750.31Ref.[10]
    下载: 导出CSV
  • [1] 张琪敏, 张旭, 赵康, 等. TATB基钝感炸药JB-9014的冲击起爆反应增长规律 [J]. 爆炸与冲击, 2019, 39(4): 041405.
    ZHANG Q M, ZHANG X, ZHAO K, et al. Law of reaction growth of shock initiation on the TATB based insensitive explosive JB-9014 [J]. Explosion and Shock Waves, 2019, 39(4): 041405.
    [2] HILL L G, GUSTAVSEN R L, ALCON R R, et al. Shock initiation of new and aged PBX 9501 measured with embedded electromagnetic particle velocity gauges: LA-13634-MS [R]. New Mexico, US: Los Alamos National Laboratory, 1999.
    [3] ZHANG X, WANG Y F, HUANG W B, et al. Reaction buildup of PBX explosives JOB-9003 under different initiation pressures [J]. Journal of Energetic Materials, 2017, 35(2): 197–212. doi: 10.1080/07370652.2016.1250841
    [4] GUSTAVSEN R L, GEHR R J, BUCHOLTZ S M, et al. Shock initiation of the tri-amino-tri-nitro-benzene explosive PBX9502 cooled to -55 ℃ [J]. Journal of Applied Physics, 2012, 112(7): 074909. doi: 10.1063/1.4757599
    [5] 张涛, 赵继波, 伍星, 等. 未反应JBO-9021炸药冲击雨贡纽曲线的研究 [J]. 高压物理学报, 2016, 30(6): 457–462. doi: 10.11858/gywlxb.2016.06.004
    ZHANG T, ZHAO J B, WU X, et al. Hugoniot curve of unreacted JBO-9021 explosive [J]. Chinese Journal of High Pressure Physics, 2016, 30(6): 457–462. doi: 10.11858/gywlxb.2016.06.004
    [6] GUSTAVSEN R L, SHEFFIELD S A, ALCON R R. Measurements of shock initiation in the tri-amino-tri-nitro-benzene based explosive PBX 9502: wave forms from embedded gauges and comparison of four different material lots [J]. Journal of Applied Physics, 2006, 99(11): 114907. doi: 10.1063/1.2195191
    [7] 张涛, 谷岩, 赵继波, 等. JBO-9021炸药的化学反应区宽度 [J]. 爆炸与冲击, 2017, 37(3): 415–421. doi: 10.11883/1001-1455(2017)03-0415-07
    ZHANG T, GU Y, ZHAO J B, et al. Chemical reaction zone length of JBO-9021 [J]. Explosion and Shock Waves, 2017, 37(3): 415–421. doi: 10.11883/1001-1455(2017)03-0415-07
    [8] SHEFFIELD S A, BLOOMQUIST D D, TARVER C M. Subnanosecond measurements of detonation fronts in solid high explosives [J]. The Journal of Chemical Physics, 1984, 80(8): 3831–3844. doi: 10.1063/1.447164
    [9] 裴红波, 黄文斌, 覃锦程, 等. 基于多普勒测速技术的JB-9014炸药反应区结构研究 [J]. 爆炸与冲击, 2018, 38(3): 485–490.
    PEI H B, HUANG W B, QIN J C, et al. Reaction zone structure of JB-9014 explosive measured by PDV [J]. Explosion and Shock Waves, 2018, 38(3): 485–490.
    [10] 赵同虎, 张新彦, 李斌, 等. 用光电法研究钝感炸药JB-9014反应区结构 [J]. 高压物理学报, 2002, 16(2): 111–119. doi: 10.3969/j.issn.1000-5773.2002.02.005
    ZHAO T H, ZHANG X Y, LI B, et al. Detonation reaction zones tructure of JB-9014 [J]. Chinese Journal of High Pressure Physics, 2002, 16(2): 111–119. doi: 10.3969/j.issn.1000-5773.2002.02.005
    [11] LOBOIKO B G, LUBYATINSKY S N. Reaction zones of detonating solid explosives [J]. Combustion, Explosion, and Shock Waves, 2000, 36(6): 716–733. doi: 10.1023/A:1002898505288
  • [1] 李金河訾攀登张旭张蓉 . 用组合式电磁粒子速度计研究一种活性材料的反应特性. 高压物理学报, 2017, 31(3): 309-314. doi: 10.11858/gywlxb.2017.03.013
    [2] 刘福生经福谦 . 沿等压路径求解疏松材料Hugoniot关系的微分方程组及其求解. 高压物理学报, 2004, 18(1): 10-16 . doi: 10.11858/gywlxb.2004.01.003
    [3] 杨舒棋张旭彭文杨舒俊翔刘寿先覃双钟斌 . 钝感炸药冲击起爆反应过程的PDV技术. 高压物理学报, 2020, 34(2): 023402-1-023402-9. doi: 10.11858/gywlxb.20190856
    [4] 彭克锋潘昊赵凯郑志军虞吉林 . 铜粉末动态压缩行为的多颗粒有限元分析. 高压物理学报, 2019, 33(4): 044102-1-044102-8. doi: 10.11858/gywlxb.20180665
    [5] 吴强经福谦李欣竹 . 零温物态方程输入参数B0K、B0K和0K的确定. 高压物理学报, 2005, 19(2): 97-104 . doi: 10.11858/gywlxb.2005.02.001
    [6] 傅华谭多望李金河李涛 . 未反应JOB-9003炸药冲击Hugoniot关系测试. 高压物理学报, 2009, 23(6): 427-432 . doi: 10.11858/gywlxb.2009.06.005
    [7] 刘俊明张旭裴红波舒俊翔覃双钟斌张蓉 . JB-9014钝感炸药冲击Hugoniot关系测量. 高压物理学报, 2018, 32(3): 033202-1-033202-7. doi: 10.11858/gywlxb.20170669
    [8] 于川池家春刘文翰李良忠杨淑英 . JB-9001钝感炸药冲击Hugoniot关系测试. 高压物理学报, 1998, 12(1): 72-77 . doi: 10.11858/gywlxb.1998.01.012
    [9] 刘晓旸苏文辉孟长功胡壮麒 . 岩盐结构固体的硬度与化学键特征的关系. 高压物理学报, 1994, 8(4): 241-247 . doi: 10.11858/gywlxb.1994.04.001
    [10] 范宝春崔东民陈启峰 . 驻定激波诱导的化学反应. 高压物理学报, 1997, 11(3): 182-188 . doi: 10.11858/gywlxb.1997.03.004
    [11] 戴诚达王翔谭华 . Hugoniot实验的粒子速度测量不确定度分析. 高压物理学报, 2005, 19(2): 113-119 . doi: 10.11858/gywlxb.2005.02.003
    [12] 洪时明 . 高压相变与时间的关系. 高压物理学报, 2013, 27(2): 162-167. doi: 10.11858/gywlxb.2013.02.002
    [13] 王刚张德良刘凯欣 . 气相爆轰波数值模拟中化学反应模型研究. 高压物理学报, 2008, 22(4): 350-356 . doi: 10.11858/gywlxb.2008.04.003
    [14] 田雨刘雪廷何运鸿赵慧芳姜峰谭大勇肖万生 . NaCl-O2体系高温高压化学反应的拉曼光谱证据. 高压物理学报, 2017, 31(6): 692-697. doi: 10.11858/gywlxb.2017.06.003
    [15] 胡金彪谭华经福谦 . 溴仿在冲击压缩下的光辐射及化学反应. 高压物理学报, 1996, 10(3): 220-226 . doi: 10.11858/gywlxb.1996.03.009
    [16] 唐旭东张增明赵智王中平丁泽军 . 高压诱致的反式联苯乙烯酮结构相变与化学反应. 高压物理学报, 2009, 23(6): 401-406 . doi: 10.11858/gywlxb.2009.06.001
    [17] 赵慧芳谭大勇姜峰谢亚飞姜昌国罗兴丽肖万生 . Re-H2O体系高温高压化学反应的拉曼证据. 高压物理学报, 2020, 34(4): 040102-1-040102-10. doi: 10.11858/gywlxb.20200518
    [18] 王武孝介万奇袁森 . 制备Al2O3/Cu复合材料中CuO-Al体系的化学反应. 高压物理学报, 2004, 18(4): 353-358 . doi: 10.11858/gywlxb.2004.04.011
    [19] 严祖同 . 碱卤晶体热导率与压强的关系. 高压物理学报, 1993, 7(1): 66-70 . doi: 10.11858/gywlxb.1993.01.010
    [20] 龙期威 . 分形图形周界和面积的关系. 高压物理学报, 1990, 4(4): 259-262 . doi: 10.11858/gywlxb.1990.04.004
  • 加载中
图(9)表(5)
计量
  • 文章访问数:  1953
  • 阅读全文浏览量:  1088
  • PDF下载量:  12
出版历程
  • 收稿日期:  2019-11-01
  • 录用日期:  2019-11-12
  • 网络出版日期:  2020-04-09
  • 刊出日期:  2020-06-01

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

    作者简介:杨舒棋(1993–),女,硕士研究生,主要从事冲击起爆研究. E-mail:yangshuqi77@126.com
    通讯作者: 张旭, caepzx@sohu.com
  • 中国工程物理研究院流体物理研究所,四川 绵阳 621999

摘要: 为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。

English Abstract

  • 钝感炸药(Insensitive high explosive, IHE)的安全性能是弹药设计的重点。以三氨基三硝基苯(TATB)为基的钝感炸药在国防领域有着广泛应用,PBX-3炸药为我国研制的含有少量奥克托金(HMX)的TATB基新型高能钝感塑性黏结炸药,开展其性能研究对炸药的应用和改进具有重要意义。

    PBX-3炸药含有少量HMX,冲击起爆感度更高,当入射压力为6.754 GPa时,炸药样品在9.439 mm处就开始稳定爆轰,较TATB基PBX9502炸药和JB-9014炸药具有更优异的爆轰性能和冲击感度。张琪敏等[1]研究了纯TATB基钝感炸药在入射压力为11.330~14.180 GPa范围内的冲击起爆反应增长规律,结果表明,当入射压力达到11.334 GPa时,炸药样品在9.590 mm处才开始稳定爆轰。同时,PBX-3炸药的钝感性能优于HMX基PBX9501炸药和JOB-9003炸药,因此制作加工和运输使用更加安全。

    国内外研究炸药性能常用的测试方法有锰铜压力计法、PVDF压电传感器法、组合式电磁粒子速度计测试技术和激光干涉测试技术等。锰铜压力计的响应时间长,测量精度受到压阻系数、绝缘保护膜等因素影响。PVDF压电传感器的测试上限为35 GPa,当压力大于10 GPa时不能忽略温度的影响。组合式电磁粒子速度计技术和激光干涉技术是近些年国内外广泛使用的测试技术[2-4]。激光干涉测试技术只能测量冲击加载界面的粒子速度,对于普通圆柱形炸药样品,无法在一次冲击实验中得到不同深度处样品的粒子速度,为此炸药样品的形状需特殊加工(如台阶形或楔形),并且炸药到爆轰的距离和时间需要采用其他测试方法。张涛等[5]采用激光干涉测速系统得到了冲击波进入炸药特定深度位置处的波后粒子速度剖面,同时结合高速扫描相机得到了炸药样品到爆轰的距离。组合式电磁粒子速度计具有实验操作简便、响应灵敏的特点,采用的铝金属丝更加接近炸药样品的阻抗,从而减小对流场的破坏。本研究采用铝基组合式电磁粒子速度计技术测量楔形炸药样品的撞击表面及内部不同深度处的冲击波后粒子速度,利用示踪器得到前沿冲击波位置随时间的变化关系,通过一次样品冲击实验即可为研究炸药样品的Hugoniot关系、Pop关系甚至化学反应区提供丰富的实验数据。

    • 实验采用火炮加载系统驱动蓝宝石飞片撞击炸药,进行一维平面冲击实验。实验装置主要包括火炮驱动系统、铝弹托、蓝宝石飞片、激光测速环等,实验测试系统如图1所示。通过改变火药装药量控制蓝宝石飞片速度,进而控制初始入射压力。

      图  1  组合式电磁粒子速度计实验测试系统

      Figure 1.  Measurement system of aluminum-based multiple electromagnetic particle velocity gauge

      PBX-3炸药样品是$\varnothing $42 mm × 30 mm的圆柱形组合药柱。为测量不同深度处炸药样品的冲击波后粒子速度,炸药样品被切割成两块30°角的楔形,用环氧树脂把组合式电磁粒子速度计嵌入两块楔形炸药中。利用蓝宝石飞片通过激光测速环时遮挡3组激光束的时间差来计算飞片速度,同时在炸药表面用环氧树脂粘贴3个电磁粒子速度计,用于测量蓝宝石飞片撞击炸药表面时的瞬时界面粒子速度。在冲击波传播方向,利用铝基组合式电磁粒子速度计(8个电磁粒子速度计)测量炸药样品的波后粒子速度,其深度间隔为1 mm,3个冲击波示踪器元件的间距为0.5 mm,组合式电磁粒子速度计的安装如图2所示。

      图  2  铝基组合式电磁粒子速度计的安装

      Figure 2.  Installation of aluminum-based multiple electromagnetic particle velocity gauge

    • 在冲击波作用下,组合式电磁粒子速度计在均匀磁场中切割磁感线产生感应电动势。根据法拉第电磁感应定律,粒子速度与感应电动势的关系为

      $u = \frac{E}{{Bl}}$

      式中:u为粒子速度,m/s;E为电动势,V;B为磁场强度,T;l为切割磁感线长度,m。

      因为感应电动势由电磁粒子速度计电阻、连接的导线以及示波器自带电阻以串联形式共同分担表示,所以在计算炸药不同深度处电磁粒子速度计测量得到的粒子速度时,需要乘以补偿系数

      $\frac{E}{{{E_{{\rm{oscilloscope}}}}}} = \frac{{{R_{{\rm{oscilloscope}}}} + {{\left( {R - \Delta R} \right)}}}}{{{R_{{\rm{oscilloscope}}}}}}$

      式中:Eoscilloscope为实验中由示波器测量得到的感应电动势,Roscilloscope为示波器的电阻(50 Ω),R为实验中电磁粒子速度计和连接导线本身电阻之和(实验使用万用表测量电阻,故计算时需减去万用表本底电阻),$\Delta R$ 为万用表本身的电阻。

      结合式(1)和式(2),得到粒子速度计算公式为

      $u = \frac{{{E_{{\rm{oscilloscope}}}}}}{{Bl}} \left( {1 + \frac{{{{R - \Delta R}}}}{{{R_{{\rm{oscilloscope}}}}}}} \right)$

    • 通过改变火炮加载系统的火药装药量,得到递增的蓝宝石飞片撞击速度。6发实验的火药装药量为500、600、700、800、900和1 100 g,分别命名为Shot 01、Shot 02、Shot 03、Shot 04、Shot 05、Shot 06。利用阻抗匹配法分析计算蓝宝石飞片撞击炸药过程产生的入射压力

      $ {p_{0{\rm{L}}}} = {\rho _{0{\rm{L}}}}{D_{\rm{L}}}({u_{{\rm{imp}}}} - {u_{\rm{p}}}) $

      $ {p_{0{\rm{S}}}} = {\rho _{0{\rm{S}}}}{D_{\rm{S}}}{u_{\rm{p}}} $

      式中:p0L为蓝宝石初始入射冲击压力,GPa;p0S为炸药样品初始入射冲击压力,GPa;ρ0L为蓝宝石初始密度,g/cm3ρ0S为炸药样品初始密度,g/cm3DL为蓝宝石中的冲击波速度,km/s;DS为炸药样品的入射冲击波速度,km/s;uimp为蓝宝石飞片的撞击速度,km/s;up为波后粒子速度,km/s。

      根据蓝宝石材料和PBX-3炸药的Hugoniot关系,将式(4)写为

      $ {p_{0{\rm{L}}}} = {\rho _{0{\rm{L}}}}[{C_{0{\rm{L}}}} + {\lambda _{\rm{L}}}({u_{{\rm{imp}}}} - {u_{\rm{p}}})]({u_{{\rm{imp}}}} - {u_{\rm{p}}}) $

      式中:蓝宝石的初始密度ρ0L = 3.985 g/cm3;Hugoniot系数C0L = 11.19 km/s,λL = 1.0。

      由界面连续性可知p0L = p0S,结合式(5)和式(6),得到未反应炸药的波后粒子速度up,并计算得出实验的初始入射冲击压力p0,结果如表1所示。

      Shot No.m/gρ0/(g·cm–3) DS/(km·s–1)up/(km·s–1)p0/GPa
      01 5001.9003.9660.724 5.456
      02 6001.9004.2830.830 6.754
      03 7001.9004.5161.019 8.743
      04 8001.9004.4921.096 9.354
      05 9001.9004.8081.13210.341
      061 1001.9005.2811.29212.964

      表 1  PBX-3炸药平面冲击实验参数

      Table 1.  Parameters of plane impact experiments on PBX-3 explosive

    • 利用式(3)对实验数据进行处理,得到炸药样品的波后粒子速度,根据示踪器数据得到炸药样品到爆轰的位置与时间的关系,结果如图3所示。图3(a)所示为Shot 01实验,入射压力为5.456 GPa。可见,随着样品发生化学反应释放能量,粒子速度上升后缓慢增长,粒子速度达到最大值后,随着化学反应结束又逐渐减小,速度信号呈“驼峰”形状。随着冲击波强度的增加,炸药反应速率不断加快,“驼峰”形状中速度上升部分加速减小,但并没有转变为爆轰。图3(b)所示为Shot 02实验,入射压力为6.745 GPa,粒子速度达到最大值的时间间隔逐渐减小,“驼峰”形状逐渐尖锐,说明此时化学反应释放的能量使冲击波逐渐加强。当入射压力达到8.743 GPa时,如图3(c)所示的Shot 03实验,距离炸药样品8 mm处的粒子速度曲线出现迅速增长,达到最大速度后转变为逐渐减小的三角波形,说明炸药样品形成稳定的爆轰波,完成了冲击转爆轰过程(SDT)。

    • 实际应用中,炸药样品的入射冲击波速度DS与粒子速度up之间近似遵循线性关系

      $ {D_{\rm{S}}} = {C_0} + \lambda {u_{\rm{p}}} $

      式中:DS为炸药样品中的入射冲击波速度,km/s;C0λ为常数;up为波后粒子速度,km/s。

      由于相邻两个速度计间隔是已知的,从粒子速度曲线上读取冲击波到达相邻两个测试点的时间即可计算得到炸药样品不同深度处的入射冲击波速度 DS。如表2所示,随着冲击波在炸药样品内传播,冲击波波速逐渐增大。以图3(d)中3 mm处粒子速度曲线为例,利用图4所示的方法[6],将粒子速度分成两部分进行线性拟合,两条拟合直线的交点则为Hugoniot关系中的波后粒子速度up表2为6发实验所测PBX-3炸药的 DS-up实验数据,选取炸药样品达到爆轰之前的部分实验数据,绘制成如图5所示的冲击Hugoniot关系图,并得到PBX-3炸药的冲击Hugoniot关系式

      Depth/mmShot 01Shot 02Shot 03
      up/(km·s–1)DS/(km·s–1)up/(km·s–1)DS/(km·s–1)up/(km·s–1)DS/(km·s–1)
      00.7244.1890.8304.2241.0193.489
      11.1043.489
      21.0494.259
      30.7764.1890.8944.2241.0334.686
      40.7373.9560.9114.5211.1184.726
      50.7323.7960.9384.2291.2824.546
      60.7663.8641.0344.3181.4744.799
      70.8333.9501.1874.5791.6825.252
      80.8734.0071.2794.6041.9575.165
      90.9374.1901.4505.470
      100.9834.2801.5984.960
      Depth/mmShot 04Shot 05Shot 06
      up/(km·s–1)DS/(km·s–1)up/(km·s–1)DS /(km·s–1)up/(km·s–1)DS/(km·s–1)
      01.0964.1641.1314.9551.2925.152
      11.1134.164
      21.1224.119
      31.1374.6801.3474.9551.6345.152
      41.2324.6821.5454.6991.9985.774
      51.4074.6641.6544.9472.2306.187
      61.6374.9901.8675.1442.0667.353
      71.9625.2132.1045.6562.2827.082
      82.5105.8682.1377.0832.1207.043
      92.3086.746
      102.3186.6992.2837.374

      表 2  PBX-3炸药的 DS-up实验数据

      Table 2.  Experimental results of DS-up for PBX-3 explosive

      图  4  两条拟合直线交点细节

      Figure 4.  Details of two fitted straight lines’ intersection

      图  5  未反应 PBX-3炸药的Hugoniot关系

      Figure 5.  Hugoniot relationship of uncreated PBX-3 explosive

      $ {D_{\rm{S}}} = (3.102 \pm 0.140) + (1.234 \pm 0.115){u_{\rm{p}}} $

    • 炸药的初始入射压力与到爆轰距离xD(到爆轰时间tD)的关系称为Pop关系,反映了炸药的冲击起爆性能。组合式电磁粒子速度计的冲击波示踪器呈“梯子锯齿”状,示踪器元件切割磁感线产生的电压会发生正负周期性变化。在冲击波传播方向,冲击波示踪器中每个测量元件的间距为0.5 mm,结合正负电压信号转变的时间间隔,可以得到冲击波速度。根据示踪器信号,可以绘制出显示前沿冲击波位置随时间的变化关系(x-t图),如图3所示。由图3中的实验数据拟合出两条直线,其斜率分别表示未反应炸药的平均冲击速度和爆轰速度,两条拟合线的交点表示到爆轰的位置和时间,结果如表3所示。根据表3数据可以得到初始入射压力在5.456 ~ 12.964 GPa范围内的PBX-3炸药的Pop关系拟合线,如图6所示。Pop关系表达式为

      Shot No.ρ0/(g·cm–3)p0/GPaxD/mmtD/μs
      011.900 5.456
      021.900 6.7549.4392.224
      031.900 8.7437.3291.674
      041.900 9.3546.2571.421
      051.90010.3416.2341.292
      061.90012.9644.1430.785

      表 3  PBX-3炸药Pop关系相关参数

      Table 3.  Related parameters of PBX-3 explosive’s Pop-plot

      图  6  PBX-3炸药的Pop关系

      Figure 6.  Pop-plot of PBX-3 explosive

      $ {\rm{lg}}\;{x_{\rm{D}}} = (2.013 \pm 0.133) - (1.236 \pm 0.136){\rm{lg}}\;{p_0} $

      为更好地研究PBX-3炸药的冲击起爆性能,将纯TATB基炸药TATB-1的实验数据和PBX-3进行对比,如图7所示。由图7可知,PBX-3的Pop关系曲线在TATB-1的下方,在相同入射压力条件下,PBX-3到爆轰的时间更短。PBX-3是少量HMX和TATB混合的钝感炸药,由于HMX的爆轰压力和爆轰速度比TATB高,故PBX-3炸药的冲击起爆性能比纯TATB基炸药TATB-1更敏感。

      图  7  PBX-3和TATB-1的Pop关系

      Figure 7.  Pop-plots of PBX-3 and TATB-1 explosive

    • ZND(Zeldovich-von Neumann-Doring)模型是一种描述炸药爆轰波剖面的经典模型。模型假设爆轰波阵面是由前导冲击波和具有一定厚度的化学反应区组成。对于稳定传播的爆轰波,其终态点为Chapman-Jouguet(C-J)点,化学反应终止时压力急剧下降至C-J点,C-J面后为等熵膨胀流动区,即泰勒波稀疏区,如图8所示。炸药样品发生稳定爆轰时,冲击波后粒子速度突跃至最大速度,随着化学反应的进行,从粒子速度曲线斜率的变化可知,样品粒子加速度由开始的快速下降逐渐变成缓慢下降。化学反应结束后,样品粒子受到泰勒波稀疏影响,在终态C-J点处发生分散,最终粒子加速度趋于零。

      图  8  爆轰反应区结构

      Figure 8.  Structure of the detonation reaction zone

      图3(f)所示,组合式电磁粒子速度计测得的粒子速度曲线存在信号噪声,但PBX-3炸药界面粒子速度有较明显的拐点,可以通过读取炸药界面粒子速度的拐点(即受到泰勒波稀疏影响的粒子速度曲线分散点)粗略地确定炸药的C-J点,如图9所示。

      图  9  6个速度曲线修正

      Figure 9.  Correction of 6 speed curves

      为了更准确地测量反应区时间,选择了爆轰波形一致性较好的实验数据。在Shot 06实验(入射压力为12.964 GPa)中,炸药样品在4.143 mm处开始发生稳定爆轰,将达到爆轰的6条速度曲线(即5~10 mm深度处的粒子速度数据)修整成相同零点。在0.226 μs处6条实验曲线分离,分离点即为反应区末端相应的C-J点,由此可得化学反应时间为0.226 μs。确定炸药化学反应区宽度的关系式[7]

      $x = \int_0^\tau {\left( {{D_{{\rm{CJ}}}} - {{\bar u}_{\rm{P}}}} \right)} {\rm{d}}t$

      式中:τ为化学反应区持续时间,为0.226 μs;DCJ为PBX-3炸药样品的C-J爆速,为7.86 km/s;${\bar u}_{\rm{P}}$为化学反应区持续时间内反应产物粒子平均速度(利用软件将速度曲线积分,得到反应区持续时间内反应产物粒子运动的位移,结合化学反应区持续时间,得到反应产物的粒子平均速度)。表4所列数据为计算得到的化学反应区宽度,均值为1.449 mm。国内外学者[8-10]主要采用速度和时间分辨率较高的激光干涉技术,结合状态方程、数值模拟和粒子速度求导等方法进行TATB基炸药反应区的研究。表5列出了国内外学者得到的TATB基炸药的反应区宽度和时间。由于PBX-3炸药含有少量HMX,因而其化学反应区的持续时间和反应区宽度较纯TATB基炸药[11]更短、更窄。

      Depth/mmup/(km·s–1)t/μsx/mmDepth/mmup/(km·s–1)t/μsx/mm
      51.5390.2261.429 81.4250.2261.454
      61.5400.2261.428 91.3700.2261.467
      71.4360.2261.452101.3860.2261.463

      表 4  PBX-3炸药化学反应区实验参数

      Table 4.  Parameters of chemical reaction zone of PBX-3 explosive

      Shot No.x/mmt/μsSource
      PBX-31.449 ± 0.2000.226 ± 0.030This work
      PBX95022.10.28Ref.[8]
      JB-90141.5 ± 0.20.26 ± 0.02Ref.[9]
      JB-90141.750.31Ref.[10]

      表 5  TATB基炸药反应区相关参数

      Table 5.  Related parameters of reaction zone of TATB-based explosive

    • 采用火炮驱动蓝宝石飞片方法和铝基组合式电磁粒子速度计技术,对高能钝感炸药PBX-3进行了一维平面冲击实验。利用11个电磁粒子速度计测量得到炸药撞击表面及内部不同深度处的冲击波后粒子速度,通过3个冲击波示踪器绘制了炸药到爆轰距离与时间的x-t关系曲线。通过实验撞击表面及内部不同深度处冲击波后粒子速度up和入射冲击波速度DS,得到PBX-3炸药的Hugoniot关系为 DS =(3.102 ± 0.140) + (1.234 ± 0.115)up。根据炸药到爆轰时间与距离的x-t关系曲线和初始入射压力,得到反映炸药冲击起爆性能的Pop关系为lg xD = (2.013 ± 0.133)-(1.236 ± 0.136)lg p0。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点,得到化学反应区时间为(0.226 ± 0.030) μs,化学反应区宽度为(1.449 ± 0.200) mm。

参考文献 (11)

目录

    /

    返回文章
    返回