Muhammad Aqeel,高延伟,李永强,刘 青,张照煌
(1.华北电力大学能源动力与机械工程学院,北京 102206;2.中国建筑出版社,北京 100037;3.山西省忻州市滹沱河水利管理局,山西 忻州 034000)
风力发电叶片是捕获风能并将之转化为机械能的核心部件,其气动外形对实现这一功能具有至关重要的作用,因此,一直以来,国内外学者对其气动外形设计理论的研究都非常重视。文献[1]根据风力发电叶片优化设计目标的不同将风力发电叶片优化设计分为以风能利用系数(CP)最大或者年能量输出(Annual Energy Production,AEP)最大为目标的优化设计;以能量成本(Cost of Energy,COE,风力机总成本C与年能量输出AEP的比值)最小为目标及多学科、多目标的优化设计。风力发电叶片优化目标确定后,一般采用动量-叶素理论[2]、Schmitz方法[3]、Glauert方法[4]或 Wilson方法[5]等对之实施优化设计,而后辅之以叶片外形修正以实现设计要求[6]。随着计算机技术的发展,以Na⁃vier-Stokes(N-S)方程为基础的计算流体力学(Computational Fluid Dynamics,简称CFD)法等应用于风力发电叶片优化已逐渐受到人们的重视[7-8]。研究发现,对风力发电叶片进行优化设计,无论建立什么模型、采用哪种优化方法,一般步骤是:首先选定叶片总体参数(如:额定功率P,叶片数B,风能利用系数CP,叶尖速比λ,设计风速v0,机电效率η,空气密度ρ,风轮直径D,风轮转速n等);然后将叶片分段,设定叶素,选择对应翼型,即确定翼型沿叶片展向的布置;最后,根据优化目标,选定的优化方法进行优化设计[9]。优化设计变量一般是①弦长、扭角(或桨矩角[10]);②弦长、扭角、桨矩角[11];③弦长、扭角、叶尖线速度和叶尖速比[12];④外形参数、翼型特征、风轮转速和叶尖桨矩角[13];⑤更多[14]等。从中可看出对风力发电叶片进行优化设计,一般不包括叶片截面轮廓曲线特征参数的优化,因为叶片截面轮廓曲线是通过选择翼型确定的,国外常用风力发电叶片专用翼型有[15-16]:NREL的S系列翼型、DU系列翼型、RisØ系列翼型和FFA系列翼型等。因商业保密等原因,国内叶片厂商很难得到这些翼型的相关数据,一般靠测绘仿制。因此,对流动空气与风力发电叶片作用理论进行研究就显得非常重要。文中通过对流经风力发电风轮空气建模,应用伯努利方程推导出环状流动空气流过风力发电风轮时传递给风轮的能量计算解析式;应用动量矩原理推导出环状流动空气流过风力发电风轮时作用于风轮上的力矩解析式;通过对导出的解析式的分析,提出了流过风轮的环状空气向风轮释放的功率的解析式,从而发现风轮正常工作时,流过风轮的环状空气向风轮释放的功率等于环状空气流量、空气密度、环缘处风轮的牵连运动速度及环缘处流出空气的绝对运动速度在此点牵连运动速度方向上的投影的连乘积。从而为叶片截面轮廓曲线后缘特征参数设置提供了依据。
以叶片半径为r、在半径方向长为dr的圆环状流体为研究对象(见图1所示)。设空气流入叶片的初始绝对运动速度、初始相对运动速度、初始牵连运动速度、初始压力能头和空气相对于零位能面的高度分别为和z1,相应的,空气流出叶片的初始绝对运动速度、初始相对运动速度、
初始牵连运动速度初始压力能头、和空气相对于零位能面的高度分别为和z2,此过程中的能量损失为∑h。空气刚流入叶片间的能量E1为,
空气刚从叶片间流出的能量E2为,
考虑此过程中的能量损失∑h,则圆环状流体输送给风轮的能量E为:
将式(1)和(2)代入式(3)并整理,得:
在叶片间流动的空气,遵循伯努力方程,即:
现在假设给风轮一个大小与其转速相等、旋转方向与其相反的转速,则相当于在刚流入叶片间的单位重量的流体上加了一负能量在刚要流出叶片间的单位重量的流体上加了一负能量则式(5)的能量转换还可写为:
将式(6)代入式(4)并整理,得
考虑环域内有:u1=u2,所以式(7)可表示为:
图1 流动空气与风力发电叶片相互作用分析示意图
此即所考察的环状流动空气流过风力发电风轮时传递给风轮的能量。由此可看出,环状流动空气流过风力发电风轮时传递给风轮的能量只与此环状流动空气刚流出风轮的状态有关,即,只与此环状流动空气刚流出风轮的绝对运动速度、牵连运动速度和相对运动速度有关,且随空气相对于风力发电叶片运动的运动速度的平方的增加而增加,随该点的牵连运动速度平方和相对运动速度平方的增加而减小。
流动空气流过风轮过程中,其与风力发电叶片的相互作用满足动量矩定理。设考察环域内空气刚流入风轮时对风轮产生的动量矩为:
式中:ρ为空气密度;l为r处叶片宽度。
式中各符号意义同前。
假设环域空气流过叶片宽度l的时间为t,根据动量矩原理,环域内空气对风轮产生的力矩M为:
式中:q为环域内空气流量;其它符号意义同前。
此即所考察的环状流动空气流过风力发电风轮时作用于风轮上的力矩。由此可看出,环状流动空气流过风力发电风轮时作用于风轮上的力矩除与空气自身的密度、空气流量和考察点半径有关外,就只与此环状流动空气刚流出风轮的状态有关,即,只与此环状流动空气刚流出风轮的绝对运动速度、牵连运动速度和相对运动速度有关。
式(9)给出的能量是以水头表示的,所以其对应的功率NE为:
式(13)给出的功率NM可表示为:
式中ω为风轮旋转角速度;其它符号意义同前。
比较式(14)、(15)发现,流过风力发电风轮的空气,无论从能量的角度进行分析,还是从动量矩定理进行分析,流动空气传递给风轮的功率的表达式都完全一致。即,只与此环状流动空气刚流出风轮的状态有关,具体来说,也就是只与此环状流动空气刚流出风轮的绝对运动速度、牵连运动速度和相对运动速度有关。
根据图1中ΔEFG,由余弦定理,得
将式(16)代入式(14)或式(15),得
式(17)表示,风轮正常工作,流动空气流出风轮叶片的绝对运动速度与该点牵连运动速度的夹角必为钝角;流过风轮的环状空气向风轮释放的功率等于环状空气流量、空气密度、环缘处风轮的牵连运动速度及环缘处流出空气的绝对运动速度在此点牵连运动速度方向上投影的连乘积。
通过分析流经风电叶片间的能量流和动量矩的传递发现,流动空气与风力发电叶片间的相互作用如下:(1)环状流动空气流过风轮过程中,若向风轮传递功率,其流出风轮叶片的绝对运动速度与该点牵连运动速度的夹角必为钝角;(2)流过风轮的环状空气向风轮释放的功率等于环状空气流量、空气密度、环状空气从叶片后缘流出的牵连运动速度及其绝对运动速度在此点牵连运动速度方向上投影的连乘积;(3)改变风力发电叶片迎风面后缘曲线切线角(相当于改变空气从风力发电叶片后缘流出角),可改变流过风轮的环状空气传递给风轮的功率。
[1 ]DANIIL Perfiliev.Methodology for wind turbine blade geome-tryoptimization[D].Finland:Lappeenranta Univer⁃sity of Tech-nology,2013.
[2]段宏,赵明,梁俊宇,等 .基于动量叶素理论的风力机叶轮建模[J].云南电力技术,2016,44(2):130-133.
[3]李连波,陈涛,王凡,等 .基于Schmitz理论的风叶气动设计研究[J].应用力学学报,2012,29(2):225-229.
[4]方程.基于涡流理论的立轴风力机气动特性分析[D].北京:华北电力大学,2012.
[5]钟友富.基于Wilson理论的大型风力机叶片三维实体建模[J].科技创新,2015(10):39-40.
[6]李强,金哲岩,杨志刚.风力机叶片气动外形优化设计综述[J].机械设计,2016,33(7):1-7.
[7]姚世刚,戴丽萍,康顺 .风力机叶片气动性能及流固耦合分析[J].工程热物理学报,2016,37(5):988-992.
[8]徐宇,廖猜猜,张淑丽,等.大型风电叶片设计制造技术发展趋势[J].中国科学:物理学力学天文学,2016,46(12):124702-1-10 .
[9]GUO Tingting,WU Dianwen,XU Jihui,et al.The method oflarge-scale wind turbine blades design based on MATLAB programming[C]//Nanjing:Sustainable Power Generation and Supply,2009:1-5 .
[10]WOOD D H.Dual purpose design of small wind turbine blades[J].Wind Engineering,2004,28(5):511-528 .
[11]吕小静.大型水平轴风力机风轮气动性能计算与优化设计[D].兰州:兰州理工大学,2011.
[12]ERNESTO Benini,ANDREA Toffolo.Optimal design of horizontal-axis wind turbines using blade-element theo⁃ry and evolutionary computation[J].Journal of Solar Energy Engineering,2002(124):357-363 .
[13]FUGLSANG P,MADSEN H A.Optimization method for wind turbine rotors[J].Journal of Wind Engineering and Industrial Aerodynamics,1999(80):191-206 .
[14]WANG Tongguang,WANG Long,ZHANG Wei,et al.Large-scale wind turbine blade design and aerodynamic analysis[J].Fluid Mechanics,2012,57(5):466-472 .
[15]NAMBU T,HASHIMOTO A,UENO M,et al.Evaluation of wall-interference correction method using numerical analysis and porous wall model[J].Aircraft,2015,52:226-234.
[16]KANG S H,SHIN E S,RYU K W,et al.Separation blockage-correction method for the airfoil of a wind turbine blade[J].Mech.Sci.Tech.,2013,27:1321-1327 .