基于物质点法的轨道炮刨削机理三维数值研究

吴金国; 林庆华; 弯港; 金涌; 李海元; 栗保明

doi:10.11883/1001-1455(2017)02-0307-08

基于物质点法的轨道炮刨削机理三维数值研究

doi: 10.11883/1001-1455(2017)02-0307-08

南京理工大学瞬态物理国家重点实验室, 江苏南京 210094

基金项目:

国家自然科学基金 11402266

中央高校基本科研业务费专项资金项目 1151210420

详细信息

作者简介:
吴金国(1989-)，男，博士研究生，wujg8848@163.com

中图分类号: O347
计量
- 文章访问数: 4510
- HTML全文浏览量: 1293
- PDF下载量: 400
- 被引次数: 0
出版历程
- 收稿日期: 2015-07-07
- 修回日期: 2016-02-29
- 刊出日期: 2017-03-25

3D numerical research of railgun gouging mechanism based on material point method

National Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China

摘要

摘要: 基于轨道炮结构特点以及冲击热力学理论，采用物质点法建立了轨道微颗粒诱发刨削的三维模型，模拟了轨道刨削的形成过程，并对其产生机理与影响因素进行了分析。结果表明：电枢与轨道的局域高速冲击产生瞬时的能量交换，形成的高热高压金属流对轨道表面的斜侵彻作用形成了刨坑；刨削的产生存在速度阈值，超过速度阈值，随着电枢速度增加，刨削越严重；低于速度阈值，可产生轨道擦伤；减小轨道表面微颗粒尺寸、增加电枢头部倾角均可降低刨削损伤。
- 冲击动力学 /
- 物质点法 /
- 轨道炮 /
- 刨削 /
- 超高速
Abstract: Based on the structural features of the railgun and the theories of impact thermodynamics, a 3D gouging model containing a micro particle was established using the material point method to simulate the formation process of the rail gouging, and the gouging mechanism and some influencing factors were also analyzed. The results show that the local impact between the armature and the rails at a high velocity produces transient energy exchanges, thus simultaneously forming metal flows with high energy and high pressure that penetrate obliquely into the rail and cause the formation of gouging, for which there is a threshold velocity to produce. As the armature velocity increases, the gouging damage gets more serious. On the other hand, the galling damage occurs when the armature velocity is below the threshold velocity. Both controlling the particle size within a certain range and increasing the head angle of the armature will help to suppress the formation of gouging.
- impact dynamics /
- material point method /
- railgun /
- gouging /
- hypervelocity

HTML全文

爆炸成型弹丸（explosively-formed projectile，EFP）在侵彻装甲后形成的破片会对装甲内部人员和仪器造成毁伤^[1-2]，因此通过靶板厚度和EFP着靶速度预测靶后破片质量、数量以及速度分布特性对EFP靶后毁伤效应评估、工程防护与设计具有重要意义。而准确计算EFP垂直侵彻有限厚装甲钢板时靶板和EFP产生的靶后破片质量，是预测上述分布特性的重要基础。

在侵彻问题中，Alekseevskii^[3]和Tate^[4-7]所建立的A-T模型取得了很好的效果，很多学者对此模型进行了改进和发展^[8-13]，但所涉及侵彻体的主体部分都是细长等截面圆柱杆，没有考虑小长径比和变截面的情况。文献[3-9]在研究长杆弹侵彻问题时忽略了蘑菇头轴向长度对侵彻过程的影响；文献[14-15]在研究EFP侵彻问题时忽略了蘑菇头轴向长度对侵彻过程的影响；文献[15]在研究EFP垂直侵彻RHA（rolled homogenous armour）靶板时，忽略了EFP的强度，直接套用了射流的扩孔模型：这些都是与实际情况不符的。目前对于靶后效应的理论研究均是基于侵彻体的主体部分是细长等截面圆柱杆的情况^[2,16]，对于EFP靶后效应的研究多集中在仿真和试验方面^[17-19]，而在EFP靶后效应的理论研究中均未考虑EFP变截面的特性^[1,14-15]。EFP主要用于末敏弹的战斗部，这类EFP具有变截面的特性，因此不能直接套用主体部分是细长等截面圆柱杆的侵彻理论；由于其长径比小，因此不能忽略蘑菇头轴向长度对侵彻过程的影响；由于其头部形状不规则，因此在侵彻有限厚钢靶板后产生的冲塞体形状不确定；由于其炸高很大（可达1 000倍装药直径），因此不能像射流一样忽略其强度。基于上述原因，使得利用现有的理论计算EFP垂直侵彻有限厚装甲钢板时靶板和EFP产生的靶后破片质量将会带来较大的偏差。

因此，本文将在考虑了EFP蘑菇头轴向长度、EFP强度的基础上，对冲塞体形状进行假设，并运用理论的方法分别分析考虑和不考虑EFP变截面的特性时，靶板和EFP产生的靶后破片质量。本文将借助文献[20]中的的试验数据验证理论分析中的关键参量，借助文献[18-19]中的的试验数据验证仿真方法的可信性，并利用此仿真方法对某典型EFP垂直侵彻装甲钢板时靶板和EFP产生的靶后破片质量进行计算，利用所得的仿真结果分析和判断EFP变截面的特性对靶板和EFP产生的靶后破片质量的影响。

