Numerical simulation on penetration of concrete target by shaped charge jet with SPH method
-
摘要: 在完全变光滑长度SPH(smoothed particle hydrodynamics)方法的基础上，利用F.Ott等提出的修正SPH方法处理在求解多介质大密度问题时的数值不稳定性问题，运用Holmquist-Johnson-Cook本构模型处理混凝土在冲击载荷下的变形和损伤问题，对聚能装药射流侵彻混凝土靶板的过程进行了数值模拟，同时利用LS-DYNA非线性有限元程序进行对比，分析了2种方法得到的混凝土von Mises应力变化、射流头部特定节点处的速度变化及裂纹演变，验证了SPH方法的准确性。分析了另外2种不同尺寸的靶板在射流侵彻作用下的破坏形式，结果符合射流侵彻物理规律，表明该方法适合模拟聚爆炸与冲击等大变形破坏等问题。Abstract: On the basis of the smoothed particle hydrodynamics (SPH) method with fully variable smoothing lengths, the modified Ott-Schnetter SPH method was used to cope with the computational instability in the large density gradient problems. The Holmquist-Johnson-Cook constitutive model was used to cope with the deformation and damage of concrete under impact. By comparing with the corresponding ones simulated by the non-linear finite element program LS-DYNA, the changing courses of the von-Mises stress and cracks in the concrete targets, and the velocities at the some special points of the jet head were analyzed, which proved the feasibility and accuracy of the SPH method. Other two concrete plates with different sizes penetrated by shaped charge jets were also simulated. The results are in good agreement with the physical principle. So this method can be used to deal with the multi-material and large deformation problems such as detonation and impact.
-
Key words:
- mechanics of explosion /
- penetration /
- smoothed particle hydrodynamics /
- concrete /
- damage
-
电热化学发射技术使用等离子体发生器代替常规点火源，将电能转化为等离子体内能来点燃固体发射药。等离子体与固体发射药相互作用不仅能显著缩短发射药点火延迟时间[1-3]，还能消除发射药的温度效应[4]，提高燃速，从而有效提高炮口初速和动能[5]。常规发射通常使用指数燃速公式来描述发射药燃烧特性[6-7]，C.R.Woodley等[8-9]假定电热化学发射时注入的电能只用于增强发射药的内能，引入了电功率增强因子，给出了与输入电功率相关的发射药燃速公式。M.J.Taylor等[10]进一步研究了辐射、发射药侵蚀和电热裂解对燃速的影响，考虑了引起燃速增强的电功率阈值的影响。李海元等[11-12]、A.Brik等[13]利用改进的密闭爆发器对等离子体增强发射药燃速的特性进行了实验研究; 实验表明等离子体作用下，发射药燃气压力变化迅速，燃气压力梯度增大，点火和燃烧特性发生显著变化。因此，等离子体作用下燃气压力梯度对发射药燃速的影响不可忽略。
本文中，拟通过密闭爆发器实验得到燃气压力曲线，分析等离子体作用下固体发射药的燃烧特性。基于指数燃速公式，结合燃气压力梯度和电功率增强对发射药燃速的影响，得到等离子体作用下固体发射药的瞬态燃速公式。
1. 基本原理
1.1 实验装置
图 1所示为等离子体点火密闭爆发器实验装置结构示意图，主要由脉冲成形网络、底喷式等离子体发生器和密闭爆发器等组成。密闭爆发器为耐高压厚壁圆筒，其中一端装有底喷式等离子体发生器，其主体由聚乙烯毛细管组成，一端采用杆状封闭电极，另一端装有敞开式环状电极，通过爆炸丝将两端电极导通。脉冲成形网络由多个模块并联组成，能通过时序设置产生幅值、脉宽可调的脉冲电流。固体发射药试样由硝化棉纸包裹放置于密闭爆发器内。
图 1 等离子体点火密闭爆发器实验装置结构示意图Figure 1. Experimental setup of closed bomb with plasma igniter实验测试系统主要由传感器和数据采集设备组成。燃气压力采用Kistler6215压力传感器测量，发生器两端的电压和电流分别采用电阻分压器和Rogowski线圈测量。采用JV5200瞬态记录仪对脉冲电源和发生器的电参数以及密闭爆发器内的压力信号进行同步采集记录。
1.2 瞬态燃速理论模型
1.2.1 基本假设
忽略密闭爆发器内固体发射药试样燃烧过程中气相压力和温度随空间的变化，即假设任意时刻密闭爆发器内各位置的气相压力和温度均相同。
固体发射药试样燃烧过程满足几何燃烧定律假设：(1)装药的所有颗粒具有均一的理化性质以及完全相同的几何形状和尺寸；(2)所有药粒表面同时着火；(3)所有药粒具有相同的燃烧环境。
由于气相与密闭爆发器壁面存在换热，需考虑试样燃烧过程中的热损失。本文中采用热损失因数修正发射药火药力的方式来计算发射药试样燃烧过程中的热损失。
1.2.2 气相状态方程
考虑点火药的影响，得到密闭爆发器内的气体状态方程[14]：
p{V0−mρp[1−ψ(t)]−αmψ(t)−αigmig(t)}=T[Rmψ(t)+Rigmig(t)] (1) 式中:p为燃气压力，Pa；V0为密闭爆发器容积，m3；m为装药质量，kg；ρp为固体发射药密度，kg/m3；ψ (t)为t时刻发射药已燃体积分数；α为余容，m3/kg；αig为点火药余容，m3/kg；mig(t)为t时刻已燃点火药质量，kg；T为混合气体温度，K；R为发射药燃气气体常数，J/(mol·K)；Rig为点火药燃气气体常数，J/(mol·K)。
1.2.3 气相能量守恒方程
气相能量守恒方程[14]为：
(1−c1)fmψ(t)+figmig(t)+(κ−1)cplEpl(t)=T[Rmψ(t)+Rigmig(t)] (2) 式中:cl为热损失因数；f为发射药火药力，J/kg；fig为点火药火药力，J/kg；κ为燃气比热比；cpl为密闭爆发器中输入等离子体电能的利用率；Epl(t)为t时刻输入等离子体发生器的电能，J。
式(2)通用于常规点火和等离子体点火。常规点火时，Epl(t)=0；等离子体点火时，mig(t)=0。
1.2.4 固体发射药瞬态燃速公式
根据Vieille定律，常规点火时固体发射药燃速公式通常采用指数形式[14]：
u=u1pn1 (3) 式中:u1为燃速系数，n1为燃速指数。
C.R.Woodley等[8]在Vieille定律的基础上引入电功率增强因子，得到等离子体增强发射药燃速公式:
u=u1pn1(1+βePe) (4) 式中:βe为电功率增强因子，MW-1; Pe为输入发生器的电功率，MW。
等离子体点火时，在Woodley燃速公式的基础上，加入燃气压力梯度对发射药燃速的影响[15]，得到等离子体作用下固体发射药瞬态燃速公式：
u=u1pn1[1+A(t)n1u21p2n1+1dpdt](1+βePe) (5) 式中:A(t)为发射药燃烧过程中与压力及火焰结构有关的时间函数。
1.2.5 相关系数处理方法
由气体状态方程和能量守恒方程可知：
p{V0−mρp[1−ψ(t)]−αmψ(t)−αigmig(t)}=(1−c1)fmψ(t)+figmig(t)+(κ−1)cplEpl(t) (6) 首先，通过常规点火实验确定热损失因数。常规点火时, Epl(t)=0。发射药试样完全燃完时刻有：ψ(t) =1, p=pm, pm为实验测得的最高燃气压力。由此可求出发射药燃烧过程中的热损失因数：
c1=1+figmig(t)fm−pmf[V0m−α−mig(t)mαig] (7) 其次，通过等离子体点火实验确定电能利用率。等离子体点火时，mig(t)=0。同样发射药试样完全燃完时刻有:ψ(t) =1, p=pm。假定等离子体点火时热损失因数未发生变化，进一步求出等离子体点火过程中的电能利用率：
cpl=pm(V0−αm)−(1−c1)fm(κ−1)Epl(t) (8) 1.2.6 拟合精度
采用均方误差σ来衡量模拟压力曲线与实验压力曲线间的误差：
σ=√n∑i=1[ps(i)−pt(i)]2n (9) 式中:n为压力曲线上选取不同时刻的点数；ps(i)、pt(i)分别为同一时刻模拟和实验测得的压力。
2. 实验结果与分析
对4/7高固体发射药试样进行了2发常规点火和3发等离子体点火密闭爆发器实验。其中密闭爆发器容积为145 cm3，装药量为36.1 g; 实验中采用单个模块放电，电容约为1 300 μF，电感为40 μH。实验参数与结果如表 1所示，其中第3、4发实验发射药置于密闭爆发器中间，第5发实验发射药置于发生器出口处。表 1中, Uc为电容器放电电压；tig为发射药点火延迟时间(定义为爆发器内压力达到20 MPa的时刻)；tend为发射药完全燃完时刻，即爆发器内压力达到峰值pm对应的时刻。
表 1 实验参数和结果Table 1. Experimental parameters and results编号 点火方式 Uc/kV tig/ms tend/ms pm/MPa 实验1 2号电点火 1.604 5.05 298 实验2 2号电点火 1.673 4.95 300 实验3 等离子体点火 8.3 0.261 2.93 318 实验4 等离子体点火 10.0 0.197 2.44 319 实验5 等离子体点火 10.1 0.167 2.16 329 实验测得的膛内压力曲线如图 2所示，由图可知等离子体点火大大缩短了发射药试样的点火延迟时间和燃烧过程时间；同时，随着等离子体的注入，密闭爆发器内燃气压力不断增大。
实验测得的膛内压力梯度曲线如图 3所示，由图 3可知在发射药燃烧初始时刻，等离子体作用下的燃气压力梯度远大于常规点火，而且等离子体明显影响了燃气的压力梯度，随着发射药的燃烧，影响程度逐渐降低。
图 4所示为常规点火时实验和数值模拟得到的压力曲线，其中数值模拟压力曲线分别由指数燃速公式和瞬态燃速公式推得。对比可知，指数燃速公式和瞬态燃速公式均有较高的精确度；但是在发射药燃烧后期瞬态燃速公式与实验结果更吻合。
图 4 常规点火时实验与数值模拟压力曲线Figure 4. Experimental and simulated pressure-time curvesafter conventional ignition通过分析常规点火和等离子体点火实验相关数据，得到4/7高发射药的电功率增强因子为0.005 MW-1。图 5~7分别给出了第3、4和5发等离子体点火实验时压力与放电功率随时间变化的曲线，其中压力曲线分别通过实验测量、Woodley燃速公式和瞬态燃速公式所得。当加入电功率增强因子后，Woodley燃速公式和瞬态燃速公式得到的压力曲线均符合实验曲线。因此，拟合得到的电功率增强因子适用于4/7高发射药试样。
图 5 第3发实验的压力与输入电功率曲线Figure 5. Pressure-time and electric power-time curves for test 3
图 6 第4发实验的压力与输入电功率曲线Figure 6. Pressure-time and electric power-time curves for test 4
图 7 第5发实验的压力与输入电功率曲线Figure 7. Pressure-time and electric power-time curves for test 5结合图 3中的压力梯度曲线，分析等离子体输入过程中压力梯度与发射药位置和输入电功率的关系。对比第3、4发实验曲线可知，发射药处于相同位置时，早期燃气压力梯度峰值与输入电功率峰值成正比，且压力梯度变化趋势与输入电功率曲线相似。对比第4、5发实验曲线可知，在相似电功率下，缩短发生器与发射药距离提高了电能利用系数，从而增大了燃气压力梯度。因此，缩短等离子体与发射药距离，提高输入电功率均能获得更理想的等离子体增强效应。
对比图 5~7中各压力曲线可知，在燃烧中间段Woodley燃速公式所得压力曲线略低于实验压力曲线。其原因可能是Woodley燃速公式中仅考虑了等离子体注入期间的电增强效应，未考虑放电结束后的增强效应。因此，脉冲放电结束后Woodley燃速公式所得压力曲线与实验有偏差，而在瞬态燃速公式中通过压力梯度项考虑了放电结束后的燃速增强效应，提高了数值模拟精度。
进一步对比各燃速公式的拟合程度，计算得拟合的压力曲线的均方误差如表 2所示。由表 2可知，随着输入电功率的增加以及发生器与发射药间距的减小，Woodley燃速公式均方误差增大；瞬态燃速公式得到的压力均方误差小于Woodley燃速公式，且受到输入电能、电功率和发生器与发射药间距的影响更小。因此，瞬态燃速公式得到的膛压曲线与实验数据符合情况优于Woodley燃速公式，能更真实地反映等离子体作用下密闭爆发器内燃气压力的变化。
表 2 模拟压力曲线与实验压力曲线间的均方误差Table 2. Mean squared errors between simulated pressure curves and test ones编号 σ/MPa Woodley燃速公式 瞬态燃速公式 实验3 4.325 4.294 实验4 9.312 4.910 实验5 13.506 5.715 3. 结论
利用密闭爆发器实验，研究了等离子体增强4/7高固体发射药的燃速特性，实验表明等离子体作用下燃气压力梯度增大，燃速明显增强。根据实验数据，综合考虑燃气压力梯度和燃速增强因子的影响，拟合了4/7高固体发射药等离子体作用下的瞬态燃速公式，计算得到等离子体注入期间的燃速增强因子为0.005 MW-1。瞬态燃速公式与实验压力曲线符合程度比C.R.Woodley等[8]提出的燃速公式更高，能更精确描述固体发射药在等离子体作用下的燃烧过程。实验表明，在低装填密度时，缩短等离子体与发射药的距离和提高电功率均能获得更明显的等离子体增强效应。
-
图 6 聚能射流侵彻过程中不同时刻，混凝土靶板中Von Mises应力的分布
Figure 6. Von Mises stress distribution in concrete target penetrated by shaped charge jet at different times
图 9 在聚能射流的侵彻下不同时刻混凝土靶板的损伤情况
Figure 9. Damage in concrete target penetrated by shaped charge jet at different times
图 10 不同方法模拟得到的裂纹长度及分布情况
Figure 10. Length and distribution of crack simulated by different methods
图 11 另外2种不同尺寸的混凝土靶板在聚能射流侵彻作用下的损伤
Figure 11. Damage in other two concrete targets with different sizes penetrated by shaped charge jets
-
[1] 段建, 杨黔龙, 周刚, 等.串联随进战斗部侵彻混凝土靶实验研究[J].爆炸与冲击, 2007, 27(4):364-369. doi: 10.3321/j.issn:1001-1455.2007.04.012Duan Jian, Yang Qianlong, Zhou Gang, et al. Experimental studies of a tandem follow-through warhead penetrating concrete target[J]. Explosion and Shock Waves, 2007, 27(4):364-369. doi: 10.3321/j.issn:1001-1455.2007.04.012 [2] Swegle J W, Attaway S W. On the feasibility of using smoothed particle hydrodynamics for underwater explosion calculations[J]. Computational Mechanics, 1995, 17(3):151-168. doi: 10.1007-BF00364078/ [3] Libersky L D, Petscheck A G. Smoothed particle hydrodynamics with strength of materials[C]//Trease H, Fritts J, Crowley W. Proceedings of the Next Free Lagrange Conference. NY: SpringerVerlag, 1991, 395: 248-257. [4] Libersky L D, Petscheck A G, Carney T C, et al. High strain Lagrangian hydrodynamics: A three-dimensional SPH code for dynamic material response[J]. Journal of Computational Physics, 1993, 109(1):67-75. http://cn.bing.com/academic/profile?id=f5b87b9c4f2fcce90394ac322bb29f60&encoded=0&v=paper_preview&mkt=zh-cn [5] Liu M B, Liu G R, Zong Z, et al. Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodology[J]. Computers & Fluids, 2003, 32(3):305-322. http://cn.bing.com/academic/profile?id=78358a2b72043bd2db9a95d83b7df716&encoded=0&v=paper_preview&mkt=zh-cn [6] Liu M B, Liu G R, Lam K Y, et al. Meshfree particle simulation of the detonation process for high explosives in shaped charge unlined cavity configurations[J]. Shock Waves, 2003, 12(6):509-520. doi: 10.1007/s00193-003-0185-2 [7] Qiang Hongfu, Wang Kunpeng, Gao Weiran. Numerical simulation of shaped charge jet using multi-phase SPH method[J]. Transactions of Tianjin University, 2008, 14(1):495-499. http://cn.bing.com/academic/profile?id=2700e225ea1da83de80808d57d0ac670&encoded=0&v=paper_preview&mkt=zh-cn [8] Ott F, Schnetter E. A modified SPH approach for fluids with large density differences[J]. Arxiv Physics E-prints, 2003:3112. doi: 10.1002-zaac.201100286/ [9] 强洪夫, 高巍然.完全变光滑长度SPH法及其实现[J].计算物理, 2008, 25(5):569-575. doi: 10.3969/j.issn.1001-246X.2008.05.008Qiang Hongfu, Gao Weiran. SPH method with fully variable smoothing lengths and implementation[J]. Chinese Journal of Computational Physics, 2008, 25(5):569-575. doi: 10.3969/j.issn.1001-246X.2008.05.008 [10] 强洪夫, 王坤鹏, 高巍然.基于完全变光滑长度SPH方法的HE爆轰过程的数值试验[J].含能材料, 2009, 17(1):27-31. doi: 10.3969/j.issn.1006-9941.2009.01.008Qiang Hongfu, Wang Kunpeng, Gao Weiran. Numerical study of high explosive detonation process using SPH method with fully variable smoothing lengths[J]. Chinese Journal of Energetic Materials, 2009, 17(1):27-31. doi: 10.3969/j.issn.1006-9941.2009.01.008 [11] Johnson G R, Cook W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures[C]//Proceedings of the Seventh International Symposium on Ballistics. Hague, Netherlands, 1983: 571-574. [12] Steinberg D J, Cochran S G, Guinan M W. A constitutive model for metals applicable at high strain rate[J]. Journal of Applied Physics, 1980, 51(3):1498-1504. doi: 10.1063/1.327799 [13] Holmquist T J, Johnson G R, Cook W H. A computational constitutive model for concrete subjected to large strains, high strain rates, and high pressures[C]//Proceedings of the 14th International Symposium on Ballistics. Quebec, Canada, 1993: 591-600. [14] Monaghan J J. Smoothed particle hydrodynamics[J]. Reports on Progress in Physics, 2005, 68(8):1703-1759. doi: 10.1088/0034-4885/68/8/R01 [15] Liu G R, Liu M B.光滑粒子流体动力学: 一种无网格粒子法[M].韩旭, 杨刚, 强洪夫, 译.长沙: 湖南大学出版社, 2005: 195-197. [16] Livermore Software Technology Corporation. LS-DYNA keyword user's manual[M]. Livermore: Livermore Software Technology Corporation, 2012:17-45. [17] Liu M B, Liu G R, Zong Z, et al. Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodology[J]. Computers and Fluids, 2003, 32(3):305-322. http://cn.bing.com/academic/profile?id=78358a2b72043bd2db9a95d83b7df716&encoded=0&v=paper_preview&mkt=zh-cn 期刊类型引用(15)
1. 薛颖杰，陈智刚，杨芮，崔晋，康彦姝，付建平. 椭圆截面弹体侵彻多层混凝土靶的数值仿真. 中北大学学报(自然科学版). 2025(01): 35-44 . 
百度学术2. 张丁山，张博，付良，徐笑，李鹏飞. 弹体尾部斜锥面形状对侵彻偏转的影响. 高压物理学报. 2024(01): 145-154 . 百度学术3. 刘泓甫，黄风雷，白志玲，段卓平. 刚性弹体带攻角斜侵彻贯穿混凝土靶板的理论模型. 兵工学报. 2023(08): 2381-2390 . 百度学术4. 朱超，张晓伟，张庆明，张陶. 弹体斜侵彻双层钢板的结构响应和失效研究. 爆炸与冲击. 2023(09): 140-154 . 
本站查看5. 杨璞，李继承，陈建良，张斌，何丽灵，陈刚. 撞击姿态对构型弹体非正侵彻多层间隔钢靶弹道特性的影响规律. 爆炸与冲击. 2023(09): 123-139 . 本站查看6. 张丁山，全嘉林，付良，张博，徐笑. 侵彻弹体尖卵形头部形状对偏转力矩的影响. 火炸药学报. 2023(09): 834-839 . 百度学术7. 姜安邦，李典，李永清，侯海量. 偏转式抗侵彻防护技术研究现状. 材料导报. 2023(24): 99-107 . 百度学术8. 陈柏翰，沈子楷，邹慧辉，王伟光，王可慧. 硬目标内弹体斜侵彻状态的基本演化特性. 兵工学报. 2023(S1): 59-66 . 百度学术9. 陈柏翰，王笠镔，邹慧辉，王伟光，王可慧. 自旋对弹体侵彻效应的影响. 兵工学报. 2023(S1): 117-124 . 百度学术10. 郭松林，高世桥，李泽章，李云彪. 弹引系统攻角侵彻混凝土仿真与试验研究. 兵器装备工程学报. 2022(01): 135-139+205 . 百度学术11. 成丽蓉，汪德武，贺元吉，赵宏伟，杨杰，邓斌. 弹体斜侵彻多层间隔混凝土薄靶姿态偏转机理研究. 北京理工大学学报. 2022(10): 1009-1016 . 百度学术12. 李鹏程，张先锋，刘闯，魏海洋，刘均伟，邓宇轩. 攻角和入射角对弹体侵彻混凝土薄靶弹道特性影响规律研究. 爆炸与冲击. 2022(11): 92-104 . 本站查看13. 廉璞，牟东，青泽，刘军. 侵彻弹药姿态测量技术研究现状及发展. 探测与控制学报. 2021(02): 1-9 . 百度学术14. 李东伟，刘俞平，王筱锋，王昭明，赵双双. 弹形参数对战斗部斜穿甲姿态偏转影响研究. 兵器装备工程学报. 2021(09): 40-44 . 百度学术15. 张健东，武海军，李伟，李金柱，皮爱国. 头部带肋板的异形结构弹体斜贯穿混凝土薄靶实验和数值模拟. 高压物理学报. 2021(06): 186-194 . 百度学术其他类型引用(5)
-
下载:
下载:
百度学术
计量
- 文章访问数: 6796
- HTML全文浏览量: 2528
- PDF下载量: 756
- 被引次数: 20