课程设计风力机.doc

《课程设计风力机.doc》由会员分享，可在线阅读，更多相关《课程设计风力机.doc（31页珍藏版）》请在三一办公上搜索。

1、课程设计风力机 作者： 日期：机械 与动 力 工 程学院风力机空气动力学课程设计设计题目： 小型三叶片风力机叶片设计 设计 人: 王 伦 班 级: 风能1101 组 号： 4 指导教师: 姚桂焕 设计时间： 2周 成 绩: 日期：21.23-204. 设计内容及要求1设计目的:风力机空气动力学课程设计是风能与动力工程专业中重要的实践性教学环节。通过该教学环节,使学生熟练掌握风力机叶片工作原理，并能够通过结合动量叶素理论相关知识与给定的环境条件设计出工作叶片，巩固和提高其风力机叶片设计及制造知识,树立其理论知识指导设计的工作思想,加深其对现场生产实际的了解，培养其对工程技术问题严肃认真、负责的态

2、度,为其以后从事实际工作打下坚实的基础。2.设计内容和要求：一、设计内容: （1）基于叶素和动量理论设计水平轴风力机叶片; （2）绘制风力机叶片弦长随叶片展向长度的变化曲线； （3)绘制风力机叶片扭角随叶片展向长度的变化曲线； （4）绘制设计风力机的性能曲线； (5）编写设计说明书,并附上必要的计算公式 （6）分组分别设计出不同翼型的风力机叶片。二、 设计要求: (1)掌握风力机叶片设计原理; （）掌握风力机叶片设计过程; (3）完成确定风力机叶片的参数； （)按时提交课程设计说明书、图纸,按时参加答辩。3.设计工作任务及工作量的要求一、对于给定的风力机工作环境以及功率按照ACA给定的参数设计

3、风力机工作叶片 每组同学数据给定情况不同可以采用手算和计算机编程序两种方式进行,在答辩时要提前说明计算方法以及在每种方法中遇到的问题，经教师校验方法正确后方可进行答辩。二、编制课程设计计算说明书设计计算说明书中应附上主要计算公式以及适用条件、工作原理、设计方法、系统构成及流程、计算成立条件,字数不少于0000字（至少要8000字),要求条理清晰，逻辑严密，字迹工整。.主要参考文献：风力机空气动力学.artin O.L anen（著）,肖劲松(译). 中国电力出版社，2010风工程与工业空气动力学. 贺德馨 国防工业出版社,06风力机设计理论及方法 赵丹平 北京大学出版社, 2015所用基础理论

4、: 根据动量理论，描述作用在风轮上的力与来流速度之间的关系。 根据叶素理论，将风轮叶片沿展向分成许多微段,即叶素，并假设在每个叶素上作用的气流相互之间没有干扰,作用在叶片上的力可分解为升力和阻力。 叶素-动量理论，假设各个叶素单元作用相互独立,各个圆环之间没有径向干扰，轴向诱导因子并不沿着径向方向改变。6.相关参数的选定:（1)调速方式：变桨距（)翼型 (AA系列)（3)叶片数目：3片(4）风力机功率：10 kW（5)设计风速:1 m/s(6)风轮转速:160rpm（）风能利用系数：05（）传动效率：09(9)发电机效率：0.9（0）叶片材料:碳纤维复合材料.常用变量参数声明;Pu-额定功率P

5、-设计计算功率-风流经风轮损失功率D-风轮直径-叶片长度B-叶片数C-叶片剖面弦长-风轮实度-风轮半径-叶素距旋转中心的距离-风轮扫过的面积1-额定风速，来流风速2-风流过风轮后的下风向风速-叶素攻角,最佳攻角C-升力系数Cd-阻力系数-叶素入流角-叶素扭角,桨距角 -rR叶素到旋转中心距离与风轮半径的比值 C-法向力系数 Cy-切向力系数 -风轮旋转角速度，额定旋转速度 -叶素旋转速度与主流风速比 。-叶尖速比 W-叶素合成流速 -叶片弦长实度 M-叶片转矩 T-风轮轴向推力 a-轴向诱导因子 b-周向诱导因子 第一章 风力机发展程 风能作为一种清洁的可再生能源，越来越受到世界各国的重视。其

