Failure modes and response characteristics of finite-thickness aluminum targets under normal penetration of elliptical cross-section projectiles
-
摘要: 基于30 mm口径弹道炮平台,开展了3种不同椭圆横截面弹体在200~600 m/s撞击速度范围内正侵彻2A12铝靶的实验,获得了2A12铝靶的破坏形貌及弹体的剩余速度。在此基础上,建立了相应的数值模型,结合实验结果验证了所建模型的有效性,并系统分析了弹体横截面长短轴长度比对靶体的破坏情况及响应特性的影响。研究结果表明:弹体最大横截面面积是影响弹体剩余速度的主要因素,而弹体横截面长短轴长度比对弹体剩余速度的影响较弱;在圆形横截面弹体侵彻下靶体背部形成的花瓣大小和形状一致,空间分布均匀,而在椭圆横截面弹体侵彻下,随着弹体横截面长短轴长度比的增大,靶体背部形成的花瓣数量增加、尺寸变小,且在短轴方向的花瓣数量和靶体表面隆起高度均大于长轴方向的;靶体在圆形横截面弹体侵彻下的径向位移、径向应力和切向应力与其在椭圆横截面弹体侵彻下的显著不同,前者沿周向方向各点的变化规律基本一致,靶体处于简单的压缩状态,切向应力为零,而后者各点的应力状态与弹体横截面长短轴长度比和周向角密切相关,靶体受到压缩和剪切应力的耦合作用。Abstract: By means of a 30-mm-caliber ballistic gun platform, a series of experiments were carried out on 2A12 aluminum targets subjected to normal penetration by three kinds of 30CrMnSi2A steel projectiles with different elliptical cross-section shapes in the striking velocity range from 200 m/s to 600 m/s. The residual velocities of the projectiles and the failure modes of the targets were experimentally obtained. Based on the experimental results, the corresponding numerical models were established and verified. And the influences of the major-to-minor axis length ratios of the projectile cross-sections on the failure modes and response characteristics of the targets were systematically analyzed. The results show as follows. The maximum cross-sectional areas of the projectiles are the main factor affecting the residual velocities of the projectiles, while the major-to-minor axis length ratios of the projectile cross-sections have little effect on the residual velocities. Therefore, in engineering applications, the engineering model for the circular cross-section projectile penetrating a target can be directly used to calculate the residual velocity of the elliptical cross-section projectile with the same maximum cross-sectional area. In addition, under normal penetration of the circular cross-section projectiles, the sizes, shapes and distribution of the petals induced at the back faces of the targets are uniform. However, under normal penetration of the elliptical cross-section projectiles, as the major-to-minor axis length ratios of the projectile cross-sections increase, the numbers of the petals induced at the back faces of the targets increase and the petal sizes decrease, and the petal numbers and the uplifted height in the minor axis direction are greater than those in the major axis direction. The radial displacement, radial stress and tangential stress of the targets under the normal penetration of the elliptical cross-section projectiles are obviously different from those of the targets under the normal penetration of the circular cross-section projectiles. Under normal penetrations of the circular cross-section projectiles, the above response characteristics of the targets change basically the same along the circumferential directions and the targets are under simple compression states with the tangential stress of zero. But, under normal penetrations of the elliptical cross-section projectiles, the stress states of different points of the targets are closely related to the major-to-minor axis length ratios and the circumferential angles of the projectiles, and the targets are subjected to the coupling effects of the compression and shear stresses.
-
可燃气体在化工过程及能源发展方面均有广泛的运用。研究可燃气体的燃烧爆炸特性对于火灾爆炸事故的评估及安全生产都具有重要的意义,尤其是在不同环境条件下的燃烧爆炸特性,例如高温高压条件下。能够维持可燃气体火焰传播的最低/最高浓度被称为可燃气体的爆炸下限(lower explosion limit, LEL)/爆炸上限(upper explosion limit, UEL)。燃料的爆炸极限是预测火灾、评估爆炸可能性和设计保护系统的重要依据[1]。尽管有很多不同的数值方法可以确定可燃气体的爆炸极限,例如一维平面模拟、极限火焰温度和极限燃烧速度的应用,但是实验测定在获得特定参数及验证数值结果可靠性方面仍然非常必要[2-4]。爆炸极限随初始条件(例如初始温度、初始压力)的不同而改变。学者们对常温常压下气体及蒸汽的爆炸极限已经进行了一些研究[5-7]。通常认为,常温常压下甲烷在空气中的爆炸下限为5%,上限为15%。然而,高温高压条件下甲烷-空气混合物爆炸极限的研究非常有限[6, 8-9],而且温度压力对爆炸极限共同影响的研究仅涉及爆炸上限[10-13]。本文中借助特殊环境20 L爆炸特性测试系统,测定不同初始温度(25~200 ℃)和初始压力(0.1~1.0 MPa)条件下甲烷-空气混合物的爆炸上限及下限,研究初始温度和压力对甲烷-空气混合物爆炸极限的耦合影响规律。以期研究结果不仅对煤层气的安全利用及煤矿的安全生产具有指导意义,也为进一步研究多因素耦合影响作用提供参考。
1. 实验系统与方法
采用特殊环境20 L爆炸特性测试系统进行高温高压条件下甲烷-空气混合物爆炸极限的实验测定,该实验系统工作原理如图 1所示。实验系统主要包括容积为20 L的爆炸罐体、配气系统、抽真空系统、点火系统、加热系统、控制系统和采集系统7部分。爆炸罐体为不锈钢双层球形结构,配气系统根据体积分压原则由电磁阀控制配制实验设定浓度的实验气体,该系统配气精度为±0.1%,实验中运用高压电点火,点火能量为10 J。实验中先对罐体进行抽真空,根据设定的初始压力及甲烷浓度,配气系统按体积分压依次充入甲烷和空气。若进行常温下试验,则静置5 min左右以使实验气体混合均匀,若进行高温试验,则开启加热系统加热实验气体至实验温度。之后借助罐体中心处的点火电极进行点火,设置200 ms的延迟时间以使点火系统在更稳定的状态下进行点火操作。
依据ASTM标准[14],以最大爆炸压力超过初始压力7%的压力升高作为爆炸极限实验的爆炸判据,若连续3次实验均未爆炸,则认为此时的甲烷浓度为此条件下的爆炸极限值。为研究初始温度和初始压力对甲烷爆炸极限的影响,罐体内实验气体温度变化范围为25~200 ℃,压力变化范围为0.1~1.0 MPa。点火前罐体内的实验气体均处于静止状态。
2. 温度、压力单因素的影响
图 2为甲烷-空气混合物爆炸极限随初始温度和初始压力的变化情况。随着初始温度的升高和初始压力的增大,甲烷-空气混合物爆炸上限升高,爆炸下限降低,爆炸极限范围扩大。
由图 2(a)可知:在常压条件下,当初始温度为25 ℃时爆炸上限为15.8%,下限为5.1%,极限范围为10.7%;而当初始温度升高至200 ℃时爆炸上限上升为17.7%,下限下降为4.6%,极限范围扩大为13.1%,爆炸极限扩大的百分率为22.4%。由图 2(b)可知:在常温条件下,当初始压力由0.1 MPa增加到1.0 MPa,甲烷-空气混合物爆炸上限上升到21.4%,下限降低到4.5%,极限范围扩到到16.9%,爆炸极限较常温常压时扩大的百分率为57.9%。在实验温度和压力范围内,甲烷-空气混合物的爆炸上限和下限与初始温度和初始压力基本呈线性关系:
y=A+Bx (1) 式中:当x=T0/℃时,25≤x≤200;当x=p0/MPa时,0.1≤x≤1.0。各参数值如表 1所示。
表 1 拟合函数的参数Table 1. Parameters for fitting functiony% x A B R UEL T0/℃ 15.464 3 0.011 4 0.995 5 UEL p0/MPa 15.659 7 6.142 4 0.987 6 LEL T0/℃ 5.142 9 -0.002 5 -0.996 6 LEL p0/MPa 5.151 2 -0.666 9 -0.998 5 通过比较各拟合直线的斜率可以看出,初始温度和压力对甲烷-空气混合物爆炸上限的影响要大于其对下限的影响。这是由于在爆炸下限浓度附近,可燃混合气体中甲烷的含量很小,过量的空气作为反应中的惰性气体,一方面阻碍了甲烷与氧气分子的有效碰撞,减少了反应发生的可能性;另一方面吸收反应放热,不利于链式反应的持续及火焰的蔓延。因此,温度的升高和压力的增大虽然使甲烷爆炸下限降低,但由于过量空气的存在,其受影响的程度要小一些。
3. 温度、压力的耦合影响
3.1 甲烷-空气混合物爆炸上限
图 3为压力温度耦合影响下甲烷-空气混合物爆炸上限的变化情况。在不同的初始压力条件下,爆炸上限均随初始温度的上升而升高,变化趋势基本相同,如图 3(a)所示。当处于相同的初始温度时,随着初始压力的升高,单位压力对爆炸上限的影响逐渐减弱。由图 3(b)可知,在不同初始温度条件下,爆炸上限均随初始压力的上升而升高。从变化趋势来看,压力的升高使单位温度升高对爆炸上限产生更大的影响。在相同的初始压力下,随着初始温度的升高,单位温度升高对爆炸上限的影响逐渐减弱。这是由于随着甲烷-空气混合物爆炸上限的升高,反应系统中的氧气含量逐渐减少,处于负氧状态的系统发展成为爆炸系统将更加困难,这就对系统提出了更高的要求,因此,单位温度和单位压力的升高对瓦斯爆炸上限的影响逐渐变小。
常温常压条件下,甲烷-空气混合物爆炸上限为15.8%,当初始温度升高到200 ℃且初始压力增大到1.0 MPa时,爆炸上限上升为25.7%,上升幅度几乎达到了62.7%。常压下,当初始温度由25 ℃升至200 ℃时,上限的升幅仅为8.2%;而常温下,当初始压力由0.1 MPa增大到1.0 MPa时,上限升幅也只为35.4%。由此可见,甲烷-空气混合物爆炸上限在初始温度和初始压力耦合影响作用下的变化幅度比单一因素影响下的变化幅度大得多,且远大于单一影响因素下变化幅度的加和。为了综合分析初始温度和初始压力对甲烷爆炸上限的耦合影响,以初始温度作为x轴、初始压力作为y轴、甲烷-空气混合物爆炸上限作为z轴得到爆炸上限随初始温度和初始压力的变化曲面,如图 4所示。该曲面更直观地反映了初始温度和初始压力对甲烷-空气混合物爆炸上限的耦合影响作用,根据拟合函数:
z=z0−12Aexp[(x−xcw1)+(y−ycw2)2] (2) 可预估实验温度压力范围内的甲烷爆炸上限, 式(2)中50≤x≤200, 0.2≤y≤1.0, 拟合函数各参数如表 2所示。
表 2 拟合函数的参数Table 2. Parameters for fitting functionz x y z0 A UEL/% T0/℃ p0/MPa 13.80 12.30 xc w1 yc w2 R2 263.81 267.86 1.13 0.75 0.998 9 3.2 甲烷-空气混合物爆炸下限
甲烷-空气混合物的爆炸下限随初始温度和初始压力耦合影响的变化情况如图 5所示。在不同初始压力条件下,爆炸下限均随初始温度的上升而降低;在不同初始温度条件下,下限均随初始压力的增大而降低,且变化趋势均具有较好的一致性。甲烷-空气混合物常温常压时的爆炸下限为5.1%,当初始条件改变为200 ℃、1.0 MPa时爆炸下限下降为4.1%,下降幅度为19.6%。而在相同实验条件下,仅初始温度升高至200 ℃时,爆炸下限下降幅度为9.8%,而仅初始压力增大至1.0 MPa时下降幅度为11.8%。由此可见,甲烷-空气混合物爆炸下限在初始温度和初始压力耦合影响作用下的变化幅度也比单一因素影响下的变化幅度要大,但其与两者单一影响下的变化幅度的加和基本一致。这是因为在爆炸下限浓度附近,空气的惰性效应是反应进行中的重要阻碍,初始温度的升高和初始压力的增大虽然在一定程度上增加了甲烷分子与氧气分子的有效碰撞,提高了反应速率,有利于链式反应的发展,但是由于甲烷浓度的限制,两者耦合的正反馈效应并不是很明显,因而温度压力对甲烷爆炸下限的影响基本上是两者单一影响效果的叠加。图 6为初始温度和初始压力对甲烷-空气混合物爆炸下限的耦合影响曲面,拟合曲面的函数表达式为:
z=z0+Aexp(−xB)exp(−yC)50≤x≤200,0.2≤y≤1.0 (3) 根据此式(3)可以求得实验温度压力范围内的甲烷-空气混合物爆炸下限值, 拟合函数各参数见表 3。
表 3 拟合函数的参数Table 3. Parameters for fitting functionz x y z0 LEL/% T0/℃ p0/MPa 3.08 A B C R2 2.34 587.71 1.81 0.995 3 初始条件为200 ℃、1.0 MPa时甲烷-空气混合物的爆炸极限范围为21.6%,比常温常压时的10.7%扩大的百分率为101.9%,扩大的危险浓度范围使甲烷-空气混合气体在高温高压条件下的危险性也有了很大程度的提高,因此,如果甲烷-空气混合气体有处于高温高压条件的可能性,则在生产流程及安全设施上需要对更宽浓度范围的可燃气体进行监管和治理。
4. 结论
随着初始温度或初始压力的升高,甲烷-空气混合物的爆炸下限降低,爆炸上限升高,爆炸极限范围扩大。在常温常压下,爆炸上限和下限与初始压力、初始温度近似呈线性关系。
在相同的初始温度条件下,随着初始压力的升高,单位压力升高对甲烷-空气混合物爆炸上限的影响逐渐减弱。当初始压力相同时,单位温度升高对爆炸上限的影响随初始温度的升高而减弱。当初始条件改变时,单位温度/压力的升高对爆炸下限的影响变化并不显著。
初始温度和初始压力对甲烷-空气混合物爆炸上限和下限的耦合影响比单因素的影响要大得多。两者对爆炸上限的耦合影响远大于两者单一影响的加和,而对爆炸下限的影响几乎等同于两者单一影响的加和。温度压力与爆炸上限和下限的关系均可用特定的曲面方程进行描述。
为了深入了解爆炸极限与初始条件之间的依赖关系,还需要继续研究更高的温度和压力条件下甲烷-空气混合物等可燃气体的爆炸极限以及其他其他影响因素之间的耦合影响作用。
-
表 1 三种弹体主要参数
Table 1. Main parameters of three projectiles
弹体类型 弹体轮廓 2a/mm 2b/mm β L/mm ψ m/g C1 23.6 23.6 1.00 43.2 3.6 360 T1 30.0 18.6 1.61 43.2 5.6 360 T2 30.0 24.0 1.25 43.2 3.5 360 表 2 弹体正侵彻铝靶的实验结果
Table 2. Experimental results for normal penetration of projectiles into aluminum targets
弹体 v0/(m·s−1) vr/(m·s−1) α/(°) γ/(°) C1-1 402.3 336.0 1.08 +0.63 C1-2 310.7 214.5 0.34 +0.19 C1-3 256.2 120.3 0.25 +0.17 C1-4 566.5 522.5 1.11 −0.68 T1-1 402.0 338.0 1.26 −0.62 T1-2 229.4 0 2.03 −1.51 T1-3 570.3 531.6 1.91 −1.49 T2-1 405.3 322.7 0.76 −0.30 T2-2 229.5 0 2.69 −2.39 T2-3 569.3 509.2 1.57 −1.13 表 3 材料参数
Table 3. Material parameters
表 4 弹体剩余速度模拟结果与实验结果的对比
Table 4. Comparison of residual velocities of projectiles between simulation and experiment
弹体类型 v0/(m·s−1) vr/(m·s−1) ɛr/% 实验 数值模拟 C1 256.2 120.3 109.2 −9.2 310.7 214.5 202.7 −5.5 402.3 336.0 323.3 −3.8 566.5 522.5 510.7 −2.3 T1 229.4 0 21.3 402.0 338.0 323.5 −4.3 570.3 531.6 514.5 −3.2 T2 229.5 0 0 0 405.3 322.7 302.5 −6.3 569.3 509.2 497.5 −2.3 表 5 靶体背部塑性应变区域范围对比
Table 5. Comparison of plastic strain ranges on the back of targets
弹体类型 2a/mm 2b/mm 长轴最大坐标/mm 短轴最大坐标/mm 长轴相对增量/% 短轴相对增量/% C1 23.60 23.60 18.55 18.55 57.20 57.20 T3 26.40 21.10 17.30 18.08 31.06 71.37 T1 30.00 18.60 18.13 20.10 20.87 116.13 T4 33.36 16.68 19.13 22.30 14.69 167.39 -
[1] VEDERNIKOV Y A, LEVIN V A, KHUDYAKOV Y S. Evolution and comparative analysis of group armor- and aeroballistics of ancient and modern ruled and poly-wedge arrows [J]. Bulletin of the Novosibirsk Computing Center. Series: Mathematical Modeling in Geophysics, 2005, 10: 93–116. [2] DAI X H, WANG K H, LI M R, et al. Rigid elliptical cross-section ogive-nose projectiles penetration into concrete targets [J]. Defence Technology, 2021, 17(3): 800–811. DOI: 10.1016/j.dt.2020.05.011. [3] MA X H, ZHANG Q M, ZHANG X W. A model for rigid asymmetric ellipsoidal projectiles penetrating into metal plates [J]. International Journal of Impact Engineering, 2022, 163: 104140. DOI: 10.1016/j.ijimpeng.2021.104140. [4] BEN-DOR G, DUBINSKY A, ELPERIN T. Optimal 3D impactors penetrating into layered targets [J]. Theoretical and Applied Fracture Mechanics, 1997, 27(3): 161–166. DOI: 10.1016/S0167-8442(97)00018-9. [5] BEN-DOR G, DUBINSKY A, ELPERIN T. A model for predicting penetration and perforation of FRP laminates by 3-D impactors [J]. Composite Structures, 2002, 56(3): 243–248. DOI: 10.1016/S0263-8223(02)00009-0. [6] YAKUNINA G Y. The construction of optimum three-dimensional shapes within the framework of a model of local interaction [J]. Journal of Applied Mathematics and Mechanics, 2000, 64(2): 289–298. DOI: 10.1016/S0021-8928(00)00051-4. [7] YAKUNINA G Y. The optimum non-conical and asymmetrical three-dimensional configurations [J]. Journal of Applied Mathematics and Mechanics, 2000, 64(4): 583–591. DOI: 10.1016/S0021-8928(00)00084-8. [8] WOO H J. Cavity expansion analysis of non-circular cross-sectional penetration problems [D]. Austin, Texas, USA: The University of Texas at Austin, 1997: 132–154. [9] BLESS S J. Penetration mechanics of non-circular rods [J]. AIP Conference Proceedings, 1996, 370(1): 1119–1122. DOI: 10.1063/1.50861. [10] 杜忠华, 曾国强, 余春祥, 等. 异型侵彻体垂直侵彻半无限靶板试验研究 [J]. 弹道学报, 2008, 20(1): 19–21.DU Z H, ZENG G Q, YU C X, et al. Experimental research of novel penetrator vertically penetrating semi-infinite target [J]. Journal of Ballistics, 2008, 20(1): 19–21. [11] 杜忠华, 朱建生, 王贤治, 等. 异型侵彻体垂直侵彻半无限靶板的分析模型 [J]. 兵工学报, 2009, 30(4): 403–407. DOI: 10.3321/j.issn:1000-1093.2009.04.005.DU Z H, ZHU J S, WANG X Z, et al. Analytical model on non-circular penetrator impacting semi-infinite target perpendicularly [J]. Acta Armamentarii, 2009, 30(4): 403–407. DOI: 10.3321/j.issn:1000-1093.2009.04.005. [12] 高光发, 李永池, 刘卫国, 等. 长杆弹截面形状对垂直侵彻深度的影响 [J]. 兵器材料科学与工程, 2011, 34(3): 5–8. DOI: 10.3969/j.issn.1004-244X.2011.03.002.GAO G F, LI Y C, LIU W G, et al. Influence of the cross-section shapes of long rod projectile on the vertical penetration depth [J]. Ordnance Material Science and Engineering, 2011, 34(3): 5–8. DOI: 10.3969/j.issn.1004-244X.2011.03.002. [13] DONG H, LIU Z H, WU H J, et al. Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete [J]. International Journal of Impact Engineering, 2019, 132: 103311. DOI: 10.1016/j.ijimpeng.2019.05.025. [14] DONG H, WU H J, LIU Z H, et al. Penetration characteristics of pyramidal projectile into concrete target [J]. International Journal of Impact Engineering, 2020, 143: 103583. DOI: 10.1016/j.ijimpeng.2020.103583. [15] GAO X D, LI Q M. Trajectory instability and convergence of the curvilinear motion of a hard projectile in deep penetration [J]. International Journal of Mechanical Sciences, 2017, 121: 123–142. DOI: 10.1016/j.ijmecsci.2016.12.021. [16] 高旭东, 李庆明. 带攻角斜侵彻混凝土的弹道偏转分析 [J]. 兵工学报, 2014, 35(Supp1 2): 33–39.GAO X D, LI Q M. Trajectory analysis of projectile obliquely penetrating into concrete target at attack angle [J]. Acta Armamentarii, 2014, 35(Supp1 2): 33–39. [17] 王浩, 武海军, 闫雷, 等. 椭圆横截面弹体斜贯穿双层间隔薄钢板失效模式 [J]. 兵工学报, 2020, 41(Suppl 2): 1–11. DOI: 10.3969/j.issn.1000-1093.2020.S2.001.WANG H, WU H J, YAN L, et al. Failure mode of oblique perforation of truncated ogive-nosed projectiles with elliptic cross-section into double-layered thin steel plate with gap space [J]. Acta Armamentarii, 2020, 41(Suppl 2): 1–11. DOI: 10.3969/j.issn.1000-1093.2020.S2.001. [18] 王浩, 潘鑫, 武海军, 等. 椭圆截面截卵形刚性弹体正贯穿加筋板能量耗散分析 [J]. 爆炸与冲击, 2019, 39(10): 103203. DOI: 10.11883/bzycj-2018-0350.WANG H, PAN X, WU H J, et al. Energy dissipation analysis of elliptical truncated oval rigid projectile penetrating stiffened plate [J]. Explosion and Shock Waves, 2019, 39(10): 103203. DOI: 10.11883/bzycj-2018-0350. [19] LANDKOF B, GOLDSMITH W. Petalling of thin, metallic plates during penetration by cylindro-conical projectiles [J]. International Journal of Solids and Structures, 1985, 21(3): 245–266. DOI: 10.1016/0020-7683(85)90021-6. [20] 张中国, 黄风雷, 段卓平, 等. 弹体侵彻带加强筋结构靶的实验研究 [J]. 爆炸与冲击, 2004, 24(5): 431–436. DOI: 10.3321/j.issn:1001-1455.2004.05.009.ZHANG Z G, HUANG F L, DUAN Z P, et al. The experimental research for projectile penetrating the structural target with rebar [J]. Explosion and Shock Waves, 2004, 24(5): 431–436. DOI: 10.3321/j.issn:1001-1455.2004.05.009. [21] CHEN Y, WANG Y, TANG P, et al. Impact characteristics of stiffened plates penetrated by sub-ordnance velocity projectiles [J]. Journal of Constructional Steel Research, 2008, 64(6): 634–643. DOI: 10.1016/j.jcsr.2007.12.006. [22] SONG W D, NING J G, WANG J. Normal impact of truncated oval-nosed projectiles on stiffened plates [J]. International Journal of Impact Engineering, 2008, 35(9): 1022–1034. DOI: 10.1016/j.ijimpeng.2007.05.008. [23] 徐双喜, 吴卫国, 李晓彬, 等. 截锥形弹穿甲单加筋板的破坏特性 [J]. 爆炸与冲击, 2011, 31(1): 62–68. DOI: 10.11883/1001-1455(2011)01-0062-07.XU S X, WU W G, LI X B, et al. Falure characteristics of a conical projectile penetrating single stiffened plate [J]. Explosion and Shock Waves, 2011, 31(1): 62–68. DOI: 10.11883/1001-1455(2011)01-0062-07. [24] LIU J W, ZHANG X F, WEI H Y, et al. Study on the penetration of elliptical cross-section projectiles into concrete targets: theory and experiment [J]. Latin American Journal of Solids and Structures, 2022, 19(3): 23. DOI: 10.1590/1679-78256939. [25] 李磊, 张先锋, 吴雪, 等. 不同硬度30CrMnSiNi2A钢的动态本构与损伤参数 [J]. 高压物理学报, 2017, 31(3): 239–248. DOI: 10.11858/gywlxb.2017.03.005.LI L, ZHANG X F, WU X, et al. Dynamic constitutive and damage parameters of 30CrMnSiNi2A steel with different hardnesses [J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 239–248. DOI: 10.11858/gywlxb.2017.03.005. [26] 米双山, 张锡恩, 陶贵明. 钨球侵彻LY-12铝合金靶板的有限元分析 [J]. 爆炸与冲击, 2005, 25(5): 477–480. DOI: 10.11883/1001-1455(2005)05-0477-04.MI S S, ZHANG X E, TAO G M. Finite element analysis of spherical fragments penetrating LY-12 aluminum alloy target [J]. Explosion and Shock Waves, 2005, 25(5): 477–480. DOI: 10.11883/1001-1455(2005)05-0477-04. -