航行体开槽包裹式缓冲头帽结构设计及其降载性能

施瑶; 刘振鹏; 潘光; 高兴甫

doi:10.11883/bzycj-2021-0426

航行体开槽包裹式缓冲头帽结构设计及其降载性能

doi: 10.11883/bzycj-2021-0426

施瑶^{1, 2,},
刘振鹏^{1, 2},
潘光^{1, 2},
高兴甫^{1, 2}

1.
西北工业大学航海学院，陕西西安 710072
2.
西北工业大学航海学院无人水下运载技术工业和信息化部重点实验室，陕西西安 710072

基金项目: 国家自然科学基金（U21B2055，52171324）；中央高校基本科研业务费（3102019JC006）

详细信息

作者简介:
施　瑶（1988－　），男，博士，副研究员，shiyao@nwpu.edu.cn

中图分类号: O352
计量
- 文章访问数: 469
- HTML全文浏览量: 171
- PDF下载量: 60
- 被引次数: 5
出版历程
- 收稿日期: 2021-10-18
- 修回日期: 2022-08-10
- 网络出版日期: 2022-10-26
- 刊出日期: 2022-12-08

Structural design of a slotted wrapping buffer head cap of vehicles and its load reduction performance

SHI Yao^{1, 2
,},
LIU Zhenpeng^{1, 2},
PAN Guang^{1, 2},
GAO Xingfu^{1, 2}

1.
School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an 710072, Shaanxi, China
2.
Key Laboratory of Unmanned Underwater Vehicle, Ministry of Industry and Information Technology, School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an 710072, Shaanxi, China

摘要

摘要: 针对空投航行体和火箭助飞航行体高速入水过程中遭受巨大的冲击载荷，可能导致的结构损坏、弹道失控等问题，提出了一种开槽包裹式缓冲头帽，用于保护航行体入水过程中的结构安全。首先，给出了缓冲头帽的详细设计参数，基于任意拉格朗日-欧拉算法，建立了航行体带缓冲头帽高速入水数值模型，并对该数值模型的正确性进行了验证。然后，在此基础上，研究了不同入水角度下，空泡流场的演变过程，分析了入水时缓冲材料的应力分布情况。最后，探究了不同入水速度和入水角度下缓冲头帽的降载性能。结果表明，数值计算所得空泡形态与实验图像基本吻合，且数值计算和实验测试所得的冲击加速度变化趋势基本一致，两者轴向加速度峰值相对误差为6.72%，径向加速度峰值相对误差为7.52%。航行体装备所设计的缓冲头帽以300 m/s的速度垂直入水时轴向降载率为22.17%；以100 m/s的速度60°入水时，轴向降载率为31.83%，径向降载率为66.80%。
- 航行体 /
- 高速入水 /
- 缓冲头帽 /
- 降载性能 /
- 任意拉格朗日-欧拉算法
Abstract: In order to solve the problems of structural damage and ballistic runaway caused by huge impact loads suffered by air-drop vehicles and rocket-assisted vehicles during high-speed water-entry, a slotted wrapping buffer head cap was proposed to guarantee the structural safety of the vehicles during water entry. Firstly, the structural composition and detailed parameters of the head cap were given, and a numerical model for the high-speed water-entry of the vehicles was established based on the arbitrary Lagrangian-Eulerian (ALE) algorithm. The Lagrangian viewpoint was used to solve the small deformation of the vehicle and the head cap, and the Eulerian viewpoint was used to capture the large deformation of the free surface such as water and air, thereby overcoming the problems that the Eulerian mesh was not accurate enough to solve the structural deformation and the numerical oscillation caused by mesh distortion in solving large deformation problems by the Lagrangian mesh. On this basis, the evolution processes of the cavity and flow field around the vehicle entering the water with the head cap at different angles were studied by numerical simulation, and the interaction process between the head cap and the water was given. Furthermore, the distribution of effective stress of the buffer was analyzed when it entering the water vertically and obliquely. Finally, the load reduction performances of the head cap when the vehicle entered the water at different velocities and angles were investigated. The results show that the cavities obtained by the simulation are basically consistent with the experimental images, and the change trends of impact acceleration are basically consistent with the experimental results. The relative error of the axial peak acceleration between the numerical simulation and experiment is 6.72%, and the relative error of the radial peak acceleration is 7.52%. The ratio of axial load reduction is 22.17% when the vehicle enters the water vertically with a head cap at 300 m/s. At the same time, the ratio of axial load reduction is 31.83% and the ratio of radial load reduction is 66.80% when the vehicle with a head cap enters the water at 100 m/s and 60°. So this research has a certain guiding role in the design of new load-reduction structure.
- vehicle /
- high-speed water-entry /
- buffer head cap /
- load reduction performance /
- arbitrary Lagrangian-Eulerian algorithm

HTML全文