1. 质量模型

1.1 基本假设

冲塞体是当EFP侵彻有限厚靶板到一定深度时，在靶板背面形成的一个内部布满裂纹的金属块^[16]，部分消蚀的EFP内部也布满裂纹。冲塞体和部分消蚀的EFP内部的裂纹会随着侵彻的进行，逐渐延伸扩展，当侵彻完成，冲塞体和部分消蚀的EFP将在靶板后形成若干碎块，即为靶板和EFP产生的靶后破片。为了简化侵彻过程，本文中做如下假设：（1）EFP开坑阶段所需的时间短、消耗的能量少，因此不计EFP开坑对侵彻过程的影响。（2）初始扩孔的轴向压力与径向压力相等，径向压力与扩孔面积成反比。（3）针对本文所涉及的EFP垂直侵彻装甲钢板的情况，假设冲塞体为圆台且母线与底面夹角45°。理由是：Backman等^[21]指出脆性靶板产生靶后破片时靶后出孔较大呈喇叭状；Yarin等^[16]在研究钨合金长杆弹侵彻装甲板时在其理论中指出靶后出孔较大呈喇叭状；Arnold等^[22]和Dalzell等^[17]分别运用仿真的方法发现射流和EFP垂直侵彻装甲板时靶后出孔较大呈喇叭状；王昕等^[18]、李睿等^[15]在运用试验的方法研究EFP垂直侵彻装甲板时发现侵彻出孔呈锥形，外侧直径较大；张先锋等^[14]在研究EFP对有限厚靶板侵彻过程及后效中认为冲塞崩落块为圆台且母线与底面夹角45°。

1.2 变截面EFP对装甲钢板的侵彻模型

邢柏阳等^[23]、Xing等^[24]在之前的研究中发现，密实结构EFP具有旋转对称的特性，可以用任意通过其旋转轴的截面表示其形状。本文中根据某典型EFP的激光高速照相测量结果（如图1上半部分所示），绘制了EFP垂直侵彻装甲钢板的过程示意图，如图2所示，H₀为靶板厚度，H为EFP未消蚀部分和消蚀部分交界面距离靶板底面的距离，R表示侵彻孔（蘑菇头）半径，u为EFP的侵彻速度，v表示EFP未消蚀部分的速度，r为EFP未消蚀部分和消蚀部分交界面处的EFP截面半径，l为EFP未消蚀部分和消蚀部分交界面距离EFP尾部的距离，P为EFP的侵彻深度。

图 1 某典型爆炸成型弹丸的形状

Figure 1. The shape of a certain typical explosively-formed projectile

下载: 全尺寸图片幻灯片

图 2 爆炸成型弹丸（EFP）垂直侵彻靶板的过程

Figure 2. The process of an explosively-formed projectile (EFP) normally penetrating into a target

下载: 全尺寸图片幻灯片

由于EFP具有变截面的特性且没有函数表达式可以描述其横截面半径与其轴向位置的对应关系，因此只能根据EFP的激光高速照相测量结果获取横截面半径（r_i）与其轴向位置（l_i）关系的离散点数据，记为r（l_i），l为EFP未消蚀部分的长度，由于离散点的数量足够大，因此 ${l_{i + 1}} - {l_i} \ll l$ ，可以认为dl=l_i+1−l_i，结合牛顿第二定律，可知EFP未消蚀部分的减速度满足以下关系，

$\quad\quad\quad\quad\quad\frac{{{\rm{d}}v}}{{{\rm{d}}t}} = - \frac{{{Y_{\rm{p}}}r_i^2}}{{\int_0^l {{\rho _{\rm{p}}}r_i^2{\rm{d}}l} }}\quad\quad\quad\quad\quad\quad\quad(1)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

式中：i为自某典型EFP尾部起的离散点序号，由于EFP具有变截面的特性，使得没有函数表达式可以描述其横截面半径与其轴向位置的对应关系，因此通过对激光高速照相图像识别而获得离散点；l_i为第i个离散点处与某典型EFP尾部的距离；r_i为某典型EFP的l_i处外轮廓横截面半径；dl为某典型EFP上任意两个离散间的距离，i足够大， ${\rm{d}}l = {l_{i + 1}} - {l_i} \ll l$ ；t为侵彻时间，从t=0时刻开始侵彻；ρ_p为EFP的密度；Y_p为EFP的强度因子。

不考虑EFP变截面的特性时，即认为EFP轴上所有位置处的横截面半径均相同，其值为前述所有离散点横截面半径（r_i）的平均值：

$r' = \frac{1}{k}\sum\limits_{i = 1}^k {{r_i}}$

(2)

式中：k为离散点的数量。此时EFP未消蚀部分的减速度关系即为A-T模型^[3-7]中的减速度关系：

$\frac{{{\rm{d}}v}}{{{\rm{d}}t}} = - \frac{{{Y_{\rm{p}}}}}{{{\rho _{\rm{p}}}l}}$

(3)

