前抛体对弹体入水载荷影响数值模拟研究

肖睿; 魏继锋; 吉耿杰; 冯宇浪

doi:10.11883/bzycj-2022-0431

前抛体对弹体入水载荷影响数值模拟研究

doi: 10.11883/bzycj-2022-0431

北京理工大学爆炸科学与技术国家重点实验室，北京 100081

详细信息

作者简介:
肖　睿（1997－　），女，硕士研究生，xr1xiaorui@163.com

通讯作者:
魏继锋（1977－　），男，博士，副教授，weijifeng@bit.edu.cn

中图分类号: O359; TJ630
计量
- 文章访问数: 591
- HTML全文浏览量: 183
- PDF下载量: 95
- 被引次数: 3
出版历程
- 收稿日期: 2022-10-08
- 修回日期: 2022-12-08
- 网络出版日期: 2023-01-14
- 刊出日期: 2023-04-05

Numerical research on the effect of front body on water-entry load of a projectile

State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 弹体穿过气-液界面时因密度突变会遭受强冲击载荷，出现结构损伤或破坏。为降低弹体入水的冲击载荷，基于Rabbi降载思想，提出了一种在主弹体前附加前抛体的结构形式。运用S-ALE (structured arbitrary Lagrange-Euler)算法和罚函数流固耦合方法进行空泡形态和弹体运动状态数值模拟，数值模拟结果与弹体入水实验图像吻合，验证了数值模拟算法的有效性；随后研究了前抛弹体的入水角度、主弹体与前抛体的无量纲入水时间间隔参数、前抛体尺寸、主弹体与前抛体的入水初速度对降载效果的影响。研究结果表明，前抛体垂直入水时存在前抛体和主弹体碰撞现象，不利于降载，甚至可能给主弹体带来更大的载荷；前抛体斜入水时可避免两弹体间的碰撞，具有良好的降载效果，最大降载率可达90%。在此基础上获得了降载效果最佳的无量纲入水时间间隔范围为0.8~0.9，在此范围内主弹体的降载效果随前抛体尺寸的增大和入水初速度增加而提高。
- 弹体入水 /
- 冲击载荷 /
- 降载技术 /
- 前抛体 /
- S-ALE
Abstract: When the projectile passes through the gas-liquid interface, the sudden change of density may cause a violent impact load and do untold damage to it. It seriously affects the working effect of the projectile. In order to reduce the water entry load of the projectile, based on Rabbi’s idea of load reduction, a kind of structure of projectile with front body was proposed. The S-ALE (structured arbitrary Lagrange-Euler) algorithm and the fluid-structure coupling method with penalty function were used to simulate the shape of the cavitation wall and the projectile motion state, which is coincident by comparing with those of experiments. The validity of the numerical method is verified. Furthermore, the influence of the water entry angle of front body, the dimensionless water entry time interval parameter between main projectile and front body, the size of front body, the initial water entry velocity of main projectile and front body on impact load were researched by numerical simulation. The simulation results show that the front body will impact the main projectile when they both entering water vertically, which increases the impact load due to the collision of them. When the front body enters the water obliquely and the main projectile still enters the water vertically, the collision between main projectile and front body could be avoided and a good load reduction effect is obtained. The maximum load reduction ratio is up to 90%. The dimensionless time interval parameter range for obtaining a good load reduction effect is from 0.8 to 0.9. Within this range, variation laws of the water entry load of main projectile with the size of the front body and the initial velocity of entering water were discussed in detail. The effect of load reduction increases with the increase of the size of the front body and the initial water entry velocity.
- projectile water entry /
- impact load /
- load reduction technology /
- front body /
- S-ALE

HTML全文

结构体入水过程涉及固-液-气三相耦合作用，具有高度非定常性和非线性，在救生艇下水、航天器水上回收、水上飞机降落、空投鱼雷入水等诸多领域都有着广泛应用。当结构体穿过气-液界面时，由于环境介质的密度突变，其头部会受到剧烈的冲击，结构极易损坏失效^[1]。因此，研究结构体入水载荷特性及其降载技术具有重要的学术意义和工程应用价值。

最早开展的入水冲击载荷实验研究的是Bottomley^[2]进行的V型浮舟模型入水实验，而最早的入水冲击载荷理论研究则是von Karman^[3]提出的楔形体和平板入水冲击理论。Wagner^[4]在von Karman入水冲击模型^[3]的基础上引入水面变形因子，提出了入水冲击问题的自相似解法。Eroshin等^[5]提出钝体结构入水时，若忽略液体可压性将导致入水冲击载荷被高估，因此应考虑液体的可压性。随着对入水载荷问题认识日益深刻，入水载荷描述愈加复杂，理论解析愈发困难。数值模拟程序的研发和数值计算能力的提升，促进了数值模拟技术在入水冲击载荷研究中的广泛应用。

在入水冲击载荷降载方面，常用的降载方法有调整弹体结构、增加缓冲罩壳、加装空化器等。改善弹体头部形状可有效减小入水冲击载荷，斜锥形头部弹体受到的载荷比半球形、锥形、和球台形弹体受到的载荷更小^[6]，椭球头型受到的最大冲击力与速度和头部长短轴比有关^[7]，头部半球角越小入水冲击载荷越小^[8]，独特的海鸟头型可有效减小入水冲击载荷^[9]。缓冲罩壳在入水时的碎裂和缓冲组件的变形吸能可降低入水冲击载荷，最早由Howard^[10]提出，之后的学者又研究了不同材料、结构的罩壳及其缓冲组件的降载效果^[11–18]。弹体头部加装喷气空化器在入水前使水面产生空泡来降低冲击载荷，空化器通气量越大、弹体速度越高，降载效果越好^[19-21]。Rabbi等^[22]提出串联入水降载思路，在进行的两小球前后入水试验中发现，前抛球可使后球的入水冲击载荷下降78%；Lyu等^[23]进行了速度为15～40 m/s的两个相同小球串联入水试验，发现后入水小球的降载率随时间间隔的增大而减小。