6、蕴量巨大,全球的风能约为209MW,其中可利用的风能为2107MW，比地球上可开发利用的水能总量还要大10倍。风很早就被人们利用-主要是通过风车来抽水、磨面等,而现在，人们感兴趣的是如何利用风来发电。把风的动能转变成机械动能，再把机械能转化为电力动能，这就是风力发电。风力发电的原理,是利用风力带动风车叶片旋转，再透过增速机将旋转的速度提升,来促使发电机发电。依据目前的风车技术,大约是每秒三米的微风速度(微风的程度），便可以开始发电。风力发电正在世界上形成一股热潮，因为风力发电不需要使用燃料,也不会产生辐射或空气污染。风力机简介 风力机，将风能转换为机械功的动力机械,又称风车。广义地说,它是一种

7、以太阳为热源,以大气为工作介质的热能利用发动机。 许多世纪以来，它同水力机械一样，作为动力源替代人力、畜力,对生产力的发展发挥过重要作用。近代机电动力的广泛应用以及20世纪5年代中东油田的发现，使风力机的发展缓慢下来。70年代初期,由于“石油危机”，出现了能源紧张的问题，人们认识到常规矿物能源供应的不稳定性和有限性,于是寻求清洁的可再生能源遂成为现代世界的一个重要课题。风能作为可再生的、无污染的自然能源又重新引起了人们重视。 1.2风力机简史 风车最早出现在波斯,起初是立轴翼板式风车,后又发明了水平轴风车。风车传入欧洲后，15世纪在欧洲已得到广泛应用。荷兰、比利时等国为排水建造了功率达66千瓦

8、(90马力）以上的风车。8世纪末期以来,随着工业技术的发展,风车的结构和性能都有了很大提高，已能采用手控和机械式自控机构改变叶片桨距来调节风轮转速。 风力机用于发电的设想始于890年丹麦的一项风力发电计划。到1918年，丹麦已拥有风力发电机120台，额定功率为5千瓦不等。第一次世界大战后，制造飞机螺旋桨的先进技术和近代气体动力学理论为风轮叶片的设计创造了条件,于是出现了现代高速风力机。1931年，苏联采用螺旋桨式叶片建造了一台大型风力发电机,风速为1.5米/秒时,输出功率达10千瓦,风能利用系数提高到0.32。在第二次世界大战前后,由于能源需求量大,欧洲一些国家和美国相继建造了一批大型风力发电

9、机。41年,美国建造了一台双叶片、风轮直径达53.3米的风力发电机,当风速为3.4米秒时输出功率达125千瓦。 英国在0年代建造了三台功率为100千瓦的风力发电机。其中一台结构颇为独特,它由一个26米高的空心塔和一个直径244米的翼尖开孔的风轮组成。风轮转动时造成的压力差迫使空气从塔底部的通气孔进入塔内,穿过塔中的空气涡轮再从翼尖通气孔溢出。法国在50年代末到6年代中期相继建造了三台功率分别为1000千瓦和800千瓦的大型风力发电机。 现代的风力机具有增强的抗风暴能力,风轮叶片广泛采用轻质材料,运用近代航空气体动力学成就,使风能利用系数提高到0.4左右,用微处理机控制,使风力机保持在最佳运行状

10、态,发展了风力机阵列系统,风轮结构形式多样化。 法国人在2年代发明的垂直轴风轮在淹没了半个多世纪之后，已成为最有希望的风力机型之一。这种结构有多种形式，它具有运转速度高、效率高和传动机构简单等优点,但需用辅助装置起动。人们还提出了许多新的设想,如旋涡集能式风力机，据估计这种系统的单机功率将001000倍于常规风力机。 中国利用风车的历史至少不晚于13世纪中叶，曾建造了各种形式的简易风车碾米磨面、提水灌溉和制盐。直到20世纪50年代仍可见到“走马灯” 式风车。 1.3风力机的特点 新一代风力机的特点是：增强抗风暴能力;风轮叶片广泛采用轻质材料,如玻璃纤维复合材料等；运用近代航空气体动力学成就使风

11、能利用系数提高到0.45左右;用微处理机控制，使风力机保持在最佳运行状态;发展风力机 阵列系统;风轮结构形式多样化。14风力机的基本原理 太阳对大气层的不均匀照射和地球表面吸热能力的不同，在大气层中引起冷热空气的强烈对流而形成风。风的动能与风速的3次方成正比。用v表示空气速度，用表示质量密度,则单位时间内流过风轮扫掠面积的空气质量(m)为Av ,于是空气动能便是。由于气体的可压缩性，气体质点穿过风轮扫掠面能量转换界面时,风速由v1降为v2,即v。因自然风速v1只能有一部分被利用,若以风能利用系数表示利用程度,则可利用风能为 ,其中C1。根据气体动量 理论推导出风能利用系数的最大可能值为0.3，

12、因此风轮输出功率与风轮的工作面积成正比。C取决于风轮和叶片的结构和工艺。旧式风车C.1,现代风力机00.4，最高可达0.5。另外，现代风力机在能 量传输过程中大约还要损失13 理论上应输出的功。 5风力机的构成和分类 风力机的主要部件是风能接收装置。一般说来,凡在气流中产生不对称力的物理构形都能成为风能接收装置,它以旋转、平移或摆动运动而发出机械功。各类风能接收装置的取舍取决于使用寿命和成本的综合效益。风力机大都按风能接收装置的结构形式和空间布置来分类,一般分为水平轴结构和垂直轴结构两类。以风轮作为风能接收装置的常规风力机，按风轮转轴相对于气流方向的布置分为水平轴风轮式（转轴平行于气流方向）、

13、侧风水平轴风轮式（转轴平行于地面、垂直于气流方向)和垂直轴风轮式(转轴同时垂直于地面和气流方向）。广义风力机还包括那些利用风力产生平移运动的装置,如风帆船和中国古代的加帆手推车等。无论何种类型的风力机,都是由风能接收装置、控制机构、传动和支承部件等组成的。近代风力机还包括发电、蓄能等配套系统。 1.6风力机存在的问题 世界上已有数万台风力机在运行，作为辅助能源正在发挥作用。但风力机仍存在若干不足之处:能量输出不稳定，特别是大型风力机的利用率低，作为独 立能源的条件还不具备；安全可靠性尚无充分保障;成本在短期内尚不足以与矿物燃料相竞争。但是,随着人类对能源需求量的日益增多和科学技术的发展,上述问

14、题终会得到解决。.本课题的背景目的及主要工作 我国可开发利用的风能资源为25亿k,新疆、内蒙至东北和东南沿海两大主风带有有效风力时间百分率在70%以上。可以说,我国开发风能具有良好的自然环境和资源条件。近几年来,随着我国电网覆盖程度的提高,在各级政府、电力部门和国外政府及金融组织的援助下,我国在新疆、内蒙、广东、福建、辽宁等地区建立了20座风力发电场,总装机容量达302M,对缓解当地电力供应矛盾，提高供电质量起到了很好的作用。风力发电场的建设,加速了我国能源结构改革的进程，风能己成为真正的补充能源和发挥规模效益的生力军。 我国风力发电起步较晚,但发展较快。自0年末引进大型风力发电机以来,经过十

15、多年的不断引进、消化、吸收，积累了一定的经验。我国并网型风力发电技术在年代中期开始进行试验、示范,经过二十多年的努力,为今后进行国产化风力发电机组的规模化生产打下了一定的基础，同时也为推动国家风电产业化进程做出了努力。但遗憾的是,作为世界上的风能大国,我国尚不具备独立开发风力机尤其是大型风力机的能力，迄今为止国内已投入运行的风力机绝大部分是进口风力机。设计水平是主要制约因素，与此相关的基础研究、实验研究和新技术应用等方面与国外存在着较大的差距,有些领域国内甚至是空白。尤其是目前主流的大型风力机，我国基本上是依靠从国外引进生产技术来仿制。这不但受到成本、运输、售前售后等方面的制约，还要消耗大量的

16、资金，而且将使我国对风力机组的研制水平日益落后于国际先进水平，从根本上来说不利于我国风电产业的发展。更何况从国外引进的风机由于在设计时针对国外的风况和有一些特殊的环保要求,并不能和国内的情况非常吻合，不能很好地达到预期的性能。因此,必须以提高我国风力机的设计和研究水平为目标来实现“国产化。1.8本课程的主要工作：1. 阅读理解“风能转换原理与技术中”叶片设计的相关原理及设计知识，主要是第五章和第六章相关知识。 2. 编制叶素轴向、周向速度诱导因子、最佳弦长及扭角的计算的界面程序；3. 根据程序计算并绘制风力机叶片弦长随叶片展向长度的变化曲线；4.根据程序计算并绘制风力机叶片扭角随叶片展向长度的

17、变化曲线； 第三章风力机叶片设计一选定翼型,确定最佳攻角下的l和Cd 1.翼型(定义=r/,其中R为叶片设计半径，为叶片剖面到旋转中心距离) （1）叶片根部选用ACA418标准翼型（08） alculated plar fr: NA 41 1 1 ey numbe fixe ac number fx xtrf = 1.0 （o) 1.00（bto) Mc = .000 e= 1.0 Ncrt = 9.000 alp L CD CD M opXr Bot_tr 4.500 0.74 0.008 0.08-0.10 0.273 1.0000 4.70 .99 0.0778 .002970.101

