导航：安全管理网>> 安全论文>> 电力安全论文>>正文

流动空气与风力发电叶片作用特性研究

文档作者： Muhammad Aqeell，高延伟2，李永强3，刘青1，张照煌1 文档来源： 1．华北电力大学能源动力与机械工程学院 2．中国建筑出版社 3．山西省忻州市滹沱河水利管理局		点击数：	更新时间： 2024年03月05日
下载地址：点击这里	文件大小： 717.45 KB 共5页	文档格式： PDF	下载点数： 1 点（VIP免费）

全屏查看

内容预览 [文件共5页]

本文件共5页，如需编辑使用，请下载。

注：预览效果可能会出现部分文字乱码（如口口口）、内容显示不全等问题，下载是正常的。

	文件大小：717.45 KB 共5页文件格式：PDF 下载点数：1 点（VIP会员免费）

下一篇：丽水通济堰灌区水利管理体系的演进与启示

上一篇：麻阳二级水电站报废重建水轮发电机组的选型

文本预览

仅提取页面文字内容，供快速阅读使用。

第15卷第6期中国水利水电科学研究院学报v。1．15 N。．6
2017年12月Journal of China Institute of Water Resources and Hvdropower Research December，2017
文章编号：1672—3031(2017)06—0425—05
流动空气与风力发电叶片作用特性研究
Muhammad Aqeell，高延伟2，李永强3，刘青1，张照煌1
(1．华北电力大学能源动力与机械工程学院，北京102206；2．中国建筑出版社，北京100037；
3．山西省忻州市滹沱河水利管理局，山西忻州034000)
摘要：流动空气与风力发电叶片作用特性是风力发电叶片设计的基础，是风力发电领域研究的难题。长期以来以
动量一叶素理论为基础，以流过风力发电叶片空气流场为对象，相继建立了叶片优化设计的Schmitz方法、Clauert
方法、wilson方法和cFD方法等。但流经叶片过程中的空气对叶片作用的特性仍有许多问题需要亟待解决和深入
研究。为此，论文应用伯努利方程推导了环状流动空气流过风力发电风轮时传递给风轮的能量计算解析式；应用
动量矩原理推导了环状流动空气流过风力发电风轮时作用在风轮上的力矩解析式；提出了流过风轮的环状空气向
风轮释放功率的解析式，从而发现风轮正常工作时，流过风轮的环状空气向风轮释放的功率等于环状空气流量、
空气密度、环缘处风轮的牵连运动速度及环缘处流出空气的绝对运动速度在此点牵连运动速度方向上投影的连乘
积。为叶片截面轮廓曲线后边缘特征参数设置提供了参考。
关键词：风力发电叶片；伯努利方程；动量矩原理
中图分类号：TK83 文献标识码：A doi：10．13244，j．cnki．jiwhr．2017．06．003
1 研究背景
风力发电叶片是捕获风能并将之转化为机械能的核心部件，其气动外形对实现这一功能具有至
关重要的作用，因此，一直以来，国内外学者对其气动外形设计理论的研究都非常重视。文献[1]根
据风力发电叶片优化设计目标的不同将风力发电叶片优化设计分为以风能利用系数(c。)最大或者年
能量输出(Annual Energy Production，AEP)最大为目标的优化设计；以能量成本(cost of Energy，
c0E，风力机总成本C与年能量输出AEP的比值)最小为目标及多学科、多目标的优化设计。风力发
电叶片优化目标确定后，一般采用动量一叶素理论。2。、Schmitz方法”1、Glauert方法”1或wilson方法”。
等对之实施优化设计，而后辅之以叶片外形修正以实现设计要求‘6。。随着计算机技术的发展，以Na—
vier-stokes(N—s)方程为基础的计算流体力学(computational F1uid Dynamics，简称cFD)法等应用于风
力发电叶片优化已逐渐受到人们的重视"11。研究发现，对风力发电叶片进行优化设计，无论建立什
么模型、采用哪种优化方法，一般步骤是：首先选定叶片总体参数(如：额定功率P，叶片数曰，风
能利用系数C，，叶尖速比入，设计风速％，机电效率叼，空气密度p，风轮直径D，风轮转速n等)；
然后将叶片分段，设定叶素，选择对应翼型，即确定翼型沿叶片展向的布置；最后，根据优化目标，
选定的优化方法进行优化设计。91。优化设计变量一般是①弦长、扭角(或桨矩角。1叫)；②弦长、扭
角、桨矩角¨“；③弦长、扭角、叶尖线速度和叶尖速比‘1“；④外形参数、翼型特征、风轮转速和
叶尖桨矩角【1“；⑤更多¨41等。从中可看出对风力发电叶片进行优化设计，一般不包括叶片截面轮廓
曲线特征参数的优化，因为叶片截面轮廓曲线是通过选择翼型确定的，国外常用风力发电叶片专用
翼型有[15。1“：NREL的s系列翼型、Du系列翼型、Ris仍系列翼型和FFA系列翼型等。因商业保密等原
因，国内叶片厂商很难得到这些翼型的相关数据，一般靠测绘仿制。因此，对流动空气与风力发电
收稿日期：2016一05—26
作者简介：Muhammad Aqeel(1985一)，男，巴基斯坦，博士生，主要从事可再生能源高效利用理论与技术等方面的研究。
E一瑚ail：maqeell@qq．(’om
一425—
万方数据
流动空气与锋利发电叶片作用特性研究Muhammad Aqeel高延伟李永强刘青张照煌
叶片作用理论进行研究就显得非常重要。文中通过对流经风力发电风轮空气建模，应用伯努利方程
推导出环状流动空气流过风力发电风轮时传递给风轮的能量计算解析式；应用动量矩原理推导出环
状流动空气流过风力发电风轮时作用于风轮上的力矩解析式；通过对导出的解析式的分析，提出了
流过风轮的环状空气向风轮释放的功率的解析式，从而发现风轮正常工作时，流过风轮的环状空气
向风轮释放的功率等于环状空气流量、空气密度、环缘处风轮的牵连运动速度及环缘处流出空气的
绝对运动速度在此点牵连运动速度方向上的投影的连乘积。从而为叶片截面轮廓曲线后缘特征参数
设置提供了依据。
2风轮的能量表达式
以叶片半径为r、在半径方向长为6fr的圆环状流体为研究对象(见图1所示)。设空气流人叶片的
初始绝对运动速度、初始相对运动速度、初始牵连运动速度、初始压力能头和空气相对于零位能面
的高度分别为瓦、瓦、i、p，和z。，相应的，空气流出叶片的初始绝对运动速度、初始相对运动速度、
初始牵连运动速度初始压力能头、和空气相对于零位能面的高度分别为”：、训：、u：、p：和z：，此过程
中的能量损失为Σ允。空气刚流入叶片间的能量E。为，
牛矿嚣+若(1)
空气刚从叶片间流出的能量E：为，
即矿等+若(2)
考虑此过程中的能量损失Σ^，则圆环状流体输送给风轮的能量E为：
E=El—E2一Σ矗(3)
将式(1)和(2)代入式(3)并整理，得：
E=(rz：)十等+等一Σ危(4)
在叶片间流动的空气，遵循伯努力方程，即：
”毒+若剐：+等+号+Σ^ ㈦
现在假设给风轮一个大小与其转速相等、旋转方向与其相反的转速，则相当于在刚流人叶片问
的单位重量的流体上加了一负能量苦，在刚要流出叶片间的单位重量的流体上加了一负能量专，则
式(5)的能量转换还可写为：
”等+等～等≈：十等+等一号+Σ^ c6，
将式(6)代入式(4)并整理，得
肚寻+警一等㈩
考虑环域内有：u。=u：，所以式(7)可表示为： EE=：lF≮三一兰—一-_半二 (8)
由图中△ABc得，阳12=H?+∞?=u；+可。2代入式(8)并整理，得： E：兰掣E=—L_：』一 ((99))
一426—
万方数据
流动空气与锋利发电叶片作用特性研究Muhammad Aqeel 高延伟李永强刘青张照煌
6
|
| 八
§裂一￡一
／
7
1．风及风速：
2．风力发电风轮；
3．轮毂；
4．风力发电叶片；
5．空气环；
6．风力发电叶片背风面曲线
7．空气流出风轮速度三角形
8．风力发电叶片迎风面
Il，J发}UM轮俯眦㈩
图1 流动空气与风力发电叶片相互作用分析示意图
此即所考察的环状流动空气流过风力发电风轮时传递给风轮的能量。由此可看出，环状流动空
气流过风力发电风轮时传递给风轮的能量只与此环状流动空气刚流出风轮的状态有关，即，只与此
环状流动空气刚流出风轮的绝对运动速度、牵连运动速度和相对运动速度有关，且随空气相对于风
力发电叶片运动的运动速度的平方的增加而增加，随该点的牵连运动速度平方和相对运动速度平方
的增加而减小。
3风轮的力矩表达式
流动空气流过风轮过程中，其与风力发电叶片的相互作用满足动量矩定理。
刚流入风轮时对风轮产生的动量矩M．为：
Ml=F×(p27crdrf"1)=o
式中：p为空气密度；Z为r处叶片宽度。
环域内空气刚流出风轮时对风轮产生的动量矩M：为：
M2=i×(P27crdrf口2)=o
由图1可得，M，的大小M，为： M：=p2冗rdrf”：rc。s(瓦，夏)=p2兀rdrf”：r!寄
设考察环域内空气
(10)
(12)
式中各符号意义同前。
假设环域空气流过叶片宽度z的时间为￡，根据动量矩原理，环域内空气对风轮产生的力矩M为
一427—
万方数据
流动空气与锋利发电叶片作用特性研究Muhammad Aqeel高延伟李永强刘青张照煌
M：生当：P掣v窖兰蔓硝，堕害王⋯)
f ￡
二“2"2 二“2
式中：q为环域内空气流量；其它符号意义同前。
此即所考察的环状流动空气流过风力发电风轮时作用于风轮上的力矩。由此可看出，环状流动
空气流过风力发电风轮时作用于风轮上的力矩除与空气自身的密度、空气流量和考察点半径有关
外，就只与此环状流动空气刚流出风轮的状态有关，即，只与此环状流动空气刚流出风轮的绝对运
动速度、牵连运动速度和相对运动速度有关。
4风轮的功率表达式
式(9)给出的能量是以水头表示的，所以其对应的功率Ⅳ，为： ⅣF：印gg：竺警Pg (14)
式(13)给出的功率Ⅳ。可表示为： Ⅳ舸：彤(tJ：胖2r竺掣∞：_pq竺掣(15) M划舭胖2r二斧帖_pq』f ⋯’
式中叫为风轮旋转角速度；其它符号意义同前。
比较式(14)、(15)发现，流过风力发电风轮的空气，无论从能量的角度进行分析，还是从动量
矩定理进行分析，流动空气传递给风轮的功率的表达式都完全一致。即，只与此环状流动空气刚流
出风轮的状态有关，具体来说，也就是只与此环状流动空气刚流出风轮的绝对运动速度、牵连运动
速度和相对运动速度有关。
根据图1中△EFG，由余弦定理，得
加加；2==““2；+移+移2一；一z“2“2秽：2秽co：scloⅡs(2乏，，秽五2)J L(1l6o，)
将式(16)代入式(14)或式(15)，得
Ⅳ=Ⅳg=Ⅳ^f=一pqu2口2coslu2，。2} (17)
式(17)表示，风轮正常工作，流动空气流出风轮叶片的绝对运动速度与该点牵连运动速度的夹
角必为钝角；流过风轮的环状空气向风轮释放的功率等于环状空气流量、空气密度、环缘处风轮的
牵连运动速度及环缘处流出空气的绝对运动速度在此点牵连运动速度方向上投影的连乘积。
5结论
通过分析流经风电叶片间的能量流和动量矩的传递发现，流动空气与风力发电叶片间的相互作
用如下：(1)环状流动空气流过风轮过程中，若向风轮传递功率，其流出风轮叶片的绝对运动速度与
该点牵连运动速度的夹角必为钝角；(2)流过风轮的环状空气向风轮释放的功率等于环状空气流量、
空气密度、环状空气从叶片后缘流出的牵连运动速度及其绝对运动速度在此点牵连运动速度方向上
投影的连乘积；(3)改变风力发电叶片迎风面后缘曲线切线角(相当于改变空气从风力发电叶片后缘
流出角)，可改变流过风轮的环状空气传递给风轮的功率。
参考文献
[ 1 ] DANIIL Per“liev，Methodology for wind turbine blade geome—t。yoptimization[D]
sity of Tech—nology，2叭3．
[2]段宏，赵明，梁俊字，等．基于动量叶素理论的风力机叶轮建模[J]．云南电力技术，2016，44(2)
一428—
万方数据
流动空气与锋利发电叶片作用特性研究Muhammad Aqeel 高延伟李永强刘青张照煌
[ 3 ]
[ 8]
[9]
10
11
12
[13]
[14]
[15]
[16]
130一133．
李连波，陈涛，王凡，等．基于schmitz理论的风叶气动设计研究[J]．应用力学学报，2012，29(2)：
225—229．
方程．基于涡流理论的立轴风力机气动特性分析[D]．北京：华北电力大学，2012．
钟友富．基于Wilson理论的大型风力机叶片三维实体建模[J]．科技创新，2015(10)：39—40．
李强，金哲岩，杨志刚．风力机叶片气动外形优化设计综述[J]．机械设计，2016，33(7)：1—7．
姚世刚，戴丽萍，康顺．风力机叶片气动性能及流固耦合分析[J]．工程热物理学报，2016，37(5)：
988—992．
徐宇，廖猜猜，张淑丽，等．大型风电叶片设计制造技术发展趋势[J]．中国科学：物理学力学天文学，
2016．46(12)：124702—1一10．
GUO Tingting，WU Dianwen，XU Jihui，et a1．The method oflarge—scale wind turbine blades design based on
MATLAB programmingl c]，／Nanjing：sustainable Power Generation and supply，2009：1—5．
w00D D H．Dual purpose design of small wind turbine blades[J]．wind Engineering，2004，28(5)：511—528，
吕小静．大型水平轴风力机风轮气动性能计算与优化设计[D]．兰州：兰州理工大学，2011．
ERNESTO Benini，ANDREA Toffolo．0ptimal design of horizontal—axis wind turbines using blade—element theo—
ry and evolutiona。y computation l J J．Journal of Solar Energy Engineering，2002(1 24)：357—363．
FUGLSANG P，MADSEN H A．0ptimization method for wind turbine rotors[j]，journal of wind Enginee“ng
and Industrial Aerodynamics，1999(80)：191—206．
WANG Tongguang，WANG Long，ZHANG Wei，et aI．Large—scale wind turbine blade design and aerodynamic
analvsis[J]．Fluid Mechanics，2叭2．57(5)：466—472．
NAMBU T，HASHIMOTO A，UEN0 M，et a1．Evaluatjon of wall—jnfeI如rence correc￡ion method using numerical
anal丫sis and porous wall modellJj．Aircraft，2015，52：226—234．
KANG S H，SHIN E S，RYU K W，et a1．Separation blockage—correction method for the airfoil of a wind turbine
bladelJ]．Mech．Sci．Tech．，2013，27：l 321一1327．
Study on characteristics of the interaction between nowing air and wind turbine blade
Muhammad Aqeel。， GAo Yanwei2， LI Yongqianf， LIu Qin91，zHANG zhaohuan91
(1．School《Energy P0wer oRd Mechonicoz Engineertng，～onh chtn8 Ekc￡rtc Po蜘er Un抽ersi姆，Be电tng 102206，Chino；
2．C^inese口rc充眺cfur口Z p阳5s， Be盯ing 100037， C矗inⅡ；
3．HufⅡo^e Ri口er Bns胁缈口ler Resource5 Mnnngemen￡曰ure口u矿5^nnxi Pro"ince， Xinz^oⅡ 034000， C^in。)
Abstract： The characteristic of interaction between nowing air and wind turbine blade is the foundation in
the wind turbine blade design and has been a dimcult problem in the wind power generation． For a long
time， taking the air now field of wind turbine blades as the objective and on the basis of momen—
tum—blade element lheorv， the Schmitz method， Glauert method， Wilson method and CFD memod have
been established to optimize the balde design． However， still a lot of problems need fu rther improvement
and study in view of the characteristics of the effect of air on blades during the air nowing through blades．
Therefore， the analytical formula of energy transferred from annular nowing air to rotor was proposed based
on Bernoulli equation， the analvtical fbrmula of torque on the rotor under the annular flow situation was
put forward with the principle of momentum， and the analytical formula of transmission power from annular
nowing air to I时“)r was presented． It is obta；ned that the transmission power is 8qual to the continuous
product of annular air now， air density， convected velocity of rotor at rim and the projection of the abso—
lute veloeity of the air outflowi“g around the rim． It provides guide for setting trailing edge feature parame—
ters of blade profile curve．
KeywOrds： wind turbine blade； bernoulli equation； momentum moment principIe
(责任编辑：李琳)
一429—
1j
1●J
4
5
6
7
『L

网友评论 more

电力安全论文最新内容

03-05