留言板

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

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

T型翼主动控制策略的数值研究

姜宜辰 白钧宇 孙一方 李永坤 张桂勇

姜宜辰, 白钧宇, 孙一方, 等. T型翼主动控制策略的数值研究[J]. 中国舰船研究, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757
引用本文: 姜宜辰, 白钧宇, 孙一方, 等. T型翼主动控制策略的数值研究[J]. 中国舰船研究, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757
JIANG Y C, BAI J Y, SUN Y F, et al. The numerical study on active control strategies of T-foil[J]. Chinese Journal of Ship Research, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757
Citation: JIANG Y C, BAI J Y, SUN Y F, et al. The numerical study on active control strategies of T-foil[J]. Chinese Journal of Ship Research, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757

T型翼主动控制策略的数值研究

doi: 10.19693/j.issn.1673-3185.01757
基金项目: 青岛海洋科学与技术国家实验室开放基金资助项目(QNLM2016ORP0402);国家自然科学基金资助项目(51639003,51279030);中央高校基本科研业务费专项资金资助项目(DUT2017TB05)
详细信息
    作者简介:

    姜宜辰,男,1984年生,副教授,硕士生导师。研究方向:高性能船舶,船舶智能控制,海上新能源。E-mail:ycjiang@dlut.edu.cn

    白钧宇,男,1996年生,硕士生。研究方向:高速船T型翼主动运动控制技术。E-mail:junyubai7788@mail.dlut.edu.cn

    孙一方,男,1987年生,博士生。研究方向:船舶与海洋结构物设计制造。E-mail:sunyf4787@foxmail.com

    张桂勇,男,1978年生,博士,教授,博士生导师。研究方向:流固耦合,极地船舶,无网格方法。E-mail:gyzhang@dlut.edu.cn

    通讯作者:

    张桂勇

  • 中图分类号: U661.3

