舰船科学技术  2020, Vol. 42 Issue (3): 67-71    DOI: 10.3404/j.issn.1672-7649.2020.03.014   PDF    
潜伏式武器水下施放后初始运动过程仿真
练永庆1,2, 宋保维1, 李宗吉2     
1. 西北工业大学 航海学院,陕西 西安 710072;
2. 海军工程大学兵器工程学院,湖北 武汉 430033
摘要: 为了研究潜伏式武器从潜艇外部施放后初始阶段的运动规律,在分析潜伏式武器水下施放过程的基础上,推导武器从潜艇从运载装置施放后的水下纵平面运动数学模型,进行施放后初始阶段运动过程的仿真计算与分析。通过仿真对潜伏式武器水下施放后的初始阶段的运动规律及关键参数的影响进行了分析,仿真结果证明所建立的武器施放后初始阶段运动模型基本正确,能反映武器施放后初始阶段的基本运动规律。
关键词: 潜伏式武器     水下施放     初始运动     仿真    
Simulation on the initial motion of the latent weapon deployed underwater
LIAN Yong-qing1,2, SONG Bao-wei1, LI Zong-ji2     
1. School of Marine Science and Technology, Northwestern Polytechnical University, Xi′an 710072, China;
2. Department of Weaponry Engineering, Naval University of Engineering, Wuhan 430033, China
Abstract: In order to study the initial motion of the latent weapon deployed from the submarine, the mathematic models are deduced on the base of the analysis of the motion process after the latent weapon being deployed. The characteristics of the initial motion and the effect of the critical parameters are analysed by simulation, and the simulation result shows that the mathematic models of the initial motion are correct, and can be used in the study of the latent weapon's initial motion.
Key words: latent weapon     underwater deployment     initial motion     simulation    
0 引 言

随着海上斗争区域逐渐转向深海,深海潜伏式武器[1]等将逐渐登上战场。水下平台如潜艇作为水下隐蔽作战平台,将是运送和搭载该类武器的最佳平台之一。由于所研究的潜伏式武器口径大于现有鱼雷发射管口径,因此无法通过鱼雷发射装置[2]进行发射或布放。因此采用外部挂载的方式是潜艇等水下平台携带及施放潜伏式武器的有效方法之一。潜艇外搭载武器方面国外早已有一定的应用[3-5],国内也开展了部分研究[6]。在潜伏式武器以及水下平台外挂装置的总体论证研究中,武器施放后的初始运动规律是其所关注的重点之一,为此本文通过建立潜伏式武器施放后水下运动过程相关数学模型并进行仿真计算,对潜伏式武器施放后初始阶段的运动特性进行进行定量分析,从而为潜伏式武器及潜艇外挂装置的总体方案设计和论证提供理论依据。

1 潜伏式武器水下施放方案简介

本文研究的武器水下施放的方案如见图1所示。该方案的运载装置外挂架通过连接紧固机构固定在潜艇底部,在潜艇航行到预定区域需投放武器时,首先通过艇内的仪器对武器进行初始发射参数设定。完成设定后,则通过运载装置的操纵控制设备首先对模块储放箱的底部挡板进行解锁,再通过液压开闭结构将底部挡板打开,这时武器在其负浮力的作用下离开外挂架,开始初始运动。本文主要对武器施放后的水下纵平面初始运动进行建模与仿真研究。

图 1 施放方案原理图 Fig. 1 Schematic diagram of the weapon underwater deployment
2 潜伏式武器水下施放后纵平面运动数学模型 2.1 建模坐标系简介

在武器水下纵平面运动建模与仿真中,采用2种坐标系:固定坐标系 $o - xyz$ (简称“定系”)和运动坐标系 ${o_1} - {x_1}{y_1}{z_1}$ (简称“动系”),如图2所示。

图 2 固定坐标系和运动坐标系 Fig. 2 Fixed coordinate system & moving coordinate system

