鲁春佳，侯 健，魏 平

（海军工程大学兵器工程学院，湖北 武汉 430032）

全浸水带间隙发射作为一种新的水下发射方式[1]，通过间隙燃气在膛内贴壁运动，卷吸回流后逐渐汇聚成弹前气幕，排出弹前水柱，将射弹在膛口的发射环境由水介质转化为气体介质。当气体射流流出枪口后，膛口处射流迅速膨胀成球形气体空腔，射弹穿过气体空腔与水介质接触，产生强烈的冲击载荷。水下高速射弹的弹体主要由硬铝合金尾杆和钨合金头部组成，两者镶嵌连接，连接强度有限，入水瞬间弹头会承受巨大冲击，因此射弹的入水冲击载荷成为水下射弹设计中的一个重要问题。

Karman[2]最先开始对入水冲击现象进行研究，采用动量定理并引入附加质量的概念，推导出入水冲击载荷的计算公式。Wagner[3]将Karman 的方法理论化，提出了近似平板理论及自相似解法，得出了冲击压力在结构沾湿面的分布情况，使理论分析更加符合实际情况，为后来学者的理论研究奠定了基础。在国内，秦洪德等[4]、王永虎等[5]对入水冲击问题的现状和进展进行了详细的分析。宋保维等[6]基于不可压缩的非定常势流理论，建立了空投水雷入水冲击计算的数学模型。卢炽华等[7]利用不同浸深的附加质量，对刚性细长体斜姿态落水冲击进行建模，得出其入水角很小，会使弹体处在最危险的状态。王永虎等[8-9]先后对刚性尖拱体垂直姿态高速入水和斜入水的冲击理论进行了建模和仿真。魏卓慧等[10]建立了刚性截锥弹体垂直入水冲击载荷的数学模型，并对其进行了数值计算。陈诚等[11]对超空泡航行器倾斜入水冲击载荷进行了试验研究，得出了峰值时刻的阻力系数。朱珠等[12]利用商业软件FLUENT，建立了柱体回转体高速入水冲击的数值模拟模型，得到了速度对入水冲击载荷的影响规律。然而以上研究主要针对由空中入水的冲击载荷分析，对于水下入水冲击问题研究较少。本工作在此前提条件下，计算分析全水下发射高速射弹入水的冲击载荷，这对于水下发射武器研究具有一定的现实意义。

本研究拟建立锥形弹体水平及斜入水的冲击载荷理论模型，模型中考虑弹体重力、弹体浮力、附加质量、弹头锥角及入水攻角的影响，对不同头部结构参数的锥形弹体以不同入水速度入水的冲击载荷进行计算，分析入水速度、弹头锥角和入水攻角对冲击载荷的影响。研究结果对于弹体入水冲击载荷的预测及全水下发射方式发射的射弹头部结构设计具有参考价值。

1 数学模型

1.1 射弹水平入水

假设射弹为刚体，不考虑射弹入水时空泡的影响，根据动量定理，射弹高速入水冲击时的动量方程为

式中：M 为射弹的质量，m 为射弹的附加质量， Fd为 射弹入水时受到的阻力， Fb为射弹所受的浮力。

对式（1）等号两边进行微分，得到射弹入水冲击时的动力学方程

式中：h 为射弹侵入的距离，A 为阻力面积， Cd为 阻力系数， ρ为水的密度。

射弹沿 x 轴水平入水时，由于入水冲击过程的瞬时性，入水初期其运动方向基本保持不变，流体动力主要作用于射弹轴线方向，在射弹轴线方向上重力和浮力对射弹影响很小，基本可忽略不计。

射弹入水冲击过程中，射弹浸没在水中的体积会排挤液面流体产生隆起现象，如图1 所示，有效液面决定了自由水平液面的抬高程度，这取决于射弹锥头的外形和入水角等初始状态。沾湿因子定义为有效液面与实际液面的比值。

利用轴长体假设和Tayler 关于在不同浸深时附加质量的表达式[13]，参考垂直入水相关文献[10,14]，求出锥形射弹入水的附加质量为

图1 射弹水平入水冲击示意图Fig. 1 Schematic diagram of horizontal impact of the projectile

式中： Cw为 沾湿因子， β为射弹头部的半锥角，如图1 所示。

将式（3）～式（5）代入式（2），得到总方程为

采用MATLAB 软件，利用龙格-库塔方法进行求解，可以得到弹体入水时的冲击载荷。

1.2 带攻角入水

如图2 所示，采用全新的发射方式时，燃气排出并在膛口形成球形气体空腔。以地面为坐标系，枪口斜向上（即 y 轴正方向）、斜向下（即 y轴负方向）发射时，根据圆切线定理，穿过球形气腔仍可看作垂直于液面入水，然而射弹入水在有攻角的情况下，轴向上会受到重力和浮力的作用分力影响，攻角正负值相反时，重力与浮力在射弹轴向上的作用分力方向也相反。