目前，对串联入水体降载问题的研究刚刚起步，仅仅分析了两个小球垂直入水的情形，对于其他结构体、不同入水速度、不同入水姿态的入水特性和降载规律尚未见报道。本文据此开展平头圆柱弹体串联入水冲击数值模拟研究，采用S-ALE (structured arbitrary Lagrange-Euler)算法和罚函数流固耦合方法，深入探讨前抛体入水角度、时间间隔、入水速度等参量对入水空泡和入水冲击载荷的影响，探索串联降载机制和降载影响规律，以期为弹体降载新技术的运用提供重要的理论和数据支撑。

1. 含前抛体弹体入水过程及理论分析

根据von Karman的平板入水冲击模型^[3]，平头圆柱形弹体入水时受到的冲击载荷为

$p = {\rho _{\text{w}}}{v_{\text{p}}}{c_{\text{w}}}$

(1)

式中：p为平板入水所受压力；ρ_w为水的密度；v_p为弹体与水自由面间的相对速度；c_w为水中声速。

弹体受到的加速度为

$a = \frac{F}{m} = \frac{{p{A_0}}}{m} = \frac{{{\rho _{\text{w}}}{v_{\text{p}}}{c_{\text{w}}}{A_0}}}{m}$

(2)

式中：a为弹体加速度；F为弹体受到的冲击力；A₀为平头圆柱与水的初始接触面积；m为弹体质量。由式(2)知，决定弹体冲击载荷的因素有水的密度与声速、弹体与水自由面间的相对速度和弹体入水接触面积，因此可通过减小弹体与水面之间相对速度、降低水体密度以及减小初始接触面积的方式来降低弹体入水冲击载荷。

含前抛体的弹体入水时，前抛体改变水体的状态，使原本静止的水体具有了一定的速度，减小了弹体与水体之间的相对速度，降低了主弹体撞击水体的砰击载荷。水体在前抛体的作用下会经历空泡形成、空泡闭合（包括面闭合、浅闭合、深闭合）、射流形成等阶段。与之相对应的，主弹体可能在这之中的任意阶段入水。当然，前抛体和主弹体顺序入水时可能出现相撞的现象。如图1所示，将含前抛体弹体入水过程中会出现的情况归结为以下四种：Ⅰ，撞击前抛体；Ⅱ, 撞击空泡壁面；Ⅲ，撞击飞溅水体；Ⅳ，撞击射流。

图 1 含前抛体弹体入水情况示意图

Figure 1. Projectile with front projectile entering the water

下载: 全尺寸图片幻灯片

2. 数值模拟算法及设置

2.1 S-ALE与流固耦合算法

S-ALE (structured arbitraty Lagrange-Euler)算法于2015年提出，在LS-DYNA中得到运用^[24]，与耦合了Lagrange算法和Euler算法的ALE方法接近^[25]：计算时先执行Lagrange时间步，单元网格随材料产生变形；然后保持变形后的边界条件，对内部单元网格重新划分，网格拓扑关系保持不变；最后将变形网格中的单元变量和节点速度矢量输送到重分后的新网格中^[26]。S-ALE与ALE方法的区别在于，前者的流场采用正交的结构化网格，只需在设置中输入i、j、k三个方向的节点数量及坐标即可由计算机自动生成网格，便于设置且更容易收敛，并行计算效率更高。

S-ALE算法的流场控制方程为

$\frac{{\partial \rho }}{{\partial t}} = - \rho \frac{{\partial {v_i}}}{{\partial {x_i}}} - {w_i}\frac{{\partial \rho }}{{\partial {x_j}}}$

(3)

$\rho \frac{{\partial {v_i}}}{{\partial t}} = \frac{{\partial {\sigma _{ij}}}}{{\partial {x_j}}} + \rho {b_i} - \rho {w_j}\frac{{\partial {v_i}}}{{\partial {x_j}}}$

(4)

$\rho \frac{{\partial e}}{{\partial t}} = {\sigma _{ij}}\frac{{\partial {v_i}}}{{\partial {x_j}}} + \rho {b_i}{v_i} - \rho {w_j}\frac{{\partial e}}{{\partial {x_i}}}$

(5)

式中：ρ为流体密度；v为流体速度；w为流体与网格间的相对速度，w=v−u，u为网格移动速度；σ_i为应力张量，σ_ij=−p_fδ_ij+μ(v_ij+v_ji)，其中p_f为流场中的压力，μ为动力黏性系数，δ_ij为Kronecker函数，当i=j时为1，否则为0；b为单位体积力；e为能量；下标i、j表示不同方向的欧拉坐标。

网格控制方程组为

$\frac{{\partial f\left( {{X_i},t} \right)}}{{\partial t}} = \frac{{\partial \left( {{x_i},t} \right)}}{{\partial t}} + {w_i}\frac{{\partial \left( {{x_i},t} \right)}}{{\partial t}}$

(6)

式中：X_i为拉格朗日坐标。

计算入水问题时气-液界面捕捉是研究的重要内容。S-ALE采用的界面捕捉方法是根据网格单元中不同物质的体积分数梯度方向和大小来确定切割面或分界面的位置，其计算方式如图2所示，图中η为物质的体积分数，上标“−”表示节点的平均体积分数；V为单元体积。方向n由公式(7)确定。

图 2 S-ALE的界面捕捉法

Figure 2. Surface capturing method of S-ALE

下载: 全尺寸图片幻灯片

${\boldsymbol{n}}={\left\|\frac{\partial {\bar{\eta }}_{1}}{\partial \boldsymbol{x}}\right\|}^{-1}\cdot \frac{\partial {\eta }_{1}}{\partial \boldsymbol{x}}$

(7)

