[网站公告] [设置首页] []

首 页

找产品

找企业

资讯中心

采购中心

项目招标

风能技术

专访报道

展会会议

人才招聘

风能书屋

百科资料

欢迎关键词竞价推广,热门搜索: nordex  叶片  基础环  株洲    请选好搜索类型!
当前位置:首页 » 资讯 » 技术市场 » 风能技术 » 正文
风电机组叶片流固耦合的数值模拟方法
2016-10-13 来源:中国风能设备网
 随着风电机组单机容量的不断增大,风电机组叶片尺寸也越来越长。对于超长风电机组叶片的气动计算及结构计算,常规的叶片刚性体假设所引起的误差越来越大,因此必须考虑叶片在风载条件下结构变形对气动力的影响。流固耦合(fluid-structure interaction, FSI)是研究风电机组叶片气弹特性及结构特性的重要手段,国内外通过流固耦合对叶片的结构响应、流场变化、叶片载荷分布等特性进行了一系列研究。流固耦合按照求解方法可分为完全耦合和交替求解。李媛等以MPCCI为数据交换平台连接商用CFD,CSD软件,模拟考虑风切变时轮毂高度额定风速下的风轮全三维流固耦合,分析了风剪切对叶片绕流、变形以及耦合作用对载荷分布的影响规律。胡丹梅等采用k-SST紊流模型和滑移网格技术,对海上风电机组叶片进行了流固藕合计算分析,发现叶片流固耦合作用使叶片气动攻角、扭矩增大。Zahle等采用不可压缩重叠网格技术,分析了水平轴风电机组风轮与塔架间的耦合作用。Hsu等通过对5 MW三叶片水平轴风电机组风轮在有无塔架情况下进行了流固耦合研究,发现叶片转过塔架时由于塔影效应使单个叶片气动扭矩降低10%-12%。 Zhang等通过对比研究,发现风切变对风电机组叶片变形及应力的影响明显大于流固耦合所产生的作用。任年鑫等基于三维黏性不可压缩Navier-Stokes方程和重整化群k-ε(RNG)湍流模型,数值模拟美国可再生能源实验室(NREL)5 MW海上风电机组的气动性能,并研究了浮式平台不同运动幅值及运动周期对风电机组叶片气动性能的影响规律。陈海萍等选用k-ε湍流模型,将风作为黏性不卡压缩流体,计算出流固耦合作用下风电机组叶片的应变分布。李德源等采用计算多体系统动力学理论,基于R-W(Roberson-Wittenburg)建模方法,建立柔性叶片的气弹耦合方程,结合相应的风电机组气动载荷分析模块,实现了叶片的气弹耦合分析。王旭东等将叶片的启动与结构力学模型进行耦合,研究叶片变形对来流的影响,提出一种旋转叶片结构性能分析方法。Dong等应用CFD-CSD祸合方法对风电机组叶片的载荷及气弹特性进行研究,发现叶片气弹变形显著降低气动载荷。MEXICO (model experiment incontrolled conditions)实验是由欧盟资助的一项大型风电机组实验。梁明轩通过对叶片模态分析得到叶片各阶模态,发现叶片扭转、挥舞、摆振三者之间的耦合,并对叶片气动弹性问题进行研究,建立叶片颤振方程,确定了叶片颤振发散条件。杨华等应用风电机组旋转叶片表面压力的三维PIV(particle image velocimetry)流场测量方法,测试得到无偏航工况下叶片上5个测试断面的表面压力系数分布,并计算了该工况下叶轮所受的轴向力和切向力。潘旭运用Ansys、CFX基于Ansys Workbench软件平台对某MW级水平轴风电机组叶片进行流固耦合数值模拟,分析了藕合过程中叶片的变形及应力变化。陆洋等利用Newmark数值积分法获得叶片气弹响应的稳态周期解,分别以NREL Phase VI非定常空气动力学实验及其公开的1.5 MW风电机组叶片为算例计算了有/无预弯叶片的气弹响应。以上研究均通过建立风轮旋转的计算模型对风轮流固祸合进行研究,网格划分困难,涉及到移动域与静止域网格匹配等复杂过程,计算成本高。为此,本文采用风切变形式模拟风轮旋转及来流风速的综合效应,研究风电机组叶片流固耦合数值模拟的一种新方法。

