海冰作用下近海桩柱式风力机动力学响应研究

叶柯华+李春+叶舟+杨阳+缪维跑

摘 要： 高纬度低温海域海平面存在大量运动状态的海冰，位于此处的近海樁柱式风力机容易受到不规则的海冰载荷作用，风力机平台、塔架和叶片等结构部件的动力学响应均受其影响.为定性及定量分析海冰载荷对叶片和塔架的结构动力学响应的影响程度，以NREL 5 MW 近海4桩柱式风力机为研究对象，耦合风载荷、波浪载荷及海冰载荷，通过Kane方法建立风力机动力学模型，其中海冰载荷通过冰力函数定义.对比分析了在IEC Lock in冰力函数、Mttnen海冰模型和无海冰作用三种工况下叶片和塔架的结构动力学响应，结果表明：海冰载荷使塔顶位移增加，在Mttnen海冰模型作用下塔顶位移增加了24.1%，在IEC Lock in模型作用下则增加了16%；两种海冰模型均不同程度地使叶片挥舞振动的频率增大，其中在Mttnen海冰模型中变化更加剧烈，这极大地增加了叶片的疲劳载荷.

关键词： 桩柱式风力机； 海冰载荷； 冰力函数； 动力学响应

中图分类号： TP 392 文献标志码： A

Abstract： There exists a lot of floating ice at low temperature and high latitude sea.The dynamic response of offshore piling wind turbines platform，tower and blades is affected by the irregular ice load.In order to analyze the effect of ice load on the dynamic response of the blade and tower qualitatively and quantitatively，the Kane method was applied to build the dynamic model based on NREL 5 MW offshore 4piling wind turbine with the loads of wind，wave and ice.The ice load was defined by the ice force function.The comparative analysis was performed for the blade and tower dynamic response among IEC Lock in the ice force function，Mttnen model and no ice load.The results showed that the displacement of tower top increased significantly with the ice load.It increased by 24.1% and 16% for Mttnen model and IEC Lock in model，respectively.Both models increased the waved vibration frequency of the blade with different levels.In particular，Mttnen model increased the fatigue loading of the blade.

Keywords： piling wind turbine； ice load； ice force function； dynamic response

随着石化能源逐渐枯竭，人们对环境保护要求日益提高，风能因其技术成熟、储量丰富和汲取方便等优点，受到各国重视[1].由于海上风能具有能量密度高、风速稳定和蕴藏量巨大等优点，发展海上风电已成为近年来的研究热点[2].太阳能辐射至地球的能量的2%转化为风能，约为2.74万亿kW[3].截至2013年，全球海上装机容量已达6 837 MW.我国“十二五”规划提出，2015年投入运行海上风电装机容量5 000 MW，2020年将达到30 000 MW[1].风能的发展，一方面满足了人们对可再生能源的需求，另一方面也迫使风力机面对更多样、更复杂的运行环境.近海桩柱式风力机主要受到大湍动度的海风载荷、不规则的海浪载荷、海流载荷及高纬度低温海域存在的海冰载荷作用[4]，分析风力机结构在如此复杂的运行工况下的动态响应对风力机安全运行显得愈加重要.

目前，对近海桩柱式风力机的研究主要考虑风、波、流单激励作用或多激励耦合作用，而忽略了海冰对风力机平台的动静冰力作用[5-7].作用在垂直圆柱桩柱上的冰力可达5 MN，即使采用锥形抗冰结构，冰力仍可达800 kN[8].海冰对风力机动态响应影响显著，尤其体现在塔顶和叶片位移与变形[9].

历史事实和现有研究均表明海冰足以对海工平台造成破坏甚至使其报废.1969年海冰造成渤海海域“海井一号”石油平台倾覆，“海井二号”受到破坏[10]；1986年，库页岛附近的Molikpaq钻井平台因海冰作用导致的“频锁”振动，使砂土严重液化，平台损毁[11]；文献[12]研究表明，当海冰冰激振动频率和海工平台固有频率接近时，海工平台发生“频锁”振动，平台振幅和应力急剧增加.以上研究主要针对高宽比较小的海工平台，而风力机为大高宽比悬臂梁结构，结构稳定性更低，对海冰作用更敏感，受海冰作用影响更明显，所以研究海面海冰对风力机时域动态响应的影响十分必要.

