Dynamic Mechanical Behavior of Graded Metallic Foams Based on Lagrangian Analysis Method
-
摘要: 采用Lagrangian分析法，对梯度泡沫金属在高速冲击下的变形机理和应力响应进行研究。基于3D-Voronoi技术，构建了5种不同密度梯度的泡沫金属细观有限元模型，并进行了高速冲击下的Taylor数值实验，得到不同密度梯度泡沫金属的质点速度分布规律。采用Lagrangian分析法并结合数值实验结果，研究了高速冲击下密度梯度参数对泡沫金属的局部应变分布、应力分布以及冲击波传播与衰减规律的影响。结果表明：负密度梯度泡沫金属比正密度梯度泡沫金属具有更强的抵抗变形能力，且密度梯度参数越小，变形程度越小；负密度梯度泡沫金属的局部压实应力呈线性减小，最大局部压实应力随着密度梯度参数的减小而增大，在冲击端附近可以承受更大的载荷；正密度梯度泡沫金属的局部压实应力分布呈平台状，其最大局部压实应力小于负密度梯度泡沫金属。
-
关键词:
- Lagrangian分析法 /
- 梯度泡沫金属 /
- 3D-Voronoi技术 /
- 密度梯度参数
Abstract: The Lagrangian analysis method was employed to investigate the deformation mechanism and stress response of graded metallic foams. The mesoscopic finite element models of the graded metallic foams with five different density gradient parameters were constructed by the 3D-Voronoi technique, and the corresponding Taylor numerical tests were performed under high-speed impact, and the particle velocity distributions of different graded foams were obtained. By combining the Lagrangian analysis method with the results of Taylor numerical tests, the effects of density gradient parameters on the local strain distribution, stress distribution, shock wave propagation and attenuation of metallic foams under high-speed impact were investigated. The results show that the metallic foams with negative density gradient have better resistance to vertical deformation than those with positive density gradient, and the deformation degree decreases with the decrease of the density gradient parameter. The local densification stress distribution of the metallic foams with negative density gradient decreases linearly, and the maximum local densification stress increases with the decrease of the density gradient parameter. The metallic foams with negative density gradient have high load bearing capability near the impact end. The local densification stress distribution of the metallic foams with positive density gradient has a plateau stage, and the maximum local densification stress is less than metallic foams with negative density gradient. -
为应对日益复杂的战场环境，需要不断创新毁伤模式。横向效应增强体(Penetrator with Enhanced Lateral Effect, PELE)是一种新型侵彻弹，最早由法德圣路易斯研究所于1996年提出[1]。PELE弹的典型结构特点是采用密度不同的外壳和内芯，无需引信和装药，其中：外壳一般为高密度材料，如合金钢、钨合金；弹芯为低密度材料，如铝、尼龙、聚乙烯等。弹丸侵彻目标时，由于低密度弹芯的侵彻能力弱，弹芯被压缩产生径向膨胀，弹体内部压力迅速上升；弹丸穿过目标后，内外压力差使壳体断裂，形成破片。由此可见，PELE弹兼具穿甲和破片杀伤效果，具有结构简单、低成本、弹体使用安全等优点，可出色地完成高效毁伤[2-3]。
传统侵彻弹的毁伤机理和作用特点研究已经相当成熟[4]，而PELE弹不同于一般的侵彻弹，需要在靶后形成有效破片，其弹体结构的侵彻规律也有自己的特点，因此学者们开展了大量的研究工作。张谋等[5]通过数值仿真技术分析了PELE弹内芯与横向效应的关联；朱建生等[6]对PELE弹的破碎机理和弹体结构进行了理论分析；杜忠华等[7]通过理论分析与数值仿真，认为PELE弹撞击金属薄板时，壳体的横向速度随着弹丸着靶速度和装填材料声阻抗的增大而增大。
在现代战场上，随着目标防护性能的增强，单一功能的弹丸很难达到理想的破坏效果，这就要求进攻模式多元化。为此，本研究提出分段式PELE弹的概念，期望通过分段式PELE弹体结构增强毁伤能力，尤其是应对多层目标。然而，分段式PELE弹的穿靶过程受多种因素影响，本研究采用数值模拟方法重点探讨弹体的侵彻速度和靶板厚度对其终点效应的影响规律，以实现分段式PELE弹应对不同目标时发挥出最强的优势性能。
1. 数值方法与模型
利用非线性动力学软件LS-DYNA进行数值模拟。弹体采用两段PELE结构，用螺纹连接，其结构如图 1所示，有限元模型见图 2。第1段和第2段弹体的长度分别为32和40 mm，外径10 mm，弹芯尺寸为Ø6 mm×30 mm。靶板为4层，每层靶的材料和厚度均相同，设靠近弹丸一侧为第1层。弹壳均采用钨合金，内芯材料为铝，弹尾采用钨合金。弹体材料参数列于表 1，其中ρ0为密度，E为弹性模量，G为剪切模量，ν为泊松比。
表 1 材料参数Table 1. Material parametersMaterial ρ0/(g·cm-3) E/GPa G/GPa ν Constitutive model Tungsten alloy 17.67 354 138 0.28 Johnson-Cook Aluminum 2.70 69 0.33 Plastic_Kinematic 921 steel 7.85 210 78 0.30 Johnson-Cook 为了提高计算效率，采用1/4模型模拟，约束对称面的位移和转动，设置为轴对称状态。靶板周围边界采用固定约束，限制所有方向的运动。网格尺寸为0.1~0.5 mm，弹体经过区域网格加密，其他区域采用渐变网格。网格单元均采用SOLID164单元，为八节点六面体[8]。靶板和壳体均采用Johnson-Cook模型和Grüneisen状态方程描述，弹芯采用Plastic_Kinematic模型描述。同时添加Add_Erosion，以控制材料的失效。弹体内部采用Contact_Automatic_Surface_To_Surface接触算法，弹体与靶板之间采用Contact_Eroding_Surface_To_Surface接触算法。
2. 模拟结果与分析
2.1 分段PELE弹与一般PELE弹对比
设计一个尺寸与分段PELE弹一致的普通PELE弹进行模拟，入射速度为1.4 km/s，单层靶板厚度为8 mm。两类弹体的轴向速度(va)变化曲线如图 3所示。可见：穿过第1层靶时，两种弹的轴向速度变化基本一致；180 μs左右两类弹体穿过第2层靶，轴向速度产生分离；此后，普通PELE弹的轴向速度降幅明显大于分段PELE弹。两类弹穿过4层靶板后，弹壳未破碎长度Lr如表 2所示。结果显示，分段PELE弹主要在第1层与第2层靶板之间以及第3层与第4层靶板之间形成破片。在侵彻多层靶板过程中，由于分段PELE弹在弹芯之间设置壳体间隔保护，能够在一定程度上限制每次穿靶后弹体的破碎长度，保留其后续侵彻和产生破片的能力，因此分段PELE弹在侵彻多层靶过程中的破片分配更合理。相比之下，普通结构的PELE弹在穿过前两层靶板后，弹体基本全部破碎，难以对后续目标产生有效毁伤。对比表明：分段PELE弹相较普通PELE弹，具有更强的侵彻能力。
表 2 弹壳未破裂长度Table 2. Length of unbroken PELE shell after penetrationProjectile type Lr/mm 1st layer 2nd layer 3rd layer 4th layer Segmented PELE 39.75 22.82 5.12 0 Normal PELE 31.49 8.71 0 0 2.2 靶板厚度对分段PELE弹终点效应的影响
一般来说，PELE弹主要打击较薄的防护壳体及其后方的目标，而靶板厚度对弹丸的终点效应有重要影响[9]。设弹丸的初始侵彻速度为1.4 km/s，改变单层靶板厚度H(4、5、6和8 mm)，进行数值模拟。不同H下，弹丸穿过各层靶板后壳体径向速度峰值vr以及壳体破裂长度L如图 4所示。
图 4 不同靶板厚度条件下壳体的径向速度峰值和破裂长度Figure 4. Peak radial velocity and crushing length of shell under different thicknesses of target plate如图 4(a)所示，随着H的增加，分段PELE弹穿过各层靶板后壳体的径向速度峰值呈现增加趋势；只有当靶板厚度H为8 mm时，弹丸在侵彻第4层靶板过程中完全被侵蚀消耗，无破片产生，速度为零。在H不同的条件下，弹丸穿过多层靶板后壳体径向速度峰值基本逐层递减，但穿过第2层靶时例外，其壳体径向速度峰值最低。由图 4(b)可知，弹丸穿过第1层靶板后第1段弹体几乎完全破碎，表明该阶段产生的破片最多。故而弹体侵彻第2层靶板时，由于破片提前对靶板进行破坏，导致穿靶后弹壳的径向速度峰值相对其他3次较低。
H对弹壳破碎的影响主要体现在后3层靶板的侵彻过程中。增加H会使弹丸在穿过第2层靶板后的破碎长度L增加；穿过第3层靶板后，L多数情况大于穿过第2层靶后，且此时L先随H的增加而增加，当H=8 mm时，因前两次穿靶时L较大，导致第3次穿靶后L反而减少；穿过第4层靶后，L随着H的增加而逐渐变小。由此可见，当H处于5~6 mm区间时，弹丸穿透各层靶后L的分配比较均匀合理。
虽然增加H可使穿靶后壳体的径向速度峰值更大，但也会使弹体在前几次穿靶过程中壳体破碎过多，不利于后续侵彻和毁伤目标。总体而言，H在5~6 mm时，本弹体结构每次穿靶后的壳体径向速度峰值较高，破碎长度均匀，效果最好。
设计一个与分段PELE弹相同尺寸的钨合金杆，以相同的速度侵彻多层靶板。得到两种弹体侵彻不同厚度靶板后的弹孔直径(D)，如图 5所示。从图 5(b)中可以看到，分段PELE弹侵彻不同厚度靶板时，在第2层靶上的弹孔最大。随着H的增加，分段PELE弹在每层靶板上的开孔规律并不一致：对于第1层靶，D随着H的增加略微增大；第2层靶中，H在5~6 mm时D有极小值；第3层靶板中，靶板厚度为6 mm时D达到极大值；而第4层靶板中，D随着H的增加而减小。
图 5 两种弹体侵彻不同厚度靶板的弹孔直径Figure 5. Layer thickness effect on diameter of crater on multi-layer targets图 6中e表示分段PELE弹与钨合金杆侵彻靶板后弹孔直径的相对偏差，e＞0表示分段PELE弹的弹孔直径更大。数据显示：第1层靶板中的e值均较低，H=4 mm时达到最大值，但此时分段PELE弹的D相对钨合金杆仅有8.67%的提升；对于不同厚度的靶板，分段PELE弹侵彻第2层靶板的D值均远大于钨合金杆，即使对于e最低的6 mm厚靶板，D也增加27.68%，而H=8 mm时，D的增幅达到54.1%；第3层靶板的e值在H=6 mm时取极大值；第4层靶板的e值随着H的增加而减小。
图 6 钨合金杆和分段PELE杆在多层靶中造成的开孔大小差异Figure 6. Relative difference in diameters of craters caused by tungsten alloy rod and segmented PELE penetration综合来看，当H为5~6 mm时，分段PELE弹侵彻各层靶板时的弹孔均大于钨合金杆，即就相对开孔能力而言，本研究中的分段PELE弹适应的最佳靶板厚度依然是5~6 mm。
2.3 侵彻速度对分段PELE弹终点效应的影响
根据以上分析结果，选定靶板厚度为5 mm，改变弹丸初始速度v0(0.8、1.1、1.4和1.7 km/s)，进行数值模拟，所得vr和L随v0的变化曲线如图 7所示。
图 7 不同侵彻速度下壳体的径向速度峰值和破裂长度Figure 7. Peak radial velocity (a) and length of broken PELE shell (b) after different velocity penetration图 7(a)显示，弹丸以不同的v0侵彻靶板时，穿过第1层和第4层靶板后，壳体径向速度峰值随着v0的增加而增大；而穿过第2层靶后，壳体径向速度峰值相对较低；穿过第2层和第3层靶板后，壳体径向速度峰值在v0=1.4 km/s时出现极大值，在v0=1.7 km/s时反而降低。从图 7(b)中也可以看到，当v0达到1.7 km/s时，弹丸穿过第2层和第3层靶板后壳体破裂很少，此时弹丸的横向效应不显著。这表明在一定的靶板厚度条件下，分段PELE弹的v0过高反而不利于其在侵彻多层薄靶后产生适量破片。
当H=5 mm、v0在0.8~1.7 km/s范围内时：弹丸穿过第1层靶板后前段弹体完全破裂，L均为30 mm，体现分段PELE弹可在一定程度上控制破碎范围；之后的几层靶板侵彻过程中，随着侵彻速度的变化，L与vr的变化趋势一致，径向速度峰值越高，破碎长度越长。综合考虑后认为弹丸以1.4 km/s的初始速度侵彻最佳。
3. 结论
(1) 相对普通PELE弹，分段PELE弹侵彻多层靶板时，壳体破裂产生的破片分布更合理，侵彻能力更强。
(2) 分段PELE弹在侵彻多层金属薄靶后，各层靶的弹孔直径普遍大于相同金属杆侵彻形成的弹孔直径，特别是第2层靶板，其弹孔直径增大超过26%；分段PELE弹在侵彻不同厚度靶板时，均在第2层靶上的弹孔直径最大。
(3) 随着靶板厚度的增加，弹丸贯穿各层靶板后的壳体径向速度峰值一同增加，弹丸侵彻第2层至第4层靶时壳体随靶板厚度的增加而产生更多破碎。随着弹丸初速度的增加，弹丸在穿过第2层和第3层靶板后壳体径向速度峰值和破碎长度均在初始速度为1.4 km/s时达到极大值。一定范围内改变靶板厚度和弹丸初始速度，弹丸贯穿第2层靶后的壳体径向速度峰值最低，而穿过第1层靶板后的弹体破碎长度基本保持不变。
(4) 此分段PELE弹在单层靶板厚度接近一半弹丸口径、初始速度在1.4 km/s附近时有较好的终点效应。
-
图 2 Taylor数值实验(a)以及理论和实际密度分布(b)
Figure 2. Virtual Taylor impact test scenario (a) and theoretical and practical density distributions (b)
图 3 γ=–0.8的泡沫金属在Taylor撞击过程中的变形
Figure 3. Deformation of metallic foams with γ=–0.8 in the virtual Taylor impact test
图 4 质点速度时程曲线(a)和X=40 mm处的速度时程曲线(b)
Figure 4. Time histories of particle velocity (a) and velocity at X=40 mm (b)
图 6 局部动态应变时程曲线(a)和X=40 mm的动态应变时程曲线(b)
Figure 6. Local dynamic strain-time curves (a) and time histories of dynamic strain at X=40 mm (b)
图 8 局部动态应力时程曲线((a)～(c))和X=40 mm处的动态应力时程曲线(d)
Figure 8. Local dynamic stress-time curves ((a)–(c)) and time histories of dynamic stress atX=40 mm (d)
表 1 局部最大压实应变、压实应变、最大局部压实应力及其位置
Table 1. Maximum local locking strain, locking strain, and local maximum densification stress and its location
γ Maximum local locking strain Locking strain Maximum local densification stress Value/MPa Lagrangian location/mm –0.8 0.88 0.68 47.44 58 –0.4 0.91 0.70 39.72 56 0 0.92 0.72 34.77 58 0.4 0.93 0.74 30.08 35 0.8 0.96 0.76 30.75 37 
下载: 导出CSV
-
[1] WANG S, DING Y, WANG C, et al. Dynamic material parameters of closed-cell foams under high-velocity impact [J]. International Journal of Impact Engineering, 2017, 99: 111–121. doi: 10.1016/j.ijimpeng.2016.09.013 [2] WANG P, WANG X, ZHENG Z, et al. Stress distribution in graded cellular materials under dynamic compression [J]. Latin American Journal of Solids and Structures, 2017, 14(7): 1251–1272. doi: 10.1590/1679-78253428 [3] 王礼立. 应力波基础 [M]. 北京: 国防工业出版社, 2005. [4] DESHPANDE V S, FLECK N A. High strain rate compressive behaviour of aluminium alloy foams [J]. International Journal of Impact Engineering, 2000, 24(3): 277–298. doi: 10.1016/S0734-743X(99)00153-0 [5] KOLSKY H. An investigation of the mechanical properties of materials at very high rates of loading [J]. Proceedings of the Physical Society Section B, 1949, 62(11): 676. doi: 10.1088/0370-1301/62/11/302 [6] FOWLES R. Conservation relations for spherical and cylindrical stress waves [J]. Journal of Applied Physics, 1970, 41(6): 2740–2741. doi: 10.1063/1.1659298 [7] GRADY D E. Experimental analysis of spherical wave propagation [J]. Journal of Geophysical Research, 1973, 78(8): 1299–1307. doi: 10.1029/JB078i008p01299 [8] SEAMAN L. Lagrangian analysis for multiple stress or velocity gages in attenuating waves [J]. Journal of Applied Physics, 1974, 45(10): 4303–4314. doi: 10.1063/1.1663050 [9] GUPTA Y M. High strain-rate shear deformation of a polyurethane elastomer subjected to impact loading [J]. Polymer Engineering & Science, 1984, 24(11): 851–861. [10] FOREST C A, WACKERLE J, DICK J J, et al. Lagrangian analysis of MIV (Magnetic Impulse-Velocity) gauge experiments on PBX 9502 using the mass-displacement moment function: LA-UR-89-2898 [R]. Los Alamos: Los Alamos National Lab, 1989. [11] WANG L, ZHU J, LAI H. A new method combining Lagrangian analysis with Hopkinson pressure bar technique [J]. Strain, 2011, 47(2): 173–182. doi: 10.1111/str.2011.47.issue-2 [12] 朱珏. 混凝土类材料冲击本构特性的SHPB技术及Lagrange反解法的研究 [D]. 合肥: 中国科学技术大学, 2006. [13] WANG L, DING Y, YANG L. Experimental investigation on dynamic constitutive behavior of aluminum foams by new inverse methods from wave propagation measurements [J]. International Journal of Impact Engineering, 2013, 62: 48–59. doi: 10.1016/j.ijimpeng.2013.06.002 [14] DING Y, WANG S, ZHENG Z, et al. Dynamic crushing of cellular materials: a unique dynamic stress–strain state curve [J]. Mechanics of Materials, 2016, 100: 219–231. doi: 10.1016/j.mechmat.2016.07.001 [15] 张建军. 冲击载荷下梯度多孔金属力学行为研究 [D]. 太原: 太原理工大学, 2016. [16] ZHENG Z, WANG C, YU J, et al. Dynamic stress–strain states for metal foams using a 3D cellular model [J]. Journal of the Mechanics and Physics of Solids, 2014, 72: 93–114. doi: 10.1016/j.jmps.2014.07.013 -
下载:
下载:
点击查看大图
计量
- 文章访问数: 7034
- HTML全文浏览量: 3145
- PDF下载量: 40
