李慧, 赵琳, 毛英

(哈尔滨工程大学 自动化学院，黑龙江 哈尔滨 150001)

海况干扰下潜艇六自由度运动分析

李慧, 赵琳, 毛英

(哈尔滨工程大学 自动化学院，黑龙江 哈尔滨 150001)

为了弥补实际海试条件的局限性，辅助实现陆上鉴定敏感潜艇运动的导航系统设备的性能，本文对潜艇六自由度运动规律进行了研究，包括静水以及海洋环境干扰下潜艇的直航和回转运动。通过分析潜艇运动时所受的力及力矩，以及海浪、海流的基本模型，建立海洋环境干扰下的潜艇六自由度运动方程。针对潜艇的线运动、角运动相关信息，结合已有的经过实际试验修正的某潜艇的部分水动力系数，对其直航、回转运动进行了研究分析，仿真结果能够较为真实地反应潜艇在海洋环境干扰下的运动情况，满足实际工程中辅助测试的需求。

潜艇；运动模型；六自由度；海洋环境；海浪干扰；海流干扰

对于敏感潜艇运动的导航系统而言，在装备前要对其进行性能鉴定，需要尽可能地遍历潜艇出海时的各种运动，包括各种海洋环境干扰。实艇试验成本高、风险大[1]，而且由于复杂海况的不可预见性，造成潜艇的短周期运动具有不可控性，海况(尤其是复杂海况)干扰造成的潜艇受迫运动无法实现人为控制，无法全面考察设备性能。这使得全面分析潜艇在海洋环境干扰下的六自由度运动十分有意义。利用陆地上建立的真实的潜艇模拟仿真环境，一定程度上弥补实际海试时条件的局限性，是一种很好的解决办法。

本文研究了海况干扰下的潜艇运动，首先介绍了描述潜艇运动时所用到的坐标系，然后以格特勒的潜艇标准运动方程为基础，进行适当简化，结合海洋环境干扰模型，建立海况干扰下的潜艇六自由度运动方程。以某潜艇为研究对象，进行静水和海况干扰下的直航、回转运动仿真分析，最后得出潜艇六自由度运动的速度和姿态的变化特点。

1 坐标系

描述潜艇六自由度运动时采用了固定坐标系(或称地面坐标系)及运动坐标系(或称艇体坐标系)，按研究惯例，一律采用右手坐标系，参数、符号体系与国际水池会议(ITTC)及造船和轮机工程学会(SNAME)术语公报推荐的一致[2]。

固定坐标系的原点E可以任意选取海面或海中的某一点，本文将t=0时刻潜艇重心所在位置定为坐标原点，Eξ轴沿水平方向，以潜艇的主航向作为正向，Eη轴位于Eξ轴所在的水平面，且由Eξ轴顺时针旋转90°得到，Eζ轴垂直于Eη轴和Eξ轴构成的平面，其正向符合右手定则。

运动坐标系的坐标原点O通常选取在艇体的重心G处，Ox轴位于艇体纵中剖面内，沿水平线指向船艏方向为正，Oy轴与纵中剖面垂直，以指向右舷方向为正，Oz轴垂直于Ox轴与Oy轴所在的平面，以指向船底方向为正，固定坐标系及运动坐标系如图1所示。

图1 固定坐标系及运动坐标系Fig.1 Fixed coordinate system and moving coordinate system

设两个坐标系原点重合，则固定坐标系可以通过三次正交旋转与动系完全重合。表1中列出了本文在描述潜艇运动时常用参量的符号表示。

表1 运动坐标系的各参量符号

2 潜艇空间运动模型

2.1 潜艇六自由度运动模型

潜艇六自由度空间运动方程[3-4]以1967年泰勒海军潜艇研究和发展中心发表的格特勒《用于潜艇模拟研究的标准运动方程》为基础，该方程含108个水动力系数，实际中由于试验的困难很难全部得到，在不影响潜艇六自由度运动的主要特性下，可对标准方程进行简化处理，假设条件主要有[5-6]：1)潜艇左右对称；2)潜艇重心与原点重合；3)潜艇运动状态对舵力没有影响，且不考虑舵之间的相互耦合；4)忽略对潜艇运动影响非常小的二阶水动力系数及其耦合系数；5)忽略潜艇垂直面的速度对水平面的运动耦合。