本文以NREL 5 MW近海4桩柱式风力机为仿真模型，通过Kane方法建立风力机动力学模型，使用IEC Kaimal风谱模型生成相应的时域风场以计算风载荷，选用JONSWAP波浪谱生成的波浪载荷，分别采用IEC Lock in海冰模型[13]和Mttnen海冰模型[14]，对比分析两种海冰模型下和无海冰作用下风力机结构部件的时域动态响应.endprint

1 Kane动力学模型

模拟风力机整机动力学响应是多体结构耦合的复杂过程，包括气动弹性、水动弹性和风轮塔架结构耦合等.叶片、塔架、低速轴等为柔体，机舱、轮毂、高速轴等为刚体[15].本文以Kane方法建立动力模型，在保证计算精度的前提下，避免采用有限元计算产生的巨大计算量[16].

求解时，在每个时间步中，方程的数值解的第一步采用四阶AdamsBeshforth预测修正算法的预测方法确定低阶项的值，并以此构成方程的右边项，然后采用Gauss消元法求解系统自由度的加速度，计算得到的加速度用于修正预测值，并给出该时间步的最终解.由于该预测修正算法不是自发的，前四个时间步的解需用四阶RungeKutta法确定.

2 海冰作用

海冰运动由地球动量系统和能量系统决定，实际海冰运动则受风、波、流等驱动.由于驱动力的不确定性及非连续性，海冰作用力（简称冰力）属于随机载荷.当冰力变化较小时，才可视冰力为定值，载荷作用方式为静冰力作用，其他为动冰力[19].

静冰力对结构的作用主要包括挤压破坏、压屈破坏、纵向剪切破坏、弯曲破坏、摩擦破坏和混合破坏等[19].动冰力作用不仅造成上述破坏，还将引起风力机结构振动.由动冰力造成的结构振动称为冰激振动.目前对于冰激振动形成的原因学界存在两种不同观点：强迫振动和自激振动[20-21].完全取决于冰本身条件和环境因素，与海工结构自身特点无关的振动为强迫振动；而自激振动不仅取决与冰本身条件和环境因素，还受到海工结构本身特性影响.

当海冰对海工平台间歇性的碰撞频率与结构固有频率接近时，结构振动将迫使海冰碰撞频率保持在结构固有频率附近，此类海冰结构碰撞称为“锁定”碰撞[22].

2.3 Mttnen冰力函数

Mttnen模型为典型的自激振动模型.该模型解释冰激振动的依据是冰体连续破坏假设和海冰抗压强度与应力速率关系曲线.海冰抗压强度与应力速率关系曲线中，抗压强度随加载速率的增加先增大后减小，在曲线的下降段由于破坏应力随加载速率增大而降低，因此产生负阻尼效应.当负阻尼大于结构阻尼时会导致结构振动失稳，振幅增大，而由于阻尼非線性，振幅将趋于在一个较大且稳定的范围[24].

3 算例模型

3.1 风力机模型

本文以美国可再生能源实验室（National Renewable Energy Laboratory，NREL）的5 WM近海4桩柱式风力机模型作为研究对象.塔架底端连接着桁架，桁架将固定在海床的桩柱与风力机连接为一体.四个桩柱以塔架为中心对称排列，桩柱之间关于海冰载荷的相互影响通过桩遮蔽系数表示.桁架桩柱直接受到海冰冰力作用，它是保证风力机支撑稳定的重要部分.

为方便计算，定义多个坐标系，分别为叶片坐标系{O，X，Y，Z}、轮毂坐标系{O′，X′，Y′，Z′}、塔架坐标系{O″，X″，Y″，Z″}和地基坐标系{O，X，Y，Z}.风力机坐标示意图如图2所示.风力机模型和桩柱参数分别如表1、2所示.

3.2 波浪载荷