流固耦合中固体对流体的作用通过添加在流场控制方程中的“源项”来实现，流场对固体结构的作用通过结构附近的流体质点速度的插值实现^[27]。其控制方程组为

$\rho \frac{{{\text{d}}{\boldsymbol{v}}}}{{{\text{d}}t}} = - \nabla p + \mu \Delta {\boldsymbol{v}} + {\boldsymbol{F}}\left( {{\boldsymbol{x}},t} \right)$

(8)

$\nabla \cdot {{\boldsymbol{v}}}=0$

(9)

${\boldsymbol{F}}\left( {{\boldsymbol{x}},t} \right) = \int {{\boldsymbol{f}}\left( {{\boldsymbol{s}},t} \right)\delta \left[ {{\boldsymbol{x}} - {\boldsymbol{X}}\left( {{\boldsymbol{s}},t} \right)} \right]{\text{d}}{\boldsymbol{s}}}$

(10)

$\frac{{\partial {\boldsymbol{X}}}}{{\partial t}}\left( {{\boldsymbol{s}},t} \right) = {\boldsymbol{u}}\left[ {{\boldsymbol{X}}\left( {{\boldsymbol{s}},t} \right),t} \right]$

(11)

式中：F(x,t)为边界传递给流场的体积力；x为流固耦合边界的位移；f为流固耦合边界产生的单位力；X(s,t)为流体与固体耦合边界的位移；s为流体与固体的边界坐标；δ为Dirac delta函数，该函数除0坐标以外函数值为0，在其整个定义域上的积分为1。S-ALE的流固耦合算法可以有效解决ALE“罚函数”耦合方法产生的“渗漏”问题，获得更准确的计算结果。

2.2 计算模型

主弹体和前抛体均为平头圆柱，主弹体尺寸为 $\varnothing$ 6 mm×24 mm，前抛体尺寸为 $\varnothing$ 2 mm×2 mm。流体计算域宽200 mm，水深300 mm，对弹体经过区域的网格进行加密。弹体采用Lagrange网格划分，流体域用ALE网格剖分。由于计算模型的对称性，采用1/2模型进行计算，以主弹体入水点为原点O，水平向右方向为x轴正向，竖直向下为z轴正向，将xOz平面设为对称平面。弹体模型与流体域网格如图3所示。进行计算时，首先设置主弹体速度为0，前抛体以预定速度入水，计算至所需时间间隔时设置主弹体速度，通过LS-DYNA的重启动计算主弹体入水载荷。

图 3 弹体模型与流体域网格

Figure 3. Projectile models and fluid domain mesh

下载: 全尺寸图片幻灯片

2.3 材料模型

计算模型包括流体（水和空气）、主弹体和前抛体四种材料。

水采用NULL材料模型和Grüneisen状态方程：

${p_{\text{w}}} = \dfrac{{{\rho _0}{c_{\text{w}}}^2{\mu _{\text{w}}}\left[ {1 + \left( {1 - \dfrac{{{\gamma _0}}}{2}} \right){\mu _{\text{w}}} - \dfrac{{{a_1}}}{2}\mu _{\text{w}}^{\text{2}}} \right]}}{{{{\left[ {1 - \left( {{S_1} - 1} \right){\mu _{\text{w}}} - {S_2}\dfrac{{\mu _{\text{w}}^2}}{{{\mu _{\text{w}}} + {\text{1}}}}- {S_3}\dfrac{{\mu _{\text{w}}^{\text{3}}}}{{{{\left( {{\mu _{\text{w}}} + 1} \right)}^2}}}} \right]}^2}}} + \left( {{\gamma _0} + {a_1}{\mu _{\text{w}}}} \right){E_{\text{w}}}$

(12)

式中：p_w为水的压力；μ_w为水的体积变化率，μ_w = ρ_w/ρ_w0−1；γ₀为Grüneisen系数；a₁为γ₀的一阶修正；S₁～S₃为u_s-u_p曲线斜率的无量纲系数，u_s为冲击波速度，u_p为水的质点速度；E_w为水的体积内能。表1为水的Grüneisen状态方程参数值。

表 1 水的Grüneisen状态方程参数

Table 1. Grüneisen equation of state parameters of water

ρ_w/(kg·m⁻³)	c_w/(m·s⁻¹)	S₁	S₂	S₂	γ₀	a₁	E_w/(J·m⁻³)	μ_w0
1000	1480	1.979	0	0	0.11	3	3.07×10⁵	1

下载: 导出CSV

| 显示表格

空气采用NULL材料模型和线性多项式状态方程：

${p_{\text{a}}} = {C_0} + {C_1}{\mu _{\text{a}}} + {C_2}\mu _{\text{a}}^{\text{2}} + {C_3}\mu _{\text{a}}^{\text{3}} + \left( {{C_4} + {C_5}{\mu _{\text{a}}} + {C_6}\mu _{\text{a}}^{\text{2}}} \right){E_{\text{a}}}$

(13)

式中：p_a为空气压力；μ_a为空气的体积变化率；C₀～C₆为状态方程系数；E_a为空气的体积内能。用Gamma定律描述气体状态方程时

$C_0=C_1=C_2=C_3=C_6=0$

(14)

$C_4=C_5=\gamma -1$

(15)

式中：γ为绝热常数，空气为1.4。表2为空气的多项式状态方程参数值，其中μ_a0为μ_a的初始值。

表 2 空气的多项式状态方程参数

Table 2. Polynomial equation of state parameters of air

ρ_a/(kg·m⁻³)	C₀	C₁	C₂	C₃	C₄	C₅	C₆	E_a/(J·m⁻³)	μ_a0
1.225	0	0	0	0	0.4	0.4	0	2.5×10⁵	1

下载: 导出CSV

| 显示表格

主弹体和前抛体由钢材料制成，采用刚体材料模型，密度为7 800 kg/m³，杨氏模量为2.1×10¹¹ Pa和泊松比为0.3。

