中等水深三立柱式浮式风机结构的风浪联合动力响应_王朝辉.pdf
《中等水深三立柱式浮式风机结构的风浪联合动力响应_王朝辉.pdf》由会员分享,可在线阅读,更多相关《中等水深三立柱式浮式风机结构的风浪联合动力响应_王朝辉.pdf(10页珍藏版)》请在咨信网上搜索。
1、SHIP ENGINEERING 船 舶 工 程 Vol.45 No.1 2023 总第 45 卷,2023 年第 1 期 24 中等水深三立柱式浮式风机结构的 风浪联合动力响应中等水深三立柱式浮式风机结构的 风浪联合动力响应 王朝辉,梁 峰(中国广核新能源控股有限公司,北京 100070)摘 要:摘 要:海上风资源丰富的风场大多处于水深大于 50 m 的海域,此时传统的固定式风机不再适用,因此,漂浮式风机的发展对于海上风电的开发利用有着重要的意义。然而,对于漂浮式风机系统在空气动力及水动力共同作用下的耦合运动响应特点仍需进一步研究。以 WindFloat 风机基础为例,基于势流理论以及叶素动
2、量理论,在考虑系泊缆作用下,求解漂浮式风机系统的时域耦合方程,模拟在中等水深下三立柱式风机结构的空气-水动力耦合响应,对漂浮式风机的耦合响应特点进行分析。研究表明:在工作工况下,受湍流风与风机偏航影响,横荡、横摇和艏摇具有明显的非零特性;纵摇响应的功率谱出现明显的多峰特征;在停机工况下,平台动力响应改由波浪荷载主控。关键词:关键词:中等水深;三立柱式;浮式风机结构;风浪耦合响应 中图分类号:中图分类号:P752 文献标志码:文献标志码:A 【DOI】10.13788/ki.cbgc.2023.01.04 Wind-Wave Coupled Response of Three Column Fl
3、oating Wind Turbine in Medium Water Depth WANG Zhaohui,LIANG Feng(CGN New Energy Holding Co.,Ltd.,Beijing 100070,China)Abstract:The offshore wind source is mainly located at the site with more than 50 m water depth where the fixed offshore wind turbine is not suitable.The floating offshore wind turb
4、ine(FOWT)is the solution which is a milestone on offshore renewable energy development.However,the characteristics of the FOWT response under the wind-wave coupled environment need further investigation.Taking the WindFloat floating wind turbine foundation as an example,based on the potential flow t
5、heory and the blade element momentum theory,the time-domain coupling equation of the floating fan system is solved under the consideration of the mooring line,and the air-hydrodynamic coupling response of the three-column floating wind turbine structure under the medium water depth is simulated,and
6、the coupling response characteristics of the floating wind turbine structure are analyzed.The results show that:under operational conditions,sway,roll and yaw have obvious non-zero characteristics due to the influence of turbulent wind and the yaw operation of the turbine;the power spectrum of the p
7、itch response has obvious multi-peak characteristics;under parked condition,the platform dynamic response is controlled by wave load.Key words:medium water depth;three columns;floating wind turbine;wind-wave coupled response 0 引言引言 近年来,全球对清洁能源的需求不断增加,海上风电这一项丰富的可再生能源正在各国迅速发展。海上风机主要分为固定式和漂浮式 2 种。在浅水区域
8、,主要以单桩和导管架式风机基础为主。随着水深的加深,固定式风机的建造成本大幅上升,因此采用漂浮式风机。漂浮式风机可摆脱海床及水深的束缚,同时其漂浮式基础的随浪运动特点使其具备大大降低在海 收稿日期:2022-01-25;修回日期:2022-03-21 作者简介:王朝辉(1967),男,高级工程师。研究方向:新能源规划设计、海上风电规划及接入系统、新型电力系统接网技术,以及海上风电结构基础设计。王朝辉等,中等水深三立柱式浮式风机结构的风浪联合动力响应 25 上风浪下的结构破坏可能性。漂浮式风机的概念首先由 HERONEMUS1在1972 年提出,但由于当时海洋工程技术的限制,其构想未能实现。然而
9、,随着海上油气发展带动的海洋工程技术进展以及全球对于清洁能源的需求,漂浮式风机相关研究在近年来取得了巨大进展。挪威石油公司与挪威科技大学合作开发了世界第一台漂浮式风机样机 Hywind2,NIELSEN 等3对 Hywind 进行了 140的模型试验以及数值分析。BROMMUNDT 等4对半潜式漂浮式基础 Tri-floater 的系泊系统进行了分析。BACHYNSIKI 等5对张力腿式、半潜式以及单筒式漂浮基础风机在故障下的动力特性响应做了研究。美国国家能源实验室(National Renewable Energy Laboratory,NREL)开展了对漂浮式风机基础DeepCWind 的
10、一系列研究(OC3-4-5-6),其中包含了势流数值模拟,模型试验以及计算流体力学验证6。NAVA 等7对单筒式以及半潜式漂浮式风机基础在典型风浪作用下进行了长期预报。赵永生等8对张力腿平台在极端海况情况下的动力响应特征进行了研究,对张力腿平台的极限承载能力进行了研究。唐友刚等9设计了一种新型单筒式漂浮风机基础,并针对该风机基础进行了运动与受力特性的研究。张亮等10针对三立柱式海上风机半潜台进行了对考虑风载荷情况下的整体稳性与破舱稳性进行了分析。任年鑫等11设计了一种新型张力腿及锚缆组合式系泊系统并将其应用于多种漂浮式风机结构。ZHENG 等12对提出了三峰谱分析方法对漂浮式风机结构疲劳进行分
11、析。REN 等13对多立柱间的浮体水动力耦合效应进行了分析,认为应充分考虑立柱间距对浮体间波面高程的影响。LIU 等14采用OpenFOAM 编写了针对 OC4 的风浪联合数值计算模型,并获得了良好仿真效果。为了进一步揭示漂浮式风机在风浪耦合作用下的运动响应特点,本文以 WindFloat 风机基础为例,对三立柱半潜式平台在中等水深情况下的系泊系统进行设计及优化并对其耦合运动响应特点进行分析。1 计算原理计算原理 1.1 理想流体运动方程理想流体运动方程 流体在流场中运动时需遵循质量守恒(连续方程)和动量守恒(运动方程)。在研究海水波动时又通常假定海水是无黏性不可压缩流体。波动海域中小尺度直立
12、柱波浪力的计算坐标系见图 1,海水的连续方程和运动方程可分别写为 0yxzuuuxyz+=(1)d1()duuuupFtt=+=-+(2)式(1)和式(2)中:ux、uy和 uz为水质点速度 u 在3个坐标轴方向上的分量;t 为时间;为海水密度;p为压强;(1/)p 和F分别为单位质量流体所受到的压强梯度力和重力。无黏性不可压缩流体的运动是无旋的,因此存在单值速度势(x,y,z,t),满足:(,)ux y z t=(3)图1 小尺度直立柱波浪力计算坐标系 则连续方程和运动方程可分别写为 2222220 xyz+=(4)1()()02pgzt+=(5)式(5)中:g为重力加速度。式(4)即是拉普
13、拉斯方程,一旦得到速度势(x,y,z,t),便可根据式(3)和式(5)求得波动场各点的速度u和压强p,进而得到整个波动场的波面形状。对于波动场速度势的求解需满足以下边界条件15:n00012(,)(,)(,)(,)zdzdzxyzzzzzttunuuutxygtx y tx yx y z tv x y z=-=-=|=|=+|+|=|=海底运动边界条件:自由表面运动边界条件:自由表面动力边界条件:()=0 无穷远处边界条件:(6)式中:d为水深;n表为海底的法向;un为水质点速度的法向分量;为自由水面的铅直位移,也即波面方程,是关于(x,y,t)的函数;(x,y)为初始水面起伏条件;v(x,y
14、,z)为初始速度场。1.2 莫里森方程莫里森方程 莫里森方程是MORISON等于1950年提出的计算小尺度结构物受到的波浪力的经验公式16。假定结构物的存在对波浪运动无明显影响,结构物受到的波专题:海洋可再生能源 26 浪力可分解为拖曳力和惯性力。如式(7)所示,单位长度直立圆柱受到的水平波浪力可写为()2Hdm2dMd1d1+24dd124dxxxxxxuDFCDu uCtuDCDu uCt=+=+(7)式中:ux和dux/dt分别为柱体中轴线高度z处水质点的水平速度和加速度;D为圆柱体的直径;Cd为垂直于圆柱体轴线方向的拖曳力系数;Cm为附加质量系数;CM为质量系数。式(7)中的第1项即为
15、拖曳力项,第2项为惯性力项。沿着柱体积分即可得柱体受到的波浪力。莫里森方程是半经验公式,其精度取决于Cm和Cd。1.3 动量动量-叶素理论叶素理论 根据叶素理论和动量理论17均可得到叶素环受到的推力和扭矩。利用2种理论下的推力和扭矩分别相等(dT叶素=dT动量,dM叶素=dM动量),即可求出轴向干扰系数b和周向干扰系数b。()ld2cossin18sinNC CCbb+=-(8)()ldsincos18sincosNC CCbb-=+(9)式(8)和式(9)中:N为风机叶片的数目;C为叶素的弦长;Cl和Cd分别为叶片的升、阻力系数;为相对风向角。得到轴向和周向干扰系数后,即可得到风轮受到的推力
16、和扭矩。1.4 塔筒风荷载塔筒风荷载 风机塔筒是1个长度近百米的高耸构件,塔筒遭受的风荷载不可忽略。单位长度塔筒上的风荷载可按式(10)计算。2wairdtower1dd2FC DUz=(10)式中:air为空气密度,取1.28 kg/m3;Cd为塔筒的气动拖曳力系数,取1;Dtower为塔筒直径;U为风速。1.5 时域耦合分析方法时域耦合分析方法 在浮式系统耦合分析过程中,风机/浮体及其系泊系统的动力响应同时求解,可以得到所有的动力相互作用。具体来说,系泊系统产生的作用力、阻尼和惯性力等的耦合效应包括6个部分:1)系泊系统产生的静态回复力,该力是浮体位移的函数;2)流载荷及其 对系泊系统回复
17、力的影响;3)系泊缆躺底段与海床的接触载荷;4)系泊缆和立管的运动(特别是横向运动)以及流载荷的作用等产生的阻尼;5)浮体与柔性构件相互接触而产生的阻尼;6)系泊缆和立管的存在而产生的附加惯性力。风机对浮体产生的耦合效应主要包括:风载荷对风机叶片及塔筒的作用力产生的力/弯矩在浮体上的作用;浮体运动导致风机机头迎风角度变化对漂浮式风机上部结构风载荷的影响。ENGSETH等18提出了浮体/系泊系统的整体耦合动力分析方法,空间离散系统的控制动态平衡方程为 IDSE(,)(,)(,)(,)R r r tRr r tRr tRr r t+=?(11)式中:RI、RD、RS和RE分别为惯性力、阻尼力、内部
18、相互作用力和外部激励力;r、r?和r?分别为结构运动的位移、速度和加速度矢量。惯性力矢量和阻尼力矢量的形式为 I(,)()R r r tr r=M?(12)?D(,)()Rr r tC r r=?(13)式(12)和式(13)中:M?为系统的质量矩阵,包括结构质量、立管内部流体对整体质量的影响和水动力质量;C?为系统的阻尼,包括结构内部阻尼、水动力阻尼和特定离散阻尼。内部相互作用力基于结构单元的瞬时应力状态得到,外力包括重力、浮力、平台运动产生的力、环境载荷和指定的节点集中力。式(11)为非线性系统的微分方程,非线性可能是由于惯性力和阻尼力中的位移相关性,以及由于外部载荷矢量与结构位移和速度之
19、间的耦合导致。此外,惯性力与柔性结构变形之间可能存在非线性关系。在计算求解时,浮体所受到的波浪激励力及势流阻尼由频域水动力分析得到的参数经过频转时的变换得到,而细长构件所受的波浪荷载及其与浮体之间的相互作用则根据二者实时的位置状态进行计算,其耦合状态可采用Newmark-方法和牛顿迭代法进行时域内的迭代求解。2 数值模型数值模型 2.1 模型主尺度及网格划分模型主尺度及网格划分 按照WindFloat的基本参数(见表1)进行数值模型的建立。本文中采用边界元法对漂浮式风机的水动力参数进行计算,首先建立平台基础结构模型,然后对基础模型进行边界元网格的划分。表 1 WindFloat-10 MW 风
20、机参数表 参数 数值 参数 数值 参数 数值 参数 数值立柱直径/m 10.7 平台质心位置/m(32.56,0,11.07)垂荡板边长/m15.629壁厚/m 0.03立柱间距/m 56.4 平台浮心位置/m(32.56,0,12.13)垂荡板厚度/m0.04 风机塔筒高度/m 125 立柱高度/m 33.6 压舱水质量/kg 3.80106 横撑直径/m 2.2 叶片长度/m 86 平台质量/kg 2.06106 吃水/m 22.9 斜撑直径/m 2.0 机头质量/t 1 200 王朝辉等,中等水深三立柱式浮式风机结构的风浪联合动力响应 27 对WindFloat平台进行内部舱室建模,将平
21、台三立柱引入舱室结构(单个立柱为1个舱室),压舱为水压载,从而计算得到不同舱室的压舱水结构,进而求得压载水的质量、质心以及转动惯量等关键参数。随后对整体结构模型进行网格剖分,生成对应的边界元模型,对应浮式平台的面元模型、结构模型和水动力模型等见图2图4。图 2 WindFloat 面元模型 图 3 WindFloat 结构模型 图4 WindFloat 水动力模型 2.2 水动力参数分析与模型验证水动力参数分析与模型验证 开展WindFloat半潜式浮式风机的频域水动力性能分析,浪向计算角度范围0360,间隔15,波浪计算频率范围05 rad/s,间隔0.05 rad/s,进而得到WindFl
22、oat半潜式浮式风机频域水动力数值结果,见图5。由图5可知:1)相同波浪入射角下:在纵荡和横荡方向,RAO随波浪入射频率的增加而逐渐减小;在垂荡、横摇和纵摇方向,RAO峰值最大处对应的频率分别为0.30 rad/s、0.25 rad/s和0.25 rad/s,分别对应垂荡、横摇和纵摇的固有频率,同时,RAO在垂荡、横摇和纵摇方向在0.4 rad/s1.0 rad/s正常波浪频率范围内出现较大,表明波浪对垂荡、横摇和纵摇方向的运动具有较大影响,垂荡和纵摇的运动幅值最大分别为2.5 m/m和1.5()/m左右;艏摇的RAO较小,峰值集中在0.8 rad/s附近,最大值约为0.5()/m。图 5 W
23、indFloat 六自由度运动幅值响应算子(RAO)2)在不同波浪入射角度下:0浪向时,纵荡、垂荡和纵摇运动较为明显,最大值分别为4.50 m/m、2.50 m/m和2.25()/m;垂荡运动对于波浪入射角的变化不敏感,在不同入射角的波浪作用下,变化趋势基本一致,RAO运动幅值较为接近;纵摇方向在不同波浪入射角的作用下运动幅值有较为显著的变化,0浪(a)纵荡(b)横荡(c)垂荡(d)横摇(e)纵摇(f)艏摇 专题:海洋可再生能源 28 向角入射时,响应最大。3 数值模拟分析数值模拟分析 3.1 计算工况计算工况 为了考虑不同工况下的浮式风机动力响应特性,开展风浪联合作用下的浮式海上风机耦合动力
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 中等 水深 立柱 式浮式 风机 结构 风浪 联合 动力 响应 王朝
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。