随着军事科技的发展，各种新型的军事武器装备不断出现，其中采用固定翼飞机、直升飞机投放或火箭助飞等新型发射方式的跨水空介质航行体具有防区外发射、突防能力强及隐蔽性高等优点，成为了军事武器装备发展的热点。攻击过程中，需要由空中弹道转变为水中弹道，该类航行体将经历一个跨介质高速入水的过程。在该过程中，由于介质密度的突变以及复杂的湍流流动等原因，航行体将承受巨大的冲击载荷，引起壳体结构的动态响应，导致壳体结构发生弹塑性变形，甚至造成断裂、破损和屈曲；或对航行体的弹道轨迹产生影响，致使航行体入水过程中发生忽扑、跳水和弹道失控等严重问题。因此，减小航行体高速入水时受到的冲击载荷，确保安全入水，是至关重要的。

比较常见的降载方法有：在航行体尾部挂载降落伞，降低入水速度并稳定空投弹道^[1-2]；利用航行体的头部喷出气体，缓和入水冲击环境^[3-5]；在航行体头部安装一个气囊^[6]或者缓冲头帽，利用材料的变形和破碎来吸收航行体高速入水过程中承受的巨大冲击能量。严忠汉^[7]讨论了入水弹头缓冲材料的动态特性，提出了切片法、临界设计状态等概念，这些概念对后续缓冲头帽的设计具有重要意义。王永虎等^[8]引进了聚氨酯泡沫塑料在冲击加载条件下的本构关系式，代替了传统的三段式经验模式，分析了航行体垂直入水情况下缓冲器的动态缓冲特性。早期，对于缓冲头帽的研究大多基于实验研究，或者对单一的缓冲材料特性进行分析。实验研究存在耗资昂贵、可重复差等问题，对单一缓冲材料的特性分析反映到实际的入水过程中存在较大的局限性，因此如今对于缓冲头帽的研究大多基于数值计算手段。其中比较常用的算法为任意的拉格朗日-欧拉（arbitrary Lagrangian-Eulerian, ALE）算法^[9]，该算法为解决作用过程复杂的水气固多介质耦合问题提供了新的求解思路。Wang等^[10]研究了轻型金字塔夹层板结构低速入水时的流固耦合动态响应，获得了冲击压力和结构变形的特性。李建阳等^[11]通过对大型返回舱入水时不同工况的数值计算，分析了入水速度和入水角度对冲击载荷特性的影响。胡明勇等^[12]建立了低亚声速射弹垂直入水的流固耦合数值计算模型，对入水空泡、流场与弹道间的多介质耦合问题进行了数值求解。颜彬等^[13]在考虑弹体内部支撑结构的情况下，对弹体入水时结构的力学响应进行了数值计算。Wu等^[14]对装备缓冲头帽的射弹的入水过程进行了数值计算，提出了计算射弹入水过程中减速度的经验公式。Li等^[15]对装备缓冲头帽的航行体的入水过程进行了数值计算，重点研究了入水过程中罩壳的变形失效形式。魏海鹏等^[16]对适用于直径为324 mm的航行体的缓冲组件进行了数值研究，发现在撞水时罩壳的头部和预置沟槽处会出现明显的应力集中现象，分层的缓冲泡沫出现了二次缓冲等现象。

对于航行体入水缓冲头帽的研究不乏公开报道，但是现有的研究主要是针对于直径为324 mm的小尺度航行体入水情形，适用于大尺度航行体高速安全入水的缓冲头帽结构设计及其降载性能分析研究较少。直径为533 mm的大尺度航行体在高速入水时，所面临的冲击环境更为极端，流体对航行体壳体结构的作用更为剧烈，小尺度航行体降载结构的相关设计理论难以直接应用于大尺度航行体。

本文中，利用LS-DYNA软件，并基于ALE算法，建立装备缓冲头帽的航行体高速入水数值模型，对适用于大尺度航行体高速入水的开槽包裹式缓冲头帽的缓冲性能进行数值研究，探究其缓冲性能的相关影响因素，以期为新型缓冲头帽的结构设计提供指导。

1. 缓冲头帽

根据鱼雷缓冲头帽设计可行域^[17]，本文中所设计的缓冲头帽主要由外部开槽罩壳和内部缓冲材料组成。罩壳可以承受航行体投放时的气动压力，保证缓冲头帽在入水之前整体结构的完整性，同时使装备缓冲头帽的航行体在空中飞行时阻力较小，提高航行体的飞行速度，增大航行体的投送距离；缓冲材料在入水过程中依靠自身变形吸收巨大的冲击能量，从而降低作用在航行体上的冲击载荷。

1.1 罩　壳

航行体的外形参考鱼雷MK48，最大直径为533 mm，总体长度为5 850 mm，如图1所示。为了保证装备缓冲头帽后的航行体具有良好的气动外形，本文中所设计的罩壳采用尖拱体外形，罩壳圆柱段外径为555 mm，厚度为5 mm，尖拱段长度为600 mm，圆柱段长度为700 mm，如图2所示。为了减小缓冲头帽对航行体入水弹道的影响，使它更容易与航行体分离，沿罩壳的周向开8个强度减弱槽，轴向开槽长度为900 mm，开槽深度为3 mm，轴向开槽角度为15°，如图3所示。

图 1 航行体外形

Figure 1. Shape of the vehicle

下载: 全尺寸图片幻灯片