The numerical study on active control strategies of T-foil

  • 摘要:   目的  主动控制T型翼对船舶有着较好的纵向运动减摇效果,但最佳的T型翼主动控制策略仍未有定论,因此需对T型翼主动控制策略进行研究计算与分析。  方法  首先,采用细长体理论对Wigley III船型在不同波长规则波中的纵向运动响应进行数值计算;然后,在相同航速与波浪工况下,对加装了主、被动T型翼的船舶纵向运动进行数值预报,分析不同主动控制策略的减摇效果。  结果  计算结果显示:T型翼被动控制在共振区具有较好的减摇效果,长波区和短波区的减摇效果相对有限;纵摇角速度信号控制对纵摇运动和艏部运动加速度的减摇效果非常显著,但对于长波区的垂荡运动会有减益效果;船艏部运动速度信号对纵摇运动和艏部运动加速度的减摇效果较纵摇角速度信号均有10%~20%的提升,并且对垂荡运动的减摇效果也有一定提升。  结论  对于T型翼的主动控制,艏部运动速度信号控制综合考虑了纵摇与垂荡这2种运动,可以实现更好的减摇效果。
  • 图  1  T型翼攻角分析图

    Figure  1.  Analysis of T-foil angle of attack

    图  2  Wigley III船型

    Figure  2.  Hull form of Wigley III

    图  3  T型翼示意图

    Figure  3.  Schematic diagram of the T-foil

    图  4  数值计算结果与试验值的对比

    Figure  4.  Comparison of numerical results and the experimental results

    图  5  T型翼摆角与纵摇角速度的对应关系(Fr=0.3,$\lambda /L = 1.75$

    Figure  5.  The relationship between T-foil angle and angle velocity of pitch(Fr=0.3,$\lambda /L = 1.75$

    图  6  T型翼摆角与艏部运动速度的对应关系(Fr=0.3,$\lambda /L $=1.75)

    Figure  6.  The relationship between T-foil angle and bow velocity (Fr=0.3,$\lambda /L$=1.75)

    图  7  Fr=0.3时3种控制方式的垂荡运动响应

    Figure  7.  Heave motion response with three types of control method at Fr=0.3

    图  8  Fr=0.5时3种控制方式的垂荡运动响应

    Figure  8.  Heave motion response with three types of control method at Fr=0.5

    图  9  Fr=0.3时3种控制方式的纵摇运动响应

    Figure  9.  Pitch motion response with three types of control method at Fr=0.3

    图  10  Fr=0.5时3种控制方式的纵摇运动响应

    Figure  10.  Pitch motion response with three types of control method at Fr=0.5

    图  11  Fr=0.3时3种控制方式的艏部运动加速度运动响应

    Figure  11.  Bow acceleration motion response with three types of control method at Fr=0.3

    图  12  Fr=0.5时3种控制方式的艏部运动加速度运动响应

    Figure  12.  Bow acceleration motion response with three types of control method at Fr=0.5

    图  13  2种T型翼主动控制策略在2种工况下对船产生的纵摇力矩

    Figure  13.  The pitching moment of two type T-foil active control strategies for the ship in two cases

    图  14  采用艏部运动速度控制方法的升力矩与纵摇角速度和垂荡速度的变化关系

    Figure  14.  The relationship between lifting moment and pitching angle velocity and vertical velocity controlled by the heading velocity control method

    图  15  采用纵摇角速度控制方法的升力矩与纵摇角速度和垂荡速度的变化关系

    Figure  15.  The relationship between lifting moment and angular velocity in pitch and heave velocity controlled by the pitch angle velocity control method

    表  1  Wigley III船型参数

    Table  1.   Parameters of Wigley hull III

    参数数值
    垂线间长/cm300
    船宽/cm30
    设计吃水/cm18.75
    排水量/g78 000
    中横剖面系数0.666 7
    重心高度/cm0.17
    下载: 导出CSV

    表  2  T型翼参数

    Table  2.   Parameters of the T-foil

    参数数值
    展长/cm9.0
    弦长/cm6.0
    垂直翼长度/cm5.0
    下载: 导出CSV

    表  4  Fr=0.3时的垂荡运动减摇效果

    Table  4.   Damping effect of heave motion at Fr=0.3

    $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
    被动
    控制
    纵摇角
    速度信号
    艏部运动
    速度信号
    0.759.7823.0096.2087.15
    1.256.787.1911.1912.55
    1.755.37−2.60−11.422.28
    2.254.54−2.17−11.361.35
    2.753.97−1.67−10.931.90
    下载: 导出CSV

    表  5  Fr=0.5时的垂荡运动减摇效果

    Table  5.   Damping effect of heave motion at Fr=0.5

    $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
    被动
    控制
    纵摇角
    速度信号
    艏部运动
    速度信号
    0.7512.81−5.7219.1140.46
    1.258.6055.1049.0273.72
    1.756.6736.4433.8526.88
    2.255.55−1.2911.3316.50
    2.754.80−3.514.6515.85
    下载: 导出CSV

    表  6  Fr=0.3时的纵摇运动减摇效果

    Table  6.   Damping effect of pitch motion at Fr=0.3

    $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
    被动
    控制
    纵摇角
    速度信号
    艏部运动
    速度信号
    0.759.783.7037.9929.44
    1.256.7827.9140.4244.98
    1.755.379.3724.0839.82
    2.254.547.2422.7147.47
    2.753.976.8827.8256.34
    下载: 导出CSV

    表  7  Fr=0.5时的纵摇运动减摇效果

    Table  7.   Damping effect of pitch motion at Fr=0.5

    $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
    被动
    控制
    纵摇角
    速度信号
    艏部运动
    速度信号
    0.7512.81−5.0568.7249.70
    1.258.6033.3879.7831.16
    1.756.6765.0482.7887.33
    2.255.5518.8662.1583.19
    2.754.8011.6662.1286.05
    下载: 导出CSV

    表  8  Fr=0.3时的艏部运动加速度减摇效果

    Table  8.   Damping effect of bow acceleration motion at Fr=0.3

    $\lambda /L$${\omega _ {\rm {e}} }$/(rad·s−1)减摇效果/%
    被动
    控制
    纵摇角
    速度信号
    艏部运动
    速度信号
    0.759.789.7157.7647.37
    1.256.7824.2334.1141.36
    1.755.373.914.4327.80
    2.254.541.38−5.0426.79
    2.753.970.85−9.4027.45
    下载: 导出CSV

    表  9  Fr=0.5时的艏部运动加速度减摇效果

    Table  9.   Damping effect of bow motion acceleration at Fr=0.5

    $\lambda /L$${\omega _ {\rm {e}} }$/(rad·s−1)减摇效果/%
    被动
    控制
    纵摇角
    速度信号
    艏部运动
    速度信号
    0.7512.813.9186.7574.00
    1.258.6047.5565.0876.32
    1.756.6757.6472.6879.27
    2.255.554.3629.0858.50
    2.754.80−0.9518.4058.45
    下载: 导出CSV
  • [1] DAVIS M R, WATSON N L, HOLLOWAY D S. Wave response of an 86 m high speed catamaran with active T-foils and stern tabs[J]. Transactions of the Royal Institution of Naval Architects, Part A: International Journal of Maritime Engineering, 2003, 145(A2): 15–34.
    [2] FANG C C, CHAN H S. An investigation on the vertical motion sickness characteristics of a high-speed catamaran ferry[J]. Ocean Engineering, 2007, 34(14/15): 1909–1917.
    [3] HAYWOOD A J, DUNCAN A J, KLAKA K P, et al. The development of a ride control system for fast ferries[J]. Control Engineering Practice, 1995, 3(5): 695–702. doi: 10.1016/0967-0661(95)00046-W
    [4] 李鹏, 姜治芳, 黄卫刚. 被动翼复合船型航行性能模型试验研究[J]. 舰船科学技术, 2016, 38(10): 25–28, 33.

    LI P, JIANG Z F, HUANG W G. Research on the experiment of navigation performance for hybrid monohull combined with passive foil[J]. Ship Science and Technology, 2016, 38(10): 25–28, 33 (in Chinese).
    [5] 杨强, 林壮, 郭志群. 带T型翼三体船纵向运动的预报研究[J]. 华中科技大学学报: 自然科学版, 2014, 42(4): 86–90.

    YANG Q, LIN Z, GUO Z Q. Research on motion simulation of a T-foil fitted trimaran[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2014, 42(4): 86–90 (in Chinese).
    [6] 杨帅, 谢伟, 朱志洁, 等. 采用自控T型翼和尾板改善穿浪双体船波浪中运动研究[J]. 中国水运, 2015, 15(9): 1–3.

    YANG S, XIE W, ZHU Z J, et al. Research on the wave motion of catamaran ship using automatic control T-foil and stern tab[J]. China Water Transport, 2015, 15(9): 1–3 (in Chinese).
    [7] 刘冰. 高速双体船纵向运动及其控制研究[D]. 哈尔滨: 哈尔滨工程大学, 2012.

    LIU B. Study of longitudinal motion and control of high speed catamaran[D]. Harbin: Harbin Engineering University, 2012 (in Chinese)
    [8] 孙一方, 姜宜辰, 宗智. T型翼控制船舶在波浪上纵向运动的数值研究[J]. 中国舰船研究, 2018, 13(2): 16–26. doi: 10.3969/j.issn.1673-3185.2018.02.003

    SUN Y F, JIANG Y C, ZONG Z. Numerical research on control of longitudinal motion of ship in waves by T-foil[J]. Chinese Journal of Ship Research, 2018, 13(2): 16–26 (in Chinese). doi: 10.3969/j.issn.1673-3185.2018.02.003
    [9] YEUNG R W, JIANG Y C. Bilge keel influence on the free decay of roll motion of a three-dimensional hull[C]//Proceedings of the 33rd International Conference on Ocean, Offshore and Arctic Engineering. San Francisco: ASME, 2014.
    [10] YEUNG R W, KIM S H. Radiation forces on ships with forward speed[C]//Proceedings of the 3rd International Conference on Numerical Ship Hydrodynamics. Paris: Palais des Congres, 1981: 499–510.
    [11] JOURNÉE J M J. Experiments and calculations on 4 Wigley hull forms in head waves[R]. Delft, Netherlands: Delft University of Technology, 1992.
  • [1] 张庆伟, 俞翔, 杨理华.  反馈式多线谱主动隔振控制算法研究 . 中国舰船研究, 2021, 16(): 1-8. doi: 10.19693/j.issn.1673-3185.02107
    [2] 胡浩俊, 张军, 刘志林, 李国胜.  考虑附体输入约束的高速多体船预测控制减纵摇方法 . 中国舰船研究, 2021, (): 1-8. doi: 10.19693/j.issn.1673-3185.02050
    [3] 王文晋, 王冠学, 李奔, 徐国华.  大尺度高速水下无人艇控制系统设计与试验验证 . 中国舰船研究, 2020, 15(2): 95-103. doi: 10.19693/j.issn.1673-3185.01759
    [4] 陆超, 崔敬玉, 孟凡华, 余克涛.  高速船型融合特型球鼻艏的多方案优化设计 . 中国舰船研究, 2020, 15(3): 54-60. doi: 10.19693/j.issn.1673-3185.01556
    [5] 姚金艺, 曾庆军, 赵强, 朱志宇, 戴文文.  全驱动型AUV三维路径跟踪控制系统设计及分析 . 中国舰船研究, 2019, 14(6): 22-29. doi: 10.19693/j.issn.1673-3185.01520
    [6] 张光义, 曾庆军, 戴晓强, 朱春磊, 凌宏杰.  水下安全检测与作业型机器人控制系统 . 中国舰船研究, 2018, 13(6): 113-119. doi: 10.19693/j.issn.1673-3185.01167
    [7] 孙一方, 姜宜辰, 宗智.  T型翼控制船舶在波浪上纵向运动的数值研究 . 中国舰船研究, 2018, 13(2): 16-26. doi: 10.3969/j.issn.1673-3185.2018.02.003
    [8] 张玮康, 王冠学, 徐国华, 刘畅, 申雄.  腹部作业型水下机器人控制系统研制 . 中国舰船研究, 2017, 12(2): 124-132. doi: 10.3969/j.issn.1673-3185.2017.02.016
    [9] 岳林.  船机电系统自主控制平台及资源分配策略 . 中国舰船研究, 2016, 11(4): 102-110. doi: 10.3969/j.issn.1673-3185.2016.04.015
    [10] 赵淑琴, 张永生.  船用蒸汽发生器给水系统的容错控制 . 中国舰船研究, 2016, 11(2): 117-120,138. doi: 10.3969/j.issn.1673-3185.2016.02.016
    [11] 张永生, 赵淑琴.  船用蒸汽发生器水位的模糊自适应PID控制 . 中国舰船研究, 2013, 8(3): 106-109. doi: 10.3969/j.issn.1673-3185.2013.03.020
    [12] 周芃, 李为, 陈韬颖.  多块结构化网格在高速水面舰船阻力计算中的应用 . 中国舰船研究, 2013, 8(2): 34-40. doi: 10.3969/j.issn.1673-3185.2013.02.007
    [13] 刘辉, 戈亮, 严军.  船用蝶阀液压控制系统的数值仿真 . 中国舰船研究, 2012, 7(2): 97-101107. doi: 10.3969/j.issn.1673-3185.2012.02.018
    [14] 张永生, 马运义.  基于PID的船用汽轮机功频控制系统 . 中国舰船研究, 2011, 6(4): 89-91. doi: 10.3969/j.issn.1673-3185.2011.04.019
    [15] 宋媚婷, 洪明, 崔洪宇.  基于舰船振动信号激励下压电智能梁振动主动控制 . 中国舰船研究, 2011, 6(6): 28-33. doi: 10.3969/j.issn.1673-3185.2011.06.006
    [16] 邱斌, 吴卫国, 刘恺.  高速船全频段舱室噪声仿真预报 . 中国舰船研究, 2011, 6(6): 49-53. doi: 10.3969/j.issn.1673-3185.2011.06.010
    [17] 董海运, 肖航.  船用锅炉汽包水位模糊控制研究 . 中国舰船研究, 2009, 4(1): 73-76. doi: 10.3969/j.issn.1673-3185.2009.01.016
    [18] 王学杰, 贺国, 陈连树.  基于二叉树的某型船柴油机控制系统故障诊断专家系统 . 中国舰船研究, 2009, 4(06): 34-37. doi: 10.3969/j.issn.1673-3185.2009.06.008
    [19] 郦茜, 吴卫国.  基于AutoSEA的高速船静噪声预报与控制 . 中国舰船研究, 2008, 3(1): 28-30. doi: 10.3969/j.issn.1673-3185.2008.01.007
    [20] 金全洲, 李天匀, 赵耀, 严谨.  基于主动力方式的圆柱壳结构振动功率流控制分析 . 中国舰船研究, 2006, 1(1): 17-20,49. doi: 10.3969/j.issn.1673-3185.2006.01.004
  • 加载中
图(15) / 表 (8)
计量
  • 文章访问数:  18
  • HTML全文浏览量:  5
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-09-04
  • 修回日期:  2019-12-22
  • 网络出版日期:  2020-08-05
  • 刊出日期:  2021-04-01

T型翼主动控制策略的数值研究

doi: 10.19693/j.issn.1673-3185.01757
    基金项目:  青岛海洋科学与技术国家实验室开放基金资助项目(QNLM2016ORP0402);国家自然科学基金资助项目(51639003,51279030);中央高校基本科研业务费专项资金资助项目(DUT2017TB05)
    作者简介:

    姜宜辰,男,1984年生,副教授,硕士生导师。研究方向:高性能船舶,船舶智能控制,海上新能源。E-mail:ycjiang@dlut.edu.cn

    白钧宇,男,1996年生,硕士生。研究方向:高速船T型翼主动运动控制技术。E-mail:junyubai7788@mail.dlut.edu.cn

    孙一方,男,1987年生,博士生。研究方向:船舶与海洋结构物设计制造。E-mail:sunyf4787@foxmail.com

    张桂勇,男,1978年生,博士,教授,博士生导师。研究方向:流固耦合,极地船舶,无网格方法。E-mail:gyzhang@dlut.edu.cn

    通讯作者: 张桂勇
  • 中图分类号: U661.3

摘要:   目的  主动控制T型翼对船舶有着较好的纵向运动减摇效果,但最佳的T型翼主动控制策略仍未有定论,因此需对T型翼主动控制策略进行研究计算与分析。  方法  首先,采用细长体理论对Wigley III船型在不同波长规则波中的纵向运动响应进行数值计算;然后,在相同航速与波浪工况下,对加装了主、被动T型翼的船舶纵向运动进行数值预报,分析不同主动控制策略的减摇效果。  结果  计算结果显示:T型翼被动控制在共振区具有较好的减摇效果,长波区和短波区的减摇效果相对有限;纵摇角速度信号控制对纵摇运动和艏部运动加速度的减摇效果非常显著,但对于长波区的垂荡运动会有减益效果;船艏部运动速度信号对纵摇运动和艏部运动加速度的减摇效果较纵摇角速度信号均有10%~20%的提升,并且对垂荡运动的减摇效果也有一定提升。  结论  对于T型翼的主动控制,艏部运动速度信号控制综合考虑了纵摇与垂荡这2种运动,可以实现更好的减摇效果。

English Abstract

姜宜辰, 白钧宇, 孙一方, 等. T型翼主动控制策略的数值研究[J]. 中国舰船研究, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757
引用本文: 姜宜辰, 白钧宇, 孙一方, 等. T型翼主动控制策略的数值研究[J]. 中国舰船研究, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757
JIANG Y C, BAI J Y, SUN Y F, et al. The numerical study on active control strategies of T-foil[J]. Chinese Journal of Ship Research, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757
Citation: JIANG Y C, BAI J Y, SUN Y F, et al. The numerical study on active control strategies of T-foil[J]. Chinese Journal of Ship Research, 2021, 16(2): 15–23 doi: 10.19693/j.issn.1673-3185.01757
    • 高速船由于其细长的船体以及高弗劳德数,相比在较低弗劳德数下航行的传统单体船,有更大的垂荡和纵摇运动响应[1]。而T型翼作为一种重要的船舶纵向运动减摇附体,具有能有效提高船舶的纵向运动阻尼和附加阻力小的优点,在高速船上得到了广泛应用[2]。因为可以随航态实时改变水平翼攻角的T型翼具有很好的减摇效果[3],所以T型翼的主动控制方法一直是近年来国内外研究的热点。

      李鹏等[4]采用模型试验方法,针对一对艏部水翼和T型翼对阻力及垂向运动的综合影响进行了分析,结果表明,被动翼复合船型在风浪中航行时,其快速性与原船型相当,而耐波性则有大幅提升。杨强等[5]对某三体船及其适配T型翼在规则波与不规则波下的运动响应进行了计算,结果显示,在高速运动下最高可使纵摇响应幅值降低50.33%。杨帅等[6]通过Simulink建立控制模型,并采用PID控制对控制效果进行了仿真,结果显示,采用主动控制系统的穿浪双体船的运动响应与被动式相比降低了30%。刘冰[7]构建了T型翼和扰流板的升力模型并对某高速双体船进行了水动力计算,计算了采用纵摇角、纵摇角速度和纵摇角加速度联合控制方式下的减摇效果,结果显示,采用主动控制策略后垂荡与纵摇的最大减摇幅值可达50%。

      孙一方等[8]采用纵摇信号对T型翼的减摇效果进行分析讨论,采用力矩控制的方式得到了T型翼攻角随船体运动变化的控制方程,通过裸船在波浪下的运动响应计算T型翼的线性控制参数C,再应用控制参数C对T型翼的摆角进行主动控制,得出在纵摇角、纵摇角速度和纵摇角加速度中,采用纵摇角速度信号控制的T型翼减摇效果最佳。其研究表明,该主动控制方式对纵摇运动的减摇效果显著,但对垂荡与艏部加速度的减摇效果有限。

      为此,本文将以Wigley船型为研究对象,基于细长体理论,分别计算2种航速(Fr=0.3,0.5)下及5种波长/船长比($\lambda /{L}$=0.75,1.25,1.75,2.25,2.75)规则波中裸船与被动式T型翼船运动的垂荡、纵摇和艏部加速度运动响应;然后,在此基础上讨论不同主动控制方式的减摇效果;最后,将线性控制计算结果与无T型翼裸船及被动式T型翼船的结果进行对比,得出不同主动控制方式下T型翼减摇效果的影响并分析其优劣性。

    • 由于高速船具有细长体的特征,船舶垂向与横向运动速度的量级远高于船舶沿船长方向速度的量级,所以可以通过边界元法将三维流体问题降阶为多个二维流体问题,以此对船舶的运动进行近似求解。在这种近似求解方法中,每个固定的空间位置都有相应的切片,以计算船舶在穿过某个切片时的运动[9-10]。将船舶穿过的一系列二维切片的运动在空间上组合起来,即可得到船舶在航行过程中各个方向的运动响应。本文计算依据的流场速度势公式与边界条件、三维船体的运动方程和船体与流体运动动态耦合的求解方程的推导过程参见文献[8]。

    • 为求解T型翼产生的力或者力矩,文献[4]给出了求解方法。船体以给定的弗劳德数迎浪航行,T型翼的水平翼的投影面积设为A,升力系数为${C_ {\rm {L}} }$,T型翼产生的垂向升力${F_ {\rm {T}}}$和升力矩${M_ {\rm {T}}}$的表达式如下:

      $${F_ {\rm {T}}} = \frac{1}{2}\rho {U^2}A\frac{{{\rm{d}}{C_ {\rm {L}}}}}{{{\rm{d}}\alpha }}\alpha $$ (1)
      $${M_ {\rm {T}}} = - \frac{1}{2}{l_ {\rm {F}}}\rho {U^2}A\frac{{{\rm{d}}{C_ {\rm {L}}}}}{{{\rm{d}}\alpha }}\alpha $$ (2)

      式中:${l_ {\rm {F}}}$为水平翼尾缘距离船舯的距离;$\rho $为流体密度;$U$为船的航速;α为水平翼的攻角。

      若要求得升力与升力矩,首先需确定水平翼的攻角$\alpha $,其表达式为

      $$\alpha = \varphi - {\gamma _ {\rm {b}}} + {\theta _ {\rm {F}}}$$ (3)

      式中:$\varphi $为T型翼摆角,本文规定摆角正方向与船纵摇方向相反;${\gamma _ {\rm {b}}}$为船纵摇角;${\theta _ {\rm {F}}}$为船体运动速度影响的附加角,如图1所示。图中:${\dot y_ {\rm {b}}}$为船体垂荡速度;$\dot \gamma $为纵摇角速度;$\dot \xi $为水平翼处流体质点垂直方向速度;${\theta _ {\rm {F}}}$为小量,其表达式为

      $${\theta _ {\rm {F}}} = \arctan \frac{{{l_ {\rm {F}}}{{\dot \gamma }_ {\rm {b}}} - {{\dot y}_ {\rm {b}}} + \dot \xi }}{U} = \frac{{{l_ {\rm {F}}}{{\dot \gamma }_ {\rm {b}}} - {{\dot y}_ {\rm {b}}} + \dot \xi }}{U}$$ (4)

      图  1  T型翼攻角分析图

      Figure 1.  Analysis of T-foil angle of attack

      由于水平翼处流体质点垂直方向速度$\dot \xi $为小量,其影响可以忽略不计,故最终得到T型翼的升力和升力矩的表达式为:

      $${F_ {\rm {T}}} = \frac{1}{2}\rho {U^2}A\frac{{{\rm{d}}{C_ {\rm {L}}}}}{{{\rm{d}}\alpha }}\left(\varphi - \gamma + \frac{{{l_ {\rm {F}}}\dot \gamma - \dot y}}{U}\right)$$ (5)
      $${M_ {\rm {T}}} = - \frac{1}{2}{l_ {\rm {F}}}\rho {U^2}A\frac{{{\rm{d}}{C_ {\rm {L}}}}}{{{\rm{d}}\alpha }}\left(\varphi - \gamma + \frac{{{l_ {\rm {F}}}\dot \gamma - \dot y}}{U}\right)$$ (6)

      在一定的弗劳德数下,确定水平翼给定摆角范围在$ \pm {\rm{1}}{{\rm{5}}^\circ }$之间,$\dfrac{1}{2}\rho {U^2}A\dfrac{{{\rm{d}}{C_ {\rm {L}}}}}{{{\rm{d}}\alpha }}$为一定值,设

      $${K_ {\rm {F}}} = \frac{1}{2}\rho {U^2}A\frac{{{\rm{d}}{C_ {\rm {L}}}}}{{{\rm{d}}\alpha }}$$ (7)
    • T型翼进行被动控制时,其摆角始终为某一固定角度${\varphi _ {\rm {fixed}}}$,可得被动控制下T型翼产生的升力与升力矩为:

      $${F_ {\rm {T}}} = {K_ {\rm {F}}}\left({\varphi _ {\rm {fixed}}} - \gamma + \frac{{{l_ {\rm {F}}}\dot \gamma - \dot y}}{U}\right)$$ (8)
      $${M_ {\rm {T}}} = - {K_ {\rm {F}}}{l_ {\rm {F}}}\left({\varphi _ {\rm {fixed}}} - \gamma + \frac{{{l_ {\rm {F}}}\dot \gamma - \dot y}}{U}\right)$$ (9)
    • T型翼的线性主动控制是指T型翼产生的垂向力和纵摇力矩与T型翼控制信号通过控制参数呈线性关系,因此其水平翼的摆角变化是连续的。采用纵摇角速度和艏部运动速度这2种信号分别对T型翼的摆角进行实时控制,进而实时控制T型翼产生的力矩。现将各种控制信号统称为$S$。船在刚开始运动时,T型翼为被动控制,此时监测$S$的最大值${S_{\max }}$、最小值${S_{\min }}$和幅值${S_ {\rm {a}}}$。最大值与最小值分别对应T型翼的摆角$ \pm {\varphi _{\max }}$。主动控制开始后,线性控制T型翼摆角为

      $$ {\varphi _ {\rm {L}}} = {\varphi _{\max }}\frac{S}{{{S_ {\rm {a}}}}}(- {\varphi _{{\rm{max}}}}\leqslant{\varphi _ {\rm {L}}}\leqslant{\varphi _{\max }}) $$ (10)

      线性控制控制力矩为

      $${M_ {\rm {T}}} = - {K_ {\rm {F}}}{l_ {\rm {F}}}\left({\varphi _ {\rm {L}}} - \gamma + \frac{{{l_ {\rm {F}}}\dot \gamma - \dot y}}{U}\right)$$ (11)
    • 本文数值计算采用的是Wigley III船型,这是一种细长的理想船型,如图2所示,船型参数如表1所示。

      T型翼(NACA 0018)纵向位置位于距船舯1.3 m的船底处。T型翼的详细尺寸如表2所示,三维图如图3所示。

      图  2  Wigley III船型

      Figure 2.  Hull form of Wigley III

      表 1  Wigley III船型参数

      Table 1.  Parameters of Wigley hull III

      参数数值
      垂线间长/cm300
      船宽/cm30
      设计吃水/cm18.75
      排水量/g78 000
      中横剖面系数0.666 7
      重心高度/cm0.17

      表 2  T型翼参数

      Table 2.  Parameters of the T-foil

      参数数值
      展长/cm9.0
      弦长/cm6.0
      垂直翼长度/cm5.0

      图  3  T型翼示意图

      Figure 3.  Schematic diagram of the T-foil

      在实际工程应用中,T型翼的摆角范围有严格的限制,一般在$ \pm 15^\circ $以内。因此,本文研究控制T型翼减摇的过程中,给定T型翼摆角范围为$ - 15^\circ < \varphi < 15^\circ $

    • Fr=0.3的航速及不同规则波($\lambda /L = $0.75,0.907,1.25,1.38,1.5,1.75,2;H=0.038 m)下对Wigley III船型的纵摇与垂荡运动的运动响应幅值(RAO)曲线进行数值计算,并将计算结果与试验值[11]进行对比,以验证本文数值方法的准确性。

      本文船的垂荡、纵摇和艏部运动加速度响应幅值均采用无量纲化计算:

      $${y^*} = \frac{y}{\zeta }$$ (12)
      $${\gamma ^*} = \frac{{\gamma \cdot L}}{{2 \cdot \ {\text{π}} \cdot \zeta }}$$ (13)
      $${b^ * } = \frac{b}{g}$$ (14)

      式中:$L$为船长;$g$为重力加速度;$b$为艏部运动加速度;${y^ * }$为无量纲垂荡幅值;${\gamma ^ * }$为无量纲纵摇幅值;${b^ * }$为无量纲艏部运动加速度;$\zeta $为波幅。通过计算以上7种工况下垂荡与纵摇运动响应,得到其运动响应幅值曲线如图4所示。

      图  4  数值计算结果与试验值的对比

      Figure 4.  Comparison of numerical results and the experimental results

      图4可见,数值计算的趋势与试验结果基本吻合,误差小于20%,说明本文数值模型计算结果能够与试验值很好地匹配,验证了本文数值模型的准确性。

    • 采用被动控制和线性主动控制(纵摇角速度信号和艏部运动速度信号)的Wigley III模型,对比不同工况下的垂荡、纵摇和艏部加速度幅值,并分析不同策略下的减摇效果。

    • 在线性主动控制过程中,T型翼的摆角是随控制信号呈线性变化的。对控制信号和T型翼摆角做归一化处理,得

      $${S^1} = \frac{S}{{{S_ {\rm {a}}}}}$$ (15)
      $${\varphi ^1} = \frac{\varphi }{{{\varphi _{{\rm{max}}}}}}$$ (16)

      式中:${S^1}$为归一化控制信号;${\varphi ^1}$为归一化T型翼摆角。

      Fr=0.3,$\lambda /L = 1.75$工况为例,对于以纵摇角速度为信号的线性控制,T型翼摆角和控制信号的对应关系如图5所示。

      图  5  T型翼摆角与纵摇角速度的对应关系(Fr=0.3,$\lambda /L = 1.75$

      Figure 5.  The relationship between T-foil angle and angle velocity of pitch(Fr=0.3,$\lambda /L = 1.75$

      在同样的工况下,对于以艏部运动速度为信号的线性控制,T型翼摆角和控制信号的对应关系如图6所示。

      图  6  T型翼摆角与艏部运动速度的对应关系(Fr=0.3,$\lambda /L $=1.75)

      Figure 6.  The relationship between T-foil angle and bow velocity (Fr=0.3,$\lambda /L$=1.75)

      船在航行过程中,T型翼始终处于船艏下方。图5图6中的纵摇角速度和艏部运动速度均是在被动控制下采集得到,符合加装T型翼船舶在运动过程中的实际情况。这也可以使更小的运动响应对应到更大的T型翼摆角,提高T型翼的减摇效果。

    • 2种航速(Fr=0.3,0.5)下的垂荡运动响应分别如图7图8所示。图中,${\omega _ {\rm {e}}}$为遭遇频率。表4表5所示为各种控制策略对垂荡运动的减摇效果。

      图  7  Fr=0.3时3种控制方式的垂荡运动响应

      Figure 7.  Heave motion response with three types of control method at Fr=0.3

      图  8  Fr=0.5时3种控制方式的垂荡运动响应

      Figure 8.  Heave motion response with three types of control method at Fr=0.5

      表 4  Fr=0.3时的垂荡运动减摇效果

      Table 4.  Damping effect of heave motion at Fr=0.3

      $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
      被动
      控制
      纵摇角
      速度信号
      艏部运动
      速度信号
      0.759.7823.0096.2087.15
      1.256.787.1911.1912.55
      1.755.37−2.60−11.422.28
      2.254.54−2.17−11.361.35
      2.753.97−1.67−10.931.90

      表 5  Fr=0.5时的垂荡运动减摇效果

      Table 5.  Damping effect of heave motion at Fr=0.5

      $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
      被动
      控制
      纵摇角
      速度信号
      艏部运动
      速度信号
      0.7512.81−5.7219.1140.46
      1.258.6055.1049.0273.72
      1.756.6736.4433.8526.88
      2.255.55−1.2911.3316.50
      2.754.80−3.514.6515.85

      在中速($Fr = 0.3$)情况下,裸船的垂荡运动响应峰值在${\omega _ {\rm {e}}}{\rm{ = }}6.78{\text{ rad/s}}$处,此时采用艏部运动速度信号进行控制具有最佳的减摇效果,可达12.55%;在${\omega _ {\rm {e}}}{\rm{ = }}9.78{\text{ rad/s}}$时,纵摇角速度信号具有最佳的减摇效果,可达96.20%;被动控制在高频区(${\omega _ {\rm {e}}} > 9{\text{ rad/s}}$)对垂荡运动的减摇效果相比另外2种均较弱,减摇效果不明显。在低频区(${\omega _ {\rm {e}}} < 5{\text{ rad/s}}$),3种控制的减摇效果均减弱了,被动控制和纵摇角速度信号的线性控制会使垂荡运动有所增加;艏部运动速度信号的线性控制在低频区减摇效果最佳。

      在高速($Fr = 0.5$)情况下,纵摇角速度信号和艏部运动速度信号的线性主动控制的减摇效果较被动控制有很大提高。裸船的垂荡运动响应峰值在${\omega _ {\rm {e}}}{\rm{ = }}8.60{\text{ rad/s}}$处,采用艏部运动速度信号进行控制具有最佳减摇效果,可达73.72%,纵摇角速度信号减摇效果为49.02%;在低频区(${\omega _ {\rm {e}}}{\rm{ = }}4.80{\text{ rad/s}}$时),艏部运动速度信号的减摇效果可达15.85%,纵摇角速度信号仅为4.65%。

      综上所述,对于垂荡运动的减摇效果,在中速和高速情况下,艏部运动速度信号具有最大的减摇效果。

    • 2种航速(Fr=0.3,0.5)下的纵摇运动响应分别如图9图10所示。表6表7所示为各种控制策略对纵摇运动的减摇效果。

      图  9  Fr=0.3时3种控制方式的纵摇运动响应

      Figure 9.  Pitch motion response with three types of control method at Fr=0.3

      图  10  Fr=0.5时3种控制方式的纵摇运动响应

      Figure 10.  Pitch motion response with three types of control method at Fr=0.5

      表 6  Fr=0.3时的纵摇运动减摇效果

      Table 6.  Damping effect of pitch motion at Fr=0.3

      $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
      被动
      控制
      纵摇角
      速度信号
      艏部运动
      速度信号
      0.759.783.7037.9929.44
      1.256.7827.9140.4244.98
      1.755.379.3724.0839.82
      2.254.547.2422.7147.47
      2.753.976.8827.8256.34

      表 7  Fr=0.5时的纵摇运动减摇效果

      Table 7.  Damping effect of pitch motion at Fr=0.5

      $\lambda /L$${\omega _ {\rm {e}}}$/(rad·s−1)减摇效果/%
      被动
      控制
      纵摇角
      速度信号
      艏部运动
      速度信号
      0.7512.81−5.0568.7249.70
      1.258.6033.3879.7831.16
      1.756.6765.0482.7887.33
      2.255.5518.8662.1583.19
      2.754.8011.6662.1286.05

      在中速($Fr = 0.3$)情况下,裸船的纵摇运动响应峰值在${\omega _ {\rm {e}}}{\rm{ = }}6.78\;{\rm rad/s}$处,艏部运动速度信号具有最佳减摇效果,可达44.98%。3种控制方式下,船纵摇运动的变化趋势大体相似,在峰值处减摇效果最大,远离峰值处减摇效果逐渐降低。在低频区(${\omega _ {\rm {e}}}{\rm{ = }}4.54,\;3.97\;{\rm rad/s}$),艏部运动速度信号的减摇效果可达47.47%和56.34%,纵摇角速度的减摇效果在20%~30%之间,被动控制的减摇效果则低于10%。

      在高速($Fr = 0.5$)情况下,裸船的纵摇运动响应峰值在${\omega _ {\rm {e}}}{\rm{ = }}6.67{\text{ rad/s}}$处,艏部运动速度信号具有最佳减摇效果,可达87.33%,纵摇角速度信号次之,可达82.78%。在高频区(${\omega _ {\rm {e}}}{\rm{ = }}12.81,\;8.60\;{\rm rad/s}$),纵摇角速度信号的减摇效果最为明显,约达70%;在低频区(${\omega _ {\rm {e}}}{\rm{ = }}5.55,\;4.80{\text{ rad/s}}$),艏部运动速度信号具有最佳减压效果,达80%以上。

      综上所述,对于纵摇运动的减摇效果,在2种航速下,在高频区时纵摇角速度的减摇效果更加显著,在低频区则是艏部运动速度信号的减摇效果更佳。

    • 2种航速(Fr=0.3,0.5)下的艏部运动加速度运动响应分别如图11图12所示。表8表9所示为各种控制策略对艏部运动加速度运动的减摇效果。

      图  11  Fr=0.3时3种控制方式的艏部运动加速度运动响应

      Figure 11.  Bow acceleration motion response with three types of control method at Fr=0.3

      图  12  Fr=0.5时3种控制方式的艏部运动加速度运动响应

      Figure 12.  Bow acceleration motion response with three types of control method at Fr=0.5

      表 8  Fr=0.3时的艏部运动加速度减摇效果

      Table 8.  Damping effect of bow acceleration motion at Fr=0.3

      $\lambda /L$${\omega _ {\rm {e}} }$/(rad·s−1)减摇效果/%
      被动
      控制
      纵摇角
      速度信号
      艏部运动
      速度信号
      0.759.789.7157.7647.37
      1.256.7824.2334.1141.36
      1.755.373.914.4327.80
      2.254.541.38−5.0426.79
      2.753.970.85−9.4027.45

      表 9  Fr=0.5时的艏部运动加速度减摇效果

      Table 9.  Damping effect of bow motion acceleration at Fr=0.5

      $\lambda /L$${\omega _ {\rm {e}} }$/(rad·s−1)减摇效果/%
      被动
      控制
      纵摇角
      速度信号
      艏部运动
      速度信号
      0.7512.813.9186.7574.00
      1.258.6047.5565.0876.32
      1.756.6757.6472.6879.27
      2.255.554.3629.0858.50
      2.754.80−0.9518.4058.45

      在中速($Fr = 0.3$)情况下,裸船的艏部运动加速度峰值出现在${\omega _ {\rm {e}}}{\rm{ = }}6.78{\text { rad/s}}$处。除高频区(${\omega _ {\rm {e}}}{\rm{ = }}9.78{\text { rad/s}}$)外,艏部运动速度信号在另外4种波长下均具有最佳减摇效果,在${\omega _ {\rm {e}}}{\rm{ = }}3.97{\text { rad/s}}$时,减摇效果可达27.45%,远高于另外2种控制方式。纵摇角速度信号在高频区(${\omega _ {\rm {e}}}{\rm{ = }}9.78{\text { rad/s}}$)减摇效果较好,但在低频区(${\omega _ {\rm {e}}}{\rm{ = }}4.54,\;3.97{\text { rad/s}}$)则有一定程度的减益影响。

      在高速($Fr = 0.5$)情况下,裸船的艏部运动加速度峰值出现在${\omega _{\rm {e}}}{\rm{ = }}6.67{\text { rad/s}}$处。与低速时情况一样,除高频区(${\omega _{\rm {e}}}{\rm{ = }}12.81{\text { rad/s}}$)外,艏部运动速度信号在另外4种波长下均具有最佳的艏部运动加速度减摇效果,在${\omega _{\rm {e}}}{\rm{ = }}4.80{\text { rad/s}}$时,减摇效果可达58.45%,而纵摇角速度信号的减摇效果则仅为18.40%。

      综上所述,对于艏部运动加速度,艏部运动速度信号的减摇效果和纵摇角速度相比有很大的提升,特别是在低频区,减摇效果十分显著。

    • 船在运动过程中,T型翼对船产生的升力和纵升力矩与船的垂荡及纵摇运动方向相反,这样的力和力矩可以抑制船的纵向运动响应。在计算过程中,针对T型翼对船产生的升力及升力矩进行了监测。升力矩是升力与水平翼尾缘距船舯距离的乘积,与升力具有同样的变化形式。本节将以Fr=0.5时,$\lambda /L = 1.75$,2.25这2工况为例进行分析。2种主动控制T型翼对船产生的升力矩如图13所示。

      图  13  2种T型翼主动控制策略在2种工况下对船产生的纵摇力矩

      Figure 13.  The pitching moment of two type T-foil active control strategies for the ship in two cases

      图13中,T型翼在第6个波浪周期后开始进行主动运动控制,在此之前,为被动控制。进行主动控制之后,T型翼对船产生的升力矩有大幅度的增长,说明主动控制相比被动控制具有更大的升力矩可用于抑制船舶的纵向运动。

      在T型翼主动控制减摇过程中,所产生的升力与升力矩是随控制信号变化的,如图14图15所示。图中,工况为Fr=0.5,$\lambda /L = 1.75$的升力矩与运动响应均做了归一化处理。

      图  14  采用艏部运动速度控制方法的升力矩与纵摇角速度和垂荡速度的变化关系

      Figure 14.  The relationship between lifting moment and pitching angle velocity and vertical velocity controlled by the heading velocity control method

      图  15  采用纵摇角速度控制方法的升力矩与纵摇角速度和垂荡速度的变化关系

      Figure 15.  The relationship between lifting moment and angular velocity in pitch and heave velocity controlled by the pitch angle velocity control method

      图13所示,T型翼产生的升力矩为正弦变化形式,有幅值、频率和相位3种主要参数。在Fr=0.5,$\lambda /L = 1.75$和2.25这2个工况下,艏部运动速度信号相比纵摇角速度信号具有更好的减摇效果,并且在同一工况下采用不同的控制策略产生的升力矩幅值与频率一致。如图14图15所示,艏部运动速度信号控制策略所产生的升力矩与纵摇速度和垂荡速度在相位上均相差约1/8个周期,纵摇角速度信号控制策略所产生的升力矩与纵摇角速度无相位差,但与垂荡速度相差约1/4个相位,导致在1/4个周期内T型翼产生的升力矩与垂荡速度方向同向,故对垂荡与艏部运动的减摇有减益效果。而艏部运动速度信号为纵摇与垂荡这2种运动的耦合信号,具有更佳的减摇效果。

    • 本文采用细长体理论方法计算了Wigley III船型在规则波上的运动,并分别采用被动控制、纵摇角速度信号线性控制和艏部运动速度线性控制这3种控制方法,计算并分析了其对船舶纵向运动的减摇效果,主要得到如下结论:

      1) 在共振区,被动控制T型翼对垂荡运动、纵摇运动和艏部运动加速度具有较好的减摇效果,但在高频区和低频区,其减摇效果相对有限。在低频区,被动控制对纵摇与艏部运动加速度的减摇效果约为5%~10%。

      2) 纵摇角速度信号线性主动控制解决了被动控制在低频区减摇效果不理想的问题。在高速情况下,纵摇运动的减摇效果最佳可达82.78%,在低频区其减摇效果也可达62.12%,艏部运动加速度的减摇效果最佳可达72.68%。但是在中速低频区,纵摇角速度信号对垂荡运动和艏部运动加速度均有一定程度的减益,为约−10%。

      3) 艏部运动速度信号线性主动控制综合考虑了纵摇与垂荡这2种纵向运动,与被动控制和纵摇角速度信号相比,在低频区,垂荡运动的减摇效果有较大改善,没有出现减益现象。高航速下,艏部速度信号线性主动控制在各个工况下都具有很好的减摇效果,在低频区的减摇效果要优于纵摇角速度信号和被动控制。在高速低频区,纵摇运动的减摇效果可达86.05%。对于艏部运动加速度,在中、高速情况下,艏部运动速度信号的减摇效果较纵摇角速度信号有很大的提升,平均提高了10%~20%。

      综上所述,以艏部运动速度信号对T型翼进行线性控制对于垂荡、纵摇与艏部运动加速度具有最佳的减摇效果。

参考文献 (11)

目录

    /

    返回文章
    返回