2024.7.8
理论基础:
横摇运动响应的研究成果,是船舶摇荡运动的基础。
一、固有周期:
不同船大小对应的阻尼系数:
不同种类船的固有周期:
惯性矩公式:杜埃尔公式
I
x
x
=
D
12
g
(
B
2
+
4
z
g
2
)
{{I}_{xx}}=\frac{D}{12g}({{B}^{2}}+4{{z}_{g}}^{2})
Ixx=12gD(B2+4zg2)
D:船舶排水量,单位:吨(t)
z
g
{{z}_{g}}
zg:以基线算起的重心高度
B:船宽
I
x
x
≈
I
x
x
′
{{I}_{xx}}\approx {{I}_{xx}}^{'}
Ixx≈Ixx′
由于上述式子形式简单并具有相当的精准性,为此,我国海船稳性规范在计算横摇固有周期时就采用这个公式。
我国常用的附加惯性矩: J x x = 0.28 I x x {{J}_{xx}}=0.28{{I}_{xx}} Jxx=0.28Ixx
固有周期:
T
ϕ
=
2
π
ω
ϕ
=
2
π
I
x
x
′
D
h
(
s
)
{{T}_{\phi }}=\frac{2\pi }{{{\omega }_{\phi }}}=2\pi \sqrt{\frac{{{I}_{xx}}^{'}}{Dh}}(s)
Tϕ=ωϕ2π=2πDhIxx′(s)
h:初稳性高度
综合上述两式, T ϕ = 0.58 B 2 + 4 z g 2 h {{T}_{\phi }}=0.58\sqrt{\frac{{{B}^{2}}+4{{z}_{g}}^{2}}{h}} Tϕ=0.58hB2+4zg2
实际运动中,阻尼会对周期和振幅产生影响,定义有阻尼的横滚周期: T ϕ ′ = T ϕ 1 1 − μ 2 {{T}_{\phi }}^{'}={{T}_{\phi }}\frac{1}{\sqrt{1-{{\mu }^{2}}}} Tϕ′=Tϕ1−μ21
贝尔登公式:
μ
=
0.001472
L
B
4
D
h
T
ϕ
2
\mu =0.001472\frac{L{{B}^{4}}}{Dh{{T}_{\phi }}^{2}}
μ=0.001472DhTϕ2LB4
L:船长
二、横摇角度计算
θ
a
=
87.5
C
1
C
2
C
3
0.216
+
Z
g
d
{{\theta }_{a}}=87.5{{C}_{1}}{{C}_{2}}{{C}_{3}}\sqrt{0.216+\frac{{{Z}_{g}}}{d}}
θa=87.5C1C2C30.216+dZg
d:船舶吃水深度
根据航区可以看出C1的值:
近似公式:
C
1
=
ξ
w
λ
≈
0.32
−
0.0313
T
ϕ
{{C}_{1}}=\sqrt{\frac{{{\xi }_{w}}}{\lambda }}\approx 0.32-0.0313\sqrt{{{T}_{\phi }}}
C1=λξw≈0.32−0.0313Tϕ
C2:
C2近似为0.8
C3:
C3近似为0.9
综上,船舶横摇角公式:
θ
a
=
(
9.36
−
0.9126
T
ϕ
)
Z
g
d
{{\theta }_{a}}=(9.36-0.9126\sqrt{{{T}_{\phi }}})\sqrt{\frac{{{Z}_{g}}}{d}}
θa=(9.36−0.9126Tϕ)dZg
三、结论和规律
通过公式可以看出,船舶固有周期与振幅是成反比的,固有周期越长,横摇角越小。因此在设计船舶时将其固有周期设计的较大,以避免与波浪周期重合,发生严重的共振。
在船舶三自由度摆动设计时,效仿横滚角的设计
在船舶震荡过程中,周期逐渐变长、横滚角逐渐减小。
四、仿真验证
1.参数:
lengths = 5:1:15; % 船舶长度从5到15变化
broads = 3;
zg = 2;
h = 2;
D = 4;
d = 1.8;