留言板

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

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

厚板结构对波浪作用下带月池船体垂荡运动的影响分析

曾智宏 邓潇潇 孙雷

曾智宏, 邓潇潇, 孙雷. 厚板结构对波浪作用下带月池船体垂荡运动的影响分析[J]. 中国舰船研究, 2020, 15(增刊): 21–34 doi:  10.19693/j.issn.1673-3185.01993
引用本文: 曾智宏, 邓潇潇, 孙雷. 厚板结构对波浪作用下带月池船体垂荡运动的影响分析[J]. 中国舰船研究, 2020, 15(增刊): 21–34 doi:  10.19693/j.issn.1673-3185.01993
ZENG Z H, DENG X X, SUN L. effect analysis of thick plate structure on heave motion of ship with moonpool in waves [J]. Chinese Journal of Ship Research, 2020, 15(Supp): 21–34 doi:  10.19693/j.issn.1673-3185.01993
Citation: ZENG Z H, DENG X X, SUN L. effect analysis of thick plate structure on heave motion of ship with moonpool in waves [J]. Chinese Journal of Ship Research, 2020, 15(Supp): 21–34 doi:  10.19693/j.issn.1673-3185.01993

厚板结构对波浪作用下带月池船体垂荡运动的影响分析

doi: 10.19693/j.issn.1673-3185.01993
基金项目: 国家自然科学基金资助项目(51679035);辽宁省兴辽英才计划资助项目(XLYC1908027);中央高校基本科研业务费资助项目(DUT20TD108);大连理工大学大学生创新创业训练计划资助项目(2019101410401011423)
详细信息
    作者简介:

    曾智宏,男,1998年生,本科生。研究方向:船舶水动力学。E-mail:qq1774481948@mail.dlut.edu.cn

    邓潇潇,男,1994年生,硕士。研究方向:带月池平台水动力分析。E-mail:3051478711@qq.com

    孙雷,男,1982年生,博士,副教授。研究方向:船舶与海洋工程水动力载荷。E-mail:sunlei@dlut.edu.cn

    通讯作者:

    孙雷

  • 中图分类号: U661.1

