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

CL-20/HMX共晶冲击起爆模拟

刘海 李毅 李俊玲 马兆侠 陈鸿

引用本文:
Citation:

CL-20/HMX共晶冲击起爆模拟

    作者简介: 刘 海(1985- ),男,博士,助理研究员,liumy2016@163.com;
  • 中图分类号: O389

Simulations of shock initiation of CL-20/HMX co-crystal

  • CLC number: O389

  • 摘要: 采用非平衡分子动力学方法模拟了CL-20/HMX共晶炸药冲击压缩和化学反应行为,获得了密度以及粒子速度的时空分布、冲击雨贡纽、冲击起爆压力、爆轰压力等数据,以及主要中间产物和稳定产物分布。模拟结果显示,共晶的初始分解路径是CL-20中N—NO2键断裂,主要稳定产物是N2、CO2和H2O。CL-20和HMX的分解速率随着冲击波速度的增加而增加,并且逐渐接近,但各冲击条件下CL-20分子的衰减速率均大于HMX。
  • 图 1  CL-20/HMX模型以及动量镜激发冲击波传播原理图

    Figure 1.  CL-20/HMX structure model and schematic diagram of the shock wave propagation induced by momentum mirror

    图 2  冲击波传播引发粒子速度和密度变化的x-t

    Figure 2.  x-t diagrams of particle velocity and density induced by shock wave propagation

    图 3  冲击波速度-粒子速度关系及p-V/V0雨贡纽

    Figure 3.  Shock wave velocity-particle velocity relation and p-V/V0 Hugoniot

    图 4  共晶中CL-20和HMX分子衰减曲线

    Figure 4.  Decay curves of CL-20 and HMX molecules in co-crystal

    图 5  诱导时间和衰减速率及冲击到达右端自由面时主要产物的数量统计

    Figure 5.  Induction time and decay rate under different shock conditions and quantitative statistics of the main products when the shock wave reaches the free surface

    图 6  冲击波传播的序列图像

    Figure 6.  Sequence images of shock wave propagation

    表 1  各冲击条件下共晶中CL-20和HMX诱导时间和衰减速率比较

    Table 1.  Comparison of induction time and decay rate of CL-20 and HMX in co-crystal under different shock conditions

    up /(km·s−1)τCL-20/psτHMX/psrd, CL-20rd, HMXrd, CL-20/rd, HMX
    1.50.800 02.100 00.004 00.001 13.636 4
    2.00.254 00.800 00.016 00.002 36.956 5
    2.50.100 00.400 00.062 00.031 02.000 0
    3.00.052 00.300 00.160 00.114 01.403 5
    3.50.049 00.160 00.234 00.205 01.141 5
    4.00.001 00.150 00.289 00.278 01.039 6
    下载: 导出CSV
  • [1] BOLTON O, SIMKE L R, PAGORIA P F, et al. High power explosive with good sensitivity: a 2∶1 cocrystal of CL-20∶HMX [J]. Crystal Growth & Design, 2012, 12(9): 4311–4314.
    [2] SUN T, XIAO J J, LIU Q, et al. Comparative study on structure, energetic and mechanical properties of a ε-CL-20/HMX cocrystal and its composite with molecular dynamics simulation [J]. Journal of Materials Chemistry A, 2014, 2(34): 13898–13904. DOI: 10.1039/C4TA01150C.
    [3] LIU Z, WU Q, ZHU W, et al. Insights into the roles of two constituents CL-20 and HMX in the CL-20: HMX cocrystal at high pressure: a DFT-D study [J]. RSC Advances, 2015, 5(43): 34216–34225. DOI: 10.1039/C5RA01829C.
    [4] XUE X, MA Y, ZENG Q, et al. Initial decay mechanism of the heated CL-20/HMX cocrystal: a case of the cocrystal mediating the thermal stability of the two pure components [J]. The Journal of Physical Chemistry C, 2017, 121(9): 4899–4908. DOI: 10.1021/acs.jpcc.7b00698.
    [5] DOBLAS D, ROSENTHAL M, BURGHAMMER M, et al. Smart energetic nanosized co-crystals: exploring fast structure formation and decomposition [J]. Crystal Growth & Design, 2015, 16(1): 432–439.
    [6] OKOVYTYY S, KHOLOD Y, QASIM M, et al. The mechanism of unimolecular decomposition of 2, 4, 6, 8, 10, 12-hexanitro-2, 4, 6, 8, 10, 12-hexaazaisowurtzitane: a computational DFT study [J]. The Journal of Physical Chemistry A, 2005, 109(12): 2964–2970. DOI: 10.1021/jp045292v.
    [7] ISAYEV O, GORB L, QASIM M, et al. Ab initio molecular dynamics study on the initial chemical events in nitramines: thermal decomposition of CL-20 [J]. The Journal of Physical Chemistry B, 2008, 112(35): 11005–11013. DOI: 10.1021/jp804765m.
    [8] WANG F, CHEN L, GENG D, et al. Effect of density on the thermal decomposition mechanism of ε-CL-20: a ReaxFF reactive molecular dynamics simulation study [J]. Physical Chemistry Chemical Physics, 2018, 20(35): 22600–22609. DOI: 10.1039/C8CP03010C.
    [9] WANG F, CHEN L, GENG D, et al. Thermal decomposition mechanism of CL-20 at different temperatures by ReaxFF reactive molecular dynamics simulations [J]. The Journal of Physical Chemistry A, 2018, 122(16): 3971–3979. DOI: 10.1021/acs.jpca.8b01256.
    [10] XUE X, WEN Y, ZHANG C. Early decay mechanism of shocked ε-CL-20: a molecular dynamics simulation study [J]. The Journal of Physical Chemistry C, 2016, 120(38): 21169–21177. DOI: 10.1021/acs.jpcc.6b05228.
    [11] WEN Y, XUE X, ZHOU X, et al. Twin induced sensitivity enhancement of HMX versus shock: a molecular reactive force field simulation [J]. The Journal of Physical Chemistry C, 2013, 117(46): 24368–24374. DOI: 10.1021/jp4072795.
    [12] GE N N, WEI Y K, JI G F, et al. Initial decomposition of the condensed-phase β-HMX under shock waves: molecular dynamics simulations [J]. The Journal of Physical Chemistry B, 2012, 116(46): 13696–13704. DOI: 10.1021/jp309120t.
    [13] GE N N, WEI Y K, SONG Z F, et al. Anisotropic responses and initial decomposition of condensed-phase β-HMX under shock loadings via molecular dynamics simulations in conjunction with multiscale shock technique [J]. The Journal of Physical Chemistry B, 2014, 118(29): 8691–8699. DOI: 10.1021/jp502432g.
    [14] LIU L, LIU Y, ZYBIN S V, et al. ReaxFF-lg: Correction of the ReaxFF reactive force field for London dispersion, with applications to the equations of state for energetic materials [J]. The Journal of Physical Chemistry A, 2011, 115(40): 11016–11022. DOI: 10.1021/jp201599t.
    [15] NOMURA K, KALIA R K, NAKANO A, et al. Dynamic transition in the structure of an energetic crystal during chemical reactions at shock front prior to detonation [J]. Physical review letters, 2007, 99(14): 148303. DOI: 10.1103/PhysRevLett.99.148303.
    [16] BUDZIEN J, THOMPSON A P, ZYBIN S V. Reactive molecular dynamics simulations of shock through a single crystal of pentaerythritoltetranitrate [J]. The Journal of Physical Chemistry B, 2009, 113(40): 13142–13151. DOI: 10.1021/jp9016695.
    [17] LI Y, KALIA R K, NAKANO A, et al. Multistage reaction pathways in detonating high explosives [J]. Applied Physics Letters, 2014, 105(20): 204103. DOI: 10.1063/1.4902128.
    [18] HE L, SEWELL T D, THOMPSON D L. Molecular dynamics simulations of shock waves in oriented nitromethane single crystals [J]. The Journal of chemical physics, 2011, 134(12): 124506. DOI: 10.1063/1.3561397.
    [19] REED E J, FRIED L E, HENSHAW W D, et al. Analysis of simulation technique for steady shock waves in materials with analytical equations of state [J]. Physical Review E, 2006, 74(5): 056706. DOI: 10.1103/PhysRevE.74.056706.
    [20] REED E J, MAITI A, FRIED L E. Anomalous sound propagation and slow kinetics in dynamically compressed amorphous carbon [J]. Physical Review E, 2010, 81(1): 016607. DOI: 10.1103/PhysRevE.81.016607.
    [21] ZHANG L, ZYBIN S V, VAN DUIN A C T, et al. Modeling high rate impact sensitivity of perfect RDX and HMX crystals by ReaxFF reactive dynamics [J]. Journal of Energetic Materials, 2010, 28(S1): 92–127.
    [22] PLIMPTON S. Fast parallel algorithms for short-range molecular dynamics [J]. Journal of computational physics, 1995, 117(1): 1–19. DOI: 10.1006/jcph.1995.1039.
    [23] GUMP J C, PEIRIS S M. Phase transitions and isothermal equations of state of epsilon hexanitrohexaazaisowurtzitane (CL-20) [J]. Journal of Applied Physics, 2008, 104(8): 083509. DOI: 10.1063/1.2990066.
    [24] BRUNDAGE A L. EOS development and numerical modeling of CL-20 compaction [J]. AIP Conference Proceedings, 2009, 1195(1): 1365−1368.
    [25] MARSH S P. LASL shock Hugoniot data [M]. California: University of California Press, 1980: 595.
    [26] 刘海, 李毅, 马兆侠, 等. 定常冲击波作用下六硝基六氮杂异伍兹烷(CL-20)/奥克托今(HMX) 含能共晶初始分解机理研究 [J]. 物理化学学报, 2019, 35(8): 858–867. DOI: 10.3866/PKU.WHXB201812011.
    LIU H, LI Y, MA Z X, et al. Study on the initial decomposition mechanism of energetic co-crystal 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaiso-wurtzitane (CL-20) /1,3,5,7-tetranitro-1,3,5,7-tetrazacy-clooctane (HMX) under a steady shock wave [J]. Acta Physico-Chimica Sinica, 2019, 35(8): 858–867. DOI: 10.3866/PKU.WHXB201812011.
    [27] 李维新. 一维不定常流与冲击波[M]. 2版. 北京: 国防工业出版社, 2003: 212−215.
    [28] ZHANG L, ZYBIN S V, VAN DUIN A C T, et al. Carbon cluster formation during thermal decomposition of octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine and 1, 3, 5-triamino-2, 4, 6-trinitrobenzene high explosives from ReaxFF reactive molecular dynamics simulations [J]. The Journal of Physical Chemistry A, 2009, 113(40): 10619–10640. DOI: 10.1021/jp901353a.
    [29] FURMAN D, KOSLOFF R, DUBNIKOVA F, et al. Decomposition of condensed phase energetic materials: Interplay between uni-and bimolecular mechanisms [J]. Journal of the American Chemical Society, 2014, 136(11): 4192–4200. DOI: 10.1021/ja410020f.
    [30] ZHANG X Q, CHEN X R, KALIAMURTHI S, et al. Initial decomposition of the co-crystal of CL-20/TNT: sensitivity decrease under shock loading [J]. The Journal of Physical Chemistry C, 2018, 122(42): 24270–24278. DOI: 10.1021/acs.jpcc.8b06953.
  • [1] 姚成宝李若田宙郭永辉 . 空气自由场中强爆炸冲击波传播二维数值模拟. 爆炸与冲击, 2015, 35(4): 585-590. doi: 10.11883/1001-1455(2015)04-0585-06
    [2] 张柱晋艳娟 . 反向起爆模型下的冲击波加载. 爆炸与冲击, 2014, 34(2): 223-228. doi: 10.11883/1001-1455(2014)02-0223-06
    [3] 师华强宗智贾敬蓓 . 水下爆炸冲击波的近场特性. 爆炸与冲击, 2009, 29(2): 125-130. doi: 10.11883/1001-1455(2009)02-0125-06
    [4] 周杰陶钢王健 . 爆炸冲击波对肺损伤的数值模拟. 爆炸与冲击, 2012, 32(4): 418-422. doi: 10.11883/1001-1455(2012)04-0418-05
    [5] 刘贵兵侯海量朱锡张国栋 . 液滴对爆炸冲击波的衰减作用. 爆炸与冲击, 2017, 37(5): 844-852. doi: 10.11883/1001-1455(2017)05-0844-09
    [6] 周佩杰王坚陶钢周杰 . 泡沫材料对冲击波的衰减特性. 爆炸与冲击, 2015, 35(5): 675-681. doi: 10.11883/1001-1455(2015)05-0675-07
    [7] 李梅蒋建伟王昕 . 复合装药空气中爆炸冲击波传播特性. 爆炸与冲击, 2018, 38(2): 367-372. doi: 10.11883/bzycj-2016-0209
    [8] 许志宇谭永华李小明 . 冲击波的小波数值计算方法. 爆炸与冲击, 2020, 40(1): 014201-1-014201-9. doi: 10.11883/bzycj-2018-0467
    [9] 侯俊亮蒋建伟门建兵王树有 . 预制孔靶板在爆炸冲击波载荷作用下的动态响应. 爆炸与冲击, 2013, 33(6): 662-666. doi: 10.11883/1001-1455(2013)06-0662-05
    [10] 侯日立周平彭建祥 . 冲击波作用下LY12铝合金结构毁伤的数值模拟. 爆炸与冲击, 2012, 32(5): 470-474. doi: 10.11883/1001-1455(2012)05-0470-05
    [11] 顾文彬郑向平刘建青李丹俊陆鸣 . 浅层水中爆炸冲击波对混凝土墩斜碰撞作用试验研究. 爆炸与冲击, 2006, 26(4): 361-366. doi: 10.11883/1001-1455(2006)04-0361-06
    [12] 余尚江李科杰 . 混凝土结构内冲击波应力传感器设计及其行为. 爆炸与冲击, 2005, 25(4): 350-354. doi: 10.11883/1001-1455(2005)04-0350-05
    [13] 赵新颖王伯良李席 . 温压炸药在野外近地空爆中的冲击波规律. 爆炸与冲击, 2016, 36(1): 38-42. doi: 10.11883/1001-1455(2016)01-0038-05
    [14] 陈文张庆明 , . 爆炸冲击波下导弹结构动力分析模型初步研究. 爆炸与冲击, 2009, 29(2): 199-204. doi: 10.11883/1001-1455(2009)02-0199-06
    [15] 彭常贤 , . 声振耦合分析在冲击波测量中的应用. 爆炸与冲击, 2008, 28(5): 427-432. doi: 10.11883/1001-1455(2008)05-0427-06
    [16] 庞伟宾李永池何翔 . 化爆冲击波在T型通道内到时规律的实验研究. 爆炸与冲击, 2007, 27(1): 63-67. doi: 10.11883/1001-1455(2007)01-0063-05
    [17] 杜红棉何志文马铁华 . 冲击波超压测试系统二次仪表频域特性. 爆炸与冲击, 2015, 35(2): 261-266. doi: 10.11883/1001-1455(2015)02-0261-06
    [18] 颜事龙王尹军 . 冲击波作用下乳化炸药压力减敏的表征方法. 爆炸与冲击, 2006, 26(5): 441-447. doi: 10.11883/1001-1455(2006)05-0441-07
    [19] 郭亚丽韩焱王黎明 . 基于广义逆算法的冲击波超压场重建方法. 爆炸与冲击, 2014, 34(6): 764-768. doi: 10.11883/1001-1455(2014)06-0764-05
    [20] 李永池姚磊董杰朱林法胡秀章 . 粘塑性变截面杆中冲击波的演化. 爆炸与冲击, 2005, 25(6): 481-486. doi: 10.11883/1001-1455(2005)06-0481-06
  • 加载中