图 2 罩壳外形

Figure 2. Shape of the nose cap

下载: 全尺寸图片幻灯片

图 3 罩壳开槽示意图

Figure 3. Schematic diagrams of the slotted nose cap

下载: 全尺寸图片幻灯片

1.2 缓冲件

本文中设计了一种开槽包裹式缓冲件，如图4所示，缓冲件内部线型与航行体头部线型一致，能够完全包裹航行体头部曲线段，二者接触面积大，在航行体入水角度较小时，可以有效地保护航行体头部壳体设备。同时，为了避免缓冲件被压实后不容易破碎的问题，沿缓冲件周向均匀开8个槽，开槽底部距离缓冲件前端100 mm，缓冲件开槽前端深度为50 mm，缓冲件包裹段后端厚度为6 mm，开槽深度为3 mm。罩壳主要发生撑进破坏^[17]，适当增大罩壳前段厚度有利于它破碎，因此缓冲件前端向后缩进一部分，罩壳前段加厚。航行体装备缓冲头帽后的整体示意图如图5所示。

图 4 开槽包裹式缓冲件

Figure 4. A slotted wrapping buffer

下载: 全尺寸图片幻灯片

图 5 航行体装配缓冲头帽后整体

Figure 5. The whole body of the vehicle assembled with the buffer head cap

下载: 全尺寸图片幻灯片

2. 计算模型

有限元算法中比较常见的算法有拉格朗日算法与欧拉算法以及任意的拉格朗日-欧拉算法（ALE算法）。拉格朗日算法的网格会随着物质进行运动，对结构边界的捕捉比较精确，但是当结构发生大变形时，可能导致网格出现畸变的现象，引起计算结果不准确。欧拉算法的网格是固定不动的，物质在网格之间进行输运，对网格的变形没有限制，但对于结构边界的捕捉精度不够。对于航行体高速入水的问题，选择ALE算法进行求解，发挥二者的优势，利用拉格朗日网格求解航行体及缓冲头帽的变形破碎，捕捉精度高，利用欧拉网格来求解自由液面的大变形，数值计算稳定。

2.1 ALE算法基本原理

ALE算法的控制方程由质量守恒方程、动量守恒方程和能量守恒方程组成^[18]。

(1)质量守恒方程：

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

(1)

(2)动量守恒方程：

$v\frac{{\partial {v_i}}}{{\partial t}} = {\sigma _{ij,j}} + \rho {b_i} - \rho {w_i}\frac{{\partial {v_i}}}{{\partial {x_j}}}$

(2)

(3)能量守恒方程：

$\rho \frac{{\partial E}}{{\partial t}} = {\sigma _{ij}}{v_{i,j}} + \rho {b_i}{v_i} - \rho {w_j}\frac{{\partial E}}{{\partial {x_j}}}$

(3)

式中：ρ为流体的密度；x为欧拉坐标；v为流体的速度；u为网格的速度；w=v−u；E为流体的内能；σ_ij为应力张量的分量；b为流体所受的体积力载荷。

ALE算法的基本实现过程如下^[19]：

(1)先进行拉格朗日步计算，单元网格随着材料运动而变形，保持变形后的物体边界条件，对内部重新划分网格，网格的拓扑关系保持不变，称为光滑步。

(2)将变形网格中的单元变量和节点速度矢量输运到新网格中，称为对流步。一般每个单元解的各种变量都要进行输运，要输运的数量取决于材料模型。

ALE算法中流固耦合力计算方法如图6所示。对于结构节点，在每一时间步通过相对速度(v_s−v_f)更新侵彻深度d，结构节点定义为从节点，流体节点定义为主节点，使用流体单元的等参坐标计算主节点的位置。当t=0时，假定d⁽⁰⁾=0，流体节点和结构节点重合；当t⁽ⁿ⁾=t⁽ⁿ⁻¹⁾+Δt，侵彻向量通过下式更新：

图 6 流固耦合算法

Figure 6. Algorithm of fluid-structure interaction

下载: 全尺寸图片幻灯片

${{\boldsymbol{d}}^{\left( n \right)}} = {{\boldsymbol{d}}^{(n - 1)}} + \left( {{{\boldsymbol{v}}_{\text{s}}} - {{\boldsymbol{v}}_{\text{f}}}} \right)\Delta t$

(4)

耦合力 ${\boldsymbol F}$ 采用下式计算：

${\boldsymbol{F}} = k{\boldsymbol{d}}$

(5)

式中：k为基于主从节点质量模型特性的刚度系数。

2.2 有限元模型

选定一个局部坐标系(Oxyz)和一个全局坐标系(O₀x₀y₀z₀)^[20]，如图7所示。其中，全局坐标系原点O₀处于航行体入水点；O₀x₀轴位于水平面内，此轴的方向指向航行体的前进方向，称为参考航向；O₀y₀轴垂直于地面指向上方，称为铅垂轴；O₀z₀轴处于水平面内，垂直于x₀O₀y₀平面，方向由右手坐标法则确定。局部坐标系原点O位于航行体重心处，Ox轴处于航行体对称面内，与航行体几何对称轴一致，指向航行体头部，称为航行体纵轴；Oy轴位于航行体纵对称面内，垂直于Ox轴，指向上方，称为航行体立轴；Oz轴垂直于xOy平面，从航行体尾部观察，Oz轴指向右方，称为航行体横轴。定义Ox轴与水平面的夹角为入水角θ。

