基于matlab软件风力机叶片数字化设计.docVIP

基于matlab软件风力机叶片数字化设计 　　摘要：本文以风力机叶片为研究对象,选用常用的NACA4412翼型，基于Glauert设计模型，利用matlab工具对一20KW风力机进行了叶片气动外形的设计，得到风力机叶片的外形参数。计算过程利用matlab能够简化繁琐计算便于修改参数，能够推广到其他产品的设计，为产品设计提供一定的依据。 　　关键词：风力机；叶片设计；matlab 　　中图分类号：TK474.8+11 文献标识码：A 　　1 引言 　　当今煤、石油、天然气等传统能源消耗的日益增加，能源危机日益严重，寻找新的可再生能源成为人类社会发展所面临的重要问题。风能作为一种无污染可再生的新能源有着巨大的发展潜力，风能的开发利用受到各国的广泛重视[1]。如今风力发电产业正逐步发展成为新兴产业，而风力机是风力发电的核心，风力发电的关键问题之一是风力机叶片气动外形的设计[2]，因此，研究风力机的设计具有重要意义。 　　目前国内外在风力机的设计方面已做了许多工作[3,4]，具有代表性的设计模型有基于圆盘理论的简化设计模型、基于涡流理论的Schmitz设计模型、Glauert设计模型和Wilson设计模型。本文选择具有代表性的Glauert设计模型利用matlab语言编制风力机叶片气动外形优化设计程序，得到风力机叶片的外形参数。 　　2 风力机叶片的设计 　　要求确定一个风力发电机的主要参数，当额定风速为10m/s[5]（本文选取的设计风速是针对我国风能资源丰富的地区而言）时输出功率为20KW，气体密度ρ为1.225kg/m3，风轮功率系数Cp为0.42，传动系统效率η1为0.92，发电机效率η2为0.95。 　　2．1 风轮直径D的确定 　　对于给定输出功率的风力机，由参考文献[2]得风轮直径D为: 　　式中:P-风力机输出功率；Cp-风能利用系数；ρ-气体密度；V-额定风速；η-传动装置和发电机总效率。 　　经计算D=10.64m，取风轮直径D为11m,半径为5.5m,除去中心轮毂直径,取风轮叶片长为5m。 　　2．2 叶片数N及叶尖速比λ的确定 　　风力发电用的风力机叶片数一般为2-3 片, 叶片数与叶尖速比有一定的匹配关系[6], 见表1。由于三叶片风力机的运行和功率输出较平稳, 一般中小型风力机都选用三叶片, 因此本文采用的叶片数N=3，叶尖速比λ=6。 　　表1 尖速比与叶片数的关系 　　2．3 风轮转速n的确定 　　风轮的额定转速为： 　　（2） 　　式中：λ-尖速比；V-额定风速；D-风轮直径 　　经计算得：n=104r/min 　　2．4 叶片翼型的选择 　　NACA翼型是具有代表性的传统风力机翼型，它是美国国家宇航局（NASA）前身国家航空咨询委员会（NACA）提出的。NACA4412翼型[7,8]是常用的风力机翼型，如图1所示。它的相对厚度较小，气动性能较好，数据较全，研究较多的翼型之一。其最大相对弯度是4%；最大弯度在弦长的40%处；最大相对厚度为12%。因此本文设计的叶片从叶根到叶尖全部采用NACA4412翼型。 　　图1 NACA4412翼型 　　2．5 叶片攻角的确定 　　一般来说，应该选取相近雷诺数附近、最佳升阻比所对应的角度作为叶片的攻角a,再由攻角a确定对应的升力系数Cl与阻力系数Cd。但在实际工程设计中，人们为了简化设计步骤，通常忽略了雷诺数对翼型数据选取的影响，整个叶片选取一个统一的角度作为攻角，实践也证明这样的做法完全能够满足设计精度的要求[8]。 　　NACA4412的数据如表2[7]，由表2知，NACA4412最佳升阻比为94.4435，此时攻角为6度，因此，整个叶片选取的攻角为6度。 　　表2 NACA4412的数据 　　2．6 安装角与弦长的确定 　　风力机叶片的叶展形状是风力机叶片设计的重要组成部分，直接影响风力机的气动性能，叶展形状的设计包括各叶素上弦长和安装角的设计。目前流行的叶展形状设计模型很多，如简化设计模型、Schmitz设计模型、Glauert设计模型、Wison设计模型。其中Glauert设计模型是较具有代表性的设计模型之一。 　　Glauert设计模型考虑风轮后有涡流流动的叶素理论(即考虑诱导因子a和b)，由于叶片翼型阻力和叶梢损失的影响对叶片外形设计的影响小，因此给予忽略不计。参考文献[2,5]，Glauert设计模型有两种计算过程。其中一种和Schmitz模型类似，未明确引进诱导因子。 　　中间参数计算： 　　式中, λ-为尖速比；R为风轮半径。 　　弦长： 　　（4） 　　式中，Z-叶片数； Cl-各半径r处的升力系数。 　　安装角：（5） 　　根据叶素理论，将叶片沿径向分为10等分，每