根据以上简化条件，得到简化后的潜艇六自由度运动力与力矩的方程如式(1)～(6)，其中前三个方程为力的方程，后三个方程为力矩方程。

轴向力方程：

(1)

横向力方程：

(2)

垂向力方程：

(3)

横倾力矩方程：

(4)

纵倾力矩方程：

(5)

偏航力矩方程：

(6)

结合潜艇六自由度运动辅助方程为[7]

(7)

若已知潜艇水动力系数，对于任何操纵信息(包括舵、桨等)，采用四阶龙格-库塔积分算法对潜艇六自由度运动模型进行求解，可快速得到舰潜艇动状态的数值解，包括潜艇的三维速度、角速度、姿态及轨迹。

2.2 海洋环境干扰模型

潜艇的战术任务主要是在近水面完成的，如水深十多米处发射导弹等，不可避免地受到海浪和海流的影响而产生摇荡[8]，为了更真实的模拟潜艇在海洋中的运动状况，本文对海浪和海流干扰下的潜艇运动进行了仿真研究。

2.2.1 作用在潜艇上的海浪干扰模型

为了提高仿真模型的精度，同时满足实时仿真的需求，本文采用十二届ITTC提供的双参数随机海浪频谱来实时仿真长峰波随机海浪。

ITTC提供的双参数海浪谱为[9]

(8)

根据能量等效原则，遭遇频率能量谱密度函数与自然频率能量密度函数之间的关系为

(9)

式中：V为潜艇航速，β为浪向角，g为重力加速度。

潜艇在海浪作用下所受到力及力矩可主要考虑为横摇力矩，纵摇力矩及垂向作用力，其受力模型分别介绍如下：

1)潜艇受海浪作用横摇力矩

作用于潜艇上的海浪横摇力矩[10]：

(10)

海浪作用下潜艇横摇的波倾角仿真模型为

式中：n为谐波个数，ωi为第i个海浪谐波的频率，Δω为频率增量，Sα(ωi)为第i个谐波波倾角能量谱密度。

2)潜艇受海浪作用纵摇力矩

作用于潜艇上的海浪纵摇力矩[10]：

(12)

本文中所用到的海浪作用下的潜艇纵摇的波倾角仿真模型为

(13)

3)潜艇垂向受海浪的作用力

本文采用Hirom近似公式来计算作用于潜艇垂向上的海浪一阶垂荡力[11]，当潜艇深度变化值和纵摇角都较小的情况下，潜艇在近水面所受的波浪力为

(14)

2.2.2 作用在潜艇上的海流干扰模型

海流的处理，采用均匀定常流的假设，对于非高速运动的潜艇，只是从运动学上考虑其对潜艇的影响。若VC为流速，ψC为流向，则流速VC在附体坐标轴上的分量为

(15)

潜艇运动对水流的相对速度的投影ur、vr：

(16)

将海浪模型中的力和力矩加入到潜艇六自由度运动方程，引入海浪干扰信息，而对于海流，由于只考虑运动学上的影响，故按上式进行处理。

3 某潜艇六自由度运动仿真试验

水下回转是潜艇转向运动的重要研究内容[12]，本文根据上小节中建立的潜艇六自由度运动模型及海洋环境干扰模型，以某潜艇为仿真对象，给出了该潜艇在静水中、海洋环境干扰下的直航及回转运动情况。

3.1 静水中潜艇回转

静水中潜艇直航相对较简单，本文仅给出静水中回转运动的仿真，设定初始航速20 kn，航向角为0°，即正北向航行，调整螺旋桨推力以固定航速，方向舵角分别给为10°、20°、30°，艏艉升降舵角为0，即不进行定深控制。

(a)运动轨迹

(b)纵向速度曲线

(c)横摇角速率曲线

(d)横摇角曲线

(e)纵摇角速率曲线