图 7 坐标系的定义

Figure 7. Definition of coordinate systems

下载: 全尺寸图片幻灯片

考虑到航行体壳体以及罩壳在入水过程中的弹塑性变形，对它们采用*MAT_PLASTIC_KINEMATIC材料模型。航行体壳体材料选用铝合金，该铝合金的密度为2700 kg/m³，弹性模量为75 GPa，泊松比为0.33，屈服应力为275 MPa，剪切模量为1.33 GPa。罩壳材料选用不饱和聚酯树脂复合材料，该材料的密度为1160 kg/m³，弹性模量为3.5 GPa，泊松比为0.34，屈服应力为101 MPa。对缓冲件采用*MAT_PLASITIC_KINEMATIC材料模型，并设置失效应变，用于模拟缓冲件的破碎过程。当某个单元的应变超过失效应变时，该单元即失效，不再参与计算。不同密度的发泡材料的缓冲性能不同^[21]，所采用缓冲件的材料为硬质聚氨酯泡沫，该泡沫的密度为90 kg/m³，弹性模量为129 MPa，泊松比为0.024，失效应变为0.5，阻尼系数为0.5。

对水和空气等流体采用*MAT_NULL材料模型，需要采用*EOS_GRÜNEISEN状态方程^[22]描述水的运动：

${p_{\text{w}}} = \dfrac{{{\rho _{\text{w}}}{c^2}{\mu _{\text{w}}}\left( {1 + \left( {1 - \dfrac{{{\gamma _0}}}{2}} \right){\mu _{\text{w}}} - \dfrac{\alpha }{2}\mu _{\text{w}}^2} \right)}}{{{{\left( {1 - ({S_1} - 1){\mu _{\text{w}}} - {S_2}\dfrac{{\mu _{\text{w}}^2}}{{{\mu _{\text{w}}} + 1}} - {S_3}\dfrac{{\mu _{\text{w}}^3}}{{{{({\mu _{\text{w}}} + 1)}^3}}}} \right)}^2}}} + ({\gamma _0} + \alpha {\mu _{\text{w}}}){e_{\text{w}}}{\text{，}}\;\;\;{\mu _{\text{w}}} = \dfrac{1}{{{\psi _{\text{w}}}}} - 1$

(6)

式中：p_w为水的压力，c为水中声速，γ₀为Grüneisen常数，α为γ₀的一阶体积修正量，e_w为水的体积内能；S₁、S₂、S₃为斜率系数，ψ_w为水的相对体积。本文中：c=1480 m/s, S₁=2.56, S₂=−1.986, S₃=0.226, S₃=0.226；水的初始体积内能e_w0=0，水的初始相对体积ψ_w0=1。

采用*EOS_LINEAR_POLYNOMIAL状态方程^[23]描述空气的运动：

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

(7)

式中：C₁、C₂、C₃、C₄、C₅、C₆为多项式系数；e_a为空气的体积内能； ${\mu _{\text{a}}} = 1/{\psi _{\text{a}}} - 1$ ，ψ_a为空气的相对体积。本文中：C₀=C₁=C₂=C₃=0, C₄=C₅=0.4, C₆=0；空气的初始体积内能e_a0=250 kJ/m³。

2.3 网格无关性验证

流体网格尺寸与结构网格尺寸之比：

$N = {{{L_{\text{E}}}} \mathord{\left/ {\vphantom {{{L_{\text{E}}}} {{L_{\text{L}}}}}} \right. } {{L_{\text{L}}}}}$

(8)

式中：L_E为流体单元网格的特征尺寸；L_L为结构单元网格的特征尺寸。

航行体的入水冲击加速度用无量纲的系数c_d表示：

${c_{\text{d}}} = \frac{{2ma}}{{{\rho _{\text{w}}}{\text{π }}{R^2}{v^2}}}$

(9)

式中：m为航行体的质量，a为航行体的入水加速度，R为航行体的最大半径，v为航行体的入水速度。

航行体网格尺寸为10 mm， $N$ 取值为1～5；水和空气单元的网格尺寸一致，分别为10、20、30、40和50 mm。航行体的入水速度为100 m/s，入水角度为90°。不同网格比例下的冲击加速度系数时程曲线、加速度峰值系数c_d,p如图8所示，可以看出，随着流-固网格比例的减小，加速度峰值和脉宽逐渐增大，当N=4时计算结果就满足了收敛性要求。综合考虑了计算精度的要求及所需要的计算时间，航行体的网格尺寸设置为10 mm，水和空气的网格尺寸均设置为40 mm，同时，为减小壁面效应对计算结果的影响，空气域和水域的尺寸均设置为10 m×6.4 m×6 m，水域底部设置为固定边界，其余边界条件设置为无反射边界条件，总体划分网格数为8 763 225，如图9所示。

