偏心起爆对战斗部装药能量分配增益的影响

邓海; 全嘉林; 梁争峰

doi:10.11883/bzycj-2021-0280

偏心起爆对战斗部装药能量分配增益的影响

doi: 10.11883/bzycj-2021-0280

西安近代化学研究所，陕西西安 710065

详细信息

作者简介:
邓　海（1990－　），男，博士研究生，hai179590@163.com

通讯作者:
梁争峰（1972－　），男，博士，研究员，2042lzf@sohu.com

中图分类号: O389
计量
- 文章访问数: 447
- HTML全文浏览量: 226
- PDF下载量: 38
- 被引次数: 3
出版历程
- 收稿日期: 2021-07-05
- 修回日期: 2021-11-11
- 网络出版日期: 2022-03-30
- 刊出日期: 2022-05-27

Influence of eccentric initiation on energy distribution gain of a warhead charge

Xi’an Institute of Modern Chemistry, Xi’an 710065, Shaanxi, China

摘要

摘要: 为研究不同方式的偏心起爆对炸药装药能量分配及增益的影响，建立了偏心起爆战斗部的计算模型，通过局部装填比这一变量，给出了偏心起爆战斗部破片的初速计算公式。采用数值模拟与试验验证结合的方法，对六分位条件下不同偏心起爆方式的破片速度增益和能量增益进行了对比，得出以中心起爆为基准，分别以邻位双线、连位三线、间位双线、偏心单线方式起爆，定向方位内破片的速度增益依次增大；邻位双线起爆时，目标方向破片速度增益达25.47%，定向区域破片动能占总能量的24.57%，能量增益超过40%。
- 偏心起爆 /
- 能量分配 /
- 能量增益 /
- 起爆方式 /
- 破片速度
Abstract: In order to study the influence of different ways of eccentric initiation on the energy distribution and the gain of explosive charge, a theoretical model of eccentric initiation warhead is established, and the concept of energy distribution center is introduced. By introducing the variable of local loading ratio, the calculation formula of initial velocity of fragments of eccentric initiation warhead is formulated. In this paper, the velocity gain of fragments and energy gain with different initiation modes under the sextile condition are compared and analyzed by using numerical simulation and experimental verification. The results show that at the directional orientation, the maximum pressure at the edge of multi-line eccentric initiation is significantly greater than that of eccentric single line initiation and central initiation, and the detonation pressure at the edge of charge increases from 23.5 GPa of central initiation to 36.2 GPa of asymmetrical two lines 60° initiation; The distribution law of fragment velocity in the direction of 0°−30° is similar to the distribution law of maximum pressure at the edge of charge. Taking the central initiation as the benchmark, the relationship of velocity gain with the directional orientation takes the following relation: asymmetrical two lines 60°>asymmetrical three lines 120°>asymmetrical two lines 120°>asymmetrical one line. When asymmetrical two lines 60° initiation, the fragment velocity gain in the target direction is 25.47%. Finally, through the verification of experiments and theoretical calculation, it is concluded that the energy proportion in the directional area of adjacent asymmetrical two lines 60° is the highest, with the energy gain in this area being 47.42%; followed by asymmetrical three lines 120einitiation, with the energy gain being 38.84%; then symmetrical two lines 1202initiation, with the energy gain being 36.98%; and finally asymmetrical one line initiation, with the energy gain in the directional area being 32.72%.
- eccentric initiation /
- energy distribution /
- energy gain /
- initiation mode /
- fragment velocity

HTML全文

对于滑翔式高超声速飞行器而言，由于其气动外形要求，有效载荷空间一般呈现为不对称椭圆、不对称菱形、舌形和梯形等异形截面形状。椭圆截面弹体以其对平台适应性好、有效载荷空间利用率高等优点，成为滑翔式高超声速武器平台的首选方案。诸多学者在椭圆截面弹体对目标的侵彻作用机理方面已开展了大量研究工作^[1-15]。王文杰等^[1]开展了椭圆截面弹体侵彻砂浆靶试验，通过分析弹靶的变形和破坏，发现弹体横截面长轴端点附近区域受力更大，且弹体长短轴参数的改变对侵彻性能影响较显著。Dong等^[2-3]开展了一系列椭圆截面弹体侵彻半无限厚混凝土靶试验，弹体的长短轴比分别为1.5和2.0。试验数据表明，椭圆截面弹体在700～1000 m/s的撞击速度范围内具有良好的侵彻性能和弹道稳定性。同时，通过数值模拟研究了椭圆截面弹体侵彻混凝土的受力特性，结果表明，在相同侵彻速度下，弹体头部的法向应力从短轴向长轴逐渐增大。王浩等^[4-5]开展了椭圆截面弹体正贯穿加筋板试验研究，获得了椭圆截面弹体剩余速度和弹道极限速度的预测公式，得出了随着椭圆截面弹体长短轴比增大，靶板的弹道极限速度近似线性增大的结论。进一步地，田泽等^[6]基于王浩的试验结果，采用能量守恒和虚功原理分阶段分析弹体受力特征，建立了弹体姿态偏转理论模型，并分析了弹体撞击速度、初始倾角、质心位置、弹体翻滚角及截面长短轴比等参数对椭圆变截面弹体姿态偏转的影响。邓希旻等^[7]开展了上下非对称结构异形弹体正/斜贯穿多层间隔921A薄钢板实验，基于Abaqus/Explicit、VUMAT和Python子程序分析了靶体的能量耗散情况和损伤机理。在弹体侵彻阻力方面，Dai等^[8]基于空腔膨胀模型，系统分析了弹体结构参数对椭圆截面弹体侵彻性能的影响，并给出了封闭形式的弹体深度预测公式。刘子豪等^[9]通过数值模拟、理论推导与试验研究相结合的方法，对椭圆截面弹体高速侵彻混凝土特性进行了研究，得出在侵彻速度相同时，弹体侵彻深度随着弹体长短轴比的增大而增大。Liu等^[10-11]开展了椭圆截面弹体侵彻半无限混凝土靶和金属靶试验研究，分析了弹体侵彻阻力与弹靶参数的关系，建立了考虑截面形状影响的阻力修正函数。Ma等^[12]基于修正的空腔膨胀模型，建立了非对称椭圆弹体斜侵彻有限厚金属靶弹道偏转模型，结果表明，通过合理设计弹体结构，非对称椭圆截面弹体可以保持与圆形截面弹体相同的侵彻能力。魏海洋等^[13]研究了椭圆截面弹体对金属厚靶的侵彻弹道规律，基于空腔膨胀理论和局部相互作用模型，建立了椭圆截面弹体侵彻弹道模型，结果表明，椭圆截面弹体侵彻弹道稳定性随长短轴比的增大而变弱，最优长短轴比为1.0。谭远深等^[14]基于152 mm口径轻气炮开展了椭圆截面弹体反弹道侵彻试验研究，并结合数值模拟结果建立了椭圆截面侵彻弹体弯曲结构响应计算模型。刘均伟等^[15]开展了3种不同椭圆截面弹体正贯穿有限厚2A12铝板的试验与数值模拟，结果表明，弹体长短轴比对靶体破坏形貌和响应特性影响显著，椭圆截面弹体侵彻作用下靶体呈现非对称破坏。

