风力机组气动特性分析与载荷计算-1

1目录1前言22风轮气动载荷...............................................................22.1动量理论.................................................................................................................................22.1.1不考虑风轮后尾流旋转.................................................................................................22.1.2考虑风轮后尾流旋转.....................................................................................................32.2叶素理论.................................................................................................................................42.3动量──叶素理论.................................................................................................................42.4叶片梢部损失和根部损失修正.............................................................................................62.5塔影效果.................................................................................................................................62.6偏斜气流修正.........................................................................................................................62.7风剪切.....................................................................................................................................63风轮气动载荷分析...........................................................73.1周期性气动负载......................................................................................错误!未定义书签。4.1载荷情况DLC1.3..................................................................................................................104.2载荷情况DLC1.5..................................................................................................................104.3载荷情况DLC1.6..................................................................................................................104.4载荷情况DLC1.7..................................................................................................................114.5载荷情况DLC1.8..................................................................................................................114.6载荷情况DLC6.1..................................................................................................................112风力发电机组气动特性分析与载荷计算1前言风力发电机是靠风轮吸取风能的，将气流动能转为机械能，再转化为电能输送电网，风力机气动力学计算是风力机设计中的一项重要工作。特别是对于大、中型风机，其意义更为重大。风力机处于自然大气环境中，大气紊流、风剪切、风向的变化（侧偏风）和塔影效应等，这些现象使叶片受到非常复杂气动载荷的作用，对风力机的气动性能和结构疲劳寿命产生很大的影响。对一台大型风力发电机组来说，除风轮叶片产生机组的气动载荷外，机舱和支撑风轮和机舱的塔筒也产生气动载荷，这些都对机组的载荷产生影响。2风轮气动载荷目前计算风力发电机的气动载荷有动量—叶素理论、CFD等方法。动量—叶素理论是将风轮叶片沿展向分成许多微段，称这些微段为叶素，在每个叶素上的流动相互之间没有干扰，叶素可以认为是二元翼型，在这些微段上运用动量理论求出作用在每个叶素上的力和力矩，然后沿叶片展向积分，进而求得作用在整个风轮上的力和力矩，算得旋翼的拉力和功率。动量—叶素理论形式比较简单，计算量小，便于工程应用，估算机组初始设计时整机的气动性能，被广泛用于风力机的设计和性能计算，而且还用来确定风力机的动态载荷，不断地被进一步改进和完善。CFD数值计算不需要对数学模型作近似处理，直接对流体运动进行数值模拟，从物理意义上说，数值求解N-S方程的CFD方法应该是最全面准确计算风力机气动特性的方法。但是，由于极大的计算工作量，数值计算的稳定性等原因，目前CFD求解N-S方程方法还远不能作为风力机气动设计和研究的日常工具。作为解决工程问题的工具还不太实际。为此在计算中应用动量—叶素理论方法来计算机组的气动载荷。2.1动量理论动量理论是经典的风力机空气动力学理论。风轮的作用是将风的动能转换成机械能，但是它究竟能够吸收多大的风的动能就是动量理论回答的问题。下面分不考虑风轮后尾流旋转和考虑风轮后尾流旋转两种情况应用动量理论。2.1.1不考虑风轮后尾流旋转首先，假设一种简单的理想情况：（1）风轮没有偏航角、倾斜角和锥度角，可简化成一个平面桨盘；（2）风轮叶片旋转时不受到摩擦阻力；（3）风轮流动模型可简化成一个单元流管；（4）风轮前未受扰动的气流静压和风轮后的气流静压相等，即p1=p2；（5）作用在风轮上的推力是均匀的；（6）不考虑风轮后的尾流旋转。将一维动量方程用于风轮流管，可得到作用在风轮上的轴向力为21VVmT（1）式中m为流过风轮的空气流量TAVm（2）于是21VVAVTT（3）而作用在风轮上的轴向力又可写成ppAT（4）由伯努利方程可得pVpVT222121（5）3pVpVT222222（6）根据假设，p1=p2，（5）式和（6）式相减可得22221VVpp（7）由（3）式、（4）式和（7）式可得221VVVT（8）（8）式表明：通过风轮的风速是风轮前的风速和风轮后的尾流速度的平均值。设定轴向诱导因子11Vuaa，ua为风轮处的轴向诱导速度，则111aVVT（9）11221aVV（10）（9）式和（10）式代入（3）式得2)1(42111AVaaT（11）)1(4)2(1121aaAVTCT（12）轴向诱导因子a1又可写成121121UVVa（13）（13）式表示，如果风轮全部吸收风的能量，即V2=0时，a1有一个最大值1/2，但实际情况不可能这样，所以a11/2。根据能量方程，风轮吸收的能量（即风轮轴功率P）等于风轮前后气流动能之差222222212221VVAVVVmPT（14）将（9）式、（10）式代入（14）式，可得2113112aaAVP（15）当dPda/10时，P出现极值，则03412211311aaAVdadP（16）a1=1和a1=1/3是（16）式的根。又因为a11/2，故只考虑a1=1/3的情况462131212aAVdaPd（17）当a1=1/3时，dPda2120/，P取极大值，由于P的连续性，因此极大值就是最大值31max212716AVP（18）相应地，功率系数CP为最大值593.027162/31maxmaxAVPCP（19）这个值被称为贝兹极限，它表明在理想情况下，风轮最大能吸收593%.的风的动能。2.1.2考虑风轮后尾流旋转实际上，风轮尾流是旋转的，这时如果风轮处气流的角速度和风轮角速度相比是个小量的话，一维动量方程仍然可用，而且假设p1=p2。风轮作用盘假设是由许多以风轮轴线为对称轴的小圆环（内半径r，外半径r+dr）构成。这时)(21VVmddT（20）而rdrVdAVmdTT2（21）假设（11）式仍然成立，则有11212VaVV（22）将（21）式、（22）式与（9）式代入（20）式可得draaVrdT)1(41121（23）作用在整个风轮上的轴向力为RrdraaVdTT01121)1(4（24）4由动量矩方程，作用在该圆环上的转矩为dMdmurt()（25）式中urt，为风轮叶片r处的周向诱导速度，为风轮叶片r处的周向诱导角速度。设定周向诱导因子a22/，为风轮转动角速度。将uart22，（20）式及（9）式代入（25）式可得draaVrdM2113)1(4（26）因此风轮轴功率为RdrraaVdMdPP031212)1(4（27）设定风轮叶尖速比1/VR，2RA，则RdrraaRAVP03124231)1(4（28）风能利用系数为CRaardrPR81242130()（29）2.2叶素理论叶素理论的基本出发点是将风轮叶片沿展向分成许多微段，称这些微段为叶素，在每个叶素上的流动相互之间没有干扰，叶素可以认为是二元翼型，将作用在每个叶素上的力和力矩沿展向积分，求得作用在风轮上的力和力矩。从动量理论可知，当考虑风轮尾流旋转后，风轮处轴向速度)1(11aVVa，周向速度)1(2arVt，实际流经风轮处的气流速度是taVVW。对每个叶素来说，是迎角，是入流角，是扭转角])1()1([211raVaarctg（30）（31）求出后，查翼型手册得到作用在叶素上的升力系数Cy和阻力系数Cx。由于dFdYdXncossin（32）dFdYdXtsincos（33）则法向力系数Cn和切向力系数Ct分别为CCCnyxcossin（34）CCCtyxsincos（35）作用在每个叶片上的叶素的轴向力为dTcdrWCn~22（36）式中c为