留言板

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

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

航行姿态对半滑行三体船型静水阻力影响的数值研究

贾敬蓓 宗智 金国庆 王海英

贾敬蓓, 宗智, 金国庆, 等. 航行姿态对半滑行三体船型静水阻力影响的数值研究[J]. 中国舰船研究, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859
引用本文: 贾敬蓓, 宗智, 金国庆, 等. 航行姿态对半滑行三体船型静水阻力影响的数值研究[J]. 中国舰船研究, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859
JIA J B, ZONG Z, JIN G Q, et al. A numerical investigation of sinkage and trim effects on the resistance of trimaran hull form in calm water[J]. Chinese Journal of Ship Research, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859
Citation: JIA J B, ZONG Z, JIN G Q, et al. A numerical investigation of sinkage and trim effects on the resistance of trimaran hull form in calm water[J]. Chinese Journal of Ship Research, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859

航行姿态对半滑行三体船型静水阻力影响的数值研究

doi: 10.19693/j.issn.1673-3185.01859
基金项目: 国家自然科学基金面上项目资助(51679037);国家自然科学基金重点项目资助(51639003);国家重点基础研究发展计划资助(2013CB036101)
详细信息
    作者简介:

    贾敬蓓,女,1984年生,博士,讲师。研究方向:船舶水动力学。E-mail:jiajingbei@dlou.edu.cn

    宗智,男,1964年生,博士,教授,博士生导师。研究方向:高性能船舶,船舶水动力学。E-mail:zongzhi@dlut.edu.cn

    金国庆,男,1994年生,博士生。研究方向:船舶水动力学,计算流体力学,离散涡方法。E-mail:jinguoqing2012@mail.dlut.edu.cn

    王海英,女,1980年生,博士,讲师。研究方向:船舶振动与噪声机理研究。E-mail:wanghaiying@dlou.edu.cn

    通讯作者:

    宗智

  • 中图分类号: U661.31+1