考虑和不考虑EFP变截面的特性时，EFP的消蚀速度、侵彻速度均满足A-T模型^[3-7]：

$\frac{{{\rm{d}}l}}{{{\rm{d}}t}} = - (v - u)$

(4)

$\frac{{\rm{1}}}{{\rm{2}}}{\rho _{\rm{p}}}{(v - u)^2} + {Y_{\rm{p}}} = \frac{{\rm{1}}}{{\rm{2}}}{\rho _{\rm{t}}}{u^2} + {\sigma _{\rm{t}}}$

(5)

式中： ${\rho _{\rm{t}}}$ 、 ${\sigma _{\rm{t}}}$ 分别为靶板的密度和强度因子。

由于本文所研究的某典型EFP的材料为紫铜，因此 ${Y_{\rm{p}}} \text{＜} {\sigma _{\rm{t}}}$ ，由式（5）可以得到：

$u = \frac{{v - \sqrt {\displaystyle\frac{{{\rho _{\rm{t}}}}}{{{\rho _{\rm{p}}}}}{v^2} + \left( {1 - \displaystyle\frac{{{\rho _{\rm{t}}}}}{{{\rho _{\rm{p}}}}}} \right)v_{{\rm{cr}}}^2} }}{{1 - \displaystyle\frac{{{\rho _{\rm{t}}}}}{{{\rho _{\rm{p}}}}}}}$

(6)

${v_{{\rm{cr}}}} = \sqrt {\frac{{2({\sigma _{\rm{t}}} - {Y_{\rm{p}}})}}{{{\rho _{\rm{p}}}}}}$

(7)

式中：v_cr为临界侵彻速度。当v＞v_cr时，EFP可以侵彻靶板；当v≤v_cr时，EFP停止侵彻靶板。

分析图2中的受力关系，得到EFP对靶板的轴向作用力满足：

${F_{\rm{t}}} = {\rm{\pi }}{R^2}\left( {\frac{1}{2}{\rho _{\rm{t}}}{u^2} + {\sigma _{\rm{t}}}} \right)$

(8)

由假设（2）可以得到初始扩孔的径向压力p₀、侵彻过程中扩孔的径向压力p与扩孔半径的关系：

${p_{\rm{0}}} = \frac{1}{2}{\rho _{\rm{p}}}{({v_{\rm{0}}} - {u_{\rm{0}}})^2} + {Y_{\rm{p}}}$

(9)

$p = \frac{{{r^2}}}{{{R^2}}}{p_{\rm{0}}}$

(10)

式中：u₀为t=0时刻的侵彻速度，v₀为EFP着靶速度，r为EFP消蚀部分和未消蚀部分交界面处的横截面半径。可以发现，r是影响EFP-靶板交界面处压力的关键因素。另外，r的取值视情况而定，当考虑EFP变截面的特性时取r_i，当不考虑EFP变截面的特性时取 $r'$ ，后文中r的取值规则与此相同。

侵彻过程中扩孔的径向压力可以使用含有靶板强度因子的流体力学Bernoulli方程得到:

$p = \frac{1}{2}{\rho _{\rm{t}}}{\left( {\frac{{{\rm{d}}R}}{{{\rm{d}}t}}} \right)^2} + {\sigma _{\rm{t}}}$

(11)

由式（9）～（11）可以得到扩孔速度：

$\frac{{{\rm{d}}R}}{{{\rm{d}}t}} = \sqrt {\frac{{{r^2}}}{{{R^2}}}\frac{{{\rho _{\rm{p}}}{{({v_{\rm{0}}} - {u_{\rm{0}}})}^2} + 2{Y_{\rm{p}}}}}{{{\rho _{\rm{t}}}}} - \frac{{2{\sigma _{\rm{t}}}}}{{{\rho _{\rm{t}}}}}}$

(12)

由式（12）可以发现，扩孔速度与R负相关，与r和v₀正相关。当考虑EFP变截面的特性时，随着侵彻进行：结合图1可以发现r先快速增大再缓慢变化，然后又快速增大，再缓慢变化；结合式（12）可知R先快速增大（bI阶段），再缓慢增大（bII阶段，R≈R_bII（常数）），然后又快速增大（bIII阶段），再缓慢增大（bIV阶段）。当不考虑EFP变截面的特性时，随着侵彻的进行，r恒为定值，由于扩孔速度恒为正值，但r/R的值逐渐减小，因此R先快速增大（dI阶段），然后一直保持缓慢增大（dII阶段，R≈R_dII（常数）），理论分析并计算可得R_bII≈R_dII。

1.3 靶板和EFP产生的靶后破片质量

冲塞体形成时刻 t_c的各参数关系如图3所示，网格表示冲塞体内部的碎块，H_c为冲塞体高度，R_c为冲塞体形成时刻的侵彻孔（蘑菇头）半径，u_c为冲塞体形成时刻的EFP的侵彻速度，v_c为冲塞体形成时刻的EFP未消蚀部分的速度，r_c为冲塞体形成时刻的EFP未消蚀部分和消蚀部分交界面处的EFP截面半径，l_c为冲塞体形成时刻的EFP未消蚀部分和消蚀部分交界面距离EFP尾部的距离，P_c为冲塞体形成时刻的EFP的侵彻深度。