定系 $o - xyz$ 固定在地球,原点 $o$ 选在施放瞬间潜艇与武器交界面上的某一点。 $ox$ 轴位于水平面上,以潜艇的主航向为正向; $oy$ 轴位于 $ox$ 轴所在的水平面,按右手法则将 $ox$ 轴顺时针旋转 ${90^ \circ }$ 即是; $oz$ 轴垂直于 $xoy$ 的坐标面,以指向地心为正。

动系 ${o_1} - {x_1}{y_1}{z_1}$ 固联于武器,随武器一起运动,原点 ${o_1}$ 取在潜伏式武器的形心处。 ${o_1}{x_1}$ 轴、 ${o_1}{y_1}$ 轴和 ${o_1}{z_1}$ 轴,分别是经过中点的水线面、横剖面和纵中剖面的交线,正向按右手系的规定,即以 ${o_1}{x_1}$ 轴指向武器的前部、 ${o_1}{y_1}$ 轴指向右侧、 ${o_1}{z_1}$ 轴指向水平龙骨为正,并认为 ${o_1}{x_1}$ ${o_1}{y_1}$ ${o_1}{z_1}$ 是武器的惯性主轴。

为了简化建模过程,建模中不考虑水下海流的影响。

2.2 潜伏式武器水下运动数学模型

当武器施放后,武器在自身负浮力作用下作下沉运动。同时由于武器重心与浮心不在同一位置,武器还将绕其形心作旋转运动。因此该阶段的运动为武器整体向下运动的同时作绕自身的旋转运动,经推导,该阶段的纵平面运动方程如下:

横向力方程( ${x_1}$ 轴向力)

