风力机叶片设计与制作课程设计

1、风力机叶片设计与制作课程设 计风力机空气动力学课程设计（综合实验）报告（2012 - 2013年度第1学期）名 题 院 班 学称：风力机空气动力学目：风力机叶片设计与制作系：可再生能源学院级：号：学生姓名：指导教师：设计周数：成 绩：日期：2014年1月11日一、目的与要求主要目的:掌握动量叶素理论设计风力机叶片的原理和方法熟悉工程中绘图软件及办公软件的操作掌握科研报告的撰写方法主要要求:要求独立完成叶片设计参数的确定，每人提供 一份课程报告每小组提供一个手工制作的风力机叶片 设计并计作翌个风力机叶片原始数据三叶片风力机功率P=6.03KW来流风速7m/s风轮转速72rpm风力机功率系数Cp=

2、0.43传动效率为0.92发电机效率为0.95空气密度为1.225kg/m35人一组，每个小组采用一种风力机翼型，翼 型的气动数据（升力系数，阻力系数，俯仰力矩 系数）已知。设计任务1）风力机叶片设计：根据动量叶素理论对各个 不同展向截面的弦长和扭角进行计算，按比例 画出弦长、扭角随叶高的分布。2）根据以上计算结果手工制作风力机叶片，给 出简单的制作说明。三一、进度计划序 号设计（实验）内容完成时间备注1叶片参数的计算3天2叶片制作1天3报告撰写1天4课程设计答辩半天四、数据计算选用翼型s8301.叶片半径的计算:由风力发电机输出功率:得，叶片直径：p = -pv3兀d2c nn8 rP 1

3、2D _土一 =：8 X 6.03 X 10 3 _ 9.863 m片 r pV 3兀Cpn 1n 2* 1.225 x 73 X 兀 x 0.43 x 0.95 x 0.92D 9.863 R _ _ 4.932m222.叶尖速比的计算： 整个叶片的叶尖速比人 _QR _ 2兀n /60 x R _ 2兀 x 72/60 x 4.932 _ 5 31 0 _ W.1设计中取9处截面，分别是叶片半径的0.15,0.20,0.30，0.90.3.各截面处翼型弦长设计（综合实验）报告确定每个剖面的形状参数N: 可根据公珞-竺9 （4求得:八 rN16兀RC =CiB 9C B人：人 2（r）2 +

4、 4i 0； 0 R 9由气动数据表查得最大升力系数C立响，取风机 叶片数B=3，不同半径处叶片弦长的计算：宙程序直接 给出结果及线性优化后修正弦长如下：弦长随叶高分布叶高 /m0.73980.98641.47961.97282.4662.95923.45243.94564.4388弦长1.00.90.60.50.40.30.20.20.2/m895025553083134477997632345修正0.80.70.60.50.40.30.20.20.2弦长/m114342214083134477997632345图如下（系列1为计算弦长；系列2为修正后弦长）弦长随叶高的分布-一修正后的弦长

5、T-未修正的弦长fj ： 0 =。一 a，a为翼型临界攻角， 系数对成的攻角为其中为各截面处的 且由气动数据表查得最 a =5即为翼型临界攻角 根据相关关系式就可以通过迭代方法求得轴向诱 导因子a和周向诱导因子b，迭代步骤如下：1）假设a和b的初值，一般可取0；2）计算入流角； 3）计算扭角0 =。-a；4）根据翼型空气劫综嘿唯曲线得到叶素的升 力系数Cl和阻力系数Cd；5）计算叶素的法向力系数Cn和切向力系数CtC = C cos + C sin 巾6）计算a和b的新值十C，洲弋顽Bcc =2 7Ta _ cCn1 - a 4 F sin2b cCT4 1 + b 4F sin cos7）比

6、较新计算的a和b值与上一次的a和b值, 如果误差小于设定的误差值（一般可取0.001）， 则迭代终止；否则，再回到（2）继续迭代。其中，F为普朗特叶尖损失修正因子，且F - arccos兀由以上迭代公式得到计算程序（见附录），并由程序计 算得到个界面处的入流角：因为根部有轮毂，故从15%处开始计算扭角叶 高 /m0.73 980.98641.47961.97282.46 62.95 923.45243.94 564.43 88扭 角 n/28.825.119.214.511.28.766.875.263.4920132334659728579017869计算直接由程序给出扭角沿叶高的分布图如下

7、扭角随叶高的分布五、1.将打印好的翼型贴在硬纸板上，沿线剪下；在各截面1/4处确定为转轴的位置；扭角的固定：将翼型截面转至相应的扭角位画一条直线，在直线上选取两点挖点并以两根铁叶片制作选取比例1： 5；2.利用绘图工具分别绘出不同界面的翼型，打印;置，丝将翼型穿起，形成骨架；6.外形生成：将纸裁成条后，沿支架贴于表面， 形成叶片。课程设计心得体会10月8号，叶片厂实习。从学校到保定华翼风电 叶片厂有超过3个小时的路程，初下车，迎面而来的寒风让人清醒，闯入视野的庞大的风力机叶片让我感 到十分震撼。厂方的工作人员热情地接待了我们，耐心而细致的给我们讲解风力机叶片设计制造工艺。我提问讲解人员三个问题

