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

高温下带金属壳PBX炸药低速撞击敏感性数值模拟

胡偲 吴艳青 黄风雷

引用本文:
Citation:

高温下带金属壳PBX炸药低速撞击敏感性数值模拟

    作者简介: 胡 偲(1994- ),女,博士,3120160068@bit.edu.cn;
    通讯作者: 吴艳青, wuyqing@bit.edu.cn
  • 中图分类号: O383

Numerical simulation of confined PBX charge under low velocity impact at high temperature

    Corresponding author: WU Yanqing, wuyqing@bit.edu.cn ;
  • CLC number: O383

  • 摘要: 炸药撞击感度和热安全性是评价炸药安全性能的重要指标。为了对高温下炸药撞击敏感性变化规律进行可靠预测,本文中通过数值模拟,研究不同预加热温度下带壳PBX炸药装药在小弹丸低速撞击下的热力学响应,得到炸药点火前至点火阶段局部高温区的位置、形态、温度和应变随时间在炸药中分布的变化。结果显示,炸药发生点火的撞击阈值速度与烤燃温度的关系并非单一随温度升高而降低,而是在加热至348.15 K时达到最高;根据温度和应力应变云图分析可得,随着烤燃温度的提高,炸药强度下降,PBX炸药装药局部高温区快速升温的主导因素由局部剪切变为压缩。热软化对炸药的撞击敏感性起重要作用。
  • 图 1  使用表1参数的应变率-温度-最大流动应力曲面图

    Figure 1.  Maximum flow stress as a function of temperature and strain rate at selected parameters

    图 2  二维轴对称有限元模型及其网格划分示意图

    Figure 2.  Meshing of two-dimensional axisymmetric finite-element model

    图 3  预加热温度为301.15 K时,不同撞击速度下炸药装药的反应过程分数云图

    Figure 3.  Reaction fractions at different impact velocities and a preheating temperature of 301.15 K

    图 4  预加热温度为348.15 K时,不同撞击速度下炸药装药的反应过程分数云图

    Figure 4.  Reaction fractions at different impact velocities and a preheating temperature of 348.15 K

    图 5  预加热温度为378.15 K时,不同撞击速度下炸药装药的反应过程分数云图

    Figure 5.  Reaction fractions at different impact velocities and a preheating temperature of 378.15 K

    图 6  不同预加热温度和撞击速度下炸药装药中的最大反应过程分数

    Figure 6.  The maximum reaction process fraction at different pre-heating temperatures and impact velocities

    图 7  不同预加热温度和撞击速度下PBX炸药装药的平均反应过程分数

    Figure 7.  The average reaction fractions of the PBX charge at different temperatures and impact velocities

    图 8  预加热温度为301.15 K时炸药装药的温度云图

    Figure 8.  Temperature of the explosive charge at a preheating temperature of 301.15 K

    图 9  预加热温度为301.15 K时炸药装药的体积应变云图

    Figure 9.  Volumetric strain of the explosive charge at a preheating temperature of 301.15 K

    图 10  预加热温度为301.15 K炸药装药的等效应变云图

    Figure 10.  Equivalent strain of the explosive charge at a preheating temperature of 301.15 K

    图 11  预加热温度为301.15 K炸药装药的Mises应力云图

    Figure 11.  Mises stress of the explosive charge at a preheating temperature of 301.15 K

    图 12  预加热温度为348.15 K时炸药装药的温度云图

    Figure 12.  Temperature of the explosive charge at a preheating temperature of 348.15 K

    图 13  预加热温度为348.15 K时炸药装药的体积应变云图

    Figure 13.  Volume strain of the explosive charge at a preheating temperature of 348.15 K

    图 14  预加热温度为348.15 K炸药装药的等效应变云图

    Figure 14.  Equivalent strain of the explosive charge at a preheating temperature of 348.15 K

    图 15  预加热温度为348.15 K炸药装药的Mises应力云图

    Figure 15.  Mises stress of the explosive charge at a preheating temperature of 348.15 K

    图 16  预加热温度为378.15 K时炸药装药的温度云图

    Figure 16.  Temperature of the explosive charge at a preheating temperature of 378.15 K

    图 17  预加热温度为378.15 K时炸药装药的体积应变云图

    Figure 17.  Volumetric strain of the explosive charge at a preheating temperature of 378.15 K

    图 18  预加热温度为378.15 K炸药装药的等效应变云图

    Figure 18.  Equivalent strain of the explosive charge at a preheating temperature of 378.15 K

    图 19  预加热温度为378.15 K炸药装药的Mises应力云图

    Figure 19.  Mises stress of the explosive charge at a preheating temperature of 378.15 K

    图 20  局部高温区的位置变化

    Figure 20.  The change of the locations of localized heating at different preheating temperatures

    图 21  不同预加热温度下PBX炸药中局部高温区的温度时间曲线

    Figure 21.  Temperature as a function of time for the localized heating regions at different preheating temperatures

    图 22  不同预加热温度下PBX炸药中局部高温区的塑性功时间曲线

    Figure 22.  Plastic work as a function of time for the localized heating regions at different preheating temperatures

    图 23  不同预加热温度下PBX炸药中局部高温区内能时程曲线

    Figure 23.  Histories of internal energy for the localized heating regions at different preheating temperatures

    表 1  Johnson–Cook强度模型参数

    Table 1.  Parameters of Johnson–Cook model

    A/MPa B/MPa C m n Tref/K Tm/K
    15 0 0.2 0.6 1 301.15 540
    下载: 导出CSV

    表 2  固体反应物的状态方程参数

    Table 2.  Equation of the state parameters for solid reactants

    pcrush/MPa μcrush K1 K2 K3 Klock μlock
    0.689 5 0.125×10−3 0.102 7×1011 0.217 2×1012 0.288 8×1013 0.665 4×1011 0.09
    下载: 导出CSV

    表 3  气体产物的状态方程参数

    Table 3.  Equation of the state parameters for gaseous products

    R1/GPa R2 R3/GPa R4 ω
    166.89 5.9 5.969 2.1 0.45
    下载: 导出CSV

    表 4  PBX化学动力学参数

    Table 4.  Chemical kinetic parameters for PBX

    No. Zk/s−1 Ea k/(kJ·g−1·mol·K)
    1 1.08×1026 148.93
    2 1.41×1021 220.64
    3 2.60×1016 185.48
    下载: 导出CSV

    表 5  用于温度计算的参数

    Table 5.  Parameters for temperature calculation

    TCJ/K vCJ/cm3/g R5/GPa CV, s/(J·kg−1·K−1) CV, CJ/(J·kg−1·K−1)
    2 500 0.387 3 0.85 1 600 2 000
    下载: 导出CSV
  • [1] KIM S, BARUA A, HORIE Y, et al. Prediction of probabilistic ignition behavior of PBXs from microstructural stochasticity[C]// Biennial International Conference of the APS Topical Group on Shock Compression of Condensed Matter, 2013.
    [2] PERRY W L, ZUCKER J, DICKSON P M, et al. Interplay of explosive thermal reaction dynamics and structural confinement [J]. Journal of Applied Physics, 2007, 101(7): 074901–5. DOI: 10.1063/1.2713090.
    [3] ASANTE D O, KIM S, CHAE J, et al. CFD cook-off simulation and thermal decomposition of confined high energetic material [J]. Propellants Explosives Pyrotechnics, 2015, 40(5): 699–705. DOI: 10.1002/prep.201400296.
    [4] CHIDESTER S K, GREEN L G, LEE C G. A frictional work predictive method for the initiation of solid high explosives from low-pressure impacts[R]. California, the United States: Lawrence Livermore National Laboratory, 1993.
    [5] BERGHOUT H L, SON S F, SKIDMORE C B, et al. Combustion of damaged PBX 9501 explosive [J]. Thermochimica Acta, 2002, 384(1-2): 261–77. DOI: 10.1016/S0040-6031(01)00802-4.
    [6] 张文英, 邓晓雅, 陈思瑾, 等. PBX9501热感度、热膨胀及力学性能的分子动力学模拟 [J]. 火炸药学报, 2016, 39(3): 32–6. DOI: 10.14077/j.issn.1007-7812.2016.03.006
    ZHANG Wenying, DENG Xiaoya, CHEN Sijin, et al. Molecular dynamics simulation of thermal sensitivity, thermal expansion and mechanical properties of PBX9501 [J]. Chinese Journal of Explosives and Propellants, 2016, 39(3): 32–6. DOI: 10.14077/j.issn.1007-7812.2016.03.006
    [7] SHUKLA P, DEEPU M. Experimental and numerical investigations of thermal ignition of a phase changing energetic material [J]. Defence Science Journal, 2016, 66(3): 228–35. DOI: 10.14429/dsj.66.9357.
    [8] DICKSON P, ASAY B, HENSON B, et al. Thermal cook-off response of confined PBX 9501[C]//Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 2004.
    [9] BURNHAM K, WEESE R K, WEMHOFF A P, et al. A historical and current perspective on predicting thermal cookoff behavior [J]. Journal of Thermal Analysis and Calorimetry, 2007, 89(2): 407–15. DOI: 10.1007/s10973-006-8161-6.
    [10] GUSHANOV A R, VOLODINA N A, BELOV G V, et al. Numerical simulation of experiments on the low- velocity impact on HMX- based HE using explosive transformation initiation kinetics [J]. Journal of Energetic Materials, 2010, 28(S1): 50–65. DOI: 10.1080/07370651003720429.
    [11] DRODGE D R, WILLIAMSON D M. Understanding damage in polymer-bonded explosive composites [J]. Journal of Materials Science, 2016, 51(2): 668–79. DOI: 10.1007/s10853-013-7378-6.
    [12] WILLEY T M, LAUDERBACH L, GAGLIARDI F, et al. Mesoscale evolution of voids and microstructural changes in HMX-based explosives during heating through the beta-delta phase transition [J]. Journal of Applied Physics, 2015, 118(5): 055901. DOI: 10.1063/1.4927614.
    [13] YAN Q L, ZEMAN S, SANCHEZ-JIMENEZ P E, et al. The effect of polymer matrices on the thermal hazard properties of RDX-based PBXs by using model-free and combined kinetic analysis [J]. Journal of Hazardous Materials, 2014, 271(1): 85–95. DOI: 10.1016/j.jhazmat.2014.02:019.
    [14] FORBES J, TARVER C, URTIEW P, et al. The effects of confinement and temperature on the shock sensitivity of solid explosives[C]//The 11th International Detonation Symposium. Snowmass, Colorado, 1998.
    [15] SANDUSKY H W, ANDREWS S A, ALEXANDER K E. Hazard tests on a heated TATB-based high explosive[C]//The 11th International Detonation Symposium. Snowmass, Colorado, 1998.
    [16] DAI X G, WEN Y S, HUANG H, et al. Impact response characteristics of a cyclotetramethylene tetranitramine based polymer-bonded explosives under different temperatures [J]. Journal of Applied Physics, 2013, 114(11): 114906. DOI: 10.1063/1.4820248.
    [17] GRUAU C, PICART D, BELMAS R, et al. Ignition of a confined high explosive under low velocity impact [J]. International Journal of Impact Engineering, 2009, 36(4): 537–50. DOI: 10.1016/j.ijimpeng.2008.08.002.
    [18] APPAN A S, RENLUND A M, STACHOWIAK J C, et al. Mechanical properties determination by real-time ultrasonic characterization of thermally damaged energetic materials[C]//The 12th International Detonation Symposium. 2002.
    [19] DAVANDE H, BEDROV D, SMITH G D. Thermodynamic, transport, and viscoelastic properties of PBX-9501 binder: a molecular dynamics simulations study[J]. Journal of Energetic Materials, 2008, 26(2): 115−38. DOI:10.1080/07370650701801937.
    [20] HOLMQUIST T J, JOHNSON G R. Determination of constants and comparison of results for various constitutive models [J]. Journal De Physique Iii, 1991, 1(C3): 853–60. DOI: 10.1051/Jp4:19913119.
    [21] GRAY III G, IDAR D, BLUMENTHAL W, et al. High-and low-strain rate compression properties of several energetic material composites as a function of strain rate and temperature: LA-UR-98-3059[R]. Los Alamos National Laboratory, 1998.
    [22] PARKER G R, PETERSON P D, ASAY B W, et al. Examination of morphological changes that affect gas permeation through thermally damaged explosives [J]. Propellants Explosives Pyrotechnics, 2004, 29(5): 274–81. DOI: 10.1002/prep.200400057.
    [23] PETERSON P D, MANG J T, ASAY B W. Quantitative analysis of damage in an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazonic-based composite explosive subjected to a linear thermal gradient [J]. Journal of Applied Physics, 2005, 97(9): 7. DOI: 10.1063/1.1879072.
    [24] WANG X J, WU Y Q, HUANG F L. Thermal-mechanical-chemical responses of polymer-bonded explosives using a mesoscopic reactive model under impact loading [J]. Journal of Hazardous Materials, 2017, 321(2): 56–67. DOI: 10.1016/j.jhazmat.2016.08.061.
    [25] YOH J J, MCCLELLAND M A, MAIENSCHEIN J L, et al. Simulating thermal explosion of octahydrotetranitrotetrazine-based explosives: Model comparison with experiment [J]. Journal of Applied Physics, 2006, 100(7): 9. DOI: 10.1063/1.2357418.
    [26] TARVER C M, CHIDESTER S K. On the violence of high explosive reactions [J]. Journal of Pressure Vessel Technology-Transactions of the ASME, 2005, 127(1): 39–48. DOI: 10.1115/1.1845474.
    [27] HANDLEY C A. Numerical modelling of two HMX-based plastic-bonded explosives at the mesoscale[D]. University of St Andrews, 2011,
    [28] MENIKOFF R, SEWELL T D. Constituent properties of HMX needed for mesoscale simulations [J]. Combustion Theory and Modelling, 2002, 6(1): 103–25. DOI: 10.1088/1364-7830/6/1/306.
    [29] HORIE Y. Hot spots, high explosives ignition, and material microstructure[C]//Materials Science Forum, 2014.
    [30] DIENES J K, ZUO Q H, KERSHNER J D. Impact initiation of explosives and propellants via statistical crack mechanics [J]. Journal of the Mechanics & Physics of Solids, 2006, 54(6): 1237–75. DOI: 10.1016/j.jmps.2005.12.001.
    [31] MENIKOFF R. Detonation waves in PBX 9501 [J]. Combustion Theory and Modelling, 2006, 10(6): 1003–21. DOI: 10.1080/13647830600851754.
    [32] HENSON B, ASAY B W, DICKSON P M, et al. Measurement of explosion time as a function of impact pressure for PBX 9501 subject to weak planar shocks[R]. Office of Scientific & Technical Information Technical Reports, 1997.
  • [1] 吴飞鹏徐尔斯刘静尉雪梅蒲春生任扬 . 组合脉冲压裂加载过程耦合模拟及火药配比影响敏感性分析. 爆炸与冲击, 2018, 38(3): 683-687. doi: 10.11883/bzycj-2016-0302
    [2] 马丽英李向东周兰伟蓝肖颖宫小泽姚志军 . 高速破片撞击充液容器时容器壁面的损伤. 爆炸与冲击, 2019, 39(2): 023302-1-023302-12. doi: 10.11883/bzycj-2018-0009
    [3] 张晓敏彭向和张培源 . 热力藕合问题的交替算法. 爆炸与冲击, 2006, 26(4): 309-314. doi: 10.11883/1001-1455(2006)04-0309-06
    [4] 李俊玲傅华谭多望卢芳云 . PBX 炸药的拉伸断裂损伤分析. 爆炸与冲击, 2011, 31(6): 624-629. doi: 10.11883/1001-1455(2011)06-0624-06
    [5] 傅华李俊玲谭多望 . PBX 炸药本构关系的实验研究. 爆炸与冲击, 2012, 32(3): 231-236. doi: 10.11883/1001-1455(2012)03-0231-06
    [6] 王晨陈朗刘群皮铮迪胡晓棉 . 多组分PBX 炸药细观结构冲击点火数值模拟. 爆炸与冲击, 2014, 34(2): 167-173. doi: 10.11883/1001-1455(2014)02-0167-07
    [7] 张永康李玉龙汤忠斌杨洪许海 . 冰雹撞击下泡沫铝夹芯板的动态响应. 爆炸与冲击, 2018, 38(2): 373-380. doi: 10.11883/bzycj-2016-0232
    [8] 吴桂英赵勇刚赵隆茂杨桂通 . 圆板动力反直观行为实验研究. 爆炸与冲击, 2006, 26(2): 121-124. doi: 10.11883/1001-1455(2006)02-0121-04
    [9] 路亮龙源郭涛谢全民赵长啸高福银 . 爆破地震波作用下城市隧道结构动力响应的敏感性. 爆炸与冲击, 2014, 34(6): 701-708. doi: 10.11883/1001-1455(2014)06-0701-08
    [10] 于津平高耀林范钦文王荪源 . 多孔颗粒火药床燃烧转爆轰的敏感性研究. 爆炸与冲击, 1991, 11(2): 139-145.
    [11] 史飞飞索涛侯兵李玉龙 . YB-2航空有机玻璃的应变率和温度敏感性及其本构模型. 爆炸与冲击, 2015, 35(6): 769-776. doi: 10.11883/1001-1455(2015)06-0769-08
    [12] 松全才魏化震 . PBX颗粒力学性质对其撞击感度影响的初步研究. 爆炸与冲击, 1990, 10(3): 244-249.
    [13] 吴艳青黄风雷艾德友 . HMX 颗粒炸药低速撞击点火实验研究. 爆炸与冲击, 2011, 31(6): 592-599. doi: 10.11883/1001-1455(2011)06-0592-08
    [14] 吴艳青鲍小伟王明扬黄风雷张柱 . RDX/HMX颗粒炸药落锤撞击点火-燃烧机理. 爆炸与冲击, 2017, 37(2): 339-346. doi: 10.11883/1001-1455(2017)02-0339-08
    [15] 赵锋孙承纬文尚刚赵剑衡龙新平 . 飞片撞击下JO-9159炸药的脆性损伤. 爆炸与冲击, 2001, 21(2): 121-125.
    [16] 张小宁徐更光王廷增 . 射流撞击粉碎法制备超细炸药颗粒原理分析. 爆炸与冲击, 1999, 19(3): 262-268.
    [17] 方青卫玉章张克明张寿齐赵玉华 . 射弹倾斜撞击带盖板炸药引发爆轰的条件. 爆炸与冲击, 1997, 17(2): 153-158.
    [18] 袁方强蔡从中赵帅 . 用结构参数预测硝基类炸药的撞击感度. 爆炸与冲击, 2013, 33(1): 79-84. doi: 10.11883/1001-1455(2013)01-0079-06
    [19] 章冠人 . 关于分形和炸药撞击感度的一点记注. 爆炸与冲击, 1994, 14(2): 105-111.
    [20] 卢强王占江刘晓新郭志昀吴玉蛟 . 薄片炸药与固体靶冲量耦合的计算模型. 爆炸与冲击, 2017, 37(1): 84-91. doi: 10.11883/1001-1455(2017)01-0084-08
  • 加载中