$ \begin{split} & {{m}}\left[ {\mathop u\limits^ \bullet + wq + {z_G}\mathop q\limits^ \bullet - {x_G}{q^2}} \right] = \\ & \frac{1}{2}\rho {L^3}{X'_{\dot u}}\mathop u\limits^ \bullet + {X_u}u - \left( {W - B} \right)\sin \theta {\text{;}} \end{split} $ (1)

竖向力方程( ${z_1}$ 轴向力):

$ \begin{split} & {{m}}\left[ {\mathop w\limits^ \bullet - u{q^2} - {x_G}\mathop q\limits^ \bullet - {z_G}{q^2}} \right] = \\ & \frac{1}{2}\rho {L^4}{{\rm{Z'}}_{\dot q}}\mathop q\limits^ \bullet + \frac{1}{2}\rho {L^3}{{\rm{Z'}}_{\dot w}}\mathop w\limits^ \bullet + \frac{1}{2}\rho {L^3}{{\rm{Z'}}_q}q + \frac{1}{2}\rho {L^2}{Z'_w}w +\\ & \left( {W - B} \right)\cos \theta{\text{;}} \end{split} $ (2)

纵倾力矩方程( ${y_1}$ 轴向力矩)

$ \begin{split} & {J_y}\mathop q\limits^ \bullet + {\rm{m}}\left[ {{z_G}\left( {\mathop u\limits^ \bullet + wq} \right) - {x_G}\left( {\mathop u\limits^ \bullet - uq} \right)} \right] = \\ & \frac{1}{2}\rho {L^5}{M'_{\dot q}}\mathop q\limits^ \bullet + \frac{1}{2}\rho {L^4}{M'_{\dot w}}\mathop w\limits^ \bullet + \frac{1}{2}\rho {L^4}{M'_q}q + \frac{1}{2}\rho {L^3}{M'_w}w+\\ & \left( {{x_G}W - {x_C}B} \right)\cos \theta - \left( {{z_G}W - {z_C}B} \right)\sin \theta {\text{;}} \end{split} $ (3)

式中: $m$ 为武器质量; $L$ 为武器长度; $\rho $ 为海水密度; ${J_y}$ 为武器对动坐标系 $y$ 轴的转动惯量;WB分别为武器重力与浮力; $u$ $w$ 分别为武器速度在动坐标系 $x$ $z$ 轴上的投影; $q$ 为角速度 $\varOmega $ 在动系 ${y_1}$ 轴上的投影; $\theta $ 为武器纵倾角;xGzG分别为武器重心 $G$ 在固定坐标系坐标;xCzC分别为武器浮心C在固定坐标系坐标;其余参数为武器的无因次水动力及力矩系数,具体含义可参见文献[7]。

解以上方程可得在运动坐标系下的潜伏式武器任一时刻质心的加速度和绕质心的角加速度。

武器在固定坐标系下任一时刻的位置为:

$ \begin{split} & {x_{\rm{G}}} = \int {{v_x}} {\rm{d}}t {\text{,}}\\ & {z_{\rm{G}}} = \int {v_z} {\rm{d}}t {\text{。}} \end{split} $ (4)

式中:vxvz分别为武器在固定坐标系下x轴、z轴方向速度,可按下式计算:

$ \left[ {\begin{array}{*{20}{l}} {{v_x}}\\ 0\\ {{v_z}} \end{array}} \right] = C_o^o {\text{。}} $ (5)

式中: $C_o^{{o_1}}$ 为运动坐标系到固定坐标系的转换矩阵; $u$ $w$ 分别为武器速度在动坐标系 $x$ $z$ 轴上的投影可按下式计算:

$\begin{gathered} u = \int {\dot u{\rm{d}}t} {\text{,}} \\ w = \int {\dot w{\rm{d}}t} {\text{。}} \end{gathered} $ (6)

武器在任一时刻的姿态角为:

$ \theta = \iint {\dot q{\rm d}{t^2}} {\text{。}} $ (7)
3 仿真结果与分析 3.1 标准参数仿真结果

第一阶段的主要仿真参数为:武器初速度 $u = 0$ ;武器初始负浮力 $W - B = 1\,000\;{\rm{kg}}$ ;武器转动惯量 ${I_y} = 15\,000\;{\rm{kg}} \cdot {{\rm{m}}^{\rm{2}}}$ ;武器的重、浮心坐标: $G(0.4,0)$ C(–0.2,0)。

将参数的初始值代入到仿真模型中进行运行仿真,得到武器的潜深曲线、速度 ${v_x}$ 曲线、速度 ${v_z}$ 曲线、纵倾角变化曲线(见图3图7)。

图 3 潜深曲线图 Fig. 3 The curve of the weapon diving depth in vertical plane

图 7 0-5 s内姿态变化图 Fig. 7 The weapon's movement posture changes in 5 s

为了更加直观反映武器在下落时初始阶段的姿态变化,仿真中给出初始段(0–5 s)武器的姿态变化图(见图7),图中直线表示为武器中心轴线,圆圈表示武器底部。

从以上仿真结果可以看出:

1)武器在施放后,由于武器的重心与形心(动坐标原点)不重合,因此武器除了自身转动外(见图7),其重心还围绕施放点垂线作类似钟摆运动(见图3),并经过一段时间的摆动后,最终保持竖直下沉方向。

2)由武器在x轴方向的速度 ${v_{\rm{x}}}$ 仿真结果(见图4)可见,速度 ${v_{\rm{x}}}$ 先增大后减小,此后在 ${v_{\rm{x}}} = 0$ 附近震荡变化,其幅值逐渐减小并趋向于0;

图 4 ${v_x}$ 速度曲线 Fig. 4 The curve of the weapon's x direction velocity in vertical plane

3)由武器在z轴方向的速度 ${v_{\rm{z}}}$ 仿真结果(见图5)可见,由于有海水阻力的影响, ${v_{\rm{z}}}$ 速度开始时是波动式递增的,并最终稳定在 ${v_{\rm{z}}} = 6\;{\rm{m}}/{\rm{s}}$ 左右。

图 5 ${v_z}$ 速度曲线 Fig. 5 The curve of the weapon's z direction velocity in vertical plane