A numerical investigation of sinkage and trim effects on the resistance of trimaran hull form in calm water

  • 摘要:   目的  半滑行船型(弗劳德数Fr >0.4)的航行姿态(升沉和纵倾)对船舶阻力性能影响显著,重点探讨半滑行状态下航行姿态对三体船型阻力成分(摩擦阻力和压差阻力)的影响规律。  方法  采用k-ε湍流模型,使用重叠网格技术,对三体船型的静水阻力问题进行CFD数值模拟研究。计算三体船的2种状态(自由航行状态和固定航行姿态的状态),比较上述状态下的阻力成分,分析航行姿态对阻力成分的影响规律。  结果  计算结果表明,在半滑行状态下,由航行姿态引起的总阻力变化超过20%;航行姿态变化对压差阻力影响显著,约高达50%,而对摩擦阻力的影响则相对较小,不到8.5%;由航行姿态引起的阻力增加,90%以上来自压差阻力的变化。  结论  准确预报航行姿态及由此产生的压差阻力变化对半滑行三体船型非常重要。
  • 图  1  三体船模型横剖面示意图

    Figure  1.  Schematic diagram of transverse section of trimaran model

    图  2  三体船布局示意图

    Figure  2.  Schematic diagram of trimaran layout

    图  3  边界条件设置

    Figure  3.  The boundary conditions

    图  4  不同时间步长下的总阻力时历曲线

    Figure  4.  Time history curves of total resistance for different time steps

    图  5  不同网格模型的总阻力时历曲线

    Figure  5.  Time history curves of total resistance for different grid sizes

    图  6  约束模型总阻力与试验值的对比

    Figure  6.  Comparison between total resistance and experimental values of constrained model

    图  7  艏、艉升沉测量点位置示意图

    Figure  7.  The position diagram of sinkage measurement points at bow and stern

    图  8  三体船模型艏、艉测量点的升沉曲线

    Figure  8.  The sinkage curves of measurement points at bow and stern of trimaran model

    图  9  三体船模型的纵倾曲线

    Figure  9.  The trim curves of trimaran model

    图  10  自由模型总阻力与试验值的对比

    Figure  10.  Comparison of total resistance and experimental values of free model

    图  11  自由模型的摩擦阻力和压差阻力曲线

    Figure  11.  The frictional resistance and pressure resistance curves of free model

    图  12  自由模型摩擦阻力和压差阻力百分比曲线

    Figure  12.  The percentage curves of frictional resistance and pressure resistance of free model

    图  13  自由模型与约束模型的总阻力曲线

    Figure  13.  The total resistance curves of free model and constrained model

    图  14  自由模型的总阻力变化率曲线

    Figure  14.  The curve of total resistance change rate of free model

    图  15  自由模型和约束模型的摩擦阻力曲线

    Figure  15.  The frictional resistance curves of free model and constrained model

    图  16  自由模型的摩擦阻力变化率曲线

    Figure  16.  The curve of frictional resistance change rate of free model

    图  17  自由模型的摩擦阻力改变量曲线

    Figure  17.  The curves of frictional resistance change of free model

    图  18  高航速时的试验图片和计算图片(Fr=0.9)

    Figure  18.  The experimental picture and computational contours at high speed(Fr=0.9)

    图  19  自由模型和约束模型的压差阻力曲线

    Figure  19.  The pressure resistance curves of free model and constrained model

    图  20  自由模型的压差阻力变化率曲线

    Figure  20.  The curve of pressure resistance change rate of free model

    图  21  自由模型摩擦阻力和压差阻力的变化量

    Figure  21.  The change of friction resistance and pressure resistance of free model

    表  1  船模主尺度参数

    Table  1.   Main dimensions of ship model

    参数主船体侧体
    船长/m4.01.0
    设计水线宽/m0.358 40.085
    设计吃水/m0.170.10
    排水量/kg129.074.45
    湿表面积/m21.8990.201
    方形系数0.5300.524
    棱形系数0.7690.746
    下载: 导出CSV

    表  2  不同时间步长下的计算结果

    Table  2.   The calculation results for different time steps

    时间步长Δt/s总阻力Rt/N偏差/%
    0.0192.090
    0.0291.98−0.1
    0.0493.14 1.1
    下载: 导出CSV

    表  3  模型网格信息和计算结果

    Table  3.   The grids and calculation results

    网格数总阻力Rt/N偏差/%
    150万91.25−1.0
    240万91.98−0.2
    392万92.190
    下载: 导出CSV
  • 姜次平, 邵世明. 船舶阻力[M]. 上海: 上海交通大学出版社, 1985: 32-36.

    JIANG C P, SHAO S R. Ship resistance[M]. Shanghai: Shanghai Jiao Tong University Press, 1985: 32-36 (in Chinese).
    JACKSON J. Fast trimaran ferry for islands run "Benchijigua Express"[J]. Engineers Australia, 2015, 77(2): 42–43.
    卢晓平, 王中. 美国海军新舰艇研发与采办[M]. 北京: 电子工业出版社, 2010: 55-56.

    LU X P, WANG Z. R&D and procurement of U.S. navy ships[M]. Beijing: Publishing House of Electronics Industry, 2010: 55-56 (in Chinese).
    刘鑫旺, 詹开宇, 万德成. 势流理论评估下三体船主体船型及片体间距优化[C]//第十九届中国海洋(岸)工程学术讨论会论文集. 重庆, 2019: 203-208.

    LIU X W, ZHAN K Y, WAN D C. Optimization of the main ship form and the space between sidehulls based on potential flow theory[C]//Proceedings of the 19th China Ocean (Shore) Engineering Symposium. Chongqing, 2019: 203-208 (in Chinese).
    戴原星, 张志远, 刘建国, 等. 喷水推进三体船阻力与自航数值模拟研究[J]. 船舶, 2019, 30(1): 105–111. doi:  10.3969/j.issn.1001-9855.2019.01.019

    DAI Y X, ZHANG Z Y, LIU J G, et al. Numerical simulation of resistance and self-propulsion performance for water-jet propelled trimaran[J]. Ship & Boat, 2019, 30(1): 105–111 (in Chinese). doi:  10.3969/j.issn.1001-9855.2019.01.019
    张明霞, 韩兵兵, 卢鹏程, 等. 网格因素及湍流模型对三体船阻力计算影响分析[J]. 武汉理工大学学报(交通科学与工程版), 2019, 43(5): 876–881.

    ZHANG M X, HAN B B, LU P C, et al. Impact analysis of mesh factors and turbulence models on resistance calculation of trimaran[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2019, 43(5): 876–881 (in Chinese).
    ZONG Z, HONG Z C, WANG Y G, et al. Hull form optimization of trimaran using self-blending method[J]. Applied Ocean Research, 2018, 80: 240–247. doi:  10.1016/j.apor.2018.09.003
    崔焰, 卢晓平. 基于线性兴波阻力理论的三体船中体船型优化[J]. 船舶工程, 2018, 40(1): 32–36.

    CUI Y, LU X P. Trimaran center hull optimization based on linear wave resistance theory[J]. Ship Engineering, 2018, 40(1): 32–36 (in Chinese).
    YANG C, LÖHNER R. Calculation of ship sinkage and trim using a finite element method and unstructured grids[J]. International Journal of Computational Fluid Dynamics, 2002, 16(3): 217–227. doi:  10.1080/10618560290034690
    王中, 卢晓平, 付攀. 三体船升沉和纵倾计算及其对兴波阻力的影响[J]. 上海交通大学学报, 2010, 44(10): 1388–1392.

    WANG Z, LU X P, FU P. Calculation of the sinkage and trim of trimaran and their effect on wave making resistance[J]. Journal of Shanghai Jiaotong University, 2010, 44(10): 1388–1392 (in Chinese).
    倪崇本, 朱仁传, 缪国平, 等. 计及航行姿态变化的高速多体船阻力预报[J]. 水动力学研究与进展(A辑), 2011, 26(1): 101–107.

    NI C B, ZHU R C, MIAO G P, et al. The resistance prediction for high speed multi-hull vessels with consideration of hull gesture variation during voyage[J]. Chinese Journal of Hydrodynamics (Ser.A), 2011, 26(1): 101–107 (in Chinese).
    DENG R, LI C, HUANG D B, et al. The effect of trimming and sinkage on the trimaran resistance calculation[J]. Procedia Engineering, 2015, 126: 327–331. doi:  10.1016/j.proeng.2015.11.199
    周广利, 艾子涛, 邓锐, 等. 高速三体船的阻力预报方法研究[J]. 船舶力学, 2016, 20(7): 805–815. doi:  10.3969/j.issn.1007-7294.2016.07.003

    ZHOU G L, AI Z T, DENG R, et al. Resistance prediction method research of high speed trimaran[J]. Journal of Ship Mechanics, 2016, 20(7): 805–815 (in Chinese). doi:  10.3969/j.issn.1007-7294.2016.07.003
    YEUNG R W. Sinkage and trim in first-order thin-ship theory[J]. Journal of Ship Research, 1972, 16(1): 47–59.
    李柯, 宗智, 孙雷. 三体船粘压阻力预报方法[J]. 中国舰船研究, 2013, 8(4): 27–35. doi:  10.3969/j.issn.1673-3185.2013.04.005

    LI K, ZONG Z, SUN L. Estimation method of the viscous resistance for trimarans[J]. Chinese Journal of Ship Research, 2013, 8(4): 27–35 (in Chinese). doi:  10.3969/j.issn.1673-3185.2013.04.005
    VERSTEEG H K, MALALASEKERA W. An introduction to computational fluid dynamics: the finite volume method[M]. England: Longman Scientific & Technical, 1995.
    王福军. 计算流体动力学分析——CFD软件原理与应用[M]. 北京: 清华大学出版社, 2004: 114-122.

    WANG F J. Computational fluid dynamics analysis-principles and applications of CFD software[M]. Beijing: Tsinghua University Press, 2004: 114-122 (in Chinese).
    贾敬蓓. 三体船型阻力和运动性能试验研究[D]. 大连: 大连理工大学, 2010.

    JIA J B. Experimental study of resistance and motion characteristics of trimaran hull form[D]. Dalian: Dalian University of Technology, 2010 (in Chinese).
  • [1] 周渝航, 冯国庆, 张旭辉2020: 三体船屈服强度的可靠性研究 , 中国舰船研究, 15, 117-122. doi: 10.19693/j.issn.1673-3185.01546
    [2] 王硕, 宗智2019: 基于蒙特卡洛方法的三体船侧体运动响应优化 , 中国舰船研究, 14, 122-128. doi: 10.19693/j.issn.1673-3185.01348
    [3] 孙源, 卢晓平, 李井煜, 王中2019: 滑行艇阻力计算方法对比研究 , 中国舰船研究, 14, 27-32. doi: 10.19693/j.issn.1673-3185.01277
    [4] 张明霞, 韩兵兵, 卢鹏程2019: 小水线面三体船与细长型三体船剩余阻力对比分析 , 中国舰船研究, 14, 21-29. doi: 10.19693/j.issn.1673-3185.01233
    [5] 邓琦, 毛筱菲, 吴铭浩2016: 主体瘦长度对三体船耐波性和波浪载荷的影响 , 中国舰船研究, 11, 8-14. doi: 10.3969/j.issn.1673-3185.2016.06.002
    [6] 魏成柱, 毛立夫, 李英辉, 易宏2015: 单体半滑行穿浪船船型与静水航行性能 , 中国舰船研究, 10, 16-21. doi: 10.3969/j.issn.1673-3185.2015.05.003
    [7] 李柯, 宗智, 孙雷2013: 三体船粘压阻力预报方法 , 中国舰船研究, 8, 27-35. doi: 10.3969/j.issn.1673-3185.2013.04.005
    [8] 黄卫刚, 姜治芳, 邱辽原2012: 撞角球鼻艏形式和参数变化对阻力影响的研究 , 中国舰船研究, 7, 23-28. doi: 10.3969/j.issn.1673-3185.2012.01.005
    [9] 张乔斌, 尹成彬, 吴开峰2012: 滑行艇阻力近似计算方法对比研究 , 中国舰船研究, 7, 25-29. doi: 10.3969/j.issn.1673-3185.2012.03.005
    [10] 王毅, 卢晓平, 赵军强, 段晔鑫2012: 以Rankine源三维面元法求解三体船纵摇与升沉运动 , 中国舰船研究, 7, 29-36. doi: 10.3969/j.issn.1673-3185.2012.02.006
    [11] 王文江, 宗智, 倪少玲, 张磊, 陈志强2012: 半滑行船扰流板阻力试验研究 , 中国舰船研究, 7, 18-22. doi: 10.3969/j.issn.1673-3185.2012.01.004
    [12] 崔健, 文逸彦, 陈鹏, 杨松林2012: 一种具有倾斜侧体的三体船阻力试验研究 , 中国舰船研究, 7, 57-62. doi: 10.3969/j.issn.1673-3185.2012.06.009
    [13] 邱辽原, 谢伟, 姜治芳, 冯佰威, 刘祖源2011: 基于参数化CAD模型的船型阻力/耐波性一体化设计 , 中国舰船研究, 6, 18-21,29. doi: 10.3969/j.issn.1673-3185.2011.01.004
    [14] 卢晓平, 姚迪, 王中2010: 三体船操纵特性计算机数值仿真 , 中国舰船研究, 5, 1-7. doi: 10.3969/j.issn.1673-3185.2010.03.001
    [15] 黄晓琼, 陈立, 杨雄辉, 周心桃2009: 三体船连接桥结构波浪载荷研究 , 中国舰船研究, 4, 42-46. doi: 10.3969/j.issn.1673-3185.2009.04.009
    [16] 黄晓琼, 陈立, 杨雄辉, 周心桃2009: 三体船耐波性预报 , 中国舰船研究, 4, 42-47. doi: 10.3969/j.issn.1673-3185.2009.06.010
    [17] 岳国强, 姚朝帮, 董文才2009: 深V型滑行艇静水阻力性能影响因素研究 , 中国舰船研究, 4, 24-27,37. doi: 10.3969/j.issn.1673-3185.2009.03.005
    [18] 刘昌明, 李云波, 李裕龙2009: 三体船兴波阻力计算方法比较及兴波干扰研究 , 中国舰船研究, 4, 28-32. doi: 10.3969/j.issn.1673-3185.2009.05.006
    [19] 黄晓琼, 陈立, 杨雄辉, 周心桃2008: 三体船波浪诱导载荷计算研究 , 中国舰船研究, 3, 16-20. doi: 10.3969/j.issn.1673-3185.2008.05.004
    [20] 孙荣, 吴晓光, 姜治芳, 吴启锐, 卢晓晖2008: 三体船阻力数值计算及方案优选 , 中国舰船研究, 3, 21-24. doi: 10.3969/j.issn.1673-3185.2008.05.005
  • 加载中
图(21) / 表 (3)
计量
  • 文章访问数:  29
  • HTML全文浏览量:  10
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-20
  • 修回日期:  2020-04-22
  • 网络出版日期:  2020-11-04
  • 刊出日期:  2020-12-30

航行姿态对半滑行三体船型静水阻力影响的数值研究

doi: 10.19693/j.issn.1673-3185.01859
    基金项目:  国家自然科学基金面上项目资助(51679037);国家自然科学基金重点项目资助(51639003);国家重点基础研究发展计划资助(2013CB036101)
    作者简介:

    贾敬蓓,女,1984年生,博士,讲师。研究方向:船舶水动力学。E-mail:jiajingbei@dlou.edu.cn

    宗智,男,1964年生,博士,教授,博士生导师。研究方向:高性能船舶,船舶水动力学。E-mail:zongzhi@dlut.edu.cn

    金国庆,男,1994年生,博士生。研究方向:船舶水动力学,计算流体力学,离散涡方法。E-mail:jinguoqing2012@mail.dlut.edu.cn

    王海英,女,1980年生,博士,讲师。研究方向:船舶振动与噪声机理研究。E-mail:wanghaiying@dlou.edu.cn

    通讯作者: 宗智
  • 中图分类号: U661.31+1

摘要:   目的  半滑行船型(弗劳德数Fr >0.4)的航行姿态(升沉和纵倾)对船舶阻力性能影响显著,重点探讨半滑行状态下航行姿态对三体船型阻力成分(摩擦阻力和压差阻力)的影响规律。  方法  采用k-ε湍流模型,使用重叠网格技术,对三体船型的静水阻力问题进行CFD数值模拟研究。计算三体船的2种状态(自由航行状态和固定航行姿态的状态),比较上述状态下的阻力成分,分析航行姿态对阻力成分的影响规律。  结果  计算结果表明,在半滑行状态下,由航行姿态引起的总阻力变化超过20%;航行姿态变化对压差阻力影响显著,约高达50%,而对摩擦阻力的影响则相对较小,不到8.5%;由航行姿态引起的阻力增加,90%以上来自压差阻力的变化。  结论  准确预报航行姿态及由此产生的压差阻力变化对半滑行三体船型非常重要。

English Abstract

贾敬蓓, 宗智, 金国庆, 等. 航行姿态对半滑行三体船型静水阻力影响的数值研究[J]. 中国舰船研究, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859
引用本文: 贾敬蓓, 宗智, 金国庆, 等. 航行姿态对半滑行三体船型静水阻力影响的数值研究[J]. 中国舰船研究, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859
JIA J B, ZONG Z, JIN G Q, et al. A numerical investigation of sinkage and trim effects on the resistance of trimaran hull form in calm water[J]. Chinese Journal of Ship Research, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859
Citation: JIA J B, ZONG Z, JIN G Q, et al. A numerical investigation of sinkage and trim effects on the resistance of trimaran hull form in calm water[J]. Chinese Journal of Ship Research, 2020, 15(6): 106–114 doi:  10.19693/j.issn.1673-3185.01859
    • 对于压差阻力占总阻力比例较大的高速船舶,细长的船型可以显著降低阻力[1]。但是,若船体过细,其横稳性会变差。而在主船体两侧各设置一个小的侧体,则可在提高航速的同时大大改进横稳性。近20多年来,备受关注的三体船概念[2-3]就充分考虑了这一因素[4-8],三体船的应用带来了船型设计上的变革。澳大利亚Austal公司的Benchijigua Express三体渡船(船长L=126.7 m)和美国海军“独立”号三体濒海战斗舰(船长L=127.6 m)的投入使用,激起了人们对大型高速三体船研究的兴趣[2-3]。Benchijigua Express三体渡船和“独立”号濒海战斗舰的速度均高达40 kn以上,长度超过100 m,其显著特征是弗劳德数Fr > 0.6,处于半滑行状态。与排水型船型(一般Fr < 0.3)或滑行艇型(一般Fr > 1.0)相比,半滑行状态由于动升力的作用,静浮力和浮力分布发生了变化,致使船体的航行姿态(升沉和纵倾)发生了显著变化,进而湿表面积和兴波状况也产生了相应的变化,带来的效果就是阻力性能发生显著变化[9-13]。半滑行船型的阻力性能不同于排水型船型和滑行艇型,对于高速半滑行三体船型,其航行姿态对阻力的影响是阻力预报方面必不可少的一个重要研究内容。

      造船界很早就意识到了航行姿态对阻力的影响。Yeung[14]最早从船舶兴波阻力薄船理论的角度开展了升沉和纵倾对阻力影响的研究。航行姿态对三体船阻力的影响,始于Yang和Löhner[9]、王中等[10]及李柯等[15]的研究,不同于Yeung的早期薄船理论,他们分别使用更加准确的Rankine源法或细长体理论考察了这一问题。数值结果表明,在0.3<Fr<0.6范围内,航行姿态的影响显著。同时,结果也表明,势流理论在更高速时与试验结果存在较大偏差,其原因是势流理论在数学模型和数值模拟等方面均存在误差,因而会产生一些不真实的现象。

      随着计算机技术和船舶水动力学的迅猛发展,近年来,较多学者开始倾向于采用CFD方法考察这方面的问题[11-13]。本文将以一艘方艉折角线型的前三体船为例,采用CFD方法重点探讨航行姿态对阻力成分(摩擦阻力和压差阻力)及不同成分阻力大小的影响规律。

    • 浮力的本质是压力沿船体表面的积分。根据伯努利方程,当航速较大时,压力就会降低,从而沿船体表面积分的浮力也会发生变化。对于排水型船,当航速较低时,航行中压力改变会很小,故浮力的变化通常可忽略不计。因此,排水型船(体积弗劳德数满足${{F}}{{{r}}_\nabla }{\rm{ < 1}}{\rm{.0}}$,一般对应的长度弗劳德数满足Fr < 0.3)的航行姿态影响通常忽略不计;而对于滑行船型(体积弗劳德数满足${{F}}{{{r}}_\nabla }{\rm{ > 3}}{\rm{.0}}$,对应的长度弗劳德数Fr >1.0),因升力是主要的托举力,故航行姿态由升力决定。对于半滑行船型(体积弗劳德数满足$ 1.0 < {Fr}_\nabla < 3.0 $,一般对应的长度弗劳德数满足0.3~0.4< Fr <1.0),由于航速变化带来的压力变化、进而浮力以及航行姿态的变化相当显著,因此必须予以考虑,这也是航行姿态对半滑行船型(或者过渡船型)的阻力预报影响较大的原因。三体船型属高速船型,处于半滑行状态。为分析航行姿态对阻力的影响,本文以某三体船为例,考虑三体船的2种升沉和纵倾状态,即自由航行状态(释放三体船的升沉和纵倾运动,简称“自由模型”)和固定航行姿态的状态(约束三体船的升沉和纵倾运动,简称“约束模型”),采用CFD方法对各自的阻力进行了计算。

    • 绕三体船的流动应满足连续性方程和动量方程(N-S方程)。假定流体为不可压缩流体,流体域内时均形式的连续性方程和动量方程用张量形式表示为[16-17]

      $$\frac{{\partial {u_i}}}{{\partial {x_i}}} = 0$$ (1)
      $$\rho \frac{{\partial {u_i}}}{{\partial t}} + \rho \frac{{\partial ({u_i}{u_j})}}{{\partial {x_j}}} = - \frac{{\partial p}}{{\partial {x_i}}} + \frac{\partial }{{\partial {x_j}}}\left[ {\mu \frac{{\partial {u_i}}}{{\partial {x_j}}} - \rho \overline {u_i'u_j'} } \right] + \rho {g_i}$$ (2)

      式中:${u_i},{u_j}$为速度分量的时均值($i,j = 1,2,3$);$t$为时间;${x_i},{x_j}$为三维直角坐标系的3个坐标轴($i,j = 1,2,3$);$p$为压力时均值;$\rho $为流体密度;$\mu $为流体动力黏性系数;${g_i}$为重力加速度分量($i = 1,2,3$);$\rho \overline {u_i'u_j'} $为雷诺应力项。

    • 本文采用标准$k - \varepsilon $湍流模型来封闭RANS方程。湍动能$k$、湍动耗散率$\varepsilon $的输运方程分别如式(3)和式(4)所示[16-17]

      $$\begin{split} & \frac{{\partial k}}{{\partial t}} + {u_j}\frac{{\partial k}}{{\partial {x_j}}} = \frac{\partial }{{\partial {x_j}}}\left[ {\left( {\frac{{{\nu _{\rm{t}}}}}{{{\sigma _k}}} + \nu } \right)\frac{{\partial k}}{{\partial {x_j}}}} \right] +\\& \qquad {\nu _{\rm{t}}}\left( {\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{{\partial {u_j}}}{{\partial {x_i}}}} \right)\frac{{\partial {u_i}}}{{\partial {x_j}}} - \varepsilon \end{split} $$ (3)
      $$\begin{split} & \frac{{\partial \varepsilon }}{{\partial t}} + {u_j}\frac{{\partial \varepsilon }}{{\partial {x_j}}} = \frac{\partial }{{\partial {x_j}}}\left[ {\left( {\frac{{{\nu _{\rm{t}}}}}{{{\sigma _\varepsilon }}} + \nu } \right)\frac{{\partial \varepsilon }}{{\partial {x_j}}}} \right] + \\&{C_{\varepsilon 1}}{\nu _{\rm{t}}}\frac{\varepsilon }{k}\left( {\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{{\partial {u_j}}}{{\partial {x_i}}}} \right)\frac{{\partial {u_i}}}{{\partial {x_j}}} - {C_{\varepsilon 2}}\frac{{{\varepsilon ^2}}}{k} \end{split}$$ (4)

      式中:${\nu _{\rm{t}}}$为湍动黏度,定义为${\nu _{\rm{t}}} = \rho {C_\mu }\dfrac{{{k^2}}}{\varepsilon }$$\nu $为流体运动黏性系数。通常情况下,标准$k - \varepsilon $湍流模型的系数取值如下[16-17]

      $${C_\mu } = 0.09,\;{C_{_{\varepsilon 1}}} = 1.44,\;{C_{_{\varepsilon 2}}} = 1.92,\;{\sigma _k} = 1.0,{\sigma _\varepsilon } = 1.3$$ (5)
    • 本次计算采用的是一个前三体构型的三体船模型[18]。该船模的主尺度如表1所示,横剖面示意图如图1所示。

      表 1  船模主尺度参数

      Table 1.  Main dimensions of ship model

      参数主船体侧体
      船长/m4.01.0
      设计水线宽/m0.358 40.085
      设计吃水/m0.170.10
      排水量/kg129.074.45
      湿表面积/m21.8990.201
      方形系数0.5300.524
      棱形系数0.7690.746

      图  1  三体船模型横剖面示意图

      Figure 1.  Schematic diagram of transverse section of trimaran model

      三体船的构型布局如图2所示。图中:$a$为侧体中心位置与主体中心位置的纵向距离,侧体中心位于主体舯前为正;$b$为侧体中心位置与主体中心位置的横向距离,$b$值始终为正;L为主体长度;Ls为侧体长度。本次计算的三体船模型对应的构型参数a=1.0 m,b=0.7 m。

      图  2  三体船布局示意图

      Figure 2.  Schematic diagram of trimaran layout

    • 借助CFD商业软件STAR-CCM+完成数值模拟,采用有限体积法进行空间离散,流体方程则采用非定常求解器求解。所选$k - \varepsilon $湍流模型选择适用于复杂几何外部的流动问题求解。自由液面的捕捉采用流体体积(volume of fluid,VOF)法。均匀来流采用静水VOF法处理。三体船的六自由度运动模拟利用DFBI (dynamic fluid body interaction) 求解完成,自由运动模型释放了纵摇和垂荡这2个方向的自由度[16-17]

      本文的数值计算采用重叠网格技术,将计算域分成背景区域和重叠网格区域两部分。使用2套网格进行网格划分,计算总网格数240万。其中,在主体和2个侧体周围以及自由表面处均进行网格加密。

      本次计算的流体域尺度以主体长度L为基准选取。流体域长4.5L,宽3L,高2L,其中水域深度为L,空气域高度为L。计算流体域的边界条件由5个速度入口和1个压力出口组成。主体的船艏距离速度入口的距离为L

      数值计算模型的具体边界条件设置如图3所示。

      图  3  边界条件设置

      Figure 3.  The boundary conditions

    • 本次验证的模拟工况为:Fr = 0.546(航速V = 3.416 m/s)。分别求时间步长Δt = 0.01,0.02和0.04 s这3种条件下,三体船在静水中的总阻力时历曲线。取40~70 s时间段内的结果曲线进行比较,如图4(图中,Rt为总阻力)所示。由图中曲线可见,3种时间步长的总阻力计算平均值与试验值的偏差约为6%,Δt=0.02 s时的总阻力时历曲线与Δt = 0.01 s时的非常接近。

      图  4  不同时间步长下的总阻力时历曲线

      Figure 4.  Time history curves of total resistance for different time steps

      取40 ~70 s时间段内的总阻力计算平均值作为偏差衡量标准,以Δt = 0.01 s的计算结果为基准,比较其他时间步长计算结果与其之间的偏差,如表2所示。从表中可看出,Δt = 0.02 s的总阻力计算平均值与Δt = 0.01 s的计算平均值仅相差−0.1%,满足时间步无关性要求。本文数值模拟计算采用的时间步长为Δt = 0.02 s。

      表 2  不同时间步长下的计算结果

      Table 2.  The calculation results for different time steps

      时间步长Δt/s总阻力Rt/N偏差/%
      0.0192.090
      0.0291.98−0.1
      0.0493.14 1.1
    • 通过改变自动网格模块下基础尺度的值来调整网格数量。研究了3种不同网格数下的CFD计算模型,3种计算模型的网格数分别为150万、240万和392万。模拟工况为Fr =0.546(V =3.416 m/s),Δt = 0.02 s。

      图5所示为3种网格模型的总阻力时历曲线比较。由图可见,网格数为240万和392万的计算结果曲线非常接近。

      图  5  不同网格模型的总阻力时历曲线

      Figure 5.  Time history curves of total resistance for different grid sizes

      取40~70 s时间段内的总阻力计算平均值作为偏差衡量标准,以网格数392万的计算结果为基准,比较其他网格数计算结果与其之间的偏差,结果如表3所示。由表3中可以看出,网格数为240万和392万的计算结果平均值仅相差−0.2%,这表明该网格数对计算结果几乎无影响,满足网格收敛性要求。本文的数值模拟采用网格数为240万的计算模型。

      表 3  模型网格信息和计算结果

      Table 3.  The grids and calculation results

      网格数总阻力Rt/N偏差/%
      150万91.25−1.0
      240万91.98−0.2
      392万92.190

      根据以上验证结果,选择时间步长Δt = 0.02 s、网格数量为240万的数值模型进行计算。数值计算范围为:Fr = 0.1~1.0(${{F}}{{{r}}_\nabla } = 0.28\sim 2.80$),共选取21个速度点。

    • 计算前三体船约束模型的静水总阻力,并将其与自由状态下的静水阻力试验值[18]进行比较,如图6所示。

      图  6  约束模型总阻力与试验值的对比

      Figure 6.  Comparison between total resistance and experimental values of constrained model

      由图可见,Fr < 0.3时,约束模型的计算总阻力与自由状态下的试验值非常接近;而在Fr > 0.3的中、高速段,两者相差较大,其中在Fr = 0.5时试验值比约束模型计算值高出了34.4%。

    • 在自由状态下的静水阻力试验中,测量该前三体船模艏、艉测量点的升沉值[16]。模型试验中艏、艉测量点的位置如图7所示。其中,升沉值沿垂直方向上升为正,下降为负;纵倾值指船体纵向绕质心旋转的角度,艏倾为正,艉倾为负。

      图  7  艏、艉升沉测量点位置示意图

      Figure 7.  The position diagram of sinkage measurement points at bow and stern

      对前三体船自由模型的升沉和纵倾进行数值计算,监测三体船的升沉和纵倾随Fr的变化,并与自由状态下的模型试验值进行比较,结果分别如图8图9所示。

      图  8  三体船模型艏、艉测量点的升沉曲线

      Figure 8.  The sinkage curves of measurement points at bow and stern of trimaran model

      图  9  三体船模型的纵倾曲线

      Figure 9.  The trim curves of trimaran model

      图8可见,自由模型艏、艉测量点的升沉计算值与试验数据吻合良好。船体艏部测量点在Fr > 0.4时出现了明显的抬升,而艉部测量点则是从Fr = 0.3开始出现明显的下沉。随着航速的不断增加,艏、艉测量点的升沉值也随之增大,至Fr > 0.6时艏、艉测量点的升沉值均趋于一个比较稳定的值,随航速的变化开始变缓。

      图9可见,自由模型的纵倾计算值与试验数据吻合较好。在整个航速区间内,三体船的纵倾角均为负值,即船体在不同航速下航行时,始终处于艉倾状态。当Fr < 0.3时,三体船的纵倾并不明显;当0.3 < Fr < 0.6时,纵倾角随着航速的增加快速增大;在Fr = 0.6时,纵倾角达到最大值,约1.44°;当0.6 < Fr < 0.8时,纵倾角几乎保持不变,其大小稳定在1.4°左右;而当Fr > 0.8时,纵倾角则随航速的增加缓慢减小,在Fr = 1.0时,其值约为1.25°。

      由以上数值计算和模型试验结果可以看出,在Fr > 0.3的中、低速范围内,三体船的航行姿态变化并不明显,而在Fr > 0.3的中、高速段,三体船出现了非常显著的航行姿态变化。

    • 将数值计算得到的自由模型静水总阻力与自由状态下的模型试验总阻力进行对比,如图10所示。由图可见,当Fr < 0.3时,自由模型的计算总阻力与试验值非常接近,这与约束模型的计算结果一致,说明当Fr < 0.3时,航行姿态变化较小,对总阻力的影响甚微,是否考虑航行姿态的影响几乎不影响总阻力计算结果。而当Fr > 0.3时,约束模型阻力计算值与试验值的偏差幅度增大,这说明考虑了航行姿态变化的自由模型计算值精度显著提高。当Fr = 0.5时,约束模型的总阻力计算值与试验值相差34.4%,而自由模型阻力计算值与试验值的偏差值则降低了7.2%,即计算中考虑航行姿态变化后计算结果精度显著提高。

      图  10  自由模型总阻力与试验值的对比

      Figure 10.  Comparison of total resistance and experimental values of free model

    • 自由模型的摩擦阻力Rf和压差阻力Rp随航速变化的曲线如图11所示。

      图  11  自由模型的摩擦阻力和压差阻力曲线

      Figure 11.  The frictional resistance and pressure resistance curves of free model

      图11可见,在Fr > 0.65的高速段,压差阻力并没有随着航速的增大而迅速增加,反而呈现出非常平缓的增长趋势。这主要与船体自身的型线特征有关,说明该船型在高航速下具备优良的兴波阻力性能。

      为了定量分析摩擦阻力和压差阻力在自由模型总阻力中所占的比例,绘制了两种阻力成分占总阻力的比例曲线,如图12所示。由图可见,当0.3 < Fr < 0.65时,压差阻力明显大于摩擦阻力;而在Fr > 0.65的高速段,摩擦阻力成为总阻力的主要成分,并且其比例随航速的增加而迅速增大,其中在Fr = 1.0时,摩擦阻力是压差阻力的2倍。

      图  12  自由模型摩擦阻力和压差阻力百分比曲线

      Figure 12.  The percentage curves of frictional resistance and pressure resistance of free model

    • 为便于比较,定义了以下参数:δRtδRfδRpεtεfεp。其中:δRt=Rt自由Rt约束δRf=Rf自由Rf约束δRp=Rp自由Rp约束δRtδRfδRp分别为自由模型的总阻力、摩擦阻力和压差阻力变化量,Rt自由Rt约束Rf自由Rf约束Rp自由Rp约束分别为自由模型和约束模型的总阻力、摩擦阻力和压差阻力;εt =δRt/Rt约束εf =δRf/Rf约束εp =δRp/Rp约束εtεfεp分别为自由模型的总阻力变化率、摩擦阻力变化率和压差阻力变化率。

    • 自由模型与约束模型的数值计算总阻力曲线如图13所示。由图可见,2种模型的总阻力曲线在0.3 < Fr < 0.88时差别较明显,且自由模型的总阻力明显高于约束模型;在Fr > 0.88的甚高速段,自由模型的总阻力略低于约束模型。

      图  13  自由模型与约束模型的总阻力曲线

      Figure 13.  The total resistance curves of free model and constrained model

      为定量分析不同航速下航行姿态对总阻力的影响幅度,绘制了自由模型的总阻力变化率εt曲线,如图14所示。

      图  14  自由模型的总阻力变化率曲线

      Figure 14.  The curve of total resistance change rate of free model

      图14可见,在整个航速范围内,自由模型的总阻力变化率εt变化显著。其中,0.3 < Fr < 0.85时的总阻力变化率εt明显高于其他航速段,并在Fr = 0.5,0.7时出现了2个峰值,总阻力变化率分别为24.8%和14.4%;在Fr>0.88的高速段,自由模型的总阻力反而略低于约束模型,总阻力变化率约为−2%。

      由此可见,在半滑行高速段,航行过程中的航行姿态变化对三体船总阻力的影响非常显著。

    • 自由模型和约束模型的摩擦阻力比较曲线如图15所示。由图可见,当Fr < 0.55时,两者相差较很小;而当Fr > 0.8时,两者的差距较为显著。

      图  15  自由模型和约束模型的摩擦阻力曲线

      Figure 15.  The frictional resistance curves of free model and constrained model

      自由模型的摩擦阻力变化率εf曲线如图16所示。由图可见,当0.3 < Fr < 0.5时,自由模型的摩擦阻力变化率约为2%~3%;在Fr > 0.55的高速段,自由模型的摩擦阻力变化率为负值,其中在Fr = 0.9时达到最低值,约为−8.4%。

      图  16  自由模型的摩擦阻力变化率曲线

      Figure 16.  The curve of frictional resistance change rate of free model

      与约束模型相比,自由模型的摩擦阻力在Fr < 0.55时是增加的,但在Fr > 0.55时不增反降。

      为分析主体和双侧体对三体船摩擦阻力改变量的影响效果,绘制了主体、双侧体以及三体船的摩擦阻力改变量曲线,如图17所示。

      图  17  自由模型的摩擦阻力改变量曲线

      Figure 17.  The curves of frictional resistance change of free model

      图17可见,当Fr > 0.4时,双侧体的摩擦阻力随Fr的增大而不断减小。这主要是由于自Fr = 0.4之后,三体船出现了非常显著的升沉和纵倾现象,双侧体的部分船体抬出水面,从而导致湿表面积不断减小,进而引起摩擦阻力的减小。主体的摩擦阻力增加量在Fr < 0.85时始终为正,而在Fr > 0.85时为负。这是因为主体艏部尖瘦而艉部宽肥,纵倾角较小时的艉倾现象能够带来湿表面积的增大,从而导致摩擦阻力增加;而当速度更高时,方艉船型会使得船体艉部小幅抬升,艉部横剖面不再接触水面,并出现“虚长度”,此时湿表面积减小,从而导致摩擦阻力减小。图18所示的图片中也显示出了这一现象。

      图  18  高航速时的试验图片和计算图片(Fr=0.9)

      Figure 18.  The experimental picture and computational contours at high speed(Fr=0.9)

    • 自由模型和约束模型的压差阻力曲线比较如图19所示。由图可见,自由模型和约束模型的压差阻力曲线在Fr > 0.3时差别非常显著。相应的压差阻力变化率εp曲线如图20所示。

      图  19  自由模型和约束模型的压差阻力曲线

      Figure 19.  The pressure resistance curves of free model and constrained model

      图  20  自由模型的压差阻力变化率曲线

      Figure 20.  The curve of pressure resistance change rate of free model

      图20可见,在整个航速区间内,自由模型的压差阻力变化率均为正值。当Fr > 0.3时,压差阻力变化率均大于10%,并出现了2个波峰,分别对应Fr = 0.5和0.7,变化率分别为48.7%和40.7%。对比自由模型的纵倾曲线可以看到,在整个航速范围内,纵倾值始终为负,即船舶始终处于艉倾状态航行。由于该船型在船舯之后的型线偏肥大,因此这样的航行姿态会导致更大的兴波阻力和黏压阻力。

    • 为了分析自由模型状态下不同阻力成分对总阻力变化的贡献,绘制了如图21所示的曲线,纵坐标为摩擦阻力和压差阻力的变化量。

      图  21  自由模型摩擦阻力和压差阻力的变化量

      Figure 21.  The change of friction resistance and pressure resistance of free model

      图21可见,在整个航速区间内,压差阻力的变化量始终为正;摩擦阻力的变化量在Fr < 0.55时为正,在Fr > 0.55时为负;压差阻力的变化量曲线明显高于摩擦阻力。

      结合总阻力变化率分析(图14)可知,在Fr < 0.88范围内,总阻力的变化量为正值。当Fr=0.5时,压差阻力变化量δRp =15.7 N,摩擦阻力变化量δRf =0.5 N,总阻力增量的97%来自于压差阻力。当Fr =0.7时,压差阻力变化量δRp =16.76 N,而摩擦阻力变化量δRf = −1.17 N。因航行姿态变化引起的总阻力增量90%以上来自于压差阻力的变化。

    • 本文选用标准$k - \varepsilon $湍流模型,采用静水VOF法,使用重叠网格技术,针对某前三体船型的自由模型和约束模型进行了静水阻力计算,分析了该三体船型的阻力特性,并重点分析了该三体船航行姿态对阻力性能的影响。主要得到以下结论:

      1) 当弗劳德数(Fr < 0.3)较低时,航行姿态对阻力的影响可以忽略不计,这也是排水型船进行阻力预报时不考虑航行姿态影响的原因。一旦进入半滑行阶段,将发生显著的升沉和纵倾现象(纵倾可高达${1.44^ \circ }$),进而阻力也会发生显著变化。

      2) 在半滑行状态,航行姿态对摩擦阻力和压差阻力这2个阻力成分的影响不同,其对压差阻力的影响可高达近50%,而对摩擦阻力的影响则不到10%。在半滑行状态,由航行姿态变化引起的阻力增加90%以上来自压差阻力的增加。

参考文献 (18)

目录

    /

    返回文章
    返回