3. 无前抛体时主弹体的入水载荷

当主弹体以106.8 m/s的速度垂直入水时，将弹体头部接触水面的时刻记为时间零点，图4和图5给出了主弹体入水4 ms内的空泡形态及入水深度，与文献[28]的试验结果吻合。虽然文献[28]中并未测得弹体入水时的冲击载荷，但是通过弹体的运动状态可以反映其动力学特征,由此可以判断数值模拟结果可信。图6为主弹体无前抛体入水过程中的受载加速度数值模拟曲线，主弹体入水瞬间的砰击载荷高达806.3 km/s²，与式(1)和式(2)的理论计算结果844.3 km/s²的误差为4.507%。

图 4 弹体入水运动状态和空泡形态

Figure 4. Projectile entering water motion state and cavitation shape

下载: 全尺寸图片幻灯片

图 5 弹体入水深度(z)

Figure 5. Projectile water entry depth (z)

下载: 全尺寸图片幻灯片

图 6 主弹体的加速度(a)

Figure 6. The accelerationof the main projectile (a)

下载: 全尺寸图片幻灯片

入水过程中，弹体头部受到水面冲击的同时，附近的水介质也受到冲击并产生一道压缩波，以声速在水中传播。当压缩波到达水面时在自由表面发生卸载并产生稀疏波，使主弹体受到的冲击载荷卸载，因此从图6中可以看到在极高的初始冲击峰值之后主弹体的冲击加速度快速降低至零。在极短的一段时间后弹体头部压力恢复为水中动压并趋于稳定，弹体在水中开始稳定航行，弹体速度在流体阻力的作用下逐渐衰减。图6中加速度值的波动为流固耦合时耦合节点变化产生的数值波动。

4. 有前抛体时主弹体的入水载荷

4.1 前抛体角度变化对主弹体载荷的影响

4.1.1 前抛体角度变化时的空泡演化特点

将前抛体的初始入水角定义为圆柱轴线与水平面夹角，如图7所示，记为θ；主弹体的载荷变化用降载率α表征，降载率定义为无前抛体时主弹体最大加速度与含前抛体入水时主弹体最大加速度之差与无前抛体时主弹体最大加速度的比值，即

图 7 前抛体初始入水角示意

Figure 7. Schematic of the initial water entry angle of front body

下载: 全尺寸图片幻灯片