effect analysis of thick plate structure on heave motion of ship with moonpool in waves

  • 摘要:   目的  含月池的船舶在波浪作用下月池会对船舶运动产生不利影响,引起船体大幅运动,尤其是垂荡方向,导致相关功能设备受到损坏。为此,尝试通过在月池中布置厚板来解决此问题。  方法  基于势流理论和含月池钻井船频域数值模型,计算和分析不同长宽比和吃水的月池以及对称设置的不同高度厚板对船体垂荡运动的影响规律。选择某浮式生产储卸油装置(FPSO),利用AQWA软件对波浪作用下FPSO船体垂荡运动展开频域数值计算,并与文献所使用的方法及实验结果进行对比,验证所提数值模型及方法的可靠性和有效性。  结果  若设置的厚板上表面正好与静水面齐平,此时可有效抑制船体因月池内流体共振而产生的剧烈垂荡运动。  结论  研究结果可为含月池船体抑制垂荡运动提供参考。
  • 图  1  船舶数学模型示意图

    Figure  1.  Diagram of mathematical ship model

    图  2  不同长宽比月池示意图

    Figure  2.  Diagram of the moonpool with different length-width ratios

    图  3  厚板示意图和月池横剖面图

    Figure  3.  Diagrams of thick plate and cross section of moonpool

    图  4  横浪及迎浪下钻井船运动响应(工况2)

    Figure  4.  RAO of the drillship in beam sea and head sea (case 2)

    图  5  横浪及迎浪下钻井船运动响应(工况5)

    Figure  5.  RAO of the drillship in beam sea and head sea (case 5)

    图  6  标模1 FPSO船体平面视图[9](单位:m)

    Figure  6.  Plane views of the No.1 FPSO ship model [9] (Unit in meter)

    图  7  标模1 FPSO船体网格剖分

    Figure  7.  Mesh generation of the No.1 FPSO ship hull

    图  8  90°浪向下FPSO船运动响应

    Figure  8.  RAO of the FPSO ship when α0=90°

    图  9  120°浪向下FPSO船运动响应

    Figure  9.  RAO of the FPSO ship when α0=120°

    图  10  180°浪向下FPSO船运动响应

    Figure  10.  RAO of the FPSO ship when α0=180°

    图  11  标模2平面视图(单位:m)

    Figure  11.  The plane view of No.2 FPSO ship model (Unit in meter)

    图  12  标模2 FPSO船体网格剖分

    Figure  12.  Mesh generation of the No.2 FPSO ship hull model

    图  13  迎浪下FPSO船垂荡运动响应(MP1)

    Figure  13.  Heaving RAO of the FPSO ship in head sea ( case MP1)

    图  14  横浪下FPSO船垂荡运动响应 (MP1工况)

    Figure  14.  Heaving RAO of the FPSO ship in beam sea (case MP1)

    图  15  迎浪下FPSO船垂荡运动响应(MP2工况)

    Figure  15.  Heaving RAO of the FPSO ship in head sea (case MP2)

    图  16  横浪下FPSO船垂荡运动响应(MP2工况)

    Figure  16.  Heaving RAO of the FPSO ship in beam sea(case MP2)

    图  17  迎浪下FPSO船垂荡运动响应(MP3工况)

    Figure  17.  Heaving RAO of the FPSO ship in head sea(case MP3)

    图  18  横浪下FPSO船垂荡运动响应(MP3工况)

    Figure  18.  Heaving RAO of the FPSO ship in beam sea(case MP3)

    图  19  月池固有频率对比

    Figure  19.  Comparison on the nature frequency of moonpool

    图  20  横浪及不同厚板高度下钻井船垂荡运动响应(工况1)

    Figure  20.  Heaving RAO of the drillship in beam sea with different heights of thick plate(case 1)

    图  21  横浪及不同厚板高度下钻井船垂荡运动响应(工况2)

    Figure  21.  Heaving RAO of the drillship in beam sea with different heights of thick plate(case 2)

    图  22  横浪及不同厚板高度下钻井船垂荡运动响应(工况3)

    Figure  22.  Heaving RAO of the drillship in beam sea with different heights of thick plate(case 3)

    图  23  横浪及不同厚板高度下钻井船垂荡运动响应(工况4)

    Figure  23.  Heaving RAO of the drillship in beam sea with different heights of thick plate (case 4)

    图  24  横浪及不同厚板高度下钻井船垂荡运动响应(工况5)

    Figure  24.  Heaving RAO of the drillship in beam sea with different heights of thick plate (case 5)

    图  25  0°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

    Figure  25.  Heaving RAO of the drillship with different heights of thick plate when α0= 0° (case 2)

    图  26  45°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

    Figure  26.  Heaving RAO of the drillship with different heights of thick plate when α0= 45° (case 2)

    图  27  135°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

    Figure  27.  Heaving RAO of the drillship with different heights of thick plate when α0=135° (case 2)

    图  28  180°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

    Figure  28.  Heaving RAO of the drillship with different heights of thick plate when α0=180° (case 2)

    图  29  横浪下厚板上表面高于静水面时钻井船垂荡运动响应(工况2)

    Figure  29.  Heaving RAO of the drillship with thick plate above the waterline in beam sea(case 2)

    表  1  含月池钻井船的主要参数

    Table  1.   Main parameters of the drillship with moonpool

    参数工况1工况2工况3工况4工况5
    总长/m 112 112 112 112 112
    型宽/m 28.1 28.1 28.1 28.1 28.1
    型深/m 98 98 98 98 98
    吃水/m 6 7.5 9 7.5 7.5
    月池长/m 22.6 22.6 22.6 16 27.7
    月池宽/m 11.3 11.3 11.3 16 9.2
    排水量/t 14 474 18 379 22 248.4 18 379 18 379
    重心垂向位置/m
    (相对船基线)
    9.2 8.2 7.2 8.2 8.2
    厚板厚度/m 1 1 1 1 1
    厚板宽度/m 1.65 1.65 1.65 2.00 1.00
    横摇惯性半径/m 7.9 7.9 8.037 7.9 7.9
    纵摇惯性半径/m 27.864 27.954 28.149 27.954 27.954
    艏摇惯性半径/m 28.899 28.905 29.272 28.905 28.905
    下载: 导出CSV

    表  2  压载状态下标模1 FPSO船主要参数

    Table  2.   Main parameters of the No.1 FPSO ship in ballast condition

    参数数值
    垂线间长/m285
    船宽/m63
    吃水/m13
    排水量/t225 518
    重心垂向位置/m16.5
    重心纵向位置/m142.5
    横摇惯性半径/m19.45
    纵摇惯性半径/m71.25
    艏摇惯性半径/m71.25
    下载: 导出CSV

    表  3  标模2 FPSO船主尺度参数

    Table  3.   Main parameters of the No.2 FPSO ship model

    参数船模尺寸实船尺寸
    船长/m4.0138
    船宽/m0.827.6
    吃水/m0.26.9
    MP1尺度 (VMP/VS = 1.5%)
    月池长度/m0.26.9
    月池宽度/m0.26.9
    月池中心/m(0.65, 0, 0)(22.4, 0, 0)
    船舶质量/t0.56721 627.8
    船舶重心/m(−0.11, 0, 0.055)(−3.8, 0, 1.9)
    惯性半径/m(0.28, 1.02, 0.76)(9.7, 35.2, 26.2)
    MP2尺度 (VMP/VS = 6.4%)
    月池长度/m0.413.8
    月池宽度/m0.413.8
    月池中心/m(0.65, 0, 0)(22.4, 0, 0)
    船舶质量/t0.502 6920 642.3
    船舶重心/m(−0.15, 0, 0.057)(−5.2, 0, 2.0)
    惯性半径/m(0.29, 1.024, 0.75)(10.0, 35.2, 25.9)
    MP3尺度 (VMP/VS = 42.7%)
    月池长度/m2.069
    月池宽度/m0.413.8
    月池中心/m(−0.15, 0, 0)(−5.2, 0, 0)
    船舶质量/t0.374 6915 386.1
    船舶重心/m(−0.08, 0, 0.083)(−2.76, 0, 2.9)
    惯性半径/m(0.30, 1.11, 1.13)(10.4, 38.3, 39.0)
    下载: 导出CSV

    表  4  月池固有频率

    Table  4.   Natural frequency of moonpool

    工况数值模拟的固有频率/Hz经验公式计算的固有频率/Hz
    10.1570.160
    20.1470.149
    30.1380.140
    40.1450.148
    50.1480.150
    下载: 导出CSV

    表  5  横浪下钻井船垂荡运动响应(工况1)

    Table  5.   Heaving RAO of the drillship in beam sea (case 1)

    h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m−1)
    无厚板0.1571.201
    1.00.1471.406
    2.00.1481.363
    3.00.1481.335
    4.00.1491.309
    5.00.1481.302
    5.20.1471.334
    5.40.1461.339
    5.60.1461.347
    5.70.1451.375
    5.80.1401.502
    5.90.1281.626
    6.00.1800.865
    下载: 导出CSV

    表  6  横浪下钻井船船体垂荡运动响应(工况2)

    Table  6.   Heaving RAO of drillship in beam sea (case 2)

    h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m−1)
    无厚板 0.147 1.396
    1.0 0.138 1.727
    2.0 0.139 1.670
    3.0 0.139 1.627
    4.0 0.140 1.582
    5.0 0.140 1.569
    6.0 0.140 1.539
    6.5 0.140 1.539
    6.7 0.139 1.557
    6.9 0.138 1.592
    7.1 0.138 1.591
    7.2 0.137 1.625
    7.3 0.133 1.790
    7.4 0.124 1.895
    7.5 0.202 0.405
    下载: 导出CSV

    表  7  横浪下钻进船垂荡运动响应(工况3)

    Table  7.   Heaving RAO of the drillship in beam sea (case 3)

    h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m−1)
    无厚板0.1381.701
    2.00.1322.011
    4.00.1321.898
    6.00.1331.860
    8.00.1331.798
    8.20.1321.833
    8.40.1311.851
    8.60.1321.806
    8.70.1311.852
    8.80.1282.105
    8.90.1202.268
    9.00.1600.981
    下载: 导出CSV

    表  8  横浪下钻井船垂荡运动响应(工况4)

    Table  8.   Heaving RAO of the drillship in beam sea (case 4)

    h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m-1)
    无厚板0.1451.301
    1.50.1381.614
    3.00.1391.912
    4.50.1391.476
    6.00.1391.429
    7.50.1800.613
    下载: 导出CSV

    表  9  横浪下钻井船船体垂荡运动响应(工况5)

    Table  9.   Heaving RAO of the drillship in beam sea (case 5)

    h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m-1)
    无厚板0.1481.452
    1.50.1441.583
    3.00.1441.547
    4.50.1451.505
    6.00.1451.469
    7.50.2010.543
    下载: 导出CSV

    表  10  横浪下厚板上表面高于静水面时钻井船垂荡运动响应(工况2)

    Table  10.   Heaving RAO of the drillship with thick plate above the waterline in beam sea (case 2)

    h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m-1)
    无厚板0.1471.396
    7.50.2020.405
    7.60.1690.970
    7.80.1680.934
    8.00.1690.969
    下载: 导出CSV
  • 孙采微, 杨建民, 吕海宁. 波浪作用下带月池结构船体运动数值预报[J]. 海洋工程, 2013, 31(4): 21–29.

    SUN C W, YANG J M, LYU H N. Numerical investigation on motions of vessel with moonpool in wave conditions[J]. The Ocean Engineering, 2013, 31(4): 21–29 (in Chinese).
    黄磊. 开口式月池水动力特性研究[D]. 天津: 天津大学, 2014.

    HUANG L. Study on hydrodynamic characteristics of moonpool with bottom open[D]. Tianjin: Tianjin University, 2014 (in Chinese).
    姚震球, 郭琦, 凌宏杰, 等. 带月池超深水钻井船的水动力特性[J]. 船舶工程, 2017, 39(12): 22–26.

    YAO Z Q, GUO Q, LING H J, et al. Hydrodynamic performance of ultra-deep water drilling ship with moonpool[J]. Ship Engineering, 2017, 39(12): 22–26 (in Chinese).
    张永恒. 深海Truss Spar平台与月池内流体的耦合运动及参数敏感性研究[D]. 天津: 天津大学, 2014.

    ZHANG Y H. Research on coupled motion and parametric sensitivity between truss spar platform and moonpool water in deep sea[D]. Tianjin: Tianjin University, 2014 (in Chinese).
    张利军, 杨朕, 赵志坚, 等. 月池共振对超深水钻井船运动性能的影响[J]. 船舶工程, 2018, 40(5): 14–18, 34.

    ZHANG L J, YANG Z, ZHAO Z J, et al. Effect of moon pool resonance on motion performance of ultra-deepwater drillship[J]. Ship Engineering, 2018, 40(5): 14–18, 34 (in Chinese).
    FUKUDA K. Behavior of water in vertical well with bottom opening of ship, and its effects on ship-motion[J]. Journal of the Society of Naval Architects of Japan, 1977, 1977(141): 107–122. doi:  10.2534/jjasnaoe1968.1977.107
    FALTINSEN O M, ROGNEBAKKE O F, TIMOKHA A N. Two-dimensional resonant piston-like sloshing in a moonpool[J]. Journal of Fluid Mechanics, 2007, 575: 359–397. doi:  10.1017/S002211200600440X
    MOLIN B. On the piston and sloshing modes in moonpools[J]. Journal of Fluid Mechanics, 2001, 430: 27–50. doi:  10.1017/S0022112000002871
    NAM B W, KIM Y, KIM D W, et al. Experimental and numerical studies on ship motion responses coupled with sloshing in waves[J]. Journal of Ship Research, 2009, 53(2): 68–82.
    RAVINTHRAKUMAR S, KRISTIANSEN T, MOLIN B, et al. Coupled vessel and moonpool responses in regular and irregular waves[J]. Applied Ocean Research, 2020, 96: 102010. doi:  10.1016/j.apor.2019.102010
  • [1] 王路才, 周其斗2020: 桨叶振动对螺旋桨垂向激励下潜艇结构振动与声辐射的影响 , 中国舰船研究, 15, 161-166. doi: 10.19693/j.issn.1673-3185.01610
    [2] 游晶越, 赵耀, 张赣波, 储炜2020: 艇体弹性耦合边界条件下轴系纵振反共振隔振分析 , 中国舰船研究, 15, 137-142. doi: 10.19693/j.issn.1673-3185.01802
    [3] 詹星宇, 毛筱菲2020: 船舶月池的阻力与流场特性及其改进型式 , 中国舰船研究, 15, 45-53,101. doi: 10.19693/j.issn.1673-3185.01594
    [4] 孙雷, 邓潇潇, 曾智宏, 刘昌凤2020: 波浪作用下钻井船的共振数值计算及分析 , 中国舰船研究, 15, 90-105, 114. doi: 10.19693/j.issn.1673-3185.01807
    [5] 胡泽超, 何琳, 徐伟, 李正民, 赵兴乾2019: 船舶推进轴系纵向振动共振转换器的优化设计 , 中国舰船研究, 14, 107-113. doi: 10.19693/j.issn.1673-3185.01077
    [6] 李良伟, 赵耀2019: 基于波动理论的共振转换器减振特性分析 , 中国舰船研究, 14, 131-137. doi: 10.19693/j.issn.1673-3185.01445
    [7] 王硕, 宗智2019: 基于蒙特卡洛方法的三体船侧体运动响应优化 , 中国舰船研究, 14, 122-128. doi: 10.19693/j.issn.1673-3185.01348
    [8] 刘梦超, 刘延俊, 薛钢, 吴瀚崚2018: 参数化单元边界元法解势流速度场问题 , 中国舰船研究, 13, 77-84,90. doi: 10.19693/j.issn.1673-3185.01158
    [9] 姚绪梁, 孟令卫, 牛小丽2017: 基于航行阻力优化的近水面机器人减纵摇控制 , 中国舰船研究, 12, 111-119. doi: 10.3969/j.issn.1673-3185.2017.03.016
    [10] 陈长海, 徐文献, 朱锡, 侯海量2015: 超高分子量聚乙烯纤维增强层合厚板抗高速钝头弹侵彻的理论模型 , 中国舰船研究, 10, 63-69,83. doi: 10.3969/j.issn.1673-3185.2015.03.011
    [11] 李井煜, 卢晓平, 赵鹏伟2015: 直接边界元法解势流速度场问题 , 中国舰船研究, 10, 68-75. doi: 10.3969/j.issn.1673-3185.2015.01.010
    [12] 韩玉超, 卢晓平, 王中2015: 无限区域二维势流直接边界元法精度分析 , 中国舰船研究, 10, 39-45,54. doi: 10.3969/j.issn.1673-3185.2015.04.006
    [13] 张弩2014: 水下爆炸气泡作用下船体总纵强度估算方法 , 中国舰船研究, 9, 14-18,25. doi: 10.3969/j.issn.1673-3185.2014.06.003
    [14] 黄毅, 张弩2013: 水下爆炸气泡脉动载荷影响因素分析 , 中国舰船研究, 8, 33-39. doi: 10.3969/j.issn.1673-3185.2013.06.006
    [15] 黄昆仑, 朱军, 葛义军2012: 频率散射对船舶垂荡与纵摇耦合运动的影响 , 中国舰船研究, 7, 19-24. doi: 10.3969/j.issn.1673-3185.2012.03.004
    [16] 詹金林, 卢晓平, 崔焰, 王家林2011: 基于势流理论的船舶操纵水动力分析 , 中国舰船研究, 6, 41-45. doi: 10.3969/j.issn.16733185.2011.05.009
    [17] 晋文菊, 王志东, 李力军, 杨爽2011: 滑行艇升沉纵摇运动的二维数值预报 , 中国舰船研究, 6, 66-70,74. doi: 10.3969/j.issn.1673-3185.2011.06.013
    [18] 莫立新, 涂三, 周心桃2009: 一种结构垂向冲击载荷测量方法分析研究 , 中国舰船研究, 4, 38-41. doi: 10.3969/j.issn.1673-3185.2009.04.008
    [19] 黄晓琼, 陈立, 杨雄辉, 周心桃2009: 三体船耐波性预报 , 中国舰船研究, 4, 42-47. doi: 10.3969/j.issn.1673-3185.2009.06.010
    [20] 张强, 洪明2007: 损伤引起加筋板结构屈曲性态变异研究 , 中国舰船研究, 2, 24-29. doi: 10.3969/j.issn.1673-3185.2007.01.005
  • 加载中
