水平轴风力机气动性能计算模型的修正.docxVIP

水平轴风力机气动性能计算模型的修正 0 气动性能计算 计算动力性能是风力机设计和修订的重要组成部分。设计出桨叶的气动外形之后,计算其气动性能,可以作为对设计结果的评价;反之气动性能计算结果可以作为反馈,为修正桨叶气动外形提供数据,因此气动性能计算结果的准确度直接影响风力机桨叶外形优化设计的优劣;风力机桨叶强度、刚度以及稳定性的校核,需要气动性能的计算结果作为其原始输入参数,准确的气动性能计算能够使校核结果更可信。 目前流行的风力机气动性能计算方法是采用片条理论(Strip Theory),在此基础上考虑影响风力机气动性能的各种因素并对模型进行修正。风力机在实际运行时还受很多因素影响,如风剪切、偏航、风轮的结构参数和风力机安装参数等。本文基于片条理论及其模型修正,并综合考虑影响风力机气动性能的各种因素,建立了严格的气动性能计算模型,以期获得更准确的计算结果,完善所发展的风力机气动与结构设计软件Blade Design,以便更好地为风力机的设计和校核提供支持。 1 基础研究 1.1 风轮桨叶平面推力的基本概念 动量理论(Momentum Theory)主要用来估算风力机的理想出功、效率和流速。经典的动量理论定义了一个通过桨叶平面的理想流管,如图1所示,U、V1、V2分别表示来流风速、通过风轮桨叶平面的风速和桨叶平面后部尾流内的风速。 考虑风力机轴向的动量变化,可得桨叶平面上产生的推力T为: T=m(U-V2) 式中,m=ρV1A;ρ——空气密度;A——桨叶平面的面积。 应用伯努利方程可以得到: V1=U+V22 定义轴向速度诱导因子a=1-V1/U,可推导出推力的表达式为: T=2ρAU2a(1-a) (1) 桨叶平面吸收的风能转化为风轮的角速度Ω和作用在风轮上的转矩Q。定义切向速度诱导因子a′=ω/(2Ω),其中ω——桨叶平面风的角速度,那么可以得到转矩Q为: Q=πR4ρU(1-a)a′Ω(2) 式中,R——桨叶平面的半径。 1.2 叶素流动曲线 叶素理论(Blade Element Theory)将风力机桨叶简化为由有限个叶素沿径向叠加而成,因而风轮的三维气动特性可以由叶素的气动特性沿径向积分得到。相对于动量理论,叶素理论从叶素附近的空气流动来分析叶片上的受力和能量交换,从而更多的应用于风力机的设计中。 图2所示为一个叶素上的风速分量,假设该叶素在桨叶上处于半径r的位置。通过叶素的气流轴向速度为U(1-a),在风轮旋转平面内,气流相对于桨叶的角速度为Ω(1+a′),所以在该叶素处气流的切向速度为Ωr(1+a′),轴向速度和切向速度合成了相对速度W,从而该叶素的来流角为: ?=tg-1[1-aU1+a′Ωr](3) 设叶素弦长为c,则其升力和阻力分为: FL=12ρW2CLcdrFD=12ρW2CDcdr 其中,CL、CD——分别为叶素翼型的升力系数和阻力系数。 考虑桨叶数量B,半径r处叶素的轴向推力为: dT=B(FLcos?+FDsin?) =12ρW2(CLcos?+CDsin?)Bcdr(4) 转矩为: dQ=B(FLsin?-FDcos?)r =12ρW2r(CLsin?-CDcos?)Bcdr(5) 1.3 a1-ar3a1-a及型bbbdr的求解 片条理论综合了动量理论和叶素理论。考虑式(1)和式(4),式(2)和式(5),令dT动量理论=dT叶素理论及dQ动量理论=dQ叶素理论,可得: 4ρπU2a(1-a)rdr=12ρW2(CLcos?+CDsin?)Bcdr4πr3ρU(1-a)a′Ωdr=12ρW2r(CLsin?-CDcos?)? Bcdr 由以上两式可以得到a,a′的求解公式: a=Bc(CLcos?+CDsin?)8πrsin2?+Bc(CLcos?+CDsin?)(6) a′=Bc(CLsin?-CDcos?)8πrsin?cos?-Bc(CLsin?-CDcos?)(7) 1.4 气动损失系数 片条理论并没有考虑桨叶尖端和根部气流沿桨叶径向的流动情况,事实上由于桨叶两面的气压差异,在桨叶尖端和根部气流会沿叶片产生二次流动,从而引起力矩的减小。由于叶尖处的叶素受力对整个风力机性能的影响很大,所以叶尖损失不容忽视。Prandtl对叶尖区的流动作了研究,定义了叶尖损失系数为: Ft=2cos-1(eB(R-r)2Rsin?)π 轮毂部分虽然动压较小,气动阻力大而升力小,但由于弦长较大,占桨叶的比例较大,因而在计算气动性能时也有必要考虑轮毂的损失。与叶尖损失系数类似,Prandtl定义了轮毂损失系数为: Fh=2cos-1(eB(r-rhub)2rhubsin?)π 其中,rhub——轮毂半径。 那么总的损失系数为F=FtFh。应用损失系数,动量理论的式(1)和式(2)需修正成为: T