4)由纵倾角仿真曲线(见图6)可见,因武器绕自身转动,则纵倾角上则是在 $\theta = {90^ \circ }$ 附近震荡,并最终经过一段时间后,最终纵倾角趋向于 ${90^ \circ }$

图 6 纵倾角变化曲线 Fig. 6 The curves of the weapon's pitch angles
3.2 变参数仿真结果

为了进一步分析潜艇平台及武器的主要参数对其水下运动的影响,在标准参数仿真模型的基础上,逐一改变标准参数中的武器的负浮力、初始速度以及转动惯量,分别进行仿真以研究这些参数对武器的施放过程以及运动的影响规律。

1)改变武器负浮力的仿真

将武器的负浮力在标准参数基础上分别设定为减少20%( $800\;{\rm{kg}}$ )、不变( $1\;000\;{\rm{kg}}$ )、增加20%( $1\;200\;{\rm{kg}}$ )进行仿真,仿真的结果如图8图11所示。

图 8 潜深曲线(变负浮力) Fig. 8 The curve of the weapon diving depth in vertical plane according to the negative buoyancy

图 11 纵倾角变化曲线(变负浮力) Fig. 11 The curves of the weapon's pitch angles according to the negative buoyancy

从仿真的结果可以看出:

武器的负浮力对武器的稳定后的z方向的运动速度 ${v_z}$ 影响较大,随着负浮力的增大而 ${v_z}$ 相应增大(见图10)。当 $W - B = 800\;{\rm{kg}}$ 时, ${v_z} = 5.1\;{\rm{m}}/{\rm{s}}$ ;当 $W - B = 1\;000\;{\rm{kg}}$ 时, ${v_z} = 6.0\;{\rm{m}}/{\rm{s}}$ ;当 $W - B = 1\;200\;{\rm{kg}}$ 时, ${v_z} = 6.6\;{\rm{m}}/{\rm{s}}$

图 10 ${v_z}$ 速度曲线(变负浮力) Fig. 10 The curve of the weapon's z direction velocity in vertical plane according to the negative buoyancy

武器的负浮力对 ${v_x}$ 以及纵倾角 $\theta $ 的随时间震荡的幅值有一定的影响,随着负浮力的增加, ${v_x}$ 以及纵倾角 $\theta $ 的在相同时刻的振幅相应增大(见图9图11),但对稳定后的速度 ${v_x}$ 无影响。

图 9 ${v_x}$ 速度曲线(变负浮力) Fig. 9 The curve of the weapon's x direction velocity in vertical plane according to the negative buoyancy

2)改变武器初始速度的仿真

仿真中的武器初始速度可以认为是潜艇平台在运动中进行武器施放时武器在离开潜艇瞬间将带上的初速。将武器的初始速度分别设定为 $u = 0\;{\rm{m}}/{\rm{s}}$ $1\;{\rm{m}}/{\rm{s}}$ $2\;{\rm{m}}/{\rm{s}}$ 进行仿真,仿真的结果如图12图15所示。

图 12 潜深曲线(变初速) Fig. 12 The curve of the weapon diving depth in vertical plane according to its initial velocity

图 15 垂直面的纵倾角变化曲线(变初速) Fig. 15 The curves of the weapon's pitch angles according to its initial velocity

图 13 ${v_x}$ 速度曲线(变初速) Fig. 13 The curve of the weapon's x direction velocity in vertical plane according to its initial velocity

从仿真的结果可见,改变武器的初始速度,主要对武器落点位置有影响,初始速度越大,武器落点离施放开始位置的距离越大,如初始速度 $u = 1\;{\rm{m}}/{\rm{s}}$ 时,落点将距离施放点约15 m左右、当 $u = 2\;{\rm{m}}/{\rm{s}}$ ,落点将距离施放点约30 m左右(见图12)。而且由于海水阻力的作用,武器 $x$ 轴方向的速度终将趋于0。武器的初始速度对 ${v_z}$ 的初始增加段有一定影响,对稳定后的速度值无影响(见图14),武器的初始速度纵倾角 $\theta $ 基本没有影响(见图15)。