综上所述，针对椭圆截面弹体侵彻各类介质靶体，主要研究集中于弹体剩余速度、侵彻深度和弹道稳定性方面，取得了一定的研究成果。但关于圆形和椭圆截面弹体侵彻性能孰优孰劣这一基本问题尚不明确，亟待分析、明确影响椭圆截面弹体侵彻性能的主要因素。

本文中，首先，设计并开展圆锥和椭圆锥压头的准静态深压痕试验，以评估在不考虑动态效应时截面形状对侵彻阻力的影响；然后，开展椭圆截面弹体正侵彻2A12厚铝靶试验，分析靶体的破坏形貌及弹体侵彻深度；最后，结合团队前期研究得到的阻力函数修正系数^[10]和空腔膨胀模型，建立弹体侵彻动力学模型，并系统讨论截面形状和头部曲径比对弹体侵彻性能的影响。

1. 椭圆锥压头压痕试验

圆形和椭圆截面弹体对半无限靶的侵彻性能优劣尚未有明确定论，弹体侵彻时，靶体阻力和弹体结构综合决定了弹体的侵彻性能。由于两者的耦合作用，在动态侵彻试验中很难对其中的单一因素进行定量分析。众所周知，深压痕试验^[16-20]可以反映不同材料对刚性弹体的侵彻阻力。因此，为评估在不考虑动态效应时截面长短轴比对侵彻阻力的影响，开展了圆锥和椭圆锥压头深压痕试验，通过压缩试验机记录的力-位移曲线定量分析不同截面形状压头之间的阻力差异。该部分工作可认为是获得给定目标对刚性侵彻弹体施加的侵彻阻力的基础。

1.1 试验方案

3种类型压头（CC1、EC1和EC2）的长短轴比分别为1.00、1.50和2.00，几何参数如表1所示，其中a和b分别为压头截面的长半轴长和短半轴长，β为弹体长短轴比，定义为β=a/b。圆锥压头的直径为8.00 mm，锥角为21°。椭圆锥压头保持与圆锥压头相同的最大横截面面积、长度和质量，但其截面形状为椭圆。需要说明的是，为了避免压入过程中侧壁摩擦对试验结果的影响，压头直径比压头杆部直径大0.10 mm，如图1所示。压头材料为30CrMnSiNi2A，屈服强度约为1.4 GPa，远高于2A12铝合金试样的屈服强度。2A12铝合金试样直径约为100 mm，高度为80 mm。利用压缩试验机进行静态深压痕试验，记录机器工作时间、压头受力及其位移，控制压缩试验机的速率为0.5 mm/min。压痕试验中压头安装方式及布置如图2所示，通过特定的套筒保持压头的同轴度，确保垂直压入试样。

表 1 圆锥和椭圆锥压头的几何参数

Table 1. Main geometric parameters of circular and elliptical cross-section conical indenters

类型	横截面轮廓	2a/mm	2b/mm	β	最大横截面面积/mm²	长度/mm
CC1		8.00	8.00	1.00	50.27	60.00
EC1		9.80	6.53	1.50	50.27	60.00
EC2		11.31	5.66	2.00	50.27	60.00

下载: 导出CSV

| 显示表格

图 1 CC1压头的设计

Figure 1. CC1 indenter drawing

下载: 全尺寸图片幻灯片

图 2 压头安装方式及现场布局

Figure 2. Installation and layout of the indenter

下载: 全尺寸图片幻灯片

1.2 试验结果与分析

共进行6次压痕试验，得到的力-位移时间曲线如图3所示。可以看出，压痕试验的所有曲线都遵循3个明显的阶段，即图中采用蓝色虚线划分的区域。第1阶段是压头头部接触目标直至完全嵌入目标的过程，该过程中曲线类似于二次函数，此时压入深度约为15 mm。第2阶段则由于压头持续地贯入目标，压头的阻力仍在缓慢增加，力与位移的关系更接近线性。第3阶段为稳定阶段，随着位移的增加，力基本不再变化。这与文献[19-20]报道的曲线略有不同，产生差异的原因是本文中的压头直径比压头杆部直径大0.10 mm，减弱了压入过程中侧壁摩擦对试验结果的影响。

图 3 压痕试验测量的力-位移曲线

Figure 3. Force-distance record in the deep indentation test

下载: 全尺寸图片幻灯片

此外，从图3中可以看出，第2条蓝色虚线位于32 mm处，约为圆锥压头直径的4倍。意味着对于圆锥压头，当压入深度超过32 mm时，力逐渐趋于稳定。而椭圆锥压头的边界线则超过32 mm，这表明椭圆锥压头需要贯入更大的深度才能达到稳定的阻力。

从图3(a)～(c)可以看出，3种压头的深压痕试验数据均具有良好的重复性，微小的差异可以归因于试验中2A12铝样品特性的差异。图3(d)对比了不同长短轴比的椭圆锥压头的力-位移曲线。长短轴比β=2.00的EC2压头始终保持在曲线顶端，其次是β=1.50的EC1压头，最后是β=1.00的CC1压头。此外，2A12铝合金试样对3种压头的作用力分别达到约41.24、40.32和37.45 KN的恒定水平。这意味着β=1.50的椭圆锥压头所需阻力比圆锥压头大6.8%，而β=2.00的椭圆锥压头阻力则大10.1%。