1流体计算模型

1.1流体域求解方法

流场中的控制体既不能产生质量也不能消灭质量,因此它满足质量守恒定律。该定律在流场中的数学表达式为连续性方程

式中,p,t,V分别为流体密度、时间和流体速度。

动量守恒定律可表述为:控制体中流体的动量对时间的变化率等于外界作用在该微元体上的各种力之和。其在流场中的数学表达式为运动方程

能量守恒定律是流动系统必须满足的基本定理,在不考虑温度变化条件下,其数学表达式为

式中,F为质量力,P 为表面压强。

1.2叶素速度

风电机组风轮运行过程中实际参考风速为来流风速和旋转风速的共同作用,通过改变叶片不同叶展处的风速可实现正常来流风速和旋转风速的叠加效果。不同叶展处的参考风速Vref的计算公式为

式中,Vo为自然来流风速;w 为风轮旋转角速度;r 为叶展位置的旋转半径;a为轴向诱导系数;a’为切向诱导系数。

叶素速度矢量的三角形关系如图1所示,图中α、β、Ø分

2结构计算模型

2.1结构计算方法

风电机组叶片为多自由度有阻尼的非线性结构,其结构动力学运动方程为

2.2叶片结构模型

为了用静止叶片模拟叶片旋转效应,只通过将叶素处的速度按照速度矢量关系修正还不能达到应有效果,需要对叶片不同截面处翼型的气动扭角进行修正。如图1所示,以旋转导致的风速方向作为参考方向,将该截面处的翼型进行顺时针旋转,旋转角θ等于人流角Ø。新翼型各点坐标相对于原坐标的计算公式为

 

ewindpower.cn版权及免责声明:
①凡粘转自“中国风能设备网”的文章,转载时请注明来源“中国风能设备网”;
②部分文章转载自其它媒体,转载目的在于传递更多信息,其信息仅为原作者个人观点,不作任何投资参考依据;
③如有作品内容、版权和其它问题需要同本网联系的,请在15个工作日内进行;

看完了,返回首页

同类资讯
• 基于风洞实验的浮动式海上风电机组载荷控制 • 谷歌的“能源风筝”招徕壳牌合作 预计今年
• 风电零件特种修复方法之冷熔 • 风机维修中常见的九大问题
• 风电行业中几种新技术的简介 • 谷歌能量风筝 要去风电无法到达的地方
• 新技术 | 高空风力发电机 未来可期! • 早期延寿规划可将风电项目利润增加2%
• 风力发电机机舱罩轮毂罩行业现状分析及新工 • 聚氨酯叶片时代即将到来

  图片新闻:
2017年风电整机商吊装容量数据公布
2017年风电整机商吊装容量数据公布
国家能源局:新疆、甘肃、辽宁、吉林、宁夏弃风率下降超过10个百分点
国家能源局:新疆、甘肃、辽宁、吉林、宁
25家风电企业登上《财富》世界500榜单
25家风电企业登上《财富》世界500榜单
一二审股民皆胜诉 华锐风电诉讼时效倒数百日
一二审股民皆胜诉 华锐风电诉讼时效倒
东方风电自主开发直驱电控系统完成低电压穿越测试
东方风电自主开发直驱电控系统完成低电
中国大唐连续五年荣获中央企业经营业绩考核A级
中国大唐连续五年荣获中央企业经营业绩
国神集团内蒙古风电公司白云40MW风电工程顺利移交生产
国神集团内蒙古风电公司白云40MW风电工
中电联4月底风电并网1.5亿千瓦 完成投资91亿元
中电联4月底风电并网1.5亿千瓦 完成投

推荐文章

点击排行

网站首页 | 会员升级 | 诚聘英才 | 关于我们 | 付款方法 | 联系方式 | 使用帮助 | 网站留言 | 广告服务 | 竞价推广
版权所有 中国风能设备网 © 2015-2019  电话:130-7219-8820  
E-mailL:ewindpower@163.com   QQ:810978040