其他条件与水平入水时保持不变，攻角为正时，根据动量方程，得到动力学方程

图2 射弹攻角示意图Fig. 2 Angle of attack of the projectile

式中：α 为射弹攻角，射弹斜向上发射时为正值，斜向下时为负值。

射弹所受浮力为

当攻角为正值时，得到的总方程为

同理，攻角为负值时的总方程为

对式（9）、式（10）进行求解，可以得到不同攻角下的入水冲击载荷。

2 计算结果与分析

由于水下射弹的质量较轻，入水速度较大，水下超空泡射弹入水平均速度约为600 m/s，因此入水速度对射弹入水冲击载荷的影响很大。本计算中，设射弹质量为0.14 kg，射弹头部半锥角为6°，计算得到不同速度时射弹的入水冲击载荷曲线，如图3 所示，其中用射弹加速度反映入水冲击载荷。可以看出：射弹的入水冲击载荷先增大后减小，载荷变化主要发生在射弹入水前1 ms，最后渐渐趋于稳定；入水速度越大，冲击载荷峰值越大，入水后达到峰值的时间越短。此外，计算了不同速度下的入水冲击载荷峰值，如图4 所示。入水速度在400～700 m/s 范围时，射弹入水冲击载荷峰值一般为103g 量级（g 为重力加速度）。从冲击载荷峰值与速度的关系可以得出入水冲击载荷的峰值与速度基本呈线性关系。

图3 不同速度时的入水冲击载荷Fig. 3 Driving impact loads at different velocities

图4 不同速度时的冲击载荷峰值Fig. 4 Peak impact loads at different velocities

设射弹质量为0.14 kg，通过计算获得了不同锥角的锥形射弹以600 m/s 入水时的冲击载荷曲线，如图5 所示。可以看出：射弹的入水冲击载荷先增大后减小，最后趋于稳定；半锥角 β越大，冲击载荷峰值越大，并且入水后达到峰值的时间越短。改变半锥角，计算出不同半锥角情况下射弹入水冲击载荷峰值，如图6 所示。从冲击载荷峰值与半锥角的关系可以得出入水冲击载荷峰值与半锥角基本呈线性关系。

图5 不同锥角时的入水冲击载荷Fig. 5 Driving impact loads at different cone angles

图6 不同锥角时的冲击载荷峰值Fig. 6 Peak impact loads at different cone angles

当射弹质量为0.14 kg，入水速度为600 m/s，攻角 α分别取45°、-45 °和0°（即水平入水）时，计算得到的射弹入水冲击载荷曲线如图7 所示。可见，3 条曲线基本重叠，差值在1g 量级，相比于速度和半锥角对冲击载荷的影响，攻角对冲击载荷的影响几乎可以忽略不计。

此外，射弹在高速入水状态下，从液面进入水中的实际深度与理想深度存在一定的偏差，为了表达该误差，引入沾湿因子，沾湿因子的取值对冲击载荷的结果也存在影响。沾湿因子不同时射弹的入水冲击载荷如图8 所示。从图8 可以看出：沾湿因子越大，入水冲击载荷峰值越大；但沾湿因子对入水冲击 载荷的影响较小，当沾湿因子变化值为0.1 时，入水冲击载荷的变化在10%以内。

图7 不同攻角时的入水冲击载荷Fig. 7 Inlet impact loads at different angles of attack

图8 沾湿因子不同时的入水冲击载荷Fig. 8 Inlet impact loads with different wet factors

3 结 论

（1）射弹头部锥角相同时，入水冲击载荷峰值与速度呈正线性相关，入水速度越大，冲击载荷达到峰值的时间越短；射弹入水速度相同时，入水冲击载荷的大小与锥角呈正线性相关，锥角越大，冲击载荷达到峰值的时间越短。

（2）射弹锥角和入水速度相同、入水攻角不同时，入水冲击载荷曲线与水平入水曲线基本重合，说明射弹重力和浮力在轴向上对入水冲击载荷的影响相对入水阻力几乎可以忽略不计。当射弹带有攻角入水后，重力和浮力更多的是对射弹产生径向力矩影响。

（3）沾湿因子越大，入水冲击载荷峰值越大，但对入水冲击载荷影响较小，当沾湿因子变化为0.1 时，入水冲击载荷的变化在10%以内。沾湿因子作为一个变量，其大小反过来也取决于入水冲击载荷，入水冲击载荷越大，沾湿因子也越大。

（4）理论模型借鉴了高速弹体垂直入水的理论模型，数值模拟计算结果与已报道的高速射弹垂直入水冲击载荷数值计算和仿真结果高度一致，验证了本数学建模和数值计算的准确性。

（5）探讨了全水下发射高速射弹入水瞬间的冲击载荷，高速射弹入水后形成超空泡，冲击载荷迅速减小。本工作对射弹未形成超空泡的情况进行了模型推导和数值模拟计算，对形成超空泡之前的理论研究具有重要意义。