Rosenberg等^[21-22]通过分析刚性弹体侵彻铝靶的试验数据及其数值模拟结果，指出在阈值速度下弹体侵彻阻力表现为靶体强度项，并与弹体结构等若干因素相关联。深压痕试验可以反映不同材料对刚性弹体的侵彻阻力，因此利用深压痕试验获得的数据来评价刚性弹体在阈值速度下的侵彻性能是可行的。可以合理推测，当质量、长度等参数一致时，长短轴比较大的椭圆截面弹体在阈值速度下侵彻性能较差。当然，静态压痕试验和侵彻试验下材料应变率差异很大，且由于铝靶强度是率相关的，即使不考虑惯性项阻力，两种情况下材料侵彻阻力一定会存在部分差异，但仍可对认识截面形状对侵彻阻力的影响起到一定帮助。另一方面，截面形状所影响的靶体侵彻阻力只是决定弹体侵彻性能的部分因素，弹体结构对侵彻性能的影响不可忽略。因此，将通过开展椭圆截面弹体正侵彻2A12厚铝靶试验，综合分析弹体截面长短轴比和头部曲径比对侵彻性能的影响。

2. 椭圆截面弹体正侵彻2A12厚铝靶试验与理论模型

2.1 试验弹体与靶体

为分析弹体撞击速度和截面形状对弹体侵彻性能的影响，开展了椭圆截面弹体正侵彻2A12厚铝靶试验。弹体材料为30CrMnSiNi2A高强度钢，经淬火处理后洛氏硬度HRC为42，屈服强度约为1.2 GPa。试验中使用了2种椭圆截面弹体（E1和E2弹体）和1种圆形截面弹体（C1弹体），3种弹体的结构参数如图4和表2所示。表中：L为弹体头部长度，m为弹体质量，r_CRH为头部曲径比。应当指出，对于椭圆截面弹体，严格来讲并没有头部曲径比的概念。椭圆截面弹体的r_CRH与弹体方位角呈正相关，沿圆周从短轴方向到长轴方向逐渐减小。因此，在表2中列出了弹体在短轴方向上的r_CRH，其他方向上的r_CRH则可以根据弹体的几何方程计算。为保证侵彻后弹体的完整性，靶体选用塑性好、强度较低的2A12铝合金。同时，为减小边界效应对侵彻过程的影响，靶体直径选取为300 mm，厚度为400 mm，为弹体直径的13倍，如图5所示。一般而言，在该尺寸下可以忽略靶体侧面和背面边界效应对弹体侵彻过程的影响。

图 4 试验弹体

Figure 4. Projectiles used in the test

下载: 全尺寸图片幻灯片

表 2 圆形和椭圆截面弹体的主要结构参数

Table 2. Main parameters of circular and elliptical cross-section projectiles

类型	弹体轮廓	β	最大横截面面积/mm²	m/g	L/mm	r_CRH
C1		1.00	556	360	43.2	3.60
E1		1.25	720	360	43.2	3.49
E2		1.61	556	360	43.2	5.64

下载: 导出CSV

| 显示表格

图 5 试验靶体

Figure 5. Targets used in the test

下载: 全尺寸图片幻灯片

图6为试验现场布置，采用高速摄像机观察弹体飞行姿态，同时2个平面镜分别放置于靶前和靶后弹体运动轨迹的正下方，与平面呈45°夹角，通过测量并分析弹体和平面镜中弹体的位置变化，得到弹体的速度、运动轨迹和姿态。

图 6 试验示意图和现场布局

Figure 6. Layout of penetration test

下载: 全尺寸图片幻灯片

2.2 侵彻试验结果

本次试验共计8发，其中C1弹体2发，E1弹体3发，E2弹体3发，通过调整发射药量控制弹体速度在400～800 m/s范围内。图7给出了E2-2弹体的飞行姿态，根据高速摄影以及平面镜可以核定出弹体的着靶姿态和初始撞击速度，试验结果统计如表3所示，表中：v₀为弹体初始撞击速度，P_e为弹体侵彻深度。从表3可以看出，3组对照试验中各弹体初始撞击速度差异较小，具有一定的对比性；弹体俯仰角最大值为0.53°，偏航角最大为2.79°，基本满足正侵彻条件。图8给出了弹体侵彻前后的对比，可以看出，弹体长度和形状无明显变化，仅表面有一定的摩擦痕迹。因此，椭圆截面弹体正侵彻2A12厚铝靶过程可视为刚体侵彻。

图 7 E2-2弹体飞行姿态

Figure 7. Flying attitude of the E2-2 projectile

下载: 全尺寸图片幻灯片

表 3 弹体侵彻2A12厚铝靶试验结果统计

Table 3. Test results of projectiles penetrating 2A12 thick aluminum targets

编号	类型	着角/(°)	俯仰角/(°)	偏航角/(°)	v₀/(m·s⁻¹)	P_e/mm
1	C1-1	0.39	0.41	0.08	773.1	252.77
2	E1-1	1.18	0.38	0.67	810.1	217.18
3	E2-1	1.49	0.41	1.74	816.2	274.55
4	C1-2	1.44	0.30	0.86	608.0	175.73
5	E1-2	0.95	0.26	1.87	604.2	143.93
6	E2-2	0.43	0.20	1.03	563.8	159.17
7	E1-3	0.41	0.53	2.79	398.3	72.92
8	E2-3	0.25	0.39	0	402.2	99.30

下载: 导出CSV

| 显示表格

图 8 试验前后弹体对比

Figure 8. Comparison of projectiles before and after tests

下载: 全尺寸图片幻灯片

2.3 椭圆截面弹体侵彻动力学模型

球形空腔膨胀（sphere cavity expansion，SCE）模型常用于计算刚性弹体侵彻各类介质靶体的侵彻深度/剩余速度^[23-32]。因此，结合Liu等^[10]研究得到的阻力函数修正系数k_r和空腔膨胀模型，建立考虑截面形状影响的弹体侵彻动力学模型。图9给出了椭圆截面弹体的头部几何形状，其中a₁和b₁分别为弹体的长半轴和短半轴长，x₀、y₀和z₀为弹体在3个轴线的坐标，r₁为弹体最大椭圆横截面的极径，r₀为任意椭圆横截面的极径，R_s为头部任意点处的曲率半径，θ₀为任意卵形曲线与x₀Oy₀平面之间的方位角，χ₀为卵形头部表面的外法线向量与弹轴之间的夹角，η₀为弹体表面微元投影至x₀Oy₀平面上与x₀轴的夹角。

图 9 椭圆截面弹体的头部几何形状示意图

Figure 9. Nose geometry diagram of the elliptical cross-section projectile

下载: 全尺寸图片幻灯片