图(29) / 表 (10)
计量
  • 文章访问数:  13
  • HTML全文浏览量:  7
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-09
  • 修回日期:  2020-09-07
  • 网络出版日期:  2020-12-10

厚板结构对波浪作用下带月池船体垂荡运动的影响分析

doi: 10.19693/j.issn.1673-3185.01993
    基金项目:  国家自然科学基金资助项目(51679035);辽宁省兴辽英才计划资助项目(XLYC1908027);中央高校基本科研业务费资助项目(DUT20TD108);大连理工大学大学生创新创业训练计划资助项目(2019101410401011423)
    作者简介:

    曾智宏,男,1998年生,本科生。研究方向:船舶水动力学。E-mail:qq1774481948@mail.dlut.edu.cn

    邓潇潇,男,1994年生,硕士。研究方向:带月池平台水动力分析。E-mail:3051478711@qq.com

    孙雷,男,1982年生,博士,副教授。研究方向:船舶与海洋工程水动力载荷。E-mail:sunlei@dlut.edu.cn

    通讯作者: 孙雷
  • 中图分类号: U661.1

摘要:   目的  含月池的船舶在波浪作用下月池会对船舶运动产生不利影响,引起船体大幅运动,尤其是垂荡方向,导致相关功能设备受到损坏。为此,尝试通过在月池中布置厚板来解决此问题。  方法  基于势流理论和含月池钻井船频域数值模型,计算和分析不同长宽比和吃水的月池以及对称设置的不同高度厚板对船体垂荡运动的影响规律。选择某浮式生产储卸油装置(FPSO),利用AQWA软件对波浪作用下FPSO船体垂荡运动展开频域数值计算,并与文献所使用的方法及实验结果进行对比,验证所提数值模型及方法的可靠性和有效性。  结果  若设置的厚板上表面正好与静水面齐平,此时可有效抑制船体因月池内流体共振而产生的剧烈垂荡运动。  结论  研究结果可为含月池船体抑制垂荡运动提供参考。