图 8 不同网格尺寸下的加速度系数及其峰值

Figure 8. Time history curves of acceleration coefficient and its peaks under different mesh sizes

下载: 全尺寸图片幻灯片

图 9 计算域

Figure 9. Computational domain

下载: 全尺寸图片幻灯片

2.4 数值计算结果验证

为了验证数值计算模型的正确性与准确性，开展了航行体缩比模型头段带缓冲头帽高速入水实验。实验现场布局如图10所示，采用高压空气作为模型的动力来源；在模型内部安装测试装置，记录模型的加速度信息；利用高速摄像机捕捉模型的位置姿态及入水空泡演变信息；模型的入水速度为101 m/s，入水角度为60°。

图 10 实验现场布局

Figure 10. Experimental layout

下载: 全尺寸图片幻灯片

利用图像处理技术，将数值计算的空泡演变宏观轮廓进行提取，并与实验中高速摄像机拍摄所得的航行体入水空泡轮廓进行叠加对比，如图11所示。从两者的空泡轮廓和演变规律来看，本文中所建立的装备缓冲头帽的航行体高速入水数值计算模型，能较真实地反映实际情况中航行体高速入水的空泡演变宏观特性。

图 11 数值计算与实验空泡的对比

Figure 11. Comparison of cavities between simulation and experiment

下载: 全尺寸图片幻灯片

图12为航行体的入水冲击加速度数值计算结果和实验测试结果的对比。可以看出，数值计算和实验测得的加速度的变化趋势基本一致，两者轴向加速度峰值相对误差为6.72%，径向加速度峰值相对误差为7.52%，两者冲击加速度峰值和峰值脉宽均吻合良好，满足实际工程中数值计算的精度需求，能够支撑后文的数值计算。

图 12 数值计算与实验测试加速度的对比

Figure 12. Comparison of accelerations between simulation and experiment

下载: 全尺寸图片幻灯片

3. 结果与分析

3.1 入水流场演化特性和头帽破碎过程

图13～14为航行体以100 m/s的速度垂直入水和以60°斜入水时的空泡演化图像和缓冲材料破坏图像，为了观察缓冲材料的变形破碎，将罩壳隐藏。初始时刻时，整体结构位于水面上方一定的距离，5 ms时缓冲头帽已经接触到水面，水体被排开，空泡轮廓扩张，缓冲头帽破裂后刺破了空泡壁面，使空泡轮廓变得粗糙。航行体垂直入水时，空泡轮廓基本对称，水面上方液面隆起，有液体飞溅。当航行体以60°斜入水时，由于迎水面与背水面的速度不一致（见图15），导致了两侧的水压不同，压力低的一面空泡扩张速度快，压力高的一面空泡扩张速度慢。因此，航行体斜入水时产生的空泡是非对称的，同时，航行体除了受到一个较大的轴向载荷外，其径向也将承受一个不可忽视的冲击载荷。装备缓冲头帽的航行体斜入水时，缓冲头帽的破碎不对称，相较于垂直入水时，缓冲头帽破碎的范围更大，从图13～14可以看出，不管是垂直入水还是斜入水，缓冲头帽都能发生较大程度的破损，并在入水后期与航行体完全脱离，从而不影响航行体后续的入水弹道。

图 13 航行体以100 m/s、90°入水时流场演化和缓冲件的破坏过程（隐藏罩壳）

Figure 13. Flow field evolution and failure process of buffer when the vehicle enters water at 100 m/s and 90° (hide the nose cap)

下载: 全尺寸图片幻灯片

图 14 航行体以100 m/s、60°入水时流场演化和缓冲件的破坏过程（隐藏罩壳）

Figure 14. Flow field evolution and failure process of buffer when the vehicle enters water at 100 m/s and 60° (hide the nose cap)

下载: 全尺寸图片幻灯片

图 15 航行体以100 m/s的速度在不同入水角度下20 ms时水体的速度矢量

Figure 15. Vectors of velocity of the water when the vehicle enters water at 100 m/s and different angles at 20 ms

下载: 全尺寸图片幻灯片

3.2 缓冲材料的破坏特性

图16～17分别为航行体以100 m/s的速度垂直入水和60°入水角斜入水过程中缓冲材料在不同时刻的等效应力云图。

图 16 航行体以100 m/s、90°入水角入水时缓冲材料的等效应力分布

Figure 16. Distribution of effective stress of the buffer when the vehicle enters water at 100 m/s and 90°

下载: 全尺寸图片幻灯片

图 17 航行体以100 m/s、60°入水角入水时缓冲材料的等效应力分布

Figure 17. Distribution of effective stress of the buffer when the vehicle enters water at 100 m/s and 60°

下载: 全尺寸图片幻灯片

可以明显看出，在缓冲头帽撞水时，缓冲材料内部出现了明显的应力集中现象；应力首先出现在缓冲材料的前端，并以应力波的形式向缓冲材料的后端传播；在航行体垂直入水初期，应力在缓冲材料内部呈现较均匀的分层现象。而在航行体以60°斜入水时，由于缓冲材料两边触碰到水的时刻不同，下部先受到水的挤压作用，先出现应力集中现象，因而导致入水初期应力在缓冲材料内部分布不均匀。从图中也可以看出，7 ms左右缓冲材料被航行体击穿，此时航行体头部直接接触到水，航行体受到的冲击载荷到达峰值。