从图9可知，椭圆截面弹体的几何函数可以表示为：

$f\left( {{x_0},{y_0},{{\textit{z}}_0}} \right) = x_0^2 + y_0^2 + {\textit{z}}_0^2 + \frac{{{L^2}\sqrt {{{\left( {{b_1}{x_0}} \right)}^2} + {{\left( {{a_1}{y_0}} \right)}^2}} }}{{{a_1}{b_1}}} - \frac{{{a_1}{b_1}\left( {x_0^2 + y_0^2} \right)}}{{\sqrt {{{\left( {{b_1}{x_0}} \right)}^2} + {{\left( {{a_1}{y_0}} \right)}^2}} }} - {L^2} = 0$

(1)

由几何关系可知，弹体微元表面的法向速度v_n与撞击速度v₀之间的关系为：

${v_{\rm{n}}} = {v_0}\cos {\chi _0}$

(2)

其中， ${\ \chi _0}$ 可以根据几何函数求出：

$\cos {\chi _0} = \frac{{{f_{{{\textit{z}}_0}}}\left( {{x_0},{y_0},{{\textit{z}}_0}} \right)}}{{\sqrt {f_{{x_0}}^2\left( {{x_0},{y_0},{{\textit{z}}_0}} \right) + f_{{y_0}}^2\left( {{x_0},{y_0},{{\textit{z}}_0}} \right) + f_{{{\textit{z}}_0}}^2\left( {{x_0},{y_0},{{\textit{z}}_0}} \right)} }}$

(3)

式中： ${f}_{{x}_{0}}\left({x}_{0},{y}_{0},{{\textit{z}}}_{0}\right)$ 、 ${f}_{{y}_{0}}\left({x}_{0},{y}_{0},{{\textit{z}}}_{0}\right)$ 和 ${f}_{{{\textit{z}}}_{0}}\left({x}_{0},{y}_{0},{{\textit{z}}}_{0}\right)$ 分别为 $f\left({x}_{0},{y}_{0},{{\textit{z}}}_{0}\right)$ 对x₀、y₀和z₀的偏导数。

基于球形空腔膨胀理论提出的描述弹塑性应变硬化不可压缩材料的动态空腔边界应力 $\sigma _{\rm{r}}$ ^[33]可表示为：

$\frac{{{\sigma _{\rm{r}}}}}{{{\sigma _{\text{Y}}}}} = {k_{\rm{r}}}_1 + {k_{\rm{r}}}_2 {\left( {\frac{{{v_{\rm{n}}}}}{{\sqrt {{{{\sigma _{\rm{Y}}}} \mathord{\left/ {\vphantom {{{\sigma _{\rm{Y}}}} {{\rho _{\rm{t}}}}}} \right. } {{\rho _{\rm{t}}}}}} }}} \right)^2}$

(4)

${k_{{\rm{r}}1}} = \frac{2}{3}\left[ {1 + {{\left( {\frac{{2{E_{\rm{t}}}}}{{3{\sigma _{\rm{Y}}}}}} \right)}^{{n_{\rm{t}}}}}K} \right]$

(5)

$K = \int_0^{1 - {{3{\sigma _{\rm{Y}}}} \mathord{\left/ {\vphantom {{3{\sigma _{\rm{Y}}}} ({2{E_{\rm{t}})}}}} \right. } ({2{E_{\rm{t}})}}}} {\frac{{{{\left( { - \ln \delta } \right)}^{{n_{\rm{t}}}}}}}{{1 - \delta }}{\rm{d}} \delta }$

(6)

式中：σ_Y、ρ_t、E_t和n_t分别为靶体的压缩屈服强度、密度、杨氏模量和应变硬化指数；k_r1和k_r2为无量纲系数，k_r2=1.5；δ为积分变量；K为积分常数。

进一步，椭圆截面弹体头部的微元所受应力需要考虑应力分布的不均匀性，因此椭圆截面弹体头部表面应力σ_n可由传统球形空腔膨胀所得应力σ_r，即式(4)，并结合阻力函数修正系数k_r表示：

${\sigma _{\rm{n}}} = {k_{\rm{r}}}{\sigma _{\rm{r}}}$

(7)

阻力函数修正系数k_r的函数表达式为：

${k_{\rm{r}}} = {A_{\text{r}}} {\sin ^2}{\theta _0} + {B_{\text{r}}} {\cos ^2}{\theta _0}$

(8)

式中：A_r和B_r为常数，受椭圆空腔长短轴比、材料特性等影响。

图10给出了圆形和椭圆空腔的阻力函数修正系数k_r随角度的变化情况。可以看出，圆形空腔的修正系数恒为1，而椭圆空腔的修正系数沿着短轴方向至长轴方向逐渐增大，类似于正弦分布。

图 10 阻力函数修正系数

Figure 10. Correction coefficient for resistance function

下载: 全尺寸图片幻灯片

弹体表面径向应力确定后，弹体轴向所受的阻力可以通过式(7)对弹体头部微元S_r积分获得：

${F_{\textit{z}}} = \mathop {\bigcirc \iint}\limits_\varSigma {{\sigma _{\rm{n}}}{\text{d}}} {S_{\rm{r}}}$

(9)

当曲面积分转化为二重积分时，式(9)可表示为：

${F_{\textit{z}}} = \iint\limits_{\tfrac{{x_0^2}}{{a_1^2}} + \tfrac{{y_0^2}}{{b_1^2}} {\text{≤}} 1} {{\sigma _{\rm{n}}}\sqrt {1 + {\textit{z}}_{{x_0}}^2\left( {{x_0},{y_0}} \right) + {\textit{z}}_{{y_0}}^2\left( {{x_0},{y_0}} \right)} }\cos {\chi _0}{\text{d}}{x_0}{\text{d}}{y_0}$

(10)

式中： ${{\textit{z}}}_{{x}_{0}}\left({x}_{0},{y}_{0}\right)$ 和 ${{\textit{z}}}_{{y}_{0}}\left({x}_{0},{y}_{0}\right)$ 分别为式(1)对x₀和y₀所求的偏导数。

联立式(2)～(10)，可以得到弹体侵彻过程中的轴向阻力：

${F_{\textit{z}}} = {C_1} + {C_2}v_{\rm{n}}^2$

(11)

式中：C₁和C₂为积分常数。

根据牛顿第二定律，弹体运动方程可以表示为：

