船舶摆动模型

2024.7.8
理论基础：
横摇运动响应的研究成果，是船舶摇荡运动的基础。

一、固有周期：

不同船大小对应的阻尼系数：
不同种类船的固有周期：

惯性矩公式：杜埃尔公式
$$I_{xx}=\frac{D}{12g}(B^2+4{z_g}^2)$$
D：船舶排水量,单位：吨（t）
$z_g$:以基线算起的重心高度
B:船宽
$I_{xx}\approx {I_{xx}}^{{}^{\prime }}$
由于上述式子形式简单并具有相当的精准性，为此，我国海船稳性规范在计算横摇固有周期时就采用这个公式。

我国常用的附加惯性矩：$J_{xx}=0.28I_{xx}$

固有周期：$$T_{\varphi }=\frac{2\pi }{{\omega }_{\varphi }}=2\pi \sqrt{\frac{{I_{xx}}^{{}^{\prime }}}{Dh}}(s)$$
h:初稳性高度

综合上述两式，$$T_{\varphi }=0.58\sqrt{\frac{B^2+4{z_g}^2}{h}}$$

实际运动中，阻尼会对周期和振幅产生影响，定义有阻尼的横滚周期：$${T_{\varphi }}^{{}^{\prime }}=T_{\varphi }\frac{1}{\sqrt{1-{\mu }^2}}$$

贝尔登公式：$$\mu =0.001472\frac{L{B}^4}{Dh{T_{\varphi }}^2}$$
L：船长

二、横摇角度计算
$${\theta }_a=87.5C_1{C}_2{C}_3\sqrt{0.216+\frac{Z_g}{d}}$$
d:船舶吃水深度

根据航区可以看出C1的值：

近似公式：
$$C_1=\sqrt{\frac{{\xi }_w}{\lambda }}\approx 0.32-0.0313\sqrt{T_{\varphi }}$$

C2:

C2近似为0.8

C3：

C3近似为0.9

综上，船舶横摇角公式：
$${\theta }_a=(9.36-0.9126\sqrt{T_{\varphi }})\sqrt{\frac{Z_g}{d}}$$

三、结论和规律
通过公式可以看出，船舶固有周期与振幅是成反比的，固有周期越长，横摇角越小。因此在设计船舶时将其固有周期设计的较大，以避免与波浪周期重合，发生严重的共振。
在船舶三自由度摆动设计时，效仿横滚角的设计
在船舶震荡过程中，周期逐渐变长、横滚角逐渐减小。

四、仿真验证
1.参数：
lengths = 5:1:15; % 船舶长度从5到15变化
broads = 3;
zg = 2;
h = 2;
D = 4;
d = 1.8;