图(23)表(5)
计量
  • 文章访问数:  171
  • HTML全文浏览量:  774
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-07-10
  • 录用日期:  2017-12-21
  • 刊出日期:  2019-04-01

高温下带金属壳PBX炸药低速撞击敏感性数值模拟

    作者简介:胡 偲(1994- ),女,博士,3120160068@bit.edu.cn
    通讯作者: 吴艳青, wuyqing@bit.edu.cn
  • 北京理工大学爆炸科学与技术国家重点实验室,北京 100081

摘要: 炸药撞击感度和热安全性是评价炸药安全性能的重要指标。为了对高温下炸药撞击敏感性变化规律进行可靠预测,本文中通过数值模拟,研究不同预加热温度下带壳PBX炸药装药在小弹丸低速撞击下的热力学响应,得到炸药点火前至点火阶段局部高温区的位置、形态、温度和应变随时间在炸药中分布的变化。结果显示,炸药发生点火的撞击阈值速度与烤燃温度的关系并非单一随温度升高而降低,而是在加热至348.15 K时达到最高;根据温度和应力应变云图分析可得,随着烤燃温度的提高,炸药强度下降,PBX炸药装药局部高温区快速升温的主导因素由局部剪切变为压缩。热软化对炸药的撞击敏感性起重要作用。

