②下翼面平缓,V3≈V1,使其几乎保持原来的大气压,即:
P3≈P1。
结论:
由于机翼上下表面所受的压力差,使得机翼得到向上的作用力——升力。
当风流经翼型时,不仅会受到上述升力,还会受到气流对它的推力,这样,叶片翼型将受到一个合力Q。
Q在垂直于来流方向的分量称为升力,而在平行于来流方向的分量D称为阻力。
合力Q对于前缘点的力矩M,称为气动俯仰力矩。
翼型上的升力表示为
L1w2cClz(4-1)
2
翼型上的阻力表示为
12
Dw2cCdz(4-2)
2d
气动俯仰力矩表示为
122
Mw2c2Cmz(4-3)
2
式中——空气的密度,单位kg/m3;——相对速度,单位m/s;——几何弦长,单位;
Cl——升力系数,无量纲;
Cd——阻力系数,无量纲;
Cm——俯仰力矩系数,无量纲;——翼型在叶片长度方向上的长度。
(其中1w2的意义,对风功率密度1w3进行分析,22
1w31w2wFw,即F1w2相当于推动空气流动的力)
222
对于一个特定的功角,翼型上总对应地有一点(上图中C点),
空气动力Q对这个点的力矩为零,将该点称为压力中心点。
空气动
力在翼型上的作用力可由单独作用于该点的升力和阻力来表示和分析计算
升力系数、阻力系数与翼型的形状和攻角有关,其关系如下曲线所示。
可见有一段直线区,为实际运行的攻角范围,在这段区内升力系数随攻角直线增加,但在较大攻角时略有向下弯曲。
当攻角增大到
cr时,Cl达到最大值Clmax,其后则突然下降,这一现象叫做“失速”
力发电机的输出功率显著减小发电功率不超过额定值。
阻力系数在攻角不大时有最小值Cdmin,后随攻角的增加而增加。
临界攻角以后,增加更为显著。
极曲线:
以Cd为横坐标,Cl为纵
坐标,将对应于每个攻角,有确定的Cd、Cl,将点(Cd,Cl)绘制
在坐标图上,并在旁边标注攻角。
将这些连接组成的曲线叫极曲线。
极曲线放映了在攻角变化时,升力系数与阻力系数的关系。
如下图:
从极曲线原点到曲线上任何一点的矢径表示了对应攻角下总气动力系数的大小和方向。
该矢径的斜率为对应攻角下升力和阻力之比,称为升阻比,又称气动力效率。
过原点做极曲线的切线,原点与切点之间的连线的斜率为所有矢径中斜率最大者,即为最大升阻比。
切点对应风力发电机的最佳运行状态。
第三节风轮的空气动力特性
1、风轮的几何定义和参数
1)风轮直径:
叶尖旋转圆的直径,用D表示。
2)风轮扫掠面:
风轮旋转时,叶片的回转面积。
3)风轮偏角:
风轮轴线与气流方向的夹角在水平面的的投影。
4)风轮锥角:
风轮叶片与垂直与轴线的平面的夹角。
5)风轮仰角:
风轮轴线与水平面的夹角。
6)桨距角:
在叶片径向指定位置(通常叶尖)处叶片弦线与风轮旋转面间的夹角。
(故名思义,就是桨叶距离上的夹角)
7)风轮额定转速:
输出额定功率时风轮的转速。
8)风轮最高转速:
风力机处于正常状态下(空载或负载),风轮允许的最大转速。
9)风轮实度:
风轮叶片投影面积的总和与风轮扫掠面积的比值。
10)叶尖速比:
叶尖切向速度与风轮前的风速之比,用表示。
2、作用在风轮上的空气动力在叶片的径向处,取长度为的微元,称为叶素,如下图。
在风轮旋转时,叶素扫掠出一个圆环。
令风轮旋转角速度为,上游无穷远处的风速为v。
当风轮旋转时,如果在风轮旋转切线方向空气静止,则叶素与空气在旋转切向上有速度为r的相对运动,实践上由于风轮对空气的反作用力,使得空气向风轮旋转方向相反的方向旋转,这样叶素在切向与空气有相对
运动速度(1a)r,为切向气流诱导因子。
上游无穷远处的风到达风轮前,由于受到风轮阻挡而减速,实际到达风轮的速度为v(1a),为轴向气流诱导因子。
叶素上风速与受力如下图:
叶素上的合成风速
合成风速与旋转平面之间的夹角是,则
4-8)
4-9)
攻角为
4-7)
为叶素处桨距角
叶素上的升力,垂直与合成风速
12
Lw2cClr
2l
叶素上的阻力,平行与合成风速
12
Dw2cCdr
2d
作用在轴向上的力和切向上的力分别为
4-10)
LcosDsin
LsinDcos(4-11)
第四节贝兹极限讨论水平轴风力机能从风中汲取的最大功率问题。
基本风轮是理想的:
1)风轮简化为一个平面桨盘,叶片数无穷多。
平面桨盘称为致动盘。
2)风轮叶片旋转时不受摩擦阻力,不考虑能量损耗。
3)风轮前、风轮扫掠面、风轮后气流都是均匀的定常流。
风机采能简化为单元流管模型。
4)风轮前未受扰动的气流静压力和风轮后远方的气流静压力相
5)作用在风轮上的推力是均匀的。
6)不考虑风轮后的尾流旋转。
7)不考虑空气的压缩。
风轮才能物理过程用下图表示:
v——风轮前无穷远处的风速;(m/s)p——风轮前无穷远处的气流静压;(Pa)vd——风轮处的风速;(m/s)pd——风轮处,风轮前的气流静压;(Pa)pd——风轮处,风轮后的气流静压;(Pa)vw——风轮后远端风速;(m/s)p——风轮后远端气流静压。
(Pa)根据流量连续性有:
AvAdvdAwvw(4-12)
——空气密度(kg/m3);
——横截面面积(m2)。
风轮使得气流速度发生变化,风轮处的风速vd(1a)v(4-13)
——为轴向诱导因子
4-14)
应用动量定律
Fm(vvw)Advdt(vvw)
Ftt
(vvw)Advd(4-15)
(vvw)Adv(1a)
——气流所受到的作用力(N)。
气流所受到的作用力是气流对风轮作用力的反作用力。
气流对风轮作用力由风轮前后气压不同而产生。
即
F(pdpd)Ad(4-16)
从而有
(pdpd)Ad(vvw)Adv(1a)(4-17)根据伯努利方程(对于重力场中的不可压缩均质流体,满足方
1
程pgh1v2c式中、、分别为流体的压强、密度和速度;为
2
铅垂高度;为重力加速度;为常量。
)风轮前,有:
11
pghv2pddghddvd2(4-18)
22
由于气体不可压缩,d水平高度一致,hhd,(4-18)变为
1212
pv2pdvd2(4-19)
22
同样,风轮后有:
1212
pwvw2pddvd2(4-20)
22
上两式相减得:
pdpd1(v2vw2)(4-21)
2
将(4-21)代入(4-17)得
12(v2vw2)Ad(vvw)Adv(1a)(4-22)
整理得
vw(12a)v(4-23)
(4-23)表达了在风轮后静压力等于风轮前无穷远处静压力的地方风速与风轮前无穷远处风速的关系。
由(4-13)和(4-23)还可推得:
vdvvw
d2
风轮吸收的功率等于到达风轮的气流单位时间内动能的变化11
Pm(v2vw2)Advd(v2vw2)(4-24)
2w2ddw
将(4-13)和(4-23)代入:
P1Advd(v2vw2)
2
1Ad(1a)v(v2(12a)2v2)(4-25)
2d
2Adv3a(1a)2
可见风轮汲取的功率是轴向诱导因子的函数,为求取最大功率,对(4-25)式求对的导数:
dP2Adv3(1a)(13a)(4-26)
da
11
令(4-26)等于0,得解1、1,当a1时,风轮汲取功率为最
33
大
定义风能利用系数
其分母为截面积为Ad的自由流所具有的功率
将(4-25)代入(4-27)得