$m\frac{{{\text{d}}v(t)}}{{{\text{d}}t}} = - {F_{\textit{z}}} = - \left( {{C_1} + {C_2}v_{\rm{n}}^2} \right)$

(12)

式中：v(t)为弹体瞬时速度。

对式(12)进行积分，得到弹体瞬时速度与时间的关系：

$v\left( t \right)=\sqrt{\frac{{{C}_{1}}}{{{C}_{2}}}}\tan \left[ {{\arctan }}\left( \sqrt{\frac{{{C}_{2}}}{{{C}_{1}}}}{{v}_{\rm{n}}} \right)-\frac{\sqrt{{{C}_{1}}{{C}_{2}}}}{m}t \right]$

(13)

给定弹体初始速度v₀，弹体侵彻深度P_e可以表示为：

${P_{\rm{e}}} = \frac{m}{{2{C_2}}}\ln \left( {1 + \frac{{{C_2}v_0^2}}{{{C_1}}}} \right)$

(14)

3. 椭圆截面弹体侵彻动力学模型验证及影响因素分析

3.1 模型可靠性验证

表4为弹体侵彻深度计算结果与试验结果的对比，可以看出，理论计算结果与试验结果吻合较好，误差在15%以内。对比分析传统球形空腔膨胀模型和本文中模型计算结果与试验结果的误差可知，本文中考虑椭圆截面形状的理论模型计算结果误差更小，更能反映弹靶接触时真实的受力状况。图11给出了弹体侵彻深度计算结果与试验结果的对比情况。可以看出，无论是传统球形空腔膨胀模型，还是本文中考虑椭圆截面形状的理论模型计算结果，从侵彻性能来看，均为E2弹体的侵彻性能最优，其次为C1弹体，E1弹体的侵彻性能最差。对于圆形截面弹体，由于阻力函数修正系数k_r恒为1，本文中模型与传统球形空腔膨胀模型的计算结果一致，2条曲线完全重合；而对于E1和E2两种椭圆截面弹体，本文中模型计算结果均在传统球形空腔膨胀模型计算结果之下，表明考虑截面形状的影响后，椭圆截面弹体的侵彻性能有所下降。

表 4 弹体侵彻深度计算结果与试验结果的误差

Table 4. Deviation of penetration depth between calculation results and test results

弹体编号	v₀/(m·s⁻¹)	P_e/mm			相对试验结果的误差/%
弹体编号	v₀/(m·s⁻¹)	试验结果	SCE模型	本文模型	SCE模型	本文模型
C1-1	773.1	252.77	285.98	285.98	13.14	13.14
E1-1	810.1	217.18	239.32	229.05	10.19	5.47
E2-1	816.2	274.55	319.24	289.24	16.28	5.35
C1-2	608.0	175.73	191.95	191.95	9.23	9.23
E1-2	604.2	143.93	150.33	144.25	4.45	0.22
E2-2	563.8	159.17	172.56	157.41	8.41	−1.11
E1-3	398.3	72.92	79.61	76.86	9.17	5.40
E2-3	402.2	99.30	100.17	92.25	0.88	−7.10

下载: 导出CSV

| 显示表格

图 11 理论模型计算结果与试验结果的对比

Figure 11. Comparison of the theoretical model calculation results with test results

下载: 全尺寸图片幻灯片

进一步对比传统球形空腔膨胀模型计算的C1和E2弹体的侵彻深度随弹体撞击速度的曲线可知，由于具有相同的截面面积，因此采用传统球形空腔膨胀模型计算的两者曲线基本一致，而具有较大横截面面积的E1弹体的侵彻性能低于C1和E2弹体，这表明最大横截面面积是衡量弹体侵彻性能的一大因素。然而，传统球形空腔膨胀模型计算的E1和E2弹体的侵彻深度相较于试验结果偏大，而考虑椭圆截面对应力分布的影响情况，即在理论模型中加入k_r后，侵彻深度的计算结果更接近于试验结果，这意味着截面形状对弹体侵彻性的影响同样不可忽略。

3.2 弹体长短轴比对侵彻性能的影响

椭圆与圆形截面弹体之间最显著的差异是其截面长短轴比不同。为排除曲径比对侵彻性能的影响，设计了长短轴比分别为1.00、1.25、1.50、1.75和2.00的5种椭圆锥弹体，弹体结构参数如表5所示。除长短轴比外，5种弹体的质量、长度和最大横截面面积均一致。图12给出了撞击速度在200～1000 m/s时5种弹体的侵彻深度。可以看出，在相同撞击速度下，长短轴比越大的弹体侵彻深度越低。具体而言，长短轴比为2.00的椭圆锥弹体的侵彻深度比长短轴比为1.00的圆锥弹体低约13%。

表 5 5种椭圆锥弹体结构的质量参数

Table 5. Main parameters of five typical elliptical cone projectiles

类型	弹体轮廓	2a/mm	2b/mm	β	m/g	L/mm
E_C-1		23.60	23.60	1.00	360	144
E_C-2		26.38	21.11	1.25	360	144
E_C-3		28.90	19.26	1.50	360	144
E_C-4		31.22	17.84	1.75	360	144
E_C-5		33.36	16.68	2.00	360	144

下载: 导出CSV

| 显示表格

图 12 不同长短轴比下侵彻深度与撞击速度的关系

Figure 12. Relationships between penetration depth and impact velocity at different major-to-minor axis length ratios

下载: 全尺寸图片幻灯片

为分析侵彻性能差异的原因，提取撞击速度为800 m/s时5种弹体典型截面沿周向的径向应力分布，截面位置分别为距弹尖L/3和2L/3处。阻力函数修正系数k_r如图13所示，可以看出，椭圆锥弹体的径向应力分布与圆锥弹体明显不同，呈现出从短轴方向到长轴方向逐渐增大的趋势，类似于正弦分布。然而，椭圆锥的径向应力分布曲线并非严格按照基准线（k_r=1）对称。图13中曲线所包围的面积即该截面所提供的弹体阻力。5种弹体的曲线积分面积与圆锥弹体曲线积分面积的比值分别为1.00、1.02、1.06、1.12和1.19。这说明椭圆截面影响了弹体的受力分布，导致侵彻阻力增大，且增加比例随着截面长短轴比的增大而增大。

图 13 不同长短轴比弹体头部的应力分布特性

Figure 13. Stress distribution characteristics of the projectile noses at different major-to-minor axis length ratios

下载: 全尺寸图片幻灯片

