给我们留言
91学术服务平台

您好,欢迎来到91学术官网!业务合作:91xueshu@sina.com,站长邮箱:91xszz@sina.com

发布论文

论文咨询

海上风电机组三脚架支撑结构疲劳分析方法

  2024-08-07    48  上传者:管理员

摘要:为实现海上风电机组支撑结构抗疲劳设计,提出基于整体-局部方式建模的疲劳分析方法。利用该方法对热点局部区域进行详细有限元建模,基于IIW规范采用高精度六面体实体单元离散,其余非重点区域采用梁单元离散,可克服以往运用应力影响矩阵研究中忽略整体结构位移、边界条件不易施加等缺点。基于IEC 61400规范进行整机动力学和载荷计算,对某海上风电机组三脚架支撑结构进行疲劳分析。对比分析不同载荷工况、波高和水深下的累积损伤值,结果验证了所提出的整体-局部建模方法在海上风电机组三脚架支撑结构疲劳分析中的可行性和优越性。

  • 关键词:
  • 抗疲劳性
  • 海上风电机组
  • 疲劳损伤
  • 疲劳载荷
  • 网格划分
  • 加入收藏

近年来陆上风电机组倒塔事件频发,结构强度问题受到广泛关注[1-2]。由于海上风电机组的设计、安装与运维成本更高,支撑结构强度尤其是疲劳强度成为研究热点。应力集中因子(stress concentration factor,SCF)概念简单,在支撑结构设计中应用广泛[3],其最初研究集中在较简单的K型、X型管节点上,后来发展至复杂结构形式、复合材料加强的管节点上[4]。在疲劳预测方面,董文斌等[5]提出采用两参数威布尔函数和广义γ函数描述典型工况导管架应力变程分布;Yeter等[6]提出基于时频转换三脚架的频域疲劳计算方法;徐传恒等[7]和Larsen等[8]分别提出K型管节点的多轴疲劳预测方法和疲劳分析的随机有限元法;刘伯紘等[9]在应力影响矩阵中增加了重力项以得到更精确的时序应力数值;Shittu等[10]对比了疲劳理论和断裂力学评估海上风电基础强度结果的异同;朱航等[11]探究三筒导管架基础的浮运稳性以及极限强度和疲劳强度。有关风电机组支撑结构抗疲劳研究可参考文献[12]。

三脚架是海工结构中常见的支撑型式,可用于超过30 m水深且单机容量大于5 MW的风电机组。王滨等[13]分别考察了单方向载荷与全方位载荷对三脚架基础疲劳寿命的影响;马宏旺等[14]提出考虑局部和全局冲刷效应的冲刷模型;梁发云等[15]考察了冲刷作用下的三脚架桩基的拉力作用;袁春光等[16]开展了深水域三脚架冲刷实验;张力伟等[17]建立考虑桩土和流固耦合作用的三脚架支撑风电机组动力学模型;Arcigni等[18]考察了海生物对三脚架动态性能的影响;郝二通等[19]对比了不同基础支撑下抗船体碰撞能力;Luczak等[20]对比了两种基于振动信号的三脚架损伤检测方法;康海贵等[21]验证了粒子群算法在三脚架基础结构优化设计中的有效性;杨和振等[22]提出动力响应下的三脚架可靠性优化设计方法;李炜等[23]对比了固定式海上风力发电机组不同组成部分对基频的影响;鲍石榴等[24]开展了由面外弯矩引起的三脚架热点应力数值仿真与实验研究。

在三脚架支撑结构中,主筒和3个斜支撑相贯处的外形复杂,且受到面内外弯矩、推力等复杂载荷作用。目前应用较为成熟的SCF方法仍存在诸多不足,诸如:通常选取极限强度的薄弱点或根据经验来假设疲劳失效点,无法得到全场累积疲劳损伤分布;基于单一方向载荷作用的SCF评估方法忽略了多方向载荷的叠加效应;在风浪流耦合作用下,有关波浪载荷等对结构强度的影响缺乏定量研究。鉴于上述问题,本文提出整体-局部的三脚架有限元建模方法,克服以往应力影响矩阵方法无法考虑整体结构边界条件的不足,考察关键部位疲劳损伤分布特点,证明所提方法在三脚架疲劳分析中的必要性,考察波高、水深环境参数对疲劳损伤的影响规律。


1、基于应力影响矩阵三脚架疲劳分析


1.1应力影响矩阵与应力计算

假设外界载荷通道数为n,各分量方向对应的单位大小载荷记为f1,f2,…,fn,则应力影响矩阵(stress influence matrix,SIM)定义为:

式中:第j列代表由fj引起的各应力分量。值得注意的是,在实际应用中,通过n次有限元计算可一次性获得SIM。

在外载荷作用下,时序应力数值基于线性叠加原理计算,某点应力分量为:

式中:[F1F2F3F4⋯Fn]T——对应f1,f2,…,fn方向上的载荷大小。

由式(2)可知,在获取SIM矩阵的基础上,通过线性叠加原理,可避免外载荷庞大数据下的有限元计算,大幅提高疲劳分析的效率。由式(2)可计算得等效应力表达式为:

1.2焊缝S-N曲线

根据DNV-RP-C203标准确定关键部位焊缝处的S-N曲

式中:N——循环次数;lg a——S-N曲线在横轴截距;m——S-N曲线斜率的负倒数,这里选取阴极保护海水环境,在循环次数1×106前后,m取值分别为3和5;t——预计裂纹产生处的结构厚度,mm;k——指数参数,这里取值0.2。

1.3雨流计数法与累积疲劳损伤计算

常见疲劳分析计数包括穿级类型计数法、峰值计数法、幅值计数法等,这里采用普遍应用的雨流计数法对时序应力统计处理[25]。

由式(3)可知,等效应力σSVM恒为正值,因此无法反映拉压特性,这里采用带符号的等效应力作为雨流计数对象:

式中:W——表征主应力方向的矩阵。

设在第i个应力循环实际次数为ϖi,对应S-N曲线中该应力变程允许次数为,基于Miner法则判断疲劳破坏发生条件:


2、三脚架结构疲劳强度分析


2.1三脚架结构的载荷计算

某5 MW海上风电机组轮毂高度为80 m,风轮直径为118 m,基于DNV BladedTM建立整机动力学模型,并采用隐式Newmark积分求解。设参考水深42 m,1倍浪高为1.1 m(对应4 m/s风速)。根据IEC 61400-3标准选取设计载荷工况中的正常发电工况[26]。风速范围4~26 m/s,步长2 m/s,风切变指数0.14,偏航误差取-10°、0°和10°,共包含216个疲劳工况,各工况在20 a发生次数根据风速威布尔分布曲线计算得到。动力学方程求解得到整机各点的时序位移值,由此计算等效静力载荷,即在该静力载荷下产生与动力响应相同的位移。将等效静力载荷中的载荷分量施加到详细有限元模型上。每一疲劳工况共含12000个时刻点对应的载荷数据量,若进行有限元求解,将面临着庞大的计算量。以下将求解3个方向单位力和单位力矩作用下的有限元方程,形成SIM,利用线性叠加原理获取时序应力。

2.2基于整体-局部建模的三脚架疲劳分析

在应用SIM时,以往研究将整体分析的载荷或位移施加到局部模型上作为边界条件,这种处理方式忽略了整体与局部关联性,且当施加载荷或位移的通道数量较多时,不易施加边界条件。为克服上述缺陷,提出整体-局部建模方法,选取三脚架主筒与斜支撑相贯处作为重点考察部位。如图1所示,焊缝处由被焊接件和焊趾组成,根据IIW标准这里不再对焊趾结构建模,其力学性能在焊缝S-N曲线中体现。对焊缝及周边区域采用高精度六面体单元离散,非重点考察区域采用一维梁单元模拟,网格模型包含58.8万个单元和89.4万个节点。由此可知,该处理方式在获取局部结构精确应力值的同时,可避免整体离散对应的庞大有限元规模,同时弱化了施加在局部结构边界条件对整体分析结果的依赖性。

图1两种不同三脚架疲劳分析方法对比

全约束三脚架底部三节点,在顶部梁单元节点上施加6个通道载荷,即Fx、Fy、Fz、Mx、My和Mz,累积损伤分布如图2所示。由图2可知,最大累积疲劳损伤值为0.6069,寿命设计满足疲劳强度要求。在SCF方法中,通常选取如图2所示的8个均布点作为考察热点,但图2中的最大疲劳损伤点位于斜支撑与主筒相贯处的某点,并非预设热点。其次,相对于SCF方法逐点计算累积损伤值,所提方法在SIM基础上,快速计算全场累积损伤值,更加全面了解疲劳损伤分布情况。由此可证明所提方法在三脚架结构疲劳分析的可行性和优越性。

图2三脚架累积疲劳分布与常见热点分布

2.3不同载荷作用的累积损伤结果对比

