风力发电机组 叶片设计理论

12一、风力发电机风轮技术的发展一、风力发电机风轮技术的发展风轮是风力发电机吸收风能的关键部件，把风能转化成机械化能是由风力机的风轮完成的，它直接决定风力发电机的重要性能指标—风能利用系可靠的质量和优越的性能是保证机组正常稳定运行的决定因素。1888年，美国的CharlesF．Brush制成了世界上第一个用于发电的风力发电机，叶片采取平板设计，所以效率较低。31891年，丹麦的PoulLaCour教授首先将气体动力学引入风力发电机的研究，是世界上第一个利用风洞实验研究风力发电机的科学家，为设计和制造性能良好的风力发电机开辟了新途径。20世纪初，空气动力学的蓬勃发展和飞机的发明使人们对叶片的气动设计更为重视。1925年，美国的雅克布斯风力发电机生产厂发明了三叶片螺旋桨型叶轮，使水平轴叶轮的效率及性能均得到了较大提高。4随着风电技术的发展，风力发电机叶轮技术也在迅速发展。在空气动力学方面最重要的发展是进瑞典和丹麦等风能技术发达国家都发展各自的翼型系列，这些翼型系列各具特点。美国的SERI系列翼型具有较高的升阻比和一定的最大升力系数,失速时对翼型表面的粗糙度敏感性低,而NREL系列翼型能有效减小由于桨叶表面粗糙度增加而造成的风轮性能下降,并且能增加能量输出和改善功率控制；5丹麦的RISФA系列翼型在接近失速时具有较大的升阻比,并且具有对前缘粗糙度的不敏感性；瑞典的FFA-W翼型系列具有较高的最大升力系数和升阻比,并且在失速工况下具有良好的气动性同时,采用柔性叶片也是一个发展方向,利用新型材料进行设计制造,使其根据风况的变化能够改变空气动力型面,从而改进空气动力性能和叶片受力状况,增加可靠性和对风能的捕获量。6另外在开发新的空气动力装置上也进行了大量尝试,如在水平轴风力机叶端加一小翼以增加风力机的功率输出。对于叶片的材料,随着叶片尺寸的不断增大,主要朝大型化和轻量化的方向发展,由最初的木质品逐步过渡到玻璃钢,现在采用碳纤维复合材料(CFRP)的大型叶片制造技术是大型风力机制造7风轮的组成部件，主要是叶片。风力发电机的风轮，一般是由2-3片叶片组成的。为了理解叶片的功能，即它们是怎样将风能转变成机械能的，必须懂得有关翼型的空气动力学知识。8风力机空气动力学的基本概念叶片翼型的几何形状与空气动力特性风轮的组成部件，主要是叶片。风力发电机的风轮，一般是由2-3片叶片组成的。为了理解叶片的功能，即它们是怎样将风能转变成机械能的，必须懂得有关翼型的空气动力学知识。9现代风力发电机风轮叶片的剖面形状如图1所示。先考虑一个不动的翼型受到风吹时的情况。风的速度为矢量，方向与翼型平面平行，有关翼型几何形状定义如下：VLVLMθAiNB翼型的尖尾点B为后缘，圆头上A点称为前缘。连接前、后缘的直线AB=l,称为翼弦。AMB为翼型上表面，ANB称为翼型下表面。从前缘到后缘的虚线叫做翼型的中张线。迎角i是翼型与来流速度矢量之间的夹角。下面考虑风吹过叶片时所受的空气动力。翼剖面上的压力分部如图2所示。上表面压力为正，下表面压力为负，下表面压力为正。合力如图3所示。低压图2翼剖面上的压力图3翼剖面上的合力合力可用下式表达：式中：pS2F=2pCrSV—空气密度，kg/m3；—总的气动力系数。这个力可以分解为两个分力：垂直于气流速度的分力——升力平行于气流速度的分力——阻力FdFL和Fd可用下式表示：CL和Cd分别别为翼型的升力系数和阻力系数。由于这两个力互相垂直，所以图4翼剖面的升力系数随攻角变化的曲线图4翼剖面的升力系数随攻角变化的曲线图5翼剖面的阻力系数随攻角变化的曲线假设将受影响的空气与那些没有经过风轮圆盘、没有减速的空气分离出来，那么就可以画出一个包含受到影响的空气团的边界面，该边界面分别向上游和下游延伸，从而形成一个横截面为圆形的长的气流管（如图6所示）。如果没有空气横穿边界面，那么对于所有的沿气流管流向位置的空气质量流量都相等。但是因为流管内的空气减速，而没有被压缩，所以流管的横截面积就要膨胀以适应减速的空气。风力机的存在导致上游剖面接近风轮的空气逐渐减速以至于当空气到达风轮圆盘时其速度已经低于自由流风速了。风速的降低导致了流管膨胀，因为其速度没有对气体或通过气体来做功，所以气体的静压将上升以吸收其动能的减少。当空气经过风轮圆盘时显然有静压降存在，以至于空气离开风轮时其压力会小于大气压力。空气流就会以减小的速度和静压向下游前进——这个气流域被称为尾流。最终，为了保持平衡，下游远端尾流的静压要与大气压保持一致。动能的消耗使静压增加，从而导致风速进一步降低。因此在上游剖面远端和尾流远端之间，静压没有发生变化，但是有动能减一、贝茨（Betz）理论二、葛劳渥旋涡理论三、动量理论四、叶素理论五、schmitz理论六、叶栅理论贝茨理论的建立，是假定风轮是“理想”的，全部接受风能（没有轮毂），叶片无限多，对空气流空气流是连续的，不可压缩的，叶片扫掠面上的气流是均匀的，气流速度的方向不论在叶片前或叶片后都是垂直叶片扫掠面的（或称平行风轮轴线的），这时的风轮称为“理想风轮”。ss2v1s1vv2vv分析一个放置在移动的空气中的“理想风轮”叶片上所有受到的力及移动空气对风轮叶片所做的风轮扫掠后的风速面积为s2。风吹到叶片上所做的功是将风的动能转化为叶片转动的机械能，则必有v于是于是风单位时间作用在叶片上的力由动量定理求得式中：PF=Psv(v2-v1)—空气当时的密度，kg/m3；—单位时间流过的空气质量，kg。风轮所接受的功率为N=Fv=psv2(v2-v1)经过风轮叶片的风的动能转化N=ΔT2因此，风作用在风轮叶片上的力F和风轮输出的功率N分别为2N=4ps(v1-v2)(v1+v2)令其等于0，求解方程，得1求Nmax得令=0.593为Cp,称作贝茨功率系数（或称作理想风能利用系数）,有而psv1正是风速为v1的风能,故N=CTmaxppC=0.593说明风吹在叶片上,叶片所能获得的最p大功率为风吹过叶片扫掠面积风能的59.3%.贝茨理论说明,理想的风能对风轮叶片做功的最高效率通常风力机风轮叶片接受风能的效率达不到59.3%，一般设计时根据叶片的数量、叶片翼型、功率等情况，取0.25～0.45。风轮空气动力学的几何定义风轮轴：风轮旋转运动的轴线。旋转平面：与风轮轴垂直，叶片在其旋转的平面。风轮直径：风轮扫掠面的直径。叶片轴：叶片纵向轴，绕此轴可以改变叶片相对于旋转平面的偏转角（安装角）。在半径r处的叶片截面：叶片与半径r并以风轮轴为轴线的圆柱相交的截面安装角或桨距角：在半径r0处翼型的弦线与旋转面二、葛劳渥旋涡理论旋涡理论的优点在于考虑了通过叶轮的气流诱导转动。叶轮旋转工作时，流场并不是简单的一维定常流动，而是一个三维流场，在流场中会形成一种是由于气流流经旋转的叶轮，通过叶片尖部的气流迹线为螺旋线，在流场中形成螺旋涡流；同样在轮毂附近有同样的旋涡形成中心涡流；另外，气流通过叶片时，由于叶片表面上下压力不同，也形成涡流，这个涡流叫边界涡流。正因即引入诱导因子。轴向方向上，由于气流旋涡运动，在高速时形成一个平面阻碍气流在轴向方向上运动，引入轴向诱导因子a。周向方向上，由于气流旋涡运动，气流在下游周周向方向上，由于气流旋涡运动，气流在下游周向方向上产生一个旋转角速度z，上游周向方向的角速度为0，引入周向诱导因子a’。而假定叶轮以论的思想可得出：气流在叶轮处的角速度为(z+0)/2，则在叶轮平面内气流对叶轮的角速度为ω′w=w+z2Z=(a'-1)w气流对叶轮的角速度(2)三、动量理论-P+PVV∞V2 -P+PVV∞V2 图9动量理论简图在叶轮上r--r+dr的环域内应用动量定理（如图9），则风作用在该环域上的轴向推力为dF=dm(V∞-V2)=2πrdrV(V∞-V2)(3)式中：V∞—来流风速，m/s；V—叶轮处的风速，m/s；V2—叶轮后很远处的风速，m/s。由于叶轮前后有压力差ΔP=P+-P-，则轴向力又可用下式表示：dF=2πrdr(P+-P-)(4)利用伯努利方程得P∞+1/2ρV∞2=P++1/2ρV2P∞+1/2ρV22=P-+1/2ρV2以上二式相减并代入式（4），则dF=ρπrdr(V∞2-V22)(5)联立(3)、(5)两式，可得V=（V∞+V2）/2(6)因为引入轴向干扰因子，则将上式分别代入式(6)，式(5)得2dF=2)作用在环形面上的转距作用在环形面上的转距dM利用角冲量定理，由于dm—气流通过叶轮时质量的变化d将式(7)代入式(9)得(10)将式(2)代入式(10)得∞(11)dr(1+a)(a'-1)叶素理论是将叶片沿展向分成几个微段（一般化分为十个微段），每个微段称为一个叶素。这里假设作用在每个叶素上的力相互之间没有干扰，叶素本身可以视为一个二元翼型。研究叶片的受力情况，一般以叶素为研究对象，分析叶素上所受的力和力矩，然后沿翼展方向上积分，即可求得叶片上所受的力和力矩。V图10在半径r处的叶片截面图11安装角流中的受力分根据二元翼型理论，作用在翼型上的升力和阻力为式中：式中：dFL—翼型升力，方向与气流相对叶轮的相对速度dFD—翼型阻力，方向与气流相对叶轮的相对速度W—气流对叶素的相对速度，m/s；ρ—空气密度，kg/m3；C—距转轴r处的翼型弦长，m；Cl、Cd—分别为升力系数和阻力系数，Cl、Cd值由所选翼型决定。dFy=dFLsinI-dFDcosI=pCW2dr(ClsinI-CdcosI)(4)将上式代入式(3)、式(4)并消去Cd整理得，所以Cd=Cltgε那么，作用在（r，r+dr）环域上的叶素所受的轴式中：B—叶片数。转矩dM=BdFyr=）式中：令从流经环形面积（r,r+dr）的气流中获得的最大3+a)(a'1)（13）（14）ε=0，代入式得由式(9)由式(9)得将以上两式联立可推出把式(15)代入式(14)得对于给定的λ值，当，dCp/da0时d2Cp/da2<0令，cosθ式中：θ-所采用的中间变量(17)33将式(17)代入式(16)，并在等式两边同除以(λ2+1)24cos3θ-3cosθ=cos3θ所以cos3θ=-可见当θ取到一定值时,叶片不同截面处满足风能系数达到最大值,利用式(17)、式(15)、式(11)、式(7)可相应得出在不同截面处的周向诱导因子,轴向诱导因子,叶弦弦长,气相角等几何参数。λ,λ=λ0；(3)由式(18)求