实际上，尖卵形弹体的侵彻性能受到多个因素的影响。当弹体最大横截面面积和头部长度保持一致时，长短轴比和头部曲径比会协同变化。长短轴比影响弹体不同截面的径向应力分布，但这仅是侵彻阻力的基础。在计算侵彻阻力时，需要对弹体表面微元进行曲面积分，其中弹体头部表面的外法线向量和弹轴之间的角度对积分值的影响是不可忽略的，而这与弹体头部曲径比息息相关。

3.3 弹体头部曲径比对侵彻性能的影响

圆形截面弹体的侵彻性能很大程度上取决于头部曲径比。Frew等^[32]的研究结果表明，头部曲径比较大的圆形截面弹体具有更好的侵彻性能。为研究不同头部曲径比下椭圆截面弹体的侵彻性能，设计了5种典型椭圆截面弹体，其结构参数如表6所示。为控制头部曲径比的大小，E_CS-4和E_CS-5弹体的头部长度略短于前3种弹体，但5种弹体的质量、最大横截面面积和长度均保持一致。同时，需要注意的是，表中标注的为尖卵形弹体短轴方向的头部曲径比。

表 6 5种椭圆截面弹体结构的质量参数

Table 6. Main parameters of five typical elliptical cross-section projectiles

类型	弹体轮廓示意图	2a/mm	2b/mm	β	r_CRH	m/g	L/mm
E_CS-1		23.60	23.60	1.00	3.60	360	144
E_CS-2		28.90	19.26	1.50	5.28	360	144
E_CS-3		33.36	16.68	2.00	6.95	360	144
E_CS-4		28.90	19.26	1.50	3.60	360	144
E_CS-5		33.36	16.68	2.00	3.60	360	144

下载: 导出CSV

| 显示表格

图14给出了5种弹体侵彻深度随撞击速度的变化曲线。从图14可以看出，2个对照组（E_CS-2和E_CS-4、E_CS-3和E_CS-5）表现出相似的规律，即具有相同长短轴比但头部曲径比更小的E_CS-4和E_CS-5弹体的侵彻性能分别略低于E_CS-2和E_CS-3弹体，这意味着对于尖卵形椭圆截面弹体而言，头部曲径比减小降低了弹体的侵彻性能。

图 14 头部曲径比和长短轴比对侵彻深度的影响

Figure 14. Influence of r_CRH and β on penetration depth

下载: 全尺寸图片幻灯片

此外，相对于E_CS-3和E_CS-5弹体，E_CS-1弹体具有相同的头部曲径比和更小的阻力函数修正系数，因此其侵彻性能优于E_CS-3和E_CS-5弹体。然而，相对于E_CS-2和E_CS-3弹体，E_CS-1弹体在头部曲径比方面处于劣势，而从图14中可以看到，E_CS-1弹体侵彻性能更优异。这说明在长短轴比1～2的范围内，截面形状对侵彻性能的减弱可能强于头部曲径比对侵彻性能的增益。综合上述分析可知，当弹体截面面积相当且各横截面保持长短轴比不变时，圆形截面弹体的侵彻性能最优，椭圆截面弹体侵彻性能随着截面长短轴比的增大和头部曲径比的减小而降低。

4. 结　论

开展了圆锥和椭圆锥压头的准静态深压痕试验和椭圆截面弹体以400～800 m/s撞击速度正侵彻2A12厚铝靶试验，获得了椭圆锥压头的静阻力和椭圆截面弹体侵彻性能的对比结果。结合阻力函数修正系数k_r和球形空腔膨胀模型，建立了弹体侵彻动力学模型，并系统讨论了截面形状和头部曲径比对弹体侵彻性能的影响，得出以下主要结论。

(1) 准静态深压痕试验表明，长短轴比更大的椭圆锥压头贯入材料时所受的阻力更大。具体而言，β=1.50和β=2.00的椭圆锥压头所需阻力比圆形截面压头分别大6.8%和10.1%。

(2) 结合阻力函数修正系数k_r建立的考虑截面形状影响的椭圆截面弹体侵彻动力学模型的计算结果与试验结果的最大相对误差为15%，在侵彻深度预测精度方面优于传统球形空腔膨胀模型。

(3) 当弹体截面面积相当，且各横截面保持长短轴比不变时，圆形截面弹体具有最优的侵彻性能，椭圆截面弹体的侵彻性能随着截面长短轴比的增大和头部曲径比的减小而降低。

图 1 邻位两点偏心起爆后爆轰波相互作用

Figure 1. Interaction between two detonation waves produced by two-point eccentric initiation with an interval of 60°

下载: 全尺寸图片幻灯片

图 2 偏心起爆弹轴中心截面示意图

Figure 2. Schematic diagram of the eccentric detonator shaft center section

下载: 全尺寸图片幻灯片

图 3 战斗部模型

Figure 3. The warhead model

下载: 全尺寸图片幻灯片

图 4 六分圆起爆方式

Figure 4. Initiation modes of sextant circle

下载: 全尺寸图片幻灯片

图 5 不同起爆方式下爆轰波的传播过程

Figure 5. Propagation processes of detonation waves with different initiation modes

下载: 全尺寸图片幻灯片

图 6 不同起爆方式下，爆轰波压力峰值随方位角分布曲线

Figure 6. Detonation wave pressure peak varying with azimuth angle under different initiation modes

下载: 全尺寸图片幻灯片

图 7 不同起爆方式下破片速度随方位角的分布

Figure 7. Fragment velocity varying with azimuth angle under different initiation modes

下载: 全尺寸图片幻灯片

图 8 不同起爆方式下动能分配比随破片方位角的变化

Figure 8. Kinetic energy distribution ratio varying with azimuth angle under different initiation modes

下载: 全尺寸图片幻灯片

图 9 试验样弹的结构设计

Figure 9. Structure design of test sample projectile

下载: 全尺寸图片幻灯片

图 10 靶场布局

Figure 10. Shooting range layout

下载: 全尺寸图片幻灯片

图 11 样弹1爆轰过程的高速摄影

Figure 11. High-speed photography of detonation process of test bomb 1

下载: 全尺寸图片幻灯片

图 12 邻位双线起爆过程的状态

Figure 12. The states of the initiation process of two wires in adjacent positions

下载: 全尺寸图片幻灯片

表 1 定向区不同方位破片速度

Table 1. Fragment velocity in different directions of orientation area

