风力机的基础理论.docx
- 文档编号:24024877
- 上传时间:2023-05-23
- 格式:DOCX
- 页数:21
- 大小:231.60KB
风力机的基础理论.docx
《风力机的基础理论.docx》由会员分享,可在线阅读,更多相关《风力机的基础理论.docx(21页珍藏版)》请在冰豆网上搜索。
风力机的基础理论
第二章风力机的基础理论[3、4]
第一节风力机的能量转换过程
一、风能的计算
由流体力学可知,气流的动能为
(2-1)
式中m──气体的质量;
v──气体的速度。
设单位时间内气流流过截面积为S的气体的体积为L,则
L=Sv
如果以ρ表示空气密度,该体积的空气质量为
m=ρL=ρSv
这时气流所具有的动能为
(2-2)
上式即为风能的表达式。
在国际单位制中,ρ的单位是kg/m3;L的单位是m3;v的单位是m/s;E的单位是W。
从风能公式可以看出,风能的大小与气流密度和通过的面积成正比,与气流速度的立方成正比。
其中ρ和v随地理位置、海拔高度、地形等因素而变。
二、自由流场中的风轮
风力机的第一个气动理论是由德国的Betz于1926年建立的。
Betz假定风轮是理想的,即它没有轮毂,具有无限多的叶片,气流通过风轮时没有阻力;此外,假定气流经过整个叶轮扫掠面时是均匀的;并且,气流通过风轮前后的速度为轴向方向。
现研究一理想风轮在流动的大气中的情况(见图2-1),并规定:
v1──距离风力机一定距离的上游风速;
v──通过风轮时的实际风速;
v2──离风轮远处的下游风速。
设通过风轮的气流其上游截面为S1,下游截面为S2。
由于风轮的机械能图2-1叶轮的气流图
量仅由空气的动能降低所致,因而v2必然低于v1,所以通过风轮的气流截面积从上游至下游是增加的,即S2大于S1。
如果假定空气是不可压缩的,由连续条件可得:
S1v1=Sv=S2v2
风作用在风轮上的力可由Euler理论写出:
F=ρSv(v1-v2)(2-3)
故风轮吸收的功率为
(2-4)
此功率是由动能转换而来的。
从上游至下游动能的变化为
(2-5)
令式(2-4)与式(2-5)相等,得到
(2-6)
作用在风轮上的力和提供的功率可写为:
(2-7)
(2-8)
对于给定的上游速度v1,可写出以v2为函数的功率变化关系,将式(2-8)微分得
式
有两个解:
v2=-v1,没有物理意义;
v2=v1/3,对应于最大功率。
以
代入P的表达式,得到最大功率为:
(2-9)
将上式除以气流通过扫掠面S时风所具有的动能,可推得风力机的理论最大效率(或称理论风能利用系数)
(2-10)
式(2-10)即为有名的贝兹(Betz)理论的极限值。
它说明,风力机从自然风中所能索取的能量是有限的,其功率损失部分可以解释为留在尾流中的旋转动能。
能量的转换将导致功率的下降,它随所采用的风力机和发电机的型式而异,因此,风力机的实际风能利用系数Cp<0.593。
风力机实际能得到的有用功率输出是
(2-11)
对于每m2扫风面积则有:
(2-12)
三、风力机的特性系数
在讨论风力机的能量转换与控制时,以下特性系数具有特别重要的意义。
1.风能利用系数Cp
风力机从自然风能中吸取能量的大小程度用风能利用率系数Cp表示,由(2-11)知
(2-13)
式中P——风力机实际获得的轴功率,W;
ρ——空气密度,kg/m3
S——风轮的扫风面积,m2
v——上游风速,m/s
2.叶尖速比(TipSpeedRatio)
为了表示风轮在不同风速中的状态,用叶片的叶尖圆周速度与风速之比来衡量,称为叶尖速比λ
(2-14)
式中n——风轮的转速,r/s;
ω——风轮角频率,rad/s
R——风轮半径,m;
v——上游风速,m/s
3.扭矩系数CT和推力系数CF
为了便于把气流作用下风力机所产生的扭矩和推力进行比较,常以λ为变量作成扭矩和推理的变化曲线。
因此,扭矩和推力也要无因次化。
式中,T——扭矩,N·m;
F——推力,N;
第二节桨叶的几何参数和空气动力特性
无论风力机的型式如何,桨叶是其至关重要的部件。
为了很好地理解它在控制能量转换中的作用,必须知道某些空气动力学的基本知识。
先研究一静止的叶片,其承受的风速为v,假定风速方向与叶片横截面平行。
一、叶型的几何参数和气流角
图2-2叶型的几何参数和气流角
叶型的几何定义:
B点:
后缘(Trailingedge);
A点:
前缘(Leadingedge),它是距后缘最远的点;
l:
叶型的弦长,是两端点A、B连线方向上叶型的最大长度;
C:
最大厚度,即弦长法线方向之叶型最大厚度;
:
叶型相对厚度,
,通常为10%~15%;
叶型中线:
从前缘点开始,与上、下表面相切的诸圆之圆心的连线,一般为曲线;
f:
叶型中线最大弯度;
:
叶型相对弯度,
;
i:
攻角,是来流速度
与弦线间的夹角;
θ0:
零升力角,它是弦线与零升力线间的夹角;
θ:
升力角,来流速度方向与零升力线间的夹角。
i=θ+θ0(2-13)
此处θ0是负值,θ和i是正值。
二、作用在运动桨叶上的气动力
假定桨叶处于静止状态,令空气以相同的相对速度吹向叶片时,作用在桨叶上的气动力将不改变其大小。
气动力只取决于相对速度和攻角的大小。
因此,为便于研究,均假定桨叶静止处于均匀来流速度v中。
此时,作用在桨叶表面上的空气压力是不均匀的,上表面压力减少,下表面压力增加。
按照伯努里理论,桨叶上表面的气流速度较高,下表面的气流速度则比来流低。
因此,围绕桨叶的流动可看成由两个不同的流动组合而成;一个是将叶型置于均匀流场中时围绕桨叶的零升力流动,另一个是空气环绕桨叶表面的流动。
而桨叶升力则由于在桨叶表面上存在一速度环量如图(2-3)所示。
图2-3气流绕翼叶的流动
为了表示压力沿表面的变化,可作桨叶表面的垂线,用垂线的长度KP表示各部分压力的大小,
(2-14)
式中p──桨叶表面上的静压;
(ρ、po、v)──无限远处的来流条件。
连接各垂直线段长度KP的端点,得到图2-4a,其中上表面KP为负,下表面KP为正。
图2-4作用在翼叶上的力
使用在桨叶上的力F与相对速度的方向有关,并可用下式表示
(2-15)
式中S──桨叶面积,等于弦长×桨叶长度;
Cr──总的气动系数。
该力可分为两部分:
分量Fd与速度v平行,称为阻力;分量Fl与速度v垂直,称为升力。
Fd与Fl可分别表示为:
(2-16)
(2-17)
式中Cd──阻力系数;
Cl──升力系数。
因两个分量是垂直的,故可写成
若令M为相对于前缘点的由F力引起的力矩,则可求得变距力矩系数
。
(2-18)
式中
──弦长。
因此,作用在桨叶截面上的气动力可表示为升力、阻力和变距力矩三部分。
由图2-4b可看出,对于各个攻角值,存在某一特别的点C,该点的气动力矩为零,称为压力中心。
于是,作用在叶型截面上的气动力可表示为作用在压力中心上的升力和阻力。
压力中心与前缘点之间的位置可用比值CP确定。
(2-19)
一般CP=(25~30)%。
三、升力和阻力系数的变化曲线
(一)C1和Cd随攻角的变化
首先研究升力系数的变化,它由
直线和曲线两部分组成。
与
对应的
点称为失速点,超过失
速点后,升力系数下降,阻力系
数迅速增加。
负攻角时,
也呈
曲线形,
通过一最低点
。
图2-5桨叶的升力和阻力系数
阻力系数曲线的变化则不同,它的最小值对应一确定的攻角值。
阻力系数曲线的变化则不同,它的最小值对应一确定的攻角值。
不同的桨型截面形状对升力和阻力的影响很大,现分述如下:
1、弯度的影响
叶型的弯度加大后,导致上、下弧流速差加大,从而使压力差加大,故升力增加;与此同时,上弧流速加大,摩擦阻力上升,并且由于迎流面积加大,故压差阻力也加大,导致阻力上升。
因此,同一攻角时,随着弯度增加,其升、阻力都将显著增加,但阻力比升力的增加更快,使升阻比有所下降。
2、厚度的影响
叶型厚度增加后,其影响与弯度类似。
同一弯度的叶型,采用较厚的桨型时,对应于同一攻角的升力有所提高,但对应于同一升力的阻力也较大,使升、阻比有所下降。
3、前缘的影响
试验表明,当桨型的前缘抬高时,在负攻角情况下阻力变化不大。
前缘低垂时,则在负攻角时会导致阻力迅速增加。
4、表面粗糙度和雷诺数的影响
表面粗糙度和雷诺数对桨叶空气动力特性有着重要影响。
图2-6表示雷诺数和表面粗糙度对几种桨型(NACA0012,NACA23012、23015、NACA4412、4415)的气动力特性的影响曲线。
当叶片在运行中出现失速以后,噪声常常会突然增加,引起风力机的振动和运行不稳等现象。
因此,在选取Cl值时,以失速点作为设计点是不好的。
对于水平轴型风力机而言,为了使风力机在稍向设计点右侧偏移时仍能很好的工作,所取的Cl值,最大不超过(0.8-0.9)Clmax。
图2-6雷诺数和表面粗糙度对气动特性的影响
(二)埃菲尔极线(EiffelPolar)
为了便于研究问题,可将Cl和Cd表示成对应的变化关系,称为埃菲尔极线,见图2-7。
其中直线OM的斜率是:
tgθ=Cl/Cd
图2-7埃菲尔极线
四、有限叶展的影响
上述结果仅适用于桨叶无限长时,对于有限长度的叶片,其结果必须修正。
由于升力桨的下表表面压力大于大气,上表面低于大气,因此叶片两端气流企图从高压侧向低压侧流动,结果在两端形成涡流。
实际上,由于叶尖的影响,两端形成一系列的小涡流,这些小涡流又汇合成两个大涡流,卷向桨尖内侧(图2-8)
涡流形成的后果,造成阻力增加,引起一诱导阻力:
(2-20)
因此,上述阻力系数变为:
(2-21)
式中,Cdo──无限叶展的阻力系数。
图2-8有限叶展形成的涡流
为此,要想得到同样的升力,攻角必须增加一个量φ,故获得同样升力的新攻角为:
i=i0+φ(2-22)
由流体力学知,当环量呈椭圆分布时,Cdi和φ可由下列关系给出:
(2-23)
式中S──桨叶面积;
L──桨叶长度;
α──是展弦比,α=L2/S。
五、弦线和法线方向的气动力
如果将F力分解为弦线方向和垂直于弦线方向的两个分量(图2-9,a)则有:
弦线方向:
垂线方向:
上式可进一步写成:
(2-24)
式中
Cn=Clcosi﹢Cdsini(2-25)
Ct=Cdcosi﹣Clsini
Cn与Ct对应的曲线示于图2-9b,称为利兰热尔极线(LilienthalPolar)。
图2-9弦线和法线方向的气动力
第三节风轮的气动力学
一、几何定义
为了研究风力机的叶轮,先给出一些定义:
转轴:
风轮的旋转轴;
回转平面:
垂直于转轴线的平面,叶片在该平面内旋转;
风轮直径:
风轮扫掠面直径;
叶片轴线:
叶片纵轴线,围绕它,可使叶片一部分或全部相对于回转平面的倾斜变化;
安装角或节距角β:
半径r处回转平面与叶型截面弦长之间的夹角。
二、叶素特性分析
取一长度为dr的叶素,在半径r处的弦长为l,节距角为β。
则叶素在旋转平面内具有一圆周速度U=2πrn,。
如果取v为吹过风轮的轴向风速,气流相对于叶片的速度为w(图2-10b),则
,
(2-26)
而攻角为
。
其中,I为w与旋转平面间的夹角,称为倾斜角。
图2-10叶素特性分析
因此,叶素受到相对速度w的气流作用,进而受到一气动力dF作用。
dF可分为一个升力dL和一个阻力dD,分别与相对速度w垂直或平行(图2-10b),并对应于某一攻角i。
Cl和Cd的值可按相应的攻角查取所选叶型的气动特性曲线得到。
现在来计算由气动力dF产生的作用在风轮上的轴向推力以及作用在转轴上的力矩。
dFa为dF在转轴上的投影,dT为dF在回转平面上的投影对转轴的力矩,ω为风轮角速度,则有
dFa=dLcosI+dDsinI
dT=r(dLsinI-dDcosI)
代入以前的有关关系式得:
w2=v2+u2=v2+ω2r2,ωr=vctgI
dP=ωdT
于是得到dFa、dT和dP的下列表达式:
(2-27)
(2-28)
(2-29)
三、推力、扭矩和功率的一般关系式
风作用在风轮上引起的总推力Fa和作用在转轴上的总扭矩T可由所有作用在叶素上的dFa和dT求和得到。
推力Fa、扭矩T、功率P和效率η的关系式为
P=ΣdFa·v=Fav(2-30)
轴功率Pu=Tω(2-31)
效率为
(2-32)
第四节简化的风力机理论
一、基本关系的确定
为了确定叶片弦长,需要计算从转轴算起的r,r+dr一段截面上所受到的轴向推力。
可采用两种方法,该方法假定风力机按照贝兹(Betz)理论处于最佳状态下运转。
(一)第一种方法
按照贝兹理论,作用在整个风轮上的轴向推力[见公式(2-7)]为
通过风轮的风速[见公式(2-6)]是
此处v1和v2是离开风力机前后一定距离的风速。
当v2=v1/3时,功率输出达到最大值。
此时轴向推力Fa和通过扫风面的风速v是:
,
(2-33)
假设各单元扫风面产生的轴向推力正比于它的对应面积,则作用在间隔r,r+dr叶素的轴向力(图2-11)为:
(2-34)
(二)第二种方法
设角速度为ω,则半径为r的叶素圆周速度为u=ωr,此时风通过风轮的绝对速度v、相对于叶片的速度w和叶轮的圆周速度u三者的关系为v=w+u,并可写成w=v-u。
于是作用在叶片上长度为dr的升力和阻力为:
(2-35)
(2-36)
其合力为:
dF=dL/cosε
式中ε──dF和dL之间的夹角;
—半径r处的叶片弦长。
因w=v/sinI
则
(2-37)
图2-11作用在叶素上的力
将dF投影到转轴上,则r,r+dr段产生的轴向力dFa为:
(2-38)
式中b──叶片数。
令上式与式(2-24)相等,得到
(2-39)
二、上述关系的转换和简化
展开cos(I-ε)项,上述关系可写成
(2-40)
在最佳运转条件下,通过风轮的风速是
因此倾斜角可由下式确定
(2-41)
于是式(2-40)可进一步写成
(2-42)
在正常运行情况下,tgε=dD/dL=Cd/Cl一般是很小的,对于攻角在最佳值附近的普通桨型,其tgε约为0.02,则上式可简化为:
(2-43)
叶尖和半径r处的速度比分别为
及
,消去ω和v1后得到:
将λ值代入上式,有
(2-44)
叶尖速比λ0和风轮直径确定后,可由下式计算不同半径r处的倾斜角I:
(2-45)
三、与设计有关的几个问题
1、如果叶片安装角β确定了,则攻角也即确定(i=I-β),然后由叶型气动力特性曲线即可确定Cl。
当叶片数给定后,Clbl的表达式可用来确定以r为变量的各叶片截面的弦长;
2、公式(2-44)说明,给定半径r处的弦长,随叶尖速比λ0的增加而减少,因而处于高额定转速下工作的风轮,其重量也就较轻;
3、对于给定的叶尖速比λ0,弦长从叶尖向叶根增加,这一规律使叶片形成曲线边缘。
但在某些风力机中Cl沿叶片长度并不保持为常数,因而弦长也不需要从叶尖向叶根递增。
四、叶素的理论气动效率和最佳攻角
r,r+dr段叶片的气动效率可由下式确定:
式中dFu──气动力dF在旋转平面上的投影值;
dFa──气动力dF在转轴上的投影值;
dPu──dr段叶片产生的风轮功率;
dPt──流过dr段叶片的风的功率。
图2-12最佳攻角位置的确定
因dFu=dLsinI-dDcosI
dFa=dLcosI﹢dDsinI
ctgI=u/v
则得到
(2-46)
若以
代入上式,则
(2-47)
当tgε较低时,效率是较高的。
如果在tgε等于零的极限情况下,气动效率将等于1。
实际的tgε值取决于攻角的大小。
当直线OM与埃菲尔极线相切时(图2-12),与该点对应的攻角使得tgε成为最小,在这个特定的攻角时,气动效率达到最大值。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 风力机 基础理论