(f)纵摇角曲线

(g)艏摇角速率曲线

(h)艏摇角曲线

图2为该潜艇在静水中固定航速改变方向舵角进行的回转运动，图中列出了某潜艇的轨迹、轴向速率、三维角速率以及三维姿态角图。由图中可以看出，潜艇方向舵角由10°、20°、30°依次增大时，潜艇在回转过程中的漂角也相应增大，其回转圈的半径越小；潜艇纵向阻力依次增大，其速降越严重；潜艇定常回转时的向心力也依次增大，其横向力应等于向心力在横向的投影，横向力对重心的力矩与扶正力矩平衡，潜艇定常回转时的横倾角也相应的增大。

如果方向舵角固定为20°，调节螺旋桨推力使潜艇的速度分别固定在9、15、20 kn，其运动状态如图3。

(a)运动轨迹

(b)纵向速度曲线

(c)横摇角速率曲线

(d)横摇角曲线

(e)纵摇角速率曲线

(f)纵摇角曲线

(g)艏摇角速率曲线

(h)艏摇角曲线

图3为潜艇在静水中固定舵角改变航速进行的回转运动，从三种速度下的航迹图可以看出，航速对潜艇定常回转直径的影响不大；潜艇定常回转时其横倾角与航速的平方存在正比关系，因而随着航速增大，横倾角也相应增大，稳态时的横倾角也越大，当航速为20 kn时，稳定横倾角达3.8°。纵倾角也随着航速的增大而增大，当航速为20 kn时，稳定纵倾角达到5°。

3.2 海浪对潜艇运动的干扰作用

海浪的大小决定了海面的粗糙程度，进一步影响了潜艇的运动，本节在静水的基础上加上海浪干扰对潜艇运动影响进行仿真。设潜艇航速20 kn，航向角为0°，即正北向航行，初始深度为5 m，方向舵角为0°，首尾升降舵角为0°，不进行定深控制。按实际海域统计得到的海浪参数对仿真模型进行修正，其中海浪浪高为1.28 m，浪向角为30°，潜艇的轨迹及横摇和纵摇曲线如图4所示。

(a)运动轨迹

(b)纵向速度曲线

(c)横摇角速率曲线

(d)横摇角曲线

(e)纵摇角速率曲线

(f)纵摇角曲线

(g)艏摇角速率曲线

(h)艏摇角曲线

图4中给出了在海浪干扰下潜艇直航运动的状态。受海浪干扰力的作用，潜艇在水下直航运动时存在纵横摇，其中横摇角最大达到6°，纵倾角达7°。

保持上述仿真条件不变，将方向舵角改为20°，进行海浪干扰下，潜艇回转运动的仿真，其运动轨迹及横纵摇曲线如图5所示。对比潜艇在静水中回转运动时的角速率及姿态角参数，在海浪干扰下，潜艇的横向、纵向以及垂向角速度出现高频的震荡，相应的，其三维姿态角随之变化，横摇、纵摇较为明显。

与图4中的直航运动状态比较，潜艇受回转力和力矩的作用，其三维角速率变大，横摇角不再以0°为中心上下震荡，而是以-4°为中心上下震荡，纵倾角出现高频变化量。

(a)运动轨迹

(b)纵向速度曲线

(c)横摇角速率曲线

(d)横摇角曲线

(e)纵摇角速率曲线

(f)纵摇角曲线

(g)艏摇角速率曲线

(h)艏摇角曲线

为了更进一步体现海浪对潜艇的作用效果，加强海浪干扰，将波高设置为3.25 m，潜艇的轨迹及横纵摇曲线如图6所示。

(a)运动轨迹

(b)纵向速度曲线

(c)横摇角速率曲线

(d)横摇角曲线

(e)纵摇角速率曲线

(f)纵摇角曲线

(g)艏摇角速率曲线

(h)艏摇角曲线

对比图5，海浪的浪向保持不变，波高由1.28 m变为3.25 m，潜艇的三维角速率以及三维姿态角明显变大，横摇角及纵倾角高达11°。

3.3 海流对潜艇运动的干扰作用

