High speed compression behaviour of metallic cellular materials under impact loading
-
摘要: 基于微CT扫描影像信息，建立泡沫金属材料二维细观有限元模型，考虑不规则胞孔的不均匀分布，根据实验结果拟合孔壁材料的弹塑性本构参数。研究了泡沫金属在不同加载速度下的压缩变形机理，重点讨论泡沫金属中弹塑性波的传播、惯性效应和从冲击端传递到静止端的应力变化特征。对于相对密度为0.3的泡沫铝，弹性波速约为5 km/s，与孔壁材料的弹性波速相当，塑性波速表现为随着加载速度的增大而增大。在加载速度为50~100 m/s间变形模式从准静态模式转变为动态模式，未发现明显的临界速度，动态锁死应变随着加载速度的增大而增大。由于塑性波发生反射，试件会发生二次压缩过程，相应地，静止端产生二次应力平台。受惯性作用的影响，二次应力平台也随着加载速度的增大而提高。Abstract: A two-dimensional finite element model was created from a tomographic image of the aluminum foams, which represents the cell shape and geometric distribution of real foams.To determine the mechanical properties of cell wall material, the uniaxial stress versus strain curve, predicted numerically for aluminum foam, was fitted to that measured experimentally.We mainly discuss the shock wave propagation, the inertial effect and the strength of the stress on the stationary end of metallic cellular materials under high speed compression.As for aluminum foams with relative density 0.3, the elastic wave speed is calculated to be 5 km/s, whilst the plastic wave speed increases from 83 to 294 m/s, with the compression velocity increasing from 50 to 200 m/s.Within the compression velocity range of 50-100 m/s, the deformation modes change from random mode to progressive mode.However, no distinct critical velocity are observed.The dynamic locking strain increases with the increasing compression velocity.Second compression process occurs in metallic cellular materials when the plastic wave reflects on the stationary end.Accordingly, the second stress plateau appears on the stationary end, which increases with the increasing compression velocity due to inertia effect.
-
在自然灾害救援和非战争行动中, 经常面临对大块度障碍物快速破除的难题。普通民用爆破中通常采取先机械钻孔、后装药起爆的方法进行施工, 但在缺少时间、机械设备、电力和人员等条件下, 一般方法无法短时内快速清除这些大块度障碍, 常常造成救援目标生命及财产的重大损失。考虑到应急保障的高时效性, 采用串联装药技术更有效[1]。张先锋等[2]、王树有等[3]、曾必强等[4]和涂候杰等[5]分别对串联战斗部前级爆轰对后级的影响进行了数值模拟分析和实验研究, 王成等[6]对同口径串联装药前后级成型关系进行了实验研究。传统的破-爆型串联战斗部前级聚能开孔装药能力有限, 在大块度障碍物破除应用中效果不理想。因此, 本文中根据在坚硬目标上开孔的需要, 提出一种新式破-破型串联爆炸成型弹丸(explosively formed projectile, EFP)聚能装药结构作为多级串联战斗部的前级开孔装药, 并对串联EFP装药隔爆结构和装药前后级延时匹配进行优化设计; 利用有限元程序LS-DYNA对不同起爆延时条件下串联EFP装药侵彻进行数值计算, 并进行实验验证。
1. 串联装药结构
1.1 单级EFP装药选择
本文中选择前期优化设计[7]得到的Ø65 mm球缺型变壁厚EFP装药结构方案。装药结构参数为:炸药采用JH-2, 其密度为1 700 kg/m3, 装药长径比为1.0;药型罩采用紫铜材料, 罩顶厚2.1 mm, 罩内表面曲率半径为67 mm, 外表面曲率半径为62 mm。
1.2 串联EFP关键技术
为确保精确控制两级EFP装药的起爆时间, 采用本课题专门设计的精确延时起爆控制器, 设定延时间隔精度为0.1 μs, 可通过专用的应用软件在0~200 μs区间任意设定延时, 图 1是2路起爆控制信号输出时示波器采集到的触发信号波形。
1.3 隔爆体设计原则
本研究串联EFP采用逆序起爆方式, 前级装药爆炸后产生的爆轰产物和空气冲击波不可避免地作用在后级装药上, 引起后级EFP形状和性能的改变, 从而改变后级EFP的侵彻性能。因此要设计一种简单的隔爆结构, 降低前后级之间的影响。本研究的隔爆结构中选用的聚氨酯泡沫材料是一种密度小、隔爆性能好的抗冲击波材料, 广泛应用于国防、军事领域[8]。
2. 数值模拟分析
2.1 串联EFP数值模型
为了分析隔爆体及起爆延时对串联EFP装药成型性能的影响, 利用有限元软件LS-DYNA3D对其成型过程进行数值模拟, 串联装药模型如图 2所示。模型涉及聚能装药成型和侵彻2个部分, 采用流-固耦合算法来分析此类情况更贴近实际。模型中, 炸药、药型罩、隔爆体和空气等介质均采用Euler算法, 靶板则采用Lagrange算法。串联EFP装药为8701炸药, 采用高能炸药材料模型和JWL状态方程描述; 药型罩材料为军用紫铜, 用Grüneisen状态方程和Johnson-Cook本构模型描述[9-10]; 隔爆体为聚氨酯泡沫, 用Grüneisen状态方程和Elastic_Plastic_Hydro本构模型描述, 各材料状态方程参数[11-12]如表 1~3所示, 其中:ρ为密度, v为爆速, pCJ为炸药的C-J爆压, E0为材料的初始内能, A、B、R1、R2和ω为炸药的材料常数, C为材料的us-up截距, S1、S2和S3为斜度系数, γ0为材料的Grüneisen系数, a为γ0的一阶体积修正。
表 1 JH-2炸药计算参数Table 1. Computational parameters for JH-2ρ/(g·cm-3) v/(m·s-1) pCJ/GPa A/GPa E0/(J·m-3) B R1 R2 ω 1.70 8 400 30 56.4 10.0 6.801 4.1 1.3 0.36 表 2 紫铜和聚氨酯泡沫计算参数Table 2. Computational parameters for copper linerρ/(g·cm-3) C/(m·s-1) S1 S2 S3 γ0 a E0 8.96 4 750 3.8 2.74 0.125 1.346 0.34 0.0 表 3 聚氨酯泡沫计算参数Table 3. Computational parameters for polyurethane foamρ/(g·cm-3) C/(m·s-1) S1 S2 S3 γ0 a E0 0.05 886 0.78 0.0 0.0 1.55 0.00 0.0 2.2 不同炸高条件下单级EFP侵彻过程分析
要充分发挥两级串联EFP装药的侵彻能力, 首先需要确定装药的有利炸高。如果炸高太小, EFP弹丸还未完全成型, 其速度和长径比仍在快速变化中, 侵彻深度和穿孔直径等侵彻效果随炸高变化而发生较大改变。但炸高达到一定程度后, EFP成型性能已经基本稳定时, 侵彻效果变化逐渐减小。文中选取了6种不同炸高进行了数值模拟和实验, 结果如图 3所示, 其中:H为炸高, D为装药直径。
对比图 3中6种炸高条件下钢靶剖面图得出:在炸高与装药直径比小于2.8时, EFP侵彻钢靶孔径较大, 但EFP未得到有效拉伸, 侵彻深度不大; 随着炸高的增加, 侵彻深度逐渐增加, 同时孔径逐渐减小; 在炸高与装药直径比大于3.2时, 炸高增加, 侵彻深度与孔径基本不变。表 4为不同炸高下EPF侵彻钢靶的模拟和实验结果, 其中h为侵彻深度, d为侵彻孔径。分析表 4中的数据发现, 数值模拟与实验结果在侵彻深度和孔径大小上吻合较好, 误差在5%以内, 说明数值模拟结果较贴近实际, 为后面串联装药模拟提供了一种可行的方法。
表 4 EFP侵彻钢靶数值模拟与实验数据Table 4. Simulational and experimental results of EFP penetrating steel targetsH/mm h/D d/D 模拟值 实验值 模拟值 实验值 150 0.73 0.74 0.557 0.563 180 0.80 0.78 0.525 0.530 210 0.85 0.86 0.498 0.495 240 0.89 0.88 0.477 0.472 270 0.90 0.89 0.460 0.454 300 0.90 0.90 0.455 0.451 2.3 不同隔爆体形状串联EFP成型过程数值模拟
隔爆效果不仅与材料性能有关, 隔爆体的形状对后级装药成型影的响也十分巨大。文中选取了3种不同结构的隔爆体进行分析。图 4为圆柱形隔爆体、锥口向上隔爆体和锥口向下3种条件下串联EFP装药200 μs时刻后级EFP成型状态。通过对比图 4中的3种后级EFP形态, 发现隔爆体结构的不同对减少前级爆轰对后级成型的影响差异很大。采用圆柱形隔爆体时后级EFP基本成型, 但头尾连续性差, EFP侵彻能力大大降低; 采用圆柱底部挖出一个锥形空腔且锥口向上的隔爆体, 隔爆效果差, 后级弹丸基本无翻转, 后级装药侵彻能力基本消失; 锥口向下的隔爆体, 隔爆效果较理想, 可作为串联EFP的隔爆结构基本形状, 后级EFP成型较好, 但在前级爆轰场影响下, EFP长径比增大, 需要进一步优化前后级起爆时间, 来提高后级EFP的成型效果。
图 4 不同隔爆体200 μs时刻后级EFP形态Figure 4. Shapes of the postpositive EFP charge at 200 μs with different explosion-proof bodies2.4 不同延时串联EFP侵彻过程数值模拟分析
串联EFP前后级延时起爆时间Δt对后级EFP的成型和侵彻性能影响很大[8]。在两级装药间距一定的情况下, 合理控制延时起爆时间Δt, 使后级EFP受前级装药爆轰场影响最小, 才能最大限度地保持后级EFP装药的成型性能。分别对5种不同延时条件下串联EFP侵彻靶板进行数值模拟, 结果见表 5。
表 5 串联EFP装药侵彻数值模拟结果Table 5. Simulation results of tandem EFP penetrationΔt/μs h/mm h/D d/mm d/D 0 92 1.415 22.2 0.341 10 103 1.584 21.6 0.332 20 115 1.769 21.4 0.329 25 104 1.600 21.7 0.327 30 98 1.507 21.3 0.334 图 5给出了侵彻深度和侵彻孔径随起爆延时的变化关系。分析图 5中孔径和孔深的变化规律可以发现, 随着前后级延时起爆时间的增加, 后级EFP的侵彻深度先增大后快速减小, 说明前级爆轰场到达后级装药时, 会严重影响后级EFP的成型。所以, 该串联结构较为合理的延时为20 μs。图 6为同时起爆和延时20 μs的数值模拟侵彻结果。
图 5 开孔深度和直径随起爆延时的变化Figure 5. Variation of penetration depth and aperture with delay time3. 串联EFP隔爆结构实验研究
3.1 不同隔爆形体状串联EFP侵彻实验
根据前期研究[1]得到的隔爆体形状对串联EFP装药侵彻性能的影响规律, 进行6发串联EFP侵彻的验证性实验。实验分3组, 分别采用圆柱形、锥口向上和锥口向下3种隔爆体, 每组进行2发实验。图 7为3种不同条件下EFP的侵彻性能实验结果。
图 7 不同隔爆体串联EFP侵彻结果Figure 7. Penetration results of tandem EFPs with different explosion-proof bodies表 6给出了串联装药对靶板的侵彻结果, 通过对比可以发现, 同时起爆条件下, 锥口向下的隔爆体侵彻深度明显优于其他2种结构。但后级侵彻的开孔孔径较小, 难以满足后级爆破子弹随进要求, 需进一步优化。
表 6 串联EFP装药侵彻钢靶实验结果Table 6. Experimental results of tandem EFPs penetrating steel targets隔爆结构 h/mm d/mm 圆柱形 82 15.2 锥口向上 60 16.6 锥口向下 98 16.3 3.2 不同延时条件下串联EFP侵彻钢靶实验
图 8为实验设置图。前级装药炸高取210 mm, 装药间距150 mm条件下, 采用中心起爆方式, 使用8#电雷管同时起爆两级装药, 分别对延时0、10、20、25和30 μs等5种情况进行侵彻靶板分别进行2发实验, 结果依次记录为A1~A10。为对比优化后的侵彻性能, 同时进行2组分2次单独侵彻实验作为参照, 结果记录为B1和B2。侵彻结果显示串联EFP开孔形状前后基本一致。后级EFP受到前级装药爆轰场的影响, 速度和长径比都不可避免有所下降, 开孔直径比前级EFP小。
表 7给出了不同延时条件下的侵彻深度和最小侵彻孔径的实验结果。由表 7可以看出:随着延时的增加, 后级EFP受到前级爆轰场影响, 形状变得更加细长, 侵彻的孔径逐渐变小; 而穿孔深度先逐渐增加, 随后逐渐减小。这说明起爆时间间隔太大, 前级爆轰场到达后级装药后, 会严重影响后级EFP的成型, 大大降低了后级EFP的速度。将表 7与表 5进行比较, 可以看出两者最大仅相差3.7%, 可见实验结果与数值模拟结果吻合较好。
表 7 串联EFP装药侵彻钢靶结果Table 7. Experimental results of tandem EFPs penetrating steel targets编号 Δt/μs h/mm d/mm A1 0 94 21.4 A2 0 95 21.3 A3 10 105 21.0 A4 10 102 21.2 A5 20 117 20.6 A6 20 115 20.8 A7 25 107 20.4 A8 25 108 20.3 A9 30 100 20.3 A10 30 102 20.2 B1 121 21.6 B2 119 21.8 图 9所示靶板, 依次为分2次侵彻、延时0和延时20 μs侵彻实验结果。可以看出:同时起爆(Δt=0)时, 后级EFP长径比和速度受影响较大, 侵彻孔径较分2次侵彻时减小0.4 mm, 后级EFP侵彻深度只有分2次侵彻时第2次侵彻深度的58.9%, 侵彻效果大大降低。延时起爆20 μs时, 侵彻孔径较分2次侵彻时减小1 mm, 两级整体侵彻深度为分2次侵彻时的96.7%, 这体现了延时起爆对后级EFP侵彻性能的重要性。
图 9 不同起爆延时下串联EFP装药侵彻钢靶实验结果Figure 9. Tandem EFP charge penetrating steel targets with different delay time4. 结论
(1) 在已有的研究成果基础上, 分析了隔爆体形状和前后级装药延时对串联EFP侵彻能力的影响, 得到了较合理的串联装药结构。
(2) 通过数值模拟和实验研究, 对串联EFP装药隔爆结构形状和延时匹配进行了比较分析, 延时起爆20 μs时, 串联侵彻深度为分2次侵彻的96.7%, 比同时起爆侵彻深度提高了约22.8%, 大大提高了后级装药的利用效率。这可为下一步多级串联装药研究提供依据。
-
图 1 泡沫金属材料二维细观有限元模型
Figure 1. Two-dimensional mesoscale finite element models of metallic cellular materials
图 2 泡沫金属单轴压缩的名义应力应变曲线
Figure 2. Numerically predicted and experimentally measured uniaxial compressive stress versus strain curves for closed-cell aluminum foams
图 3 泡沫金属在不同速度压缩下的变形图
Figure 3. Simulated deformation modes of aluminum foam (relative density of 0.3)under different impact velocities
图 4 泡沫金属动态压缩的名义应力应变曲线
Figure 4. Numerically predicted uniaxial compressive stress versus strain curves for aluminum foams
图 5 冲击波传播示意图
Figure 5. Schematic diagram of shock wave propagation in metallic cellular materials
表 1 泡沫金属高速冲击加载计算结果
Table 1. Summary of FE predictions for aluminum foam (relative density of 0.3) under impact loading
vi/(m·s-1) vp/(m·s-1) εlock ρlock/(kg·m-3) σy(εlock)/MPa 50 83 0.60 2 025 34.5 100 167 0.60 2 025 34.5 150 234 0.64 2 314 40.1 200 294 0.68 2 700 52.7 
下载: 导出CSV
-
[1] Lopatnikov S L, Gama B A, Haque M J, et al. Dynamics of metal foam deformation during Taylor cylinder-Hopkinson bar impact experiment[J]. Composite Structures, 2003, 61(1/2): 61-71. [2] Tan P, Harrigan J, Reid S. Inertia effects in uniaxial dynamic compression of a closed cell aluminium alloy foam[J]. Materials Science and Technology, 2002, 8: 480-488. [3] Tan P J, Reid S R, Harrigan J J, et al. Dynamic compressive strength properties of aluminium foams. Part I-Experimental data and observations[J]. Journal of the Mechanics and Physics of Solids, 2005, 53(10): 2174-2205. doi: 10.1016/j.jmps.2005.05.007 [4] Elnasri I, Pattofatto S, Zhao H, et al. Shock enhancement of cellular structures under impact loading: PartⅠ-Experiments[J]. Journal of the Mechanics and Physics of Solids, 2007, 55(12): 2652-2671. doi: 10.1016/j.jmps.2007.04.005 [5] Merrett R P, Langdon G S, Theobald M D. The blast and impact loading of aluminium foam[J]. Materials & Design, 2013, 44: 311-319. [6] Pattofatto S, Elnasri I, Zhao H, et al. Shock enhancement of cellular structures under impact loading: PartⅡ-Analysis[J]. Journal of the Mechanics and Physics of Solids, 2007, 55(12): 2672-2686. doi: 10.1016/j.jmps.2007.04.004 [7] Karagiozova D, Langdon G S, Nurick G N. Propagation of compaction waves in metal foams exhibiting strain hardening[J]. International Journal of Solids and Structures, 2012, 49(19/20): 2763-2777. [8] Liu Y D, Yu J L, Zheng Z J, et al. A numerical study on the rate sensitivity of cellular metals[J]. International Journal of Solids and Structures, 2009, 46: 3988-3998. doi: 10.1016/j.ijsolstr.2009.07.024 [9] Ma G W, Ye Z Q, Shao Z S. Modeling loading rate effect on crushing stress of metallic cellular materials[J]. International Journal of Impact Engineering, 2009, 36: 775-782. [10] Reid S R, Peng C. Dynamic uniaxial crushing of wood[J]. International Journal of Impact Engineering, 1997, 19(5/6): 531-570. [11] Lopatnikov S L, Gama B A, Gillespie J W. Modeling the progressive collapse behavior of metal foams[J]. International Journal of Impact Engineering, 2007, 34(3): 587-595. doi: 10.1016/j.ijimpeng.2005.12.004 [12] Lopatnikov S L, Gama B A, Haque M J, et al. High-velocity plate impact of metal foams[J]. International Journal of Impact Engineering, 2004, 30(4): 421-445. doi: 10.1016/S0734-743X(03)00066-6 [13] Harrigan J J, Reid S R, Tan P J, et al. High rate crushing of wood along the grain[J]. International Journal of Mechanical Sciences, 2005, 47(4/5): 521-544. [14] 张健, 赵桂平, 卢天健.闭孔泡沫铝应变率效应的试验和有限元分析[J].西安交通大学学报, 2010, 44(5): 97-101.Zhang Jian, Zhao Gui-ping, Lu Tian-jian. Experimental and numerical study on strain rate effects of close-celled aluminum foams[J]. Journal of Xi'an Jiaotong University, 2010, 44(5): 97-101. [15] Hallquist J O. LSTC LS-DYNA user's manual[Z]. Livermore, CA, US: Livermore Software Technology Corporation, 2007. [16] Wang Li-li. Foundation of stress waves[M]. Beijing: National Defense Industry Press, 2005. 期刊类型引用(1)
1. Shuai Qin，Hao Liu，Jianhui Wang，Qiang Zhao，Lei Zhang. Physics-data coupling-driven method to predict the penetration depth into concrete targets. Theoretical & Applied Mechanics Letters. 2024(03): 184-192 . 
必应学术其他类型引用(0)
-
下载:
必应学术
计量
- 文章访问数: 2950
- HTML全文浏览量: 341
- PDF下载量: 576
- 被引次数: 1