空气动力学课程设计--2KW风力发电机组低风速叶片设计和分析.doc

课程：空气动力学2KW风力发电机组低风速叶片设计与分析 姓 名：余辉 学 号：201580812008 指导教师：傅彩明2016年01月05日2KW风力发电机组低风速叶片设计与分析一. 课题研究的背景和意义1941年，美国把蒙特研制的第一台风力发电机开启了风力发电;此后，大型风力涡轮发电机促进了风力发电;如今,世界许多国家都安装了超大型风力发电机进行风力发电,促进风力发电长远发展。近年来全球的风力发电发展很快，装机容量的年平均增长率超过了30%。风能是一种技术比较成熟、很有开发利用前景的可再生能源之一。开发利用风能对世界各国科技工作者具有极强的魅力，从而唤起了世界众多学者致力于风能利用方面的研究。本文将对风力发电基本原理和具体2KW风机叶片设计进行论述。目前，全球都面临着能源枯竭、环境恶化、气温升高等问题，日益增长的能源需求、能源安全问题受到世界各国广泛关注。风能具有可再生、资源广、安全、清洁、无燃料风险等优势，因此，世界各国都在加快风力发电技术的研究，以缓解越来越重的能源和环境压力。中国是世界上最大的煤炭生产国和消费国，提供电力的能源消费是以煤炭为主，燃煤发电量占总发电量的80%。但是，能为人类所用的石化资源是有限的，据第二届环太平洋煤炭会议资料介绍，若不趁早调整以石化能源为主体的能源结构，终将导致有限的石化能源趋于枯竭，人类生态环境质量下降的恶性循环，不利于经济、能源、环境的协调发展。二. 风力发电机的设计理论风力发电是通过捕风装置的风轮将风能装换成机械能，再将机械能转换成电能的过程，因此构成风轮的翼型的结构性能直接影响着分风能的转换效率。本章介绍风力机翼型的几何结构、空气动力学基础概念及基础理论，为下文的叶片分析设计奠基础。2.1风力机的基本概念（1）风力机的基本概念和参数风轮叶片的几何形状不同，则空气动力特性也不同。为了设计风机，必须对风机的有关的概念和术语加以理解，例如，风轮、叶片、叶片旋转平面、风轮直径、叶尖速比等，而翼型外形由翼的前后缘、弦、中弧线、翼的上下表面、叶片安装角、攻角、来流角、最大厚度及最大相对厚度、弯度与弯度分布等参数决定。（2）叶片空气动力学相关概念风力发电机的叶片是细而长的结构，相对于流动方向的速度分量，其叶展方向的速度分量通常很小，因此，通常假定在给定径向位置处的流动是二维的，这样就可以使用二维翼型数据对叶片的气动性能进行分析。假定翼型处于静止状态，令空气以相同的速度吹向叶片时，作用在翼型上的空气动力将不改变其大小。空气动力只取决于相对速度和攻角的大小。由于受翼型表面形状的影响，作用在翼型表面上的空气压力是不均匀的。翼型的上表面压力低于周围的气压，称为吸力面；下表面压力则高于周围气压，称为压力面。由伯努利理论，翼型的上表面气流速度较高，下表面的气流速度则比来流低。因此翼型的周围可以看做是两类气流的合成。第一类是当翼型置于均匀的气流中，在零升力条件下流过翼型的气流；另一类是围绕翼型的环流，从下表面流回到上表面。翼型的升力就是由后者产生的。作用于翼型截面上的空气动力可以由升力、阻力和俯仰力矩来表示。对于迎角的各个值总有一特殊点，使得空气动力对这一点的力矩为零，称为压力中心（也称气动中心）。空气动力在翼型截面上的影响可以由单独作用于压力中心点的升力和阻力来表示。当气流流过翼型时，叶片下表面的压力大于周围空气的压力，而上表面的压力则要小与周围的空气压力。因此，在叶尖空气要从下表面流向上表面，结果在叶稍产生涡旋。在叶片中部的对称面两边的旋涡具有不同的旋转方向，并且在离开叶片后面不远的地方翻卷长两个孤立的大旋涡。旋涡的不断形成以及叶片运动参数的变化，他们所需的能量供给必然减少气流对叶片所做的功，所以这些旋涡引起的后果就是使得阻力增加，由此产生的这部分阻力被称为诱导阻力。2.2 风力机的基本气动理论风力发电是一门包含流体力学、空气动力学、材料力学、机械加工等各方面知识的复杂学科，随着对风力发电研究的深入，特别是流体力学和空气动力学的发展及其在风力发电机设计上的应用，形成了许多关于风力发电机设计的理论，其中最重要的有贝兹理论、叶素理论和葛劳渥旋涡理论。三.叶片设计过程3.1 设计基础较好的风力机必须具有良好的气动性能，以获得较高的风能利用系数和较大的经济效益。风力机的气动性能主要表现为叶片的气动性能。因为叶片的设计对风力机性能有重大影响。叶片的核心设计包括：计算风轮直径D，确定叶片数B，选取叶素翼型，计算各叶素弦长C和安装角。叶片分析设计流程为：（1）确定风机特征风速及电机功率P；（2）计算风轮直径D；（3）确定风机尖速比；（4）确定风机叶片数B；（5）选取叶片翼型；（6）确定各叶素弦长C和安装角；（7）计算叶片性能参数；（8）计算叶片动力学特征。3.2 风力机特征风速的确定风速变化大，很难用一种数学模型准确地描述。对风速描述的方法：以实测典型风的各种风速的频率曲线为基础，模拟风速的频率曲线，用概率统计的方法建立估算风速的数学表达式。世界常用的估算风速的函数有两种：（1）瑞利（Ray leigh）函数分布；（2）威布尔（Weibull）函数分布。3.3 叶片基本设计方法3.3.1 风轮直径计算的五个模型（1）模型一 对给定的风力机，其风轮直径D为： （3-1）式中，风力机输出功率；风能利用系数；机械效率； 空气密度，1.225；额定风速；（2）模型二 对模型一进行简化得： （3-2）（3）模型三 根据以往实践，对于现代高速风力机风轮直径D可粗略的估算为： （3-3）（4）模型四 对于大型风力发电机，其风轮直径D可粗略的估算为： （3-4）模型二、三、四是在一的基础上做的某些简化，相比，模型一计算精度较高。（5）模型五 考虑温度、高度对空气密度的影响，水平轴风力发电机的风轮直径： （3-5）式中，单位换算系数； 空气高度密度换算系数，不同海拔高度空气密度的修正系数；空气温度密度换算系数，不同温度时空气密度的修正系数；风力机总效率，风力机的总效率一般取；低速风力机取小值，高速风力机取大值；一般设计时高速风力机取3050。3.3.2 叶片相关设计参数确定（1）尖速比风轮的尖速比是风轮的叶尖线速度和设计风速之比。尖速比与风轮效率密切相关，在风力机没有超速的条件下，运转于高尖速比状态下的风力机具有较高的风轮效率。尖速比-风能利用系数的影响情况如图3.1所示。从性能曲线可知，不管叶尖速比高或低，风能利用系数都不是最优，只有在某个中间状态，才可达到最佳。若风力发电机在整个运行区域内，都可保持在这个最佳叶尖速比状态，则风能利用效率就是最好的。通常，高速风力机尖速比在6之间时，风力机具有较高的风能利用系数。图3.1 尖速比对风能利用系数的影响（） 叶片数B风轮的叶片数取决于叶片的尖速比，风力发电的高速风力机一般取5 ，性能更为优越的三叶片风力机的应用较为广泛。贝茨理论和涡流理论基于无限叶片数，有限风力机叶片数B对风力机效率存在影响机理做了阐述，其方法应用于风力机时，在正常负载情况下，其风能利用系数与风洞试验结果接近。（3）翼型翼型的选取对风力机的效率十分重要，叶片通常由翼型系列组成。较好的翼型应该是在某一攻角范围内升力系数 较高，而阻力系数较小；它所适应的雷诺数与风力机实际运行情况的雷诺数相近；且具有较高的结构强度和良好的制造工艺性。由于叶片根部各翼型力臂较小，对风力机风轮输出扭矩贡献不大，所以叶根对风力机性能影响较小，主要考虑加工方便和强度问题。在尖部采用薄翼型以满足高升阻比的要求；根部采用相同翼型或较大升力系数翼型的较厚形式，以满足结构强度的需要。翼型数据选取步骤如下：(1)选取雷诺数R，选取与风力机实际运行时的雷诺数相近的值；(2)选取最佳攻角及升阻系数，选取相近雷诺数附近的最大升阻比所对应的角度作为攻角，再由攻角确定升力系数和阻力系数等特征值。3.3.3 叶展的葛劳渥（Glauert）设计模型葛劳渥（Glauert）设计模型是考虑了风轮后涡流流动的叶素理论。其设计模型有两种，一种未引入干扰系数，另一种则引入了干扰系数。（1）模型一中间参数计算： （3-6）弦长C： （3-7）安装角： （3-8）（2）模型二风轮半径r处的叶素对风轮轴功率的贡献量为： （3-9）风能利用系数： （3-10）求最大风能利用系数，即求式（3-13）的条件极值，通过运算可得到上式的极值条件为： （3-11）这样对应一个值就可以利用式（3-11）求得相应的轴向干扰系数a及切向干扰系数b的值。通过以上各式可得弦长C： （3-12）安装角： （3-13）3.4 具体2KW型风机叶片设计3.4.1 确定风轮直径D本文选择风轮直径计算模型一，即式（3-1）计算风轮直径D：式中，风能利用系数； 发电机的机械效率；空气密度； 设计风速。取风轮直径。3.4.2 确定尖速比由于通常，高速风力机尖速比在6时，风力机具有较高的风能利用系数，因此本文选取尖速比3.4.3 确定叶片数B三叶片风机的运行和输出功率较为平稳，目前风机多采用三叶片，因此本文也将采用三叶片。3.4.4 确定翼型本文的翼型选取与翼型基本气动性能计算借助于Profili软件，Profili软件是专业进行翼型设计和分析的空气动力学分析设计软件。该软件翼型库量大，且可根据需要设计新翼型；其气动性分析以专业气动性分析软件XFOIL为基础，可针对不同需要从不同角度对现有翼型进行气动性分析，计算全面且精度高。该软件的应用不仅简化了叶片设计的翼型选取环节，而且提高了翼型外形数据和气动性能数据的精度，进而提高了叶片设计精度。由于各种翼型具有不同的安装角、升阻比、尖速比和叶片扭曲，因此，各种翼型的捕获风能的能力，抗弯强度，降噪能力等等各不尽同，而复合叶片可以综合几种翼型的优点，获得一种综合性能好的叶片。所以本文采用风力发电机专用新翼型S822和S823翼型，将两者进行组合使用。此类翼型具有更大升力系数以及更小的外形阻力，还具有对表面粗糙度不敏感的特性。S822翼型形状如图3.2所示，其最大弯度16.01在39.2翼弦处，最大曲面1.89在59.6翼弦处，翼型前缘半径0.6452，翼型后缘厚度为0。图3.2 S822翼型形状 图3.3 S823翼型形状S823翼型形状如图3.3所示，其最大弯度21.14在24.3翼弦处，最大曲面2.51在70.5翼弦处，翼型前缘半径1.0749，翼型后缘厚度为0。（1）雷诺数的选取由于叶片运行于10米左右的空间环境内，所以由Profili软件计算10米高空空气的雷诺数得=500000。（2）设计攻角及升阻系数1、S822翼型的攻角及升阻系数应用Profili软件对S822翼型数据进行分析得到图3.4和图3.5以及表3.1所示的结果：图3.4 S822升阻系数随攻角的变化情况图3.5 S822升阻比及力矩系数随攻角的变化情况表3.1 S822不同攻角下的各系数值表S822= 500000AlfaClCdCl/CdCm-8-0.51870.0128-40.5234-0.0863-7.5-0.48180.0123-39.1707-0.0835-7-0.44470.0117-38.0085-0.0807-6.5-0.40420.0112-36.0893-0.0786-6-0.35860.0107-33.514-0.0772-5-0.26590.0099-26.8586-0.0747-4.5-0.21870.0095-23.0211-0.0735-4-0.16950.0092-18.4239-0.0727-3.5-0.11720.0089-13.1685-0.0723-3-0.06360.0087-7.3103-0.0721-20.04290.00864.9884-0.0716-1.50.09740.008611.3256-0.0715-10.15420.008518.1412-0.0718-0.50.21170.008524.9059-0.072200.2670.008730.6897-0.07230.50.32040.008736.8276-0.07191.50.43530.008551.2118-0.072420.48760.008557.3647-0.07172.50.54110.008464.4167-0.071230.59790.008570.3412-0.071440.69680.008581.9765-0.068650.78420.008691.186-0.06355.50.82130.008992.2809-0.059760.84040.009192.3516-0.05216.50.83550.009885.2551-0.040170.82170.010975.3853-0.02727.50.82690.012168.3388-0.0198.50.88960.015457.7662-0.015790.90940.01753.4941-0.01159.50.93330.018550.4486-0.008100.95680.0247.84-0.004510.50.97920.021745.1244-0.0011111.00590.023343.17170.001511.51.0180.025839.45740.0056121.05010.027138.74910.007312.51.06760.029436.31290.0102131.0760.032633.00610.0137从如图3.5所示的曲线中可以大致知道，在攻角为6附近升阻比为最大。由表3.1不同攻角下的升、阻力系数可知，S822翼型在攻角为6时具有最大的升阻比/=92.3516，此时升力系数=0.8404，阻力系数=0.0091。2、S823翼型的攻角及升阻系数应用Profili软件对S823翼型数据进行分析得到图3.6和图3.7以及表3.2的结果。从如图3.7所示的曲线中可以大致知道，在攻角为7附近升阻比为最大。参见不同攻角下的升、阻力系数表3.2可知，S823翼型在攻角为7时具有最大的升阻比/=85.5366，此时升力系数=1.0521，阻力系数=0.0123。图3.6 S823升阻系数随攻角的变化情况图3.7 S823升阻比及力矩系数随攻角的变化情况表3.2 S823不同攻角下的各系数值表S823= 500000AlfaClCdCl/CdCm-6-0.33250.014-23.75-0.1276-5.5-0.28420.0134-21.209-0.1274-5-0.23110.0132-17.5076-0.1276-4.5-0.17650.0133-13.2707-0.1278-4-0.12360.013-9.5077-0.1274-3.5-0.06930.0125-5.544-0.1277-3-0.0130.0121-1.0744-0.1279-2.50.04150.01173.547-0.1276-1.50.15160.010714.1682-0.127-10.20740.010719.3832-0.1267-0.50.26440.010724.7103-0.126900.32130.010530.6-0.12690.50.37770.010735.2991-0.126710.4350.010541.4286-0.12671.50.49160.010546.819-0.12662.50.6050.010657.0755-0.126130.65970.010861.0833-0.12553.50.71490.010866.1944-0.124940.76790.010970.4495-0.12384.50.82120.011273.3214-0.122850.87330.011277.9732-0.12165.50.92270.011580.2348-0.119860.970.011782.906-0.11776.51.01570.011985.3529-0.115271.05210.012385.5366-0.11117.51.08210.012983.8837-0.105981.1070.013780.8029-0.18.51.12270.014975.349-0.093291.13570.016469.25-0.08639.51.14640.018163.337-0.0798101.1590.019958.2412-0.074110.51.17370.021953.5936-0.0692111.18940.02449.5583-0.065111.51.2010.026645.1504-0.0611121.21410.029241.5788-0.057512.51.22950.031838.6635-0.05443.4.5 叶展形状设计计算考虑到叶片需要与轮毂配合，同时要考虑翼型与连接部分的过度以及强度要求，叶片小于0.1倍直径R部分采用过度设计而不采用翼型。叶片设计主要考虑0.2RR采用翼型的部分，主要计算各个叶素的弦长C和安装角。本文选择葛劳渥（Glauert）设计模型二来计算各个叶素的弦长C以及安装角。为满足叶尖部分高升阻比，根部具有较高结构强度的要求，将叶片0.2RR段沿叶展方向分成40段每段57mm，共计41个截面。其中127截面采用较厚的S823翼型，2841截面采用升阻比较大的S822翼型。根据葛劳渥（Glauert）设计模型二，并通过运用工程计算软件Matlab进行编程计算得到各个叶素的弦长和安装角，并通过函数拟合来对其所得结果进行修正。修正前后所得到的弦长和安装角如下表所示：表3.3 叶片设计参数列表叶素编号半径（mm）轴向干扰系数径向干扰系数弦长C（mm）安装角（）所选翼型修正前修正后修正前修正后10.25700.32070.1335479473.835819.53719.7539S82320.226270.32240.1127460452.893417.764518.0325S82340.267410.3250.0832421414.337914.773915.0037S82350.287980.3260.0726403396.658213.508513.679S82360.38550.32680.0638386379.999312.369712.4695S82390.3610260.32860.0452339335.81329.56169.4554S823100.3810830.3290.0408326322.90348.78818.631S823120.4211970.32980.0337301299.6127.42967.2112S823150.4813680.33050.026269270.4945.76545.5621S823170.5214820.33090.0223252254.57184.84754.7175S823180.5415390.33110.0207243247.55764.43494.3567S823190.5615960.33120.0193236241.12994.04934.0308S823210.617100.33150.0168222229.90113.34943.4676S823230.6418240.33170.0148209220.61812.73112.9982S823240.6618810.33180.0139203216.62292.44822.7907S823250.6819380.33190.0131198213.01392.18112.5977S823260.719950.3320.0124193209.75781.92832.4168S823300.7822230.33220.01218199.59472.04091.7765S822310.822800.33230.0095213197.60221.84551.6304S822340.8624510.33240.0083199192.40631.31191.2143S822350.8825080.33250.0079194190.82361.14961.0828S822360.925650.33250.0076190189.25980.99430.9556S822380.9426790.33260.0069183186.05550.70290.7176S822400.9827930.33260.0064175182.52630.43460.5102S82241128500.33270.0061172180.55640.30820.4225S822其中在对叶片的弦长和安装角进行拟合修正时，以为自变量，分别以（，C）和（，）作为样点进行函数模拟。通过运用Matlab中的polyfit函数分别对弦长C和安装角采用三次、四次多项式进行修正。最后得到的拟合函数如下：为将修正前后的结果更加的形象化，将其修正前后的弦长C和安装角随半径r的变化分布情况做成线图，如图3.8和图3.9所示。 图3.8 弦长C随半径r的变化曲线 图3.9 安装角随半径r的变化情况3.5 叶片的几何建模3.5.1 叶片各叶素三维坐标处理图形变换的实质是对组成图形的各顶点进行坐标变换。本文将采用基于点的坐标几何变换理论求解叶片各截面在空间实际位置的三维坐标。其思路为：获得翼型上下弦数据以翼型气动中心为原点的二维坐标叶素各离散点空间实际坐标(x,y,z)。建立三维坐标系：设叶片根部r=0处的叶素平面为XOY平面；叶展方向为Z轴正向；原点设在r=0的叶素平面的气动中心。叶片轮廓线各离散点坐标(x,y,z)。具体求解步骤如下:(1)获取翼型数据。通过Profili软件导出S822翼型和S823翼型的数据。(2)求解以翼型气动中心为原点、翼型前后缘连线为x轴的二维坐标。设气动中心坐标为。则 （3-19）(3)求解各叶素各离散点空间实际坐标(x,y,z)。结合弦长计算各叶素坐标： （3-20）旋转叶素得实际空间坐标： （3-21）将式（3-21）代入式（3-22）有：本文将翼型气动中心确定在翼弦上距离前缘30处，并通过Matlab编程计算得到各个叶素的三维坐标数据。 3.5.2叶片几何建模将上一步计算得到的各个叶素的三维坐标参数(x,y,z)导入EXCEL中，然后直接通过solidedge中的“按表创建曲线”命令将EXCEL中的各个叶素的三维坐标分别导入以获得各叶素的轮廓线。导入solidedge中的各个叶素轮廓线如图3.10所示。图3.10 各个叶素轮廓线图3.11 叶片的几何模型通过对所获得的各个叶素轮廓线并结合叶柄数据进行放样拉伸以及除料等处理，进一步生成叶片的整体三维模型，如图3.11所示。四. 叶片静态气动性能计算4.1 风力机气动性能指标在描述风力机的气动性能时，通常包括风轮