风力机空气动力学-第四章.ppt
《风力机空气动力学-第四章.ppt》由会员分享,可在线阅读,更多相关《风力机空气动力学-第四章.ppt(124页珍藏版)》请在三一办公上搜索。
1、华北电力大学,风力机空气动力学,1,风力机空气动力学,华北电力大学,风力机空气动力学,2,主要内容,概述基本理论风力机的空气动力设计风力机性能风电场中的空气动力学问题计算流体力学在风力机和风电场分析和设计中的应用,华北电力大学,风力机空气动力学,3,4:风力机空气动力学 4-1:概述,研究背景能源问题风能:人类最古老能源新能源、可再生能源我国丰富的风资源与 政府的大力支持风能是有很强综合性的技术学科,华北电力大学,风力机空气动力学,4,4-1:概述,研究内容风力机空气动力模型;风力机翼型空气动力特性;风力机叶片空气动力设计;风力机风轮性能计算;风力机空气动力载荷计算;风力机气动弹性稳定性和动力
2、响应;风力机空气动力噪声和风力机在风电场中的布置等。研究方法理论计算:工程计算方法和数值计算方法风洞试验风场测试,华北电力大学,风力机空气动力学,5,4-1:概述,华北电力大学,风力机空气动力学,6,4-1:概述,叶片,塔架,机舱,轮毂,尾舵调向/风向标,低速轴增速器高速轴发电机停车制动器,华北电力大学,风力机空气动力学,7,4-1:概述,华北电力大学,风力机空气动力学,8,4:风力机空气动力学 4-2:基础理论,动量理论尾流不旋转的动量理论风轮尾流旋转时的动量理论叶素理论动量叶素理论涡流理论,华北电力大学,风力机空气动力学,9,4-2:基础理论,动量理论用来描述作用在风轮上的力与来流速度之间
3、的关系,回答风轮究竟能从风的动能中转换成多少机械能。风轮尾流不旋转的动量理论 研究不考虑风轮尾流旋转时的理想情况假设:气流是不可压缩的均匀定常流;风轮简化成一个桨盘;桨盘上没有摩擦力;风轮流动模型简化成一个单元流管;风轮前后远方的气流静压相等;轴向力(推力)沿桨盘均匀分布。图-风轮流动的单元流管模型,华北电力大学,风力机空气动力学,10,3-2:基础理论,假设来自远前方的流管,在叶轮激盘处恰与激盘外径相切,并伸向下游,如此建立控制体。应用一维动量方程得 单位时间流经风轮的空气 风轮处的质量流量:,那么,激盘前后区域应用伯努利方程由假设知,则根据动量理论得,华北电力大学,风力机空气动力学,11,
4、4-2:基础理论,由于受风轮的阻挡,风流向风轮时速度减小。定义轴向诱导因子 和轴向诱导速度 风轮尾流处的轴向诱导速度是风轮处的二倍。如果风轮吸收风的全部能量,即 而实际情况下,风轮仅能吸收部分能量,因此。,华北电力大学,风力机空气动力学,12,4-2:基础理论,根据能量方程,风轮吸收的能量(风轮轴功率P)等于风轮前后气流动能之差(据假设流动前后静压不变)代入 得出,华北电力大学,风力机空气动力学,13,4-2:基础理论,定义风轮功率系数/风轮风能利用系数为因此,当 时,风轮功率系数最大:贝兹(Betz)极限 即在理想情况下,风轮最多能吸收59.3%的风的动能。对应于最大值,有,华北电力大学,风
5、力机空气动力学,14,4-2:基础理论,风轮附近速度和压力的变化规律风力平面处的风速总比来流小(风轮吸收了功率)本模型假设尾迹不旋转,意味着在转动尾迹的动能中没有能量损失。实际上肯定是有损失的。即使对最佳设计的风轮也不可能系数60的风动能。,华北电力大学,风力机空气动力学,15,4-2:基础理论,为什么风力机尾迹的流管是扩张的?由质量守恒方程以及可得对于最大功率情况,有,华北电力大学,风力机空气动力学,16,4-2:基础理论,实验表明,由前面理想风力机假设所得到的功率系数和推力系数只在约a0.4时是正确的。大于0.4,风轮前后的速度差变大,需要从外部把动量输入到尾迹中,使尾迹边沿的自由剪切层不
6、稳定,形成湍流尾迹状态。,华北电力大学,风力机空气动力学,17,4-2:基础理论,风轮尾流旋转时的动量理论由于风施加在风轮上的力矩的反作用,使转子后面的流动以与转子相反的方向旋转尾迹的旋转将减少风轮对能量的吸收。一般旋转尾迹的这部分动能将随转子力矩的增大而增加。所以,低转速风轮(小转速、大转矩)要比高转速(低转矩)产生大的尾迹旋转损失。(功率不变),华北电力大学,风力机空气动力学,18,4-2:基础理论,与前面比较,本节考虑风轮尾流的旋转。气流在风轮上产生转矩时,也受到风轮的反作用力,由此气流产生了一个反向的角速度,使尾流以相反的方向转动。即、由于流体的粘性,激盘诱导了流动的旋转,导致激盘诱导
7、的速度沿激盘径向不是常数,或诱导因子a是变化的。同时,由于激盘的转动,还会对流体产生周向的诱导速度,以及转动力矩,。如果在风轮处气流的角速度和风轮的角速度相比是个小量的话,那么一维动量方程仍可应用,仍假设风轮前后远方的气流静压相等。取控制体如图,华北电力大学,风力机空气动力学,19,4-2:基础理论,应用动量方程,作用在风轮平面圆环上的轴向力/推力为:单位时间流经风轮平面 圆环上的空气质量流量,即:风轮平面 圆环的面积,由轴向诱导因子,诱导因子a是随半径变化的!,华北电力大学,风力机空气动力学,20,4-2:基础理论,另一个求推力的方法如果采用一个控制体,它以激盘(叶轮)的角速度旋转W,激盘后
8、面气流相对与叶片的角速度为W w。把伯努利方程用于激盘的前后截面,来推导穿过叶轮的压力差:设叶轮前后的截面分别为a和b,有 并采用了前面的假设,即穿过叶轮的轴向速度V相同。简化处理后,可得:如果引入下面的诱导因子bw/2W,上式成为:,华北电力大学,风力机空气动力学,21,4-2:基础理论,对控制体应用动量矩方程,则作用在风轮平面 圆环上的转矩可以表示为 其中:风轮叶片 处的周向诱导速度:风轮叶片 处轴向诱导角速度,定义周向诱导因子,其中 为风轮转动角速度结合,华北电力大学,风力机空气动力学,22,4-2:基础理论,风轮轴功率是转矩与风轮角速度的乘积,因此引入风轮叶尖速比,风轮扫掠面积,得则风
9、轮功率系数 可表示为或其中 为当地的速比为了求解这个式子,需要知道b,a随r或lr的变化关系,华北电力大学,风力机空气动力学,23,4-2:基础理论,通过使前述的由两种方法推导的推力相等,可以求得或由前面功率的表达式可知,当b(1-a)取最大值时,是产生最大功率的空气动力条件。把上式代入并消去b,然后对所得到的式子进行求导,置零。就可得到在最大功率条件下,轴向诱导因子与当地速比的关系:把该式代入到第一个是在,得,华北电力大学,风力机空气动力学,24,4-2:基础理论,对前页第三个式子进行微分,得把该式代入到功率系数得表达式,得其中下限对应于lr0时得轴向诱导因子,上限对应于lr l时得诱导因子
10、。由下式知,a的最小值为0.25,最大值为1/3.积分后得(x1-3a),华北电力大学,风力机空气动力学,25,4-2:基础理论,尖速比越大,理论最大功率系数越大当考虑风轮后尾流旋转时,风轮轴功率有损失,风轮功率系数要减小。在轮毂附近以外,a就接近其理想值1/3在转子中叶高以上,b接近于0。,b,b,b,a和b随半径的变化,尖速比7.5,华北电力大学,风力机空气动力学,26,人们已经习惯地认为一个风力机只能处于我们所期望的状态:吸收风的动能并转换成轴功率,伴随着风速的减小。实际上其它状态也可能出现,设计者应该从物理上理解其它状态,以及它们对风轮载荷、动力学、和速度控制的重要影响。水平轴风力机的
11、叶片几乎都采用带弯度的翼型,其它状态的发生或转变与翼型的零升力线有关,而不是弦线。,零升力几何定义,华北电力大学,风力机空气动力学,27,考虑来流速度和转子的转速恒定,叶片浆矩角变化:如果浆矩角q比气流角f大,(a为负),转子处于螺旋桨状态,使气流加速,必须给转子提供能量,以维持转速,a小于零。如果浆矩角q与气流角f相对,(a为零),转子处于零滑移状态,a为零如果浆矩角q小于气流角f,(a大于零),转子处于风轮状态,使气流加速,转子吸收风中的能量,01.0),华北电力大学,风力机空气动力学,28,4-2:基础理论,前面讲述了如何应用物理学定律确定绕流风力机的流动特性,以及从风中所能获得的最大功
12、率的极限。下面介绍应用翼型获得逼近这个理论上可吸收的功率。叶素理论基本出发点叶素:风轮叶片沿展向分成许多微段,假设微段间流动相互没有干扰,即可以视为二维翼型。叶片上的力只有升力和阻力将作用在每个叶素上的力和力矩沿展向积分,即可求得作用在风轮上的力和力矩。,华北电力大学,风力机空气动力学,29,叶素上气流速度三角形和空气动力分量图:入流角(气流角):迎角:几何扭角(浆矩角)Vx0:垂直于风轮旋转平面速度分量 Vy0:垂直于风轮旋转平面速度分量,华北电力大学,风力机空气动力学,30,4-2:基础理论,叶素上的升力和阻力必须与有效的相对速度垂直或平行。有效的相对速度为风力机前的轴向有效风速V1(1-
13、a),旋转速度分量Wr(1+b)构成的矢量。推导过程根据动量理论,考虑尾流旋转时 即叶素处的入流角和迎角可表示为如此,求出迎角后,即可根据翼型空气动力特性曲线得到叶素的升力系数Cl 和阻力系数Cd。,华北电力大学,风力机空气动力学,31,4-2:基础理论,合成气流速度V0引起的作用在长度为dr 叶素上的空气动力dFa可以分解为法向力dFn和切向力dFt,则 其中c:叶素剖面弦长 Cn、Ct:法向力系数和切向力系数 则 这时,作用在风轮平面dr圆环上的轴向力(推力)可表示为 其中B:叶片数,作用在风轮平面dr圆环上的转矩为,阻力使切向力,即力矩减小,而使推力增加,华北电力大学,风力机空气动力学,
14、32,4-2:基础理论,动量叶素理论需要计算风轮旋转面中的轴向诱导因子a和周向诱导因子b,以便计算作用在风轮叶片上的力和力矩。这就是所谓的动量叶素理论(条带理论或叶素动量理论(BEM)推导过程 根据叶素理论知 得出,华北电力大学,风力机空气动力学,33,4-2:基础理论,动量理论 叶素理论,结合,结合,当地实度,华北电力大学,风力机空气动力学,34,4-2:基础理论,如果考虑普朗特叶尖损失修正因子(Prandtl L.1957)上结论可表示为根据上面的关系式可以通过迭代方法求得轴向诱导因子a和周向诱导因子b,从而可得到气流角,然后就可计算处对应的力和力矩系数,沿展向积分得到总的力和力矩系数,以
15、及功率。计算a和b的迭代步骤为:,华北电力大学,风力机空气动力学,35,4-2:基础理论,假设a和b的初值,一般可取0;计算入流角计算迎角=-根据翼型空气动力特性曲线得到叶素的升力系数Cl和阻力系数Cd计算叶素的法向力系数Cn和切向力系数Ct:计算a和b的新值:比较a、b新值与初值,如果误差小于设定的误差值(一般可取0.001),则;若否,则重新假设初值,返回继续迭代。迭代终止。,华北电力大学,风力机空气动力学,36,4-2:基础理论,前面所有结果的条件是:风沿着风轮的轴向且均匀。如果风向偏斜、风剪切、垂直风分量和叶片有锥角等时,动量方程不再适用;就要采用一些修正。对于风向偏斜的情况,可用下面
16、的经验公式对动量叶素理论进行修正。威尔森(Wilson)修正方法(Wilson R.E.1976)当a0.38时,第步中由 代替 葛劳渥特(Glarert)修正方法(Glauert H.1935)当a0.2时,第步中由 代替,其中,华北电力大学,风力机空气动力学,37,4-2:基础理论,这里计算上述轴向诱导因子a和周向诱导因子b时,都假设风轮的锥角 为零。当风轮的锥角不为零时 需化为如果考虑普朗特叶尖损失修正因子,则用 其中,华北电力大学,风力机空气动力学,38,4-2:基础理论,涡流理论 叶片静止时,据赫姆霍兹定理,叶片附着涡和后缘尾涡组成马蹄涡系。简化后,将叶片分成无限多沿展向宽度很小的微
17、段。若假设每个微段上的环量沿展向是个常量,则可用在每个微段上布置的马蹄涡系来代替风轮叶片。若考虑环量沿弦向的变化,认为每个微段上的马蹄涡系由许多个等强度马蹄涡组成。沿叶片展向,每个微段马蹄涡系的附着涡总强度等于绕该微段叶片的环量;而从后缘拖出的尾涡强度是相邻两微段叶片环量之差(每个微段的尾涡都与相邻微段的尾涡重合,且方向相反)。,华北电力大学,风力机空气动力学,39,4-2:基础理论,对长度较大的风轮叶片,则可以简化成用一个位于1/4弦线变强度 的附着涡线和从附着涡向下游拖出的尾涡系来代替。尾涡系由许多个与轴线平行的直涡线所组成。当叶片旋转时,从后缘拖出的尾涡系将变成一个由螺旋形涡面组成的复杂
18、涡系。而且随着涡与涡之间的相互干扰,该涡系不断变形。图为一个两叶片水平轴风力机叶片旋转时沿展向等环量分布的涡系。它由附着涡、叶尖螺旋形自由涡和叶根中心涡三部份组成。,华北电力大学,风力机空气动力学,40,4-2:基础理论,为了预测风力机风轮的性能,针对尾涡系又发展了不同的尾涡模型。1)刚性尾涡模型 模型假设:叶片数无限多,且实度一定,从而叶片尖部后缘拖出的尾涡形成一个管状的螺旋形涡面。假设该管状涡直径不变,即形成一个圆柱状的螺旋形涡面,又称为柱涡。叶片根部接近风轮旋 转轴,从旋转叶片根 部后缘拖出的尾涡认 为形成一个绕风轮旋 转轴旋转的中心涡。,华北电力大学,风力机空气动力学,41,4-2:基
19、础理论,模型结果:转矩功率功率系数与不考虑尾流的动量理论 相比,由于风轮尾流旋转需要消耗一部分能量,来平衡旋转流动产生的离心力所引起的压力梯度而造成的静压损失。由于风轮旋转时,流经风轮旋转面各点的速度是不均匀的,因此,风力机叶片附着涡的强度沿叶片展向和叶片方位角都有变化。于是,刚性尾涡模型不足以完全描述尾涡的几何形态,需要发展半刚性尾涡模型和自由尾涡模型。,华北电力大学,风力机空气动力学,42,4-2:基础理论,2)半刚性尾涡模型模型的两个组成部分:附着环量沿展向变化产生的尾随涡;附着环量沿方位角变化产生的脱落尾涡。模型的三个组成区域:近尾涡区、中间尾涡区和远尾涡区,它们随旋转频率呈周期性变化
20、。,华北电力大学,风力机空气动力学,43,4-2:基础理论,3)两种尾涡模型比较 图为采用两种尾涡模型计算的轴向诱导因子沿展向的分布曲线。可知半刚性尾涡模型计算的轴向诱导因子值沿叶片展向是变化的,趋于合理。,华北电力大学,风力机空气动力学,44,4:风力机空气动力学 4-3:风力机空气动力设计,风力机几何参数叶片几何参数风轮几何参数风力机空气动力设计参数风力机翼型风力机叶片气动外形设计,华北电力大学,风力机空气动力学,45,4-3:风力机空气动力设计,风力机几何参数 叶片几何参数 风轮叶片的平面形状一般为梯形,叶片主要几何参数为:叶片长度叶片展向方向上的最大长度,用L表示。叶片弦长叶片各剖面处
21、翼型的弦长,用c来表示。叶片弦长沿展向变化,叶片根部剖面的翼弦称翼根弦,用cr表示,叶片梢部剖面的翼弦称翼梢弦,用ct表示。,华北电力大学,风力机空气动力学,46,4-3:风力机空气动力设计,叶片面积通常指叶片无扭角时在风轮旋转平面上的投影面积叶片平均几何弦长叶片面积与叶片长度的比值,即叶片桨距角通常指叶片尖部剖面的翼弦与旋转平面之间的夹角,用表示;叶片各剖面的桨距角是叶片尖部剖面的桨距角与叶片各剖面的几何扭角之和。叶片扭角通常指叶片的几何扭角。它是叶片尖部桨矩角为零情况下,叶片各剖面的翼弦与风轮旋转平面之间的夹角。叶片扭角沿展向变化,叶片梢部的扭角比根部小。叶片转轴通常位于叶片各剖面的0.2
22、50.35翼弦处,与各剖面气动中心的连线重合或尽量接近,以减少作用在转轴上的转矩。,华北电力大学,风力机空气动力学,47,4-3:风力机空气动力设计,风轮几何参数 风轮由叶片和轮毂组成,主要几何参数为:风轮叶片数组成风轮的叶片个数,用B表示。风轮直径风轮旋转时的风轮外圆直径,用D表示。风轮面积 一般指风轮扫掠面积。风轮锥角指叶片与旋转轴垂直的平面的夹角,用表示。其作用是在风轮运行状态下,减少离心力引起的叶片弯曲应力以及防止叶片梢部与塔架碰撞。,华北电力大学,风力机空气动力学,48,4-3:风力机空气动力设计,风轮仰角风轮旋转轴与水平面的夹角,用表示。其作用是防止叶片梢部与塔架碰撞。风轮偏航角来
23、流速度矢在水平面上的分量与通过风轮旋转轴的铅垂面的夹角,用表示。风轮实度 风轮叶片面积与风轮扫掠面积的比值,风轮高度风轮高度是指风轮轮毂中心的离地高度,用Hh表示。,华北电力大学,风力机空气动力学,49,4-3:风力机空气动力设计,风力机的特性曲线风力机几何参数风力机空气动力设计参数叶片数风轮直径额定风速叶尖速比风轮转速塔架高度风力机翼型风力机叶片气动外形设计,华北电力大学,风力机空气动力学,50,风力机的特性曲线一个风力机可主要由三个量:功率、扭矩和推力表示功率确定风轮所能吸收的能量;扭矩确定了齿轮箱的尺度,并与驱动发电机所需的力矩相吻合;推力要影响塔架的结构设计。通常把它们表示成无量纲的形
24、式如果假设转子叶片的空气动力性能不变,那么风轮的空气动力性能就只与尖速比有关。,华北电力大学,风力机空气动力学,51,Cpl曲线例如,一个三叶片风轮最大的Cp值只有0.47(尖速比7).这是由于阻力和叶尖损失等在小尖速比时,小的Cp是由于失速,损失对性能的影响,华北电力大学,风力机空气动力学,52,叶轮实度的影响小的实度在大的尖速比区间内产生宽扩平坦 Cp,但最大Cp值大的实度产生窄的性能曲线,尖的峰值区使得风力机对尖速比的变化非常敏感。如果实度太大,最大Cp值反而会减少。最优的实度显然是三叶片风轮,华北电力大学,风力机空气动力学,53,CQl曲线(扭矩)用Cp除以转速可以得到;对带齿轮箱的机
25、组有用。扭矩随着实度的增加而增加对于先进的高速风力机,为了减少齿轮箱的花费,期望尽可能地降低扭矩.由失速引起的峰值处的尖速比,比功率曲线的峰值处的尖速比较小。,华北电力大学,风力机空气动力学,54,CTl曲线(推力)一般地,作用在风轮上的的推力随着实度的增加而增加.,华北电力大学,风力机空气动力学,55,恒定转速运行与电网连接,需要控制转速恒定。CTl曲线表示了但风速不变时,功率随转速的变化。而KP1/l 表示了当强迫转速恒定时,功率随风速的变化。恒速运行的缺点:风力机的效率随着风速增加而下降。,华北电力大学,风力机空气动力学,56,风力机总效率机械效率,包括电效率,华北电力大学,风力机空气动
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 风力机 空气动力学 第四
链接地址:https://www.31ppt.com/p-6492091.html