-
船舶稳性是船舶航行安全的重要保障,一直以来受到国际海事组织(IMO)的高度关注。IMO基于大量的事故统计资料和静力学理论,制定了《2008年国际完整稳性规则》[1]。然而,许多满足现行完整稳性衡准的船舶在航行时仍然发生了事故,为此,IMO又有针对性地开展了第二代完整稳性衡准的制定工作。作为第二代完整稳性的5种失效模式之一,参数横摇衡准包含3层:第1层和第2层为薄弱性衡准,第3层为直接计算[2- 3]。在2019年2月SDC第6次会议上,IMO通过投票初步确定了参数横摇薄弱性衡准规则,后续的工作除了进行大量样船计算并验证第1层和第2层的一致性外,还需要研究船型对参数横摇薄弱性的影响,以期在新一代完整稳性规范实施后为船舶设计提供指导。
目前,针对船型对参数横摇薄弱性影响的研究主要集中在船舶主尺度和船型系数等方面,关于船体型线对参数横摇薄弱性影响的研究还很少。马坤等[2]选取6艘母型船进行了大量的船型变换,研究了船舶主尺度对参数横摇薄弱性的影响;付超等[4]针对C11集装箱船的敏感性因素进行了分析,并以散货船型线为参照对其型线进行了修改,使其第2层计算结果明显下降。参数横摇的发生机理与船在波浪中复原力臂的变化密切相关,而船体型线则直接影响到船舶的复原力臂曲线,因此有必要研究船体型线对参数横摇薄弱性的影响。
本文拟选取C11集装箱船、大连400 t渔政船和秋刀鱼渔船作为母型船,首先应用自由曲面变换(free-form deformation,FFD)法对船艏UV度进行变换,然后利用本课题组前期开发的参数横摇薄弱性衡准评估软件对母型船和变换型线后的系列船型进行第1,2层薄弱性计算(该程序的可靠性已得到验证[5]),最后通过对比计算结果得出船体型线对参数横摇薄弱性的影响规律。
-
参数横摇薄弱性衡准评估分为2层,其中第1层评估通过计算纵向规则波经过船体时,船舶初稳性高的变化幅值来判断是否满足衡准要求:
$$\frac{{\Delta GM}}{{GM}} < {R_{{\rm{PR}}}}$$ (1) 式中:
$\Delta GM$ 为纵向规则波通过船体时初稳性高变化幅值,该纵向规则波的波长为船长,波陡为0.0167;$GM$ 为船舶在计算载况下,在静水中考虑自由液面修正后的初稳性高;${R_{{\rm{PR}}}}$ 为衡准值,其与船是否为尖舭船,以及舭龙骨宽度和其他船型要素相关。当船舶不满足第1层衡准时,需进行第2层校核。参数横摇薄弱性衡准的第2层校核分C1和C2两种不同的计算方法。其中C1用于计算船舶在给定的16种波浪条件下的失效概率,C2用于计算船舶在北大西洋海浪谱(或船舶航行海域的海浪谱)所给定的海况下,以不同航速和航向角航行时的失效概率。计算公式[6-9]为
$$C1 = \sum\limits_{i = 1}^{16} {{W_i}{C_i}} $$ (2) $$C2 = \left[ {\sum\limits_{i = 1}^{12} {C2(Fn{ _i},{\beta _{\rm {h}}})} + C2(0,{\beta _h}) + \sum\limits_{i = 1}^{12} {C2(Fn{_i},{\beta _{\rm {f}}})} } \right]/25$$ (3) 式中:
${W_i}$ 为权重值;${C_i}$ 为0或1,通过计算认为在某海况下容易发生参数横时摇取1,否则为0;$Fn{_i}$ 为计算航速的弗劳德数;${\beta _{\rm {h}}}$ ,${\beta _{\rm {f}}}$ 分别代表迎浪和随浪。当满足
$C1 < 0.06$ 或$C2 < 0.025$ 的条件时,认为通过了第2层校核。具体计算方法详见参考文献[6]。C2方法因考虑海况数更多,且需要通过求解横摇运动方程来计算横摇幅值,故其计算更为精确,因此本文在进行参数横摇薄弱性衡准第2层校核时仅采用C2方法。 -
选取C11集装箱船、大连400 t渔政船和秋刀鱼渔船作为母型船进行船型变换,其中C11集装箱船在实际航行中发生过参数横摇。3艘母型船的主要船型参数如表1所示,3艘母型船的型线图如图1~图3所示。
表 1 母型船主要船型参数
Table 1. Main parameters of parent ships
参数 C11集装箱船 大连400 t渔政船 秋刀鱼渔船 垂线间长/m 262 48.6 65.6 船宽/m 40 8 11.1 型深/m 24.45 4.3 7 吃水/m 11.5 2.85 4.447 排水体积/m3 67384 543 2293 初稳性高/m 1.902 1.862 0.548 中横剖面系数 0.958 0.752 0.978 方形系数 0.559 0.516 0.695 服务航速/kn 25 16 14 -
本文采用FFD法实现母型船的型线变换。FFD法是船型设计优化过程中常用的方法,该方法可应用于任何曲面,不受曲面表示方法的限制,还可对物体的局部或整体进行变形,变形后的曲面能够保持几何连续性。
FFD法的基本步骤如下:
1) 选定一个长方体形状的控制体包围待变形曲面;
2) 在控制体内选择一系列控制点;
3) 调节控制点位置改变控制体形状,将控制体的变形传递给待变形曲面,实现曲面变形[10-11]。
图4所示为FFD法控制体及坐标系示意图。其中,O-XYZ为全局坐标系,O′-RST为控制体坐标系,控制体内任意一点的全局坐标与控制点的全局坐标满足伯恩斯坦(Bernstein)多项式:
$$X(r,s,t) = \sum\limits_{i = 0}^l {\sum\limits_{j = 0}^m {\sum\limits_{k = 0}^n {{B_{i,l}}(r){B_{j,m}}(s){B_{k,n}}(t)} } } {Q_{i,j,k}}$$ (4) 式中:l, m, n分别为控制体在R轴,S轴和T轴方向的控制点个数;r, s, t为控制体内任意点在控制体坐标系下的坐标,满足
$0 \leqslant r \leqslant 1$ ,$0 \leqslant s \leqslant 1$ ,$0 \leqslant t \leqslant 1$ ;${Q_{i,j,k}}$ 为控制点坐标;B为伯恩斯坦多项式基函数[11]。 -
使用控制体包围船体舯部至艏部(前半体),将水线以下部分控制点沿船宽方向向外侧拉伸,水线以上部分控制点则沿船宽方向向内压缩,使母型船船艏由V 型变化为U型,由此产生的新船型命名为5号船。假设
5号船第i个控制点的变化量为 ${d_i}$ (i=1, 2,···, N,其中N为控制点总数)。将控制点变化量依次设定为$0.2{d_i}$ ,$0.4{d_i}$ ,$0.6{d_i}$ ,$0.8{d_i}$ ,从而生成4个新船型,分别命名为1,2,3,4号船,母型船编号为0号。应用上述方法变换3艘母型船艏部UV度,在变换时,只保证主尺度不变,并未约束排水量,各种变形均保证线条的光顺性,生成3个系列样船。即每个系列共6艘样船,由0号到5号船,其中母型船编号为0;船艏型线由V型向U型逐渐变化。系列船型中0号和5号船的型线图如图5~图7所示,系列船型的主要船型参数如表2~表4所示。图 6 大连400 t 渔政船系列0号和5号样船型线图
Figure 6. No.0 and No.5 Ships' lines generated from Dalian 400 t fishery administration ship
表 2 C11集装箱船系列样船主要船型参数
Table 2. Main parameters of series ships generated from C11 container ship
参数 0号船 1号船 2号船 3号船 4号船 5号船 垂线间长/m 262 262 262 262 262 262 船宽/m 40 40 40 40 40 40 型深/m 24.45 24.45 24.45 24.45 24.45 24.45 吃水/m 11.5 11.5 11.5 11.5 11.5 11.5 排水体积/m3 67 384 67918 68 452 68 986 69 517 70 047 初稳性高/m 1.902 1.793 1.685 1.581 1.482 1.386 前半体方形系数 0.505 0.512 0.519 0.527 0.534 0.542 水线面系数 0.774 0.775 0.776 0.777 0.778 0.779 表 3 大连400 t渔政船系列样船主要船型参数
Table 3. Main parameters of series ships generated from Dalian 400 t fishery administration ship
参数 0号船 1号船 2号船 3号船 4号船 5号船 垂线间长/m 48.6 48.6 48.6 48.6 48.6 48.6 船宽/m 8 8 8 8 8 8 型深/m 4.3 4.3 4.3 4.3 4.3 4.3 吃水/m 2.85 2.85 2.85 2.85 2.85 2.85 排水体积/m3 543 545 547 549 551 553 初稳性高/m 1.862 1.817 1.772 1.729 1.688 1.648 前半体方形系数 0.460 0.463 0.465 0.467 0.470 0.472 水线面系数 0.830 0.825 0.820 0.815 0.810 0.805 表 4 秋刀鱼渔船系列样船主要船型参数
Table 4. Main parameters of series ships generated from a fishing ship
参数 0号船 1号船 2号船 3号船 4号船 5号船 垂线间长/m 65.6 65.6 65.6 65.6 65.6 65.6 船宽/m 11.1 11.1 11.1 11.1 11.1 11.1 型深/m 7 7 7 7 7 7 吃水/m 4.447 4.447 4.447 4.447 4.447 4.447 排水体积/m3 2293 2294 2294 2295 2296 2297 初稳性高/m 0.548 0.530 0.512 0.494 0.477 0.460 前半体方形系数 0.675 0.675 0.675 0.675 0.675 0.675 水线面系数 0.908 0.905 0.903 0.900 0.898 0.895 -
研究船体艏部UV度对参数横摇薄弱性第1层校核(以下简称“第1层”)结果的影响。选取3个系列共18艘样船作为研究对象,采用本课题组前期开发的参数横摇薄弱性校核程序进行计算。计算时,假设同一系列样船的重心、吃水等参数相同,且未设置舭龙骨。具体计算结果如表5~表7所示。
表 5 C11集装箱船系列样船第1层计算结果
Table 5. Calculation results of level 1 of series ships generated from C11 container ship
船型编号 $\Delta GM$/m $GM$/m $\Delta GM/GM$ 衡准值 0 1.559 1.902 0.819 0.17 1 1.500 1.793 0.837 0.17 2 1.442 1.685 0.856 0.17 3 1.386 1.581 0.876 0.17 4 1.331 1.482 0.898 0.17 5 1.276 1.386 0.921 0.17 表 6 大连400 t 渔政船系列样船第1层计算结果
Table 6. Calculation results of level 1 of series ships generated from Dalian 400 t fishery administration ship
船型编号 $\Delta GM$/m $GM$/m $\Delta GM/GM$ 衡准值 0 0.042 1.862 0.022 0.17 1 0.037 1.817 0.020 0.17 2 0.032 1.772 0.019 0.17 3 0.027 1.730 0.016 0.17 4 0.023 1.688 0.014 0.17 5 0.019 1.648 0.012 0.17 表 7 秋刀鱼渔船系列样船第1层计算结果
Table 7. Calculation results of level 1 of series ships generated from Saury fishing ship
船型编号 $\Delta GM$/m $GM$/m $\Delta GM/GM$ 衡准值 0 0.120 0.548 0.218 0.17 1 0.117 0.530 0.221 0.17 2 0.115 0.512 0.224 0.17 3 0.112 0.494 0.228 0.17 4 0.110 0.477 0.231 0.17 5 0.108 0.460 0.235 0.17 由表5~表7可知,C11集装箱船系列样船和秋刀鱼渔船系列样船未通过参数横摇薄弱性第1层衡准校核,大连400 t渔政船系列船型通过了校核。3个系列样船,随着艏部UV度由V型向U型变化,
$\Delta GM$ 与GM同时下降。图8为3个系列样船的第1层计算结果。由图可知,随着艏部型线由V型向U型变化,$\Delta GM/GM$ 的变化规律并不一致,C11集装箱船和秋刀鱼渔船是上升的,而大连400 t渔政船则是降低的,因为型线变化后$\Delta GM$ 与GM同时变化,$\Delta GM/GM$ 的变化趋势与$\Delta GM$ ,GM的变化程度有关。第1层评估只是简单的经验公式,仅反映了船舶静水中初稳性高和波浪中初稳性高的变化,没有反映阻尼和遭遇频率因素,所以无法从$\Delta GM/GM$ 的角度判断船艏UV度对参数横摇薄弱性的影响大小。 -
研究船体艏部UV度对系列船型参数横摇薄弱性第2层(以下简称“第2层”)C2计算结果的影响。选取3个系列共18艘样船作为研究对象,采用本课题组前期开发的参数横摇薄弱性校核程序进行计算。在计算时,假设同一系列样船的航速、重心、吃水等参数相同,计算航速为25个,阻尼系数采用池田法进行计算,自摇周期应用《2008年国际完整稳性规则》中的公式计算。若考虑舭龙骨,大连400 t渔政船与秋刀鱼渔船系列样船第2层C2的计算结果将均为0,为了观察C2的变化趋势,未设置舭龙骨。3个系列样船第2层C2的计算结果如表8所示。
表 8 3个系列样船第2层C2计算结果
Table 8. Calcution results of C2 of three series ships
船型编号 C11
集装箱船大连400 t
渔政船秋刀鱼
渔船衡准值 0 0.0857 1.20×10−6 2.38×10−3 0.025 1 0.0801 1.60×10−6 2.45×10−3 0.025 2 0.0750 1.60×10−6 9.86×10−4 0.025 3 0.0728 3.36×10−6 7.21×10−4 0.025 4 0.0719 4.72×10−6 7.00×10−4 0.025 5 0.0649 6.48×10−6 6.81×10−4 0.025 校核结果表明,C11集装箱船系列样船未能通过参数横摇薄弱性第2层衡准校核,大连400 t渔政船和秋刀鱼渔船通过了校核。随着船舶艏部UV度由V型向U型变化,C11集装箱船与秋刀鱼渔船的参数横摇薄弱性C2整体呈下降趋势,大连400 t渔政船的参数横摇薄弱性C2小幅上升,如图9所示。
由表8可见,大连400 t渔政船系列船型和秋刀鱼渔船系列船型的0号到1号船第2层C2的计算结果略有增加。将以大连400 t渔政船为例分析C2计算结果略微上升的原因。表9所示为大连400 t渔政船迎浪航行时不同航速、不同波陡条件下横摇幅值的计算结果。表中,横摇幅值的单位为(°),波陡航速序号0,1,2,3分别代表服务航速为0,
${V_{\rm{S}}}$ ,${{\sqrt 3 }}{V_{\rm{S}}} /{2}$ ,${V_{\rm{S}}}/{2}$ 。由表9可知,在航速为${{\sqrt 3 }}{V_{\rm{S}}}/{2}$ ,波陡为0.1时,随着船艏型线由V型向U型变化,横摇幅值逐渐变小;在航速为${{\sqrt 3 }}{V_{\rm{S}}}/{2}$ ,波陡为0.09以及航速为${V_{\rm{S}}}/{2}$ ,波陡为0.1时,由2号船到3号船横摇幅值由0°分别增加到38.5°和44°,使得接下来插值过程中部分海况的横摇幅值大于25°,最终导致3号船第2层C2的校核结果比2号船大。这2种情况下横摇幅值突变的原因是0~2号船满足式(5),按照草案的规定不易发生参数横摇,横摇幅值取0°,而3~5号船因不满足该条件,故需要计算横摇幅值。表 9 大连400 t渔政船迎浪航行不同航速和不同波陡下横摇幅值
Table 9. Roll amplitude of Dalian 400 t fishery administration ship heading with the sea at different speeds and wave slops
船型编号 波陡航速序号 波陡 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 38.7 2 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 37.8 2 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 36.7 2 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 38.5 35.6 2 0 0 0 0 0 0 0 0 0 44 3 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 38.5 35.8 2 0 0 0 0 0 0 0 0 0 43.5 3 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 37.7 34.7 2 0 0 0 0 0 0 0 0 0 42.8 3 0 0 0 0 0 0 0 0 0 0 $$\frac{{G{M_{{\rm{amp}}}}}}{{GM}} < 2\sqrt {{{\left\{ {\left( {1 + \frac{{G{M_{{\rm{mean}}}}}}{{GM}}} \right) - \frac{1}{4}{{\left( {\frac{{{\omega _{\rm{e}}}}}{{{\omega _\phi }}}} \right)}^2}} \right\}}^2} + {{\left( {\frac{\alpha }{{{\omega _\phi }}}} \right)}^2}{{\left( {\frac{{{\omega _{\rm{e}}}}}{{{\omega _\phi }}}} \right)}^2}} $$ (6) 式中:
$G{M_{{\rm{amp}}}}$ 为波浪中初稳性高变化幅值,m;$GM$ 为静水中初稳性高,m;$G{M_{{\rm{mean}}}}$ 为波浪中初稳性高变化的平均值,m;${\omega _{\rm{e}}}$ 为遭遇频率;${\omega _\phi }$ 为自摇频率;$\alpha $ 为一次横摇阻尼系数。若不考虑该条件,对0到5号船应用平均法计算横摇幅值,由表10~表11可知,2种计算条件下,0~2号船 的横摇幅值小于25°且逐渐减小,3号船的横摇幅值产生了突变。3号船产生突变的原因为:从0号到5号船,船舶自摇频率
${\omega _\phi }$ 逐渐接近遭遇频率${\omega _{\rm{e}}}$ 的一半,3号船进入参数横摇激励频率范围,因此使得横摇幅值产生突变,最终导致C2计算结果略有增加。表 10 大连400 t渔政船迎浪航行时的横摇幅值计算结果(航速
${{\sqrt 3 }}{V_S}/{2}$ ,波陡0.09)Table 10. Roll amplitude of Dalian 400 t fishery administration heading with the sea at the speed of
${{\sqrt 3 }}{V_{\rm{S}}}/{2}$ and wave slop is 0.09船型编号 ${\omega _{\rm{e}}}$ ${\omega _\phi }$ 横摇幅值/(°) 0 2.179 1.288 20.2 1 2.179 1.272 16.4 2 2.179 1.257 11.1 3 2.179 1.241 38.5 4 2.179 1.226 38.5 5 2.179 1.212 37.7 表 11 大连400 t渔政船迎浪航行时的横摇幅值计算结果(航速
${V_S}/{2}$ ,波陡0.1)Table 11. Roll amplitude of Dalian 400 t fishery administration heading with the sea at the speed of
${V_{\rm{S}}}/{2}$ and wave slop is 0.1船型编号 ${\omega _{\rm{e}}}$ ${\omega _\phi }$ 横摇幅值/(°) 0 2.038 1.288 17.1 1 2.038 1.272 12.7 2 2.038 1.257 4.9 3 2.038 1.241 44.0 4 2.038 1.226 43.5 5 2.038 1.212 42.8 按式(6)计算3艘样船的横摇阻尼力矩。
$$R( {\dot \phi } ) = ( {{I_{xx}} + {J_{xx}}} )\{ {2\alpha \dot \phi + \gamma {{\dot \phi }^3}} \}$$ (7) 式中:
$R$ 为横摇阻尼力矩,N·m;${I_{xx}}$ ,${J_{xx}}$ 分别为横摇惯性矩与附加惯性矩,m4;$\alpha $ ,$\gamma $ 分别为线性阻尼系数和三次阻尼系数;$\dot \phi $ 为横摇角速度,rad/s。航速为服务航速,3个系列船型不同横摇角速度下的横摇阻尼力矩(单位:N·m)计算结果如表12~表14所示。从表中可看出,随着船艏型线由V型向U型变换,C11集装箱船的横摇阻尼力矩增加了;秋刀鱼渔船在横摇角速度大于0.14 rad/s时,阻尼力矩增加;而大连400 t渔政船的横摇阻尼力矩减小,该系列船型C2的变化趋势与其他系列船型不同。
表 12 C11集装箱船横摇阻尼力矩
Table 12. Roll damping moment of series ships generated from C11 container ship
横摇角速度/(rad·s-1) 横摇阻尼/(N·m) 0号船 1号船 2号船 3号船 4号船 5号船 0.08 2.614×107 2.623×107 2.643×107 2.670×107 2.705×107 2.749×107 0.10 4.880×107 4.923×107 4.982×107 5.054×107 5.137×107 5.235×107 0.12 8.220×107 8.318×107 8.440×107 8.581×107 8.739×107 8.920×107 0.14 1.285×108 1.303×108 1.324×108 1.348×108 1.375×108 1.404×108 0.16 1.899×108 1.927×108 1.961×108 1.998×108 2.039×108 2.085×108 0.18 2.684×108 2.727×108 2.777×108 2.831×108 2.891×108 2.957×108 0.20 3.663×108 3.724×108 3.794×108 3.871×108 3.953×108 4.045×108 0.22 4.857×108 4.940×108 5.035×108 5.138×108 5.250×108 5.373×108 0.24 6.288×108 6.398×108 6.522×108 6.658×108 6.804×108 6.964×108 0.26 7.976×108 8.118×108 8.278×108 8.452×108 8.638×108 8.844×108 表 13 大连400 t渔政船横摇阻尼力矩
Table 13. Roll damping moment of series ships generated from Dalian 400 t fishery administration ship
横摇角速度/(rad·s-1) 横摇阻尼/(N·m) 0号船 1号船 2号船 3号船 4号船 5号船 0.08 4.481×104 4.370×104 4.258×104 4.145×104 4.033×104 3.921×104 0.10 5.636×104 5.497×104 5.358×104 5.217×104 5.078×104 4.938×104 0.12 6.814×104 6.648×104 6.481×104 6.314×104 6.147×104 5.979×104 0.14 8.020×104 7.827×104 7.633×104 7.438×104 7.245×104 7.050×104 0.16 9.258×104 9.039×104 8.818×104 8.597×104 8.377×104 8.156×104 0.18 1.053×105 1.029×105 1.004×105 9.794×104 9.547×104 9.300×104 0.20 1.185×105 1.158×105 1.131×105 1.103×105 1.076×105 1.049×105 0.22 1.321×105 1.292×105 1.262×105 1.232×105 1.202×105 1.173×105 0.24 1.463×105 1.431×105 1.398×105 1.366×105 1.334×105 1.302×105 0.26 1.609×105 1.575×105 1.541×105 1.506×105 1.471×105 1.437×105 表 14 秋刀鱼渔船横摇阻尼力矩
Table 14. Roll damping moment of series ships generated from a Saury fishing ship
横摇角速度/(rad·s-1) 横摇阻尼/(N·m) 0号船 1号船 2号船 3号船 4号船 5号船 0.08 2.126×104 2.094×104 2.066×104 2.042×104 2.022×104 2.006×104 0.10 3.567×104 3.541×104 3.519×104 3.504×104 3.494×104 3.489×104 0.12 5.616×104 5.603×104 5.598×104 5.600×104 5.610×104 5.627×104 0.14 8.392×104 8.405×104 8.427×104 8.459×104 8.500×104 8.551×104 0.16 1.202×105 1.207×105 1.213×105 1.220×104 1.229×105 1.239×105 0.18 1.662×105 1.672×105 1.683×105 1.697×105 1.711×105 1.728×105 0.20 2.230×105 2.247×105 2.266×105 2.287×105 2.310×105 2.334×105 0.22 2.921×105 2.946×105 2.974×105 3.004×105 3.037×105 3.072×105 0.24 3.744×105 3.780×105 3.819×105 3.861×105 3.906×105 3.954×105 0.26 4.714×105 4.763×105 4.815×105 4.870×105 4.929×105 4.992×105 -
本文选取C11集装箱船、大连400 t渔政船、秋刀鱼渔船为母型船,应用FFD法对3艘母型船艏部UV度进行变换,生成了3个系列共15艘新船型,并对3个系列共18艘样船进行了参数横摇薄弱性校核,探究了船体型线对参数横摇薄弱性的影响规律,可为参数横摇薄弱性衡准在国内的推广应用提供技术依据。通过本文的研究,主要得到以下结论:
1) 对于参数横摇薄弱性第1层衡准,船舶艏部型线由V型向U型变化时,船舶在静水中的初稳性高减小,波浪中初稳性高的变化幅值逐渐减小,3个系列样船第1层
$\Delta GM/GM$ 的变化并不一致,且仅从$\Delta GM/GM$ 不能充分反映船体形状对参数横摇薄弱性的影响规律。2) 对于参数横摇薄弱性第2层衡准,随着船舶艏部型线由V型向U型变化,C11集装箱船与秋刀鱼渔船的参数横摇薄弱性第2层C2整体呈现下降的趋势,船舶参数横摇失效概率降低,船舶越不容易发生参数横摇,而大连400 t渔政船的C2则有小幅的上升。
3) 大连400 t渔政船的C2有小幅的上升是因为自摇频率逐渐接近于遭遇频率的一半,进入了参数横摇的激励频率范围,在某些航速与波陡组合下,横摇幅值还出现了突然上升情况,同时随着型线由V型变换为U型,大连400 t渔政船的横摇阻尼力矩减小,最终导致第2层C2的计算结果上升。
Influence of U or V degree of stem on vulnerability of parametric rolling
-
摘要:
目的 以船艏UV度为对象探究其对参数横摇薄弱性的影响。 方法 选取集装箱船、渔政船和渔船这3种不同类型的船舶为母型船,应用自由曲面变换(free-form deformation,FFD)法,改变母型船艏部UV度,生成3个系列共15艘新船型,采用本课题组前期开发的参数横摇薄弱性校核程序对该系列船型进行参数横摇薄弱性计算,并对计算结果进行分析。 结果 结果显示,船艏型线由V型向U型变化时船舶参数横摇薄弱性变好。 结论 所做研究可为参数横摇薄弱性衡准在国内的推广应用提供技术依据,也可为船体型线设计提供参考。 Abstract:Objectives This paper selects the U or V degree of ships' stems as the research target for identifying its influence on parametric rolling vulnerability. Methods Three kinds of ships were chosen as the parent ships: a container ship, fishing vessel and fishery administration ship. Three series of ships were generated by changing the U or V degree of these parent ships' stems using the free-form deformation method. Each series contained five newly created ships for a total of 15 new ships. The parametric rolling vulnerability of the ships was assessed and the rule of the influence of the U or V degree of the stem on vulnerability to parametric rolling was obtained from the assessment results. Results The assessment results show that changing the form of the stem from a V section to a U section has an advantageous impact on vulnerability to parametric rolling. Conclusions The research in this paper provides a technical basis for the application of parametric rolling vulnerability criterion in China, and can provide valuable references for hull form design. -
Key words:
- parametric rolling /
- vulnerability criteria /
- U or V degree /
- free-form deformation (FFD)
-
表 1 母型船主要船型参数
Table 1. Main parameters of parent ships
参数 C11集装箱船 大连400 t渔政船 秋刀鱼渔船 垂线间长/m 262 48.6 65.6 船宽/m 40 8 11.1 型深/m 24.45 4.3 7 吃水/m 11.5 2.85 4.447 排水体积/m3 67384 543 2293 初稳性高/m 1.902 1.862 0.548 中横剖面系数 0.958 0.752 0.978 方形系数 0.559 0.516 0.695 服务航速/kn 25 16 14 表 2 C11集装箱船系列样船主要船型参数
Table 2. Main parameters of series ships generated from C11 container ship
参数 0号船 1号船 2号船 3号船 4号船 5号船 垂线间长/m 262 262 262 262 262 262 船宽/m 40 40 40 40 40 40 型深/m 24.45 24.45 24.45 24.45 24.45 24.45 吃水/m 11.5 11.5 11.5 11.5 11.5 11.5 排水体积/m3 67 384 67918 68 452 68 986 69 517 70 047 初稳性高/m 1.902 1.793 1.685 1.581 1.482 1.386 前半体方形系数 0.505 0.512 0.519 0.527 0.534 0.542 水线面系数 0.774 0.775 0.776 0.777 0.778 0.779 表 3 大连400 t渔政船系列样船主要船型参数
Table 3. Main parameters of series ships generated from Dalian 400 t fishery administration ship
参数 0号船 1号船 2号船 3号船 4号船 5号船 垂线间长/m 48.6 48.6 48.6 48.6 48.6 48.6 船宽/m 8 8 8 8 8 8 型深/m 4.3 4.3 4.3 4.3 4.3 4.3 吃水/m 2.85 2.85 2.85 2.85 2.85 2.85 排水体积/m3 543 545 547 549 551 553 初稳性高/m 1.862 1.817 1.772 1.729 1.688 1.648 前半体方形系数 0.460 0.463 0.465 0.467 0.470 0.472 水线面系数 0.830 0.825 0.820 0.815 0.810 0.805 表 4 秋刀鱼渔船系列样船主要船型参数
Table 4. Main parameters of series ships generated from a fishing ship
参数 0号船 1号船 2号船 3号船 4号船 5号船 垂线间长/m 65.6 65.6 65.6 65.6 65.6 65.6 船宽/m 11.1 11.1 11.1 11.1 11.1 11.1 型深/m 7 7 7 7 7 7 吃水/m 4.447 4.447 4.447 4.447 4.447 4.447 排水体积/m3 2293 2294 2294 2295 2296 2297 初稳性高/m 0.548 0.530 0.512 0.494 0.477 0.460 前半体方形系数 0.675 0.675 0.675 0.675 0.675 0.675 水线面系数 0.908 0.905 0.903 0.900 0.898 0.895 表 5 C11集装箱船系列样船第1层计算结果
Table 5. Calculation results of level 1 of series ships generated from C11 container ship
船型编号 $\Delta GM$ /m$GM$ /m$\Delta GM/GM$ 衡准值 0 1.559 1.902 0.819 0.17 1 1.500 1.793 0.837 0.17 2 1.442 1.685 0.856 0.17 3 1.386 1.581 0.876 0.17 4 1.331 1.482 0.898 0.17 5 1.276 1.386 0.921 0.17 表 6 大连400 t 渔政船系列样船第1层计算结果
Table 6. Calculation results of level 1 of series ships generated from Dalian 400 t fishery administration ship
船型编号 $\Delta GM$ /m$GM$ /m$\Delta GM/GM$ 衡准值 0 0.042 1.862 0.022 0.17 1 0.037 1.817 0.020 0.17 2 0.032 1.772 0.019 0.17 3 0.027 1.730 0.016 0.17 4 0.023 1.688 0.014 0.17 5 0.019 1.648 0.012 0.17 表 7 秋刀鱼渔船系列样船第1层计算结果
Table 7. Calculation results of level 1 of series ships generated from Saury fishing ship
船型编号 $\Delta GM$ /m$GM$ /m$\Delta GM/GM$ 衡准值 0 0.120 0.548 0.218 0.17 1 0.117 0.530 0.221 0.17 2 0.115 0.512 0.224 0.17 3 0.112 0.494 0.228 0.17 4 0.110 0.477 0.231 0.17 5 0.108 0.460 0.235 0.17 表 8 3个系列样船第2层C2计算结果
Table 8. Calcution results of C2 of three series ships
船型编号 C11
集装箱船大连400 t
渔政船秋刀鱼
渔船衡准值 0 0.0857 1.20×10−6 2.38×10−3 0.025 1 0.0801 1.60×10−6 2.45×10−3 0.025 2 0.0750 1.60×10−6 9.86×10−4 0.025 3 0.0728 3.36×10−6 7.21×10−4 0.025 4 0.0719 4.72×10−6 7.00×10−4 0.025 5 0.0649 6.48×10−6 6.81×10−4 0.025 表 9 大连400 t渔政船迎浪航行不同航速和不同波陡下横摇幅值
Table 9. Roll amplitude of Dalian 400 t fishery administration ship heading with the sea at different speeds and wave slops
船型编号 波陡航速序号 波陡 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 38.7 2 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 37.8 2 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 36.7 2 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 38.5 35.6 2 0 0 0 0 0 0 0 0 0 44 3 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 38.5 35.8 2 0 0 0 0 0 0 0 0 0 43.5 3 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 37.7 34.7 2 0 0 0 0 0 0 0 0 0 42.8 3 0 0 0 0 0 0 0 0 0 0 表 10 大连400 t渔政船迎浪航行时的横摇幅值计算结果(航速
${{\sqrt 3 }}{V_S}/{2}$ ,波陡0.09)Table 10. Roll amplitude of Dalian 400 t fishery administration heading with the sea at the speed of
${{\sqrt 3 }}{V_{\rm{S}}}/{2}$ and wave slop is 0.09船型编号 ${\omega _{\rm{e}}}$ ${\omega _\phi }$ 横摇幅值/(°) 0 2.179 1.288 20.2 1 2.179 1.272 16.4 2 2.179 1.257 11.1 3 2.179 1.241 38.5 4 2.179 1.226 38.5 5 2.179 1.212 37.7 表 11 大连400 t渔政船迎浪航行时的横摇幅值计算结果(航速
${V_S}/{2}$ ,波陡0.1)Table 11. Roll amplitude of Dalian 400 t fishery administration heading with the sea at the speed of
${V_{\rm{S}}}/{2}$ and wave slop is 0.1船型编号 ${\omega _{\rm{e}}}$ ${\omega _\phi }$ 横摇幅值/(°) 0 2.038 1.288 17.1 1 2.038 1.272 12.7 2 2.038 1.257 4.9 3 2.038 1.241 44.0 4 2.038 1.226 43.5 5 2.038 1.212 42.8 表 12 C11集装箱船横摇阻尼力矩
Table 12. Roll damping moment of series ships generated from C11 container ship
横摇角速度/(rad·s-1) 横摇阻尼/(N·m) 0号船 1号船 2号船 3号船 4号船 5号船 0.08 2.614×107 2.623×107 2.643×107 2.670×107 2.705×107 2.749×107 0.10 4.880×107 4.923×107 4.982×107 5.054×107 5.137×107 5.235×107 0.12 8.220×107 8.318×107 8.440×107 8.581×107 8.739×107 8.920×107 0.14 1.285×108 1.303×108 1.324×108 1.348×108 1.375×108 1.404×108 0.16 1.899×108 1.927×108 1.961×108 1.998×108 2.039×108 2.085×108 0.18 2.684×108 2.727×108 2.777×108 2.831×108 2.891×108 2.957×108 0.20 3.663×108 3.724×108 3.794×108 3.871×108 3.953×108 4.045×108 0.22 4.857×108 4.940×108 5.035×108 5.138×108 5.250×108 5.373×108 0.24 6.288×108 6.398×108 6.522×108 6.658×108 6.804×108 6.964×108 0.26 7.976×108 8.118×108 8.278×108 8.452×108 8.638×108 8.844×108 表 13 大连400 t渔政船横摇阻尼力矩
Table 13. Roll damping moment of series ships generated from Dalian 400 t fishery administration ship
横摇角速度/(rad·s-1) 横摇阻尼/(N·m) 0号船 1号船 2号船 3号船 4号船 5号船 0.08 4.481×104 4.370×104 4.258×104 4.145×104 4.033×104 3.921×104 0.10 5.636×104 5.497×104 5.358×104 5.217×104 5.078×104 4.938×104 0.12 6.814×104 6.648×104 6.481×104 6.314×104 6.147×104 5.979×104 0.14 8.020×104 7.827×104 7.633×104 7.438×104 7.245×104 7.050×104 0.16 9.258×104 9.039×104 8.818×104 8.597×104 8.377×104 8.156×104 0.18 1.053×105 1.029×105 1.004×105 9.794×104 9.547×104 9.300×104 0.20 1.185×105 1.158×105 1.131×105 1.103×105 1.076×105 1.049×105 0.22 1.321×105 1.292×105 1.262×105 1.232×105 1.202×105 1.173×105 0.24 1.463×105 1.431×105 1.398×105 1.366×105 1.334×105 1.302×105 0.26 1.609×105 1.575×105 1.541×105 1.506×105 1.471×105 1.437×105 表 14 秋刀鱼渔船横摇阻尼力矩
Table 14. Roll damping moment of series ships generated from a Saury fishing ship
横摇角速度/(rad·s-1) 横摇阻尼/(N·m) 0号船 1号船 2号船 3号船 4号船 5号船 0.08 2.126×104 2.094×104 2.066×104 2.042×104 2.022×104 2.006×104 0.10 3.567×104 3.541×104 3.519×104 3.504×104 3.494×104 3.489×104 0.12 5.616×104 5.603×104 5.598×104 5.600×104 5.610×104 5.627×104 0.14 8.392×104 8.405×104 8.427×104 8.459×104 8.500×104 8.551×104 0.16 1.202×105 1.207×105 1.213×105 1.220×104 1.229×105 1.239×105 0.18 1.662×105 1.672×105 1.683×105 1.697×105 1.711×105 1.728×105 0.20 2.230×105 2.247×105 2.266×105 2.287×105 2.310×105 2.334×105 0.22 2.921×105 2.946×105 2.974×105 3.004×105 3.037×105 3.072×105 0.24 3.744×105 3.780×105 3.819×105 3.861×105 3.906×105 3.954×105 0.26 4.714×105 4.763×105 4.815×105 4.870×105 4.929×105 4.992×105 -
[1] 卜淑霞, 顾民, 鲁江, 等. 参数横摇薄弱性衡准研究[J]. 船舶力学, 2015, 19(1/2): 62–68. BU S X, GU M, LU J, et al. Study on vulnerability criteria of parametric rolling[J]. Journal of Ship Mechanics, 2015, 19(1/2): 62–68 (in Chinese). [2] 马坤, 刘旭, 梅树猛. 船型参数对参数横摇薄弱性衡准的影响[J]. 武汉理工大学学报(交通科学与工程版), 2019, 43(1): 16–19. MA K, LIU X, MEI S M. Influence of ship type parameters on parametric rolling of vulnerability criteria[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2019, 43(1): 16–19 (in Chinese). [3] 王志荣, 周耀华, 户艳宏. 第二代完整稳性衡准研究最新进展[J]. 造船技术, 2019(2): 77–85. doi: 10.3969/j.issn.1000-3878.2019.02.016 WANG Z R, ZHOU Y H, HU Y H. Recent progress in research of 2nd generation intact stability criteria[J]. Marine Technology, 2019(2): 77–85 (in Chinese). doi: 10.3969/j.issn.1000-3878.2019.02.016 [4] 付超, 马山, 段文洋, 等. C11集装箱船参数横摇发生敏感性因素分析和改善措施研究[J]. 中国造船, 2015, 56(S1): 72–80. FU C, MA S, DUAN W Y, et al. Sensitivity factors in parametric rolling of containership C11 and improvement measure[J]. Shipbuilding of China, 2015, 56(S1): 72–80 (in Chinese). [5] 徐光龙. 船舶参数横摇薄弱性衡准研究[D]. 大连: 大连理工大学, 2014. XU G L. A Study on vulnerability criteria for ship parametric roll[D]. Dalian: Dalian University of Technology, 2014 (in Chinese). [6] IMO. Draft explanatory notes on the vulnerability of ships to the parametric rolling stability failure mode[C]//IMO SDC 1/INF. 8. Japan: IMO, 2013. [7] IMO. Draft interim guidelines on vulnerability criteria for the second generation intact stability criteria [C]//IMO SDC 6, 2019. [8] IMO. Proposed amendments to part B of the 2008 IS code to assess the vulnerability of ships to the parametric rolling stability failure mode[C]//IMO SDC 1/INF. 8, 2013. [9] IMO. Proposal of revision of updated draft vulnerability criteria of levels 1 and 2 for the failure modes of pure loss of stability and parametric roll[C]//Submission to the ISCG. Japan: IMO, 2012. [10] 吴建威. 船体型线数值优化软件OPTShip-SJTU船型变换模块开发[D]. 上海: 上海交通大学, 2017. WU J W. Development of hull form deformation module in numerical optimization software OPTShip-SJTU[D]. Shanghai: Shanghai Jiao Tong University, 2017 (in Chinese). [11] COPPEDÉ A, VERNENGO G, VILLA D. A combined approach based on Subdivision surface and free form deformation for smart ship hull form design and variation[J]. Ships and Offshore Structures, 2018, 13(7): 769–778. doi: 10.1080/17445302.2018.1457235 [12] 李冬琴, 李国焕, 戴晶晶, 等. 基于NFFD算法的船体几何变形技术[J]. 船舶工程, 2018, 40(6): 17–22, 35. LI D Q, LI G H, DAI J J, et al. Hull geometric deformation technique based on NFFD algorithm[J]. Ship Engineering, 2018, 40(6): 17–22, 35 (in Chinese). -