本小节仿真了潜艇在不同等级海流下的直航和回转运动。潜艇初始航速20 kn，调整螺旋桨推力使潜艇航速保持不变，航向角0°，正北航向，方向舵角为0°，首尾升降舵角设置为平衡，取首舵角为0.231 8°，尾舵角0.446 6°。其轨迹与速度曲线如图7。

(a)运动轨迹

(b)纵向速度曲线

(c)横向速度曲线

(d)垂向速度曲线

对于匀速运动的潜艇，海流的作用只是运动学上的。潜艇的初始速度是潜艇相对于海流的速度，潜艇相对固定坐标系的速度为潜艇相对于海流的速度加上海流的速度。根据海流方向及大小，流速1.20 kn，流向角45°时，可计算流速在船体纵向分量为0.849 kn，横向分量为0.849 kn。潜艇纵向速度和横向速度增量为流速的分量，由图7可看出，潜艇横向和纵向速度分量与流速在船体坐标系投影相等。对于潜艇的角速度及姿态角与静水匀速直航时相同，海流不会引起它们的变化。

在上述仿真条件下，若将方向舵角设置为20°，对潜艇在不同海流下的回转运动进行仿真，则其速度及运动轨迹如图8所示。

(a)运动轨迹

(b)纵向速度曲线

(c)横向速度曲线

(d)垂向速度曲线

潜艇做回转运动时，海流速度在艇体坐标系上的分量随航向角的变化而做周期性的变化，故潜艇的纵向速度与横向速度在速降后稳定在某一值处不断周期变化。此外，由于潜艇相对于固定坐标系具有沿海流方向的速度，大小等于流速，故潜艇在做回转运动时的回转圈会沿海流方向漂移。

4 结束语

本文结合已有的经过实际试验修正过的部分水动力系数，对潜艇在静水及海况干扰下的回转运动进行了仿真，通过定量设定仿真条件，定性分析了其运动轨迹、三维速率、角速率及三维姿态角的变化，将各种条件设置下的潜艇运动结果应用到敏感潜艇运动的导航系统的陆上性能鉴定中，从而为全面分析导航系性能提供了条件，一定程度上补充了实际海试条件的局限性。文中所用到的仿真对象某潜艇，因水动力系数有限，建立在许多假设条件上，是基于简化的潜艇六自由度运动模型进行的仿真研究；对于海浪干扰的处理，采用海浪等级表中相关海浪参数，定性描述海浪大小；对于海流干扰的处理，因为缺少数据，只从运动学上考虑其对潜艇运动造成的影响，并不适用于潜艇高速运动，有待进一步研究改进。

[1]李玉东, 胡大斌, 陈勇, 等. 潜艇操纵控制系统仿真平台设计与实现[J]. 舰船科学技术, 2013, 35(5): 123-127. LI Yudong, HU Dabin, CHEN Yong, et al. Design and realization of the submarine maneuver controlling system simulation platform[J]. Ship science and technology, 2013, 35(5): 123-127.

[2]施生达. 潜艇操纵性[M]. 北京: 国防工业出版社, 1995: 48-59. SHI Shengda. Submarine maneuverability[M]. Beijing: National Defence Industry Press, 1995: 48-59.

[3]吴秀恒, 刘祖源, 施生达, 等. 船舶操纵性[M]. 北京: 国防工业出版社, 2005: 60-62. WU Xiuheng, LIU Zuyuan, SHI Shengda, et al. Ship maneuverability[M]. Beijing: national defence industry press, 2005: 60-62.

[4]林雄伟, 胡大斌, 戴余良, 等. 潜艇操纵控制方法[J]. 舰船科学技术, 2013, 35(5): 7-10, 42. LIN Xiongwei, HU Dabin, DAI Yuliang, et al. Survey for maneuver control of submarine[J]. Ship science and technology, 2013, 35(5): 7-10, 42.

[5]时俊宇. 基于潜器空间运动的惯导系统适应性问题研究[D]. 哈尔滨: 哈尔滨工程大学, 2011: 23-26. SHI Junyu. The research of inertial navigation system adaptability based on underwater vehicle space motion[D]. Harbin: Harbin Engineering University, 2011: 23-26.