$a=1-\frac{{a}_{\mathrm{max}}'}{{a}_{\mathrm{max}}}$

(16)

式中：a_max为主弹体降载前的载荷峰值； ${a}_{\mathrm{max}}'$ 为降载后的载荷峰值。

图8展示了θ为70°～90°时前抛体以100 m/s入水产生的空泡在0～20 ms内的演化过程。前抛体入水后与周围水体产生了剧烈的动量交换，水体被排开形成逐渐扩大的空泡，水面升高并形成了飞溅水体。随着空泡的扩大，为补充空泡内的气体，大气内的空气从空泡开口处流入空泡，开口处因为Bernoulli效应而形成低压区，飞溅水体由于压力差的作用向中心聚拢闭合，并形成了上下两个方向相反的射流。其中，前抛体以θ=70°入水时未产生明显的射流，而是形成了一片较大的水幕，水幕形态不规则且散乱。

图 8 前抛体以不同角度入水形成的空泡演化过程

Figure 8. The evolution process of the cavitation formed by the front body entering the water at different angles

下载: 全尺寸图片幻灯片

此时引入一个入水时间间隔的无量纲参数τ，其表达式为

$\tau =\Delta t/t_{\rm p}$

(17)

式中：Δt为主弹体与前抛体入水的时间间隔，t_p为空泡的闭合时间。由图8可得当θ=90°时t_p=13 ms；当θ=80°时t_p=13 ms；当θ=70°时t_p=10 ms。

4.1.2 前抛体90°入水时的主弹体载荷

图9为前抛体和主弹体以100 m/s垂直入水时的加速度载荷峰值随时间间隔的变化曲线。图10给出了三个不同τ时的主弹体载荷曲线，可见在不同τ值范围内主弹体载荷具有不同的特性。当τ<0.77时，主弹体直接进入至空泡内部，主弹体不与空泡壁面发生任何接触，与前方受到水阻力而减速的前抛体相撞，出现情况Ⅰ。此时，主弹体受到的冲击加速度远远超过无前抛体时的冲击加速度。图11为τ=0.08时主弹体与前抛体的速度曲线，可以看到主弹体与前抛体相撞之后主弹体减速，前抛体加速，二者再次分离。前抛体受到水阻力快速减速的同时排开周围的水体，使主弹体受到的阻力减小，主弹体的加速度较小，二者的速度差使它们再次碰撞。以上过程多次循环后二者的速度差逐渐减小，最终接触后不再分离，由主弹体推动前抛体一起运动。

图 9 前抛体垂直入水时的主弹体载荷峰值曲线

Figure 9. Load peak curves of the main projectile when the front body enters water vertically

下载: 全尺寸图片幻灯片

图 10 不同τ时前抛体垂直入水时的主弹体载荷

Figure 10. Load of the main projectile when the front body vertically enters the water while τ is different

下载: 全尺寸图片幻灯片

图 11 当τ=0.08时主弹体与前抛体的速度曲线

Figure 11. Velocity curves of main projectile and front body when τ=0.08

下载: 全尺寸图片幻灯片

当0.77≤τ<1时，主弹体冲击加速度曲线中出现两组峰值。主弹体在进入空泡前先撞击水面以上的飞溅水体，主弹体受到第一次冲击，出现情况Ⅲ；之后主弹体与空泡内的前抛体相撞，出现情况Ⅰ。对比两次载荷峰值，可见情况Ⅰ产生的峰值远远大于情况Ⅲ。

当τ≥1时，在主弹体和前抛体的撞击载荷峰值前存在一次较小的载荷，这个载荷由空泡壁面产生。主弹体入水过程中先撞击射流，然后撞击空泡壁面，最后撞击前抛体，此过程中依次出现情况Ⅳ、情况Ⅱ、情况Ⅰ。对比三次载荷，峰值由大到小依次为情况Ⅰ、情况Ⅳ、情况Ⅱ。

前抛体90°入水时，虽然τ不同时主弹体受到的载荷特性不同，但是主弹体和前抛体的碰撞都是主弹体所受冲击载荷的最主要原因。在目前的计算时间内主弹体的入水载荷随时间间隔的增大而减小。

4.1.3 前抛体80°入水时的主弹体载荷

图12为前抛体80°、主弹体垂直，两弹体以100 m/s入水时主弹体的加速度载荷峰值随时间间隔的变化曲线。图13给出了三个不同τ时的主弹体载荷曲线，在不同τ值范围内主弹体载荷具有不同的特性。图14给出了τ=1时主弹体入水过程中与水体的砰击过程。当τ<0.7时，主弹体直接进入空泡撞击空泡壁面，出现情况Ⅱ；当0.7≤τ<1时，冲击加速度曲线出现两组峰值，分别由飞溅水体和空泡壁面产生，出现情况Ⅲ、Ⅱ；当τ≥1时，两组峰值分别由射流和空泡壁面产生，出现情况Ⅳ、Ⅱ。

图 12 前抛体80°入水时的主弹体载荷峰值曲线

Figure 12. Load peak curves of the main projectile when the front body enters the water at 80°

下载: 全尺寸图片幻灯片

图 13 不同τ时前抛体80°入水时的主弹体载荷

Figure 13. Load of main projectile when the front body enters the water at 80° while τ is different

下载: 全尺寸图片幻灯片

图 14 前抛体80°入水，τ=1时的两次砰击过程

Figure 14. Two slamming processes when the front body enters the water at 80 ° and τ=1

下载: 全尺寸图片幻灯片

前抛体80°入水时，对主弹体的载荷由空泡壁面和飞溅水体或射流产生，其中射流产生的载荷最大，其次为空泡壁面，飞溅水体产生的载荷最小，即产生的载荷峰值由大到小依次为情况Ⅳ、情况Ⅱ、情况Ⅲ。在射流生成前主弹体入水载荷随时间间隔的增大而减小。

4.1.4 前抛体70°入水时的主弹体载荷

图15为前抛体以70°、主弹体垂直以100 m/s入水时主弹体冲击加速度载荷随时间间隔的变化曲线。图16给出了两个不同τ时的主弹体载荷曲线，在不同τ值范围内主弹体载荷具有不同的特性。图17给出了τ=0.8时主弹体入水过程中与水体的砰击过程。当τ<0.6时，主弹体不与飞溅水体相撞，只产生一次载荷，为情况Ⅱ；当τ≥0.6时，加速度曲线中出现两组峰值，分别由飞溅水体和空泡壁面产生，出现情况Ⅲ、Ⅱ。

图 15 前抛体70°入水时的主弹体载荷峰值曲线

Figure 15. The load peak curves of the main projectile when the front body enters the water at 70°

下载: 全尺寸图片幻灯片

图 16 不同τ时前抛体70°入水时的主弹体载荷

Figure 16. The load of the main projectile when the front body enters the water at 70° while τ is different

下载: 全尺寸图片幻灯片

图 17 前抛体70°入水，τ=0.8时的两次砰击过程

Figure 17. Two slamming processes when the front body enters the water at 70° and τ=0.8

下载: 全尺寸图片幻灯片

前抛体以70°入水时形成的空泡和飞溅水体的形态都不规律，因此主弹体的载荷随时间的变化也不规律。情况Ⅱ产生的载荷峰值大于情况Ⅲ。

4.1.5 总结

如图18所示，三种不同的前抛体入水角度分别产生了不同的入水情况，每种入水情况下主弹体的入水载荷特性各不相同，最终的降载效果也不同。

图 18 前抛体入水角度变化时主弹体载荷峰值随τ的变化曲线

Figure 18. Curves of the main projectile load peak with τ when the water-entry angle of the front body changes

下载: 全尺寸图片幻灯片

前抛体90°入水产生了情况Ⅰ～Ⅳ。入水情况Ⅰ会产生较大的主弹体载荷，而且该载荷有可能远超无前抛体时的载荷，应尽量避免。载荷随τ的增大而减小。当τ=1.46时主弹体入水载荷最小，为575.8 km/s²，降载率为28.588%，降载效果并不理想。

前抛体80°入水产生了情况Ⅱ～Ⅳ。在形成射流前主弹体入水载荷由撞击空泡壁面（情况Ⅱ）产生，载荷随时间间隔的增大而减小。射流形成后主弹体载荷主要由撞击射流（情况Ⅲ）产生，由于射流与主弹体速度相反，它们之间的相对速度较大，产生了较大的冲击，但是射流的密度较小，冲击产生的主弹体载荷依然小于无前抛体时的载荷。射流形成前主弹体的降载率在80%～90%之间，射流形成后在50%～60%之间，其中τ =0.85时主弹体入水载荷最小，为80.9 km/s²，降载率为89.967%。

前抛体70°入水时产生情况Ⅱ、Ⅲ，这一条件下的空泡、飞溅水体的形态都不规律，主弹体入水载荷以空泡壁面载荷为主，随τ的变化不规律。总体降载率在70%～85%之间，其中τ=0.9时主弹体入水载荷最小，为131.8 km/s²，降载率为83.670%。

从以上分析可知，前抛体以80°入水、τ处于0.8～0.9时降载效果最好。此时主弹体先经过结构松散的飞溅水体经历第一次减速，然后再砰击扩展中的空泡壁面，主弹体和水体自由表面间的相对速度降低，主弹体受到的冲击减小。

4.2 前抛体尺寸变化对主弹体载荷的影响

根据前述分析，选取前抛体以80°、100 m/s入水，τ=0.8作为入水条件，研究直径为2～6 mm，长径比为1的前抛体对主弹体入水载荷的影响。主弹体依然以100 m/s垂直入水。图19给出了不同前抛体尺寸下的主弹体载荷曲线，可见随着前抛体尺寸增大，其入水动能增大，产生的空泡壁面扩展速度更快，主弹体与空泡壁面间的相对速度减小，降低了主弹体的入水载荷。当前抛体直径为6 mm时，主弹体入水冲击的降载率为93.905%。此外，由于空泡扩展速度更快，开口处的Bernoulli效应更加明显，空泡闭合时间t_p也随之减小。

图 19 τ=0.8时前抛体不同尺寸下的主弹体入水载荷曲线

Figure 19. Water-entry load curve of the main projectile under different sizes of the front body while τ=0.8

下载: 全尺寸图片幻灯片

4.3 速度变化对主弹体载荷的影响

为研究在不同入水速度对降载效果的影响，记主弹体和前抛体的初速度为v₀，设置v₀=40～100 m/s，τ=0.8。前抛体直径为6 mm，以80°入水；主弹体垂直入水。前抛体入水速度越大，空泡壁面的扩展速度越快，主弹体和空泡壁面的相对速度越小，如图20所示主弹体载荷随入水初速度的增加而增大。如图21所示，随入水初速度的增加，主弹体的降载率增大，所有计算工况的降载率均在90%以上。

图 20 不同入水初速度下的主弹体载荷

Figure 20. The main projectile load at different initial velocities

下载: 全尺寸图片幻灯片

图 21 主弹体降载率随入水初速度的变化

Figure 21. Load reduction ratio of the main projectile varied with initial water-entry velocity

下载: 全尺寸图片幻灯片

5. 结　论

通过LS-DYNA数值模拟研究了前抛体产生空泡为随后入水的主弹体降载的方法。研究了前抛体入水角度、两弹体入水时间间隔、前抛体尺寸、两弹体入水初速度对主弹体入水载荷的影响，分析了不同入水条件下主弹体的降载规律，得出如下结论：

(1) 前抛体垂直入水时，主弹体会因与前抛体相撞产生更大的冲击载荷，随二者入水时间间隔的增大，载荷有所降低；前抛体倾斜入水时可以避免主弹体与前抛体相撞的情况出现，主弹体撞击飞溅水体后再撞击空泡壁面入水，起到了很好的降载效果，降载率在70%～90%之间；

(2) 当主弹体与前抛体之间的无量纲入水时间间隔为0.8～0.9时，主弹体的降载效果最佳；

(3) 主弹体入水载荷随前抛体尺寸增大而减小，前抛体尺寸越大，降载效果越好；

(4) 有前抛体时的主弹体入水降载率随两弹体入水初速度的增大而增大，前抛体对高速入水的主弹体的降载效果更好。

图 1 含前抛体弹体入水情况示意图

Figure 1. Projectile with front projectile entering the water

下载: 全尺寸图片幻灯片

图 2 S-ALE的界面捕捉法

Figure 2. Surface capturing method of S-ALE

下载: 全尺寸图片幻灯片

图 3 弹体模型与流体域网格

Figure 3. Projectile models and fluid domain mesh

下载: 全尺寸图片幻灯片

图 4 弹体入水运动状态和空泡形态

Figure 4. Projectile entering water motion state and cavitation shape

下载: 全尺寸图片幻灯片

图 5 弹体入水深度(z)

Figure 5. Projectile water entry depth (z)

下载: 全尺寸图片幻灯片

图 6 主弹体的加速度(a)

Figure 6. The accelerationof the main projectile (a)

下载: 全尺寸图片幻灯片

图 7 前抛体初始入水角示意

Figure 7. Schematic of the initial water entry angle of front body

下载: 全尺寸图片幻灯片

图 8 前抛体以不同角度入水形成的空泡演化过程

Figure 8. The evolution process of the cavitation formed by the front body entering the water at different angles

下载: 全尺寸图片幻灯片

图 9 前抛体垂直入水时的主弹体载荷峰值曲线

Figure 9. Load peak curves of the main projectile when the front body enters water vertically

下载: 全尺寸图片幻灯片

图 10 不同τ时前抛体垂直入水时的主弹体载荷

Figure 10. Load of the main projectile when the front body vertically enters the water while τ is different

下载: 全尺寸图片幻灯片

图 11 当τ=0.08时主弹体与前抛体的速度曲线

Figure 11. Velocity curves of main projectile and front body when τ=0.08

下载: 全尺寸图片幻灯片

图 12 前抛体80°入水时的主弹体载荷峰值曲线

Figure 12. Load peak curves of the main projectile when the front body enters the water at 80°

下载: 全尺寸图片幻灯片

图 13 不同τ时前抛体80°入水时的主弹体载荷

Figure 13. Load of main projectile when the front body enters the water at 80° while τ is different

下载: 全尺寸图片幻灯片

图 14 前抛体80°入水，τ=1时的两次砰击过程

Figure 14. Two slamming processes when the front body enters the water at 80 ° and τ=1

下载: 全尺寸图片幻灯片

图 15 前抛体70°入水时的主弹体载荷峰值曲线

Figure 15. The load peak curves of the main projectile when the front body enters the water at 70°

下载: 全尺寸图片幻灯片

图 16 不同τ时前抛体70°入水时的主弹体载荷

Figure 16. The load of the main projectile when the front body enters the water at 70° while τ is different

下载: 全尺寸图片幻灯片

图 17 前抛体70°入水，τ=0.8时的两次砰击过程

Figure 17. Two slamming processes when the front body enters the water at 70° and τ=0.8

下载: 全尺寸图片幻灯片

图 18 前抛体入水角度变化时主弹体载荷峰值随τ的变化曲线

Figure 18. Curves of the main projectile load peak with τ when the water-entry angle of the front body changes

下载: 全尺寸图片幻灯片

图 19 τ=0.8时前抛体不同尺寸下的主弹体入水载荷曲线

Figure 19. Water-entry load curve of the main projectile under different sizes of the front body while τ=0.8

下载: 全尺寸图片幻灯片

图 20 不同入水初速度下的主弹体载荷

Figure 20. The main projectile load at different initial velocities

下载: 全尺寸图片幻灯片

图 21 主弹体降载率随入水初速度的变化

Figure 21. Load reduction ratio of the main projectile varied with initial water-entry velocity

下载: 全尺寸图片幻灯片

表 1 水的Grüneisen状态方程参数

Table 1. Grüneisen equation of state parameters of water

ρ_w/(kg·m⁻³) c_w/(m·s⁻¹) S₁ S₂ S₂ γ₀ a₁ E_w/(J·m⁻³) μ_w0

1000 1480 1.979 0 0 0.11 3 3.07×10⁵ 1

下载: 导出CSV

表 2 空气的多项式状态方程参数

Table 2. Polynomial equation of state parameters of air

ρ_a/(kg·m⁻³) C₀ C₁ C₂ C₃ C₄ C₅ C₆ E_a/(J·m⁻³) μ_a0

1.225 0 0 0 0 0.4 0.4 0 2.5×10⁵ 1

下载: 导出CSV

参考文献(28)

[1]	路龙龙. 空投鱼雷入水技术研究 [D]. 西安: 西北工业大学, 2006. LU L L. Research on theory and application of torpedo’s water entry [D]. Xi’an: Northwestern Polytechnical University, 2006.
[2]	BOTTOMLEY G H. The impact of a model seaplane float on water: No 583 [R]. 1919.
[3]	VON KARMAN T. The impact on seaplane floats during landing: NACA TN321 [R]. Washington, DC: National Advisory Committee on Aeronautics, 1929.
[4]	WAGNER V H. Phenomena associated with impacts and sliding on liquid surfaces [J]. Zeitschrift fur Angewandet Mathematik und Mechanik, 1932, 12(4): 193–215. DOI: 10.1002/zamm.19320120402.
[5]	EROSHIN V A, ROMANENKOV N I, SEREBRYAKOV I V, et al. Hydrodynamic forces produced when blunt bodies strike the surface of a compressible fluid [J]. Fluid Dynamics, 1981, 15(6): 829–835. DOI: 10.1007/BF01096631.
[6]	石汉成, 蒋培, 程锦房. 头部形状对水雷入水载荷及水下弹道影响的数值仿真分析 [J]. 舰船科学技术, 2010, 32(10): 104–107. DOI: 10.3404/j.issn.1672-7649.2010.10.027. SHI H C, JIANG P, CHENG J F. Research on numerical simulation of mine water-entry impact acceleration and underwater ballistic trajectory under the different mine’s head shape [J]. Ship Science and Technology, 2010, 32(10): 104–107. DOI: 10.3404/j.issn.1672-7649.2010.10.027.
[7]	刘华坪, 余飞鹏, 张岳青, 等. 不同头型鱼雷入水冲击载荷研究 [J]. 水下无人系统学报, 2018, 26(6): 527–532. DOI: 10.11993/j.issn.2096-3920.2018.06.003. LIU H P, YU F P, ZHANG Y Q, et al. Analyzing water-entry impact load on torpedo with different head types [J]. Journal of Unmanned Undersea Systems, 2018, 26(6): 527–532. DOI: 10.11993/j.issn.2096-3920.2018.06.003.
[8]	SHI Y, PAN G, YIM S C, et al. Numerical investigation of hydroelastic water-entry impact dynamics of AUVs [J]. Journal of Fluids and Structures, 2019, 91: 102760. DOI: 10.1016/j.jfluidstructs.2019.102760.
[9]	SHARKER S I, HOLEKAMP S, MANSOOR M M, et al. Water entry impact dynamics of diving birds [J]. Bioinspiration & Biomimetics, 2019, 14(5): 056013. DOI: 10.1088/1748-3190/ab38cc.
[10]	HOWARD E A. Protective nose cap for torpedoes: U. S. Patent 2889772 [P]. 1959-06-09.
[11]	宣建明, 宋志平, 严忠汉. 鱼雷入水缓冲保护头帽解体试验研究 [J]. 鱼雷技术, 1999, 7(2): 41–46.
[12]	权晓波, 包健, 孙龙泉, 等. 基于耦合欧拉-拉格朗日算法的航行体缓冲头帽冲击性能 [J]. 兵工学报, 2022, 43(4): 851–860. DOI: 10.12382/bgxb.2021.0168. QUAN X B, BAO J, SUN L Q, et al. Impact performance of cushion nose cap of underwater vehicle based on CEL method [J]. Acta Armamentarii, 2022, 43(4): 851–860. DOI: 10.12382/bgxb.2021.0168.
[13]	孙龙泉, 王都亮, 李志鹏, 等. 基于CEL方法的航行体高速入水泡沫铝缓冲装置降载性能分析 [J]. 振动与冲击, 2021, 40(20): 80–88. DOI: 10.13465/j.cnki.jvs.2021.20.011. SUN L Q, WANG D L, LI Z P, et al. Analysis on load reduction performance of foamed aluminum buffer device for high speed water entry of vehicle based on a CEL method [J]. Journal of Vibration and Shock, 2021, 40(20): 80–88. DOI: 10.13465/j.cnki.jvs.2021.20.011.
[14]	魏海鹏, 史崇镔, 孙铁志, 等. 基于ALE方法的航行体高速入水缓冲降载性能数值研究 [J]. 爆炸与冲击, 2021, 41(10): 115–126. DOI: 10.11883/bzycj-2020-0461. WEI H P, SHI C B, SUN T Z, et al. Numerical study on load-shedding performance of a high-speed water-entry vehicle based on an ALE method [J]. Explosion and Shock Waves, 2021, 41(10): 115–126. DOI: 10.11883/bzycj-2020-0461.
[15]	LI Y, ZONG Z, SUN T Z. Crushing behavior and load-reducing performance of a composite structural buffer during water entry at high vertical velocity [J]. Composite Structures, 2021, 255: 112883. DOI: 10.1016/j.compstruct.2020.112883.
[16]	李向阳. 反潜导弹缓冲头帽入水冲击特性的数值研究 [D]. 哈尔滨: 哈尔滨工程大学, 2021. DOI: 10.27060/d.cnki.ghbcu.2021.000904. LI X Y. Numerical study on the impact characteristics of anti-submarine missile's cushion nose cap into water [D]. Harbin: Harbin Engineering University, 2021. DOI: 10.27060/d.cnki.ghbcu.2021.000904.
[17]	刘登科. 大尺度航行体高速入水缓冲降载特性研究 [D]. 哈尔滨: 哈尔滨工程大学, 2020. DOI: 10.27060/d.cnki.ghbcu.2020.002187. LIU D K. Research on load-reducing characteristics of water-entry vehicle with large scale and high speed [D]. Harbin: Harbin Engineering University, 2020. DOI: 10.27060/d.cnki.ghbcu.2020.002187.
[18]	SHI Y, GAO X F, PAN G. Design and load reduction performance analysis of mitigator of AUV during high speed water entry [J]. Ocean Engineering, 2019, 181: 314–329. DOI: 10.1016/j.oceaneng.2019.03.062.
[19]	潘龙, 王焕然, 姚尔人, 等. 头部喷气平头圆柱体人水缓冲机制研究 [J]. 工程热物理学报, 2015, 36(8): 1691–1695. PAN L, WANG H R, YAO E R, et al. Mechanism research on the water-entry impact of the head-jetting flat cylinder [J]. Journal of Engineering Thermophysics, 2015, 36(8): 1691–1695.
[20]	刘华坪, 余飞鹏, 韩冰, 等. 头部喷气影响航行体入水载荷的数值模拟 [J]. 工程热物理学报, 2019, 40(2): 300–305. LIU H P, YU F P, HAN B, et al. Numerical simulation study on influence of top jet in object water entering impact [J]. Journal of Engineering Thermophysics, 2019, 40(2): 300–305.
[21]	赵海瑞, 施瑶, 潘光. 头部喷气航行器高速入水空泡特性数值分析 [J]. 西北工业大学学报, 2021, 39(4): 810–817. DOI: 10.1051/jnwpu/20213940810. ZHAO H R, SHI Y, PAN G. Numerical simulation of cavitation characteristics in high speed water entry of head-jetting underwater vehicle [J]. Journal of Northwestern Polytechnical University, 2021, 39(4): 810–817. DOI: 10.1051/jnwpu/20213940810.
[22]	RABBI R, SPEIRS N B, KIYAMA A, et al. Impact force reduction by consecutive water entry of spheres [J]. Journal of Fluid Mechanics, 2021, 915: A55. DOI: 10.1017/jfm.2020.1165.
[23]	LYU X J, YUN H L, WEI Z Y. Influence of time interval on the water entry of two spheres in tandem configuration [J]. Experiments in Fluids, 2021, 62(11): 222. DOI: 10.1007/s00348-021-03300-w.
[24]	辛春亮, 涂建, 王俊林, 等. 由浅入深精通LS-DYNA [M]. 北京: 中国水利水电出版社, 2019.
[25]	桂蜀旺. 高速弹体入水冲击过程流固耦合分析 [D]. 南京: 南京航空航天大学, 2020. DOI: 10.27239/d.cnki.gnhhu.2020.000018. GUI S W. Analysis of high speed projectile water-entry process based on fluid-solid coupling [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2020. DOI: 10.27239/d.cnki.gnhhu.2020.000018.
[26]	黄兆铭. 单立柱三桩式海上风机受冰风荷载作用的损伤研究 [D]. 哈尔滨: 哈尔滨工业大学, 2021. DOI: 10.27061/d.cnki.ghgdu.2021.004343. HUANG Z M. Research on the damage of single-column three-piles offshore wind turbines subjected to ice and wind loads [D]. Harbin: Harbin Institute of Technology, 2021. DOI: 10.27061/d.cnki.ghgdu.2021.004343.
[27]	王明振, 曹东风, 吴彬, 等. 基于S-ALE流固耦合方法的飞机水上迫降动力学数值分析 [J]. 重庆大学学报, 2020, 43(6): 21–29. DOI: 10.11835/j.issn.1000-582X.2020.06.003. WANG M Z, CAO D F, WU B, et al. Numerical analysis of aircraft dynamic behavior in ditching based on S-ALE fluid-structure interaction method [J]. Journal of Chongqing University, 2020, 43(6): 21–29. DOI: 10.11835/j.issn.1000-582X.2020.06.003.
[28]	CHEN T, HUANG W, ZHANG W, et al. Experimental investigation on trajectory stability of high-speed water entry projectiles [J]. Ocean Engineering, 2019, 175: 16–24. DOI: 10.1016/j.oceaneng.2019.02.021.

施引文献

期刊类型引用(3)

1.	职明洋，燕国军，孙龙泉，王鹏霄. 带气囊结构航行体入水回收动力学特性研究. 力学学报. 2024(04): 943-959 . 百度学术
2.	郑伟，李强，范旭东，吕续舰. 跨介质航行器高速入水降载方法研究综述. 水下无人系统学报. 2024(03): 411-425 . 百度学术
3.	迟铁. 流固耦合分析下的船体高速入水冲击数值模拟. 舰船科学技术. 2023(19): 60-63 . 百度学术