18、0.4110 1.000 5.00 1.0254 .079 0.0011 -998 3979 1.0000 .250 1018 0.008 0.0326 -00995 0.381 10000 .50 1.077 0.0834 0.034 .092 0.731 1.0 5.50 1.131 0.00857 0.0035 -0.88 075 1.000 6.00 .120 0.0088 .0037 -.0983 0.398 1.00 根据以上资料,NC412的最佳攻角为5.25，此时对应的升力系数l=1018,阻力系数Cd=0.813二.设定叶尖速比,U,Pu,计算U 1.设定叶尖速比： 根据设定

19、参数,可求得=R/V=6.072. 设定=10m，Pu=10kW 3.计算： 根据风力机的额定功率计算公式:=.125Cpik 计算得：Dmin=22 根据风轮动量理论,功率关系有以下结论: U=.(1+2)- /（Cp*k)- 在单位时间内，有P=0.5A(-）- 1U- A=- 将已知设定参数：U=1m/s，Pu=0W,p=0.45,=.92, k=0.95, =.25k/m3代入以上各式得： U2=580m/ U=70ms三 计算叶片半径:R2=.65四 计算角速度： =2n/6=6.747rad五 将R分为1份，分别取=0,.2，1情况下的a和 由于风力机采用变桨距调节，当风速为额定风

20、速U1=m时,控制系统会调节桨距角,使得叶片处于最佳攻角5.75,此时对应升力系数Cl=1.061, 阻力系数d=.00882 入流角arctan r=R(.1,0.2,.0) =-计算流程:1.叶素弦长： 2.假设，b的初值，一般可取3.计算来流角:=rtan.计算各截面扭角：=-5.根据空气动力特性曲线得到叶素得升力系数Cl和阻力系数Cd; 6. 计算法向系数Cx和切向系数Cy C=cs+Cdsin Csin+Cdcos7. 计算新的a和b值: a b8. 比较新的a和b与上一次的和b，如果误差小于0.,则迭代 终止；否则回到3继续迭代。 由于需要重复计算,因此用计算机编程计算较为方便。

21、将以上步骤及公式编写V程序代码,计算结果如下:R弦长/扭角轴向诱导轴向诱导b0.1.73.7205120.4947.20.7674860.20820.0870.54217.330.520.059040.16012.470.1890.050.50.3370920186900150.60.22960860090.24355.250.846.070.0233.960.140.000090.191320.3600471.00.1742.0.1830038弦长曲线图:扭角曲线图：六.功率校核: （1）假设值，计算值；(以前面的,b为参考) (2)将假设值代入公式计算W,Sin,cos，Cx, Cy, （

22、3)通过以上值计算a和b (4)计算值与假设值比较，得到各自的误差, 若误差大于%，回到（）重新假设 (）计算每个下的dM,dP,将d相加得到P ()计算P与Pu的误差,若误差大于15%,返回（)重新假设 通过代码计算： 风轮叶片尖部的切向诱导速度可表示为: 转矩可表示为： 功率可表示为: 功率系数为： 校正后计算结果：010.20.3040.5功率19C0.50.52005.960.0.80.91.0功率1846108331823Cp0.40490480.40.8七绘制p曲线:（1） 由以上计算确定,R, 设为已知（2） 分别取=11，在每个下假设a的值（3） 通过迭代法求出诱导因子a,b，

