Calculation on multi-step thermal decomposition of HMX-and TATB-based composite explosive under cook-off conditions
-
摘要: 采用多步热分解反应动力学模型,描述单质炸药热分解反应,提出了多组分网格单元计算方法,对以HMX/TATB为基的多元混合炸药在烤燃条件下的热反应过程进行了计算。通过炸药烤燃实验测量了炸药内部温度,获得了炸药点火时间,验证了计算的准确性。分析了混合炸药组成比例的变化对炸药热反应性能的影响。在HMX/TATB混合炸药热反应阶段,主要是HMX发生分解反应释放热量,TATB的反应量很少,随着混合炸药中TATB含量的增多,炸药的点火时间逐渐增长,点火温度逐渐增高,炸药热安全性增强。Abstract: The thermal decomposition of each component in the composite explosive was described by the multi-step chemical kinetics model.A multi-component grid unit calculation method was put forward to calculate the thermal decomposition reaction processes of the HMX and TATB based composite explosives under the cook-off conditions.And the cook-off tests were carried out to measure the temperatures and ignition times of the explosives and verify the accuracies of the calculations.The thermal reaction performances of the composite explosives were analyzed by changing the mass fraction of each component in them.The results show that the exothermic decomposition of the HMX predominates in the thermal decomposition process of the HMX/TATB composite explosive.Moreover, the ignition time and temperature of the composite explosives increase with the increasing of the TATB contents in the composite explosives.
-
Key words:
- mechanics of explosion /
- multistep reaction /
- multicomponent grid unit /
- composite explosive /
- cook-off
-
炸药装药在热刺激作用下的点火时间、点火温度和反应情况, 可以通过炸药烤燃实验进行观测。而根据炸药的结构组成, 对炸药热反应过程进行预测计算, 一直是人们关心的问题。
早在1928年N.N.Semenov[1]首次提出采用单步反应动力学模型描述炸药放热分解反应, 通过解析计算得到了均温系统中炸药的热爆炸临界温度和绝热爆炸延滞期。但炸药热反应解析计算理论, 由于受计算求解方法的制约, 对复杂几何结构, 非线性加热条件的炸药装药不能求解, 难以在工程实际中广泛应用。20世纪60年代以后, 随着数值计算的进步, 数值计算成为炸药热反应理论分析的主要手段[2-3]。但炸药热反应模型仍采用单步反应模型, 即认为炸药由反应物直接生成最终产物, 忽略了中间产物对反应过程的影响。由于通过烤燃实验标定炸药反应动力学参数后, 采用单步反应模型能够较好地计算炸药装药的点火时间和温度, 因此单步反应模型被大量使用[4-5]。但炸药烤燃过程中的热反应与其组分密切相关, 单步反应模型不能描述炸药中各组分分解反应历程, 无法实现混合炸药热反应过程的预测计算。1981年, R.R.McGuire等[6]提出了单质炸药多步热分解反应动力学模型, 准确计算了HMX、RDX和TNT等炸药的一维热爆炸实验的点火时间。多步热分解反应动力学模型将中间产物考虑到反应过程中, 更细致地描述了炸药热反应过程, 为混合炸药热反应计算奠定了基础。C.M.Tarver等[7-8]建立了交叉反应计算模型, 模型中不仅包含炸药多步分解反应, 还考虑了粘结剂与炸药分解产物发生的交叉反应, 分析了吸热粘结剂和放热粘结剂对PBX炸药点火时间的影响。J.J.Yoh等[9-10]采用多步热分解反应动力学模型描述了RDX和HMX炸药在慢速烤燃实验中的热分解过程, 并利用ALE3D软件对炸药点火后的燃烧爆炸现象进行了模拟计算。W.L.Prerry等[11]对HMX多步热分解反应动力学模型进行改进, 采用过渡态理论描述了HMX发生β相→δ相晶型转变的吸热过程, 改进后的模型能够获得HMX晶型转变过程中的温度平台。W.L.Perry等[12]建立了四步可逆反应动力学模型, 对慢速烤燃实验和摩擦引起的PBX9501炸药点火进行数值模拟。从现有研究情况看, 炸药多步热分解反应动力学模型主要被用于单质炸药或单组元的混合炸药计算, 而对于多组元混合炸药, 如何处理各炸药组分对热反应的影响是需要解决的主要问题。
本文中, 通过烤燃实验获得以HMX和TATB为基的混合炸药的基本热反应特征, 采用多步热分解反应动力学模型描述单质炸药热分解反应, 提出多组分网格单元计算方法, 计算各组分对炸药热量的贡献, 从而实现混合炸药热反应的预测计算, 分析炸药组分变化对其热反应性能的影响。
1. 炸药烤燃实验
图 1是炸药烤燃实验装置简图, 实验装置结构和加热方法与传统烤燃弹实验[2]类似。为了更全面地观测炸药热反应特征, 测量了炸药内部不同位置的温度。烤燃装置主要由炸药、壳体和端盖组成。受试炸药为HMX/TATB/Kel-F混合炸药。烤燃弹壳体是厚4mm的钢壳, 药柱直径为50mm, 长88mm。实验时, 在烤燃弹内部中心横截面处放置3个测温热电偶, 电偶1位于药柱中心, 电偶2位于药柱1/2半径处, 电偶3位于药柱边缘。壳体外侧包裹加热套, 加热速率为0.5K/min。根据热电偶测量的温度随时间的变化曲线, 可以获得炸药点火时间和点火时此位置的炸药温度, 从而获得炸药基本的热反应特征。
2. 混合炸药多步热分解反应计算
混合炸药受热后, 外部环境的热量以热传导的形式通过炸药外部的壳体传递到炸药内部, 炸药中各个组分, 即每种单质炸药和粘结剂, 在受热后均要发生分解反应, 并且组分之间也有可能发生化学反应。每种单质炸药在相同温度条件下都有各自的分解规律, 它们并不是由反应物直接生成产物, 而是要经过复杂的多步反应才能完成。炸药发生自热分解反应后, 将炸药的化学能转变为热能, 随着热量积累, 炸药温度逐渐升高, 同时又加速了自热分解反应, 放出更多的热量, 最终导致炸药点火。对混合炸药烤燃过程进行数值模拟, 关键是要准确描述混合炸药的热分解反应。采用多步热分解反应动力学模型描述各个组分的分解反应。其中, 各个组分在不同热分解阶段时, 既发生吸热反应又发生放热反应, 它们的吸热量和放热量之和即为炸药的总反应热。
HMX热分解历程可用4步化学反应进行描述[5]:
反应 1A⟶B,r1=Z1exp(−E1RT)ρA (1) 反应 2 B⟶C,r2=Z2exp(−E2RT)ρB (2) 反应 3C⟶D,r3=Z3exp(−E3RT)ρC (3) 反应 4D⟶E,r4=Z4exp(−E4RT)ρ2D (4) 式中:rm为质量反应速率, Zm为指前因子, Em为反应的活化能, m为反应序号, m=1, 2, …; ρn为质量浓度, n=A, B, …; R为普适气体常数; A为固体β-HMX, B为固体δ-HMX, C为固体中间产物, D和E分别为中间气体产物和最终气体产物。反应1是HMX发生β相→δ相晶型转变的吸热过程; 反应2是δ-HMX吸收热量生成固体中间产物; 反应3是轻微放热反应, 固体中间产物生成气体中间产物CH2O、N2O、HCN和HNO2; 反应4是二级气相反应, 气体中间产物生成最终气体反应产物, 如CO2、CO、N2和H2O, 同时释放出大量化学能, 最终导致炸药点火。
HMX热分解过程为连续反应过程, 即前一步的生成物就是下一步的反应物, 如此依次连续进行, 直至反应结束。连续反应可表示为:
Ak2⟶Bk2⟶Ck4⟶2Dk5⟶E (5) 各步反应速率方程可表示为:
−dρn dt=kmρxn (6) 式中:km为反应速率常数, x为反应级数, x=1, 2;且
km=Zmexp(−EmRT) (7) 对应每一步反应, 热量生成速率可表示为:
Sm=QmZmexp(−EmRT)ρxn (8) 式中:Qm为反应热。
因此, HMX热分解过程中单位时间内生成的热量为:
S(HMX)=ρAQ1k1+ρBQ2k2+ρCQ3k3+ρ2DQ4k4 (9) TATB热分解历程可用3步化学反应进行描述[5]:
反应 5 F⟶G+H2O,r5=Z5exp(−E5RT)ρF (10) 反应 6G⟶H+ Gases ,r6=Z6exp(−E6RT)ρG (11) H⟶I,r7=Z7exp(−E7RT)ρ2H (12) 式中:F为固体TATB, G和H均为固体中间产物, I为最终气体产物。
TATB的分解反应主要发生在固相中, 反应5和反应6均是吸热反应, H2O分子先分解出来, 生成中间产物。反应7是二级放热反应, 由固体中间产物生成最终气体产物, 释放出大量化学能。计算中, 忽略H2O分子和气体中间产物引起的固体反应物的质量消耗。TATB分解过程中单位时间内生成的热量为:
S(TATB)=ρFQ5k5+ρGQ6k6+ρ2HQ7k7 (13) Kel-F是一种吸热粘结剂, 其分解温度低于HMX和TATB的分解温度, 因此, 在相同温度下, Kel-F的热稳定性最高。Kel-F的分解反应如下[5]:
反应 8 J⟶K,r8=Z8exp(−E8RT)ρJ (14) 式中:J为Kel-F, K为分解产物。反应8是一级吸热化学反应。
同时, Kel-F中含有C-H、C-Cl和C-F键, 这些键易与TATB分解的气体产物如HCN、CO、CO2、HNCO和NO等发生吸热反应[11], 最终导致加入吸热粘结剂Kel-F后混合炸药的点火时间比单质炸药的点火时间长。反应机制如下:
反应 9I+J⟶L,r9=Z9exp(−E9RT)ρIρJ (15) 式中:I为TATB的最终气体产物, L为反应产物。反应9是吸热反应, 对于每一个反应物反应级数均为1。Kel-F分解过程中单位时间内生成的热量为:
S(Kel−F)=ρJQ8k8+ρIρJQ9k9 (16) 298K时, HMX、TATB和Kel-F的物性参数及其热分解反应动力学参数[7]见表 1。
表 1 HMX、TATB和Kel-F的物性参数及其热分解反应动力学参数[7]Table 1. Physical parameter for HMX, TATB and Kel-F as well as their thermal and chemical kinetic parameters[7]材料 ρ/
(kg·m-3)cV/
(J·g-1·K-1)λ/
(W·m-1·K-1)m Z/s-1 E/
(J·mol-1)x Q/
(J·g-1)HMX 1 850 1 004.62 0.535 8 1 7.990 9×1020 203 574 1 +42.0 2 1.413 0×1021 221 340 1 +252.0 3 2.608 1×1016 186 060 1 -558.6 4 1.598 4×1012 143 220 2 -5 615.4 TATB 1 835 1 088.33 0.879 0 5 7.016 7×1023 252 000 1 +210.0 6 8.749 3×1012 176 400 1 +210.0 7 4.356 0×1011 141 960 2 -2 940.0 Kel-F 2 020 1 000.43 0.052 7 8 9.934 9×1017 272 328 1 +5 871.6 9 7.694 8×1023 272 328 1 +2 520.0 根据炸药烤燃实验, 建立炸药烤燃三维数值计算模型。模型中主要考虑炸药、壳体和空气槽, 其中壳体外壁为加热边界, 以替代外部加热器的作用。炸药烤燃过程涉及的能量变化, 主要是外部环境的热量以热传导的形式通过炸药外部的壳体传递到炸药内部, 炸药受热后各组分分别进行自热分解反应。在模型中, 将计算网格单元设为多种单质炸药及粘结剂组成的混合热分解单元, 各组分按自身的多步热分解反应释放或者吸收热量, 单元总热量为每个组分吸收和释放热量之和, 计算思路如图 2所示。
在直角坐标系中, 固体炸药烤燃过程中的能量守恒方程为:
ρcp∂T∂t=λ∂∂x(∂T∂x)+λ∂∂y(∂T∂y)+λ∂∂z(∂T∂z)+S (17) 式(17)左边为控制单元内能的增量, 右边前3项分别为从x、y和z方向进入控制单元的净热量, S为炸药自热反应源项。HMX/TATB/Kel-F炸药热分解过程中单位时间内生成的总热量为:
Sa=ρAQ1k1+ρBQ2k2+ρCQ3k3+ρ2DQ4k4+ρFQ5k5+ρGQ6k6+ρ2HQ7k7+ρJQ8k8+ρIρJQ9k9 (18) 根据单质炸药和粘结剂的热分解反应动力学参数, 采用多组分网格单元计算方法, 即可对不同组分配比的混合炸药的烤燃过程进行计算, 不需要通过实验标定混合炸药的反应动力学参数。
3. 结果与讨论
对HMX、TATB和Kel-F的质量分数w分别为0.38、0.57和0.05的HMX/TATB/Kel-F混合炸药的烤燃实验进行数值计算, 炸药内部3个测点的温度曲线如图 3所示。测点1和2处计算曲线的升温速率略高于实验曲线, 这可能是由于计算中未考虑壳体与炸药接触面的接触热阻, 导致炸药温度的计算值略高于实验值。测点3位于炸药表面, 与钢壳接触, 计算得到的温度曲线与实验曲线基本吻合。
加热速率为0.5K/min的条件下, 炸药点火时间的实验值和计算值分别为24 857和21 641s, 误差为12.9%, 两者基本一致, 表明利用建立的混合炸药多步热分解反应动力学模型可以较好地计算烤燃炸药热反应过程。由于实验中不能直接测量到炸药的点火温度, 炸药的点火温度只能由计算获得, 计算得到炸药点火时点火点的温度为527K。
图 4为计算得到的炸药点火点的温度-时间曲线和热量-时间曲线。当炸药放热量急剧增加时, 炸药最高温度区域的温度发生突变, 此时炸药点火, 对应的温度为炸药点火温度。通过炸药点火区域的温度变化, 可以分析炸药热反应的全过程。
图 5为不同时刻炸药点火区域各组分的质量分数分布。当温度高于414K时, 约1%的β-HMX开始发生晶型转变, 生成δ-HMX; 当温度达到452K时, β-HMX基本反应完毕, δ-HMX含量达到最大值, 同时有少量固体中间产物产生; 炸药点火前, TATB和Kel-F的反应量很少; 当温度达到527K时, 炸药点火; 此种配比的混合炸药点火主要是由HMX分解放热引起的。
图 6给出了不同时刻炸药内部的温度分布。10 000s时, 外部钢壳温度最高, 炸药内部温度最低, 热量由钢壳向炸药内部传递。21 641s时, 炸药发生点火, 点火区位于药柱底部以炸药柱中线为轴线的环状区域。这是由于炸药的侧面和底面与壳体紧密接触直接受热, 而炸药上表面为空气层, 空气的导热系数远小于钢的导热系数, 所以炸药上部的温度低于下部, 点火位置位于炸药的底部。
利用混合炸药热反应计算模型, 分别对组分质量分数不同的6种混合炸药在1K/min加热速率下的烤燃过程进行了数值计算, 计算结果如表 2所示。从表 2中可以看出, 加热速率为1K/min时, 随着混合炸药中HMX含量的逐渐减少、TATB含量的逐渐增加, 炸药的点火时间t逐渐延长, 点火温度T逐渐增高。HMX/TATB/Kel-F炸药烤燃过程中, HMX反应活化能低, 热稳定性差, 首先发生分解反应释放热量; HMX发生快速放热反应4之前, TATB主要吸收热量, 随着混合炸药中TATB含量的增多, 吸收的热量增多, 从而延长了混合炸药的点火时间, 提高了炸药的热安全性。HMX发生快速放热反应4时, 释放的热量远多于热传导所消耗的热量, 导致炸药温度迅速升高, 热分解反应加速, 最终导致炸药点火。炸药点火前, HMX对炸药热量积累的贡献最大。
表 2 不同组分含量混合炸药烤燃计算结果Table 2. Calculation results of composite explosives at various content of component炸药组分 t/s T/K w(HMX):w(Kel-F)=95:5 10 965 508 w(HMX):w(TATB):w(Kel-F)=88:7:5 10 974 510 w(HMX):w(TATB):w(Kel-F)=78:17:5 10 988 510 w(HMX):w(TATB):w(Kel-F)=55:40:5 11 036 519 w(HMX):w(TATB):w(Kel-F)=38:57:5 11 095 523 w(HMX):w(TATB):w(Kel-F)=15:80:5 11 289 552 w(TATB):w(Kel-F)=95:5 15 588 620 图 7给出了不同时刻w(HMX):w(Kel-F)=95:5的HMX/Kel-F混合炸药点火区域各组分的质量分数分布。当温度达到413K时, 约有1%的β-HMX开始发生晶型转变, 生成δ-HMX; 当温度达到457K时, β-HMX基本反应完毕, δ-HMX含量达到最大值, 同时有少量固体中间产物生成; 炸药点火前, Kel-F未反应; 当温度达到509K时, 该HMX/Kel-F混合炸药点火。
图 8给出了不同时刻w(TATB):w(Kel-F)=95:5的TATB/Kel-F混合炸药点火点处各组分的质量分数分布。当温度达到482K时, 少量TATB开始发生分解反应, 生成固体中间产物G; 随着温度的升高, 反应加剧; 当温度达到515K时, 少量固体中间产物H和气体产物生成, 同时Kel-F开始发生吸热反应; 当温度达到627K时, TATB/Kel-F混合炸药点火, 此时TATB和Kel-F反应完毕。Kel-F与TATB的分解产物发生吸热反应, 同时自身分解吸热, 延长了TATB/Kel-F混合炸药的点火时间。
4. 结论
在单质炸药多步热分解反应动力学模型的基础上, 采用多组分网格单元计算方法, 实现了混合炸药热反应过程的计算, 可根据炸药的配方预测炸药的热安全性。
HMX/TATB/Kel-F混合炸药烤燃过程中, 炸药点火前, 主要是HMX发生分解反应释放热量, TATB吸收热量, 随着混合炸药中TATB含量的增多, TATB吸收的热量增多, 炸药的点火时间逐渐延长, 点火温度逐渐升高, 炸药的热安全性增强; Kel-F与TATB的分解产物发生吸热反应会延长炸药的点火时间。
-
表 1 HMX、TATB和Kel-F的物性参数及其热分解反应动力学参数[7]
Table 1. Physical parameter for HMX, TATB and Kel-F as well as their thermal and chemical kinetic parameters[7]
材料 ρ/
(kg·m-3)cV/
(J·g-1·K-1)λ/
(W·m-1·K-1)m Z/s-1 E/
(J·mol-1)x Q/
(J·g-1)HMX 1 850 1 004.62 0.535 8 1 7.990 9×1020 203 574 1 +42.0 2 1.413 0×1021 221 340 1 +252.0 3 2.608 1×1016 186 060 1 -558.6 4 1.598 4×1012 143 220 2 -5 615.4 TATB 1 835 1 088.33 0.879 0 5 7.016 7×1023 252 000 1 +210.0 6 8.749 3×1012 176 400 1 +210.0 7 4.356 0×1011 141 960 2 -2 940.0 Kel-F 2 020 1 000.43 0.052 7 8 9.934 9×1017 272 328 1 +5 871.6 9 7.694 8×1023 272 328 1 +2 520.0 表 2 不同组分含量混合炸药烤燃计算结果
Table 2. Calculation results of composite explosives at various content of component
炸药组分 t/s T/K w(HMX):w(Kel-F)=95:5 10 965 508 w(HMX):w(TATB):w(Kel-F)=88:7:5 10 974 510 w(HMX):w(TATB):w(Kel-F)=78:17:5 10 988 510 w(HMX):w(TATB):w(Kel-F)=55:40:5 11 036 519 w(HMX):w(TATB):w(Kel-F)=38:57:5 11 095 523 w(HMX):w(TATB):w(Kel-F)=15:80:5 11 289 552 w(TATB):w(Kel-F)=95:5 15 588 620 -
[1] Semenov N N. Theories of combustion process[J]. Physics: Z, 1928, 48: 571-582. doi: 10.1007/BF01340021 [2] Jones D A, Parker R P. Heat flow calculations for the small-scale cook-off bomb test[R]. AD-A236829, 1991. [3] Wang Pei, Chen Lang, Wang Xiao-feng, et al. Cook-off test and numerical simulation for explosive heated by fire[J]. Journal of Beijing Institute of Technology, 2009, 18(2): 146-151. [4] Williams M R, Matei M V. The decomposition of some RDX and HMX based materials in the one-dimensional time to explosion apparatus. Part 1: Time to explosion and apparent activation energy[J]. Propellants, Explosives, Pyrotechnics, 2006, 31(6): 435-441. doi: 10.1002/prep.200600058 [5] 陈朗, 马欣, 黄毅民, 等.炸药多点测温烤燃实验和数值模拟[J].兵工学报, 2011, 32(10): 1230-1236.Chen Lang, Ma Xin, Huang Yi-min, et al. Multi-point temperature measuring cook-off test and numerical simulation of explosive[J]. Acta Armamentarii, 2011, 32(10): 1230-1236. [6] McGuire R R, Tarver C M. Chemical decomposition models for the thermal explosion of confined HMX, TATB, RDX and TNT explosives[C]//Proceedings of the 7th International Symposium on Detonation. Anapolis, MD, USA, 1981: 1-8. [7] Tarver C M, Koerner J G. Effects of endothermic binders on times to explosion of HMX-and TATB-based plastic bonded explosives[J]. Journal of Energetic Materials, 2008, 26(1): 1-28. [8] Tarver C M. Effects of exothermic binders on times to explosion of HMX-plastic bonded explosives[C]//Proceedings of the 14th International Symposium on Detonation. Coeur d'Alene Resort, Idaho, 2010: 861-870.. [9] Yoh J J, McClelland M A, Maienschein J L, et al. Simulating thermal explosion of cyclotrimethylenetrinitraminebased explosives: Model comparison with experiment[J]. Journal of Applied Physics, 2005, 97, 083504: 1-11. [10] Yoh J J, McClelland M A, Maienschein J L, et al. Simulating thermal explosion of octahydrotetranitrotetrazinebased explosives: Model comparison with experiment[J]. Journal of Applied Physics, 2006, 100, 073515: 1-9. [11] Dickson P M, Asay B W, Henson B F, et al. Measurement of phase change and thermal decomposition kinetics during cookoff of PBX 9501[C]//Shock Compression of Condensed Matter-1999.2000: 837-840. [12] Perry W L, Gunderson J A, Dickson P M. Application of a reversible four-step HMX kinetic model to an impactinduced friction ignition problem[C]//The 14th International Detonation Symposium. Coeur d'Alene Resort, Idaho, 2010. [13] Land T A, Siekhaus W J, Foltz M F, et al. Condensed-phase thermal decomposition of TATB investigated by AFM and STMBMS[C]//The 10th International Detonation Symposium. Boston, MA, 1993. 期刊类型引用(14)
1. Kebin Zhang,Wenbin Li,Changfang Zhao,Zhifang Wei,Shuxia Zhang,Jin Li. Study on cook-off characteristics and thermal safety venting area of RBOE charge. Defence Technology. 2025(01): 271-287 . 必应学术
2. 郭璐,智小琦,屈可朋,柳星河,贾杰,李劲. HMX基含铝压装炸药慢速烤燃点火时刻的压力参量计算. 爆炸与冲击. 2024(06): 50-61 . 本站查看
3. 胡平超,李涛,刘仓理,傅华. 初始空腔体积率对约束HMX基PBX-3炸药慢烤相变的影响. 爆炸与冲击. 2023(06): 54-61 . 本站查看
4. 邓海,李刚,任新联,梁争峰,李吉祯,张亚俊. 某含能黏结剂PBX炸药烤燃特性的数值模拟与实验研究. 兵器装备工程学报. 2022(08): 7-12 . 百度学术
5. 刘瑞鹏,贾宪振,郭洪卫,李鸿宾. FOX-7/HMX混合炸药烤燃试验的数值计算. 科学技术与工程. 2021(09): 3606-3611 . 百度学术
6. LIU JiPing,LIU LiLi,LIU XiaoBo. Development of high-energy-density materials. Science China(Technological Sciences). 2020(02): 195-213 . 必应学术
7. LIU JiPing,YANG WeiWei,LIU Ying,LIU XiaoBo. Mechanism and characteristics of thermal action of HMX explosive mixture containing high-efficiency fuel. Science China(Technological Sciences). 2019(04): 578-586 . 必应学术
8. 蒋超,闻泉,王雨时,张志彪,王光宇. 不敏感弹药烤燃试验技术综述. 探测与控制学报. 2019(02): 1-9 . 百度学术
9. 常天笑,闻泉,王雨时,王光宇,蒋超,闫丽. 引信内部装药对烤燃试验响应的影响. 探测与控制学报. 2019(05): 17-24 . 百度学术
10. 余一,张蕾,姜胜利,陈军. TATB在高温下分解的动力学模拟. 高压物理学报. 2018(01): 57-64 . 百度学术
11. 袁俊明,张林炎,唐鑫,刘玉存,张喜亮,陈银刚,李硕. 小型传爆装置慢燃实验及数值计算. 兵工学报. 2017(08): 1541-1546 . 百度学术
12. 李文凤,余永刚,叶锐,杨后文. 不同升温速率下AP/HTPB底排装置慢速烤燃的数值模拟. 爆炸与冲击. 2017(01): 46-52 . 本站查看
13. 姚奎光,钟敏,代晓淦,吕子剑. 缓慢热作用下PBX-9炸药的响应特性. 火炸药学报. 2015(06): 56-60 . 百度学术
14. 刘文杰,李小东,王晶禹,黄业令,宋原. 升温速率对烤燃弹温度影响的数值模拟. 中北大学学报(自然科学版). 2015(04): 440-445 . 百度学术
其他类型引用(23)
-