Effect of Near-Field Overpressure Enhancement of Reactive Material on Low Collateral Damage Ammunition
-
摘要: 为了实现低附带毁伤弹药的近场爆炸威力增强效应,提出在分装式低附带毁伤弹药的重金属颗粒嵌层中加入活性元成分,以期增强近场超压与比冲量。开展了不同含量活性元的静爆实验,利用自由场压力测试系统测得爆炸后近场及中远场的冲击波压力曲线。结果表明:在重金属颗粒嵌层中加入一定含量的活性元后,冲击波超压峰值和比冲量在37.5倍装药直径处分别提高31.6%和21.3%。根据实验结果,利用数值模拟确定了Miller反应速率模型参数,讨论了活性元后燃反应能量释放规律以及活性元组分反应度随时间的变化关系,在充分燃烧的理想情况下,活性元二次燃烧持续时间可达300 ms,且活性元含量区间极有可能存在最优配比。研究结果可为分装式低附带毁伤武器的近场冲击波区域增强效应及其工程化设计提供参考。Abstract: To realize the effect of the near-field shock wave enhancement of low collateral damage ammunition, this study proposes to add reactive material into the embedding layer of heavy metal particle in the sub-package low collateral damage ammunition to enhance the near-field overpressure and specific impulse. Static explosion experiments with different contents of reactive material were carried out, and the pressure curves of the near field and the mid-far field shock wave after the explosion were measured by the free field pressure test system. The results show that: after adding a certain content of reactive material into the heavy metal particle intercalation, the peak value of the overpressure and the specific impulse of shock wave at 37.5 times charge diameter are increased by 31.6% and 21.3%, respectively. According to the experimental results, the parameters of the Miller reaction rate model were determined by numerical simulation, and the law of energy release after-combustion reaction of the active element and the time-dependent change of the reactivity of the reactive material components were discussed. The duration of the secondary combustion can reach 300 ms under the ideal situation of sufficient combustion, and an optimal proportion of the range of reactive material content is likely exist. This research provides considerable development direction and application prospects for the regional enhancement effect of near-field shock wave and its engineering design of sub-packaged low collateral damage weapons.
-
夹芯结构由于具有较好的吸能缓冲性能以及较高的比强度和比刚度,被广泛应用于航空航天、交通运输、航海等领域。近年来随着对夹芯板研究的不断深入,提出了各种桁架夹芯和褶皱型夹芯等形式。其中,层级结构思想引起了众多学者的关注。Lakes[1]通过研究发现,层级结构与非层级结构相比具有更加显著的刚度和强度,同时也具有足够的稳定性及可靠性。Zhang等[2-3]通过梯形波纹板夹芯结构的静态压缩和低速冲击实验,研究了不同芯层材料、不同载荷下结构的能量吸收和变形过程,结果表明不同芯层材料和冲击能量下结构的能量吸收性能不同。Hou等[4]对不同层数的梯形夹芯板结构进行了实验研究,发现不同层数结构的能量吸收性能不同。Meza等[5]对多层级桁架结构的演变模式进行了分析,提出了几种不同结构的破坏形式和变形模式。Kooistra等[6]研究了层级波纹板结构的6种失效模式,并给出了结构在不同相对密度及角度下的破坏机理图。Velea等[7]提出了一种制备二级夹芯结构的方法,将一级波纹板结构芯层替换为多孔填充芯层,通过实验和理论分析发现,二级结构有着良好的刚度和强度。Wu等[8]将金字塔型夹芯结构的芯层替换为层级夹芯结构,理论推导出了结构的等效刚度,通过实验研究了芯层几何参数对结构变形模式的影响,并总结出结构的6种不同变形模式。Wu等[9]设计了二级金字塔形夹芯结构,进行了理论分析和数值模拟,结果表明二级芯层结构具有较好的强度和稳定性。Liu等[10]将传统格栅夹芯结构芯层板替换为蜂窝夹芯结构,通过实验研究和数值模拟研究发现,在压缩载荷下二级芯层结构有着较好的能量吸收性能和稳定的变形性能,在轻量化和能量吸收方面具有更好的工程应用前景。
多孔结构及蜂窝结构有着良好的力学性能和能量吸收性能[11]。Jing等[12]通过实验研究了3种不同多孔芯层夹芯板结构在脉冲载荷下夹芯板的变形模式,结果表明不同芯层夹芯结构的变形模式不同,芯层和面板尺寸不同的夹芯结构的力学性能不同。Chen等[13-14]根据层级蜂窝结构的设计理论设计了不同芯层形式的层级蜂窝结构,并由理论分析得到结构变形模式,利用有限元对结构变形模式进行了验证,与传统蜂窝结构相比,层级蜂窝结构具有更好的能量吸收性能。Yin[15]、Qiao[16]等利用有限元软件对层级蜂窝结构进行了数值模拟,研究结果表明层级蜂窝结构的力学性能远优于传统蜂窝。
已有的研究主要集中于传统波纹板夹芯结构和层级蜂窝结构,关于多层级夹芯结构的力学性能与能量吸收的研究较少。本研究基于波纹板夹芯结构,设计了3种二级波纹板夹芯结构,通过理论分析和数值模拟相结合的方法研究一级和二级夹芯结构的变形模式与能量吸收性能,理论预测结构的临界失效载荷;建立不同芯层层数的二级波纹板夹芯结构的有限元模型,研究不同芯层厚度对多层级波纹板夹芯结构变形模式和能量吸收性能的影响,并与传统的一级夹芯波纹板结构进行对比,分析其比吸能与结构效率。
1. 层级夹芯结构的芯层设计和理论分析
1.1 芯层结构设计
根据层级蜂窝结构的设计理论,将一级梯形波纹板夹芯结构的芯层替换为三角形多孔夹芯结构,共设计了3种芯层结构,如图1所示。图1(a)为一级波纹板夹芯结构的芯层,斜板长度为la,顶板长度为l,斜板角度为
ω ,芯层板厚度为ta;二级结构的芯层如图1(b)、图1(c)和图1(d)所示,其中图1(b)为二级单层结构,图1(c)和图1(d)分别为二级双层和三层结构,二级结构芯层由大支撑面板和小支撑构成,水平方向小支撑为横向小支撑。二级结构中不同层数结构的尺寸参数相同,斜向大支撑面板长度为la,厚度为ta,小支撑长度为lb,厚度为tb,小支撑与外板角度均为θ ,外板与水平方向夹角为ω 。结构垂直于纸面方向的厚度为30 mm。1.2 理论分析
小支撑厚度与大支撑面板厚度相同的结构模型如图2所示,在外部载荷Fx和Fy作用下考虑横向小支撑的变形,基于小变形假设,由受力分析可以得到结构变形和力之间的关系,通过几何分析可以得到合力与变形之间的关系
{Fx2=Fncosω+FtsinωFy2=Fnsinω−Ftcosω (1) {δx=δncosω+δtsinωδy=δnsinω−δtcosω (2) 式中:Fx和Fy分别为沿x轴和y轴的外部载荷,
δx 和δy 为沿x轴和y轴的变形,Fn和Ft分别为夹芯夹层在外部载荷下的轴向压缩分量和剪切分量,δn 和δt 分别为夹芯夹层的轴向变形和剪切变形。当只有Fx作用时,
δy =0,δx ≠0,此时有Fn=EAδnl=Fx2cosω (3) δn=FnlEA (4) 式中:EA为夹芯夹层截面的抗拉刚度。由结构几何关系可知
δx=δncosω (5) 当只有Fy作用时,
δy ≠0,δx =0,此时有Fn=EAδnl=Fy2sinω (6) 同样有
δy=δnsinω (7) 当结构芯层为二级结构时,如图2(d)所示,(1)式~(7)式仍然成立。根据文献[17]给出的不同尺寸比下结构的失效机制,当
lb/la>√εy/πsinθ (其中εy 为材料屈服应变)时层级结构出现大支撑面板塑性屈服的失效模式。本研究中结构模型的尺寸满足上述条件,因此失效模式为大支撑面板发生塑性屈服,即大支撑面板达到材料屈服应力时结构失效。由(1)式和(2)式可得单胞结构大支撑面板在y方向上发生塑性屈服时结构的合力Fy′=(2n+2)σybtasinω (8) 式中:n为二级结构层数,b为结构垂直于纸面方向的厚度,σy为屈服应力。
根据结构变形模式,考虑横向小支撑的变形,如图3所示,Fb为横向小支撑在y方向的力,
δb 为δn 在y方向的分量,由此可知δb=Fbl3b3D (9) 弯曲刚度
D=Et3b12(1−ν2) (10) 式中:E为弹性模量,
ν 为泊松比。由此可以得到Fb=2nσyt3bsinω4(1−ν2)l2b (11) 式中:tb和lb分别为二级结构芯层小支撑厚度和长度。二级结构单胞在y方向的合力为
Fy=Fy′+Fb (12) 2. 有限元模型
有限元模型主要包括3部分:刚性压板、刚性支撑底板和波纹板夹芯结构。波纹板夹芯结构由芯层和上、下面板构成,上、下面板设置为刚体,厚度为2 mm。图4为二级双层结构有限元模型。波纹板夹芯结构和下部底板采用壳单元,刚性压板采用实体单元。刚性压板速度V=1 m/s,下部刚性底板固定,中部波纹板夹芯结构芯层与上、下两面板之间采用绑定约束。结构材料为铝合金6061-T6,材料参数如表1所示,其中
ρ 为密度。图5为材料应力-应变曲线[6]。结构单元类型采用四节点线性缩减积分壳单元(S4R),单元特征长度L=1 mm。芯层和上、下面板为自接触,刚性压板与上面板、底板与下面板为面面接触,静摩擦系数和动摩擦系数分别为0.20和0.15[18]。表 1 6061-T6铝合金材料性能Table 1. Material properties of aluminum 6061-T6Material ρ/(g·cm–3) σy/MPa E/GPa ν Al 6061-T6 2.700 251 69 0.33 2.1 准静态分析
利用LS-DYNA进行准静态分析时[19]:首先在模拟过程中结构动能要相对于内能较小,其次结构的力-位移曲线相对于速度的变化是稳定的。如图6(a)所示,结构在压板速度为0.5 m/s和1.0 m/s时动能相对于内能较小,结构整体内能受速度变化的影响较小;图6(b)为结构在压板速度为0.5 m/s和1.0 m/s时的力-位移曲线,可以看出两条曲线的重合度较高,速度变化对结构力-位移曲线的影响较小。通过对结构进行网格敏感性验证,并考虑到计算时间等因素,计算时取压板的速度为1.0 m/s,网格单元尺寸为1 mm。
2.2 能量吸收指标
为了评估结构的能量吸收性能[20],选取总能量吸收(EA)、平均载荷(Pm)、比吸能(
ξEA )、载荷效率(AE)4个评价指标。总能量吸收EA为结构在变形过程中吸收的能量
EA=∫δ0F(x)dx (13) 式中:
δ 为结构的变形位移(取δ =70 mm),F(x)为压缩过程中的瞬时载荷。平均载荷Pm定义为单位变形长度的能量吸收
Pm=EA/δ (14) 比吸能
ξEA 定义为结构单位质量的能量吸收ξEA=EA/m (15) 式中:m为结构总质量。比吸能越大,结构的能量吸收性能越好。
载荷效率定义为平均载荷与峰值载荷的比值
AE=Pm/Pk (16) 式中:Pk为峰值载荷。较好的吸能结构应该具有较大的载荷效率,并且力-位移曲线波动较小。
3. 模拟结果与分析
3.1 芯层板厚度的影响
建立一级和二级结构模型,其中一级结构芯层壁厚为2 mm,二级结构模型中ta=tb,即二级结构内部小支撑厚度与外部大支撑面板厚度相同。层级结构芯层板厚度分别为0.565、0.75、1.00、1.25 mm,结构几何参数如表2所示。图7给出了一级和二级层级结构在不同芯层板厚度下的力-位移曲线。由图7可知,二级结构平台力明显高于一级结构。一级结构在靠近芯层上部产生塑性铰,随着变形的增大未产生新的塑性铰,力-位移曲线变化平稳。二级结构随着变形增大发生折叠塑性变形,力-位移曲线会随之发生波动;在芯层厚度为1.25 mm时二级单层结构由于变形模式发生变化,力-位移曲线发生波动而随之上升;随着二级结构层数的增加,结构平台力增大;随着二级结构厚度的增大,其平台力增大。
表 2 模型尺寸Table 2. Size parameters of the modelType l/mm la/mm lb/mm ta/mm tb/mm θ/(°) ω/(°) First order 50 100 2 60 Second order-1 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 Second order-2 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 Second order-3 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 图8所示为不同芯层厚度下结构失效载荷理论值与模拟值的对比。对于数值模拟结果,当结构的承载力达到峰值时认为结构失效,即取力-位移曲线上的峰值载荷为结构失效载荷,可以看出理论值与模拟值吻合较好。当结构厚度为0.565 mm时,由于其厚度较小,在压缩载荷下的变形模式受屈曲影响较大,故理论值较模拟值偏大;当结构厚度为1.25 mm时,小支撑端部变形受边界条件的影响较大,理论值相对于模拟值偏小。
3.2 变形模式
图9所示为芯层板厚度t=2 mm时一级结构在不同应变下的变形模式,应变为结构在受载荷方向上的变形位移与结构上、下两面板之间原始距离之比。由图9可知:一级结构首先发生屈曲,在靠近芯层斜支撑上部产生塑性铰,随着变形的增大,结构发生失稳,从而产生非对称变形。图10所示为二级结构在应变为0.3时不同芯层厚度下的变形模式,可以看出:当厚度较小时,结构的变形发生在整体斜支撑与整体水平支撑的交界处,整体斜支撑内部的小支撑发生逐层折叠的塑性变形;当芯层厚度为1.25 mm时,结构斜支撑整体产生塑性铰而发生非对称变形;厚度变化对二级双层和三层结构的变形模式影响较小。在芯层厚度较小时,二级结构主要为小支撑整体逐层发生塑性变形,从而有更多能量转换为非弹性能,提高了结构整体的能量吸收能力。
3.3 能量吸收性能分析
二级结构变形模式主要为小支撑逐层折叠的塑性变形,相对于一级结构,其在变形过程中小支撑不断发生塑性变形产生塑性铰,从而有更多的能量转换为非弹性能,提高了结构的能量吸收性能。图11为不同芯层厚度结构的比吸能柱状图。由图11可知:层级波纹板夹芯结构的比吸能高于传统波纹板夹芯结构;随着结构芯层厚度增大,结构的比吸能不断提高;由于二级单层结构变形模式的变化,芯层厚度较小时二级单层结构的比吸能高于二级双层及三层结构,厚度较大时二级单层结构的比吸能低于二级双层及三层结构;二级双层结构的比吸能略高于二级三层结构。
图12为结构载荷效率随芯层厚度的变化关系。随着芯层厚度的增大,结构的载荷效率提高,二级单层结构由于变形模式的变化,载荷效率发生波动。由于二级单层结构在芯层厚度为1.25 mm时结构变形模式为斜支撑整体发生失稳,相对于小支撑发生折叠变形时其能量吸收性能降低,从力-位移曲线可以看出,结构平台力减小导致结构的载荷效率下降。二级双层和三层结构的载荷效率不断增加,二级三层结构的载荷效率高于二级双层结构。
4. 结 论
设计了不同层数的层级波纹板夹芯结构,利用数值模拟方法研究了在压缩载荷下层级夹芯结构的变形规律与能量吸收性能,理论推导了结构临界失效载荷公式;分析了结构参数对其变形模式和能量吸收性能的影响,并与一级结构进行了对比分析,得到如下研究结果。
(1)理论分析得到的结构失效载荷与数值模拟结果吻合较好。
(2)在准静态压缩载荷作用下,一级结构首先发生屈曲,随着变形增大,结构发生失稳而产生非对称变形;二级单层结构在厚度为1.25 mm时由于芯层整体产生塑性铰而失稳;二级双层和二级三层结构在整体斜支撑与整体水平支撑的交界处,整体斜支撑内部的小支撑发生逐层折叠的塑性变形。
(3)二级波纹板夹芯结构的比吸能显著大于一级结构;二级结构芯层小支撑发生逐层折叠的塑性变形时,结构能量吸收性能较好;随着芯层厚度的增大,二级结构的比吸能和载荷效率增加,芯层厚度较小时二级单层结构的比吸能高于二级双层和三层结构,二级双层结构的比吸能略大于二级三层结构。
-
表 1 活性重金属颗粒嵌层配方
Table 1. Embedded formulation of reactive heavy metal particle
Serial number Mass fraction of heavy metal particle embedded/% WC diameter/µm WC RM Additive LCD-1 90 0 10 150–250 LCD-2 82 10 8 150–250 表 2 JWL状态方程参数
Table 2. Parameters of JWL equation of state
Explosives and RM A/GPa B/GPa R1 R2 ω Q/(kJ·g−1) E/(kJ·g−1) Composition B/0% Al 524.63 7.678 4.2 1.1 0.34 0 4.95 Composition B/10% Al 524.63 7.678 4.2 1.1 0.34 15.67 4.95 表 3 冲击波超压峰值数值模拟结果与实验对比
Table 3. Comparison of numerical simulation results of the peak overpressure of shock wave with experiments
R/m Δp/kPa δ/% Δpe Δps 1.5 169.2 163.6 3.3 2.0 86.0 89.5 4.4 3.0 40.0 42.8 3.3 -
[1] SIRAK M. Air Force focuses on new low-collateral-damage warhead for small bomb [N]. Defense Daily, 2006, 230(10): 1. [2] 刘意. 高密度惰性金属炸药爆轰与粒子流形成过程研究 [D]. 北京: 北京理工大学, 2015: 7–37.LIU Y. Detonation of dense inert metal explosive and the formation of particles flow [D]. Beijing: Beijing Institute of Technology, 2015: 7–37. [3] 申超. 重金属粉末嵌层CFRP壳体内爆下低附带毁伤特性表征 [D]. 北京: 北京理工大学, 2015: 9–22.SHEN C. The low collateral damage characterization of CFRP shell structure with heavy mental powder embedded as a layer under implosion [D]. Beijing: Beijing Institute of Technology, 2015: 9–22. [4] FROST D L, ORNTHANALAI C, ZAREI Z, et al. Particle momentum effects from the detonation of heterogeneous explosives [J]. Journal of Applied Physics, 2007, 101(11): 113529. doi: 10.1063/1.2743912 [5] FROST D L, GRÉGOIRE Y, PETEL O, et al. Particle jet formation during explosive dispersal of solid particles [J]. Physics of Fluids, 2012, 24(9): 091109. doi: 10.1063/1.4751876 [6] XUE K, YU Q Q, BAI C H. Dual fragmentation modes of the explosively dispersed granular materials [J]. The European Physical Journal E, 2014, 37(9): 88. doi: 10.1140/epje/i2014-14088-y [7] 黄德雨, 王坚茹, 陈智刚, 等. 炸药配比对陶瓷低附带毁伤战斗部能量输出的影响 [J]. 弹箭与制导学报, 2012, 32(1): 108–110, 130. doi: 10.3969/j.issn.1673-9728.2012.01.032HUANG D Y, WANG J R, CHEN Z G, et al. Influence of explosive radio on energy output of ceramic low collateral damage warhead [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2012, 32(1): 108–110, 130. doi: 10.3969/j.issn.1673-9728.2012.01.032 [8] 霍奕宇, 王坚茹. 某新型破片材料在低附带毁伤弹药中的应用研究 [J]. 科学技术与工程, 2015, 15(35): 175–178. doi: 10.3969/j.issn.1671-1815.2015.35.031HUO Y Y, WANG J R. Application of the research on the low collateral damage ammunition about new materials fragment [J]. Science Technology and Engineering, 2015, 15(35): 175–178. doi: 10.3969/j.issn.1671-1815.2015.35.031 [9] 霍奕宇, 王坚茹, 陈智刚, 等. 碳纤维壳体壁厚对陶瓷球初速及性能的影响 [J]. 爆破器材, 2016, 45(1): 30–33, 38. doi: 10.3969/j.issn.1001-8352.2016.01.007HUO Y Y, WANG J R, CHEN Z G, et al. Influence of thickness of carbon fiber shell on initial velocity and capability of ceramic ball [J]. Explosive Materials, 2016, 45(1): 30–33, 38. doi: 10.3969/j.issn.1001-8352.2016.01.007 [10] 李俊承, 樊壮卿, 梁斌, 等. 一种低附带弹药金属颗粒定向加载技术 [J]. 爆炸与冲击, 2018, 38(4): 869–875. doi: 10.11883/bzycj-2016-0376LI J C, FAN Z Q, LIANG B, et al. Experimental study on directed loading metal particles of low collateral damage ammunition [J]. Explosion and Shock Waves, 2018, 38(4): 869–875. doi: 10.11883/bzycj-2016-0376 [11] 刘俊, 姚文进, 郑宇, 等. 低附带毁伤弹药的炸药/钨粉质量比对钨粉抛撒特性的影响 [J]. 含能材料, 2015, 23(3): 258–264. doi: 10.11943/j.issn.1006-9941.2015.03.011LIU J, YAO W J, ZHENG Y, et al. Effect of explosive/tungsten powder mass ratio for LCD ammunition on dispersal characteristics of tungsten powder [J]. Chinese Journal of Energetic Materials, 2015, 23(3): 258–264. doi: 10.11943/j.issn.1006-9941.2015.03.011 [12] 左腾, 皮爱国. 亚毫米级金属微粒的爆炸驱动: 模型与数值计算 [J]. 兵工学报, 2016, 37(Suppl 2): 165–170.ZUO T, PI A G. Detonation driving of submillimeter-sized metal particles: model and numerical calculation [J]. Acta Armamentarii, 2016, 37(Suppl 2): 165–170. [13] 杨世全, 孙传杰, 钱立新, 等. 非金属壳体低附带战斗部实验与破片飞散分析 [J]. 高压物理学报, 2018, 32(4): 045103. doi: 10.11858/gywlxb.20170573YANG S Q, SUN C J, QIAN L X, et al. Experimentation and fragment flight analysis of low-collateral-damage warhead with nonmetal shell [J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 045103. doi: 10.11858/gywlxb.20170573 [14] 冯吉奎, 皮爱国, 刘源, 等. 爆炸驱动亚毫米级金属颗粒群的飞散特性 [J]. 高压物理学报, 2019, 33(6): 065104. doi: 10.11858/gywlxb.20190741FENG J K, PI A G, LIU Y, et al. Scattering characteristics of sub-millimeter metal particle group driven by explosion [J]. Chinese Journal of High Pressure Physics, 2019, 33(6): 065104. doi: 10.11858/gywlxb.20190741 [15] 仲倩, 王伯良, 邓金榜. 化学活性材料对温压炸药冲击波参数影响规律研究 [C]//中国力学大会2011暨钱学森诞辰100周年纪念大会论文集. 哈尔滨: 中国力学学会, 2011.ZHONG Q, WANG B L, DENG J B. Study on effect rule of chemically active material on blast wave of thermobaric explosive[C]//Proceedings of the Chinese Society of Theoretical and Applied Mechanics. Harbin: Chinese Society of Mechanics, 2011. [16] 王晓峰, 冯晓军. 温压炸药设计原则探讨 [J]. 含能材料, 2016, 24(5): 418–420. doi: 10.11943/j.issn.1006-9941.2016.05.00XWANG X F, FENG X J. Discussion on the design principle for thermobaric explosives [J]. Chinese Journal of Energetic Materials, 2016, 24(5): 418–420. doi: 10.11943/j.issn.1006-9941.2016.05.00X [17] JIANG C, LU G, MAO L, et al. Effects of aluminum content on the energy output characteristics of CL-20-based aluminized explosives in a closed vessel [J]. Shock Waves, 2021, 31(2): 141–151. doi: 10.1007/s00193-021-01001-1 [18] 隋树元, 王树山. 终点效应学 [M]. 北京: 国防工业出版社, 2000.SUI S Y, WANG S S. Terminal effects [M]. Beijing: National Defense Industry Press, 2000. [19] 田少康, 李席, 刘波, 等. 一种RDX基温压炸药的JWL-Miller状态方程研究 [J]. 含能材料, 2017, 25(3): 226–231. doi: 10.11943/j.issn.1006-9941.2017.03.009TIAN S K, LI X, LIU B, et al. Study on JWL-miller equation of state of RDX-based thermobaric explosive [J]. Chinese Journal of Energetic Materials, 2017, 25(3): 226–231. doi: 10.11943/j.issn.1006-9941.2017.03.009 [20] 张宝銔, 张庆明, 黄风雷. 爆轰物理学 [M]. 北京: 兵器工业出版社, 2001.ZHANG B P, ZHANG Q M, HUANG F L. Detonation physics [M]. Beijing: Weapon Industry Press, 2001. [21] MILLER P J, GUIRGUIS R H. Experimental study and model calculations of metal combustion in Al/Ap underwater explosives [J]. MRS Online Proceedings Library, 1992, 296(1): 299–304. doi: 10.1557/PROC-296-299 [22] 张奇, 白春华, 刘庆明, 等. 一次引爆燃料空气炸药及其爆炸效应研究 [J]. 实验力学, 2000, 15(4): 448–453. doi: 10.3969/j.issn.1001-4888.2000.04.014ZHANG Q, BAI C H, LIU Q M, et al. Investigation on single igniting fuel-air explosive and its explosion effect [J]. Journal of Experimental Mechanics, 2000, 15(4): 448–453. doi: 10.3969/j.issn.1001-4888.2000.04.014 [23] 黄菊, 王伯良, 仲倩, 等. 温压炸药能量输出结构的初步研究 [J]. 爆炸与冲击, 2012, 32(2): 164–168. doi: 10.11883/1001-1455(2012)02-0164-05HUANG J, WANG B L, ZHONG Q, et al. A preliminary investigation on energy output structure of a thermobaric explosive [J]. Explosion and Shock Waves , 2012, 32(2): 164–168. doi: 10.11883/1001-1455(2012)02-0164-05 -