Dynamic Response Characteristics of Aluminum Foam Sandwich Structure under Explosion Load in Cabin
-
摘要: 为探索爆炸载荷下舱内夹芯复合结构的动态响应特性与防护效能,采用小尺度舱室结构模型实验,结合有限元数值分析,开展了不同爆炸距离下舱内双层泡沫铝夹芯结构的动响应特性和变形模式研究。分析了不同爆距下舱内爆炸载荷的作用过程和时空分布特性,讨论了在初始冲击波、初始冲击波叠加各壁面二次反射波和舱内爆炸准静态压力3种载荷下泡沫铝夹芯结构的变形模式。爆炸载荷下舱室壁板承受的载荷依次为初始冲击波、各壁面二次反射波和准静态气压。炸药在靠近舱室一端处起爆时,初始冲击波在近端壁的局部效应明显,在远端壁的作用范围更大,与舱室中心爆炸相比,其爆轰产物波动次数更少。泡沫铝夹芯结构的变形过程可分为泡沫芯层压缩、局部凸起变形和整体挠曲变形3个阶段,对应迎爆面板局部凸起叠加整体挠曲大变形、局部凸起叠加整体挠曲大变形和整体挠曲大变形3种变形模式。Abstract: To study the dynamic response characteristics and protective effectiveness of sandwich composite structure subjected to internal explosion load in a cabin, small scale structure model experiments and finite element numerical simulations were performed, and the propagation and distribution characteristics of the explosion load in the cabin with different blasting distances were analyzed. The dynamic response and deformation mode of the aluminum foam sandwich structure under the initial shock wave, the superposition of reflected shock waves and the quasi-static pressure of the internal explosion were discussed. The experimental and numerical results showed that when the explosive was detonated near one end of the cabin, the localized deformation of the specimen caused by the initial shock wave is significant, while the deformation area of the specimen on the far wall is larger and smoother. Compared with the load conditions of explosions in the center of the cabin, the fluctuation of the shock wave are gentler. In addition, the deformation process of aluminum foam sandwich structure can be divided into three stages: foam core compression, local bulge deformation and overall deflection. The corresponding deformation modes include the facing blast panel with local rising on integral deformation, double panels with local rising on integral and large deformation of integral flexure.
-
Key words:
- mechanics of explosion /
- explosion in cabin /
- sandwich structure /
- aluminum foam /
- dynamic response /
- deformation mode
-
冲击波加载多孔固体介质过程中形成的极端高温高压条件可激活常规温压条件难以发生的化学反应,为材料合成开辟了一条新途径,尤其是对合成条件苛刻的材料具有显著优势。然而,冲击波作用导致粉末从初始状态至孔隙湮灭的过程中涉及多种机械、物理和化学变化,因而冲击波压缩下的反应过程十分复杂[1]。冲击诱导化学反应取决于外部加载条件和粉末混合物的初始状态等多种因素[2-3],反应物粉末的初始孔隙率是控制冲击反应热力学和动力学的一个非常重要因素[4]。Cooper 等[5]应用理论模拟证实,冲击波和孔隙表面相互作用导致孔隙崩塌产生的大量能量集聚有利于激活化学反应,因此探讨粉末初始孔隙率对诱发超快冲击化学反应的影响极为重要。
铌硅二元金属间化合物因具有低密度和优良的高温强度等物理特性,可作为潜在的高温结构材料和金属含能材料[6-7]。根据铌-硅二元相图,铌-硅系存在Nb5Si3、Nb3Si和 NbSi2共3种金属间化合物[8]。其中Nb5Si3的熔点(2 515 ℃)最高,密度 (7.16 g/cm3)相对较低,因其具有优良的高温性能受到关注[9]。这类材料的熔点高、晶体结构复杂,采用常规方法难以合成,但该类材料的合成反应能够放热,理论上讲,一旦被激活便可自维持发生反应。由于该激活势垒高,传统方法难以达到激活条件,而冲击波加载产生的极端高温高压对于激活该类反应具有一定优势。目前,已经开展了铌硅化合物的冲击合成工作。Vecchio等[10]应用冲击波加载合成了NbSi2,发现其冲击诱导反应机制与传统固态反应机制不同;Meyers等[11]也通过铌硅粉末混合物的冲击回收,发现塑性变形可以影响铌硅粉末间的反应。此外,还有一些理论模拟研究对铌硅金属间化合物的冲击合成进行探讨[12]。但是,所有铌硅粉末混合物的冲击反应研究,无论是实验结果还是理论分析都仅局限于较低冲击强度条件(飞片速度在2.0 km/s以下),冲击产物也仅获得了较低熔点的NbSi2,而且这些结果对于冲击化学反应的程度和机理解释存在分歧。而且,高熔点金属间化合物(如Nb5Si3)这类传统方法较难合成的材料在高冲击强度范围(飞片速度在2.0 km/s以上)的冲击合成及反应机制未有报道。
课题组前期开展了铌硅粉末混合物的冲击回收实验,通过高冲击强度下的冲击回收获得了高熔点Nb5Si3化合物[13],发现相同组分配比的铌硅粉末混合物在不同的冲击强度下合成了不同组成的铌硅化合物,且反应特征不同。为更清楚地理解高冲击强度下铌硅粉末混合物的冲击合成及反应机制,在前期工作基础上,本研究进一步对不同孔隙率的铌硅粉末混合物样品进行高冲击强度下回收,通过对冲击回收产物进行表征分析,探讨不同初始孔隙率对铌硅混合物冲击反应的影响。
1. 实验材料与方法
1.1 冲击样品制备
实验使用高纯Nb(约325目,Alfa Aesar公司)和Si(约325目,Alfa Aesar公司)粉末原料。将铌粉和硅粉按原子比为5∶3称量,V型混料器中氩气气氛下将粉末混合2 h,通过不同压力将混合好的粉末冷压成直径约为12 mm、厚度约为3 mm的圆片状试样。冲击样品的准备过程(原料称量、混合前的真空封装、压片和铜回收盒的真空封装)均在水氧监测氩气气氛保护的手套箱中完成。将装有样品的铜制回收盒密封。样品孔隙率的计算方法参考文献[14]。
1.2 冲击回收实验
应用二级轻气炮加载装置对样片进行冲击回收,装置示意图和实物如图1所示。将铜飞片(
24 mm × 3.0 mm)粘在发射弹丸上,由磁测速系统测量飞片速度,在飞片速度约2.35 km/s的最高冲击速度下进行冲击回收实验。利用质量平均法计算混合物样品的雨贡纽参数[15-16],采用阻抗匹配法和疏松粉末状态方程计算样品的冲击压力和温度[15,17],表1列出了不同孔隙率粉末混合物的冲击回收实验参数和计算结果。表 1 不同孔隙度的粉末混合物的冲击回收实验参数Table 1. Shock loading conditions of Nb-Si powder mixtures with different porositySample Flyer velocity/
(km·s−1)Density/
(g·cm−3)Porosity/
%Shock
pressure/GPaSecond shock
pressure/GPaShock
temperature/KNb-Si-P1 5.505 10 45 60 1 173 Nb-Si-P2 2.35 ± 0.02 4.895 20 39 60 1 625 Nb-Si-P3 3.980 35 29 60 2 256 1.3 冲击回收产物表征
图2显示了完整回收的铜回收盒在冲击前、后的外形变化。用机械切割方式从铜制回收盒中将回收的冲击产物分离取出,切掉产物样品外围边缘的含铜部分,留下中心主体部分用于表征分析。将产物样块逐级打磨,并清洗、干燥,储存在真空干燥箱中待表征。
使用X 射线衍射仪(丹东浩元仪器有限公司,DX-2700型)对样品进行物相表征,采用CuKa 靶(λ = 0.154 06 nm) 射线,衍射角2θ为20°~80°,步长为0.03°,每步计数时间1 s。使用日本电子公司JSM-7001F型场发射扫描电子显微镜 (SEM)进行微观形貌观察,并用SEM配备的能谱仪分析各物相的成分。使用差示扫描量热计(DSC)测量冲击回收产物的热响应曲线(Ar气氛下,加热速率为10 ℃/min)。
2. 结果和讨论
2.1 冲击回收产物的XRD结果分析
图3为不同孔隙率铌硅粉末冲击产物的XRD图谱。从图3(a)中10%孔隙率冲击样品的XRD结果可见,产物主要为Nb和Si的特征衍射峰及较弱的NbSi2衍射峰,表明仍有大部分Nb和Si未发生反应,仅有极少部分NbSi2相生成。由图3(b)中20%孔隙率铌硅粉末冲击产物的XRD结果可见,该孔隙率产物中虽仍有部分未反应的Nb和Si,但生成NbSi2相的程度有所提高,初始孔隙率提高后,冲击化学反应更易进行。此外,图3(c)中35%孔隙率铌硅粉末混合物的冲击产物XRD图谱显示,其衍射峰仅由Nb5Si3(包括α-Nb5Si3相和β-Nb5Si3相)的特征峰组成。与低孔隙率样品相比,在相同冲击强度下高孔隙率粉末样品更容易获得高熔点Nb5Si3金属间化合物。因此,初始孔隙率对粉末混合物在冲击压缩下物相生成、固相化学反应的影响非常明显。参考表1估算的冲击平均温度也可看出,随着样品初始孔隙率增加,较高的冲击温度(2 256 K)更有利于铌硅金属间化合物发生反应。
2.2 冲击回收产物的SEM结果分析
图4为10%初始孔隙率铌硅粉末冲击产物的SEM结果,其中图4(b)为图4(a)的局部放大图,图4(c)为样品EDS能谱分析结果,w为质量分数。由图4(a)可以看出,样品在该冲击压力(约为45 GPa)下铌硅颗粒被明显压实,没有明显物相生成或发生化学反应。图4(b)中局部区域放大图显示,与初态未冲击样品相比,冲击后的粉末颗粒呈细粒化、等轴状,说明冲击加载作用导致颗粒发生崩裂和相互嵌合。EDS能谱显示图4(b)中亮灰色细粒子的为Nb,周围相对较暗细粒为Si,可能由于NbSi2量极少,在图中未能明确分辨。在较小孔隙率即相对密实的混合反应物样品中,冲击波更多作用于粉末颗粒的机械和物理变形。由于孔隙率低,孔隙崩塌所致的能量贮集弱,极少能够激活并促发冲击化学反应。
图5为20%初始孔隙率铌硅粉末冲击回收产物的SEM结果,其中图5(b)为图5(a)的局部放大图,图5(c)为样品EDS能谱分析结果。图5(a)中,经冲击加载后20%孔隙率回收产物与10%孔隙率回收产物明显不同,随着冲击作用下温度不断升高,由于Si的熔点较低,而Nb的屈服强度和熔点较高,Si会发生熔融,并分布于变形的Nb周围,与基体特征类似,并且其中弥散着新相。如图5(b)所示,对富Si区域高倍放大观察可见,Si基体中分布着球形、结节状粒子,表现出发生反应的微观特征。对图5(b)方框内粒子进行EDS元素分析发现,其化学组成与NbSi2化学式的原子个数比相符,与XRD结果中出现NbSi2衍射峰一致。由此可见,与10%孔隙率样品相比,在相同飞片速度冲击下,随着冲击温度升高,高初始孔隙率反应产物中的铌硅颗粒发生严重变形并伴随硅发生熔融,并且NbSi2产物增多,表明发生了明显的冲击化学反应。
图6为高孔隙率(35%)铌硅粉末混合物冲击回收样品的SEM结果。由图6(a)可知,样品已转化为相对均匀的单一微观组织结构,表明该孔隙率样品在相同冲击速度下已发生完全反应形成新的相结构。图6(b)为图6(a)中方框内区域局部放大图,图中的新相为枝条状结构,其间有微纳孔(孔隙率约为4.7%,孔隙尺寸介于0.2~4 μm之间)。EDS能谱分析结果显示,铌和硅的原子比约为63.02:36.98,与Nb5Si3化学式的原子个数比接近,表明冲击反应生成物为Nb5Si3化合物。35%孔隙率样品中,冲击加载创造了极高的温度条件(估算冲击样品中的平均温度达2 256 K),相同冲击速度作用于铌硅粉末在物相生成和反应程度上与前述低孔隙率样品相比显著提高,并且其物相形态有明显不同。
值得注意的是,高孔隙率样品在飞片速度约为2.35 km/s时就发生了完全反应,而前期对不同冲击速度下反应行为的研究发现,较低孔隙率(20%)样品[13]需要更高的冲击速度(飞片速度约为2.63 km/s)才能发生完全反应。因此,高孔隙率样品崩塌所致的高能量贮集产生的高温是发生完全反应的重要条件。
2.3 冲击回收产物的DSC分析
图7为不同孔隙率铌硅粉末冲击产物的DSC曲线。由图7可知,孔隙率为10%和20%冲击产物的DSC曲线分别在720 ℃和810 ℃附近出现了两个放热峰。结合冲击产物的XRD结果分析可知,720 ℃放热峰可能是冲击产物中未反应的Nb和Si之间进一步发生反应放热所致。可能因为未发生完全化学反应的冲击产物中,未反应的铌和硅粒子在比自蔓延燃烧反应开始温度低约200 ℃时就开始发生化学反应[18]。该异常是由于冲击作用引起粒子粉末发生塑性流动和混合,表面清洁以及样品自身缺陷导致在冲击压缩加载期间粒子粉末被活化。而810 ℃放热峰可能是由于冲击塑变的Nb颗粒发生再结晶。35%孔隙度冲击回收产物的DSC结果呈现光滑的曲线,在测试加热范围内未出现放热峰,表明冲击回收产物为稳定的化合物。这与前述的物相表征和形貌分析结果一致,也证实了35%孔隙度样品完全反应生成了Nb5Si3金属间化合物。
粉末颗粒混合物的初始状态涵盖颗粒自身组成结构及混合颗粒间的相互接触等条件[19-22]。冲击加载时,在冲击波的高压作用下,样品颗粒的变化不是孤立发生的,而受相互间直接作用和颗粒自身物理性质的限制,表现为样品内颗粒的整体协调变形及交互作用。这些作用对于颗粒的变形方式、变形程度和融合都会产生直接影响。更重要的是,不同变化方式和程度会导致不同的能量集聚,从而产生温度升高和局域热效应[23],这些变化又反过来影响物质的扩散、相互渗透甚至激活化学反应。研究表明[24-25],合适的孔隙率对于粉末颗粒反应物的冲击化学合成非常重要。
由表征分析结果可知,孔隙率对于铌硅粉末混合物冲击反应的影响,是导致其反应发生及反应机制变化的重要因素。20%孔隙率冲击反应伴随着铌硅颗粒间强烈的相互作用、较低的相互约束以及大面积的孔隙塌陷或湮灭,冲击加热效果显著提升,不仅导致硅粒子发生熔化,变形的铌粒子嵌入到熔融的硅中,而且以铌粒子所在位置为形核中心,铌与周围的硅粒子发生半“固-液”化学反应生成相应的铌硅化合物,不过该激活程度仍处于仅能激发NbSi2反应物生成的条件。在更高孔隙率(35%)下,铌硅颗粒交互作用、变形及孔隙塌陷引起的能量集聚很高,使得样品内部急剧升温,引发高熔点的铌熔化,熔融的铌硅原子发生自由迁移,在极高温度激活下,铌硅间发生“液相”化学反应获得相应的铌硅化合物,即全面激活了铌硅间的自维持反应,并获得高熔点的Nb5Si3金属间化合物且反应完全。值得注意的是,过大孔隙率会在样品中留下微纳尺度孔隙,可能会对样品的块体性能产生影响。此外,实验研究发现,若孔隙率过高,孔隙崩塌产生的过高温度可能使冲击回收样品盒局部熔化,出现破漏或样品内气体相较多导致回收盒爆裂,从而导致冲击回收样品失败。
3. 结 论
通过3种不同初始孔隙率(10%、20%和35%)铌硅粉末混合物的冲击化学反应实验,对冲击回收产物物相、微观形貌表征分析及冲击产物热响应行为进行研究分析,发现不同孔隙率铌硅粉末混合物冲击反应的发生情况、产物特征及反应程度有所不同。低孔隙率(10%)样品几乎不发生反应;当孔隙率提高至20%时,铌硅粉末混合物发生不完全化学反应,生成NbSi2化合物;当孔隙率增大到35%时,相同冲击速度下铌硅粉末混合物发生完全反应,并获得与初始配比一致的Nb5Si3化合物。实验结果证实,适当增大铌硅粉末的初始孔隙率有利于发生冲击化学反应。
-
表 1 泡沫铝的主要力学性能参数
Table 1. Main mechanical properties of aluminum foam
Density/(g·cm−3) Elastic
modulus/MPaPlateau
stress/MPa0.37 70 4.7 0.54 145 9.2 0.75 250 19.0 表 2 钢板的主要力学性能参数
Table 2. Main mechanical parameters of steel
Density/(g·cm−3) Elastic modulus/GPa Yield strength/MPa Ultimate tensile strength/MPa 7.8 194 300 390 表 3 PLASTIC_KINEMATIC本构模型参数
Table 3. Constitutive model parameters of PLASTIC_KINEMATIC
Density/(g·cm−3) E/GPa 0/MPa Eh/MPa D/s−1 7.8 194 0.3 300 250 40.5 5 0.22 表 4 TNT炸药的状态方程参数
Table 4. Equation of state parameters of TNT
A/GPa B/GPa R1 R2 E01/(kJ·m−3) 371.2 3.231 4.15 0.95 0.3 7×106 表 5 研究工况
Table 5. Working conditions
Condition No. Explosive environment SoD/mm Method Core thickness/mm 1 Confined 50 Experiment/FEM 20 2 Confined 100 Experiment/FEM 20 3 Confined 150 Experiment/FEM 20 4 Confined 350 FEM 20 5 Confined 550 FEM 20 6 Vented 50 FEM 20 7 Vented 550 FEM 20 8 Free air burst 50 FEM 20 9 Free air burst 550 FEM 20 表 6 实验与数值模拟的背爆面板最大挠度对比
Table 6. Experimental and numerical comparison of the maximum deflection of the rear plate
SoD/mm Maximum deflection/mm Relative error/% Experiment Simulation 50 40 31 22.5 100 25 24 4.0 150 26 24 7.7 -
[1] LANGDON G S, YUEN S C K, NURICK G N. Experimental and numerical studies on the response of quadrangular stiffened plates. Part Ⅱ: localised blast loading [J]. International Journal of Impact Engineering, 2005, 31(1): 85–111. doi: 10.1016/j.ijimpeng.2003.09.050 [2] YUEN S C K, NURICK G N. Experimental and numerical studies on the response of quadrangular stiffened plates. Part Ⅰ: subjected to uniform blast load [J]. International Journal of Impact Engineering, 2005, 31(1): 55–83. doi: 10.1016/j.ijimpeng.2003.09.048 [3] 侯海量, 朱锡, 古美邦. 爆炸载荷作用下加筋板的失效模式分析及结构优化设计 [J]. 爆炸与冲击, 2007, 27(1): 26–33. doi: 10.11883/1001-1455(2007)01-0026-08HOU H L, ZHU X, GU M B. Study on failure mode of stiffened plate and optimized design of structure subjected to blast load [J]. Explosion and Shock Waves, 2007, 27(1): 26–33. doi: 10.11883/1001-1455(2007)01-0026-08 [4] 陈鹏宇. 舱内爆炸载荷及夹芯复合舱壁结构动响应特性研究 [D]. 武汉: 海军工程大学, 2019: 85−90.CHEN P Y. Study on characteristics of confined explosion load and dynamic response of sandwich composite bulkhead structure [D]. Wuhan: Naval University of Engineering, 2019: 85−90. [5] 焦立启, 侯海量, 陈鹏宇, 等. 爆炸冲击载荷下固支单向加筋板的动响应及破损特性研究 [J]. 兵工学报, 2019, 40(3): 592–600. doi: 10.3969/j.issn.1000-1093.2019.03.019JIAO L Q, HOU H L, CHEN P Y, et al. Research on dynamic response and damage characteristics of fixed supported one-way stiffened plates under blast loading [J]. Acta Armamentarii, 2019, 40(3): 592–600. doi: 10.3969/j.issn.1000-1093.2019.03.019 [6] 肖锋, 谌勇, 章振华, 等. 三明治结构爆炸冲击动力学研究综述 [J]. 噪声与振动控制, 2012, 32(6): 1–7.XIAO F, CHEN Y, ZHANG Z H, et al. Review of research of blast impact of sandwich structures [J]. Noise and Vibration Control, 2012, 32(6): 1–7. [7] 邓旭辉, 李亚斌. 双层泡沫铝夹芯板抗冲击性能数值研究 [J]. 铁道科学与工程学报, 2019, 16(10): 2603–2611. doi: 10.19713/j.cnki.43-1423/u.2019.10.029DENG X H, LI Y B. Numerical study on impact resistance of double-layer foam aluminum sandwich panels [J]. Journal of Railway Science and Engineering, 2019, 16(10): 2603–2611. doi: 10.19713/j.cnki.43-1423/u.2019.10.029 [8] 王涛, 余文力, 秦庆华, 等. 爆炸载荷下泡沫铝夹芯板变形与破坏模式的实验研究 [J]. 兵工学报, 2016, 37(8): 1456–1463. doi: 10.3969/j.issn.1000-1093.2016.08.017WANG T, YU W L, QIN Q H, et al. Experimental investigation into deformation and damage patterns of sandwich plates with aluminum foam core subjected to blast loading [J]. Acta Armamentarii, 2016, 37(8): 1456–1463. doi: 10.3969/j.issn.1000-1093.2016.08.017 [9] 韩守红, 吕振华. 铝泡沫夹层结构抗爆炸性能仿真分析及优化 [J]. 兵工学报, 2010, 31(11): 1468–1474.HAN S H, LÜ Z H. Numerical simulation of blast-resistant performance of aluminum foam sandwich structures and optimization [J]. Acta Armamentarii, 2010, 31(11): 1468–1474. [10] 亓昌, 杨丽君, 杨姝. 梯度铝泡沫夹层结构抗爆性能仿真与优化 [J]. 振动与冲击, 2013, 32(13): 70–75. doi: 10.3969/j.issn.1000-3835.2013.13.014QI C, YANG L J, YANG S. Simulation and optimization for blast-resistant performances of a graded aluminum foam sandwich structure [J]. Journal of Vibration and Shock, 2013, 32(13): 70–75. doi: 10.3969/j.issn.1000-3835.2013.13.014 [11] LI S Q, LI X, WANG Z H, et al. Finite element analysis of sandwich panels with stepwise graded aluminum honeycomb cores under blast loading [J]. Composites Part A: Applied Science and Manufacturing, 2016, 80: 1–12. doi: 10.1016/j.compositesa.2015.09.025 [12] 李春鹏, 张攀, 刘均, 等. 空爆载荷下功能梯度泡沫铝夹层板动响应数值仿真 [J]. 中国舰船研究, 2018, 13(3): 77–84. doi: 10.19693/j.issn.1673-3185.01172LI C P, ZHANG P, LIU J, et al. Numerical simulation of dynamic response of functionally graded aluminum foam sandwich panels under air blast loading [J]. Chinese Journal of Ship Research, 2018, 13(3): 77–84. doi: 10.19693/j.issn.1673-3185.01172 [13] Unified Facilities Criteria. Structures to resist the effects of accidental explosions: UFC 3−340−02 [S]. Washington D C, USA: Department of the Army, the Navy and the Air Force, 2008. [14] 侯海量, 朱锡, 李伟, 等. 舱内爆炸冲击载荷特性实验研究 [J]. 船舶力学, 2010, 14(8): 901–907. doi: 10.3969/j.issn.1007-7294.2010.08.011HOU H L, ZHU X, LI W, et al. Experimental studies on characteristics of blast loading when exploded inside ship cabin [J]. Journal of Ship Mechanics, 2010, 14(8): 901–907. doi: 10.3969/j.issn.1007-7294.2010.08.011 [15] 侯海量, 朱锡, 梅志远. 舱内爆炸载荷及舱室板架结构的失效模式分析 [J]. 爆炸与冲击, 2007, 27(2): 151–158. doi: 10.11883/1001-1455(2007)02-0151-08HOU H L, ZHU X, MEI Z Y. Study on the blast load and failure mode of ship structure subject to internal explosion [J]. Explosion and Shock Waves, 2007, 27(2): 151–158. doi: 10.11883/1001-1455(2007)02-0151-08 [16] ZHANG J, ZHAO G P, LU T J. Dynamic responses of sandwich beams with gradient-density aluminum foam cores [J]. International Journal of Protective Structures, 2011, 2(4): 439–451. doi: 10.1260/2041-4196.2.4.439 [17] SYMONDS P S, YU T X. Counterintuitive behavior in a problem of elastic-plastic beam dynamics [J]. Journal of Applied Mechanics, 1985, 52(3): 517–522. doi: 10.1115/1.3169093 [18] SHANLEY F R. Inelastic column theory [J]. Journal of the Aeronautical Sciences, 1947, 14(5): 261–268. doi: 10.2514/8.1346 [19] YU T X, XU Y. The anomalous response of an elastic-plastic structural model to impulsive loading [J]. Journal of Applied Mechanics, 1989, 56(4): 868–873. doi: 10.1115/1.3176184 -