图 3 冲塞体形成时刻

Figure 3. Forming time of a plug

下载: 全尺寸图片幻灯片

根据假设（3），冲塞体为圆台且母线与底面夹角45°，因此冲塞体形成时受到的轴向剪切力如式（13）左边所示，此时其值与EFP对靶板轴向的作用力相等，如式（13）右边所示，即：

$\sqrt 2 {\text{π}}{H_{\rm{c}}}({H_{\rm{c}}} + 2{R_{\rm{c}}}){\tau _{\max }} = {\text{π}}R_{\rm{c}} ^2\left( {\frac{1}{2}{\rho _{\rm{t}}}u_{\rm{c}}^2 + {\sigma _{\rm{t}}}} \right)$

(13)

式中：τ_max为靶板的最大剪切强度。τ_max的计算方法如下：

${\tau _{\max }} = {\tau _0}\gamma _i^n\left[ {1 + {{\left( {\frac{{{{\dot \gamma }_{\rm{m}}}}}{{q'}}} \right)}^{1/q}}} \right]\exp \left\{ { - \frac{{\alpha \beta }}{{{\rho _{\rm{t}}}{C_V}}}\left[ {1 + {{\left( {\frac{{{{\dot \gamma }_{\rm{m}}}}}{{q'}}} \right)}^{1/q}}} \right]\frac{{{\tau _0}}}{{n + 1}}\gamma _i^{n + 1}} \right\}$

(14)

${\dot \gamma _{\rm{m}}} = \frac{{1 - n}}{n}\frac{u}{R}$

(15)

