留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

船艏UV度对参数横摇薄弱性的影响

马坤 梁舒凡 庄智博

马坤, 梁舒凡, 庄智博. 船艏UV度对参数横摇薄弱性的影响[J]. 中国舰船研究, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910
引用本文: 马坤, 梁舒凡, 庄智博. 船艏UV度对参数横摇薄弱性的影响[J]. 中国舰船研究, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910
MA K, LIANG S F, ZHUANG Z B. Influence of U or V degree of stem on vulnerability of parametric rolling[J]. Chinese Journal of Ship Research, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910
Citation: MA K, LIANG S F, ZHUANG Z B. Influence of U or V degree of stem on vulnerability of parametric rolling[J]. Chinese Journal of Ship Research, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910

船艏UV度对参数横摇薄弱性的影响

doi: 10.19693/j.issn.1673-3185.01910
详细信息
    作者简介:

    马坤,女,1961年生,教授,博士生导师。研究方向:船舶稳性研究。E-mail: makun@dlut.edu.cn

    梁舒凡,男,1995年生,硕士生。研究方向:第二代完整稳性研究

    庄智博,男,1994年生,硕士生。研究方向:第二代完整稳性研究

    通讯作者:

    马坤

  • 中图分类号: U662.2

Influence of U or V degree of stem on vulnerability of parametric rolling

  • 摘要:   目的  以船艏UV度为对象探究其对参数横摇薄弱性的影响。   方法  选取集装箱船、渔政船和渔船这3种不同类型的船舶为母型船,应用自由曲面变换(free-form deformation,FFD)法,改变母型船艏部UV度,生成3个系列共15艘新船型,采用本课题组前期开发的参数横摇薄弱性校核程序对该系列船型进行参数横摇薄弱性计算,并对计算结果进行分析。   结果  结果显示,船艏型线由V型向U型变化时船舶参数横摇薄弱性变好。   结论  所做研究可为参数横摇薄弱性衡准在国内的推广应用提供技术依据,也可为船体型线设计提供参考。
  • 图  1  C11集装箱船型线图

    Figure  1.  Lines of C11 container ship

    图  2  大连400 t渔政船型线图

    Figure  2.  Lines of Dalian 400 t fishery administration ship

    图  3  秋刀鱼渔船型线图

    Figure  3.  Lines of a Saury fishing ship

    图  4  FFD法控制体及坐标系示意图

    Figure  4.  Diagram of the control body and the coordinate system of FFD method

    图  5  C11集装箱船系列0号和5号样船型线图

    Figure  5.  No.0 and No.5 ships' lines generated from C11 container ship

    图  6  大连400 t 渔政船系列0号和5号样船型线图

    Figure  6.  No.0 and No.5 Ships' lines generated from Dalian 400 t fishery administration ship

    图  7  秋刀鱼渔船系列0号和5号样船型线图

    Figure  7.  No.0 and No.5 Ships' lines generated from Saury fishing ship

    图  8  3个系列样船第1层计算结果

    Figure  8.  Calculation results of level 1 of three series ships

    图  9  3个系列样船第2层C2计算结果

    Figure  9.  Calculation results of level 2 C2 of three series ships

    表  1  母型船主要船型参数

    Table  1.   Main parameters of parent ships

    参数C11集装箱船大连400 t渔政船秋刀鱼渔船
    垂线间长/m26248.665.6
    船宽/m40811.1
    型深/m24.454.37
    吃水/m11.52.854.447
    排水体积/m3673845432293
    初稳性高/m1.9021.8620.548
    中横剖面系数0.9580.7520.978
    方形系数0.5590.5160.695
    服务航速/kn251614
    下载: 导出CSV

    表  2  C11集装箱船系列样船主要船型参数

    Table  2.   Main parameters of series ships generated from C11 container ship

    参数0号船1号船2号船3号船4号船5号船
    垂线间长/m262262262262262262
    船宽/m404040404040
    型深/m24.4524.4524.4524.4524.4524.45
    吃水/m11.511.511.511.511.511.5
    排水体积/m367 384 6791868 45268 98669 51770 047
    初稳性高/m1.9021.7931.6851.5811.4821.386
    前半体方形系数0.5050.5120.5190.5270.5340.542
    水线面系数0.7740.7750.7760.7770.7780.779
    下载: 导出CSV

    表  3  大连400 t渔政船系列样船主要船型参数

    Table  3.   Main parameters of series ships generated from Dalian 400 t fishery administration ship

    参数0号船1号船2号船3号船4号船5号船
    垂线间长/m48.648.648.648.648.648.6
    船宽/m888888
    型深/m4.34.34.34.34.34.3
    吃水/m2.852.852.852.852.852.85
    排水体积/m3543545547549551553
    初稳性高/m1.8621.8171.7721.7291.6881.648
    前半体方形系数0.4600.4630.4650.4670.4700.472
    水线面系数0.8300.8250.8200.8150.8100.805
    下载: 导出CSV

    表  4  秋刀鱼渔船系列样船主要船型参数

    Table  4.   Main parameters of series ships generated from a fishing ship

    参数0号船1号船2号船3号船4号船5号船
    垂线间长/m65.665.665.665.665.665.6
    船宽/m11.111.111.111.111.111.1
    型深/m777777
    吃水/m4.4474.4474.4474.4474.4474.447
    排水体积/m3229322942294229522962297
    初稳性高/m0.5480.5300.5120.4940.4770.460
    前半体方形系数0.6750.6750.6750.6750.6750.675
    水线面系数0.9080.9050.9030.9000.8980.895
    下载: 导出CSV

    表  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$衡准值
    01.5591.9020.8190.17
    11.5001.7930.8370.17
    21.4421.6850.8560.17
    31.3861.5810.8760.17
    41.3311.4820.8980.17
    51.2761.3860.9210.17
    下载: 导出CSV

    表  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$衡准值
    00.0421.8620.0220.17
    10.0371.8170.0200.17
    20.0321.7720.0190.17
    30.0271.7300.0160.17
    40.0231.6880.0140.17
    50.0191.6480.0120.17
    下载: 导出CSV

    表  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$衡准值
    00.1200.5480.2180.17
    10.1170.5300.2210.17
    20.1150.5120.2240.17
    30.1120.4940.2280.17
    40.1100.4770.2310.17
    50.1080.4600.2350.17
    下载: 导出CSV

    表  8  3个系列样船第2层C2计算结果

    Table  8.   Calcution results of C2 of three series ships

    船型编号C11
    集装箱船
    大连400 t
    渔政船
    秋刀鱼
    渔船
    衡准值
    00.08571.20×10−62.38×10−30.025
    10.08011.60×10−62.45×10−30.025
    20.07501.60×10−69.86×10−40.025
    30.07283.36×10−67.21×10−40.025
    40.07194.72×10−67.00×10−40.025
    50.06496.48×10−66.81×10−40.025
    下载: 导出CSV

    表  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.010.020.030.040.050.060.070.080.090.10
    000000000000
    100000000038.7
    20000000000
    30000000000
    100000000000
    100000000037.8
    20000000000
    30000000000
    200000000000
    100000000036.7
    20000000000
    30000000000
    300000000000
    10000000038.535.6
    200000000044
    30000000000
    400000000000
    10000000038.535.8
    200000000043.5
    30000000000
    500000000000
    10000000037.734.7
    200000000042.8
    30000000000
    下载: 导出CSV

    表  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 }$横摇幅值/(°)
    02.1791.28820.2
    12.1791.27216.4
    22.1791.25711.1
    32.1791.24138.5
    42.1791.22638.5
    52.1791.21237.7
    下载: 导出CSV

    表  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 }$横摇幅值/(°)
    02.0381.28817.1
    12.0381.27212.7
    22.0381.2574.9
    32.0381.24144.0
    42.0381.22643.5
    52.0381.21242.8
    下载: 导出CSV

    表  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.082.614×1072.623×1072.643×1072.670×1072.705×1072.749×107
    0.104.880×1074.923×1074.982×1075.054×1075.137×1075.235×107
    0.128.220×1078.318×1078.440×1078.581×1078.739×1078.920×107
    0.141.285×1081.303×1081.324×1081.348×1081.375×1081.404×108
    0.161.899×1081.927×1081.961×1081.998×1082.039×1082.085×108
    0.182.684×1082.727×1082.777×1082.831×1082.891×1082.957×108
    0.203.663×1083.724×1083.794×1083.871×1083.953×1084.045×108
    0.224.857×1084.940×1085.035×1085.138×1085.250×1085.373×108
    0.246.288×1086.398×1086.522×1086.658×1086.804×1086.964×108
    0.267.976×1088.118×1088.278×1088.452×1088.638×1088.844×108
    下载: 导出CSV

    表  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.084.481×1044.370×1044.258×1044.145×1044.033×1043.921×104
    0.105.636×1045.497×1045.358×1045.217×1045.078×1044.938×104
    0.126.814×1046.648×1046.481×1046.314×1046.147×1045.979×104
    0.148.020×1047.827×1047.633×1047.438×1047.245×1047.050×104
    0.169.258×1049.039×1048.818×1048.597×1048.377×1048.156×104
    0.181.053×1051.029×1051.004×1059.794×1049.547×1049.300×104
    0.201.185×1051.158×1051.131×1051.103×1051.076×1051.049×105
    0.221.321×1051.292×1051.262×1051.232×1051.202×1051.173×105
    0.241.463×1051.431×1051.398×1051.366×1051.334×1051.302×105
    0.261.609×1051.575×1051.541×1051.506×1051.471×1051.437×105
    下载: 导出CSV

    表  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.082.126×1042.094×1042.066×1042.042×1042.022×1042.006×104
    0.103.567×1043.541×1043.519×1043.504×1043.494×1043.489×104
    0.125.616×1045.603×1045.598×1045.600×1045.610×1045.627×104
    0.148.392×1048.405×1048.427×1048.459×1048.500×1048.551×104
    0.161.202×1051.207×1051.213×1051.220×1041.229×1051.239×105
    0.181.662×1051.672×1051.683×1051.697×1051.711×1051.728×105
    0.202.230×1052.247×1052.266×1052.287×1052.310×1052.334×105
    0.222.921×1052.946×1052.974×1053.004×1053.037×1053.072×105
    0.243.744×1053.780×1053.819×1053.861×1053.906×1053.954×105
    0.264.714×1054.763×1054.815×1054.870×1054.929×1054.992×105
    下载: 导出CSV
  • [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).
  • [1] 何洋洋, 王馨怡, 董晶.  基于经验小波变换与谱峭度的船舶轴系故障特征提取方法 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01771
    [2] 李宗威, 纪刚, 周其斗.  多自由度系统基于单通带近似法的振动局域化分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01811
    [3] 马坤, 梅树猛, 于津钊.  瘫船稳性横摇周期计算方法讨论分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01574
    [4] 朱军, 朱韬, 王智宇, 夏齐强, 黄昆仑, 葛义军.  船舶参数横摇数值计算与力学机理分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01618
    [5] 马坤, 胡高源, 梁舒凡.  极地船舶积冰计算及其对参数横摇的影响分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01572
    [6] 张彦儒, 林焰, 陆丛红, 纪卓尚.  轻量化NURBS船体曲面自行设计垂向参数化方法 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.05.004
    [7] 马建文, 张安西, 周兆欣, 郭绍义.  小水线面双体船五自由度运动建模与仿真 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.02.017
    [8] 马骋远, 马宁, 王廷昊, 顾解忡.  超大型集装箱船参数横摇全面校核与安全评估 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.03.002
    [9] 董润鹏, 董文才, 姚朝帮.  深V型船零速横摇改善措施 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2015.06.003
    [10] 张晓, 杨和振.  不规则波中船舶参数横摇的概率分析 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2015.03.006
    [11] 姜治芳, 黄卫刚, 邱辽原, 李亚军, 李鹏.  内倾船型的非线性横摇仿真分析 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2014.01.006
    [12] 徐静, 顾解忡, 马宁.  规则波六自由度回转运动预报 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2014.03.003
    [13] 侯岳, 鄢凯, 浦金云.  舰艇横摇进水过程姿态参数数值求解与实验研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2014.04.007
    [14] 袁东红.  多自由度刚体隔振系统冲击响应计算方法 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2013.04.011
    [15] 蒋颉, 丁勇.  减摇水舱挡板开度对稳定矩影响的研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2011.03.008
    [16] 韩小溪, 胡光振, 杨元.  基于BP神经网络的维修性设计参数灵敏度分析 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2008.06.016
    [17] 唐友刚, 李红霞.  随机斜浪中船舶参数——强迫激励横摇运动计算 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2008.01.002
    [18] 王中, 卢晓平, 付攀.  基于细分曲面的船体曲面造型及静水力参数计算 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2007.01.001
    [19] 黎庆芬, 曾广武, 肖伟.  潜器圆柱形耐压壳体强度和稳性计算及衡准方法研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2006.02.011
    [20] 张寿桂.  船舶在大风浪中横摇稳性变化的理论研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2006.01.010
  • 加载中