3.3 入水速度和入水角度对降载性能的影响

探究了不同入水速度和不同入水角度下缓冲头帽的降载性能。当入水角度为90°时，入水速度分别为100、150、200、250和300 m/s；当入水速度为100 m/s时，入水角度分别为60°、70°和80°。

图18(a)为装备缓冲头帽的航行体在不同入水速度下垂直入水时的加速度时程曲线，图18(b)为未装备缓冲头帽的航行体在不同入水速度下垂直入水时的加速度时程曲线。从图18可以看出，随着航行体入水速度的不断升高，冲击加速度的峰值不断增大，脉宽不断增大，缓冲材料被击穿的时刻不断提前，缓冲头帽的缓冲性能减弱，入水速度为300 m/s时，航行体加速度峰值减小22.17%。

图 18 不同的航行体以不同入水速度垂直入水时的加速度时程曲线

Figure 18. Time-history curves of acceleration when different vehicles enter water vertically at different velocities

下载: 全尺寸图片幻灯片

图19～20为航行体以100 m/s的速度在不同入水角度下的加速度时程曲线。可以看出，在入水速度相同时，随着入水角的增大，轴向加速度峰值逐渐增大，径向加速度峰值逐渐减小，缓冲头帽对航行体轴向和径向加速度峰值的减小率越大，当入水角度为60°时，对轴向加速度峰值的减小率为31.83%，对径向加速度峰值的减小率为66.80%。

图 19 不同的航行体以100 m/s的速度在不同入水角度下的轴向加速度时程曲线

Figure 19. Time-history curves of axial acceleration when different vehicles enter water at 100 m/s and different angles

下载: 全尺寸图片幻灯片

图 20 不同的航行体以100 m/s的速度在不同入水角度下的径向加速度时程曲线

Figure 20. Time-history curves of radial acceleration when different vehicles enter water at 100 m/s and different angles

下载: 全尺寸图片幻灯片

4. 结　论

设计了适用于大尺度航行体高速安全入水的尖拱体外形开槽包裹式缓冲头帽，并给出了详细的设计参数。基于ALE算法，对装备开槽包裹式缓冲头帽的航行体高速入水过程进行了数值模拟，并用航行体头段缩比模型高速入水实验验证了所建立数值模型的正确性与准确性。从不同角度，对所设计的开槽包裹式缓冲头帽缓冲性能进行了分析，主要得到以下结论。

(1)装备缓冲头帽的航行体在高速垂直入水时，受缓冲头帽破碎形态的影响，产生了壁面较粗糙的对称空泡；而斜入水时，受两侧水流速度不一致的影响，产生了非对称空泡，对航行体的受力分布产生了不可忽视的影响。

(2)航行体垂直入水时，应力集中区域首先在缓冲材料头部出现，并以应力波的形式向后快速传递；垂直入水时应力分布呈现较均匀的分层现象，斜入水时应力分布也呈现分层现象，但不均匀；应力区域的分布可为缓冲材料的结构设计提供指导。

(3)随着航行体入水速度的提高，冲击加速度的峰值也随之增大，缓冲头帽被击穿的时刻提前，降载性能逐渐下降，以300 m/s的速度垂直入水时，冲击加速度峰值减小22.17%；当航行体斜入水时，入水角度越大，缓冲头帽的降载效果越好，入水速度为100 m/s、入水角度为60°时，轴向冲击加速度峰值减小31.83%，径向冲击加速度峰值减小66.80%，缓冲头帽的降载隔冲作用给航行体高速安全入水提供了保证。

图 1 航行体外形

Figure 1. Shape of the vehicle

下载: 全尺寸图片幻灯片

图 2 罩壳外形

Figure 2. Shape of the nose cap

下载: 全尺寸图片幻灯片

图 3 罩壳开槽示意图

Figure 3. Schematic diagrams of the slotted nose cap

下载: 全尺寸图片幻灯片

图 4 开槽包裹式缓冲件

Figure 4. A slotted wrapping buffer

下载: 全尺寸图片幻灯片

图 5 航行体装配缓冲头帽后整体

Figure 5. The whole body of the vehicle assembled with the buffer head cap

下载: 全尺寸图片幻灯片

图 6 流固耦合算法

Figure 6. Algorithm of fluid-structure interaction

下载: 全尺寸图片幻灯片

图 7 坐标系的定义

Figure 7. Definition of coordinate systems

下载: 全尺寸图片幻灯片

图 8 不同网格尺寸下的加速度系数及其峰值

Figure 8. Time history curves of acceleration coefficient and its peaks under different mesh sizes

下载: 全尺寸图片幻灯片

图 9 计算域

Figure 9. Computational domain

下载: 全尺寸图片幻灯片

图 10 实验现场布局

Figure 10. Experimental layout

下载: 全尺寸图片幻灯片

图 11 数值计算与实验空泡的对比