波浪在时域上呈现不规则性，可用无数个简单正弦波叠加表示，进行频谱变换得到连续频域分布，进一步利用有限个离散频率近似表示.波浪谱是随机海浪统计性质，表示各波浪成分不同频率中能量分布.常用波浪谱有PM谱、JONSWAP谱、Neumann谱和文氏谱等.本文采用JONSWAP谱作为波浪谱模型.谱模型为[25]

3.4 风载荷

风既是风力机汲取风能的源泉，又是最主要的环境载荷.风载荷主要作用于叶轮和塔架，不仅作用力巨大，而且大小随时间波动，沿空间变化.风载荷在利用统计学原理表达时，其受到脉动风速谱和风速廓线模型影响[26].

脉动风速谱是湍流中各种不同频率下漩涡对湍动能贡献大小的表征.Kaimal谱考虑湍流功率随高度变化，能较好地模拟近海桩柱风力机风场环境[27].

风的三个分量（K=u，v，w）中，u为来流方向风速；v为侧向风速，w为垂直方向风速.计算式为

通过风谱得到轮毂中心处时域风速，采用幂律风廓线描述垂向风速，对数风廓线定义纵向风速变化.轮毂中心处三个方向风速分量时域大小如图3所示.

4 结果分析

近海桩柱式风力是大型复杂刚柔多体结构，主要的柔性部件有叶片和塔架.叶片作为风力机唯一汲取风能的部件，其安全性尤为重要.塔架为风力机主要支撑结构，在湍流风场、不规则波浪和海冰载荷作用下，塔架承受巨大载荷，容易出现剧烈摆动.首先，对比分析了塔顶前后位移的动态响应，结果如图4所示.

由图4可知，在三种工况下塔顶前后位移均在风力机刚启动时出现巨大波动.随着时间推进，在无海冰载荷工况下，其位移逐渐趋于稳定，在30 s后平均值为0.205 m.而其他两种海冰载荷工况则出现剧烈波动，其中Mttnen海冰模型的影响更大，振动频率高，振幅更大，对塔架的疲劳载荷影响巨大.

图5为塔顶侧向位移的时域响应.相较于前后位移，塔顶侧向位移可忽略不计，三种工况下塔顶侧向位移相差不大.主要原因是海冰碰撞方向与风向相同，侧向作用力较小，对塔架作用不明显.

综合比较三种工况下塔顶的绝对位移发现，与无海冰载荷工况相比，Mttnen海冰使塔顶位移最大值增大24.1%，IEC Lock in海冰模型的增幅为16.0%.

挥舞和摆振分别为叶尖部分在旋转平面外和旋转平面内的位移，挥舞过大可能会与塔架发生碰撞造成叶片损坏，摆振过大容易使叶片发生断裂.图6为叶片在三种工况下的挥舞特性.从图中可看出，海冰载荷对叶片挥舞的振动幅值影响不大，这间接说明叶片挥舞主要受到湍流风的影响.这两种海冰模型下叶片挥舞特性比较类似.由于海冰与风力机支撑结构的碰撞，使得风力机的晃动相对更为剧烈，导致叶片振动加剧，增加了叶片的疲劳载荷.endprint

图7为三种工况下叶片摆振时域响应.从图中可看出，三种工况下叶片摆振特性大体趋势、幅值均相差很小，但在Mttnen海冰作用下叶片摆振方向有小范围的剧烈波动.相较于海冰对叶片挥舞的影响，海冰载荷对叶片摆振的影响可忽略不计.

5 结 论

通过对NREL 5 MW近海4桩柱风力机在Mttnen模型海冰作用下、IEC模型海冰作用下和无海冰作用下对比，得出以下结论：

（1） 相较无海冰作用，在Mttnen海冰作用和IEC模型海冰作用下，近海桩柱式风力机动力学响应明显增大.

（2） 在三种不同海冰作用下，塔架前后位移、叶片挥舞波动趋势不同；塔架侧向位移、叶片摆振波动趋势相同.

（3） Mttnen模型海冰作用导致风力机动力学响应最为剧烈，IEC模型海冰作用次之，无海冰作用相较平缓.