图(9) / 表 (14)
计量
  • 文章访问数:  21
  • HTML全文浏览量:  13
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-26
  • 修回日期:  2020-06-20
  • 网络出版日期:  2020-12-10

船艏UV度对参数横摇薄弱性的影响

doi: 10.19693/j.issn.1673-3185.01910
    作者简介:

    马坤,女,1961年生,教授,博士生导师。研究方向:船舶稳性研究。E-mail: makun@dlut.edu.cn

    梁舒凡,男,1995年生,硕士生。研究方向:第二代完整稳性研究

    庄智博,男,1994年生,硕士生。研究方向:第二代完整稳性研究

    通讯作者: 马坤
  • 中图分类号: U662.2

摘要:    目的  以船艏UV度为对象探究其对参数横摇薄弱性的影响。   方法  选取集装箱船、渔政船和渔船这3种不同类型的船舶为母型船,应用自由曲面变换(free-form deformation,FFD)法,改变母型船艏部UV度,生成3个系列共15艘新船型,采用本课题组前期开发的参数横摇薄弱性校核程序对该系列船型进行参数横摇薄弱性计算,并对计算结果进行分析。   结果  结果显示,船艏型线由V型向U型变化时船舶参数横摇薄弱性变好。   结论  所做研究可为参数横摇薄弱性衡准在国内的推广应用提供技术依据,也可为船体型线设计提供参考。