图 14 ${v_z}$ 速度曲线(变初速) Fig. 14 The curve of the weapon's z direction velocity in vertical plane according to its initial velocity

3)改变武器转动惯量的仿真

将武器绕y1轴的转动惯量标准参数基础上分别设定为减少10%( $1.35 \times {10^5}\;{\rm{kg}} \cdot {{\rm{m}}^2}$ )、不变( $1.5 \times $ $ {10^5}\;{\rm{kg}} \cdot {{\rm{m}}^2}$ )、增加10%( $1.65 \times {10^5}\;{\rm{kg}} \cdot {{\rm{m}}^2}$ )进行仿真,仿真的结果如图16图19所示。

图 16 潜深曲线(变转动惯量) Fig. 16 The curve of the weapon diving depth in vertical plane according to its moment of inertia

图 19 纵倾角变化曲线(变转动惯量) Fig. 19 The curves of the weapon's pitch angles according to its moment of inertia

图 17 ${v_x}$ 速度曲线(变转动惯量) Fig. 17 The curve of the weapon's x direction velocity in vertical plane according to its moment of inertia

图 18 ${v_z}$ 速度曲线(变转动惯量) Fig. 18 The curve of the weapon's z direction velocity in vertical plane according to its moment of inertia

从仿真的结果可以看出,改变武器的转动惯量, ${v_x}$ 、纵倾角 $\theta $ 随着转动惯量的增大而发生变化,武器转动惯量大,武器自身转动幅度大,使趋于稳定的时间向后发生推移,而武器的落点位置以及 ${v_z}$ 没有明显变化。

4 结 语

为了研究潜伏式武器水下施放后初始阶段的运动规律,建立潜伏式武器水下施放后的初始阶段运动模型并进行仿真分析。从仿真结果可知,无动力的潜伏式武器施放后,由于武器的重心和浮心不在其形心处,因此在投放后的武器初期运动为形心的下沉运动与绕形心钟摆式转动,这种摆动随着下沉深度加大幅度逐渐减小直至成垂直下沉状态。在这阶段中,武器水平初始速度、负浮力以及转动惯量等对其运动有影响。水平初始速度主要对武器落点位置有影响,初始速度越大,武器落点离施放开始位置的距离越大;负浮力对武器稳定后z方向的运动速度 ${v_z}$ 影响较大,随着负浮力的增大而 ${v_z}$ 相应增大;而转动惯量则主要影响投放开始过程武器自身转动幅度及趋于稳定的时间。从仿真结果分析可见,所建立的模型基本正确,能反映武器施放后各阶段的基本运动规律。后续将在现有模型基础上增加海底涌浪影响模块,完善整个运动模型,以便更好地为潜伏式武器及水下平台外挂装置的总体设计论证提供理论支撑。

参考文献
[1]
宋保维, 朱信尧, 梁庆卫, 曹永辉. 潜伏式无人水下航行器概念设计[J]. 火力指挥与控制, 2010(8): 107-110.
[2]
练永庆, 王树宗. 鱼雷发射装置设计原理[M]. 2012: 125~134.
[3]
范乃忠. 现代潜艇的外挂布雷装置[J]. 现代舰船, 1995, 6(2): 32-34.
[4]
钱东. 美国未来的大型UUV-MANTA[J]. 鱼雷技术, 2003(1): 41-50.
[5]
SIRMALIS J E. Pursuing the MANTA vision: recent at-sea technology demonstration results[C]. UDT 2001.
[6]
杨路春. 潜艇在有外部搭载情况下操纵性水动力导数的数值计算方法研究[D]. 哈尔滨: 哈尔滨工程大学, 2009.
[7]
MARTLON G, GRANT P. H. Standard Equations of motion for submarine simulation[R]. Naval Ship Research and Development Center, AD653861 June 1967.