matlab风力机叶片仿真教程,基于MATLAB的小型风力机叶片设计

文章编号: 100027709(2007) 050142203 基于MA TLAB 的小型风力机叶片设计 王　军　周丙超 (华中科技大学 能源与动力工程学院, 湖北 武汉 430074) 摘要: 结合W ilson 设计方法流程, 采用MA TLAB 语言编程, 开发了小型风力机叶片气动设计的应用程序。利用该程序设计了一台10 kW 小型水平轴风力机, 采用CFD 方法对所设计的的模型用FLU EN T 软件进行单流道分析计算。流场分析结果表明, 编制的软件提高了初期设计效率, 具有直接形成结构实体的特点, 缩短了小型风力机叶片CFD 分析前处理时间。 关键词: 叶片; 参数化; MA TLAB 中图分类号: TK83 文献标志码: A 收稿日期: 20070425, 修回日期: 200705222 作者简介: 王军(19632) , 男, 副教授、博士后, 研究方向为流体机械优化设计及实验, E2mail:w angjhust@163. com 　　目前在石油、煤炭等传统能源紧张的环境下, 世界各国均将目光投向新能源领域, 其中风能倍受关注。风力机叶片作为将风能转换为动能的“心脏” , 其设计直接决定着风力机组的效率。在对风力机叶片设计时大多采用CADgCFD 的分析方 法, 用这些软件可加速设计进程、缩短研发周期, 在产品制造前可预测其性能, 方便地对产品进行优化设计[1]。本文以提高计算精度、节约设计时间、改善小型风力发电机的输出性能为目标, 使用MAL TLAB 软件编写小型风力机叶片设计程序。同时对一实例进行了CFD 分析, 证实了该程序的可行性, 为进一步CFD 设计建模提供参考依据。 1　气动设计原理及流程 1. 1　气动设计原理 风力机叶片的设计主要包括叶片数、叶片直径、翼型、转速、安装角和弦长等参数的确定, 本文采用叶尖损失的W ilson 方法对风力机叶片主要 参数进行计算[2, 3]。 直径D 可由关系式 P = 0. 2D 2V 3 1 得[2]: D = P 0. 2V 3 1 (1) 式中, P 为风力机轴功率,W ;D 为叶片直径,m; V 1 为来流风速, mgs。 考虑到叶片气动性能与制造等方面的因素, 设计的叶片各截面采用同一种翼型, 即NACA 23015 翼型。风轮转速确定: N = (600V 1)(ΠD ) (2)式中 , N 为风轮转速, rgm in;0 为叶轮叶尖速比。 由动量—叶素理论可得[2]: dCp = 8 Κ2 0 b(1 - a)F3d(3)a (1 - aF ) = b(1 + b)Κ2 (4) F = 2 p arccos e- f (5) f = B 2 R - r R sinf (6) tan= (1 - a)V 1 (1 + b)8r = 1 Κ 1 - a 1 + b (7) 　　同时, 　　　Η=Υ-Α (8) R e = V 1 (1 - a)C nsinf (9) B CC1cosf 8p rsin2f = (1 - aF )aF (1 - a) 2 (10) 式中, Cp 为风能利用系数; a 为轴向诱导因子; b为周向诱导因子;B 为叶片数; Ф为叶素入流角;Η为 叶素安装角;Α为叶素攻角;V 1 为来流风速,mgs; 8为风轮转速, radgs; r 为叶素剖面到叶根的距 离,m; R 为叶片半径; F 为Prantl叶尖损失修正因子; C 为叶素弦长; R e 为雷诺数;Μ为运动粘性系 数,m 2gs; C1 为叶素升力系数。 根据W ilson 理论, 对于