23、使得假设值与计算值相差小于0.0（4） 计算此时a值和对应值下的对应Cp（5） 通过计算绘制C曲线通过程序输出计算结果：4.55.6.5Cp0.20.230.0.40430.470.4800577.588.59.501112p50.50.490.480410.60.30.30.2 C-曲线; 八. 载荷计算: 叶片主要载荷集中在根部,因此只需对叶片根部载荷进行校核 从叶片根部1%处取圆柱形与轮毂连接的叶根;NACA4418叶片展向10%处厚度： H=c18%=1.12718%=0.23m 叶片转动时根部圆最小直径为: 静止时最小直径： 应力:推力: 转矩： 式中： 查找相关资料，叶片所用材料碳

24、纤维复合材料的许用应力约为： 将许用应力代入应力表达式计算得到：m=1.9762 T=164.1Nm 根部直径： （D）in=092m=0.203 为了确保叶片根部安全：取叶根圆柱直径=0.30 第四章 总结为期两周的课程设计结束了。在这段时间的学习中，我收获了很多，也找到了自己身上的不足。这次的课程设计，不仅检验了之前所学的知识,也培养了我们如何去规划一件事情,如何去做一件事情,又如何完成一件事情。在设计过程中,与同学分工合作,相互探讨、相互学习、相互监督,不仅仅学会了协作,学会了宽容,更学会了理解。 风电作为一个有富有特色的专业,风力机组是我们要深入学习的，从叶片设计的的资料收集到叶片弦长