Figure 11. Comparison of cavities between simulation and experiment

下载: 全尺寸图片幻灯片

图 12 数值计算与实验测试加速度的对比

Figure 12. Comparison of accelerations between simulation and experiment

下载: 全尺寸图片幻灯片

图 13 航行体以100 m/s、90°入水时流场演化和缓冲件的破坏过程（隐藏罩壳）

Figure 13. Flow field evolution and failure process of buffer when the vehicle enters water at 100 m/s and 90° (hide the nose cap)

下载: 全尺寸图片幻灯片

图 14 航行体以100 m/s、60°入水时流场演化和缓冲件的破坏过程（隐藏罩壳）

Figure 14. Flow field evolution and failure process of buffer when the vehicle enters water at 100 m/s and 60° (hide the nose cap)

下载: 全尺寸图片幻灯片

图 15 航行体以100 m/s的速度在不同入水角度下20 ms时水体的速度矢量

Figure 15. Vectors of velocity of the water when the vehicle enters water at 100 m/s and different angles at 20 ms

下载: 全尺寸图片幻灯片

图 16 航行体以100 m/s、90°入水角入水时缓冲材料的等效应力分布

Figure 16. Distribution of effective stress of the buffer when the vehicle enters water at 100 m/s and 90°

下载: 全尺寸图片幻灯片

图 17 航行体以100 m/s、60°入水角入水时缓冲材料的等效应力分布

Figure 17. Distribution of effective stress of the buffer when the vehicle enters water at 100 m/s and 60°

下载: 全尺寸图片幻灯片

图 18 不同的航行体以不同入水速度垂直入水时的加速度时程曲线

Figure 18. Time-history curves of acceleration when different vehicles enter water vertically at different velocities

下载: 全尺寸图片幻灯片

图 19 不同的航行体以100 m/s的速度在不同入水角度下的轴向加速度时程曲线

Figure 19. Time-history curves of axial acceleration when different vehicles enter water at 100 m/s and different angles

下载: 全尺寸图片幻灯片

图 20 不同的航行体以100 m/s的速度在不同入水角度下的径向加速度时程曲线

Figure 20. Time-history curves of radial acceleration when different vehicles enter water at 100 m/s and different angles

下载: 全尺寸图片幻灯片

参考文献(23)

