• 匿名盲审
  • 学术期刊非营利性
  • 全球免费开放获取全文
  • 最新科研成果提供绿色通道

留言板

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

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

漂浮式海上风机平台阻尼结构设计与研究

陈嘉豪 尹梓炜 朱容宽 马兆荣 景雪娇 李炎

陈嘉豪,尹梓炜,朱容宽等.漂浮式海上风机平台阻尼结构设计与研究[J].南方能源建设,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003
引用本文: 陈嘉豪,尹梓炜,朱容宽等.漂浮式海上风机平台阻尼结构设计与研究[J].南方能源建设,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003
CHEN Jiahao,YIN Ziwei,ZHU Rongkuan,et al.Design and Research on the Foundation Damping Structure of Semi-Submersible Floating Offshore Wind Turbines[J].Southern Energy Construction,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003
Citation: CHEN Jiahao,YIN Ziwei,ZHU Rongkuan,et al.Design and Research on the Foundation Damping Structure of Semi-Submersible Floating Offshore Wind Turbines[J].Southern Energy Construction,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003

漂浮式海上风机平台阻尼结构设计与研究

doi: 10.16516/j.gedi.issn2095-8676.2021.03.003
基金项目: 

中国博士后基金资助项目“台风工况下的漂浮式海上风机动力特性及抗台策略研究” 2020M682995

三亚崖州湾科技城重大研发项目“漂浮式基础平台风机关键技术研究” SKJC-2020-01-005

中国能建广东院科技项目“10 MW级半潜型漂浮式海上风机耦合动力特性与安全性评价研究” EV05241W

详细信息
    作者简介:

    陈嘉豪(通信作者)1991-,男,广东惠州人,上海交通大学工学博士,主要研究领域包括:海上漂浮式风机耦合动力理论研究、相关设计、模型试验技术以及相关数值软件开发(e-mail)chenjiahao@gedi.com.cn

    尹梓炜1993-,男,广东云浮人,华南理工大学船舶与海洋工程硕士,主要从事水工结构研究及设计(e-mail)yinziwei@gedi.com.cn

    朱容宽1994-,女,江西南昌人,南京航空航天大学结构工程硕士,主要从事水工结构研究及设计(e-mail)zhurongkuan@gedi.com.cn

  • 中图分类号: TK89