[6]朱怡. 潜艇航行训练模拟器模型简化与参数优化研究[D]. 哈尔滨: 哈尔滨工程大学, 2009: 45-47. ZHU Yi. The research on model simplification and coefficient optimization for submarine training simulator[D]. Harbin: Harbin Engineering University, 2009: 45-47.

[7]金涛, 杨枫, 王京齐, 等. 大攻角情况下破损进水潜艇水下机动性[J]. 华中科技大学学报: 自然科学版, 2008, 36(12): 79-82. JIN Tao, YANG Feng, WANG Jingqi, et al. Submerged maneuverability of flooded submarine with large angle of attack[J]. Journal of Huazhong university of science and technology: nature science edition, 2008, 36(12): 79-82.

[8]王彦磊, 袁博, 朱尚卿, 等. 海洋环境对潜艇活动的影响[J]. 舰船科学技术, 2010, 32(6): 52-55. WANG Yanlei, YUAN Bo, ZHU Shangqing, et al. The influence of marine environment on submarine activities[J]. Ship science and technology, 2010, 32(6): 52-55.

[9]OSTAFICHUK P M. AUV hydrodynamics and modelling for improved control[D]. Vancouver, Canada: The University of British Columbia, 2004: 78-81.

[10]戴余良. 潜艇在随机海浪中摇荡运动的仿真研究[J]. 计算机仿真, 2001, 18(5): 42-45. DAI Yuliang. Simulation study on swaying movement of a submarine in random ocean waves[J]. Computer simulation, 2001, 18(5): 42-45.

[11]SEO M G, KIM Y. Effects of ship motion on ship maneuvering in waves[C]//Proceedings of the 26th International Workshop on Water Waves and Floating Bodies. Athens, Greece, 2011.

[12]丰少伟, 戴余良, 陈志法. 潜艇水下回转运动稳定性分析[J]. 指挥控制与仿真, 2015, 37(1): 74-79. FENG Shaowei, DAI Yuliang, CHEN Zhifa. Analysis of motion stability of submarine s underwater rotational movement[J]. Command control & simulation, 2015, 37(1): 74-79.

Analysis of six-degree-of-freedom motion in submarines under sea disturbance

LI Hui，ZHAO Lin，MAO Ying

(College of Automation, Harbin Engineering University, Harbin 150001, China)

In order to compensate for the limitations of actual conditions in sea experiments, and to test on land the equipment that can be affected in a submarine motion, in this paper, we investigate a six-degree-of-freedom motion of submarines, including the submarine’s direct navigation and turning movements in still water and the marine environment. By analyzing the force of a moving submarine using the basic model of ocean waves and currents, we established a six-degree-of-freedom motion equation of submarines under ocean wave disturbance. Based on a submarine’s hydrodynamic coefficients, which we collected and then corrected in our experiments, we could determine the submarine’s line and angular movements, and could thus analyze its direct and turning movements. Simulation results can realistically reflect the interference in the movement of submarines in the marine environment, and thus meet the need in practical engineering to test equipment.

submarine; motion model; six-degree-of-freedom; marine environment; ocean wave disturbance; ocean current disturbance

2015-11-13.

时间：2016-12-12.

国家自然科学基金项目(61273081).

李慧(1987-), 女, 讲师,博士研究生； 赵琳(1968-), 男, 教授,博士生导师.

李慧, E-mail：lihuiheu@hotmail.com.

10.1190/jheu.1006-7043.201511027

U675.9; U666.12

A

1006-7043(2017)01-0094-07

李慧, 赵琳, 毛英.海况干扰下潜艇六自由度运动分析[J]. 哈尔滨工程大学学报， 2017, 38(1): 95-101. LI Hui，ZHAO Lin，MAO Ying.Analysis of six-degree-of-freedom motion in submarines under sea disturbance[J]. Journal of Harbin Engineering University, 2017, 38(1): 95-101.

网络出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20161212.1632.036.html