English Abstract

曾智宏, 邓潇潇, 孙雷. 厚板结构对波浪作用下带月池船体垂荡运动的影响分析[J]. 中国舰船研究, 2020, 15(增刊): 21–34 doi:  10.19693/j.issn.1673-3185.01993
引用本文: 曾智宏, 邓潇潇, 孙雷. 厚板结构对波浪作用下带月池船体垂荡运动的影响分析[J]. 中国舰船研究, 2020, 15(增刊): 21–34 doi:  10.19693/j.issn.1673-3185.01993
ZENG Z H, DENG X X, SUN L. effect analysis of thick plate structure on heave motion of ship with moonpool in waves [J]. Chinese Journal of Ship Research, 2020, 15(Supp): 21–34 doi:  10.19693/j.issn.1673-3185.01993
Citation: ZENG Z H, DENG X X, SUN L. effect analysis of thick plate structure on heave motion of ship with moonpool in waves [J]. Chinese Journal of Ship Research, 2020, 15(Supp): 21–34 doi:  10.19693/j.issn.1673-3185.01993
    • 由于海洋资源开发的需求,各种海工船舶或平台应运而生。为了方便作业及避免设备受到风浪流破坏,通常在船体舯部设置自上而下贯穿整个船体的竖井,在此底部开放的竖井内,海水自由流动,其上部为自由液面,称为月池。尽管船体可以屏蔽部分外部风浪流的影响,但是在恶劣的海况下月池内的流体仍会产生较大幅度的运动,不仅如此,月池本身也会对船舶或平台的运动产生不利影响,引起其产生大幅运动,尤其是在垂荡方向,甚至会造成相关功能设备损坏。

      通过理论分析,月池内的流体运动可以归纳为活塞运动和晃荡。活塞运动指月池内液面垂直于船体基平面周期性的上下运动,晃荡则是月池内流体在波流作用下在波流方向上产生的晃荡现象,主要为液面发生的形状变化。

      当外部波浪的激励频率与月池固有频率一致时,月池内流体会产生强烈的共振,而剧烈的活塞运动还会严重加剧船体垂荡运动。迄今,已有许多学者针对月池内流体对船体或平台运动的影响问题进行了研究。例如,孙采薇等[1]基于N-S方程和建立在欧拉网格下的自由液面捕捉法的三维数值波浪水池,对迎浪无航速下含月池船体运动及月池内流体运动进行了研究,发现月池内流体活塞运动会显著加剧船体垂荡运动,而对纵摇运动的影响较小,且在月池内流体与波浪激励发生共振时,船体垂荡运动幅值将大幅增加。黄磊[2]利用Hydrous软件对水平激励下钻井船的矩形月池水动力特性进行了数值模拟,结果表明,当月池内流体质量较大时,月池产生的附加质量极大地抵消了船舶质量,在同样幅度的激励下导致船舶产生大幅运动; 姚震球等[3]利用势流分析软件AQWA对钻井船在波浪下的运动响应进行了数值模拟,发现月池内流体的活塞运动和晃荡对钻井船的纵荡及垂荡运动产生明显影响; 张永恒[4]通过模型实验发现月池对平台垂荡运动水动力系数的影响较大,而对纵摇运动的影响很小。

      鉴于月池内流体运动对船体运动会产生不利影响,许多学者尝试利用一些附加结构或改变月池结构来抑制月池内流体运动。例如,张利军等[5]应用Hydrostar软件计算波浪下钻井船的运动及月池内自由液面的抬升,发现月池内的台阶使月池各阶晃荡产生的共振频率有一定程度地降低。Fukuda[6]通过实验方法研究月池内L型薄板对于单一方向均匀流作用下的月池内流体运动的影响,发现L型薄板靠近水面的位置对月池内流体运动抑制效果最好。然而,该研究并未考虑在月池内增设厚板或较大突起物(不改变月池整体结构)情况下的抑制效果。受文献[6]的启发,本文将含月池的钻井船视为整体来研究,设想在月池内不同高度处增设厚板,旨在抑制船体整体垂荡运动。

      为此,将基于频域势流理论,采用常数边界元方法,以某浮式生产储卸油装置(FPSO)为对象,比较不同浪向下不同长宽比及吃水的月池内各位置的厚板对船体垂荡运动产生的影响,找出厚板相对最优布置位置,以期为后续海工船舶或平台的建造提供参考。

    • 仅考虑船体在波浪作用下的垂荡运动。建立如图1所示的坐标系O-xyz。原点O位于静水面上;x轴位于中纵剖面,由船艉指向船艏;y轴位于中横剖面,由右舷指向左舷;z轴垂直向上。波浪入射角为α0(迎浪时α0 =180°)。

      图  1  船舶数学模型示意图

      Figure 1.  Diagram of mathematical ship model

    • 假设流体为理想流体,流场满足无旋、无黏、不可压缩的基本假定。由流体不可压缩性,推出其体积瞬时膨胀率为零,得到连续性方程:

      $$\nabla {{u}} = 0$$ (1)

      式中,u=(u, v, w),为流体质点速度矢量,其中u, v, w分别为x, y, z方向的流体质点速度。根据流体的无旋假定,流场必然存在势函数Φ(x, y, z, t)。使u=$\nabla \varPhi$,将其代入式(1),得到Laplace方程:

      $${\nabla ^2}\varPhi = 0$$ (2)

      引入线性化假定,自由液面的运动学和动力学统一边界条件为

      $$\frac{{{\partial ^2}\varPhi }}{{\partial {t^2}}} + g\frac{{\partial \varPhi }}{{\partial {\textit{z}}}} = 0,{\rm{ }}{\textit{z}} = 0$$ (3)

      式中:t为时间;g为重力加速度;z为水质点在坐标中的位置,z=0表示水面处。

      流场中物面条件为

      $$\frac{{\partial \varPhi }}{{\partial {{n}}}} = {{V}} \cdot {{n}}$$ (4)

      式中:n为物面法向量;V为物面单元速度矢量,针对固定边界,满足V = 0。

      有势无黏流场中的动水压强p,可由线性化伯努利方程得到:

      $$p = - \rho \frac{{\partial \varPhi }}{{\partial t}}$$ (5)

      式中,ρ为流体密度。

    • 波浪与船舶的作用和船舶运动的水动力特性,采用势流理论方法计算,通过格林函数的映像方法进行求解。为便于研究问题, 将速度势分解为入射势${\phi _0} $、绕射势${\phi _7} $和在6个自由度方向上的运动分量${\xi _j} $j=1$,{\rm{2}}, \cdots ,$6)分解的6个散射势${\phi _j} $

      $$\varPhi = {\rm{Re}} \left[\left({\phi _0} + {\phi _7} + \sum\limits_{j = 1}^6 { - {\rm{i}}\omega {\xi _j}{\phi _j}} \right){{\rm{e}}^{ - {\rm{i}}\omega t}}\right]$$ (6)

      式中,ω为波浪的入射频率。

      相应地,各个散射势的物面条件可写为

      $$ \frac{{\partial {\phi _j}}}{{\partial {{n}}}} = \left\{ {\begin{array}{*{20}{l}} {{{{n}}_j},}&{j = 1,{\rm{2}}, \cdots ,6}\\ { - \dfrac{{\partial {\phi _0}}}{{\partial {{n}}}},}&{j = 7} \end{array}} \right.$$ (7)

      式中,${{n}}_j$为船体在6个运动方向的单位矢量。

      各个散射势可通过下述积分方程计算:

      $$ \begin{split} & \quad\;\;\alpha {\phi }_{j}({{{x}}}_{0})-{ {\iint }_{{S}_{\rm{B}}}\frac{\partial G({{x}},{{{x}}}_{0})}{\partial {{n}}}}{\phi }_{j}({{x}}){\rm{d}}s=\\ & \left\{ \begin{array}{l} -{ \displaystyle{\iint }_{{S}_{\rm{B}}}{{{n}}}_{j}G({{x}},{{{x}}}_{0}){\rm{d}}s,\;\;\;\;\;\;\;\;\; j=1,\rm{2},\cdots ,6}\\ { {\displaystyle\iint }_{{S}_{\rm{B}}}\dfrac{\partial {\phi }_{0}({{x}})}{\partial {{n}}}G({{x}},{{{x}}}_{0}){\rm{d}}s,\;\;\; j=7}\end{array} \right. \end{split}$$ (8)

      式中:α为固角系数,因本文采用常数元,流体域内α=1,流体域外α =0,边界上α =0.5;${S_{\rm{B}}}$为物体在静水中的湿面积;s为湿表面面积;G为满足自由水面条件、海底条件和无穷远处辐射条件的格林函数;j=$1,{\rm{2}}, \cdots ,6$,分别表示物体在6个自由度运动方向做独立受迫运动引起的辐射势,j=7表示物体固定不动时的物面绕射势; x 为场点的坐标矢量; x0 为源点的坐标矢量。式(8)经离散后可得

      $$ {{A}} \, {{{\\text{φ}} _j}} = {{B_j}} ,\quad \quad j = 1,2, \cdots ,7$$ (9)

      式中:A为由格林函数法向导数∂G/∂n求得的系数矩阵;φj为各待求物面单元的散射势$\phi_j $的列向量;Bj为式(8)右侧积分值形成的列向量。通过式(9)可求得物体表面上的绕射势和辐射势。

      物体上的波浪激振力${f_{\rm{W}}}$可通过物体表面上压强的积分求得:

      $$ {f}_{{\rm{W}}}={\rm{i}}\omega \rho { {\iint }_{{S}_{\rm{B}}}({\phi }_{0}+{\phi }_{7}){{n}}{\rm{d}}s}$$ (10)

      辐射势贡献的部分称为水动力系数,即

      $$ {f}_{ij}={\omega }^{2}\rho { {\iint }_{{S}_{\rm{B}}}{\phi }_{j}{{{n}}}_{i}{\rm{d}}s}={\omega }^{2}{a}_{ij}+{\rm{i}}\omega {b}_{ij}$$ (11)

      式中:$f_{ij}$为水动力系数;ni 为船体六个运动方向的单位矢量,与前文计算${\phi }_{j} $ 时的 nj 为不同的法向量;aijbij分别为附加质量系数和辐射阻尼系数;$i,j=1,2, \cdots ,6$,为六自由度方向的索引编号。

      将式(11)代入式(12),可得船体各自由度上的运动响应幅值:

      $$ \sum\limits_{j = 1}^6 {[ - {\omega ^2}({m_{ij}} + {a_{ij}}) - {\rm{i}}\omega ({d_{ij}} + {b_{ij}}) + {c_{ij}}]} {\xi _j} = {f_{{\rm{W}}j}} $$ (12)

      式中,mij为系统质量系数;dij为系统阻尼系数;cij为恢复力系数;fWj为各方向的波浪激振力。

    • 将钻进船月池内流体的活塞运动类比为质量−弹簧系统,Faltinsen等[7]推导了月池活塞运动固有频率的计算公式:

      $$ {\omega _{\rm{p}}} = \sqrt {\frac{g}{h}} $$ (13)

      式中:h为船吃水;${\omega _{\rm{p}}}$为活塞运动的固有频率。

      进一步考虑月池内自由液面的升沉,Molin[8]推导了月池活塞运动固有频率的计算公式:

      $${\omega _{\rm{p}}} = \sqrt {\frac{g}{{T + b{f_0}\left( {\dfrac{b}{l}} \right)}}} $$ (14)

      其中,

      $$\begin{split} & {f_0} = \frac{1}{\text{π}}\left\{ {{\sinh }^{ - 1}}\left( {\frac{l}{b}} \right) + \frac{1}{b}{{\sinh }^{ - 1}}\left( {\frac{b}{l}} \right) + \right.\\&\left. \frac{1}{3}\left( {\frac{b}{l} + \frac{{{l^2}}}{{{b^2}}}} \right) - \frac{1}{3}\left( {1 + \frac{{{l^2}}}{{{b^2}}}} \right)\sqrt {\frac{{{b^2}}}{{{l^2}}} + 1} \right\} \end{split}$$ (15)

      式中:T为月池吃水;l为月池长度;b为月池宽度。

    • 本文以如图2所示的含月池的钻井船为研究对象。该月池面积为256 m2,改变月池长宽比及其吃水,具体工况如表1所示。月池内厚板结构如图3所示,其中厚板高度h1是指钻井船基平面至厚板上表面的距离。

      图  2  不同长宽比月池示意图

      Figure 2.  Diagram of the moonpool with different length-width ratios

      表 1  含月池钻井船的主要参数

      Table 1.  Main parameters of the drillship with moonpool

      参数工况1工况2工况3工况4工况5
      总长/m 112 112 112 112 112
      型宽/m 28.1 28.1 28.1 28.1 28.1
      型深/m 98 98 98 98 98
      吃水/m 6 7.5 9 7.5 7.5
      月池长/m 22.6 22.6 22.6 16 27.7
      月池宽/m 11.3 11.3 11.3 16 9.2
      排水量/t 14 474 18 379 22 248.4 18 379 18 379
      重心垂向位置/m
      (相对船基线)
      9.2 8.2 7.2 8.2 8.2
      厚板厚度/m 1 1 1 1 1
      厚板宽度/m 1.65 1.65 1.65 2.00 1.00
      横摇惯性半径/m 7.9 7.9 8.037 7.9 7.9
      纵摇惯性半径/m 27.864 27.954 28.149 27.954 27.954
      艏摇惯性半径/m 28.899 28.905 29.272 28.905 28.905

      图  3  厚板示意图和月池横剖面图

      Figure 3.  Diagrams of thick plate and cross section of moonpool

    • 计算中,需要验证钻井船的网格收敛性,以确定合适的计算网格。为此,基于目标船体的频域计算结果,验证船体的网格收敛性。选择工况2与工况5验证,设网格尺寸分别为1.2,1.8和2.6 m,结果如图4图5所示。

      图  4  横浪及迎浪下钻井船运动响应(工况2)

      Figure 4.  RAO of the drillship in beam sea and head sea (case 2)

      图  5  横浪及迎浪下钻井船运动响应(工况5)

      Figure 5.  RAO of the drillship in beam sea and head sea (case 5)

      图4图5可知,在选定的3种网格尺寸下,频域计算结果基本一致,尤其对于垂荡运动,计算结果的收敛性好。为了节省计算时间,同时保证足够的精度,网格最大尺寸取为1.8 m。

    • 为验证本文所建立模型及数值方法的有效性,选取2种不同船型作为验证对象,包括无月池船型和有月池船型。

    • 本文研究的钻井船垂荡运动,与FPSO船存在一定的相似性。为了验证本项目选用的AQWA频域数值模型的可靠性,进一步检验该数值模型针对方形系数较大的船体水动力计算的有效性,选取一艘FPSO船型[9]为验证标模1,其主要参数如表2所示,图6分别给出了该船型的俯视图、纵剖视图以及横剖面图。针对该FPSO船型,分别计算其在90°,120°和180°浪向下的运动响应,并与已有实验结果对比,以验证应用AQWA软件建模的可靠性。

      表 2  压载状态下标模1 FPSO船主要参数

      Table 2.  Main parameters of the No.1 FPSO ship in ballast condition

      参数数值
      垂线间长/m285
      船宽/m63
      吃水/m13
      排水量/t225 518
      重心垂向位置/m16.5
      重心纵向位置/m142.5
      横摇惯性半径/m19.45
      纵摇惯性半径/m71.25
      艏摇惯性半径/m71.25

      图  6  标模1 FPSO船体平面视图[9](单位:m)

      Figure 6.  Plane views of the No.1 FPSO ship model [9] (Unit in meter)

      船体模型网格剖分如图7所示,网格基础尺寸为3 m,网格总数为12 202。

      图  7  标模1 FPSO船体网格剖分

      Figure 7.  Mesh generation of the No.1 FPSO ship hull

      在90°,120°和180°浪向时,采用AQWA频域数值模型,计算各自由度下的运动响应幅值(基于船舶重心处)。将计算结果与Nam等[9]的实验数据对比,结果如图8~图10所示。

      图  8  90°浪向下FPSO船运动响应

      Figure 8.  RAO of the FPSO ship when α0=90°

      图  9  120°浪向下FPSO船运动响应

      Figure 9.  RAO of the FPSO ship when α0=120°

      图  10  180°浪向下FPSO船运动响应

      Figure 10.  RAO of the FPSO ship when α0=180°

      图8~图10可见:计算结果与文献[9]实验结果基本吻合,验证了本文采用AQWA软件数值建模及其计算结果的可信度。

    • 为进一步验证所采用数值模型的可靠性,以及该数值模型在含月池船体水动力计算中的有效性,再选取一艘含月池船型[10]为标模2(图11),其主要参数如表3所示。MP1工况对应船模月池长、宽皆为0.2 m,MP2工况对应船模月池长、宽皆为0.4 m,MP3工况对应船模月池长、宽分别为2.0 和0.4 m。VMP/VS为月池容积占全船型体积的百分比,各工况对应值在表3列出。

      表 3  标模2 FPSO船主尺度参数

      Table 3.  Main parameters of the No.2 FPSO ship model

      参数船模尺寸实船尺寸
      船长/m4.0138
      船宽/m0.827.6
      吃水/m0.26.9
      MP1尺度 (VMP/VS = 1.5%)
      月池长度/m0.26.9
      月池宽度/m0.26.9
      月池中心/m(0.65, 0, 0)(22.4, 0, 0)
      船舶质量/t0.56721 627.8
      船舶重心/m(−0.11, 0, 0.055)(−3.8, 0, 1.9)
      惯性半径/m(0.28, 1.02, 0.76)(9.7, 35.2, 26.2)
      MP2尺度 (VMP/VS = 6.4%)
      月池长度/m0.413.8
      月池宽度/m0.413.8
      月池中心/m(0.65, 0, 0)(22.4, 0, 0)
      船舶质量/t0.502 6920 642.3
      船舶重心/m(−0.15, 0, 0.057)(−5.2, 0, 2.0)
      惯性半径/m(0.29, 1.024, 0.75)(10.0, 35.2, 25.9)
      MP3尺度 (VMP/VS = 42.7%)
      月池长度/m2.069
      月池宽度/m0.413.8
      月池中心/m(−0.15, 0, 0)(−5.2, 0, 0)
      船舶质量/t0.374 6915 386.1
      船舶重心/m(−0.08, 0, 0.083)(−2.76, 0, 2.9)
      惯性半径/m(0.30, 1.11, 1.13)(10.4, 38.3, 39.0)

      图  11  标模2平面视图(单位:m)

      Figure 11.  The plane view of No.2 FPSO ship model (Unit in meter)

      针对该FPSO船型,分别计算其在90°和180°浪向下的垂荡运动响应,并与已有的数值模拟和实验结果对比,验证本文应用AQWA软件建模的有效性。

      船体模型网格剖分如图12所示。在90°和180°浪向时,采用AQWA频域数值模型,计算各自由度下的垂荡运动响应幅值(基于船舶重心处)。计算结果与Ravinthrakumar等的WAMIT数值模拟结果和实验数据[10]的对比如图13~图18所示,实验中的波浪条件为波高与波长之比为1/100的非规则波。

      图  12  标模2 FPSO船体网格剖分

      Figure 12.  Mesh generation of the No.2 FPSO ship hull model

      图  13  迎浪下FPSO船垂荡运动响应(MP1)

      Figure 13.  Heaving RAO of the FPSO ship in head sea ( case MP1)

      图  14  横浪下FPSO船垂荡运动响应 (MP1工况)

      Figure 14.  Heaving RAO of the FPSO ship in beam sea (case MP1)

      图  15  迎浪下FPSO船垂荡运动响应(MP2工况)

      Figure 15.  Heaving RAO of the FPSO ship in head sea (case MP2)

      图  16  横浪下FPSO船垂荡运动响应(MP2工况)

      Figure 16.  Heaving RAO of the FPSO ship in beam sea(case MP2)

      图  17  迎浪下FPSO船垂荡运动响应(MP3工况)

      Figure 17.  Heaving RAO of the FPSO ship in head sea(case MP3)

      图  18  横浪下FPSO船垂荡运动响应(MP3工况)

      Figure 18.  Heaving RAO of the FPSO ship in beam sea(case MP3)

      图13~图18可见:采用本文的AQWA模型与文献[10]中的WAMIT数值模拟结果和实验结果相比,变化趋势相同。除共振时运动响应幅值差距较大以外,共振周期相近。从而进一步验证了本文采用AQWA软件进行数值建模的有效性及计算结果的可信度。

    • 为验证本文计算所得频域峰值(由月池共振产生)非由频域数值模型引入的非规则频率,需对计算得到的共振点频率分布数据进行规律性验证。因此,采用AQWA软件计算得到无厚板情况下所有工况对应的月池共振固有频率,并将所得结果与由Molin经验公式(14)计算所得的固有频率值进行比较,以验证数值模拟结果的规律性,结果如表4图19所示。

      表 4  月池固有频率

      Table 4.  Natural frequency of moonpool

      工况数值模拟的固有频率/Hz经验公式计算的固有频率/Hz
      10.1570.160
      20.1470.149
      30.1380.140
      40.1450.148
      50.1480.150

      图  19  月池固有频率对比

      Figure 19.  Comparison on the nature frequency of moonpool

      表4图19可见:数值模拟与经验公式的计算结果相差较小,前者具有规律性,频域峰值不是由频域数值模型引入的非规则频率,进一步证明了AQWA数值模型模拟的可靠性。

    • 为了分析月池内厚板结构对横浪下钻井船垂荡运动的影响,本文运用参数分析法,分别从船体吃水、厚板设置高度、浪向等几个方面系统开展数值模拟分析。

    • 在工况1、工况2和工况3下,分别在钻井船中设置沿船长方向对称的厚板,改变厚板设置高度,对横浪下船体的垂荡运动进行数值模拟,结果如图20~图22所示,数据 汇总如表5~表7

      表 5  横浪下钻井船垂荡运动响应(工况1)

      Table 5.  Heaving RAO of the drillship in beam sea (case 1)

      h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m−1)
      无厚板0.1571.201
      1.00.1471.406
      2.00.1481.363
      3.00.1481.335
      4.00.1491.309
      5.00.1481.302
      5.20.1471.334
      5.40.1461.339
      5.60.1461.347
      5.70.1451.375
      5.80.1401.502
      5.90.1281.626
      6.00.1800.865

      表 6  横浪下钻井船船体垂荡运动响应(工况2)

      Table 6.  Heaving RAO of drillship in beam sea (case 2)

      h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m−1)
      无厚板 0.147 1.396
      1.0 0.138 1.727
      2.0 0.139 1.670
      3.0 0.139 1.627
      4.0 0.140 1.582
      5.0 0.140 1.569
      6.0 0.140 1.539
      6.5 0.140 1.539
      6.7 0.139 1.557
      6.9 0.138 1.592
      7.1 0.138 1.591
      7.2 0.137 1.625
      7.3 0.133 1.790
      7.4 0.124 1.895
      7.5 0.202 0.405

      表 7  横浪下钻进船垂荡运动响应(工况3)

      Table 7.  Heaving RAO of the drillship in beam sea (case 3)

      h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m−1)
      无厚板0.1381.701
      2.00.1322.011
      4.00.1321.898
      6.00.1331.860
      8.00.1331.798
      8.20.1321.833
      8.40.1311.851
      8.60.1321.806
      8.70.1311.852
      8.80.1282.105
      8.90.1202.268
      9.00.1600.981

      图  20  横浪及不同厚板高度下钻井船垂荡运动响应(工况1)

      Figure 20.  Heaving RAO of the drillship in beam sea with different heights of thick plate(case 1)

      图  21  横浪及不同厚板高度下钻井船垂荡运动响应(工况2)

      Figure 21.  Heaving RAO of the drillship in beam sea with different heights of thick plate(case 2)

      图  22  横浪及不同厚板高度下钻井船垂荡运动响应(工况3)

      Figure 22.  Heaving RAO of the drillship in beam sea with different heights of thick plate(case 3)

      图20~图22表5~表7可得如下结论:

      1)当厚板上表面与静水面的距离大于1 m且在此高度范围内改变高度位置时,月池的固有频率与月池共振引起的船体垂荡幅值响应算子(RAO)的变化不大;厚板会使月池垂荡的固有频率降低,但也加剧了共振引起的船体垂荡运动。例如,工况1中厚板远离静水面时,月池共振频率在0.148 Hz附近,船体垂荡幅值增加0.1 m/m以上。

      2)当厚板靠近静水面,即厚板上表面与静水面的距离不大于1 m且厚板上表面不高于静水面,在此高度范围内改变厚板设置高度时,月池升沉固有频率和船体垂荡幅值呈规律性变化:在厚板上表面与静水面齐平之前,随着厚板设置高度的增加,月池升沉的固有频率逐渐靠近船体固有频率,船体垂荡幅值也随之增加。例如,工况2中月池固有频率由0.140 Hz减至0.124 Hz,而船体垂荡幅值由1.539 m/m递增至1.895 m/m。

      3)当厚板上表面与静水面齐平时,月池升沉固有频率变大,月池共振频率范围变小,因月池共振引起的船体垂荡幅值大幅减小,而因船体共振产生的垂荡幅值略微增大。例如,工况2中当厚板上表面与静水面齐平时,月池共振频率跃迁至0.160 Hz,远大于无厚板时的0.138 Hz;而船体垂荡幅值锐减至0.981 m/m,远小于无厚板时的1.701 m/m。工况1和工况3也有同样的现象,可见将厚板布置于其上表面正好与静水面齐平时,对因月池共振引起的垂荡运动的抑制效果最好,且并不会对由船体共振产生的垂荡运动产生太大影响。

      综上所述,以上变化规律出现的原因是:月池共振引起的垂荡运动的峰(下称“月池共振峰”),随着厚板高度的上升,向船体共振引起的垂荡运动的峰(下称“船体共振峰”)靠拢,并逐渐融合。曲线与坐标轴围成的面积代表能量,月池共振峰的面积即代表月池共振时的能量,船体共振峰同理。月池共振峰向船体共振峰靠拢、融合的过程,可以视为能量转移的过程:月池共振的能量随着厚板高度越来越接近水面,部分能量被船体吸收,剩余的能量引起了较小的垂荡运动。而月池共振频率的变化,可能是因为厚板的存在,改变了月池的结构形式,也改变了月池内流体流动形式,致使月池固有频率发生规律性变化。

    • 分别在工况4和工况5的钻井船中设置沿船长方向对称的厚板,改变厚板设置的高度,分别对横浪下船体的垂荡运动进行数值模拟计算,结果如图23图24所示,数据汇总见表8表9

      表 8  横浪下钻井船垂荡运动响应(工况4)

      Table 8.  Heaving RAO of the drillship in beam sea (case 4)

      h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m-1)
      无厚板0.1451.301
      1.50.1381.614
      3.00.1391.912
      4.50.1391.476
      6.00.1391.429
      7.50.1800.613

      表 9  横浪下钻井船船体垂荡运动响应(工况5)

      Table 9.  Heaving RAO of the drillship in beam sea (case 5)

      h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m-1)
      无厚板0.1481.452
      1.50.1441.583
      3.00.1441.547
      4.50.1451.505
      6.00.1451.469
      7.50.2010.543

      图  23  横浪及不同厚板高度下钻井船垂荡运动响应(工况4)

      Figure 23.  Heaving RAO of the drillship in beam sea with different heights of thick plate (case 4)

      图  24  横浪及不同厚板高度下钻井船垂荡运动响应(工况5)

      Figure 24.  Heaving RAO of the drillship in beam sea with different heights of thick plate (case 5)

      图23~图24表8~表9可见:对于不同长宽比的月池结构,改变厚板所处位置的高度,各工况下月池升沉的固有频率与月池共振引起的船体垂荡幅值具有相近的变化规律。例如,工况5中厚板布置于水面时,月池固有频率由无厚板时的0.148 Hz跃迁至0.201 Hz;船体垂荡幅值由无厚板时的1.452 m/m锐减至0.543 m/m。将厚板布置于其上表面与静水面齐平处,对月池共振产生的垂荡运动抑制效果最好。

    • 在工况2对应月池中设置高度位置不同的厚板,改变浪向(即0°,45°,135°和180°),分别计算波浪下钻井船的垂荡运动响应,对果如图25~图28所示。

      图  25  0°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

      Figure 25.  Heaving RAO of the drillship with different heights of thick plate when α0= 0° (case 2)

      图  26  45°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

      Figure 26.  Heaving RAO of the drillship with different heights of thick plate when α0= 45° (case 2)

      图  27  135°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

      Figure 27.  Heaving RAO of the drillship with different heights of thick plate when α0=135° (case 2)

      图  28  180°浪向不同厚板高度下钻井船垂荡运动响应(工况2)

      Figure 28.  Heaving RAO of the drillship with different heights of thick plate when α0=180° (case 2)

      图25~图28可见:不同浪向下,将厚板布置于其上表面与水面齐平处对月池共振引起的船体垂荡运动仍有较好的抑制效果,即垂荡幅值明显降低,共振频率向高频移动。

    • 在工况2对应月池中设置厚板,改变厚板的高度位置,令厚板上表面高于静水面,但厚板下表面仍处于静水面以下,计算横浪下钻井船的垂荡运动响应,结果如图图29所示,数据汇总见表10

      表 10  横浪下厚板上表面高于静水面时钻井船垂荡运动响应(工况2)

      Table 10.  Heaving RAO of the drillship with thick plate above the waterline in beam sea (case 2)

      h1/m月池升沉固有频率/Hz船体垂荡幅值/(m·m-1)
      无厚板0.1471.396
      7.50.2020.405
      7.60.1690.970
      7.80.1680.934
      8.00.1690.969

      图  29  横浪下厚板上表面高于静水面时钻井船垂荡运动响应(工况2)

      Figure 29.  Heaving RAO of the drillship with thick plate above the waterline in beam sea(case 2)

      图29表10可见:当厚板上表面超过静水面并继续提高厚板位置,月池固有频率从0.202 Hz减少至0.168和0.169 Hz,船体垂荡幅值由0.405 m/m升至0.900 m/m以上。可见,当厚板上表面稍微高出静水面,厚板对于月池共振引起的船体垂荡运动的抑制作用减弱了许多,且月池共振频率重新向无厚板时的月池垂荡共振频率靠近。由此表明,将厚板布置于其上表面与静水面齐平的位置,对船船体垂荡运动的抑制效果仍是最佳的。

    • 本文基于势流理论的数值模拟和结果对比,对含厚板的月池结构船体在横浪作用下的垂荡运动进行了研究,得到如下结论:

      1) 当厚板越接近水面处,月池共振频率就越小,船体垂荡幅值呈规律性增加。

      2) 当厚板上表面与静水面齐平时,船体的垂荡幅值突变至小值,此时厚板对船体垂荡运动的抑制效果最佳。厚板上表面高于静水面后,上述抑制作用减弱。

      本文研究结果,仅为定性分析结果,还需要更为广泛的普适性以及定量分析。后期将采用CFD软件,考虑流体黏性和流场的漩涡作用,对月池内的流场进行数值模拟,以观察该流场变化的情况,进一步探究导致上述变化的原因。

参考文献 (10)

目录

    /

    返回文章
    返回