Design and Research on the Foundation Damping Structure of Semi-Submersible Floating Offshore Wind Turbines

  • CHEN Jiahao,YIN Ziwei,ZHU Rongkuan,et al.Design and Research on the Foundation Damping Structure of Semi-Submersible Floating Offshore Wind Turbines[J].Southern Energy Construction,2021,08(03):18-25.
  • 摘要:   目的  为了寻求经济有效的漂浮式风机水动力性能提升方案,研究了加装垂荡板的设计方案,从而提升浮体水动力阻尼性能和减摇效果,最终改善漂浮式海上风电机输出功率的稳定性。  方法  以10 MW半潜型浮式风机为例,通过计算流体力学的方法,探讨了垂荡板结构不同设计参数对于提升漂浮式风机基础水动力阻尼的效用,寻找出优化设计方案。  结果  经过优化设计后的垂荡板结构垂荡阻尼较原结构增大14.9%,摇摆阻尼增大19.1%,而“梅花形”边缘垂荡板垂荡运动阻尼较原始模型提升36.98% 。  结论  揭示泄涡产生与垂荡板水动力阻尼改善的形成机理,创新性地提出“梅花形”垂荡板结构,有效改善浮体垂荡运动的阻尼性能,上述结论为此类漂浮式海上风机基础结构的研究与设计工作提供重要参考。
  • 图  1  海上风机基础形式3

    Fig.  1  Foundation form of offshore wind turbines3

    图  4  模型计算域及边界条件

    Fig.  4  Computational region and boundary conditions

    图  5  计算域截面网格(422万网格)

    Fig.  5  Section mesh of the computational region

    图  6  边锋延伸宽度对基础垂荡阻尼的影响

    Fig.  6  Damping variation of heave motion from the wing extension of heave-plate

    图  7  边锋延伸宽度对基础横摇阻尼的影响

    Fig.  7  Damping variation of roll motion from the wing extension of heave-plate

    图  8  垂荡板高度变化对基础垂荡阻尼的影响

    Fig.  8  Damping variation of heave motion from the height of heave-plate

    图  9  垂荡板高变化时基础横摇阻尼的影响

    Fig.  9  Damping variation of roll motion from the height of heave-plate

    图  10  垂荡板圆盘直径对基础垂荡阻尼的影响

    Fig.  10  Damping variation of heave motion from the diameter of heave-plate

    图  11  垂荡板圆盘直径对基础横摇阻尼的变化

    Fig.  11  Damping variation of roll motion from the diameter of heave-plate

    图  12  浮式风机平台垂向运动涡量图

    Fig.  12  Vorticity contours of the foundation of the floating offshore wind turbine in heave

    图  13  “梅花形”整体模型示意图

    Fig.  13  Schematic diagram of the whole model of cinquofoil

    图  14  模型强迫垂荡运动时历曲线

    Fig.  14  Schematic diagram of the whole model of cinquefoil

    表  1  OO-Star 漂浮式风机基础参数

    Tab.  1.   Basic parameters of the OO-Star floating wind turbines

    名称数值
    基础质量/kg(不包括塔筒和系泊,含压载)2.170 9×107
    重心距水面距离/m-15.225
    横摇惯性矩/(kg·m29.43×109
    纵摇惯性矩/(kg·m29.43×109
    首摇惯性矩/(kg·m21.63×1010
    塔柱低端距水面距离/m11.0
    平衡位置吃水(包括系泊系统)/m22.0
    排水量/m32.350 9×104
    浮心距水面距离/m14.236
    下载: 导出CSV

    表  2  阻尼结构设计方案

    Tab.  2.   Design schemes of the heave-plate structure

    方案直径D/m边锋外伸S/mH/m
    原模型22.800.5
    模型122.80.20.5
    模型222.80.40.5
    模型322.80.60.5
    模型422.80.20.6
    模型522.80.20.7
    模型623.80.20.5
    模型724.80.20.5
    下载: 导出CSV

    表  3  不同网格数的计算结果

    Tab.  3.   Results with different mesh quantities

    网格数/万个附加质量阻尼系数/(N·s·m-1无因次阻尼系数变化率/%
    21210 236 877 9752 249 987 2620.114 4-2.54
    42210 036 272 2562 285 078 9000.117 40
    82110 148 250 0162 268 521 3990.115 9-1.29
    下载: 导出CSV

    表  4  测试工况参数

    Tab.  4.   Parameters in test cases

    运动形式发电工况
    周期/s幅值/m
    垂荡7.51.75
    横摇7.51.29
    下载: 导出CSV

    表  5  边锋延伸宽度对阻尼的影响

    Tab.  5.   Damping impact from the wing extension

    模型及运动无因次垂荡阻尼系数增长率/%无因次横摇阻尼系数增长率/%
    原模型0.175 300.103 30
    模型1(延伸0.2 m)0.176 50.680.106 32.95
    模型2(延伸0.4 m)0.183 74.790.106 02.61
    模型3(延伸0.6m)0.192 39.700.110 46.90
    下载: 导出CSV

    表  6  垂荡板高度变化时风机基础阻尼的影响

    Tab.  6.   Damping impact from height of heave-plate

    模型及运动无因次垂荡阻尼系数增长率/%无因次横摇阻尼系数增长率/%
    原模型(高0.5 m)0.175 300.103 30
    模型1(高0.5 m)0.176 50.680.106 32.90
    模型4(高0.6 m)0.179 92.620.104 30.93
    模型5(高0.7 m)0.187 67.010.107 74.26
    下载: 导出CSV

    表  7  圆盘直径变化时风机基础垂荡阻尼的变化

    Tab.  7.   Damping impact from the diameter of heave-plate

    模型及运动无因次垂荡阻尼系数增长率/%无因次横摇阻尼系数增长率/%
    原模型(直径22.8 m)0.175 300.103 30
    模型1(直径22.8 m)0.176 50.680.106 43.00
    模型6(直径23.8 m)0.184 75.360.117 413.64
    模型7(直径24.8 m)0.201 815.120.121 217.32
    下载: 导出CSV

    表  8  优化后模型阻尼性能对比

    Tab.  8.   Damping comparison between the original model and optimized model

    模型及运动无因次阻尼系数阻尼系数增长率/%
    有量纲无量纲
    原模型(高0.5 m)0.175 3
    模型1(高0.5 m)0.201 426.17514.9
    模型4(高0.6 m)0.103 3
    模型5(高0.7 m)0.123 028.45519.1
    下载: 导出CSV

    表  9  “梅花形”垂荡板阻尼性能计算结果

    Tab.  9.   Calculation results of damping performance of "quincunx-shaped" heave-plate

    模型及运动无因次阻尼系数阻尼系数增长率/%
    原模型垂荡0.175 3
    优化后模型垂荡0.201 414.90
    “梅花形”垂荡0.204 536.98
    下载: 导出CSV
  • [1] GWEC. GLOBAL WIND REPORT 2018 [R]. Brussels: GWEC, 2019.
    [2] JEONS H, CHOY U, SEOM W, et al. Dynamic response of floating substructure of spar-type offshore wind turbine with catenary mooring cables [J]. Ocean Engineering, 2013, 72(72), 356-364.
    [3] 陈嘉豪,裴爱国,马兆荣,等.海上漂浮式风机关键技术研究进展 [J].南方能源建设,2020,7(1): 8-20.

    CHENJ H, PEIA G, MAZ R, et al. A review of the key technologies for floating offshore wind turbines [J]. Southern Energy Construction, 2020, 7(1): 8-20.
    [4] SKAAREB, HANSONT D, YTTERVIKR, et al. Dynamic response and control of the hywind demo floating wind turbine [C]// EWEA. Proceedings of Europe Wind Energy Conference, Brussels, Belgium, 2011. Belgium: EWEA, 2011: 43.
    [5] AUBAULTA, CERMELLIC, LAHIJANIANA, et al. Wind float contraption: from conception to reproduction [C]// Ocean, Offshore and Arctic Engineering Division: the American. Proceedings of the 31st International Conference on Offshore Mechanics and Arctic Engineering, Rio de Janeiro, Brazil, 2012. Rio de Janeiro: OMAE, 2012: 847-853.
    [6] BEYERF, CHOISNETT, KRETSCHMERM, et al. Coupled MBS-CFD simulation of the IDEOL floating offshore wind turbine foundation compared to wave tank model test data [C]// ISOPE. Proceedings of ISOPE. 25th International Ocean and Polar Engineering Conference, Kona, Big Island, HI, USA, Jun. 21-26, 2015. Kona: ISOPE, 2015: 367-374.
    [7] European Technology & Innovation Platform on Wind Energy. Strategic research and innovation agenda 2018 [R]. Lisbon: Etipwind.eu, 2019.
    [8] THIAGARAJANK, TROESCHA. Effects of appendages and small currents on the hydrodynamic heave damping of TLP columns [J]. J. Offshore Mech. Arct. Eng., 1998(120): 37-42.
    [9] 吴维武, 缪泉明, 匡晓峰, 等. Spar平台垂荡板受迫振荡水动力特性研究 [J]. 船舶力学, 2009, 13(1): 27-33.

    WUW W, LIAOQ M, KUANGX F, et al. Research on hydrodynamic characteristics of forced oscillation heave damping plates of Spar platforms [J]. Journal of Ship Mechanics, 2009,13(1): 27-33.
    [10] WANGY, CHENH C, VAZ G, et al. CFD simulation of semi-submersible floating offshore wind turbine under pitch decay motion [C]// Ocean, Offshore and Arctic Engineering Division. ASME 2019 2nd International Offshore Wind Technical Conference, Malta, 2019. Malta : IOWTC, 2019: 11.
    [11] PEGALAJAR-JURADOA, BREDMOSEH, BORGM, et al. State-of-the-art model for the LIFES50+ OO-Star Wind Floater Semi 10 MW floating wind turbine [C]// SINTEF, NTNU. Journal of Physics: Conference Series, Trondheim, Norway, 2018. Norway: IOP Publishing, 2018, 1104(1): 012024.
    [12] TAOL, MOLINB, SCOLANY, et al. Spacing effects on hydrodynamics of heave plates on offshore structures [J]. J. Fluids Struct., 2007, 23(8): 1119–1136.
  • [1] 赵业彬, 任建宇, 乐丛欢.  风荷载对全潜式浮式风机拖航运动响应的影响 . 南方能源建设, 2024, 11(1): 176-184. doi: 10.16516/j.ceec.2024.1.18
    [2] QI Xiling, LE Conghuan, REN Jianyu, ZHANG Puyang, DING Hongyan, QI Xiaoliang.  Comparision of Integrated Towing Motion Performance for Multi-Column Semisubmersible Floating Offshore Wind Turbines . 南方能源建设, 2024, 11(2): 17-30. doi: 10.16516/j.ceec.2024.2.02
    [3] 王立鹤, 符妃.  海上风机基础与网箱融合系统静力计算 . 南方能源建设, 2024, 11(4): 76-87. doi: 10.16516/j.ceec.2024.4.08
    [4] 罗贤勇, 施海云, 吴家凯.  海上浮动式核电站汽轮发电机组方案研究 . 南方能源建设, 2024, 11(6): 1-8. doi: 10.16516/j.ceec.2024.6.09
    [5] 陈嘉豪, 高一帆, 尹梓炜, 郑灿, 曲晓奇.  近海深水区混凝土半潜型浮式风机一体化计算与耦合动力特性分析 . 南方能源建设, 2024, 11(2): 31-41. doi: 10.16516/j.ceec.2024.2.03
    [6] 曾雨欣, 施伟, 张礼贤, 周昳鸣.  10 MW大型单桩式海上风机桩土作用研究 . 南方能源建设, 2023, 10(1): 1-12. doi: 10.16516/j.gedi.issn2095-8676.2023.01.001
    [7] 傅禹舜, 施伟, 张松浩, 张礼贤, 任亚君, 李昕, 周昳鸣.  强非线性波作用下海上风机动力特性模型试验研究 . 南方能源建设, 2023, 10(4): 11-17. doi: 10.16516/j.gedi.issn2095-8676.2023.04.002
    [8] 郑阳.  圆形扩展式风机基础钢筋设计优化分析 . 南方能源建设, 2022, 9(S1): 76-82. doi: 10.16516/j.gedi.issn2095-8676.2022.S1.012
    [9] 周林, 施伟, 张松浩, 张礼贤.  破碎波作用下单桩式海上风机水动力学数值分析 . 南方能源建设, 2020, 7(3): 70-80. doi: 10.16516/j.gedi.issn2095-8676.2020.03.009
    [10] 陈嘉豪, 裴爱国, 马兆荣, 庞程燕.  海上漂浮式风机关键技术研究进展 . 南方能源建设, 2020, 7(1): 8-20. doi: 10.16516/j.gedi.issn2095-8676.2020.01.002
    [11] 刘飞虹, 吴迪, 苗得胜, 张敏, 黎建峰.  发电量预测误差受地形图精度影响分析 . 南方能源建设, 2019, 6(S1): 41-46. doi: 10.16516/j.gedi.issn2095-8676.2019.S1.009
    [12] 朱松晔, 朱子默, 柯世堂.  海上风机灾害分析及健康监测技术综述 . 南方能源建设, 2018, 5(2): 47-59. doi: 10.16516/j.gedi.issn2095-8676.2018.02.006
    [13] 王洪庆, 刘旭东, 毕明君, 刘晋超.  海上风机单桩基础疲劳影响因素分析 . 南方能源建设, 2018, 5(4): 92-97. doi: 10.16516/j.gedi.issn2095-8676.2018.04.014
    [14] 胡琨.  水环式真空泵工作液冷却水系统设计优化 . 南方能源建设, 2018, 5(S1): 129-132. doi: 10.16516/j.gedi.issn2095-8676.2018.S1.023
    [15] 胡传鹏.  厂房零米设备基础简化动力计算适用性研究 . 南方能源建设, 2018, 5(S1): 133-139. doi: 10.16516/j.gedi.issn2095-8676.2018.S1.024
    [16] 施伟, 郑侃, 任年鑫.  南海海况下半潜浮式风机在故障工况下的动力学响应分析 . 南方能源建设, 2018, 5(4): 12-20. doi: 10.16516/j.gedi.issn2095-8676.2018.04.002
    [17] 毕明君.  海上风机单桩基础选型设计方法 . 南方能源建设, 2017, 4(S1): 56-61,72. doi: 10.16516/j.gedi.issn2095-8676.2017.S1.011
    [18] 何小华, 张力.  海上风机基础阴极保护设计与应用 . 南方能源建设, 2015, 2(3): 76-79. doi: 10.16516/j.gedi.issn2095-8676.2015.03.015
    [19] 陈章隆.  基于Hermite插值法的板型燃料不对称冷却计算 . 南方能源建设, 2015, 2(4): 132-137. doi: 10.16516/j.gedi.issn2095-8676.2015.04.024
    [20] 张灿亨, 王鑫, 郑文棠.  基于水动力学模型的暴雨浅层滑坡风险分析 . 南方能源建设, 2015, 2(S1): 198-202. doi: 10.16516/j.gedi.issn2095-8676.2015.S1.044
  • 加载中
图(17) / 表 (9)
计量
  • 文章访问数:  966
  • HTML全文浏览量:  247
  • PDF下载量:  289
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-08-31
  • 修回日期:  2021-09-12
  • 刊出日期:  2021-09-25

漂浮式海上风机平台阻尼结构设计与研究

doi: 10.16516/j.gedi.issn2095-8676.2021.03.003
    基金项目:

    中国博士后基金资助项目“台风工况下的漂浮式海上风机动力特性及抗台策略研究” 2020M682995

    三亚崖州湾科技城重大研发项目“漂浮式基础平台风机关键技术研究” SKJC-2020-01-005

    中国能建广东院科技项目“10 MW级半潜型漂浮式海上风机耦合动力特性与安全性评价研究” EV05241W

    作者简介: 陈嘉豪(通信作者)1991-,男,广东惠州人,上海交通大学工学博士,主要研究领域包括:海上漂浮式风机耦合动力理论研究、相关设计、模型试验技术以及相关数值软件开发(e-mail)chenjiahao@gedi.com.cn

    尹梓炜1993-,男,广东云浮人,华南理工大学船舶与海洋工程硕士,主要从事水工结构研究及设计(e-mail)yinziwei@gedi.com.cn

    朱容宽1994-,女,江西南昌人,南京航空航天大学结构工程硕士,主要从事水工结构研究及设计(e-mail)zhurongkuan@gedi.com.cn

  • 中图分类号: TK89

摘要:   目的  为了寻求经济有效的漂浮式风机水动力性能提升方案,研究了加装垂荡板的设计方案,从而提升浮体水动力阻尼性能和减摇效果,最终改善漂浮式海上风电机输出功率的稳定性。  方法  以10 MW半潜型浮式风机为例,通过计算流体力学的方法,探讨了垂荡板结构不同设计参数对于提升漂浮式风机基础水动力阻尼的效用,寻找出优化设计方案。  结果  经过优化设计后的垂荡板结构垂荡阻尼较原结构增大14.9%,摇摆阻尼增大19.1%,而“梅花形”边缘垂荡板垂荡运动阻尼较原始模型提升36.98% 。  结论  揭示泄涡产生与垂荡板水动力阻尼改善的形成机理,创新性地提出“梅花形”垂荡板结构,有效改善浮体垂荡运动的阻尼性能,上述结论为此类漂浮式海上风机基础结构的研究与设计工作提供重要参考。

English Abstract

陈嘉豪,尹梓炜,朱容宽等.漂浮式海上风机平台阻尼结构设计与研究[J].南方能源建设,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003
引用本文: 陈嘉豪,尹梓炜,朱容宽等.漂浮式海上风机平台阻尼结构设计与研究[J].南方能源建设,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003
CHEN Jiahao,YIN Ziwei,ZHU Rongkuan,et al.Design and Research on the Foundation Damping Structure of Semi-Submersible Floating Offshore Wind Turbines[J].Southern Energy Construction,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003
Citation: CHEN Jiahao,YIN Ziwei,ZHU Rongkuan,et al.Design and Research on the Foundation Damping Structure of Semi-Submersible Floating Offshore Wind Turbines[J].Southern Energy Construction,2021,08(03):18-25. doi:  10.16516/j.gedi.issn2095-8676.2021.03.003
  • 由于传统化石能源的过度开采与使用导致气候变化和环境污染等问题日益严重,使得能源开发转向可再生的绿色能源。近年来,随着大容量机组技术的日渐成熟,以海上风电为代表的绿色能源在国内外蓬勃发展。GWEC市场展望,预计到2023年,海上风电新装机容量将会超过55 GW1,并将助力全球碳达峰目标的实现。

    随着近岸优质风场资源逐渐被开发殆尽,海上风电未来将不断走向深远海。在过去,近海浅水风场主要以着床式的固定式风机基础为主,包括单桩、导管架等,这些基础通常需要进行海上桩基施工作业,基础设计频域介于1倍到3倍风轮转子频率之间,随着水深的进一步增加,基础造价和施工成本将急剧上升2。相比而言,漂浮式海上风机则属于顺应式柔性结构,通过系泊与海床连接,基础造价和施工安装成本随水深变化的敏感度较传统固定式风机基础低,未来有望通过技术进步和产业链的发展,进一步提高其经济竞争力3。漂浮式风机的基础形式借鉴于过往船舶与海洋工程的相关浮体经验,可大致划分为半潜型式,立柱型式,张力腿型式和驳船型式3,如图1所示。目前,在欧洲和日本已经有部分漂浮式样机和小规模商业化漂浮式风场建设投产4-6。ETI预测,漂浮式风机基础技术将在亚洲、美国和欧洲具有广阔的市场空间和发展前景7

    然而,相比于传统的海上风机固定式基础,漂浮式风机在风浪环境下,自身运动幅度、速度和加速度都比较大,容易造成风轮入流风速的剧烈变化,湍流度增加,给功率稳定性和结构安全带来较大的挑战。因此,如何以较低的成本有效地提升漂浮式风机基础的运动阻尼,从而优化基础运动特性,是一个重要且现实的研究课题之一。

    传统船舶减摇方式,包括减摇鳍、舭龙骨、减摇水舱等,而油气平台则更普遍地采用垂荡板结构。垂荡板是一种简单而有效的浮体水动力性能调节装置,一方面可以通过改变浮体的运动附加质量,从而调节浮体的垂荡固有周期,另一方面则通过垂荡板边缘泄涡而增加系统的垂荡阻尼。漂浮式风机基础作为无动力浮体结构,可参考借鉴油气平台的垂荡板结构,探索垂荡板结构对漂浮式风机的阻尼调节作用。在以往,已有一些研究开始关注垂荡板的水动力特性。Thiagarajan和Troesch8通过对裸圆柱体进行水池强迫试验,表明试验柱体底部的圆盘能有效的产生垂荡阻尼,并且阻尼大小与震荡幅度有关。吴维武等人9采用计算流体力学软件Fluent对不同形状和开孔率的垂荡板的水动力特性进行了研究,结果表明不同的KC数范围,开孔对垂荡板阻尼的影响不尽相同。目前,一些半潜型漂浮式风机的平台设计也借鉴垂荡板,以优化平台运动性能。Wang等人10通过CFD方法对OC4-5 MW半潜型浮式风机的纵摇运动进行了水动力研究,结果表明平台纵摇运动速度最大的时候,垂荡板边缘泄涡产生的阻尼也达到极大值,有效地提升水动力阻尼。

    虽然部分研究者已经开始关注到了垂荡板在优化漂浮式风机运动性能的潜力,但是关于垂荡板的设计及漂浮式风机阻尼的改善特性研究仍然不足。因此,本文将针对国内海域普遍适用的半潜型浮式风机开展垂荡板阻尼结构的设计以及阻尼性能的定量化研究,为后续的工程设计提供有益的参考。

    • 本文采用的漂浮式风机基础为OO-Star漂浮式风机,该浮式风机基础是一个“Y”型浮筒结构,中间设计锥形立柱结构以支撑DTU 10 MW风机,外围设置三个立柱结构,具体结构形式和尺寸如图 2所示。浮式基础的原始设计水深为130 m,吃水为22 m,排水量2.350 9×104 m3。具体结构参数见表1,更多的基础细节可参考文献11

      图  2  漂浮式风机基础结构视图

      Figure 2.  Dimension of the OO-Star floating offshore wind turbines

      图  2  漂浮式风机基础结构视图

      Figure 2.  Dimension of the OO-Star floating offshore wind turbines

      图  2  漂浮式风机基础结构视图

      Figure 2.  Dimension of the OO-Star floating offshore wind turbines

      表 1  OO-Star 漂浮式风机基础参数

      Table 1.  Basic parameters of the OO-Star floating wind turbines

      名称数值
      基础质量/kg(不包括塔筒和系泊,含压载)2.170 9×107
      重心距水面距离/m-15.225
      横摇惯性矩/(kg·m29.43×109
      纵摇惯性矩/(kg·m29.43×109
      首摇惯性矩/(kg·m21.63×1010
      塔柱低端距水面距离/m11.0
      平衡位置吃水(包括系泊系统)/m22.0
      排水量/m32.350 9×104
      浮心距水面距离/m14.236
    • 本文在原有垂荡板的基础上,改进垂荡板的结构型式,包括垂荡板外伸尺寸、垂荡板的甲板和底板外伸(本文称为边锋),见图3

      图  3  阻尼结构示意图

      Figure 3.  Diagram of the heave-plate structure

      图  3  阻尼结构示意图

      Figure 3.  Diagram of the heave-plate structure

      原模型垂荡板圆盘直径22.8 m,高0.5 m,无边锋延伸。为了研究垂荡板形状参数对浮体水动力阻尼的影响,设置对照模型。模型垂荡板的圆盘直径D分别为22.8 m、23.8 m、24.8 m;边锋延伸部分宽度S分别0.2 m、0.4 m、0.6 m;垂荡板高H分别为0.5 m、0.6 m 、0.7 m,立柱参数不变。计算分析阻尼结构的运动衰减特性。阻尼结构设计方案见表2

      表 2  阻尼结构设计方案

      Table 2.  Design schemes of the heave-plate structure

      方案直径D/m边锋外伸S/mH/m
      原模型22.800.5
      模型122.80.20.5
      模型222.80.40.5
      模型322.80.60.5
      模型422.80.20.6
      模型522.80.20.7
      模型623.80.20.5
      模型724.80.20.5
    • 以垂荡运动为例,给定浮体运动形式为a=Asinωt,其中A为浮体运动幅值,ω为波浪频率。根据CFD计算所得浮体受力曲线,可整理得到浮体受力的数学表达形式为F=F0sinωt+φ,其中受力幅值F0和相位差φ可根据曲线得到。将运动和受力方程进行整理得到如下方程形式:

      F0sinωt+φ=-F0cosφAω2a¨+F0sinφAωa˙ ((1))

      基于自由运动方程:

      μa¨+λa˙+F=0 ((2))

      式中:

      μ——附加质量;

      λ——阻尼系数;

      F——浮体受力。

      因此,阻尼系数可表示为:

      -F0sin(Φ)Aω ((3))

      附加质量表示为:

      F0cos(Φ)Aω2 ((4))

      由上述过程计算得到风机基础的粘性阻尼(力矩)系数。

    • 模型的计算域划分及边界条件选取如图4所示。使用ICEM CFD进行网格划分。

      图  4  模型计算域及边界条件

      Figure 4.  Computational region and boundary conditions

      由于FLUENT计算时采用动方法,为避免结构物周围网格变形过大,将流场的网格划分为如下三个部分:

      1)结构物周围采用刚性网格,计算时随结构物一起运动,该部分选取尺寸较小的四面体网格进行划分。

      2)刚性网格外为变形区域,采用尺寸略大的四面体网格进行划分。

      3)为节约计算时间,其余区域采用非结构化大尺寸四面体网格,该区域距离物体较远,对结构的受力计算影响较小,且粗网格还可以起到人工阻尼的作用,减少侧壁面反射波影响。

      为保证计算精度,每个部分流场网格都在自由液面处对网格进行了加密。选取垂荡板圆盘直径为23.8 m,垂荡板高0.5 m,边锋延伸0.2 m的模型为例,图5为该模型计算域中纵截面的网格示意,图中模型网格总数为422万个。为了说明计算结果的正确性和可靠性,此处进行网格无关性验证,即验证计算结果对于网格密度变化的敏感性。通过改变网格的疏密,观察计算结果的变化,若其变化幅度在允许的范围之内,就可以说明网格误差在可接受精度范围内。因此,对模型重新划分网格,使得网格数量分别为212万个和821万个。

      图  5  计算域截面网格(422万网格)

      Figure 5.  Section mesh of the computational region

      各个网格模型的计算附加质量和阻尼系数结果如表3所示,计算结果对于网格密度变化的敏感性很小,网格无关性得以验证。

      表 3  不同网格数的计算结果

      Table 3.  Results with different mesh quantities

      网格数/万个附加质量阻尼系数/(N·s·m-1无因次阻尼系数变化率/%
      21210 236 877 9752 249 987 2620.114 4-2.54
      42210 036 272 2562 285 078 9000.117 40
      82110 148 250 0162 268 521 3990.115 9-1.29
    • 将发电工况下波浪的谱峰周期作为浮体强迫运动的周期,选取风机基础垂荡和横摇运动响应的有义值,由此规定简谐强迫运动形式如表4所示。

      表 4  测试工况参数

      Table 4.  Parameters in test cases

      运动形式发电工况
      周期/s幅值/m
      垂荡7.51.75
      横摇7.51.29

      将划分好的网格文件导入FLUENT商业软件,把上述强迫运动编写为用户自定义函数(UDF文件)来定义物体的简谐运动,定义简谐运动周期T=7.5 s,垂荡运动幅值A=1.75 m,横摇运动A=1.29°。计算采用Standard k-ε湍流模型,以及PISO求解格式。

    • 为了研究垂荡板边锋延伸宽度变化对浮式风机水动力阻尼系数的影响,将原模型及模型方案1、2、3的计算结果进行对比分析。各模型无因次阻尼系数计算结果见表5,边锋外伸量S变化时风机基础垂荡阻尼的变化趋势见图6,风机基础横摇阻尼的变化趋势见图7

      表 5  边锋延伸宽度对阻尼的影响

      Table 5.  Damping impact from the wing extension

      模型及运动无因次垂荡阻尼系数增长率/%无因次横摇阻尼系数增长率/%
      原模型0.175 300.103 30
      模型1(延伸0.2 m)0.176 50.680.106 32.95
      模型2(延伸0.4 m)0.183 74.790.106 02.61
      模型3(延伸0.6m)0.192 39.700.110 46.90

      图  6  边锋延伸宽度对基础垂荡阻尼的影响

      Figure 6.  Damping variation of heave motion from the wing extension of heave-plate

      图  7  边锋延伸宽度对基础横摇阻尼的影响

      Figure 7.  Damping variation of roll motion from the wing extension of heave-plate

      由图表数据可知,随着边锋外伸S的增大,基础垂荡粘性阻尼系数随之增大。由此可得,在一定范围内增大S,能增加垂荡板垂荡阻尼性能。当垂荡板边锋延伸S为0.6 m时,垂荡粘性阻尼系数增大近10%,阻尼性能增加效果较为明显。

      从图表数据来看,随着S的增大,基础横摇粘性阻尼系数总体也呈增大趋势。由此可得,在一定范围内增大S,能增加垂荡板横摇阻尼性能。但是相比于垂荡运动方向的阻尼变化值,横摇运动方向的增阻效果并不十分明显。

    • 为了研究垂荡板高度变化对浮式风机基础平台运动阻尼系数的影响,把原模型及模型方案1、4、5的计算结果进行分析。各模型无因次阻尼系数计算结果见表6,垂荡板高度H变化时风机基础垂荡阻尼的变化趋势见图8,基础横摇阻尼的变化趋势见图9

      表 6  垂荡板高度变化时风机基础阻尼的影响

      Table 6.  Damping impact from height of heave-plate

      模型及运动无因次垂荡阻尼系数增长率/%无因次横摇阻尼系数增长率/%
      原模型(高0.5 m)0.175 300.103 30
      模型1(高0.5 m)0.176 50.680.106 32.90
      模型4(高0.6 m)0.179 92.620.104 30.93
      模型5(高0.7 m)0.187 67.010.107 74.26

      图  8  垂荡板高度变化对基础垂荡阻尼的影响

      Figure 8.  Damping variation of heave motion from the height of heave-plate

      图  9  垂荡板高变化时基础横摇阻尼的影响

      Figure 9.  Damping variation of roll motion from the height of heave-plate

      由图表数据可知,随着垂荡板高度H的增大,基础垂荡方向的水动力粘性阻尼系数随之增大。由此可得,在一定范围内增大垂荡板高度,能增加垂荡板的垂荡阻尼性能。当垂荡板高度为0.7 m时,阻尼增大7%。

      从图表数据来看,随着垂荡板高度H的增大,粘性阻尼系数总体呈增大趋势。由此可得,在一定范围内增大垂荡板高度,能增加垂荡板横摇阻尼性能,但增长效果小于5%。

    • 为了研究垂荡板圆盘直径变化对垂荡阻尼系数的影响,把原模型及模型方案1、6、7的计算结果进行分析。各模型无因次阻尼系数计算结果见表7,圆盘直径D变化时风机基础垂荡阻尼的变化趋势见图10。基础横摇阻尼的变化趋势见图11

      表 7  圆盘直径变化时风机基础垂荡阻尼的变化

      Table 7.  Damping impact from the diameter of heave-plate

      模型及运动无因次垂荡阻尼系数增长率/%无因次横摇阻尼系数增长率/%
      原模型(直径22.8 m)0.175 300.103 30
      模型1(直径22.8 m)0.176 50.680.106 43.00
      模型6(直径23.8 m)0.184 75.360.117 413.64
      模型7(直径24.8 m)0.201 815.120.121 217.32

      图  10  垂荡板圆盘直径对基础垂荡阻尼的影响

      Figure 10.  Damping variation of heave motion from the diameter of heave-plate

      图  11  垂荡板圆盘直径对基础横摇阻尼的变化

      Figure 11.  Damping variation of roll motion from the diameter of heave-plate

      从图表数据来看,随着圆盘直径D的增大,浮式基础垂荡运动的粘性阻尼系数随之增大。由此可得,在一定范围内增大圆盘直径,能增加垂荡板的垂荡阻尼性能,并且随着圆盘直径增加呈非线性增长趋势,增阻效果明显,圆盘直径增加2 m(相比于原模型增加8.77%)时,垂荡阻尼增加可达15.12%。

      从图表数据来看,随着圆盘直径D的增大,浮式基础横摇运动的粘性阻尼系数增大。由此可得,在测试范围内增大圆盘直径,能增加垂荡板的横摇阻尼性能。圆盘直径增加2 m(相比于原模型增加8.77%)时,基础横摇方向的运动阻尼增加值可达17.32%,增阻效果较为显著。

      对比前述的结果,相对其他参数而言,增加垂荡板直径D,垂荡板阻尼增长效果最为显著,其次是增加边锋延伸长度,而增加垂荡板高度所带来的基础运动阻尼增加值十分有限。

    • 根据之前模型计算结果分析,调整垂荡板的尺寸,进一步提高垂荡板的阻尼性能。调整优化后垂荡板的尺寸为:直径24.8 m、高0.7 m、边锋延伸0.6 m的垂荡板。按照调整后垂荡板的尺寸,计算浮式基础的垂荡及横摇阻尼性能。将计算结果与原模型进行对比,可以看出优化后阻尼性能得到了理想的提升,具体阻尼系数增长情况如表8所示。优化后,该浮式风机平台的阻尼提升25%以上。

      表 8  优化后模型阻尼性能对比

      Table 8.  Damping comparison between the original model and optimized model

      模型及运动无因次阻尼系数阻尼系数增长率/%
      有量纲无量纲
      原模型(高0.5 m)0.175 3
      模型1(高0.5 m)0.201 426.17514.9
      模型4(高0.6 m)0.103 3
      模型5(高0.7 m)0.123 028.45519.1
    • 相关研究表明12,阻尼产生的根本原因是结构运动泄涡时对动能的耗散。对于风机基础结构来说,其在运动过程中产生、脱落的漩涡越多,能量耗散越大,阻尼性能也就越好。

      图12中,红色部分表示较大的正漩涡,黄色部分表示较小的正漩涡,深蓝色部分表示较大的负漩涡,浅蓝色部分表示较小的负漩涡。随着外伸边锋的增大,运动产生漩涡的范围会扩大,产生的漩涡量会增多,从而耗散更多能量,使得模型结构阻尼增加。

      图  12  浮式风机平台垂向运动涡量图

      Figure 12.  Vorticity contours of the foundation of the floating offshore wind turbine in heave

    • 根据上述阻尼结构尺寸敏感性分析和泄涡原理分析,为了得到阻尼性能更好的垂荡板形式,设计了一种“梅花形”的垂荡板。该垂荡板模型从优化后的模型变形而来,边锋宽度、板高及圆盘直径各参数都保持不变,在减小垂荡板圆盘面积的同时增加了垂荡板圆盘的周长。网格划分方法与之前模型相同,均采用“体网格+面网格”的方式。模型网格总数为459万个。垂荡板整体模型见图13所示,模型强迫运动时历曲线见图14所示,阻尼性能计算结果对比见表9

      图  13  “梅花形”整体模型示意图

      Figure 13.  Schematic diagram of the whole model of cinquofoil

      图  14  模型强迫垂荡运动时历曲线

      Figure 14.  Schematic diagram of the whole model of cinquefoil

      表 9  “梅花形”垂荡板阻尼性能计算结果

      Table 9.  Calculation results of damping performance of "quincunx-shaped" heave-plate

      模型及运动无因次阻尼系数阻尼系数增长率/%
      原模型垂荡0.175 3
      优化后模型垂荡0.201 414.90
      “梅花形”垂荡0.204 536.98

      在其他参数不变的情况下,“梅花形”垂荡板垂荡运动阻尼性能最佳,其垂荡阻尼系数较原始模型提升36.98%。

    • 本文选取10 MW半潜型浮式风机,利用计算流体力学的方法研究外立柱垂荡板结构几何尺寸变化对基础运动粘性阻尼的影响。数值计算结果表明,本文优化参数后的垂荡板能够明显地增加基础的附加质量与附加转动惯量,增加结构的阻尼,耗散运动能量,改善结构运动性能;带边锋的垂荡板的阻尼性能与圆盘直径、边锋延伸宽度、垂荡板高等参数均有关。一般来说,在一定范围内增大圆盘直径、边锋延伸宽度、垂荡板高,均会增加系统的附加质量和垂荡、横摇粘性阻尼系数。

      从总体情况来看,对阻尼性能影响最大的因素是垂荡板圆盘直径,尤其是对于横摇运动来说阻尼效果更明显。其次,增加边锋延伸宽度对增加垂荡阻尼性能也比较有利,而加大垂荡板的高度也可增大阻尼,但是增加缓慢。因此,增大阻尼三种措施按照有效性排序为:直径加大、边锋延伸、高度加大。本文优化方案中,直径24.8 m、高0.7 m、边锋延伸0.6 m的阻尼结构,该阻尼结构无量纲阻尼系数为:垂荡阻尼系数为20.14%,较原结构增大无量纲阻尼系数14.9%;摇摆阻尼系数12.30%,较原结构增大无量纲阻尼系数19.1%。本文设计的“梅花型”垂荡板结构相比原阻尼结构有效提升36.98%的垂荡阻尼。

      但是,在实际的设计过程中,垂荡板的设计尺寸和形式需要进一步结合考虑结构的强度、施工运输可行性以及经济性。

  • 陈嘉豪,尹梓炜,朱容宽等.漂浮式海上风机平台阻尼结构设计与研究[J].南方能源建设,2021,08(03):18-25.
  • 参考文献 (12)

    目录

      /

      返回文章
      返回