8、：一个是叶尖处有两个小孔，它们的功能是什么，她回答，是排水用的；另一个问题是我看到叶片为什么没有防雷装置，她回答说实验的叶片，没有布置，我接着问防雷装置的布置是怎样=i的，她有给我耐心的讲解；最后一个问题是叶片与轮 毂的连接，她给我们介绍几种联接方式，预埋双头螺 柱、打孔进行螺栓联接等工艺。接着又到生产车间学习叶片的生产工艺，先进的复 合材料与严密的工艺流程给我很大的感触，想想现在 自己还在时不时的在做事情的时候粗心大意，觉得自 己还有很长的路要走。后来讲解人员又带领我们到叶片测试中心参观。高 耸的测试塔，密密麻麻的数据线让人眼花缭乱，才清 楚叶片设计制造与生产检测原来是那么复杂的事情。 这也

9、让我想起了材料力学的力学性能检测试验，我问 工作人员，是不是采用应变片做为传感器，她给了肯 定的回答，我非常高兴。剩下的就是叶片的设计与模型制作了。我们在动员 大会结束后，我、xxx、xxx、xxx、xxx立即结成小组, 我组建了一个QQ讨论组，把他们都邀请进来，然后 我和小组成员讨论了翼型选择，确定后分配了任务， 大家团结协作，高效地完成了任务，我们都非常高兴。在叶片骨架制作的过程中，合实选取材料不得当，铁丝太 细，强度不够，我们有马上想办法解决，最后用多年 积攒的一次性筷子解决了问题。虽然制作的不是很精 美，但我们真心付出了自身的脑力、体力，是我们团 结协作的结果，从这个角度来看，是完美的

10、。提供的宝贵资源与大力的支持，我在课程设计中收获 很多，学习到了课本中难以把握的工程实际知识，锻 炼了实际操作软件与动手能力，提高了自身的专业素 质，对风力机空气动力学原理与风力机叶片设计制造一周的风力机空气动力学课程设计马上就要结束 了，这一周过得紧张而又充实，非常感激学校、老师有了具体、深入的认识。风力机空气动力学课程设计 中，我还体会到一些做人做事的道理，如团队意识与 组织协调，专心、细心、耐心这些只有在活动中才能 感悟到其重要性的东西。七、参考资料与应用软件风力机空气动力学贺德馨等著，华北电力大学校内教材;AutoCAD Autodesk （欧特克）公司出品;MATLAB 2013b

11、美国 MathWorks 公司出品。八、附录1.程序(mtatlab):R=4.932;B=3;a=zeros(1,9);b=zeros(1,9 );n=zeros(1,9);sig=zeros(1,9);F=sig;g=5* pi/180;Cl=1.0789;Cd=0.0216;r=0.15 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9*R;Clm=1.528;jihe0=load(s830.dat);lam=5.31;lamd=T设计(窍哥报2 0.3 0.4 0.5 0.6 0.7 0.8 0.9*lam;c=zeros(1,9);v1=7;w=72*2*pi/60;f

12、or i=1:9a0=0.5;b0=0.5;a1=0;b1=0;c(i)=16*pi*R /9/Clm/B/lam/(sqrt(4/9+lam*lam*r(i)*r(i) /R/R);sig(i)=B*c(i)/2/pi/r(i);while(abs(a1-a0)=0.001|abs(b1-b0)= 0.001)a0=a1;b0=b1;x=atan(1-a0)*v1/(1+b0)/w/r(i);F(i)=2/pi*acos(exp(-B*(R-r(i)/2/r(i )/sin(x);n(i)=x-g;Cn=Cl*cos(x)+Cd*sin(x);Ct=Cl*sin(x)-Cd*cos(x);t

13、a=sig(i)*Cn/4/F(i)/sin(x)/sin(x);tb =sig(i)*Ct/4/F(i)/sin(x)/cos(x);a1=ta/(1+ta);b1=tb/(1-tb);enda(1,i)=a1;b(1,i)=b1;jihe=jihe0*c(i);plot(jihe(:,1),jihe(:,2);axis equal;hold on;endtitle(不同截面处的翼型大小)hold offn=n*180/pi;disp(弦长分布矩阵：)，c disp(a矩阵)，a disp(b矩 阵：)，b disp(扭角矩阵：)，nfigure(3),plot程设计芝笋 0-告)，axis

14、(0 5 0 L6),title(馨长随叶高的分布)，xlabel(叶 S),ylabel(S长)figure(2),plot(r,n,o-),axis(0 5 0 35),title(扭角随叶高的分布)，xlabel(叶高 ),ylabel(扭角)figure(4),plot(jihe0(:,1),jihe0(:,2) ),title(s830翼型(单位弦长)，axis equal2.气动数据图表攻 角ClCdCmcl/cd-50.0397180.015793-0.135432.51484300.5937820.015621-0.1521638.0125951.0789440.021647-

15、0.1533349.84374101.4154420.039338-0.1413335.98126121.5191860.051391-0.1369229.56123141.5283040.073211-0.131920.87531161.4978520.103856-0.1312614.42241181.4464130.14123-0.135410.24151201.3718070.180991-0.139887.579425Cl、Cd攻角Cl CdCl/Cd攻角3.翼型图S03OW型（里位弦长）.数据总表|叶高|弦长|修正后|诱导|诱导|扭角/m/m弦长/m因子a因子bn/0.73981.08950.81140.36980.275528.82010.98640.90250.73420.29030.151525.13231.47960.65530