English Abstract

马坤, 梁舒凡, 庄智博. 船艏UV度对参数横摇薄弱性的影响[J]. 中国舰船研究, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910
引用本文: 马坤, 梁舒凡, 庄智博. 船艏UV度对参数横摇薄弱性的影响[J]. 中国舰船研究, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910
MA K, LIANG S F, ZHUANG Z B. Influence of U or V degree of stem on vulnerability of parametric rolling[J]. Chinese Journal of Ship Research, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910
Citation: MA K, LIANG S F, ZHUANG Z B. Influence of U or V degree of stem on vulnerability of parametric rolling[J]. Chinese Journal of Ship Research, 2020, 15(s1): 1–8 doi: 10.19693/j.issn.1673-3185.01910
    • 船舶稳性是船舶航行安全的重要保障,一直以来受到国际海事组织(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渔政船秋刀鱼渔船
      垂线间长/m26248.665.6
      船宽/m40811.1
      型深/m24.454.37
      吃水/m11.52.854.447
      排水体积/m3673845432293
      初稳性高/m1.9021.8620.548
      中横剖面系数0.9580.7520.978
      方形系数0.5590.5160.695
      服务航速/kn251614

      图  1  C11集装箱船型线图

      Figure 1.  Lines of C11 container ship

      图  2  大连400 t渔政船型线图

      Figure 2.  Lines of Dalian 400 t fishery administration ship

      图  3  秋刀鱼渔船型线图

      Figure 3.  Lines of a Saury fishing ship

    • 本文采用FFD法实现母型船的型线变换。FFD法是船型设计优化过程中常用的方法,该方法可应用于任何曲面,不受曲面表示方法的限制,还可对物体的局部或整体进行变形,变形后的曲面能够保持几何连续性。

      FFD法的基本步骤如下:

      1) 选定一个长方体形状的控制体包围待变形曲面;

      2) 在控制体内选择一系列控制点;

      3) 调节控制点位置改变控制体形状,将控制体的变形传递给待变形曲面,实现曲面变形[10-11]

      图4所示为FFD法控制体及坐标系示意图。其中,O-XYZ为全局坐标系,O′-RST为控制体坐标系,控制体内任意一点的全局坐标与控制点的全局坐标满足伯恩斯坦(Bernstein)多项式:

      图  4  FFD法控制体及坐标系示意图

      Figure 4.  Diagram of the control body and the coordinate system of FFD method

      $$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所示。

      图  5  C11集装箱船系列0号和5号样船型线图

      Figure 5.  No.0 and No.5 ships' lines generated from C11 container ship

      图  6  大连400 t 渔政船系列0号和5号样船型线图

      Figure 6.  No.0 and No.5 Ships' lines generated from Dalian 400 t fishery administration ship

      图  7  秋刀鱼渔船系列0号和5号样船型线图

      Figure 7.  No.0 and No.5 Ships' lines generated from Saury fishing ship

      表 2  C11集装箱船系列样船主要船型参数

      Table 2.  Main parameters of series ships generated from C11 container ship

      参数0号船1号船2号船3号船4号船5号船
      垂线间长/m262262262262262262
      船宽/m404040404040
      型深/m24.4524.4524.4524.4524.4524.45
      吃水/m11.511.511.511.511.511.5
      排水体积/m367 384 6791868 45268 98669 51770 047
      初稳性高/m1.9021.7931.6851.5811.4821.386
      前半体方形系数0.5050.5120.5190.5270.5340.542
      水线面系数0.7740.7750.7760.7770.7780.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号船
      垂线间长/m48.648.648.648.648.648.6
      船宽/m888888
      型深/m4.34.34.34.34.34.3
      吃水/m2.852.852.852.852.852.85
      排水体积/m3543545547549551553
      初稳性高/m1.8621.8171.7721.7291.6881.648
      前半体方形系数0.4600.4630.4650.4670.4700.472
      水线面系数0.8300.8250.8200.8150.8100.805

      表 4  秋刀鱼渔船系列样船主要船型参数

      Table 4.  Main parameters of series ships generated from a fishing ship

      参数0号船1号船2号船3号船4号船5号船
      垂线间长/m65.665.665.665.665.665.6
      船宽/m11.111.111.111.111.111.1
      型深/m777777
      吃水/m4.4474.4474.4474.4474.4474.447
      排水体积/m3229322942294229522962297
      初稳性高/m0.5480.5300.5120.4940.4770.460
      前半体方形系数0.6750.6750.6750.6750.6750.675
      水线面系数0.9080.9050.9030.9000.8980.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$衡准值
      01.5591.9020.8190.17
      11.5001.7930.8370.17
      21.4421.6850.8560.17
      31.3861.5810.8760.17
      41.3311.4820.8980.17
      51.2761.3860.9210.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$衡准值
      00.0421.8620.0220.17
      10.0371.8170.0200.17
      20.0321.7720.0190.17
      30.0271.7300.0160.17
      40.0231.6880.0140.17
      50.0191.6480.0120.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$衡准值
      00.1200.5480.2180.17
      10.1170.5300.2210.17
      20.1150.5120.2240.17
      30.1120.4940.2280.17
      40.1100.4770.2310.17
      50.1080.4600.2350.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度对参数横摇薄弱性的影响大小。

      图  8  3个系列样船第1层计算结果

      Figure 8.  Calculation results of level 1 of three series ships

    • 研究船体艏部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
      渔政船
      秋刀鱼
      渔船
      衡准值
      00.08571.20×10−62.38×10−30.025
      10.08011.60×10−62.45×10−30.025
      20.07501.60×10−69.86×10−40.025
      30.07283.36×10−67.21×10−40.025
      40.07194.72×10−67.00×10−40.025
      50.06496.48×10−66.81×10−40.025

      校核结果表明,C11集装箱船系列样船未能通过参数横摇薄弱性第2层衡准校核,大连400 t渔政船和秋刀鱼渔船通过了校核。随着船舶艏部UV度由V型向U型变化,C11集装箱船与秋刀鱼渔船的参数横摇薄弱性C2整体呈下降趋势,大连400 t渔政船的参数横摇薄弱性C2小幅上升,如图9所示。

      图  9  3个系列样船第2层C2计算结果

      Figure 9.  Calculation results of level 2 C2 of three series ships

      表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.010.020.030.040.050.060.070.080.090.10
      000000000000
      100000000038.7
      20000000000
      30000000000
      100000000000
      100000000037.8
      20000000000
      30000000000
      200000000000
      100000000036.7
      20000000000
      30000000000
      300000000000
      10000000038.535.6
      200000000044
      30000000000
      400000000000
      10000000038.535.8
      200000000043.5
      30000000000
      500000000000
      10000000037.734.7
      200000000042.8
      30000000000
      $$\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 }$横摇幅值/(°)
      02.1791.28820.2
      12.1791.27216.4
      22.1791.25711.1
      32.1791.24138.5
      42.1791.22638.5
      52.1791.21237.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 }$横摇幅值/(°)
      02.0381.28817.1
      12.0381.27212.7
      22.0381.2574.9
      32.0381.24144.0
      42.0381.22643.5
      52.0381.21242.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.082.614×1072.623×1072.643×1072.670×1072.705×1072.749×107
      0.104.880×1074.923×1074.982×1075.054×1075.137×1075.235×107
      0.128.220×1078.318×1078.440×1078.581×1078.739×1078.920×107
      0.141.285×1081.303×1081.324×1081.348×1081.375×1081.404×108
      0.161.899×1081.927×1081.961×1081.998×1082.039×1082.085×108
      0.182.684×1082.727×1082.777×1082.831×1082.891×1082.957×108
      0.203.663×1083.724×1083.794×1083.871×1083.953×1084.045×108
      0.224.857×1084.940×1085.035×1085.138×1085.250×1085.373×108
      0.246.288×1086.398×1086.522×1086.658×1086.804×1086.964×108
      0.267.976×1088.118×1088.278×1088.452×1088.638×1088.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.084.481×1044.370×1044.258×1044.145×1044.033×1043.921×104
      0.105.636×1045.497×1045.358×1045.217×1045.078×1044.938×104
      0.126.814×1046.648×1046.481×1046.314×1046.147×1045.979×104
      0.148.020×1047.827×1047.633×1047.438×1047.245×1047.050×104
      0.169.258×1049.039×1048.818×1048.597×1048.377×1048.156×104
      0.181.053×1051.029×1051.004×1059.794×1049.547×1049.300×104
      0.201.185×1051.158×1051.131×1051.103×1051.076×1051.049×105
      0.221.321×1051.292×1051.262×1051.232×1051.202×1051.173×105
      0.241.463×1051.431×1051.398×1051.366×1051.334×1051.302×105
      0.261.609×1051.575×1051.541×1051.506×1051.471×1051.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.082.126×1042.094×1042.066×1042.042×1042.022×1042.006×104
      0.103.567×1043.541×1043.519×1043.504×1043.494×1043.489×104
      0.125.616×1045.603×1045.598×1045.600×1045.610×1045.627×104
      0.148.392×1048.405×1048.427×1048.459×1048.500×1048.551×104
      0.161.202×1051.207×1051.213×1051.220×1041.229×1051.239×105
      0.181.662×1051.672×1051.683×1051.697×1051.711×1051.728×105
      0.202.230×1052.247×1052.266×1052.287×1052.310×1052.334×105
      0.222.921×1052.946×1052.974×1053.004×1053.037×1053.072×105
      0.243.744×1053.780×1053.819×1053.861×1053.906×1053.954×105
      0.264.714×1054.763×1054.815×1054.870×1054.929×1054.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的计算结果上升。

参考文献 (12)

目录

    /

    返回文章
    返回