本节考察不同载荷作用下的累积损伤值,以证明所提方法的必要性。假设Fx、Mx、My和Mz单独作用,得到如图3所示的累积损伤结果。由图3可知,4种情况下的累积损伤值分别为8.07×10-5、3.48×10-4、0.121和0.02293,均低于0.6069。结果表明My是引起疲劳破坏的主要因素,但其他载荷贡献也不可忽略,所以采用单一载荷作用结合SCF方法将低估累积损伤值,得到此较危险的设计。上述结果证明了所提方法在三脚架结构疲劳强度校核中的必要性。

图3不同载荷作用工况下累积损伤结果对比

2.4不同波高下的累积损伤结果对比

在以往研究中,较少考察风浪流耦合作用对累积损伤的影响。这里设置1~5倍的波高,不同波高下累积损伤结果对比,如图4所示。由图4可知,累积损伤最大值随波高倍数的增加而增加,在当超过4倍波高时急剧增加,这说明随着波高的加高,浪流载荷对疲劳损伤的影响作用加强。

图4不同波高下的累积损伤结果对比

2.5不同水深下的累积损伤结果对比

保持波高不变,设置水深变化范围为45~60 m,不同水深下的累积损伤结果对比如图5所示。

由图5可知,在一定水深范围内,最大累积损伤值随水深的增加而增加;水深超过50 m时的累积损伤值超过1,不再满足疲劳强度设计要求。

图5不同水深下的累积损伤结果对比


3、结论


为实现海上风电机组三脚架支撑结构的抗疲劳设计,提出基于整体-局部方式的建模方法,针对某海上风电机组,计算得到所选局部结构的累积损伤分布,考察了不同方向载荷引起的疲劳损伤结果,对比了不同波高、水深下的累积损伤数值,可得到如下主要结论:

1)所提出的整体-局部有限元建模方法能保持原有SIM方法精细建模、计算精度高的优点,同时克服了其忽略整体结构响应、边界条件不宜施加的缺点。相比较较为成熟的SCF法,可快速获取全场累积损伤分布。

2)三脚架结构累积损伤是不同方向载荷作用叠加的结果,单一载荷方向结合SCF方法会低估累积损伤值,得到比较危险的设计结构。

3)最大累积损伤值随着波高、水深数值的增加而增加,当波高或水深增加到一定数值时,累积损伤值会超过许用值,浪和流引起的载荷在一定程度上起到主导作用。


参考文献:

[2]龙凯.风电机组塔筒结构强度校核与优化设计[M].北京:中国水利水电出版社,2020.

[11]朱航,李孟超,郭耀华,等.深水海上风电三筒导管架基础受力及浮运分析[J].太阳能学报,2022, 43(11):269-276.

[13]王滨,李昕,明小燕.海上风力机基础结构疲劳损伤研究[J].太阳能学报,2015, 36(3):769-774.

[21]康海贵,李玉刚,郇彩云.基于可靠度的海上风机基础结构优化设计方法研究[J].太阳能学报,2009, 30(12):1602-1607

[23]李炜,张敏,刘振亚,等.三脚架式海上风电基础结构基频敏感性研究[J].太阳能学报,2015, 36(1):90-95.


基金资助:南方电网新能源联合实验室(GDXNY2024KF03);中国华能集团科技项目“超大型国产化海上风电机组关键技术研究(一期)”(HNKJ23-H13);华能集团海上风电与智慧能源系统科技专项(HNKJ20-H88-01);新能源电力系统全国重点实验室开放课题(LAPS23015);


文章来源:陆飞宇,李成功,龙凯,等.海上风电机组三脚架支撑结构疲劳分析方法[J].太阳能学报,2024,45(07):699-703.

分享:

91学术论文范文

相关论文

推荐期刊

网友评论

我要评论

太阳能

期刊名称:太阳能

期刊人气:1861

期刊详情

主管单位:中国科学技术协会

主办单位:中国可再生能源学会

出版地方:北京

专业分类:科技

国际刊号:1003-0417

国内刊号:11-1660/TK

邮发代号:2-164

创刊时间:1980年

发行周期:月刊

期刊开本:大16开

见刊时间:7-9个月

论文导航

查看更多

相关期刊

热门论文

【91学术】(www.91xueshu.com)属于综合性学术交流平台,信息来自源互联网共享,如有版权协议请告知删除,ICP备案:冀ICP备19018493号

400-069-1609

微信咨询

返回顶部

发布论文

上传文件

发布论文

上传文件

发布论文

您的论文已提交,我们会尽快联系您,请耐心等待!

知 道 了

登录

点击换一张
点击换一张
已经有账号?立即登录
已经有账号?立即登录

找回密码

找回密码

你的密码已发送到您的邮箱,请查看!

确 定