专业课程设计风力机.doc

专业课程设计风力机机 械 与 动 力 工 程 学 院风力机空气动力学课程设计 设计题目： 小型三叶片风力机叶片设计 设 计 人： 王 伦 班 级： 风能1101 组 号： 4 指导教师： 姚桂焕 设计时间： 2周 成 绩： 日期：2014.6.23-2014.7.4 - 2 -设计内容及要求1设计目的：风力机空气动力学课程设计是风能与动力工程专业中重要的实践性教学环节。通过该教学环节，使学生熟练掌握风力机叶片工作原理，并能够通过结合动量叶素理论相关知识与给定的环境条件设计出工作叶片，巩固和提高其风力机叶片设计及制造知识，树立其理论知识指导设计的工作思想，加深其对现场生产实际的了解，培养其对工程技术问题严肃认真、负责的态度，为其以后从事实际工作打下坚实的基础。2设计内容和要求：一、设计内容： （1）基于叶素和动量理论设计水平轴风力机叶片； （2）绘制风力机叶片弦长随叶片展向长度的变化曲线； （3）绘制风力机叶片扭角随叶片展向长度的变化曲线； （4）绘制设计风力机的性能曲线； （5）编写设计说明书，并附上必要的计算公式 （6）分组分别设计出不同翼型的风力机叶片。二、 设计要求： （1）掌握风力机叶片设计原理； （2）掌握风力机叶片设计过程； （3）完成确定风力机叶片的参数； （4）按时提交课程设计说明书、图纸，按时参加答辩。3设计工作任务及工作量的要求一、对于给定的风力机工作环境以及功率按照NACA给定的参数设计风力机工作叶片 每组同学数据给定情况不同可以采用手算和计算机编程序两种方式进行，在答辩时要提前说明计算方法以及在每种方法中遇到的问题，经教师校验方法正确后方可进行答辩。二、编制课程设计计算说明书设计计算说明书中应附上主要计算公式以及适用条件、工作原理、设计方法、系统构成及流程、计算成立条件，字数不少于10000字（至少要8000字），要求条理清晰，逻辑严密，字迹工整。4主要参考文献：风力机空气动力学. Martin O.L. Hansen（著），肖劲松（译）. 中国电力出版社，2010风工程与工业空气动力学. 贺德馨. 国防工业出版社，2006风力机设计理论及方法. 赵丹平. 北京大学出版社, 20125所用基础理论： 根据动量理论，描述作用在风轮上的力与来流速度之间的关系。 根据叶素理论，将风轮叶片沿展向分成许多微段，即叶素，并假设在每个叶素上作用的气流相互之间没有干扰，作用在叶片上的力可分解为升力和阻力。 叶素-动量理论，假设各个叶素单元作用相互独立，各个圆环之间没有径向干扰，轴向诱导因子a并不沿着径向方向改变。6相关参数的选定：（1）调速方式：变桨距（2）翼型 （NACA系列）（3）叶片数目：3片（4）风力机功率：10 kW（5）设计风速：10 m/s（6）风轮转速：160rpm（7）风能利用系数：0.45（8）传动效率：0.92（9）发电机效率：0.95（10）叶片材料：碳纤维复合材料7.常用变量参数声明;Pu-额定功率P-设计计算功率P-风流经风轮损失功率D-风轮直径L-叶片长度B-叶片数C-叶片剖面弦长-风轮实度R-风轮半径r-叶素距旋转中心的距离A-风轮扫过的面积U1-额定风速，来流风速U2-风流过风轮后的下风向风速-叶素攻角，最佳攻角Cl-升力系数Cd-阻力系数-叶素入流角-叶素扭角，桨距角 -r/R叶素到旋转中心距离与风轮半径的比值 Cx-法向力系数 Cy-切向力系数 -风轮旋转角速度，额定旋转速度 -叶素旋转速度与主流风速比 。-叶尖速比 W-叶素合成流速 r-叶片弦长实度 M-叶片转矩 T-风轮轴向推力 a-轴向诱导因子 b-周向诱导因子 第一章 风力机发展程 风能作为一种清洁的可再生能源，越来越受到世界各国的重视。其蕴量巨大，全球的风能约为2.74109MW，其中可利用的风能为2107MW，比地球上可开发利用的水能总量还要大10倍。风很早就被人们利用-主要是通过风车来抽水、磨面等，而现在，人们感兴趣的是如何利用风来发电。 把风的动能转变成机械动能，再把机械能转化为电力动能，这就是风力发电。风力发电的原理，是利用风力带动风车叶片旋转，再透过增速机将旋转的速度提升，来促使发电机发电。依据目前的风车技术，大约是每秒三米的微风速度（微风的程度），便可以开始发电。风力发电正在世界上形成一股热潮，因为风力发电不需要使用燃料，也不会产生辐射或空气污染。1.1风力机简介 风力机，将风能转换为机械功的动力机械，又称风车。广义地说，它是一种以太阳为热源，以大气为工作介质的热能利用发动机。 许多世纪以来，它同水力机械一样，作为动力源替代人力、畜力，对生产力的发展发挥过重要作用。近代机电动力的广泛应用以及20世纪50年代中东油田的发现，使风力机的发展缓慢下来。70年代初期，由于“石油危机”，出现了能源紧张的问题，人们认识到常规矿物能源供应的不稳定性和有限性，于是寻求清洁的可再生能源遂成为现代世界的一个重要课题。风能作为可再生的、无污染的 自然能源又重新引起了人们重视。 1.2风力机简史 风车最早出现在波斯，起初是立轴翼板式风车，后又发明了水平轴风车。风车传入欧洲后，15世纪在欧洲已得到广泛应用。荷兰、比利时等国为排水建造了功率达66千瓦（90马力）以上的风车。18世纪末期以来，随着工业技术的发展，风车的结构和性能都有了很大提高，已能采用手控和机械式自控机构改变叶片桨距来调节风轮转速。 风力机用于发电的设想始于1890年丹麦的一项风力发电计划。到1918年，丹麦已拥有风力发电机120台，额定功率为525千瓦不等。第一次世界大战后，制造飞机螺旋桨的先进技术和近代气体动力学理论为风轮叶片的设计创造了条件，于是出现了现代高速风力机。1931年，苏联采用螺旋桨式叶片建造了一台大型风力发电机，风速为13.5米/秒时,输出功率达100千瓦,风能利用系数提高到0.32。 在第二次世界大战前后，由于能源需求量大，欧洲一些国家和美国相继建造了一批大型风力发电机。1941年，美国建造了一台双叶片、风轮直径达53.3米的风力发电机，当风速为13.4米/秒时输出功率达1250千瓦。 英国在50年代建造了三台功率为100千瓦的风力发电机。其中一台结构颇为独特，它由一个26米高的空心塔和一个直径24.4米的翼尖开孔的风轮组成。风轮转动时造成的压力差迫使空气从塔底部的通气孔进入塔内，穿过塔中的空气涡轮再从翼尖通气孔溢出。法国在50年代末到60年代中期相继建造了三台功率分别为1000千瓦和800千瓦的大型风力发电机。 现代的风力机具有增强的抗风暴能力，风轮叶片广泛采用轻质材料，运用近代航空气体动力学成就，使风能利用系数提高到0.45左右，用微处理机控制，使风力机保持在最佳运行状态，发展了风力机阵列系统，风轮结构形式多样化。 法国人在20年代发明的垂直轴风轮在淹没了半个多世纪之后，已成为最有希望的风力机型之一。这种结构有多种形式，它具有运转速度高、效率高和传动机构简单等优点，但需用辅助装置起动。人们还提出了许多新的设想，如旋涡集能式风力机，据估计这种系统的单机功率将1001000倍于常规风力机。 中国利用风车的历史至少不晚于13世纪中叶，曾建造了各种形式的简易风车碾米磨面、提水灌溉和制盐。直到20世纪50年代仍可见到“走马灯” 式风车。 1.3风力机的特点 新一代风力机的特点是：增强抗风暴能力；风轮叶片广泛采用轻质材料,如玻璃纤维复合材料等；运用近代航空气体动力学成就使风能利用系数提高到0.45左 右；用微处理机控制，使风力机保持在最佳运行状态；发展风力机 阵列系统；风轮结构形式多样化。 1.4风力机的基本原理 太阳对大气层的不均匀照射和地球表面吸热能力的不同，在大气层中引起冷热空气的强烈对流而形成风。风的动能与风速的3次方成正比。用v表示空气速度，用表示质量密度，则单位时间内流过风轮扫掠面积A的空气质量(m)为 Av ，于是空气动能便是。 由于气体的可压缩性，气体质点穿过风轮扫掠面能量转换界面时,风速由v1降为v2，即v1v2。因自然风速v1只能有一部分被利用，若以风能利用系数C表示利用程度，则可利用风能为 ，其中C1。根据气体动量 理论推导出风能利用系数的最大可能值为0.593，因此风轮输出功率与风轮的工作面积成正比。C取决于风轮和叶片的结构和工艺。旧式风车 C0.10，现代风力机C0.30.4，最高可达0.5。另外，现代风力机在能 量传输过程中大约还要损失1/3 理论上应输出的功。 1.5风力机的构成和分类 风力机的主要部件是风能接收装置。一般说来，凡在气流中产生不对称力的物理构形都能成为风能接收装置，它以旋转、平移或摆动运动而发出机械功。各类风能接收装置的取舍取决于使用寿命和成本的综合效益。风力机大都按风能接收装置的结构形式和空间布置来分类，一般分为水平轴结构和垂直轴结构两类。以风轮作为风能接收装置的常规风力机，按风轮转轴相对于气流方向的布置分为水平轴风轮式（转轴平行于气流方向）、侧风水平轴风轮式（转轴平行于地面、垂直于气流方向）和垂直轴风轮式（转轴同时垂直于地面和气流方向）。广义风力机还包括那些利用风力产生平移运动的装置，如风帆船和中国古代的加帆手推车等。无论何种类型的风力机，都是由风能接收装置、控制机构、传动和支承部件等组成的。近代风力机还包括发电、蓄能等配套系统。 1.6风力机存在的问题 世界上已有数万台风力机在运行，作为辅助能源正在发挥作用。但风力机仍存在若干不足之处：能量输出不稳定，特别是大型风力机的利用率低，作为独 立能源的条件还不具备；安全可靠性尚无充分保障；成本在短期内尚不足以与矿物燃料相竞争。但是，随着人类对能源需求量的日益增多和科学技术的发展，上述问题终会得到解决。 1.7本课题的背景目的及主要工作 我国可开发利用的风能资源为253亿kW，新疆、内蒙至东北和东南沿海两大主风带有有效风力时间百分率在70以上。可以说，我国开发风能具有良好的自然环境和资源条件。近几年来，随着我国电网覆盖程度的提高，在各级政府、电力部门和国外政府及金融组织的援助下，我国在新疆、内蒙、广东、福建、辽宁等地区建立了20座风力发电场，总装机容量达302MW，对缓解当地电力供应矛盾，提高供电质量起到了很好的作用。风力发电场的建设，加速了我国能源结构改革的进程，风能己成为真正的补充能源和发挥规模效益的生力军。 我国风力发电起步较晚，但发展较快。自80年末引进大型风力发电机以来，经过十多年的不断引进、消化、吸收，积累了一定的经验。我国并网型风力发电技术在80年代中期开始进行试验、示范，经过二十多年的努力，为今后进行国产化风力发电机组的规模化生产打下了一定的基础，同时也为推动国家风电产业化进程做出了努力。 但遗憾的是，作为世界上的风能大国，我国尚不具备独立开发风力机尤其是大型风力机的能力，迄今为止国内已投入运行的风力机绝大部分是进口风力机。设计水平是主要制约因素，与此相关的基础研究、实验研究和新技术应用等方面与国外存在着较大的差距，有些领域国内甚至是空白。尤其是目前主流的大型风力机，我国基本上是依靠从国外引进生产技术来仿制。这不但受到成本、运输、售前售后等方面的制约，还要消耗大量的资金，而且将使我国对风力机组的研制水平日益落后于国际先进水平，从根本上来说不利于我国风电产业的发展。更何况从国外引进的风机由于在设计时针对国外的风况和有一些特殊的环保要求，并不能和国内的情况非常吻合，不能很好地达到预期的性能。因此，必须以提高我国风力机的设计和研究水平为目标来实现“国产化。1.8本课程的主要工作： 1. 阅读理解“风能转换原理与技术中”叶片设计的相关原理及设计知识，主要是第五章和第六章相关知识。 2. 编制叶素轴向、周向速度诱导因子、最佳弦长及扭角的计算的界面程序；3. 根据程序计算并绘制风力机叶片弦长随叶片展向长度的变化曲线；4.根据程序计算并绘制风力机叶片扭角随叶片展向长度的变化曲线；- 14 - 第三章 风力机叶片设计一选定翼型，确定最佳攻角下的Cl和Cd 1.翼型（定义=r/R，其中R为叶片设计半径，r为叶片剖面到旋转中心距离） (1)叶片根部选用NACA4418标准翼型（0.8） Calculated polar for: NACA 4412 1 1 Reynolds number fixed Mach number fixed xtrf = 1.000 (top) 1.000 (bottom) Mach = 0.000 Re = 1.000 e 6 Ncrit = 9.000 alpha CL CD CDp CM Top_Xtr Bot_Xtr 4.500 0.9734 0.00758 0.00284 -0.1004 0.4273 1.0000 4.750 0.9993 0.00778 0.00297 -0.1001 0.4110 1.0000 5.000 1.0254 0.00797 0.00311 -0.0998 0.3979 1.0000 5.250 1.0518 0.00813 0.00326 -0.0995 0.3861 1.0000 5.500 1.0777 0.00834 0.00342 -0.0992 0.3731 1.0000 5.750 1.1031 0.00857 0.00359 -0.0988 0.3575 1.0000 6.000 1.1280 0.00884 0.00379 -0.0983 0.3398 1.0000 根据以上资料，NACA4412的最佳攻角为5.25，此时对应的升力系数Cl=1.0518,阻力系数Cd=0.0813二设定叶尖速比，U,Pu,计算U 1.设定叶尖速比： 根据设定参数，可求得=R/V=6.0872. 设定U=10m/s,Pu=10kW 3.计算U: 根据风力机的额定功率计算公式：Pu=0.125Cpik 计算得：Dmin=7.272m 根据风轮动量理论，功率关系有以下结论: U=0.5(U1+U2)- P=Pu/(Cp*i*k)- 在单位时间内，有P=0.5A(-)- U1=U- A=/4- 将已知设定参数：U=10m/s,Pu=10kW，Cp=0.45，i=0.92， k=0.95, =1.225kg/m3代入以上各式得： U2=5.580m/s U=7.790m/s三 计算叶片半径： R=D/2=3.635四 计算角速度： =2n/60=16.747rad/s五 将R分为10份，分别取=0.1,0.2，1情况下的a和 a 由于风力机采用变桨距调节，当风速为额定风速U1=10m/s时，控制系统会调节桨距角，使得叶片处于最佳攻角=5.75，此时对应升力系数Cl=1.0631, 阻力系数Cd=0.00882 入流角=arctan r=R(=0.1,0.2,.1.0) =-计算流程：1.叶素弦长: 2.假设a,b的初值，一般可取03.计算来流角：=arctan4.计算各截面扭角：=-5.根据空气动力特性曲线得到叶素得升力系数Cl和阻力系数Cd; 6. 计算法向系数Cx和切向系数Cy Cx=Clcos+Cdsin Cy=Clsin+Cdcos7. 计算新的a和b值： a= b=8. 比较新的a和b与上一次的a和b，如果误差小于0.001，则迭代 终止；否则回到3继续迭代。 由于需要重复计算，因此用计算机编程计算较为方便。 将以上步骤及公式编写VB程序代码，计算结果如下：=r/R弦长/m扭角/轴向诱导a轴向诱导b0.11.127933.720.25120.49470.20.756724.860.20820.10870.30.540217.330.19520.04590.40.416012.470.18980.02510.50.33709.210.18690.01580.60.28296.940.18560.01090.70.24355.250.18460.00790.80.21363.960.18400.00600.90.19133.200.18360.00471.00.17242.380.18330.0038弦长曲线图：扭角曲线图：六功率校核： （1）假设a值，计算b值;（以前面的a,b为参考） (2)将假设值代入公式计算W,Sin，cos，Cx, Cy, (3)通过以上值计算a和b (4)计算值与假设值比较，得到各自的误差， 若误差大于8%，回到（1）重新假设 （5）计算每个下的dM,dP,将dP相加得到P (6) 计算P与Pu的误差，若误差大于15%，返回（1）重新假设 通过代码计算： 风轮叶片尖部的切向诱导速度可表示为： 转矩可表示为： 功率可表示为： P=M 功率系数为： 校正后计算结果：0.10.20.30.40.5功率1155411551111901102810938Cp0.520.520.50.50.490.60.70.80.91.0功率1089710865108461083310823Cp0.490.490.480.480.48七绘制Cp曲线：（1） 由以上计算确定，R, 设为已知（2） 分别取=112，在每个下假设a的值（3） 通过迭代法求出诱导因子a，b,使得假设值与计算值相差小于 0.001（4） 计算此时a值和对应值下的对应Cp（5） 通过计算绘制Cp曲线通过程序输出计算结果：12344.555.566.5Cp0.120.230.310.40.430.470.480.50.577.588.599.5101112Cp0.50.50.490.480.410.360.330.310.28 Cp-曲线; 八. 载荷计算： 叶片主要载荷集中在根部，因此只需对叶片根部载荷进行校核 从叶片根部10%处取圆柱形与轮毂连接的叶根； NACA4418叶片展向10%处厚度： H=c18%=1.127918%=0.203m 叶片转动时根部圆最小直径为： 静止时最小直径： 应力：推力： 转矩： 式中： 查找相关资料，叶片所用材料碳纤维复合材料的许用应力约为： 将许用应力代入应力表达式计算得到：rm=1.9762 T=1654.631N*m 根部直径： （D）min=0.0932m 0.8 Then CL = 1.0518: Cd = 0.0813 ElseIf qr = 0.2 Then CL = 1.0575: Cd = 0.00886 ElseIf qr 0.2 Then CL = 1.0909: Cd = 0.0094 End If qr = qr * r c = 16 * pi * r / (27 * CL * (pi * n * r / 30 / v) * Sqr(pi * n * r / 30 / v) 2 * (qr / r) 2 + 4 / 9) a = 0: b = 0 Do 迭代法求诱导因子 t = 8 * pi * qr / (9 * (n * pi * qr / (30 * v * (1 - a) * CL + (2 / 3) * Cd) * (9 * (n * pi * qr / (30 * v * (1 - a) 2 / 4 + 1) (0.5) = Atn(v * (1 - a) / (n * pi * qr * (1 + b) / 30) qa = (3 * t / (8 * pi * qr * (Sin() 2) * (CL * Cos() + Cd * Sin() - (3 * t / (8 * pi * qr * (Sin() 2) * (CL * Sin() - Cd * Cos() 2) * (1 - a) qb = (3 * t / (8 * pi * qr * Sin() * Cos() * (CL * Sin() - Cd * Cos() * (1 + b) dif1 = Abs(a - qa): dif2 = Abs(b - qb) a = qa: b = qb Loop Until dif1 0.001 And dif2 0.001 限制精度 Text8.Text = Format(a, 0.0000) 轴向诱导因子 Text9.Text = Format(b, 0.0000) 周向诱导因子 Text10.Text = Format(c, 0.0000) 弦长 = * (180 / pi) 计算扭角 If CL = 1.0909 Then = - 5.75 ElseIf CL = 1.0518 Then = - 5.25 ElseIf CL = 1.0575 Then = - 5.5 End If Text11.Text = Format(, 0.00)End SubPrivate Sub Command3_Click() Dim p1 As Single, a As Single, bt As Single, m As Single, cp1 As Single Dim dif3 As Single p = Text1.Text: v = Text2.Text: n = Text3.Text: = Text5.Text: cp = Text4.Text r = Sqr(8 * p / (pi * * cp * 0.92 * 0.95 * v 3) / 2 a = Text8.Text bt = 900 * v 2 * a * (1 - a) / (pi 2 * n 2 * r 2 * (1 + bt) cp1 = 4 * a * (1 - a) 2 / (1 + bt) p1 = 0.5 * * pi * r 2 * v 3 * cp1 * 0.92 * 0.95 m = p1 / (pi * n / 30) Text12.Text = p1 Text13.Text = Format(cp1, 0.00) Text14.Text = mEnd SubPrivate Sub Command4_Click() Dim qr As Single, qa As Sin