张志敏，陈 雷

(烟台职业学院 汽车与船舶工程系，山东 烟台 264670)

0 引 言

船舶稳性设计一直是船舶行业中的重要研究方向，同时也是较复杂的研究课题。稳性设计是指船舶在装载一定载荷的条件下，能够抵抗外界复杂干扰作用力，产生船体的恢复力矩，从而实现船舶漂浮于水面并保持相对平衡的设计。船舶稳性直接决定了船舶在海浪、海风条件下的航行安全性，是检验船舶设计质量的重要指标。

近年来，随着我国对海上资源开发的不断深入，海上大型船舶的需求猛增，如大型海上风电吊装船、起重船等，这些大型船舶一方面保障了海上工程的施工，产生大量的经济效益，另一方面，大型船舶的作业环境非常复杂，船舶可靠性和安全性面临着一系列的挑战。

本文重点研究大型船舶在波浪条件下的稳定问题，建立外界干扰作用力和大型船舶的运动模型，对船舶在波浪条件下的稳性问题进行详细校核，有助于改善现有船舶的稳性设计。

1 风浪条件下的船舶稳性干扰作用力分析

大型船舶在工作过程中，海浪、海风作用力是主要的外界干扰因素，针对海浪干扰作用力进行数学建模，首先建立海浪的波面方程为：

式中：为海浪的高度（波峰到波谷的距离），w为海浪的频率值，φ为海浪的初始相位角，θ为最大波面倾角，用下式计算：

海浪的波形函数为：

其中，为周期内的海浪波数，ε为 补偿相位角

船舶在与海浪的相互作用力下，海浪的波面倾角呈周期性变化，采集0.1 min内的波面函数曲线如图1所示。

图1 一段时间内的波面函数曲线图Fig. 1 Wavefront function curve over a period of time

建立大型船舶的海浪干扰力方程为：

式中：a和a分别为船舶前进方向和船舷侧方向的加速度，为船舶质量。

建立海风的干扰作用力为；

式中：为海风的频率，()为海风的能谱密度函数，θ为来风角度。

2 大型船舶装载和风浪条件下的运动建模

大型船舶的装载和风浪条件下运动建模是进行稳性设计的关键，建立大型船舶的运动坐标系如图2所示。

可知，运动坐标系的方向指向船侧方向，指向船舶前进方向，指向地心。静止坐标系为，船舶横倾角度为φ，船舶运动坐标系和静止坐标系的转换关系为：

其矩阵形式为：

图2 风浪条件下大型船舶运动坐标系Fig. 2 Motion coordinate system of large ship under wind and wave conditions

在风浪条件下考虑船舶的纵倾角θ和航向角β，可得两坐标系的转换关系为：

在该参考运动坐标系下，建立大型船舶的运动学模型为：

式中：为船舶质量；a，a，a分别为沿3个坐标轴的加速度。

船舶在运动坐标系中的位置记（x，，），可得船舶运动方程为：

式中：α为船舶的横摇角度；β为航向角，γ为横倾角，V，V，V分别为船舶沿，，轴的速度分量。w，w，w分别为船舶绕轴、轴、轴的角速度分量。

3 大型船舶装载和风浪条件下的稳性校核分析

3.1 大型船舶随浪条件下的船舵作用力分析

船舵是船舶的稳性控制关键部件，首先建立风浪条件下船舵的力学模型，如图3所示。

图中，α为船舵的攻角，为海浪作用在船舵表面的合力，可分解为

图3 风浪条件下船舵的力学模型Fig. 3 Mechanical model of rudder under wind and wave conditions

为船舵剖面的弦长，定义船舵的压力中心为x，分散系数按下式计算：

在建立船舵作用力前，根据流体动力学理论，首先建立基本控制方程分别如下：

1）连续性方程

式中：Δ为拉普拉斯算式，ρ为液体密度，v为相对运动速度。

2）动量方程

式中：为海水的体积力，σ为海水的体积域函数。

3）能量方程

式中：为海水内能，为海水的动能。

建立船舵的作用力模型为：

式中： ， ，分别为作用力沿3个坐标方向的切应力。

3.2 大型船舶随浪条件下的恢复力矩计算

船舶正常航行状态下处于相对平衡状态，当外界风浪干扰作用力突然施加时，会对船舶的稳定性产生破坏作用，使船舶产生明显的横摇角度，影响船体水线面的变化和船体浮态，进一步影响船体的稳定性。

当船舶发生大倾角稳性问题时，分析船舶倾斜后产生的复原力矩与干扰力的相互作用是关键，当船舶处于静水之中时水线面为一水平面，当在风浪条件时船体的水线面产生倾斜，建立该过程的数学模型如图4所示。

图4 船舶倾角稳性过程的数学模型Fig. 4 Mathematical model of ship inclination stability process

可知，在风浪作用下船舶的水线面产生倾角为θ，海浪作用在船体的作用力如下式：

式中：ρ为 海水的密度，A为船舶在海浪方向的投影面积，为升力系数，为航行速度。

海浪作用在船体的弯矩如下式：

式中：L为船体纵向长度的投影长度；为船体宽度；δ为船舶载重系数， δ≈0.5 ，和分别为船舶的弯矩系数，与船舶长度、吃水深度等因素有关，按下式计算：

船舶的恢复力矩由船体重心和船舵作用力共同产生，如下式：

式中，为船舶中心到剖面中线的距离。

对比不同吃水深度下的船体稳性倾角变化数据，如图5所示。

图5 不同吃水深度下的船体稳性倾角变化数据Fig. 5 Hull stability inclination variation data at different draft depths

4 结 语

船舶的稳性设计和分析是提高船舶航行安全性的关键，也是船舶设计和测试过程中需要重点考虑的因素。本文针对大型船舶在装载和风浪条件下的稳定性问题进行详细建模和研究，建立了干扰作用力模型、船舶运动模型、船舵作用力模型和恢复力矩模型，有助于提高船舶的稳性设计。