[1]	郑强, 杨日杰, 陈佳琪, 等. 直升机空投鱼雷的散布误差研究 [J]. 科学技术与工程, 2017, 17(15): 65–70. DOI: 10.3969/j.issn.1671-1815.2017.15.009. ZHENG Q, YANG R J, CHEN J Q, et al. Research on dispersion errors of helicopter’s airdrop torpedo [J]. Science Technology and Engineering, 2017, 17(15): 65–70. DOI: 10.3969/j.issn.1671-1815.2017.15.009.
[2]	温志文, 杨智栋, 王力竟. 空投鱼雷系统建模与空中弹道仿真研究 [J]. 弹箭与制导学报, 2019, 39(5): 63–66,72. DOI: 10.15892/j.cnki.djzdxb.2019.05.015. WEN Z W, YANG Z D, WANG L J. Modeling of the air-dropped torpedo system and the simulation research of the air trajectory [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2019, 39(5): 63–66,72. DOI: 10.15892/j.cnki.djzdxb.2019.05.015.
[3]	潘龙, 王焕然, 姚尔人, 等. 头部喷气平头圆柱体人水缓冲机制研究 [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.
[4]	刘华坪, 余飞鹏, 韩冰, 等. 头部喷气影响航行体入水载荷的数值模拟 [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.
[5]	赵海瑞, 施瑶, 潘光. 头部喷气航行器高速入水空泡特性数值分析 [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.
[6]	陈洋, 吴亮, 曾国伟, 等. 带环形密闭气囊弹体入水冲击过程的数值分析 [J]. 爆炸与冲击, 2018, 38(5): 1155–1164. DOI: 10.11883/bzycj-2017-0387. CHEN Y, WU L, ZENG G W, et al. Numerical analysis of the water entry process of a projectile with a circular airbag [J]. Explosion and Shock Waves, 2018, 38(5): 1155–1164. DOI: 10.11883/bzycj-2017-0387.
[7]	严忠汉. 入水弹头缓冲器特性探讨 [J]. 水动力学研究与进展, 1987(1): 112–121. DOI: 10.16076/j.cnki.cjhd.1987.01.012. YAN Z H. An approach to the behavior of water-entry missile’s mitigator [J]. Advances in Hydrodynamics, 1987(1): 112–121. DOI: 10.16076/j.cnki.cjhd.1987.01.012.
[8]	王永虎, 石秀华, 王鹏. 雷弹入水冲击动态缓冲性能分析 [J]. 西北工业大学学报, 2009, 27(5): 707–712. DOI: 10.3969/j.issn.1000-2758.2009.05.023. WANG Y H, SHI X H, WANG P. Exploring analysis of dynamic cushioning properties of water-entry missile’s shock mitigator [J]. Journal of Northwestern Polytechnical University, 2009, 27(5): 707–712. DOI: 10.3969/j.issn.1000-2758.2009.05.023.
[9]	HIRT C W, AMSDEN A A, COOK J L. An arbitrary Lagrangian-Eulerian computing method for all flow speeds [J]. Journal of Computational Physics, 1974, 14(3): 227–253. DOI: 10.1016/0021-9991(74)90051-5.
[10]	WANG H, ZHAO F, CHENG Y S, et al. Dynamic response analysis of light weight pyramidal sandwich plates subjected to water impact [J]. Polish Maritime Research, 2012, 19(4): 31–43. DOI: 10.2478/v10012-012-0038-y.
[11]	李建阳, 邢伟, 宋世鹏, 等. 舱体入水工况参数对冲击特性的影响分析 [J]. 宇航总体技术, 2019, 3(4): 49–55. LI J Y, XING W, SONG S P, et al. Analysis of effects of water entry condition parameters on impact characteristics of recovery module [J]. Astronautical Systems Engineering Technology, 2019, 3(4): 49–55.
[12]	胡明勇, 张志宏, 刘巨斌, 等. 低亚声速射弹垂直入水的流体与固体耦合数值计算研究 [J]. 兵工学报, 2018, 39(3): 560–568. DOI: 10.3969/j.issn.1000-1093.2018.03.018. HU M Y, ZHANG Z H, LIU J B, et al. Fluid-solid coupling numerical simulation on vertical water entry of projectile at low subsonic speed [J]. Acta Armamentrii, 2018, 39(3): 560–568. DOI: 10.3969/j.issn.1000-1093.2018.03.018.
[13]	颜彬, 钱韬, 马赛尔. 火箭助飞式器材高速入水冲击结构响应分析 [J]. 弹箭与制导学报, 2018, 38(5): 65–68,72. DOI: 10.15892/j.cnki.djzdxb.2018.05.016. YAN B, QIAN T, MA S E. Analysis for structure response to water entry impact for rocket-assisted equipment at high-speed [J]. Journal of Projectile, Rockets, Missiles and Guidance, 2018, 38(5): 65–68,72. DOI: 10.15892/j.cnki.djzdxb.2018.05.016.
[14]	WU S Y, SHAO Z Y, FENG S S, et al. Water-entry behavior of projectiles under the protection of polyurethane buffer head [J]. Ocean Engineering, 2020, 197: 106890. DOI: 10.1016/j.oceaneng.2019.106890.
[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]	魏海鹏, 史崇镔, 孙铁志, 等. 基于ALE方法的航行体高速入水缓冲降载性能数值研究 [J]. 爆炸与冲击, 2021, 41(10): 104201. 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): 104201. DOI: 10.11883/bzycj-2020-0461.
[17]	钱立新, 刘飞, 屈明, 等. 鱼雷头罩入水破坏模式研究 [J]. 鱼雷技术, 2015, 23(4): 257–261. QIAN L X, LIU F, QU M, et al. Failure mode of torpedo nose cap in water-entry [J]. Torpedo Technology, 2015, 23(4): 257–261.
[18]	王泽鹏, 胡仁喜, 康士廷, 等. ANSYS 13.0 LS-DYNA非线性有限元分析实例[M]. 2版. 北京: 机械工业出版社, 2011: 348–349.
[19]	尚晓江, 苏建宇, 王化峰, 等. ANSYS LS-DYNA动力分析方法与工程实例[M]. 2版. 北京: 中国水利水电出版社, 2008: 141.
[20]	潘光, 杜晓旭, 宋保维, 等. 鱼雷力学[M]. 西安: 陕西师范大学出版社, 2013: 10.
[21]	霍银磊, 张新昌. 发泡塑料缓冲设计中材料的密度选择 [J]. 塑料工业, 2007(5): 40–43. DOI: 10.3321/j.issn:1005-5770.2007.05.012. HUO Y L, ZHANG X C. Density choice of foamed plastics for cushion design [J]. China Plastics Industry, 2007(5): 40–43. DOI: 10.3321/j.issn:1005-5770.2007.05.012.
[22]	王泽鹏, 胡仁喜, 康士廷, 等. ANSYS 13.0 LS-DYNA非线性有限元分析实例[M]. 2版. 北京: 机械工业出版社, 2011: 63.
[23]	赵海鸥. LS-DYNA动力分析指南[M]. 北京: 兵器工业出版社, 2003: 167–168.

施引文献

期刊类型引用(3)

1.	赵海瑞，施瑶，潘光，黄桥高，刘晓娅. 不同厚度航行器高速入水冲击载荷及壳体变形特性. 空气动力学学报. 2024(02): 111-125 . 百度学术
2.	肖有才，邢旭阳，杨佩琮，张宏，熊言义，徐忠四，赵志颖，孙毅. 基于ALE方法的弹体入水硬质聚氨酯泡沫缓冲器降载性能分析. 船舶力学. 2024(07): 1111-1123 . 百度学术
3.	解江，潘汉源，蒋逸伦，杨祥，李漩，郭德龙，冯振宇. 带剪切销抗爆容器定向泄压特性研究. 爆炸与冲击. 2024(07): 139-156 . 本站查看