考虑气动阻尼浮式风机频域响应分析.docVIP

考虑气动阻尼浮式风机频域响应分析 　　摘 要：分析了气动阻尼对浮式风机频域响应的影响.选取美国可再生能源实验室（NREL）提出的5兆瓦（MW）浮式风机模型作为算例，利用气动阻尼计算方法建立气动阻尼矩阵，再基于三维势流理论计算浮式平台的水动力系数，并将系泊系统视为线性弹簧以考虑其刚度，最后在频域内分别建立并求解考虑与不考虑气动阻尼两种情况下的浮式风机刚体运动方程.利用求解频域方程得到的幅频响应算子（response amplitude operators，RAOs）及结合JONSWAP海浪谱得到的响应谱，在频域内分析了气动阻尼对浮式风机刚体运动的影响.结果表明：作业工况下气动阻尼能有效地降低纵荡和纵摇运动RAOs的峰值，且能在一定范围内减小对应自由度上响应谱的幅值和零阶矩. 　　关键词：海上风机；频域分析；气动阻尼；幅频响应算子；响应谱 　　中图分类号：P752 文献标志码：A 　　海上风能有着风速大、湍流??度低、储量丰富等特点，采用浮式风机是开发深水风能的首选技术手段.利用频域分析得到的RAOs及响应谱来评估浮式平台的动态性能，是其设计和优化中的重要内容：Bulder[1]等针对适用于5MW风机的三浮体式浮式平台，通过求解频域内浮式风机运动方程，得到了平台各运动自由度的RAOs及特定海况下的响应谱，并由此分析了该浮式风机在对应海洋环境下的适用性；唐友刚[2]等综合驳船式和单柱式（Spar式）平台的特点，提出了一种新型浮式平台方案，在频域内利用RAOs分析了波浪入射角和水深等因素对该平台运动的影响.而气动阻尼是多数学者在频域分析中忽略的因素. 　　Kühn [3]和Salzmann [4]指出，气动阻尼是海上风机支撑结构的振动与疲劳分析中的关键影响因素.邓露[5]等指出风机振动过大会对运行造成严重影响.Karimirad[6]等通过谱分析指出气动阻尼能有效抑制风机机舱在纵荡方向上的低频响应，且根据叶素动量理论指出浮式平台的刚体运动会受到气动阻尼的影响.但目前关于气动阻尼对浮式平台运动频域响应影响的研究较少：Roddier[7]等基于WindFloat浮式风机，通过模型实验得到浮式风机的RAOs，结果表明在研究中考虑与不考虑风荷载得到的RAOs相差不大，但该研究所涵盖的波浪周期和风速范围有限，且实验模型采用圆板代替风轮，不能准确地反映风机的气动特性；Wayman[8]等和Ramachandran[9]等联合利用分析软件WAMIT和FAST计算了多种浮式风机的RAOs并结合ISSC海浪谱计算了不同海况下各运动自由度的标准差，发现浮式风机不同于船舶和海洋平台，其RAOs的取值与海况有关，但该研究涵盖的工况较少，也没有深入探究影响气动阻尼的主要因素. 　　针对上述问题，本文在频域内分别建立并求解考虑与不考虑气动阻尼两种情况下的浮式风机运动方程.对比两种情况下RAOs和响应谱，在频域内按照不同的工况定量分析了气动阻尼对浮式风机响应的影响，并明确了影响气动阻尼的关键因素. 　　1 气动阻尼 　　1.1 气动阻尼的产生机理 　　气动阻尼源于风轮与空气的相互作用，如图1 [3]：叶素截面处于转速为Ω的风轮的径向r处，塔顶顺风向运动时，叶素处实际的轴向风速等于入流轴向风速V（1-α）减去塔顶速度top，导致叶片攻角α变小，在附着流的前提下并考虑升力系数CL，这将使得叶素受到的升、阻力变小，对应受到的推力dFx减小ΔdFx，整体上表现为风轮推力变小，故阻碍塔顶顺风向运动；同理，当塔顶逆风向运动时，风轮推力增大，进而阻碍塔顶的逆风向运动，则推力的变化始终阻碍着塔顶运动，这便是气动阻尼的机理. 　　1.2 气动阻尼的计算 　　利用NREL风机时域分析软件FAST建立风机模型，能考虑失速及小叶尖速比的情况，并能模拟具有复杂翼型的叶片[3-4].在FAST弹性动力模块中将结构阻尼设为零，并通过在软件中关闭浮式平台刚体运动自由度来限制平台运动以排除附加阻尼的影响，得到仅在气动阻尼作用下塔顶位移衰减曲线. 　　在FAST的气动模块中采用叶素动量理论计算气动荷载并考虑叶尖和叶根损失和采用Beddoes[10]等提出的动态失速模型.在FAST的控制模块中采用固定的风轮转速和浆距角，其具体值应根据不同的稳态风速确定[11]，并用指数模型来考虑竖直方向风剪切[12].通过时域模拟，得到了仅在气动阻尼作用下的塔顶振动衰减曲线，如图2所示（图中塔顶位移从对应风速下塔顶振动的平衡位置算起）.进而利用衰减曲线并结合式（1）和（2）计算对数衰减率和阻尼比 [3]： 　　式中：δ为塔顶位移的对数衰减率；ξaero为气动阻尼比；An，An+1 为相隔一个周期的两个位移峰值. 　　将风轮和机舱简化为塔架顶部的刚性质点，把塔架机舱转子体系视为有端部集中质