Structure Optimization of Square Explosive Dispersion Device
-
摘要: 为了研究方形抛撒装置壳体破坏规律及燃料分散特性,开展了抛撒装置外场实验并利用LS-DYNA仿真软件进行了数值模拟,模拟结果与实验结果相吻合。进一步分析装置倒圆角及刻槽深度对壳体破裂效果及燃料分散速度的影响规律,结果表明:倒圆角与增加刻槽深度有利于减小壳体棱边处应力集中的影响,圆角半径增至10 mm或刻槽深度增至1.2 mm时,棱边处不再破裂;同一装置壳体边部和中部位置采用不同深度的刻槽可有效减少壳体破裂不均现象。当边缘处刻槽深度为1.2 mm、中间刻槽为1.6 mm时,壳体均匀破裂;当棱边处倒10 mm圆角,边部刻槽深度为0.8 mm,中部刻槽为1.2 mm时,抛撒装置既能满足壳体均匀开裂,又可提高壳体强度,同时可将燃料分散平均速度差值降低22%,从而有效提高燃料抛撒效率。Abstract: This paper presents a numerical model of the square dispersing device for simulating the process of shell failure and fuel dispersion by LS-DYNA software. Combined with the results of the field experiments, this model reveals in detail the influence of the fillet angle and groove depth on the shell rupture process and fuel dispersion speed. The results show that the shell edge would no longer rupture when the fillet radius increases to 10 mm or the groove depth increases to 1.2 mm, since different groove depth would effectively reduce the nonuniform shell rupture. And when the depth of edge and middle groove is 1.2 mm and 1.6 mm respectively, the shell is uniformly ruptured. In addition, a special dispersing device with 10 mm fillet angle, 0.8 mm edge groove depth and 1.2 mm middle groove depth, could not only make the shell uniformly ruptured, but also increase the strength of the shell. Meanwhile, it would reduce the average velocity difference of fuel dispersion by 22%, which effectively improve the fuel dispersing efficiency.
-
Key words:
- square dispersion device /
- fuel dispersion /
- LS-DYNA /
- crack
-
碳纤维复合材料(carbon-fibre reinforced composites,CFRP)具有轻质、高强、高韧、耐腐蚀、耐高温、可设计性强等优异特性,在轨道交通、新能源、火箭导弹以及航空航天(如飞机蒙皮、高空飞艇、发动机等)等领域具有重要应用价值[1–2]。然而,在某些应用领域中,碳纤维复合材料面临极端服役环境,从而对其性能提出了极其严苛的要求。例如,作为飞机蒙皮或高空飞艇中的关键材料,碳纤维复合材料时常遭遇低温环境和结冰问题,不仅引起飞行器表面平整度、外形、气动布局的改变[3–4]以及重量的增加,导致不可预知的飞行问题,而且由于纤维与基体的热膨胀系数差异较大,在低温环境中还可能出现收缩不协调导致复合材料内部产生裂纹,进而引发灾难性破坏[5]。因此,调控碳纤维复合材料在低温环境下的力学性能,并实现其在低温环境下的除冰功能是非常必要的,这对于保障碳纤维复合材料结构件,尤其是应用于高空飞行器结构件中的碳纤维复合材料服役安全性具有重大意义。
在低温除冰方面,当前最常用的机翼除冰方法有热气除冰、化学除冰、机械除冰和电热除冰等[6–8]。热气除冰主要利用发动机的热气加热进行除冰,但是可能会降低发动机的推力。化学除冰是通过喷涂防冰液等进行除冰,例如,Xue等[9]开发出聚二甲基硅氧烷(PDMS)/TiO2纳米复合材料,将其作为牺牲层喷涂在基材上,一旦在超疏水表面形成冰,便会导致牺牲层自动脱离,达到除冰效果,但是这种方法存在涂层寿命较短、需要重复喷涂的问题。电加热除冰是将加热垫布置在需要保护的部件表面,以代替热气除冰,在其服役期间可重复利用。然而,现在的电加热除冰技术尚不成熟,在民用飞机上的应用较少,尤其是机翼表面很难布置加热垫,因此需要进一步深入研究。
不锈钢超薄带,又名手撕钢(hand-torn steel,HTS),是一种强度高、柔韧性良好、导电率高的金属材料,在航天航空等高精尖领域具有重要应用[10]。手撕钢是一种各向同性材料,将手撕钢与碳纤维材料复合,有望弥补各向异性碳纤维层合板在某些方向上的性能短板,从而有效地增强碳纤维复合材料的综合力学性能。此外,手撕钢作为即时发热材料,将其与碳纤维材料复合不仅能够实现其整体温度调控,改善材料在低温服役环境下的力学性能,而且通过发挥其致热功能,可能实现飞机机翼等结构件表面的快速除冰。将碳纤维与手撕钢相结合,形成碳纤维-金属复合材料,发挥纤维和金属的优点,不仅可以实现飞行器的轻量化设计,还能够有效避免飞行器因表面结构件结冰而造成的安全事故。
目前,学者们已经对纤维-金属层合板进行了相关研究[11–13]。例如,李艳[11]采用阳极氧化法和热压成型法制备碳纤维金属铝层合板和碳纤维金属钛层合板,结果表明碳纤维金属钛层合板相比碳纤维金属铝层合板具有更好的力学性能,这是由添加金属种类所决定的。马其华等[12]制备了金属-CFRP混合管,发现CFRP与金属管相互协调可提升整体的承载能力,改善整个混合结构的能量吸收能力。研究表明,温度对纤维复合材料的力学性能有一定程度的影响[14–15]。然而,已有研究主要针对纤维-金属复合材料在常温环境下的力学性能,对于纤维-金属复合材料的电流驱动温度调节方案,及其在电流致热效应调控下的力学性能和除冰功能尚缺乏深入研究。
为此,本研究通过设计和裁剪不同形状的手撕钢,采用铺层固化方法制备“碳纤维预浸料”、“碳纤维预浸料-未裁剪手撕钢”和“碳纤维预浸料-弓字形手撕钢”3种复合材料,通过通电致热试验、三点弯曲性能测试和通电致热除冰测试,探索电流通路和电流强度引起的致热效应对复合材料力学性能和除冰功能的影响,以期为低温服役环境下碳纤维复合材料的力学性能优化和除冰功能实现提供指导,研究结果有望在航天航空等领域得到重要应用。
1. 样品制备
1.1 主要实验原料
本研究采用的CFRP为T300单向碳纤维增强环氧树脂薄层预浸料,型号为USN 20000(碳纤维与环氧树脂的质量比为6∶4),购自威海光威集团有限责任公司;手撕钢超薄带购自上海宏伟有限公司。材料参数如表1所示。
表 1 材料的基本属性Table 1. Basic material propertiesMaterial Tensile strength/MPa Density/(g·cm−3) Thickness/mm HTS 1130 7.26 0.02 CFRP 1800 1.80 0.20 1.2 三点弯曲样品制备
为了研究电流致热效应对材料力学性能的影响,将加热层与支撑层合并,即采用2层碳纤维预浸料和1层手撕钢交替叠层铺设,共11层,如图1(a)所示,其中预浸料采用单向铺层。加热层手撕钢起支撑和通电加热作用,另外两层手撕钢则主要起支撑作用。
根据三点弯曲测试标准《纤维金属层板弯曲性能试验方法》(GB/T 41049—2021),将样品尺寸设计为15 mm×80 mm,对加热层手撕钢进行裁剪处理,形成宽度约为3 mm的“弓字形”电流加热通路,如图1(b)所示。采用铺层固化方法制备样品:(1) 将碳纤维预浸料和手撕钢裁剪成合适的尺寸并解冻6 h;(2) 按照铺层方案逐层铺设样品;(3) 将样品进行模压和抽真空消泡处理,并在干燥箱中高温固化,干燥箱中的温度先以1~2 ℃/min的速率由室温升至120 ℃并保温90 min,然后以1~2 ℃/min的速率降至60 ℃以下进行脱模处理。所有样品均采用相同的制备工艺。最终制备出“碳纤维预浸料”、“碳纤维预浸料-未裁剪手撕钢”和“碳纤维预浸料-弓字形手撕钢”3组样品,如图1(c)所示。
1.3 通电致热除冰样品制备
按照图2(a)所示形状将手撕钢裁剪成“弓字形”,以研究电流通路宽度(a)对致热效应的影响。为获取加热层温度,按照图1(a)中第3~11层方案,制备尺寸为100 mm×100 mm的层合板,如图2(b)~图2(f)所示,其中a分别为5、10、15和20 mm。此外,为了验证除冰功能,按照图1(a)中第1~11层方案制备电流通路宽度为10 mm的样品。
2. 实验测试
对制备的样品进行测试,先对碳纤维-弓字形手撕钢长方形层合板通入不同的电流,测试其力学性能,然后对正方形层合板进行通电致热和除冰功能测试。
2.1 通电致热和三点弯曲测试
为研究电流强度对材料力学性能的影响,将样品接入电路中通电加热,如图3(a)所示,通过调压器(型号TDGC2-0.5 kVA)控制输出电压,调控样品中的电流大小(见表2),利用红外热成像仪(型号FLIR C2,测量误差±2 ℃,测量范围−20~180 ℃)采集样品表面温度数据,当样品表面温度达到稳态后,采用石棉对样品进行保温处理,然后断电,并立即在万能试验机(型号STD5000)上进行三点弯曲性能测试(如图3(b)所示),并与碳纤维预浸料-未裁剪手撕钢样品进行对比。测试中,跨距为40 mm,加载速率为2 mm/min,每组样品进行3次测试,最终测试结果为多次测量的平均值。
表 2 三点弯曲样品通电方案Table 2. Three-point bending sample energization schemeSamples Current/A CFRP No power HTS-CFRP No power Cut HTS-CFRP 0, 0.5, 1.0, 1.5 2.2 电流致热除冰升温特性测试
为测试碳纤维预浸料-弓字形手撕钢复合材料的电流致热除冰功能,将样品接入如图3(a)所示的电路中,调节调压器,使样品中的电流依次增大(间隔约为0.5 A),通电方案见表3,利用红外热成像仪测量样品表面温度,记录样品在不同电流下的稳态温度。此外,对电流通路宽度为10 mm的11层除冰验证样品的升温过程和除冰功能进行测试。
表 3 电流致热样品通电方案Table 3. Energizing samples with current heatingWidth of hand-torn steel/mm Current/A 20 0, 0.5, 1.0, 1.5, 2.0, 2.5 15 0, 0.5, 1.0, 1.5, 2.0, 2.5 10 0, 0.5, 1.0, 1.5, 2.0, 2.5 5 0, 0.5, 1.0, 1.5, 1.8 3. 结果与讨论
3.1 电流通路和电流强度对复合材料力学性能的影响
对3种复合材料在未通电条件下进行三点弯曲性能测试,结果如图4所示。本研究主要考虑比强度和比吸能两方面因素,其中,比强度是材料的强度与其表观密度之比,而比吸能ESA的计算公式为
ESA=∫x00Fmdx (1) 式中:F为三点弯曲载荷;x为压头位移;x0为破坏位移(载荷下降30%对应的压头位移);m为样品的表观密度,纯碳纤维预浸料、碳纤维预浸料-未裁剪手撕钢、碳纤维预浸料-弓形手撕钢复合材料的表观密度分别为1.80、2.03、2.02 g/cm3。图4(a)表明,3种复合材料的应力-应变曲线变化趋势相同,达到弯曲强度后,载荷迅速下降,复合材料发生破坏。由图4(b)和图4(c)可知:与纯碳纤维预浸料复合材料相比,碳纤维预浸料-未裁剪手撕钢复合材料的比强度有所提升(由0.817 N·m/kg提升至0.876 N·m/kg),但是比吸能有所降低(由412.9 J/kg降低为359.2 J/kg),碳纤维预浸料-弓字形手撕钢复合材料的比强度虽然有所降低(由0.817 N·m/kg降为0.733 N·m/kg),但是破坏应变显著增大,延缓了复合材料的破坏,从而使其比吸能略有提升(由412.9 J/kg提升为414.6 J/kg)。对比碳纤维预浸料-未裁剪手撕钢复合材料与碳纤维预浸料-弓字形手撕钢复合材料,可以看出,后者的比强度降低,但比吸能提高。这说明加入未裁剪的手撕钢可显著提升碳纤维预浸料复合材料的强度,加入裁剪的手撕钢可在一定程度上提升碳纤维预浸料复合材料的韧性。下面将对该实验结果进行进一步的机理分析。
在纯碳纤维预浸料复合材料的基础上添加未裁剪手撕钢后,比强度得到提升,而添加裁剪手撕钢后比强度有所降低,这可能是由于:(1) 手撕钢具有良好的强度和柔性,有助于提升复合材料的整体强度;(2) 与碳纤维预浸料不同的是(沿纤维方向的承载能力显著高于垂直纤维方向的承载能力),手撕钢是各向同性材料,能够很好地承载并传递层间应力,从而使其整体强度得到提升;(3) 裁剪后的手撕钢由于相邻通路间含有间隙,间隙处无法承受和传递载荷,且材料质量增加,致使其比强度有所降低。
从吸能方面来看,碳纤维预浸料-弓字形手撕钢层合板的比吸能显著提升,而碳纤维预浸料-未裁剪手撕钢层合板的比吸能有所降低,其原因可能是:(1) 对于碳纤维预浸料-未裁剪手撕钢层合板,虽然强度得到提高,但是极限应变减低(手撕钢的极限应变低于碳纤维预浸料的极限应变),且质量增加导致其比吸能降低;(2) 对于碳纤维预浸料-弓字形手撕钢层合板,裁剪手撕钢留下间隙,此处的碳纤维预浸料将在固化过程中结合到一起,从而在材料破坏过程中形成更多的断裂面,提高材料的能量吸收特性。
上述分析表明,碳纤维预浸料-弓字形手撕钢层合板与纯碳纤维预浸料层合板相比,比强度降低,比吸能显著提升,即可以在一定程度上延缓层合板的破坏,同时吸收更多能量以发挥其承载和防护作用。下面对碳纤维预浸料-手撕钢复合材料在不同电流作用下的力学性能做进一步研究。
分别对碳纤维预浸料-弓字形手撕钢层合板在不同电流强度下的力学性能进行测试,并与碳纤维预浸料-未裁剪手撕钢层合板(不通电)的力学性能进行对比分析,结果如图5所示。可以看出,通电后碳纤维预浸料-弓字形手撕钢层合板的比强度和比吸能提高,但是随着电流强度增加,比强度和比吸能降低。通电不仅能够提升碳纤维预浸料-手撕钢复合材料的强度,而且增大破坏时的极限应变,使其在破坏过程中吸收更多的能量,例如: 通入0.5 A的电流后,比强度大幅提升(由0.733 N·m/kg提高至0.957 N·m/kg),甚至超过了碳纤维预浸料-未裁剪手撕钢层合板,且变形伸长量较大,通电后碳纤维预浸料-弓字形手撕钢层合板的比吸能显著提升(由414.6 J/kg降低为579.4 J/kg),相对于纯碳纤维预浸料层合板,比强度和比吸能分别提升17.2% 和40.3%。这可能是由于通电后样品温度上升,达到一定温度时预浸料中的环氧树脂基体发生链段重排并进一步固化,导致环氧树脂分子链的交联密度增加[14],增强了碳纤维与树脂之间的结合强度,从而提高了结构整体的强度和能量吸收能力。但是,随着电流强度的增加,样品温度过高会引起部分分子链断裂[15],从而降低层合板的强度和能量吸收特性。因此,需要选择合适的电流,以满足实际应用中对层合板强度和吸能特性的要求。
综上所述,在纯碳纤维预浸料中加入手撕钢会提升材料的力学性能,添加经过裁剪处理的手撕钢后,与由纯碳纤维预浸料制备的样品相比,虽然比强度有所下降,但是比吸能显著提升。通电对于碳纤维-手撕钢复合材料的力学性能也有重要影响,通过调控电流强度可实现强度和能量吸收特性的优化。
3.2 碳纤维-手撕钢层合板的电流致热除冰功能
为了实现碳纤维-手撕钢层合板电流致热除冰功能,测试得到具有不同电流通路宽度(“弓”字形手撕钢宽度)的碳纤维-手撕钢层合板在通电后的稳态温度,结果如图6所示。由图6(a)可知,随着电流通路宽度的减小,表面温度的上升趋势存在一定的差异。其中:碳纤维预浸料-未裁剪手撕钢样品通入电流后,温度几乎不变,可能是由于整片手撕钢的电阻很小,导致发热功率很小,样品表面升温过程和散热过程几乎处于动态平衡状态,因此表面温度基本维持在某一稳定值;碳纤维预浸料-弓字形手撕钢样品的电阻显著增大,通电后样品表面温度随着电流增加而升高。
此外,还发现当电流一定时,随着电流通路宽度的减小,样品表面稳态温度逐渐升高,样品表面温度场逐渐均匀。当电流达到2.5 A时,电流通路宽度为20、15、10 mm的样品的表面温度分别达到65.7、78.5、101.0 ℃,如图6(b)~图6(d)所示。需要说明的是,当电流通路宽度为5 mm时,未能完成电流为2.5 A的实验,这是因为随着电流的增大,温度迅速上升,当电流为1.6 A时,样品表面温度达到152.0 ℃,如图6(e)所示,当电流达到1.8 A时,温度超过180 ℃,且一段时间后样品发生变形损坏。图6(b)~图6(e)的结果表明:当手撕钢宽度为15和20 mm时,温度场分布均较为分散,且温度提升缓慢,不利于大面积除冰;随着电流通路宽度减小,样品表面温度场趋于均匀;当手撕钢宽度为5和10 mm时,温度场分布较为均匀,可以实现大面积除冰。结合3.1节的力学性能分析,可以发现:碳纤维预浸料-弓字形手撕钢层合板中的裂纹或缺陷增加,导致力学性能下降;电流通路为10 mm的碳纤维-手撕钢层合板相比于其他层合板能够在保证良好力学性能的情况下,达到大面积除冰效果。
在2.5 A的电流下,电流通路宽度为10 mm的样品的加热功率为2 kW/m2,其温度变化过程如图7(a)~图7(f)所示。可见,接通电路后,样品表面温度迅速提升到49 ℃,如图7(a)所示,但初始热量分布不均匀;随着时间的延长,样品表面温度变化逐渐减缓;40 s时最高温度达到86.6 ℃,如图7(d)所示,样品表面高温区域趋于均匀;120 s时,最高温度基本稳定在94.9 ℃,样品中间区域的温度相对均匀,更有利于大范围除冰。此外,将边长约为15 mm的冰块置于样品表面,大约5 min后冰块完全融化,表明手撕钢铺层材料具有很好的除冰功能。
上述针对碳纤维-手撕钢复合材料的电流致热除冰功能研究表明,不同电流通路宽度对复合材料电流致热功能具有重要影响,电流通路宽度减小,复合材料在相同电流下的致热效应提升,可用来通电致热除冰。此外,电流强度对于碳纤维-手撕钢复合材料的致热效应也有重要影响,增加电流强度会引起复合材料温度过高,可能引起材料破坏失效。通过调控手撕钢尺寸,在碳纤维-手撕钢复合材料的力学性能与除冰功能之间找到平衡点。实验表明,手撕钢宽度为10 mm的碳纤维-手撕钢复合材料可以在不降低材料性能的基础上实现通电除冰的目的。
4. 结论与展望
设计和制备了碳纤维预浸料、碳纤维预浸料-手撕钢、碳纤维预浸料-弓字形手撕钢3种复合材料,通过通电实验和三点弯曲测试,获得了手撕钢电流通路和电流强度对复合材料力学性能的影响,据此设计了具有不同电流通路宽度的碳纤维-手撕钢复合材料,研究了其电流致热除冰功能,得到如下结论。
(1) 在纯碳纤维预浸料层合板中加入未裁剪手撕钢会提升比强度,加入弓字形手撕钢会在比强度有所降低的情况下提升能量吸收特性。电流强度对碳纤维-手撕钢复合材料的力学性能也有重要影响,施加电流时,其比强度和比吸能显著提升。因此,在实际应用过程中,可通过调节电流强度优化复合材料的强度和能量吸收特性,以满足不同领域对复合材料力学性能的要求,从而为该类复合材料的应用提供积极指导。
(2) 碳纤维-手撕钢复合材料的电流致热除冰特性与其中手撕钢的电流通路宽度有密切联系。随着电流通路宽度的减小,相同电流强度下复合材料的致热温度及其均匀性提高,但是随着电流强度提升,致热温度过高会引起层合板破坏失效。通过调控手撕钢的尺寸,可以平衡其力学性能和除冰功能。相对于其他尺寸,电流通路宽度为10 mm的碳纤维-手撕钢层合板在保证材料性能和安全性的前提下,实现通电致热和快速除冰。这对于碳纤维-金属复合材料的结构设计和功能应用提供了参考依据,并有望在航天航空等领域得到应用。
本研究对于电流强度和电流通路宽度变化下复合材料的测试次数较少,实验环境为室温,且仅对比了复合材料的强度和吸能,研究结果具有一定的局限性。后续研究将细化电流通路宽度和电流强度,同时对刚度等其他重要力学性能进行分析,结合实际应用场景,增加风洞试验,研究复合材料在不同条件下的力学性能和致热效果,实现金属-碳纤维层合板力学性能和功能特性优化。
-
表 1 TNT炸药的爆轰性能及JWL状态方程参数
Table 1. Detonation properties and JWL equation-of-state parameters of TNT explosive
ρ0/(g·cm–3) D/(m·s–1) p/GPa A/GPa B/GPa R1 R2 ω E/(GJ·m–3) 1.63 6 930 21 374 7.33 4.2 0.9 0.35 7 表 2 燃料计算参数
Table 2. Parameters for the computation of fuel
ρ0/(g·cm–3) C S1 S2 S3 γ0 1.0 1.65 1.92 –0.096 0 0.35 表 3 空气计算参数
Table 3. Parameters for the computation of air
ρ0/(g·cm–3) C0 C1 C2 C3 C4 C5 C6 E/(J·m–3) V0 0.001 25 0 0 0 0 0.4 0.4 0 2.5 × 105 1.0 表 4 装置结构计算参数
Table 4. Computation parameters of device structure
Structure Material ρ0/(g·cm–3) E ν F Ce Pe Shell 5A06 Al 2.75 0.68 0.35 0.45 40 5 Plate, central tube and tube cap 2A12 Al 2.78 0.70 0.35 0.35 40 5 表 5 倒圆角对壳体棱边及刻槽的影响
Table 5. Effect of chamfer on edge and groove of shell
Shell Position σmax/MPa Crack or not tcrack/μs No chamfer Edge 376 Yes 175 Groove 1 361 No Groove 2 355 No 0.8 cm chamfer Edge 374 Yes 205 Groove 1 369 No Groove 2 355 No 1.0 cm chamfer Edge 360 No Groove 1 381 Yes 185 Groove 2 356 No 1.2 cm chamfer Edge 359 No Groove 1 385 Yes 175 Groove 2 357 No 表 6 刻槽深度对壳体棱边及刻槽的影响
Table 6. Effect of groove depth on edge and groove of shell
Depth of
groove/mmPosition σmax/MPa Crack or not tcrack/μs Edge 376 Yes 175 0.6 Groove 1 361 No Groove 2 355 No Edge 371 Yes 175 0.8 Groove 1 368 No Groove 2 358 No Edge 378 Yes 175 1.0 Groove 1 376 Yes 185 Groove 2 366 No Edge 375 No 1.2 Groove 1 378 Yes 165 Groove 2 366 No Edge 369 No 1.6 Groove 1 381 Yes 145 Groove 2 364 No 表 7 多深度刻槽条件下增加刻槽2的深度对壳体破裂产生的影响
Table 7. Effect of increasing groove 2 depth on shell rupture under multi-depth grooving conditions
Depth of
groove 2/mmPosition σmax/MPa Crack or not tcrack/μs 1.4 Groove 1 378 Yes 165 Groove 2 370 No 1.5 Groove 1 382 Yes 165 Groove 2 373 Yes 175 1.6 Groove 1 380 Yes 165 Groove 2 377 Yes 165 表 8 两因素条件下增加刻槽2的深度对壳体破裂效果的影响
Table 8. Effect of increasing groove 2 depth on shell fracture under two factors condition
Depth of
groove 2/mmPosition σmax/MPa Crack or not tcrack/μs 1.0 Groove 1 383 Yes 165 Groove 2 364 No 1.1 Groove 1 382 Yes 165 Groove 2 371 Yes 185 1.2 Groove 1 382 Yes 165 Groove 2 376 Yes 165 表 9 方案参数
Table 9. Scheme parameters
Plan Chamfer
size/cmDepth of
groove 1/mmDepth of
groove 2/mm1 0 0.6 0.6 2 0 1.2 1.6 3 1.0 0.8 1.2 -
[1] 薛田, 徐更光, 黄求安, 等. 爆炸抛撒过程的研究进展 [J]. 科学技术与工程, 2015, 15(21): 60–67. doi: 10.3969/j.issn.1671-1815.2015.21.012XUE T, XU G G, HUANG Q A, et al. Review on explosive dispersion [J]. Science Technology and Engineering, 2015, 15(21): 60–67. doi: 10.3969/j.issn.1671-1815.2015.21.012 [2] LIU G, HOU F, CAO B. Experimental study of fuel-air explosive [J]. Combustion, Explosion and Shock Waves, 2008, 44(2): 213–217. doi: 10.1007/s10573-008-0028-7 [3] 韩宝成, 雷红飞, 徐豫新, 等. 液体灭火弹动态抛撒区域工程计算方法 [J]. 含能材料, 2011, 19(3): 310–314. doi: 10.3969/j.issn.1006-9941.2011.03.015HAN B C, LEI H F, XU Y X, et al. Computational model of dynamic dispersed area of liquid fire-extinguishing ammunition [J]. Chinese Journal of Enegretic Materials, 2011, 19(3): 310–314. doi: 10.3969/j.issn.1006-9941.2011.03.015 [4] 胡涛, 刁伟, 崔正辉, 等. 水袋在爆破水雾除尘技术中的应用 [J]. 水科学与工程技术, 2012(4): 52–54. doi: 10.3969/j.issn.1672-9900.2012.04.018HU T, DIAO W, CUI Z H, et al. Computational model of dynamic dispersed area of liquid fire-extinguishing ammunition [J]. Water Sciences and Engineering Technology, 2012(4): 52–54. doi: 10.3969/j.issn.1672-9900.2012.04.018 [5] 高重阳, 施惠基, 姚振汉. 薄壁柱壳在内部爆炸载荷下膨胀断裂的研究 [J]. 爆炸与冲击, 2000, 20(2): 160–167. doi: 10.3321/j.issn:1001-1455.2000.02.012GAO C Y, SHI H J, YAO Z H. Dynamic fracture of thin cylindrical shells subject to internal explosive loading [J]. Explosion and Shock Waves, 2000, 20(2): 160–167. doi: 10.3321/j.issn:1001-1455.2000.02.012 [6] 王晔, 白春华, 李建平. 弹壳体结构对燃料装药抛撒速率影响的数值模拟研究 [J]. 兵工学报, 2017, 38(1): 43–49. doi: 10.3969/j.issn.1000-1093.2017.01.006WANG Y, BAI C H, LI J P. Influence of shell structure on dispersing velocity of fuel-air mixture [J]. Acta Armamentarii, 2017, 38(1): 43–49. doi: 10.3969/j.issn.1000-1093.2017.01.006 [7] 张奇, 白春华, 刘庆明. 壳体对燃料近区抛散速度的影响 [J]. 应用力学学报, 2000, 17(3): 102–106. doi: 10.3969/j.issn.1000-4939.2000.03.019ZHANG Q, BAI C H, LIU Q M. Influence of shell casting on fuel near-field dispersal velocity [J]. Chinese Journal of Applied Mechanics, 2000, 17(3): 102–106. doi: 10.3969/j.issn.1000-4939.2000.03.019 [8] ZHANG Q, BAI C H, LIU Q M. Study on near field dispersal of fuel air explosive [J]. Journal of Beijing Institute of Technology, 1999, 8(2): 113–118. [9] 薛社生, 刘家骢, 彭金华. 液体燃料爆炸抛撒的近场阶段研究 [J]. 南京理工大学学报(自然科学版), 1997(4): 49–52.XUE S S, LIU J C, PENG J H. Study on the near-field phase of explosive throwing of liquid fuel [J]. Journal of Nanjing University of Science and Technology (Natural Science Edition), 1997(4): 49–52. [10] 王在成, 姜春兰. 机载布撒器弹仓总体技术的发展与分析 [J]. 飞航导弹, 1999(10): 21–28.WANG Z C, JIANG C L. Development and analysis of the overall technology of airborne dispenser silo [J]. Flying Missile, 1999(10): 21–28. [11] LEE E L, HORMIG H C, KURY J W. Adiabatic expansion of high explosive detonation products: UCRL-50422 [R]. Livermore: Lawrence Livermore National Laboratory, 1968. [12] 陈明生, 白春华, 李建平. 燃料抛撒的初始速率与结构动态响应数值模拟 [J]. 含能材料, 2015, 23(4): 323–329. doi: 10.11943/j.issn.1006-9941.2015.04.003CHEN M S, BAI C H, LI J P. Simulation on initial velocity and structure dynamic response for fuel dispersion [J]. Chinese Journal of Energetic Materials, 2015, 23(4): 323–329. doi: 10.11943/j.issn.1006-9941.2015.04.003 [13] MEYERS M A. Dynamic behavior of materials [M]. John Wiley & Sons, 1994: 124. -