参考文献：

[1] 李俊峰，蔡丰波，乔黎明，等.2014中国风电发展报告[R].北京：中国循环经济协会可再生能源专业委员会，2014.

[2] 郑崇伟，潘静.全球海域风能资源评估及等级区划[J].自然资源学报，2012，27（3）：364-371.

[3] 周超，朱熀球，魏杰.我国风力发电发展现状和问题分析[J].能源研究与信息，2012，28（2）：69-75.

[4] VAN DER TEMPEL J，DIEPEVEEN N F B，DE VRIES W E，et al.15offshore environmental loads and wind turbine design：Impact of wind，wave，currents and ice[J].Wind Energy Systems：Optimising Design and Construction for Safe and Reliable Operation，2010：463-478.

[5] 高月文，李春，叶舟，等.海上风力机桩柱式结构动力响应分析[J].上海理工大学学报，2013，35（6）：591-595.

[6] 陈严，蔡安民，叶枝全，等.近海风力机在极限波浪作用下的初步计算分析[J].太阳能学报，2008，29（2）：180-187.

[7] 陈法波，李昕，周晶.近海风力涡轮机所受随机空气动力荷载模拟研究[J].太阳能学报，2011，32（10）：1528-1532.

[8] MRZ A，HOLNICKISZULC J，KRN T.Mitigation of ice loading on offshore wind turbines：Feasibility study of a semiactive solution[J].Computers & Structures：2008，86（3-5）：217-226.

[9] 李静，陈健云，陈小波，等.海上风力发电体系在风冰作用下的动力反应[J].太阳能学报，2013，34（7）：1234-1240.

[10] 柳春图，段梦兰.海冰工程中的结构力学问题[J].机械强度，1995，17（3）：7-20.

[11] 董吉武.海冰对圆桩和斜坡结构作用及其防护设施的物理模拟研究[D].大连：大连理工大学，2012.

[12] 黄焱.冰激海洋平台振动的动力模型试验研究[D].天津：天津大学，2004.

[13] International Electrotechnical Commission.IEC 614003 wind turbinepart3：Design requirements for offshore wind turbines[S].Geneva，Switzerland，2009.

[14] MTTNEN M.Numerical model for iceinduced vibration load lockin and synchronization[C]∥Proceeding of the 14th International Symposium on Ice.NY USA：Postdam，1998.

[15] 刘全，高伟，李春.兆瓦级海上风力机叶片设计及模态分析[J].能源研究与信息，2012，28（2）：112-115.

[16] 刘延柱，洪嘉振.多刚体系统动力学[M].北京：高等教育出版社，1989.

[17] 高偉.漂浮式风力机整机动力学特性分析与仿真[D].上海：上海理工大学，2013.

[18] 王磊.海上风电机组系统动力学建模及仿真分析研究[D].重庆：重庆大学，2011.

[19] 史庆增，宋安.海冰静力作用的特点及几种典型结构的冰力模型试验[J].海洋学报，1994，16（6）：133-141.

[20] 徐继祖.海冰引起的结构振动[J].海洋工程，1986，4（2）：42-47.

[21] 王翎羽，陈星，宋安，等.柔性海洋结构冰激振动问题的对策[J].中国海上油气（工程），1994，6（4）：11-17.

[22] 李明.冰激柔性直立结构振动的试验研究[D].天津：天津大学，2004.

[23] KORZHAVIN K N.Action of ice on engineering structures[M].Hanover USA：US Army CRREL Translation TL260，1971.

[24] 刘健，陈国明.基于Mttnen模型的冰激疲劳寿命分析[J].石油大学学报：自然科学版，2004，28（4）：89-93.

[25] 高月文.漂浮式风力机张力腿平台及系泊系统动力响应仿真研究[D].上海：上海理工大学，2014.

[26] 王龙文.大气风场模型研究及应用[D].长沙：国防科学技术大学，2009.

[27] 吴攀.风力机极限风况及地震载荷研究[D].上海：上海理工大学，2014.endprint