起爆方式	目标方位		方位角30°		定向区内平均
起爆方式	初速/(m·s⁻¹)	速度增益/%	初速/(m·s⁻¹)	速度增益/%	初速/(m·s⁻¹)	速度增益/%
中心起爆	2113.0		2113.0		2113.0
偏心单线	2487.6	17.72	2267.6	7.32	2397.6	13.47
邻位双线	2651.3	25.47	2389.3	13.08	2568.4	21.55
间位双线	2572.2	21.73	2313.0	9.47	2424.9	14.76
连位三线	2583.3	22.25	2319.5	9.77	2441.3	15.54

下载: 导出CSV

表 2 不同区域内破片总动能的分布情况

Table 2. Total kinetic energy distribution of fragments in different regions

起爆方式	0～30°区域		30°～90°区域		90°～150°区域		150°～180°区域
起爆方式	动能占比/%	动能增益/%	动能占比/%	动能增益/%	动能占比/%	动能增益/%	动能占比/%	动能增益/%
中心起爆	16.67		33.33		33.33		16.67
偏心单线	22.12	32.72	33.64	0.93	28.06	−15.81	16.17	−2.99
邻位双线	24.57	47.42	34.07	2.22	26.99	−19.02	14.36	−13.85
间位双线	22.83	36.98	36.65	9.96	26.16	−21.51	14.35	−13.92
连位三线	23.14	38.84	35.80	7.41	26.52	−20.43	14.54	−12.78

下载: 导出CSV

表 3 局部装填比随方位角的拟合关系式

Table 3. Fitting relationship between local loading ratio and azimuth angle

起爆方式	$\beta '$
偏心单线	$\beta (1.015 + 0.451{\text{ }}\cos \theta )$
邻位双线	$\beta ( - 0.207 + 1.923\cos \theta )$
间位双线	$\beta (0.009 + 1.511{\text{ }}\cos \theta )$
连位三线	$\beta (1.252 + 0.228{\text{ }}\cos \theta )$

下载: 导出CSV

表 4 样弹起爆方式

Table 4. Initiation modes of test bombs

样弹	起爆方式	轴向起爆点数
1	连位三线	轴向中心一点
2	间位两线
3	邻位两线
4	偏心单线
5	中心起爆

下载: 导出CSV

表 5 不同周向起爆方式下样弹壳体速度分布

Table 5. Velocity distributions of tested bomb shells under different initiation modes

样弹	起爆方式	壳体速度/(m∙s⁻¹)
		靶板1方向		靶板2、6方向	靶板3、5方向	靶板4方向
		试验	理论计算	试验	试验	试验
1	连位三线	794.6	849.9	708.0	637.5	604.2
2	间位两线	784.3	860.4	687.3	649.4	606.4
3	邻位两线	806.8	909.9	703.4	632.2	602.8
4	偏心单线	763.2	846.0	685.2	654.3	590.3
5	中心起爆	681.5	707.3		681.5

下载: 导出CSV

表 6 不同起爆方式不同方位区域动能分配对比

Table 6. Total kinetic energy distribution in different regions under different initiation modes

样弹	起爆方式	动能占比/%
		靶板1（−30°～30°）		靶板2（30°～90°）		靶板3（90°～150°）		靶板4（−150°～150°）
		试验	模拟	试验	模拟	试验	模拟	试验	模拟
1	连位三线	22.45	23.14	17.83	17.90	14.45	13.26	12.98	14.54
2	间位两线	22.20	22.83	17.05	18.33	15.22	13.08	13.27	14.35
3	邻位两线	23.22	24.57	17.65	17.03	14.26	13.49	12.96	14.36
4	偏心单线	21.37	22.12	17.22	16.82	15.70	14.03	12.78	16.17
5	中心起爆	16.67

下载: 导出CSV

表 7 不同起爆方式下不同方位区域的动能增益

Table 7. Kinetic energy gain in different regions under different initiation modes

样弹	起爆方式	动能增益/%
		靶板1（−30°～30°）		靶板2（30°～90°）		靶板3（90°～150°）		靶板4（−150°～150°）
		试验	模拟	试验	模拟	试验	模拟	试验	模拟
1	连位三线	35.9	38.8	7.9	7.4	−12.5	−20.4	−21.4	−12.7
2	间位两线	32.4	36.9	1.7	9.9	−9.2	−21.5	−20.8	−13.9
3	邻位两线	40.2	47.4	6.5	2.2	−13.9	−19.0	−21.7	−13.9
4	偏心单线	25.4	29.1	1.1	1.5	−7.8	−15.1	−24.9	−1.8

下载: 导出CSV

表 8 壳体速度

Table 8. Shell velocity

起爆方式	壳体速度/(m∙s⁻¹)
	靶板1		靶板2、6		靶板3、5		靶板4
	模拟	试验	模拟	试验	模拟	试验	模拟	试验
连位双线	849.1	794.6	742.6	708.0	703.0	637.5	654.6	604.2
间位三线	859.3	784.3	766.8	687.3	689.8	649.4	652.3	606.4
邻位双线	873.3	806.8	762.0	703.4	683.9	632.2	650.8	602.8
偏心单线	815.0	763.2	728.7	685.2	709.6	654.3	637.6	590.3

下载: 导出CSV

参考文献(18)

