风力机叶轮设计与叶片空气动力学仿真分析

1、风力机叶轮设计与叶片空气动力学仿真分析 一、概述 作用在叶轮上的空气动力是风力机最主要的动力来源，也是造成各个零部件的主要的载荷来源。要计算风力发电机组的载荷就必须先计算出空气作用在叶片上的作用力。除了气动载荷，风力机叶片在风机工作中受到的作用力主要还受到惯性力，特别应该考虑重力、离心力和陀螺力等。 风力机的叶片设计分气动设计和结构设计两大部分，气动性能计算为气动设计结果提供评价和反馈，并为叶片的结构设计提供气动载荷等原始数据。气动性能计算的准确性，直接影响叶片的气动性能和结构安全，从而影响风力机的运行效率和运行安全。二、风力机叶片几何参数 1风力机叶片翼型几何参数和气流角翼型是组成风力机叶片

2、的基本元素，因此，翼型的气动特性对风力机的性能起着决定性的作用。 以一个静止的翼型为例，其受到气流作用，风速为v，方向与翼型截面平行。 图1 翼型的几何参数和气流角 翼型的尖尾(点B)称为后缘。 圆头上的A点称为前缘，距离后缘最远。 l翼型的弦长，是两端点A、B连线方向上翼型的最大长度； C翼型最大厚度，即弦长法线方向之翼型最大厚度； 翼型相对厚度，通常为10%15%； 翼型中线从前缘点开始，与上、下表面相切诸圆之圆心的连线，一般为曲线； f翼型中线最大弯度； 翼型相对弯度，； 攻角，是来流速度方向与弦线间的夹角； o零升力角，它是弦线与零升力线间的夹角； 升力角，来流速度方向与零升力线间的夹

3、角。 =+0 （21） 此处0是负值，和是正值。 2.NACA四位数字翼型族NACA四位数字翼型分为对称翼型和有弯度翼型两种。对称翼型即为基本厚度翼型，有弯度翼型由中弧线与基本厚度翼型迭加而成。中弧线为两段抛物线，在中弧线最高点二者水平相切。四位数字翼型的表达形式为NACAXXXX第一个数字表示最大相对弯度的百倍数值;第二个数字表示最大弯度相对位置的十倍数值;最后两个数字表示最大相对厚度t的百倍数值。例如NACA2410翼型，其最大相对弯度为9%;最大弯度相对位置为30%;最大相对厚度为18%。 图2 攻角=10.5时叶片受力 图3 叶片升力系数和阻力系数三、叶轮CFD仿真 图4 叶片在气流中

4、的受力分析已知条件：风速v=8m/s ， 叶轮外径D=250 mm ， 叶轮厚度b=40 mm， 轮毂直径d=100 mm 若： 叶轮转速n=200rpm， 线速度u= Dn/60=2.6 m/s； 则气相角I=arctg(v/u)=arctg(8/2.6)=72； 选翼型为NACA2410时，攻角=10，可得：安装角=I-=72-10=62；又若：线速度u= 8m/s，则：转速n=60u/D= 611 rpm I=arctg(v/u)=arctg(8/8)=45；=I-=45-10=35 表1：fluent仿真实验数据Blade5h6h7h99 h12h18h30354062457440s7

5、474454574Pressure Differences (Pa)0.270.2550.3260.1280.2190.1300.3820.1750.0970.2830.32270.212Forces (N) X-C1.61.521.8770.661.2380.5142.230.9990.5931.5240.6470.647Moment（NM）-0.062-0.07-0.07-0.047-0.075-0.036-0.09-0.025-0.058-0.105-0.122-0.071 注：h表示d=100mm，否则d=150mm；四、理论计算 叶片升力 Fy=Cl×v2S/2 =1.1&

6、#215;1.225×82×0.003 =0.12kg五、分析结论： 1、叶片越多，扭矩越大； 2、安装角为30°45°时，扭矩较大，且无明显差异； 3、由于叶片尺寸较小，所以叶轮的气动特性差异不明显，翼型影响也较小。附录、各叶轮仿真实验详细数据1、5叶片， （1）安装角30Pressure Differences (pascal)From BoundaryTo BoundaryValueinletoutlet0.270037912556Pressure Differences (Total) (pascal)From BoundaryTo Bounda

7、ryValueinletoutlet0.266688413503Velocities (m/s)Boundary or RegionAverageMinimumMaximumoutlet8.000006616387.542238095278.03176144015inlet8.08.08.0hub0.00.00.0blade0.00.00.0flow-region7.998870189290.013.0147641438Pressures (pascal)Boundary or RegionAverageMinimumMaximumoutlet-0.509536461268-0.8370223

8、64140.0inlet-0.239498548712-0.3618087172510.113551989198hub-6.11743633104-115.75695037841.5277099609blade-10.1038917169-367.303131104129.74105835flow-region-0.474842265349-367.303131104129.74105835Pressures (Total) (pascal)Boundary or RegionAverageMinimumMaximumoutlet38.693813800734.298165131338.982

9、9790228inlet38.960502214238.838192045739.3135527521hub4.48389157665-90.802159013746.5633526369blade3.19278411713-360.95963736155.937268562flow-region38.7416665659-360.95963736155.937268562 （2）安装角40Pressure Differences (pascal)From BoundaryTo BoundaryValueinletoutlet0.326839743087Pressure Differences

10、 (Total) (pascal)From BoundaryTo BoundaryValueinletoutlet0.321891636562Velocities (m/s)Boundary or RegionAverageMinimumMaximumoutlet8.000010273617.519436261448.03983576077inlet8.08.08.0hub0.00.00.0blade10.00.00.0blade0.00.00.0flow-region7.997566485460.014.35304269Pressures (pascal)Boundary or Region

11、AverageMinimumMaximumoutlet-0.441587868739-0.6491182446480.0217580795288inlet-0.114748125652-0.2611950337890.32120808959hub-5.23889354767-129.04757690441.2422904968blade1-10.5344783064-285.36633300841.6371612549blade-9.45141871405-331.97903442441.6660194397flow-region-0.405796904955-331.97903442441.

12、6660194397Pressures (Total) (pascal)Boundary or RegionAverageMinimumMaximumoutlet38.763361000734.321014632239.1302007179inlet39.085252637338.938805729239.5212088525hub5.66634001659-112.10871334346.2597114622blade12.46701758788-242.96583283550.9565993154blade3274.5836556154.2874234729flow-region38.8087641782-274.5836556170.3868095757Mass Flow Rates (kg/s)BoundaryValueoutlet-61.456882inlet61.456909Net2.70000000029e-005F