English Abstract

  • 在炸药的存储、运输和装配过程中,常常出现炸药暴露或曾经暴露在高温条件下的部分被约束在容器中时受到低速撞击的复合场景,可能导致炸药被意外点燃,进而转化成为更剧烈的爆燃或爆轰,引发多种事故[1-3]。长久以来,炸药的撞击安全性和热安全性一般是分开研究。Steven撞击实验被广泛用于研究炸药的撞击敏感性[4],确定带金属壳约束的炸药受钢制子弹低速撞击点火的临界速度[5]。张文英等[6]对不同厚度的高能炸药进行改良的Steven实验,并提出力学因素导致的反应可能存在通用的机理。Shukla等[7]研究了不同形状的子弹撞击炸药装药导致装药压碎、穿刺和贯穿情况,并提出了对应上述情况的几种点火机制,其中包括摩擦、剪切和应变。在炸药的热安全性方面,前人设计了多种烤燃试验用以观察炸药装药的响应并预测炸药在高温下的行为[8-9],为炸药热点火的研究提供有用信息[10]。Steven实验和烤燃试验都可为建立数值模拟的初始数据库提供便利[7]。目前这两方面的研究已经阐明,炸药的内部结构可能在撞击或高温作用下发生变化[11],热损伤和力学损伤在微观和宏观尺度下同时存在,因此在实际情景即热力复合的条件下,需要考虑热力耦合效应。

    PBX高能炸药在军事方面的应用使得其特性广受关注,它的主要成分为炸药晶体和粘合剂。炸药晶体的相变[12]、吸热或放热粘结剂的影响[13]以及这两种成分之间的巨大差异使得PBX炸药的损伤机制非常复杂,而实验研究中也体现出了热损伤与撞击敏感性之间关系的复杂性。 Dickson等[8]提出,热点火期间PBX的热学性能与力学性能相互紧密耦合。Forbes等[14]测试了LX-04-01和LX-17的安全性行为,在170及250℃的高温下对两种炸药进行撞击试验,结果显示加热处理后的两种炸药的安全性均有下降倾向。Sandusky等[15]研究了PBX-9502不同温度下的撞击安全性,并与未加热的情景相互对照,结果显示,加热至240℃炸药的点火阈值速度略微降低。上述试验和计算结果均表明:加热后炸药的冲击起爆距离和发生点火临界速度降低,即高温下炸药的感度提高。Dai等[16]研究了加热至75、105和160 ℃的PBX-2炸药(HMX/TATB/olefin)的撞击敏感性,并使用扫描电镜(SEM)检验加热后炸药的热损伤情况。分析发现,在某些特定加热温度下的炸药撞击安全性将得到改善,而不是单一性地降低。

    与实验研究相比,PBX炸药热力耦合加载数值模拟研究,尤其是在数值预测与实验结果的比较方面直到现在都非常欠缺。烤燃温度对PBX敏感性的影响还有待进一步深入的研究。根据文献中的数据,PBX的模量随着烤燃温度的升高明显下降,这表明烤燃温度可能会影响受撞击的PBX的变形过程,并进一步影响撞击点火的临界速度。为研究温度对PBX炸药撞击敏感性的作用,本文中建立有限元模型模拟HMX基PBX炸药在不同预加热温度下受子弹低速撞击时的热力学性能变化,分析温度对炸药本构参数的影响。

    • Dai等[16]的实验中,PBX炸药材料的主要成分为HMX、TATB和粘结剂。模型偏应力张量由一个弹塑性模型描述,偏应力-应变关系描述如下:

      ${S_{{{ij}}}} = 2G\left( {{e_{{{ij}}}} - e_{{{ij}}}^{\rm{p}}} \right)$

      式中:${S_{ {ij}}}$为偏应力, ${e_{ {ij}}}$为偏应变,$e_{{{ij}}}^{\rm{p}}$为偏应变中的塑性部分,上标${\rm{p}}$代表塑性,下标${{i, j = 1, 2, 3}}$$G$为材料的剪切模量。若一点的von-Mises应力$\overline \sigma = \sqrt {\displaystyle\frac{3}{2}{S_{{{ij}}}}{S_{{{ij}}}}} $大于压缩强度${\sigma _{\rm{y}}}$, 则认为该点进入塑性状态。偏应力${S_{ {ij}}}$通过满足von-Mises准则抗压强度${\sigma _{\rm {y}}} = \overline \sigma $确定。

      随着温度的升高,PBX的抗拉强度下降。当加热到348.15 K以上时,PBX的拉伸强度低于0.5 MPa,而当温度超过80℃时,拉伸强度接近零[17]。此时,随着粘合剂黏度的降低主导PBX炸药的力学响应[18-19],PBX强度达最低,表现出黏性流动的性质[16]。考虑温度和应变率效应的影响,选择含有应变率项和温度项的Johnson-Cook强度模型[20],抗压强度模型表示为:

      ${\sigma _{\rm{y}}} = \left({A + B\varepsilon _{\rm{p}}^n} \right)\left({1 + C\ln {\varepsilon ^{\rm{*}}}} \right)\left({1-{T^{{\rm{*}}m}}} \right)$

      其中,${\varepsilon _{\rm{p}}}$为等效塑性应变,$ \varepsilon^{\rm{*}} =\dot \varepsilon /{\dot \varepsilon _0}$${\dot \varepsilon _0} = 1.0\;{{\rm{s}}^{ - 1}}$时的无量纲应变率,${T^*} = \left( {T - {T_{{\rm{ref}}}}} \right)/\left( {{T_{\rm{m}}} - {T_{{\rm{ref}}}}} \right)$为相对温度,${T_{{\rm{ref}}}}$为环境温度,${T_{\rm{m}}}$代表材料熔化温度,该模型仅在$0 \text{≤} {T^*} \text{≤} 1.0$时成立。$A\text{、} {\rm{ }}B\text{、} {\rm{ }}C\text{、} {\rm{ }}m\text{、} {\rm{ }}n$为材料参数,$A$为材料屈服极限,$B$为加工硬化模量,$C$为应变率常数,$n$为硬化系数,$m$为热软化常数。

      硬化系数$n$参考与PBX粘结剂成分力学性能相似的高聚物取1,环境温度取为室温301.15 K,材料熔化温度可通过PBX-2的热失重分析曲线(TG曲线)得到[16];其余参数由一系列压缩实验得到[21]。Johnson-Cook强度模型中所使用的材料参数总结至表1图1为使用这些参数时的应变率-温度-最大流应力曲面图。

      A/MPa B/MPa C m n Tref/K Tm/K
      15 0 0.2 0.6 1 301.15 540

      表 1  Johnson–Cook强度模型参数

      Table 1.  Parameters of Johnson–Cook model

      图  1  使用表1参数的应变率-温度-最大流动应力曲面图

      Figure 1.  Maximum flow stress as a function of temperature and strain rate at selected parameters

      对PBX进行热处理可能使其孔隙率发生明显变化[22-23],因此,选择孔隙率相关的压力-体积关系描述加热后固体炸药的压缩行为。固体PBX炸药的压力为:

      $ \left\{ \begin{array}{l} {p_{\rm{s}}} = \left( {\displaystyle\frac{p_{{\rm{crush}}}}{\mu _{{\rm{crush}}}}} \right)\mu \quad\quad\quad\;\;\quad\quad\quad\quad\mu \text{≤} {\mu _{{\rm{crush}}}}\\ {p_{\rm{s}}} = {p_{{\rm{crush}}}} + {K_1}\bar \mu + {K_2}{{\bar \mu }^2} + {K_3}{{\bar \mu }^3}\quad\!{\rm{ }}{\mu _{{\rm{crush}}}} \text{<} \mu \text{≤} {\mu _{{\rm{lock}}}}\\ {p_{\rm{s}}} = {K_{{\rm{lock}}}}\left( {\mu - {\mu _{{\rm{lock}}}}} \right)\quad\quad\quad\;\quad\quad\mu \text{>} {\mu _{{\rm{lock}}}} \end{array} \right.$

      式中:$\mu = {\rho \mathord{\left/ {\vphantom {\rho {{\rho _0} - 1}}} \right. \kern-\nulldelimiterspace} {{\rho _0} - 1}}$, $\bar \mu = \mu - {\mu _{{\rm{crush}}}}$, ${K_1}\text{、} {\rm{ }}{K_2}\text{、} {\rm{ }}{K_3}\text{、} {\rm{ }}{p_{{\rm{crush}}}}\text{、}{\rm{ }}{\mu _{{\rm{crsuh}}}}\text{、}{\rm{ }}{K_{{\rm{lock}}}}\text{、} {\rm{ }}{\mu _{{\rm{lock}}}}$为材料参数,下标s表示固体。参数值取自参考文献[24]. 压力模型参数列于表2。该方程即固体反应物的状态方程。

      pcrush/MPa μcrush K1 K2 K3 Klock μlock
      0.689 5 0.125×10−3 0.102 7×1011 0.217 2×1012 0.288 8×1013 0.665 4×1011 0.09

      表 2  固体反应物的状态方程参数

      Table 2.  Equation of the state parameters for solid reactants

    • 反应流模型的建立需要2个状态方程和一个反应速率定律[18],使用Jones–Wilkins–Lee状态方程来描述反应产物的压力、体积和温度之间的关系。反应气体产物的压力表示为:

      ${p_{\rm{g}}} = {R_1}\left( {{\rm{1}} - \frac{\omega }{{{R_2}\overline V}}} \right)\exp \left( { - {R_2}\overline V} \right) + {R_3}\left( {1 - \frac{\omega }{{{R_4}\overline V}}} \right)\exp \left( { - {R_4}\overline V} \right) + \frac{{\omega E}}{{\overline V}}$

      式中:$\overline V = {{{V_{\rm{g}}}} \mathord{\left/ {\vphantom {{{V_{\rm{g}}}} {{V_0}}}} \right. \kern-\nulldelimiterspace} {{V_0}}}$为相对体积,${V_{\rm g}}$为气体体积,下标g表示气体,${V_0}$为初始体积,$E$为初始体积内能,${R_1}\text{、} {\rm{ }}{R_2}\text{、} {\rm{ }}{R_3}$${R_4}$ 为材料参数,$\omega $是Grüneisen常数,根据参考文献[19]选择参数值,所使用的JWL的EOS参数列于表3

      R1/GPa R2 R3/GPa R4 ω
      166.89 5.9 5.969 2.1 0.45

      表 3  气体产物的状态方程参数

      Table 3.  Equation of the state parameters for gaseous products

    • 对于化学反应动力学,考虑到在反应的早期阶段主要是低程度热分解,只有较少部分的炸药被消耗。使用温度相关的三步Arrhenius模型来描述该时期的化学反应动力学:

      ${r_{{k}}} = {Z_{{k}}}\exp \left( { - \frac{{{E_{{{{\rm a}k}}}}}}{{R{T_{\rm{a}}}}}} \right)\left( {1 - {F_{{k}}}} \right), \;\quad k = 1, 2, 3$

      式中:${r_{{k}}}$为反应速率,${Z_{{k}}}$为频率因子,${E_{{{{\rm a}k}}}}$为反应k的活化能;$R$为理想气体常数,${T_{\rm{a}}}$为绝对温度,$F$为反应过程分数,$F = 0$表示炸药尚未反应,$F = 1.0$表示炸药反应完全。参考值的选择主要根据参考文献[25],根据Dai等[16]的实验进行了调整以提高结果的一致性,取值如表4所示。

      No. Zk/s−1 Ea k/(kJ·g−1·mol·K)
      1 1.08×1026 148.93
      2 1.41×1021 220.64
      3 2.60×1016 185.48

      表 4  PBX化学动力学参数

      Table 4.  Chemical kinetic parameters for PBX

    • 由于化学反应发展过程受到局部反应影响[26],局部温度对于形成点火核心起到重要作用。固体反应物的温度计算根据下式:

      ${T_{\rm{s}}} = {T_0} + \frac{{\Delta E}}{{{C_{V{\rm{, s}}}}\rho }}$

      式中:${T_0}$为初始温度,$\Delta E$为单位初始体积的内能增量,${C_{V{\rm{, s}}}}$为固体的比热容,$\rho $为密度。

      气体产物的温度可以使用文献[27]中给出的方程计算:

      ${T_{\rm{g}}} = {\left( {\frac{{{\upsilon _{{\rm{CJ}}}}}}{{{\upsilon _{\rm{g}}}}}} \right)^\omega }\left[{{T_{{\rm{CJ}}}} + \frac{1}{{{C_{V{\rm{, CJ}}}}}}{{\left( {\frac{{{\upsilon _{\rm{g}}}}}{{{\upsilon _{{\rm{CJ}}}}}}} \right)}^\omega }\left( {{e_{\rm{g}}}-{e_{\rm{i}}}\left( {{\upsilon _{\rm{g}}}} \right)} \right)} \right]$

      ${e}\left( {{\upsilon _{\rm{g}}}} \right) = \frac{{{\upsilon _0}}}{{{R_3}}}{R_1}\exp \left( { - {R_3}\frac{{{\upsilon _{\rm{g}}}}}{{{\upsilon _0}}}} \right) + \frac{{{\upsilon _0}}}{{{R_4}}}{R_2}\exp \left( { - {R_4}\frac{{{\upsilon _{\rm{g}}}}}{{{\upsilon _0}}}} \right) + \frac{{{R_5}{\upsilon _0}}}{\omega }{\left( {\frac{{{\upsilon _0}}}{{{\upsilon _{\rm{g}}}}}} \right)^\omega }$

      式中:${C_{V{\rm{, CJ}}}}$为CJ状态下的比热容,${\upsilon _{{\rm{CJ}}}}$为CJ状态下的比容,${T_{{\rm{CJ}}}}$为CJ状态的气体产物的温度。参数选取文献[24]的计算结果,列于表5

      TCJ/K vCJ/cm3/g R5/GPa CV, s/(J·kg−1·K−1) CV, CJ/(J·kg−1·K−1)
      2 500 0.387 3 0.85 1 600 2 000

      表 5  用于温度计算的参数

      Table 5.  Parameters for temperature calculation

    • 在使用的反应模型中所采用的混合法则假设未反应炸药和气体反应产物的压力和温度同时达平衡[27],即:

      ${p_{\rm{s}}}\left( {{\upsilon _{\rm{s}}}, {e_{\rm{s}}}} \right) = {p_{\rm{g}}}\left( {{\upsilon _{\rm{g}}}, {e_{\rm{g}}}} \right)$

      ${T_{\rm{s}}}\left( {{\upsilon _{\rm{s}}}, {e_{\rm{s}}}} \right) = {T_{\rm{g}}}\left( {{\upsilon _{\rm{g}}}, {e_{\rm{g}}}} \right)$

      混合物的比内能和比体积为:

      $e = {e_{\rm{s}}}\left( {1 - \lambda } \right) + {e_{\rm{g}}}\lambda$

      $\upsilon = {\upsilon _{\rm{s}}}\left( {1 - \lambda } \right) + {\upsilon _{\rm{g}}}\lambda$

      式中:$e, {\rm{ }}\upsilon $分别为混合物的比内能和比体积,下标${\rm{s\text{、} g}}$分别对应表示固体和气体。对${\upsilon _{\rm{s}}}$${e_{\rm{s}}}$使用Newton-Raphson迭代法来解式(9)、(10)中的平衡。

    • 图2为本文中使用的二维轴对称模型。炸药尺寸为$\varnothing ∅ $50 mm×20 mm,钢盖板厚度为3.5 mm,约束炸药的钢环厚度为5 mm,这两部分与钢制炸药基底合并成为包覆炸药的外壳,小弹丸尺寸为$\varnothing ∅ $20 mm×20 mm。计算采用四边形网格,模型节点总数为1 549,单元总数为1 370,底部加约束。顶端小弹丸以不同初始速度撞击炸药装药,撞击前炸药的初始温度均一,作为预加热温度。由于点火发生在几乎不到1 μs的时间尺度上,热传导只会影响相当小的距离,因此可忽略热传导作用[28, 29],模型绝热。60 μs内炸药将完成加载和卸载,因此模拟时长设置为60 μs。炸药钢制外壳和小弹丸的材料参数根据对低碳钢进行静态拉伸实验确定。

      图  2  二维轴对称有限元模型及其网格划分示意图

      Figure 2.  Meshing of two-dimensional axisymmetric finite-element model

    • 使用二维数值计算对HMX基PBX炸药的烤燃后撞击的热力复合试验[16]进行模拟,通过与实验结果对比,可将模拟结果中炸药的反应过程分数与实验的反应程度对应起来,判断是否发生点火,并讨论不同预加热温度对炸药反应程度的影响。

      图3为常温下(301.15 K)时不同撞击速度下炸药装药在60 μs时的反应过程分数云图。撞击速度为235 m/s时,60 μs时炸药装药中最大反应过程分数低于0.8,对比文献[16]中撞击速度为237 m/s时的实验结果,可以看出约235 m/s的撞击速度下炸药未发生点火。而当撞击速度提高到265 m/s,最大反应过程分数由0.7左右提高到0.9,与实验[16]中撞击速度269.9 m/s时PBX炸药装药发生点火且反应消耗了所有炸药相对应。同时可以注意到,该温度下反应过程分数最高的区域为子弹撞击边缘,该区域为可能的点火起始的位置。

      图  3  预加热温度为301.15 K时,不同撞击速度下炸药装药的反应过程分数云图

      Figure 3.  Reaction fractions at different impact velocities and a preheating temperature of 301.15 K

      预加热温度为348.15 K时,模拟结果显示在290 m/s的撞击速度下,60 μs时炸药中最大反应过程分数不超过0.7,大部分区域的反应过程分数低于0.5(图4);当撞击速度为360 m/s时,最大反应过程分数超过0.9,大部分区域的反应过程分数都大于0.5。对应的实验结果[16]显示,撞击速度为276 m/s和294 m/s时,炸药上表面出现压痕,但未发生明显反应,而撞击速度为367 m/s时发生点火,所有炸药都被反应消耗。

      图  4  预加热温度为348.15 K时,不同撞击速度下炸药装药的反应过程分数云图

      Figure 4.  Reaction fractions at different impact velocities and a preheating temperature of 348.15 K

      预加热温度为378.15 K时,模拟结果显示炸药在撞击速度为285 m/s时最高反应过程分数小于0.8,在撞击速度为295 m/s时最高反应过程分数大于0.9,如图5所示。根据实验[16],PBX炸药在撞击速度小于286 m/s时被压缩,但未被消耗,说明未发生点火,而在撞击速度为298.3 m/s和315.4 m/s时发生燃烧反应。不同于常温时,该温度下反应过程分数最高区域位于炸药和子弹的中轴线,靠近炸药的上表面。关于可能的点火起始位置变化将在后文中作进一步讨论。

      图  5  预加热温度为378.15 K时,不同撞击速度下炸药装药的反应过程分数云图

      Figure 5.  Reaction fractions at different impact velocities and a preheating temperature of 378.15 K

      根据上述对比可见,当炸药装药所有单元反应过程分数的最大值在60 μs内达到0.9时,局部反应可能发生失控,发展成为更高级别反应,发生点火。由此通过反应过程分数的最大值可以确定临界撞击点火速度。常温(301.15 K)时,临界撞击点火速度在260~265 m/s之间(图6(a));预加热温度为348.15 K时,临界撞击点火速度在330~360 m/s之间(图6(b));预加热温度为378.15 K时,临界撞击点火速度在285~355 m/s之间(图6(c))。模拟结果显示临界撞击点火速度不随预加热温度单调降低,而是先升高至360 m/s,此时预加热温度为348.15 K,再随温度升高而降低。也即说明炸药装药在预加热温度约为348.15 K时撞击感度最低。模拟结果与其对应的实验结果[16]相符合。

      图  6  不同预加热温度和撞击速度下炸药装药中的最大反应过程分数

      Figure 6.  The maximum reaction process fraction at different pre-heating temperatures and impact velocities

    • 为判断炸药反应的程度,计算不同预加热温度和撞击速度下PBX炸药装药的平均反应过程分数如图7所示。通过与实验结果进行对比可见,当平均反应过程分数到约0.15时,炸药装药的反应程度高。炸药平均反应过程分数随预加热温度升高而加速上升的趋势说明,在较高预加热温度下,PBX炸药装药中的反应更易传播。可以将平均反应过程分数曲线分为3个阶段,对应炸药到点火的反应过程中的3个阶段:(I)缓慢的低速热分解;(II)形成局部反应区域;(III)局部反应区域内的化学反应放热,加速这一区域的分解反应,若反应程度升高足够快,炸药装药将发生点火。

      图  7  不同预加热温度和撞击速度下PBX炸药装药的平均反应过程分数

      Figure 7.  The average reaction fractions of the PBX charge at different temperatures and impact velocities

    • 常温(301.15 K)下,撞击速度为265 m/s时,温度云图(图8)显示在55 μs左右,在子弹边缘下的炸药装药上表面出现了直径20 mm的环形局部高温区,图中为二维轴对称模型的右半部分。同时根据图911可见,这一区域体积应变和等效塑性应变最大,Mises应力最小。可以推断,撞击产生的剪切集中处的应力超过屈服强度,导致局部升温。机械能使变形能转化为大量热能,炸药材料受热点火,产生局部燃烧反应。变形可能导致热软化,材料越软,变形越大,生热越多。这一自馈过程可能形成很小的局部区域,并产生极高的局部温度[30-31]。总的反应过程释放能量,并加速分解反应。若反应释放率足够大,燃烧反应将失控成为爆轰[32]

      图  8  预加热温度为301.15 K时炸药装药的温度云图

      Figure 8.  Temperature of the explosive charge at a preheating temperature of 301.15 K

      图  9  预加热温度为301.15 K时炸药装药的体积应变云图

      Figure 9.  Volumetric strain of the explosive charge at a preheating temperature of 301.15 K

      图  10  预加热温度为301.15 K炸药装药的等效应变云图

      Figure 10.  Equivalent strain of the explosive charge at a preheating temperature of 301.15 K

      图  11  预加热温度为301.15 K炸药装药的Mises应力云图

      Figure 11.  Mises stress of the explosive charge at a preheating temperature of 301.15 K

      预加热温度为348.15 K、撞击速度为360 m/s时,约在55 μs时炸药装药中出现2个局部高温区。其一为在子弹边缘下的炸药装药上表面的环形局部高温区,直径19.30 mm;另一个是出现在炸药中轴线上的点状局部高温区,距离炸药上表面约8.5 mm(图12)。根据体积应变和等效塑性应变云图(图1314)分析可得,环形局部高温区主要由塑性剪切效应导致,而点状局部高温区由压缩功做功得到,同时两高温区位置处的Mises应力最小(图15)。

      图  12  预加热温度为348.15 K时炸药装药的温度云图

      Figure 12.  Temperature of the explosive charge at a preheating temperature of 348.15 K

      图  13  预加热温度为348.15 K时炸药装药的体积应变云图

      Figure 13.  Volume strain of the explosive charge at a preheating temperature of 348.15 K

      图  14  预加热温度为348.15 K炸药装药的等效应变云图

      Figure 14.  Equivalent strain of the explosive charge at a preheating temperature of 348.15 K

      图  15  预加热温度为348.15 K炸药装药的Mises应力云图

      Figure 15.  Mises stress of the explosive charge at a preheating temperature of 348.15 K

      预加热温度378.15 K,撞击速度295 m/s时,58 μs时,盘状局部高温区出现在中轴线上,距离炸药装药上表面约3 mm(图16)。如图17~19所示,Mises应力较小的区域主要集中在炸药装药中部,等效应变集中在子弹边缘下的炸药装药上表面,而体积应变集中区域位于炸药装药中轴线上,形状大小与温度云图中局部高温区一致,说明压缩而非剪切作用在此时对局部高温区的产生和升温起主要作用。

      图  16  预加热温度为378.15 K时炸药装药的温度云图

      Figure 16.  Temperature of the explosive charge at a preheating temperature of 378.15 K

      图  17  预加热温度为378.15 K时炸药装药的体积应变云图

      Figure 17.  Volumetric strain of the explosive charge at a preheating temperature of 378.15 K

      图  18  预加热温度为378.15 K炸药装药的等效应变云图

      Figure 18.  Equivalent strain of the explosive charge at a preheating temperature of 378.15 K

      图  19  预加热温度为378.15 K炸药装药的Mises应力云图

      Figure 19.  Mises stress of the explosive charge at a preheating temperature of 378.15 K

      可以观察到,局部高温区的位置随预加热温度升高,由炸药表面剪切集中处(图20(a))转向炸药中部受压缩最大处(图20(b),(c))。图21为不同预加热温度下局部高温区的温度时间曲线。相比于常温,高温下炸药温度突增趋于平缓,这可能是由于炸药的热软化引起的。常温下发生点火的炸药局部高温区最先进入快速升温阶段,而在预加热温度为348.15 K时,其进入快速升温阶段的时间最晚,同时其临界撞击点火速度阈值最高,塑性功做功最大(图22),说明此时炸药撞击敏感度最低。根据图23,预加热温度越高,同一时刻最高温度区内能越大,点火发生的反应将更加剧烈。

      图  20  局部高温区的位置变化

      Figure 20.  The change of the locations of localized heating at different preheating temperatures

      图  21  不同预加热温度下PBX炸药中局部高温区的温度时间曲线

      Figure 21.  Temperature as a function of time for the localized heating regions at different preheating temperatures

      图  22  不同预加热温度下PBX炸药中局部高温区的塑性功时间曲线

      Figure 22.  Plastic work as a function of time for the localized heating regions at different preheating temperatures

      图  23  不同预加热温度下PBX炸药中局部高温区内能时程曲线

      Figure 23.  Histories of internal energy for the localized heating regions at different preheating temperatures

    • 通过对HMX基PBX炸药在不同烤燃温度下低速撞击实验的数值模拟,可以得到以下结论:炸药的发生点火的撞击阈值速度与烤燃温度的关系并非单一随温度升高而降低,而是随预加热温度先升高至预加热温度为348.15 K时达最大,此时阈值速度为360 m/s,后在更高的预加热温度下下降,数值模拟复现了实验中的现象。随着预加热温度升高,局部高温区位置由炸药表面剪切集中处,转向炸药中部主要受压缩处。这样的现象可能由于受热炸药强度降低,压缩对局部高温区升温的作用超过剪切效应而产生。由此可见,热软化对炸药的撞击敏感性起重要作用。

参考文献 (32)

目录

    /

    返回文章
    返回