[1]	王宝成, 袁宝慧. 防空反导破片杀伤战斗部现状与发展 [J]. 四川兵工学报, 2013, 34(9): 20–24. DOI: 10.11809/scbgxb2013.09.007. WANG B C, YUAN B H. Research states and trend of fragment warhead for air-defense and anti-missile [J]. Journal of Sichuan Ordnance, 2013, 34(9): 20–24. DOI: 10.11809/scbgxb2013.09.007.
[2]	冯顺山, 蒋建伟, 何顺录, 等. 偏轴心起爆破片初速径向分布规律研究 [J]. 兵工学报, 1993(S1): 12–16. FENG S S, JIANG J W, HE S L, et al. Study on radial distribution of initial velocity of off axis bursting disc [J]. Acta Armamentarii, 1993(S1): 12–16.
[3]	RESNYANSKY A D, WILDEGGER-GAISSMAIER A E, KATSELIS G. Directional fragmentation warheads: a theoretical and experimental investigation [C]//Proceedings of the 18th International Symposium on Ballistics. San Antonio, TX, USA, 1999: 543–550.
[4]	KENNEDY D R. A historical review of aimable air defense warhead technology [C]//Proceedings of the 18th International Symposium on Ballistics. San Antonio, TX, USA, 1999: 618–625.
[5]	HELD M. Velocity enhanced warheads [J]. Journal of Explosives and Propellants, 2001, 17(2): 1–12.
[6]	黄静, 孔凡勋, 袁晋, 等. 多点偏心起爆对破片速度增益的影响 [J]. 现代防御技术, 2011, 39(6): 37–42. DOI: 10.3969/j.issn.1009-086x.2011.06.008. HUANG J, KONG F X, YUAN J, et al. Influence of multi-spots off-axis initiation on fragment velocity gain [J]. Modern Defence Technology, 2011, 39(6): 37–42. DOI: 10.3969/j.issn.1009-086x.2011.06.008.
[7]	王树山, 马晓飞, 隋树元, 等. 偏心多点起爆战斗部破片飞散实验研究 [J]. 北京理工大学学报, 2001, 21(2): 177–179. DOI: 10.3969/j.issn.1001-0645.2001.02.008. WANG S S, MA X F, SUI S Y, et al. Experimental research on fragments dispersion of the warhead under asymmetrical multi-spots initiation [J]. Journal of Beijing Institute of Technology, 2001, 21(2): 177–179. DOI: 10.3969/j.issn.1001-0645.2001.02.008.
[8]	刘琛, 李元, 李燕华. 偏心起爆方式对棱柱形定向战斗部破片飞散规律的影响 [J]. 含能材料, 2017, 25(1): 63–68. DOI: 10.11943/j.issn.1006-9941.2017.01.011. LIU C, LI Y, LI Y H. Influence of eccentric initiation ways on fragment dispersion rule of prismatic aimable warhead [J]. Chinese Journal of Energetic Materials, 2017, 25(1): 63–68. DOI: 10.11943/j.issn.1006-9941.2017.01.011.
[9]	LI Y, LI Y H, LIU C, et al. The initiation parameter of detonation wave aiming warhead [J]. Chinese Journal of Energetic Materials, 2016, 24(9): 915–921. DOI: 10.11943/j.issn.1006-9941.2016.09.017.
[10]	卢芳云, 李翔宇, 林玉亮. 战斗部结构与原理 [M]. 北京: 科学出版社, 2009: 112–129.
[11]	王马法, 卢芳云, 李翔宇. 爆轰波斜冲击作用下破片飞散特性研究 [J]. 国防科技大学学报, 2013, 35(1): 60–64. DOI: 10.3969/j.issn.1001-2486.2013.01.012. WANG M F, LU F Y, LI X Y. Research on the projection characteristics of fragments under the loading of the oblique shock wave [J]. Journal of National University of Defense Technology, 2013, 35(1): 60–64. DOI: 10.3969/j.issn.1001-2486.2013.01.012.
[12]	王力, 韩峰, 陈放, 等. 偏心对称起爆战斗部破片初速的增益 [J]. 爆炸与冲击, 2016, 36(1): 69–74. DOI: 10.11883/1001-1455(2016)01-0069-06. WANG L, HAN F, CHEN F, et al. Fragments’ velocity of eccentric warhead with double symmetric detonators [J]. Explosion and Shock Waves, 2016, 36(1): 69–74. DOI: 10.11883/1001-1455(2016)01-0069-06.
[13]	李元. 偏心起爆定向战斗部若干理论与技术研究 [D]. 北京: 北京理工大学, 2016. LI Y. Study on theories and technologies of asymmetrically initiated warhead [D]. Beijing, China: Beijing Institute of Technology, 2016.
[14]	沈慧铭. 多点起爆方式作用机理及其在战斗部中的应用研究 [D]. 南京: 南京理工大学, 2018. SHEN H M. Research on action mechanism of multi-point initiation way and its application in warhead [D]. Nanjing, Jiangsu, China: Nanjing University of Science and Technology, 2018.
[15]	李翔宇, 李振铎, 梁民族. D型战斗部破片飞散性及威力场快速计算 [J]. 爆炸与冲击, 2019, 39(4): 043301. DOI: 10.11883/bzycj-2017-0420. LI X Y, LI Z D, LIANG M Z. Dispersion properties and rapid calculation of fragment force field of D-shaped fragmentation warhead [J]. Explosion and Shock Waves, 2019, 39(4): 043301. DOI: 10.11883/bzycj-2017-0420.
[16]	彭军, 袁宝慧, 孙兴昀, 等. 立方形钨破片穿甲效应的数值模拟与试验 [J]. 爆破器材, 2017, 46(5): 23–28. DOI: 10.3969/j.issn.1001-8352.2017.05.005. PENG J, YUN B H, SUN X Y, et al. Numerical simulation and experiment for penetrating performance of cubic tungsten fragment [J]. Explosive Materials, 2017, 46(5): 23–28. DOI: 10.3969/j.issn.1001-8352.2017.05.005.
[17]	李刚, 刘荣忠, 郭锐. 爆炸成抛物面型弹丸战斗部的仿真研究 [J]. 计算机仿真, 2011, 28(11): 1–4; 18. DOI: 10.3969/j.issn.1006-9348.2011.11.001. LI G, LIU R Z, GUO R. Simulation on parabolic liner forming process [J]. Computer Simulation, 2011, 28(11): 1–4; 18. DOI: 10.3969/j.issn.1006-9348.2011.11.001.
[18]	周翔. 爆炸成型弹丸战斗部的相关技术研究[D]. 南京: 解放军理工大学工程兵工程学院, 2006. ZHOU X. Technological researches on explosive formed projectiles [D]. Nanjing, Jiangsu, China: Engineering Institute of Engineering Institute of Engineer Corps of PLA University of Science and Technology, 2006.

施引文献

期刊类型引用(3)

1.	薛颖杰，陈智刚，杨芮，崔晋，康彦姝，付建平. 椭圆截面弹体侵彻多层混凝土靶的数值仿真. 中北大学学报(自然科学版). 2025(01): 35-44 . 百度学术
2.	董恒，黄风雷，武海军，邓希旻，李萌，刘龙龙. 异形弹体高速侵彻/穿甲机理研究进展. 兵工学报. 2024(09): 2863-2887 . 百度学术
3.	戴湘晖，王可慧，周刚，李明，沈子楷，李鹏杰，段建，吴海军，杨慧. 刚性椭圆截面弹体侵彻性能理论研究. 现代应用物理. 2024(06): 132-141 . 百度学术