图(6)表(1)
计量
  • 文章访问数:  13
  • HTML全文浏览量:  69
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-06
  • 录用日期:  2019-11-27
  • 网络出版日期:  2020-02-25
  • 刊出日期:  2020-03-01

CL-20/HMX共晶冲击起爆模拟

    作者简介:刘 海(1985- ),男,博士,助理研究员,liumy2016@163.com
  • 中国空气动力研究与发展中心超高速碰撞研究中心,四川 绵阳 621000

摘要: 采用非平衡分子动力学方法模拟了CL-20/HMX共晶炸药冲击压缩和化学反应行为,获得了密度以及粒子速度的时空分布、冲击雨贡纽、冲击起爆压力、爆轰压力等数据,以及主要中间产物和稳定产物分布。模拟结果显示,共晶的初始分解路径是CL-20中N—NO2键断裂,主要稳定产物是N2、CO2和H2O。CL-20和HMX的分解速率随着冲击波速度的增加而增加,并且逐渐接近,但各冲击条件下CL-20分子的衰减速率均大于HMX。

English Abstract

  • CL-20是当前C/H/O/N类含能材料中能量密度最高(约2.04 g/cm3)的炸药,氧平衡接近−11%,爆速和爆压分别为9.5 km/s和43 GPa,但因较高的感度限制了在武器装备中的应用。2012年,Bolton等[1]成功制备了物质的量之比为2∶1的CL-20/HMX共晶,评估结果显示其爆速为9.484 km/s,优于β-HMX和γ-CL-20,而感度低于多种晶型的CL-20单晶。并且,CL-20/HMX共晶燃烧时特征信号较低,因此在推进剂领域具有较大的应用潜力。2014年,Sun等[2]计算得到CL-20/HMX共晶的多项物性参数,并解释了共晶感度低于ε-CL-20单晶的原因。随后,Liu等[3]分析了CL-20/HMX在0~100 GPa静水压条件下的物理、力学特性。高温热分解模拟显示CL-20/HMX共晶及其单晶组分的初始分解路径均是N—N键断裂。并且,热分解过程中CL-20反应放热将传递至HMX,加速HMX的分解[4]。Doblas等[5]采用喷雾闪蒸技术制备了纳米尺度的CL-20/HMX共晶,并研究了其结构特征以及热行为。对于单晶,单分子CL-20初始分解产物是稳定的芳香化合物1,5-二氢二咪唑[4,5-b:4′,5′-e]吡嗪[6]。Isayev等[7]则给出了不同意见,其研究结果表明单分子CL-20唯一初始反应路径是N—NO2键断裂。高温(2 000~3 500 K)条件下,CL-20热分解的初始反应路径是N—NO2键断裂,最终反应产物是N2、CO2和H2O[8-9]。在冲击条件下,当冲击波速度为8、9 km/s时,CL-20单晶的初始分解路径是N—NO2键断裂;当速度提升至10、11 km/s时,H原子脱落成为初始分解路径。N2、H2O、CO2为主要稳定产物[10]。当冲击波速度为8~10 km/s时,HMX冲击分解的主要初始产物是NO2,最终产物主要是N2、H2O、CO2、CO、H2、NH3[11]。而当冲击波速度提升至11 km/s时,初始反应路径变化为C—H键断裂[12-13]

    当前针对CL-20/HMX共晶冲击行为的研究相对较少。直接冲击作用下CL-20/HMX共晶冲击雨贡纽、爆炸参数、初始反应路径和稳定反应产物等物理、化学特征对于理解含能共晶起爆和爆轰特性具有现实意义。本文中采用非平衡分子动力学方法,同时结合在含能材料热及冲击行为研究中得到广泛应用的ReaxFF/lg势[14],模拟CL-20/HMX共晶的冲击压缩过程,旨在获得不同冲击条件下含能材料中的冲击波传播时空特征和冲击诱导化学反应行为,揭示冲击起爆和爆炸参数,并与组份单晶开展对比分析,为共晶材料设计与工程应用提供数据支撑。

    • 物质的量之比为2∶1的CL-20/HMX单晶胞中含有4个CL-20分子、2个HMX分子[1]。构建20×4×4的超晶胞,空间尺寸为32.69 nm×3.98 nm×4.86 nm。其中CL-20和HMX分子数量分别为1 280和640,共计64 000个原子。首先通过不断迭代调整原子坐标以使系统能量最小化;随后采用等温等压NPT系综将系统缓慢升至常温常压,此过程持续5 ps;接着将体系维持在常温常压下15 ps。温度和压力迭代的阻尼系数分别为300 fs和1 000 fs,整个平衡过程采用周期边界条件,平衡后共晶的密度约为1.9 g/cm3

      原子分子系统的冲击模拟中,冲击加载方法主要有非平衡加载方法(动量镜[15-17]和活塞[18]),多尺度冲击技术(MSST)[19-20]以及分阶段压缩-膨胀方法[21]等。目前,主要使用前两类方法。活塞和动量镜加载方法的关键是在近真实碰撞条件下,针对较大模型体系,从波后粒子速度up出发,获得冲击波速度us以及沿空间路径的密度、温度等冲击波传播时空特征。而MSST方法的核心是在冲击波关系式的约束下,针对较小模型体系施加冲击波条件us,获得沿时间路径的波后粒子速度up、密度、温度等动力学和热力学量变化。本文中,采用动量镜方法进行冲击加载。具体为CL-20/HMX共晶块体材料以初始速度(粒子速度)up=−0.5~−4.0 km/s穿过动量镜,动量镜反射原子动量以激发冲击波以速度us沿晶格矢量a方向(x+方向)传播。沿冲击波传播方向采用收缩边界条件,垂直于冲击波传播方向采用周期边界条件。整个模拟过程采用的时间步长均为0.1 fs。CL-20/HMX单胞、超晶胞模型以及冲击波产生和传播原理如图1所示。本文中,采用大尺度原子分子并行模拟器(large-scale atomic/molecular massivelyparallel simulator,LAMMPS)[22]开展计算。

      图  1  CL-20/HMX模型以及动量镜激发冲击波传播原理图

      Figure 1.  CL-20/HMX structure model and schematic diagram of the shock wave propagation induced by momentum mirror

    • 沿冲击波传播方向对共晶材料进行一维分层,每层的厚度为0.5 nm,约含有1 000个原子。对不同时刻每层中的密度和粒子速度进行统计分析形成x-t图,如图2所示。整体上,x-t图主要分为三个区,Ⅰ区表示未压缩材料的初始状态,Ⅱ区表示冲击压缩区,Ⅰ、Ⅱ区界面表示冲击波阵面轨迹,而Ⅲ区表示冲击波到达右端自由面后反射形成的稀疏拉伸区。随着up的增大,冲击波到达右端自由面的时间逐渐缩短。当up=0.5 km/s,冲击压缩区密度最大为2.5 g/cm3。而当up提升至2.5 km/s和4 km/s时,冲击压缩区密度最大值分别为3.5 g/cm3和4 g/cm3。对于粒子速度,冲击波到达右端自由面后,反射形成的稀疏波进入含能材料并拉伸粒子向右运动,并且在拉伸初期,右端自由面粒子速度约为2up,随后粒子加速运动。另外,当up=0.5 km/s时,Ⅰ和Ⅱ区界面有变形,这表明低速情况下,冲击波传播速度不稳定,而在高速情况下,冲击波传播速度较为稳定。

      图  2  冲击波传播引发粒子速度和密度变化的x-t

      Figure 2.  x-t diagrams of particle velocity and density induced by shock wave propagation

    • 冲击波在凝聚介质中传播时,在相当宽的速度范围内,冲击波速度us和波后粒子速度up之间具有如下形式的线性关系:us=sup+c0,其中s表示绝热体积常数,c0表示材料声速。依据冲击波到达右端自由面的时间以及含能块体的长度得到冲击波速度,并绘制CL-20/HMX的冲击波速度-粒子速度关系如图3(a)所示。空心矩形是本模拟结果,绿色线段是依据试验获得的材料声速(2.77 km/s)和绝热体积常数(2.19)绘制[23-24]。而实心矩形、空心圆和实心圆则分别表示前期获得的CL-20模拟结果[10]、HMX冲击雨贡纽试验数据[25]和模拟获得的HMX数据[11]

      图  3  冲击波速度-粒子速度关系及p-V/V0雨贡纽

      Figure 3.  Shock wave velocity-particle velocity relation and p-V/V0 Hugoniot

      从图中可以看出,当us<9.38 km/s(对应up=3 km/s)时,冲击波速度-粒子速度基本呈线性关系,而当us>9.38 km/s时,冲击波速度-粒子速度发生偏折,此时含能材料爆轰后冲击波速度-粒子速度的间断特征与CL-20/HMX共晶多尺度模拟结果[26]和动量镜冲击加载获得的RDX模拟结果一致[15]。图中黑色线表示Bolton等通过Cheetah模拟获得CL-20/HMX共晶爆速D=9.484 km/s[1]。含能材料us-up关系在发生爆轰后出现向下偏折,原因在于爆速是定值,即随着up的增加,us会略高于爆速,处于过驱爆轰状态,但随着传播而逐渐衰减为定常爆轰。对共晶未爆轰时的冲击波速度-粒子速度数据进行拟合得到us=1.724up+4.215。另外,同冲击波速度条件下,CL-20/HMX共晶波后粒子速度小于多尺度模拟结果[26],但整体上介于HMX和CL-20单晶之间。但与实验结果有一定的偏离。可能的原因是模拟对象是完美晶体,而真实含能材料多含有孔洞、杂质、颗粒边界条件等缺陷。图3(b)表示p-V/V0雨贡纽数据,实心矩形表示本模拟结果,空心圆和实心圆则分别表示前期获得的CL-20模拟结果[10]。模拟结果显示同压力条件下,CL-20的压缩程度高于CL-20/HMX共晶。依据冲击波关系式以及冲击波速度-粒子速度关系,得到冲击压缩条件下凝聚介质p-V/V0雨贡纽的一般形式[27]

      $p = \frac{{{c_0^2}{\rho _0}(1 - V/{V_0})}}{{{{(1 - s(1 - V/{V_0}))}^2}}}$

      式中:ρ0为材料初始密度。将拟合冲击波速度-粒子速度关系得到的c0s代入式(1)得到:

      $p = \frac{{33.76(1 - V/{V_0})}}{{{{(1 - 1.724(1 - V/{V_0}))}^2}}}$

      依据上述关系绘制理论解,如图3(b)中黑色曲线所示。模拟结果显示,理论解与数值模拟结果基本一致。而当共晶发生爆轰后,其p-V/V0数据发生偏离。据此计算得到CL-20/HMX爆压约为45.6 GPa,介于多尺度模拟结果获得的36.75~47.43 GPa之间[26]

    • 图4表示共晶中CL-20和HMX分子数量的衰减曲线。整体上,随着up的增加,共晶中CL-20和HMX分子数量的衰减速率加快,并且含能分子冲击分解前存在诱导期。当up=0.5,1 km/s时,共晶并未发生分解。当up提升至1.5 km/s时,共晶中仅CL-20发生分解。up≥2 km/s后,共晶中CL-20和HMX均发生分解,且分子数量基本呈线性衰减。在冲击波到达右端自由面后,反应物分子数量基本不变。这表明冲击波反射形成的稀疏波拉伸材料不会导致反应物的分解。

      图  4  共晶中CL-20和HMX分子衰减曲线

      Figure 4.  Decay curves of CL-20 and HMX molecules in co-crystal

      图5(a)中,左y轴表示CL-20和HMX诱导期τ,该诱导期表示冲击加载后至激发释能化学反应的时间[28-29],右y轴表示对冲击波到达右端自由面时,共晶中CL-20和HMX分子数量衰减曲线进行归一化后得到的平均衰减速率(decay rate)rd。整体上CL-20和HMX诱导期均随着up的增加而缩短。在各冲击条件下,HMX的诱导时间均多于CL-20,且随着冲击速度的增加,CL-20和HMX的诱导时间逐渐接近。通过对两组数据进行比较,获得各冲击条件下CL-20和HMX衰减速率的接近程度,如表1所示。CL-20和HMX分子数量的衰减速率随着up的增加而增加,并且各冲击速度条件下,共晶中CL-20分子数量的衰减速率均大于HMX,但随着冲击波速度的增加,CL-20和HMX的衰减速率逐渐接近。这主要是由于含能材料爆轰后,CL-20反应放热传递至HMX,增加了HMX的反应速率,同时使得CL-20分解速率降低,进而导致CL-20和HMX分子数量的衰减速率逐渐接近。在CL-20/HMX热分解[4]与CL-20/TNT[30]冲击诱导分解中均发现热传递现象。

      up /(km·s−1)τCL-20/psτHMX/psrd, CL-20rd, HMXrd, CL-20/rd, HMX
      1.50.800 02.100 00.004 00.001 13.636 4
      2.00.254 00.800 00.016 00.002 36.956 5
      2.50.100 00.400 00.062 00.031 02.000 0
      3.00.052 00.300 00.160 00.114 01.403 5
      3.50.049 00.160 00.234 00.205 01.141 5
      4.00.001 00.150 00.289 00.278 01.039 6

      表 1  各冲击条件下共晶中CL-20和HMX诱导时间和衰减速率比较

      Table 1.  Comparison of induction time and decay rate of CL-20 and HMX in co-crystal under different shock conditions

      图  5  诱导时间和衰减速率及冲击到达右端自由面时主要产物的数量统计

      Figure 5.  Induction time and decay rate under different shock conditions and quantitative statistics of the main products when the shock wave reaches the free surface

      图5(b)表示冲击波到达右端自由面时,共晶分解产生碎片的种类及其数量分布。随着up的增加,反应产物的类别和数量均逐渐增多。当up=1.5 km/s时,共晶开始分解并产生NO2,对比图4可知,此NO2来自于CL-20分子中N-NO2键断裂。随着冲击波速度的增加,系统中产物类别逐渐增多。up=2.5 km/s,除NO2外,其他中间小分子产物NO、OH、HONO出现,并且最终稳定产物N2形成。当up增至4 km/s时,系统中小分子产物NO2、NO、OH、HONO和H中,NO2和OH的数量较多。而稳定产物为N2、CO2和H2O。并且,N2分子数量较多,H2O分子数量居其次,CO2数量最少。另外,NO2是共晶在各冲击条件下分解的主要中间产物。CL-20和HMX单晶模拟结果显示,较高冲击速度条件下,N—NO2键断裂受到抑制,H原子首先脱落[6-7]。对比多尺度模拟结果[26]发现,二者在初始反应路径、中间小分子产物类别和分布特征上均一致,但在最终反应产物上,多尺度模拟结果显示最终稳定产物为N2、CO2、H2O、CO和H2。这主要是由于本模拟仅针对冲击压缩阶段开展分析,而多尺度模拟结果涵盖了冲击波扫过共晶材料后的压缩阶段以及部分化学反应区(膨胀阶段)中的稳定产物。因此可以推测本模拟中,压缩波到达右端自由面反射拉伸材料,体积膨胀后,压缩阶段形成的大分子团簇分解,CO和H2分子将出现。在Li[17]的模拟结果中,证明了RDX冲击压缩后拉伸膨胀阶段的新产物形成情况。

      up=1.5 km/s (us=6.8 km/s) 时,冲击诱导共晶分解。采用兰金-雨贡纽关系P=ρ0usup,其中usup分别表示冲击波激发共晶分解的冲击波速度和对应的粒子速度,计算得到冲击起爆压力为19.38 GPa,与多尺度模拟方法获得的结果(16.52 GPa)[26]基本一致,低于HMX单晶的冲击起爆压力35.1 GPa[11]

      图6直观展示了up=1.5 km/s和4.0 km/s时冲击波从左向右传播的序列图像。同一时刻,up=4.0 km/s对应的冲击波传播距离大于up=1.5 km/s的情况。当up=4.0 km/s时,在t=3.15 ps时刻,含能材料完全压缩。而后反射形成的稀疏波入射至材料中,自右向左传播并拉伸已压缩材料,如图中t=3.5,4.0 ps对应的图像所示。而当up=1.5 km/s,在4 ps内,冲击波未到达右端自由面。另外,通过压缩区的疏密程度可以直观反映出up=4.0 km/s对应的冲击压缩区的密度大于up=1.5 km/s时的压缩程度。

      图  6  冲击波传播的序列图像

      Figure 6.  Sequence images of shock wave propagation

    • 采用非平衡分子动力学方法,结合ReaxFF/lg势模拟了CL-20/HMX共晶的冲击压缩行为,分析了冲击波传播时空特征、冲击雨贡纽关系、反应物分解速率以及产物分布等,得到以下主要结论。

      (1)低速情况下,CL-20/HMX共晶中冲击波传播不稳定。随着up的增大,冲击波传播趋于稳定。当usD时,CL-20/HMX共晶中冲击波速度-粒子速度呈线性关系,usD时,冲击波速度-粒子速度出现间断。计算得到爆速约为9.38 km/s,而爆压则为45.6 GPa。

      (2)共晶中HMX冲击响应稍滞后于CL-20,随着冲击波速度的增加,二者的诱导时间逐渐接近。CL-20和HMX分子衰减速率随着冲击波速度的增加而增加,并且逐渐接近,但各冲击条件下CL-20分子的衰减速率均大于HMX。

      (3)CL-20/HMX共晶冲击压缩分解的主要中间产物和稳定产物分别是NO2和N2、CO2、H2O。冲击激发共晶起爆的冲击波速度为6.8 km/s,对应的冲击起爆压力为19.38 GPa。

参考文献 (30)

目录

    /

    返回文章
    返回