${\gamma _i} = {\left\{ {\frac{{n{\rho _{\rm{t}}}{C_V}}}{{{\tau _0}\alpha \beta \left[ {1 + {{\left( {\displaystyle\frac{{{{\dot \gamma }_{\rm{m}}}}}{{p'}}} \right)}^{1/q}}} \right]}}} \right\}^{1/(n + 1)}}$

(16)

式中：τ₀为γ=1时材料的剪切强度；n为材料的应变硬化指数；q和 $q'$ 为材料的应变率敏感指数；α为应力的温度软化系数；C_V为材料的定容比热容；β为塑形功转化为热量的比例系数（Taylor-Quinney系数^[25]）； ${\gamma _i}$ 为剪应变； ${\dot \gamma _{\rm{m}}}$ 为平均剪应变率。

式（12）～（13）表明：r可以影响冲塞体形成时刻的侵彻孔（蘑菇头）半径（R_c），并进一步影响冲塞体形成时刻的冲塞体高度（H_c）。

考虑EFP蘑菇头轴向方向的长度对侵彻过程的影响，有如下关系：

${P_{\rm{c}}} + {H_{\rm{c}}} - {H_0} - {R_{\rm{c}}} = 0$

(17)

${P_{\rm{c}}} = \int_0^{{t_{\rm{c}}}} {u{\rm{d}}t}$

(18)

式（17）～（18）表明r将会影响冲塞体形成时刻的侵彻深度（P_c），并最终决定冲塞体的形成时刻（t_c）。因此冲塞体质量，即靶板产生的靶后破片质量，满足以下关系：

${m_{\rm{t}}} = {\rho _{\rm{t}}}\left\{ {\frac{1}{3}{\text{π}}{H_{\rm{c}}}[R_{\rm{c}}^2 + {{({R_{\rm{c}}} + {H_{\rm{c}}})}^2} + {R_{\rm{c}}}({R_{\rm{c}}} + {H_{\rm{c}}})] - \frac{2}{3}{\text{π}}R_{\rm{c}}^3} \right\}$

(19)

由式（13）和（19）可以得到：

${m_{\rm{t}}} = \frac{{\text{π}}}{3}{\rho _{\rm{t}}}R_{\rm{c}}^3f({u_{\rm{c}}})$

(20)

在本文讨论的靶板厚度（30 mm≤H₀≤70 mm）和EFP着靶速度（1 650 m/s≤v₀≤1 860 m/s）范围内，f（u_c）＞0恒成立，因此m_t与R_c总是呈正相关，并且当R_c变化幅度较大时，R_c是影响m_t的主要因素；当R_c变化幅度较小时，u_c是影响m_t的主要因素。

冲塞体形成时EFP总的剩余质量即为EFP产生的靶后破片质量，即：

${m_{\rm{p}}} = {\rho _{\rm{p}}}{\text{π}}\int_0^{{l_{\rm{c}}} + {R_{\rm{c}}}} {{r^2}} {\rm{d}}l$

(21)

基于上述理论分析可以发现，r通过影响扩孔速度，进而影响冲塞体形成时刻的各个参量（R_c、H_c、P_c、t_c等），最终影响靶板和EFP产生的靶后破片质量。本文将基于某典型EFP，对不同靶板厚度和EFP着靶速度条件下，考虑和不考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，进行分析和比较，具体地验证改进后的模型可以减小与实际情况的偏差。

2. 模型验证

2.1 侵彻孔半径模型验证

Yarin等^[16]在通过试验研究钨合金长杆弹侵彻装甲板时发现，在统计靶板的质量损失时结果波动很大，接近30%，因此本文中将利用试验中可以稳定测量的参数验证力学模型的准确性。由于Held等^[20]在侵彻孔半径模型方面所做的工作非常经典，并且射流侵彻水介质和EFP侵彻装甲钢都属于侵彻体强度小于靶板强度的模式，因此将文献[20]中的初始条件代入式（12）中，得到的侵彻孔半径（R）随侵彻时间（t）的变化曲线如图4中的红线所示，黑点为Held等^[20]的理论结果。

图 4 侵彻孔半径随时间变化关系

Figure 4. The time history of the crater radius

下载: 全尺寸图片幻灯片

由图4可以发现，利用本文理论模型所得的侵彻孔半径（R）随侵彻时间（t）的变化关系与Held等^[20]的理论结果十分吻合，表明本文理论模型具有较高的准确性。

2.2 仿真模型验证

Dalzell等^[17]采用AUTODYN-3D^TM对EFP侵彻靶板形成破片云的过程进行了数值模拟，指出SPH（smoothed particle hydrodynamics）算法较Lagrange算法和Euler算法在研究靶后破片方面更具优势。在试验中所统计的靶板质量损失波动很大^[16]，并且，经常发现靶板和EFP产生的靶后破片会粘结在一起。因此，本文首先利用试验中可以稳定测量的结果^[18-19]验证本文仿真方法的可信性，进而利用得到验证的仿真方法对某典型EFP进行仿真计算，验证理论模型所得的、但在试验统计中波动较大的结果。

本文使用AUTODYN-3D^TM的SPH算法，仿真模型中EFP和靶板的材料种类、状态方程以及本构方程取自文献[19,26]。文献[18-19]的试验中，EFP战斗部由药型罩、主装药、传爆药、壳体组成，其装药直径（D）为56 mm，药型罩为壁厚3 mm的球缺形紫铜罩，装药为长径比0.86的压装JH-2，壳体厚度为0.045D。EFP垂直侵彻厚度为20 mm的45钢板，EFP的形状、尺寸如图5下半部分的X光照片所示，试验中测得EFP着靶速度为2 120 m/s，45钢板的长、宽均为200 mm。在仿真中构建与文献[18-19]的试验中完全一致的初始条件，仿真中EFP的形状、尺寸如图5上半部分所示，并对靶板边界施加固定约束，EFP侵彻45钢板前的仿真模型如图6所示。

图 5 依据文献[18-19]构建的EFP模型

Figure 5. The EFP model developed according to references [18-19]

下载: 全尺寸图片幻灯片

图 6 依据文献[18-19]构建的仿真模型

Figure 6. The simulation model developed by according to references [18-19]

下载: 全尺寸图片幻灯片

仿真所得侵彻80 µs和95 µs时的靶后破片云形状如图7（a）、（b）上半部分所示，文献[18]中给出的相应时刻的X光照片如图7（a）、（b）下半部分所示，对比发现仿真所得的靶后破片云图和文献[18]中试验所得的X光照片非常吻合。另外，仿真所得侵彻95 µs时的靶后破片云径向膨胀速度v_J=531 m/s,靶后破片云头部膨胀速度v_T=1 690 m/s与文献[18]中的试验结果v_J=525 m/s，v_T=1 678 m/s的误差不超过2%，表明本文所用的仿真方法是可信的。

图 7 本文仿真所得的靶后破片云与文献[18]的对比

Figure 7. Behind-armor debris clouds obtained by this paper compared with that in reference [18]

下载: 全尺寸图片幻灯片

3. 对比分析

图1中某典型EFP战斗部装药直径（D）为125 mm，装药长度为0.84D，药型罩为壁厚0.064D的球缺形紫铜罩，主装药、传爆药均为8701，周侧壳体厚度为0.048D，底部壳体厚度为0.12D，由西安现代控制技术研究所设计并由辽沈工业集团有限公司生产。某典型EFP的最大半径为28.5 mm，最大长度为98 mm，其形状及尺寸均在图1中标明，其激光高速照相测量结果如图1上半部分所示，其仿真中的模型如图1下半部分所示，并对靶板边界施加固定约束。使用上述仿真方法计算某典型EFP以着靶速度1 650 m/s垂直侵彻厚度为70 mm的装甲钢板，仿真得到图8所示的靶板出孔形状，其中黑色虚线表示靶板出孔侧面，黄色虚线表示冲塞体形成前EFP的侵彻路径，可以发现，靶板出孔与靶板底面夹角接近45°，与文献[14-18,21-22]比较吻合，因此假设（2）是具有可信性的。

图 8 某典型EFP垂直侵彻靶板的出孔形状(单位：mm)

Figure 8. The exit hole shape of a target normally penetrated by a certain typical EFP (unit in mm)

下载: 全尺寸图片幻灯片

质量模型中的参数全部取自文献[7,13,19,25-26]，Y_p和R_t的计算方法由文献[7]给出，主要参数为：τ₀=0.907 GPa，n=0.25， $q'$ =8.0×10⁷ s⁻¹，q=7.33，α=6.8×10⁻⁴，β=0.9，C_V=452 J/（kg·K），靶板的杨氏模量E_t=210 GPa，Tate给出的计算Y_p时的系数λ=0.7，ρ_p=8.93 g/cm³，ρ_t=7.85 g/cm³。

3.1 靶板厚度对靶后破片质量的影响

由式（19）和（21）可以发现，冲塞体的形状参数直接决定了靶板产生的靶后破片质量，间接影响了EFP产生的靶后破片质量。由式（18）可以发现，冲塞体的形状参数与靶板厚度密切相关，因此靶板厚度是影响靶后破片质量的主要因素之一。在某典型EFP着靶速度为1 650 m/s的情况下，垂直侵彻厚度为30、40、50、60、70 mm的装甲钢板，靶板和EFP产生的靶后破片质量如图9所示。

图 9 不同靶板厚度条件下的靶后破片质量

Figure 9. Mass of behind-armor debris for different thicknesses of targets

下载: 全尺寸图片幻灯片

通过编写Matlab^TM程序进行图像识别，读取某典型EFP的激光高速照相测量结果（如图7上半部分所示），辨识其边界点，从而获取横截面半径与其轴向位置关系的离散点数据，对式（20）和（21）进行数值求解，可以得到考虑EFP变截面的特性时靶板和EFP产生的靶后破片质量，分别记为m_t，bjm、m_p，bjm。利用所得的离散点数据以及式（2），求出r'=16.23 mm，对式（20）和（21）进行数值求解，可以得到不考虑EFP变截面的特性时靶板和EFP产生的靶后破片质量，分别记为m_t，djm、m_p，djm。使用本文仿真方法所得的靶板和EFP产生的靶后破片质量，分别记为m_t，sim、m_p，sim。

由图9可以发现，相比于不考虑EFP变截面特性的结果，考虑EFP变截面特性的结果与仿真结果更吻合，能够更加准确地解释靶板和EFP产生的靶后破片质量随靶板厚度的变化规律。

考虑EFP变截面的特性时，随着靶板厚度的增大，侵彻时间不断增长，在本节研究的靶板厚度范围内R_c逐渐由bII阶段过渡到bIII阶段，最后进入bIV阶段。由于bIII阶段R_c快速增大，由理论分析可知，此时R_c是影响m_t的主要因素，结合式（20）可知靶板产生的靶后破片质量先基本不变，然后快速增加，最后缓慢增加。随着靶板厚度的增大，EFP的消蚀长度不断增大，因此EFP产生的靶后破片质量不断减小。

不考虑EFP变截面的特性时，由理论分析并计算可得R_bII≈R_dII，并且在本节研究的靶板厚度范围内R_c始终处于dII阶段，因此对于厚度较小的靶板，考虑与不考虑EFP变截面所得的靶板产生的靶后破片质量接近；但随着靶板厚度的增大，考虑EFP变截面的R_c开始快速超过不考虑EFP变截面的R_c，因此考虑与不考虑EFP变截面所得的靶板产生的靶后破片质量差不断增大。由于此时R_c的变化极其缓慢，由理论分析可知，u_c是影响m_t的主要因素，随着靶板厚度的增大，u_c不断减小，但由于u_c对m_t的影响比较小，因此不考虑EFP变截面的特性时，靶板产生的靶后破片质量缓慢减小。由于不考虑EFP变截面时EFP的头部质量大于考虑EFP变截面时EFP的头部质量，这使得在侵彻开始阶段不考虑EFP变截面时EFP损失的质量较大，因此不考虑EFP变截面所得的EFP产生的靶后破片质量总是小于考虑EFP变截面所得的结果。

计算以上5种工况条件下，考虑和不考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，如表1所示。举例，考虑EFP变截面的特性时所得的靶板产生的靶后破片质量与仿真结果的偏差的具体计算方法为：ε_t，bjm=|（m_t，bjm−m_t，sim）/m_t，sim|×100%，另外3个偏差的计算方法与此类似。

表 1 不同靶板厚度条件下靶后破片质量的偏差

Table 1. Mass deviations of behind-armor debris for different thicknesses of targets

H₀/mm	ε_t，bjm/%	ε_t，djm/%	ε_p，bjm/%	ε_p，djm/%
30	4.1	13.5	3.5	32.0
40	1.8	12.3	2.9	40.3
50	0.8	25.9	4.1	45.4
60	7.4	55.4	7.4	54.8
70	6.0	70.1	10.4	64.3

下载: 导出CSV

| 显示表格

由表1可以发现，考虑EFP变截面的特性时所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，均明显小于相同靶板厚度条件下不考虑EFP变截面的特性时所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，并且，随着靶板厚度的增大，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断增强。表明，相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随靶板厚度的变化规律。

3.2 EFP着靶速度对靶后破片质量的影响

由式（12）可以发现，EFP着靶速度对扩孔速度有显著的影响，而扩孔速度与冲塞体形成的临界条件密切相关，因此EFP着靶速度是影响靶后破片质量的主要因素之一。在装甲钢板厚度为40 mm的情况下，某典型EFP分别以1 650、1 680、1 740、1 800、1 860 m/s的着靶速度进行垂直侵彻，靶板和EFP产生的靶后破片质量如图10所示。

图 10 不同EFP着靶速度条件下的靶后破片质量

Figure 10. Mass of behind-armor debris for different impact velocities of EFPs

下载: 全尺寸图片幻灯片

由图10可以发现，相比于不考虑EFP变截面特性的结果，考虑EFP变截面特性的结果与仿真结果更吻合，能够更准确地解释靶板和EFP产生的靶后破片质量随EFP着靶速度的变化规律。

考虑EFP变截面的特性时，由于靶板厚度固定为40 mm，当EFP着靶速度为1 650 m/s时，侵彻时间最长，但是由3.1节中的分析可知，此时的R_c处于bII阶段，即R_c缓慢增大，此时影响R_c的主要因素是EFP的着靶速度，次要因素是侵彻时间。因此在本节研究的EFP着靶速度范围内，随着EFP着靶速度的增高，R_c不断增大，同时u_c不断增大，结合式（20）可知靶板产生的靶后破片质量不断增大。由于侵彻时间不断减短，因此EFP的消蚀长度不断减小，因此EFP的剩余长度不断增大，使得EFP产生的靶后破片质量不断增加。

同理，不考虑EFP变截面的特性时，此时的R_c处于dII阶段，影响R_c的主要因素是EFP的着靶速度。因此在本节研究的EFP着靶速度范围内，随着EFP着靶速度的增高，靶板产生的靶后破片质量不断增大。由于侵彻时间不断减短，使得EFP产生的靶后破片质量不断增加。

计算以上5种工况条件下，考虑和不考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，如表2所示。

表 2 不同EFP着靶速度条件下靶后破片质量的偏差

Table 2. Mass deviations of behind-armor debris for different impact velocities of EFPs

v₀/（m·s⁻¹）	ε_t，bjm/%	ε_t，djm/%	ε_p，bjm/%	ε_p，djm/%
1 650	1.8	12.3	2.9	40.3
1 680	3.9	11.2	2.4	38.6
1 740	4.9	9.6	1.2	37.1
1 800	2.7	5.0	1.2	36.1
1 860	2.2	2.3	0.4	34.1

下载: 导出CSV

| 显示表格

由表2可以发现：考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，均明显小于相同EFP着靶速度条件下不考虑EFP变截面的特性时所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差；并且，随着EFP着靶速度的提高，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断减弱，当EFP着靶速度达到1 860 m/s时，考虑和不考虑EFP变截面的特性所得的靶板产生的靶后破片质量十分接近。表明，相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随EFP着靶速度的变化规律。

4. 结　论

本文中建立了变截面爆炸成型弹丸（explosively-formed projectile，EFP）垂直侵彻装甲钢板靶后破片质量模型，并得到以下结论：（1）相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随靶板厚度和EFP着靶速度的变化规律，可以减小与实际情况的偏差；（2）当EFP着靶速度为1 650 m/s时，随着靶板厚度从30 mm增大到70 mm，靶板产生的靶后破片质量先保持基本不变，然后快速增大，再缓慢增大，EFP产生的靶后破片质量不断减小，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断增强；（3）当靶板厚度为40 mm时，随着EFP着靶速度从1 650 m/s提高到1 860 m/s，靶板和EFP产生的靶后破片质量均不断增大，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断减弱。

图 1 轨道炮实验中典型的泪滴状刨坑形貌^[4]

Figure 1. Typical tear-drop gouging morphology in railgun^[4]

下载: 全尺寸图片幻灯片

图 2 轨道刨削物质点法计算模型

Figure 2. A simulation model of railgun gouging constructed with material point method

下载: 全尺寸图片幻灯片

图 3 轨道刨坑形成过程Mises应力云图

Figure 3. Mises stress distribution during the formation process of gouging

下载: 全尺寸图片幻灯片

图 4 枢轨接触面三个方向接触力变化曲线

Figure 4. Contact forces in three directions varying with time

下载: 全尺寸图片幻灯片

图 5 2个质点受冲击时的速度变化曲线

Figure 5. Velocities of two particles varying with time during impacting

下载: 全尺寸图片幻灯片

图 6 1.5 km/s冲击速度下数值模拟出的泪滴状刨坑

Figure 6. Tear-drop gouge craters from numerical simulation at 1.5 km/s

下载: 全尺寸图片幻灯片

图 7 2.0 km/s速度下数值模拟出的泪滴状刨坑

Figure 7. Gouging morphology from simulation at 2.0 km/s

下载: 全尺寸图片幻灯片

图 8 1.0 km/s速度下模拟出的轨道擦伤沟槽

Figure 8. A galling track from simulation at 1.0 km/s

下载: 全尺寸图片幻灯片

图 9 冲击峰值压强随速度的变化

Figure 9. The peak impact pressure at different velocities

下载: 全尺寸图片幻灯片

图 10 实验中轨道上出现的磨损擦伤沟槽^[6]

Figure 10. Galling tracks on the rail from an experiment^[6]

下载: 全尺寸图片幻灯片

表 1 电枢与轨道材料参数

Table 1. Material parameters of armature and rail

材料	ρ/(kg·m^-3)	E/GPa	ν	c_p/(J·kg^-1·K^-1)	A/MPa	B/MPa	n	C	m	T_m/K	T_r/K	${\dot \varepsilon }$ ₀/s^-1	χ	D₁	D₂	D₃	D₄	D₅	c₀/(m·s^-1)	s	γ₀
7075铝	2 810	71	0.33	960	369	684	0.73	0.008 3	1.7	933	293	1	0.9	0.13	0.13	-1.5	0.011	0	5 350	1.34	2
无氧铜	8 960	124	0.34	383	90	292	0.31	0.025	1.09	1 356	293	1	0.9	0.54	4	2	0.014	1.12	3 940	1.49	2

下载: 导出CSV

参考文献(15)

[1]	李军, 严萍, 袁伟群.电磁轨道炮发射技术的发展与现状[J].高电压技术, 2014, 40(4):1052-1064. http://d.old.wanfangdata.com.cn/Periodical/gdyjs201404014 Li Jun, Yan Ping, Yuan Weiqun. Electromagnetic gun technology and it's development[J]. High Votage Engineering, 2014, 40(4):1052-1064. http://d.old.wanfangdata.com.cn/Periodical/gdyjs201404014
[2]	Graff K F, Dettloff B B. The gouging phenomenon between metal surfaces at very high sliding speeds[J]. Wear, 1969, 14(2):87-97. doi: 10.1016/0043-1648(69)90339-1
[3]	Cinnamon J D, Palazotto A N, Szmerekovsky A G. Further refinement and validation of material models for hypervelocity gouging impacts[J]. American Institute of Aeronautics and Astronautics, 2008, 46(2):317-327. doi: 10.2514/1.25035
[4]	Stefani F, Parker J V. Experiments to measure gouging threshold velocity for various metals against copper[J]. IEEE Transactions on Magnetics, 1999, 35(1):312-316. doi: 10.1109/20.738423
[5]	Persad C, Prabhu G, White G, et al. Characterization of hypervelocity gouge craters in rail conductors[J]. IEEE Transactions on Magnetics, 1997, 33(1):401-405. doi: 10.1109/20.560045
[6]	Watt T J, Clay C E, Bassett P M, et al. The effect of surface indentations on gouging in railguns[J]. Wear, 2014, 310(1/2):41-50. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=266be5121f733b3a9b099eab65c32ced
[7]	Barber J P, Bauer D P, Ahmadi R. Contact phenomena at hypervelocities[J]. Wear, 1982, 78(1/2):163-169. doi: 10.1016-0043-1648(82)90028-X/
[8]	Barker L M, Trucano T G, Susoeff A R. Railgun rail gouging by hypervelocity sliding contact[J]. IEEE Transactions on Magnetics, 1989, 25(1):83-87. doi: 10.1109/20.22510
[9]	金龙文, 雷彬, 李治源, 等.轨道炮刨削形成机理分析及数值模拟[J].爆炸与冲击, 2013, 33(5):537-543. doi: 10.3969/j.issn.1001-1455.2013.05.014 Jin Longwen, Lei Bin, Li Zhiyuan, et al. Formation mechanism analysis and numerical simulation of railgun gouging[J]. Explosion and Shock Waves, 2013, 33(5):537-543. doi: 10.3969/j.issn.1001-1455.2013.05.014
[10]	Barker L M, Trucano T G, Susoeff A R. Gun barrel gouging by sliding metal contact at very high velocities[J]. IEEE Transactions on Magnetics, 1989, 1(25):83-87.
[11]	Bourell D L, Persad C. Simulation of railgun gouging[J]. IEEE Transactions on Magnetics, 1999, 35(1):274-276. doi: 10.1109/20.738416
[12]	Zhu Rengui, Zhang Qian, Li Zhiyuan, et al. Impact physics model and influencing factors of gouging for electromagnetic rail launcher[C]//17th International Symposium on Electromagnetic Launch Technology. La Jolla, California, USA, 2014.
[13]	刘峰, 赵丽曼, 张晖辉, 等.电磁轨道炮刨削的形成机理及仿真分析[J].高压物理学报, 2015, 29(3):199-205. http://www.cqvip.com/QK/96553X/201503/71898776504849534851484854.html Liu Feng, Zhao Liman, Zhang Huihui, et al. Formation mechanism and simulation analysis of railgun gouging[J]. Chinese Journal of High Pressure Physics, 2015, 29(3):199-205. http://www.cqvip.com/QK/96553X/201503/71898776504849534851484854.html
[14]	廉艳平, 张帆, 刘岩, 等.物质点法的理论和应用[J].力学进展, 2013, 43(2):237-264. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKC20132013052800030339 Lian Yanping, Zhang Xiong, Liu Yan, et al. Material point method and its applications[J]. Advances in Mechanics, 2013, 43(2):237-264. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKC20132013052800030339
[15]	张雄, 廉艳平, 刘岩, 等.物质点法[M].北京:清华大学出版社, 2013:7-13. Zhang Xiong, Lian Yanping, Liu Yan, et al. Material point method[M]. Beijing: Tsinghua University Press, 2013:7-13.