25、的计算，再到扭角计算程序的编写;从叶片叶素的绘制到叶片三维仿真的制作，我深入地了解到了叶片从计算到制作的各个过程,而原本学到的风力机空气动力学知识,也在这次的课程设计中得到实践。通过这次课程设计,让我们将理论联系到实践，开始真正地融入新能源这个行业。 附录：计算程序代码（VB） Optionxplici 通用声明im As Sngleim v s Snge, n As g, C SingeDim s Single, r As Single, h As SingleCL ASnle, CdAs SingleDima sSngl, bAs Siglim i sIteger,j A IneerDim

26、 aca As SngleDim nmAs ItegeCost = .141926-Pvte Sub Command1_Click（） 计算风轮直径 p = Text1.Txt: v=ext2.Txt:= ext3.Text: C = Text4.Tex: et5.Text r =Sqr( p /(p * Cp *0.9 05 *v 3) / 2 Txt6Tex Formt（2 r,000) 输出直径取四位小数 Sub-PivateSub Cmad2_Clk() Dim qr As Single, qAs inle,bASige, t s Single， c Sngl Dim ASigle,d

27、if Asil， if2 s nge p = Txt1ex:v =Tet2.ex： =ex3.Tex:C =ext4ext: =Text5Text r = Sr(8 p / (i * * Cp*0.92 *5 *v 3）) 2 qr = ext7.Text If r 0 Then CL =.0518： d = .013 Elsef qr = 0.8 An r 0.hn L 1.075: Cd 0.086 Elseqr 0.2 Thn = 1.009: Cd = 0.004 End I qr r c 16 * pi * r /(27 * CL *(pi n * r 30 v） * Sqr（(p*

28、 n r / 0 / v) 2 * (qr / )2+ 9) = ： b = 0 Do迭代法求诱导因子 t = 8 p *qr/(9 * ( * i qr/(3 v *(1 -a）) * C +(/ 3) * d) (9 *(n pi r (0 * v * (1 - ) 2 4 + ) (5) =At(v* （1 ） / (n*i* qr (1+ ) / ) q= ( t / (8 * iqr * (Si（)） ) （L *C(） +C* i()）- (3 * t/ （ * * qr *(Si（) 2)* （CL * Sin() - C Cs() 2) * (1- ) qb = (3 * t

29、(8 * pi qr * Si() Cos(）)* （CL Sn(）- Cd * Cos()* (1 + b) dif1 Abs( - qa): dif=As（b -q) = q:= LopUntil df1 0.1 And d2 0.001 限制精度 Tex8Text Format(,0.00) 轴向诱导因子 Tx9.Text = Fora（， 0.000) 周向诱导因子 T0.Tx =Format（, 0.00)弦长 = *（180 / pi） 计算扭角 I L = 1099 Then = - 75 EeIf CL 1.018 Then -.25 Esef L 157Th -5.5 End

30、If Text1.Txt = Fmat(, 0.00)End ubPrivat u omn3_Click（）Di1 As Snge, As Single,bt AsSile, m As ne, cp1 A Single Dm dif AsSingl p = Tex.Tex： =Tx2Tex: n = Tet3.Txt: = Text5.Text: cp= Text4.ext r S(8 p /（pi * * 0.92 * 0.9 v 3)） / 2 a Te8.et bt =0 * v 2 a* (1 a)( 2 * n * 2* (1 + bt)） c =4 a （1 - a) 2/ ( bt） p1 05 * * pi * 2 * 3*cp1 0.92 .95 p1 /(i /30) Text1.Text = p TetTet omat(cp1, 0.0) Txt14.