《风力发电机组设计与制造》课程设计任务书

1、 课程设计课程设计（综合实验综合实验）报告报告 （2012-2013年度第二学期） 名称： 题目： 院系： 班级： 学号： 学生姓名： 指导教师： 设计周数： 成绩： 日期：2013年月日风力发电机组设计与制造课程设计任务书风力发电机组设计与制造课程设计任务书 一、设计内容一、设计内容 风电机组总体技术设计 二、目的与任务二、目的与任务 主要目的： 1. 以大型水平轴风力机为研究对象，掌握系统的总体设计方法； 2. 熟悉相关的工程设计软件； 3. 掌握科研报告的撰写方法。 主要任务： 1. 确定风电机组的总体技术参数； 2. 关键零部件（齿轮箱、发电机和变流器）技术参数； 3. 计算关键零部件

2、（叶片、风轮、主轴、连轴器和塔架等）载荷和技术参数； 4. 完成叶片设计任务； 5. 确定塔架的设计方案。 ! 6. 撰写一份课程设计报告。 三、主要内容三、主要内容 选择功率范围在至6MW之间的风电机组进行设计。 1）原始参数：风力机的安装场地50米髙度年平均风速为s,60米髙度年平均风速为s,70米髙度年平均风速为m/s，当地历史最大风速为49m/s，用户希望安装MW至6MW之间的风力机。采用63418翼型，63418翼型的升力系数、阻力系数数据如表1所示。空气密度设定为m3。 2）设计内容 （1）确定整机设计的技术参数。设定几种风力机的C曲线和C曲线，风力机基本参数包括叶片 pt 数、风

3、轮直径、额定风速、切入风速、切出风速、功率控制方式、传动系统、电气系统、制动系统形式和塔架髙度等，根据标准确定风力机等级； （2）关键部件气动载荷的计算。设定几种风轮的C曲线和C曲线，计算几种关键零部件的载荷 pt （叶片载荷、风轮载荷、主轴载荷、连轴器载荷和塔架载荷等）；根据载荷和功率确定所选定机型主要部件的技术参数（齿轮箱、发电机、变流器、连轴器、偏航和变桨距电机等）和型式。以上内容建议用计算机编程实现，确定整机和各部件（系统）的主要技术参数。 （3）塔架根部截面应力计算。计算暴风工况下风轮的气动推力，参考风电机组的整体设计参数，计算塔架根部截面的应力。最后提交有关的分析计算报告。 四、进

4、度计划四、进度计划 序号序号 设计（实验）内容设计（实验）内容 完成时间完成时间 备注备注 1 风电机组整体参数设计 天 2 风电机组气动特性初步计算 2天 3 机组及部件载荷计算 | 2天 4 齿轮箱、发电机、变流器技术参数 天 4 塔架根部截面应力计算 * 1天 5 报告撰写 天 6 课程设计答辩 天 五、设计（实验）成果要求五、设计（实验）成果要求 1提供设计的风电机组的性能计算结果； 2绘制整机总体布局工程图。 六、考核方式六、考核方式 提交一份课程设计报告；准备课程设计PPT，答辩。 风电机组总体技术设计风电机组总体技术设计 一、确定风电机组的总体技术参数一、确定风电机组的总体技术参

5、数 1、叶片数、叶片数B=3 ！ 风轮的叶片数取决于风轮尖速比，一般来说，要得到很大的输出扭矩就需要较大的叶片实度；现 代风力发电机组实度较少，一般只需要13个叶片。叶片数多的风力机在低尖速比运行时有较高的风能利用系数，既有较大的转矩，而且启动风速低，适用于提水。而叶片数少的风力机则在高尖速比运行时有较高的风能利用系数，但启动风速高，因此适用于风力发电。 3叶片风轮有以下优点：平衡简单，动载荷小；3叶片使风力发电机组系统运行平稳，基本消除了系统的周期载荷，输出较稳定的转矩；3叶片风轮通常能提供较佳的效率；3叶片风轮从审美的角度来说更令人满意；3叶片风轮的受力平衡好，轮毂可以简单些；噪声相对较小

6、。 综上所述，叶片数B选取3。 2、机组额定功率及种类、机组额定功率及种类 选取3MW级别的双馈变速恒频式风力发电机组 3、切入风速、切入风速 取切入风速Vin=3m/s。 4、切出风速、切出风速 取切出风速Vout=25m/s。 5、额定风速、额定风速 取额定风速Vr=13m/s。 6、各部分效率、各部分效率 传动系统效率ni取 发电机效率n2取 变流器效率n3取 7、风能利用系数风能利用系数C p 额定功率下风能利用系数Cp取 8、风轮直径风轮直径 由风力发电机组输出功率得叶片直径： 其中： P风力发电机组额定输出功率，取3000kW；r 空气密度(一般取标准大气状态)，取m3； V-额定

7、风速，取13m/s；r ( D风轮直径； n1-传动系统效率，取； n2发电机效率，取； n3-变流器效率，取； Cp-额定功率下风能利用系数,取。 p 9、功率、功率Pw曲线，风能利用系数曲线，风能利用系数C和推力系数和推力系数C曲线的确定曲线的确定 PT 风电机组的设计需要给出功率输出特性曲线。如果确定了风电机组的切入风速，额定风速和切出 风速，可以利用下面的方程计算输出功率，查阅相关资料，得到其风能利用系数随着风速的变化值， 得到下表。(其中所用公式为:Cp=4a(l-a)2Ct=4a(l-a) 8P 风速 功率PwkW 风能利用系数 轴向诱导因 推力系数CT叶尖速比 vm/s Cp 子

8、a 3 22. 4 ！ 16. 5 13. 6 11. 】 7 8 9 10 11 12 5. r pV3兀Cqqq rp123 利用上述的公式，由excel的函数运算，计算得到随着风速的变化，Cp、和Ct的值得变化 13 5. 8P =92m D= r pV3兀C耳qq rp123 D= 14 15 16 17 18 19 20 21 22 23 24 25 风力发电机组产生的气动噪声正比于叶尖速度的5次方。通常将路基风力发电机组的叶尖速度限 制在65m/s左右。这里取叶尖速V=65m/s。 由v=60得到n=60v/，将V=65m/s带入得，r/min。 11、叶尖速比入、叶尖速比入 由公

9、式，带入额定转速n=min可以得到 参考CP-入曲线，当Cp=时对应的入的值约为5,因此所选数值合理。 12、功率控制方式、功率控制方式 功率控制方式选取当今应用较为广泛的变浆距控制。这种机组当风速过高时，通过减小叶片翼型上合成气流与翼型几何弦的夹角（攻角），改变风力发电机组获得的空气动力转矩，能使功率输出保持稳定。同时，风机在起动过程也需要通过变浆距来获得足够的起动转矩。采用变浆距技术的风力发电机组还可以使叶片和整机的受力状况大为改善，这对于大型风力发电机组十分的有利。 13、传动系统、传动系统 传动系统选定由以下四部分组成：主轴系统、增速传动系统（齿轮箱）、轴系的支撑与连接（如轴承和联轴器

10、）、制动装置。 14、制动系统、制动系统 制动系统的形式以空气动力制动和机械制动为主。 15、电气系统、电气系统 风电机组的电气系统包括以下四个主要部分：变浆距驱动，偏航驱动，发电机部分，控制和监控部分 16、塔架高度、塔架高度 由于风速与距地面高度有关，增加塔架高度可使风轮获取更多的风能，但制造更高的塔架也需要 更多的材料，使其造价相应增加。大型机组的塔架高度H可以按下式初步确定：式中D为风轮直径，单位为 m。 塔架的高度参数的选择与地形和地貌有关，如陆地与海上风电机组会有所不同。陆地的地表性对粗糙且风速随高度变化缓慢，可能需要较高的塔架；而海平面较光滑，风速沿着高度方向变化梯度大。一般而论

11、，即使相同容量的机组，在近海安装的机组塔架高度相对较低。考虑地貌因素的塔架最低设计高度一般可用下式估算： H=h+C+R 式中：h机组附近的障碍物髙度，单位m； C障碍物最髙点到风轮扫掠面最低点的距离(最小取)，单位m。 综上考虑，塔架髙度H取90m 17、设计寿命、设计寿命 按照惯例，设计寿命取20年。 18、设计标准、设计标准 设计标准为IEC61400-1 19、风力机等级、风力机等级 70m年平均风速为s。去IECIIIA,参考下表: 风电机组等级 第I级 第II级 第III级 ) 第W级 参考风速Vr=5Va(m/s) 50 30 年平均风速Va(m/s) 10 6 湍流强度等级 A

12、 e15 a 2 2 2 2 B e15 b 3 3 3 3 阶段性总结阶段性总结 叶片数 3 叶尖速比 5 切入风速 3m/s 功率控制 变浆距控制 切出风速 25m/s 制动系统形式 ) 空气动力制动和机械制动 额定风速 13m/s 塔架髙度 90m 额定功率 3MW 设计寿命 20年 风能利用系数 设计标准 IEC61400 ！ 风轮转速 min 风力机等级 IECIIIA 风轮直径 92m / 、关键部件气动载荷的计算、关键部件气动载荷的计算 包括设定几种风轮的Cp曲线和Ct曲线，计算几种关键零部件的载荷（叶片载荷、风轮载荷、主轴载荷、连轴器载荷和塔架载荷等）；根据载荷和功率确定所选定

13、机型主要部件的技术参数（齿轮箱、发电机、变流器、连轴器、偏航和变桨距电机等）和型式。以上内容建议用计算机编程实现，确定整机和各部件（系统）的主要技术参数。最后提交有关的分析报告。 根据所提供的63418翼型的升力系数、阻力系数数据，得到攻角与升力系数、阻力系数的关系曲 线。 由图可以观察出，攻角约等于13时升力系数最大，大约为左右。攻角约等于9时，升力系数约为最大升力系数的倍，为，此时的阻力系数为。当攻角为5时升阻比最大，此时升力系数为。 在叶片处选择10段截面，分别在10%R,20%R,30%R90%R,100%R。由以下公式: “16兀R/r N= 9 a、给出r,攻角i已知，C已知; 攻

14、角与升力系数攻角与升力系数攻角与阻力系数攻角与阻力系数 S（R）2+9 b、由式叫3c如壬求得角; c、由式 k=：(九2+1)COS屮得到k; d、 I1-k2 由式h=,1+得到h； e、九2 1+h f、由式I=arcctg(九)得到I； 1+k g、a=I-i,C=8nr(h-1)cosI/BC(h+1)； L h、将a-g编程进行运算得到数据-表2。（程序见附录） 通过excel计算得到以上各个参数的数值： ！ 1、主要部件功率计算、主要部件功率计算 1 1）发电机）发电机 发电机类型：双馈异步变速恒频式发电机 额定功率：P=P/n*=3000/*=3402kW(预留10%的裕度)

15、发电机r3 极对数：2 由公式60f=np可以得到发电机的额定转速为1650m/min(预留10%的裕度) 2)2) 齿轮箱齿轮箱 齿轮箱类型：增速，2级行星齿轮+1级平行齿轮 齿轮箱效率 传动比：120 低速轴转速：r/min； 高速轴转速：1650r/min 齿轮箱功率：P=3581kw GB 3)3) 联轴器联轴器 ! 低速轴联轴器动率：3400kW 髙速轴联轴器动率：3500kW 4)4) 变流器变流器变流器的功率通常为风电机组的额定功率的1/21/3，考虑到风电机组的可靠性，通常为风电机组额定功率的1/2，即变流器的功率为1500kW。 5)5) 偏航电机偏航电机 偏航电机的功率由最

16、大偏航扭矩确定。在此先初步选定4台4kW的偏航电机。 2、主要部件的载荷计算、主要部件的载荷计算 叶片绕风轮轴旋转时，有离心力作用在叶片上。离心力的方向是自旋转中心沿半径向外。作用在 叶片上的离心力F可按照下式求解: c P叶片的密度，单位是kg/m3，取m3 y A叶素处的叶片截面面积，单位是m2 r Q风轮的角速度，单位是rad/s r叶片起始处的旋转半径 0 R叶片结束处的旋转半径 其中：Q二rad/s 1)1) 2)2) 叶片的载荷计算叶片的载荷计算 由matlab计算的：F=00N。 c 3)3) 作用在叶片上的风压力作用在叶片上的风压力 风压力是作用在叶片上沿风速方向的气动力。风轮

17、转动时的风压力： 设F的作用点几粒风轮轴的距离为r,则：Vm & 由matlab计算得：Fv二755580N,rm=。 4)4) 作用在叶片上的气动力矩作用在叶片上的气动力矩 M Mb b b M是使风轮转动的力矩，由matlab计算得：M=144250Nm。bb 5)5) 作用在叶片上的陀螺力矩作用在叶片上的陀螺力矩 M是风轮对风调向时产生的惯性力矩。当风向改变时，风轮除以角速度Q绕0X轴转动外，还要以k 角速度3绕0Z轴转动。 整个叶片的转动惯量为： 由matlab计算得，J=0kgm2 对风调向时，叶片可看成以角速度3绕0Z轴转动的牵连运动、以Q绕0X轴转动的相对运动的复 合运动

18、，这时要产生柯式加速度a：=2w*v k 这里是惯性半径，B是V与3的夹角。当叶片旋转到铅垂位置时，B=0。，柯式加速度达到最大r 值： 由于柯式加速度，就产生了柯式角加速度 由动量矩定理知，叶片受到的惯性力矩M的作用，这个力矩称为陀螺力矩： k ) 由matlab计算得，M二00Nm k 3、风轮载荷计算、风轮载荷计算 轴向诱导因子，周向诱导因子 作用在整个风轮上的轴向推力 作用在整个风轮上的转矩由matlab计算得，T=718370N，M=4613600Nm。 4、 ( ( 主轴载荷计算主轴载荷计算 5 5、 低速轴角速度为：s 髙速轴角速度为：s 低速轴功率为： 髙速轴功率为： 低速轴转

19、矩为：m 髙速轴转矩为：m 假设所选择的实心钢轴的推荐最大应力为55MPa，则： 低速轴的直径： 髙速轴的直径： 6、塔架载荷计算、塔架载荷计算 除了风载荷外，风电机组的几乎所有的载荷也都将传给塔架，按风电机组的载荷源分类又要有：风轮等构件承受的空气动力载荷 重力和惯性载荷：由重力、振荡、旋转以及地震等引起的静态和动态载荷 操作载荷：在机组运行和控制过程中产生的载荷。如功率变化、偏航、变浆以及制动过程产生的载荷等 其他载荷：诸如尾迹载荷、冲击载荷、覆冰载荷等、 下面只讨论与塔架结构强度计算有关的两种载荷：即由风轮作用的最大气动推力以及塔架本身所承受最大风压产生的载荷。 1)1) 暴风工况的风轮

20、气动推力计算暴风工况的风轮气动推力计算 暴风工况所对应的最大风速v与当地的年平均风速v有关。由于当地的年平均风速为s,所以 save v取60m/s。 s 前苏联的法捷耶夫公式： F= as 式中：A叶片的投影面积，单位为m2 S V风轮中心处的暴风风速，单位为m/s s B风轮叶片数 叶片的投影面积 其中。为风轮的实度，风轮实度与叶尖速比有关，当入=5时，可以近似的认为。=由matlab计算得，F二1108400N。 as 荷兰ECN的公式 式中：C推力系数，取C= tt q动态风压，单位是N/m2 动态系数，取二 S安全系数，取s=其中的q是随着髙度变化，这里q的值为1390N/m2 由m

21、atlab计算得，F二1473800N as 德国DFVLR公式 式中：C推力系数，取C= tt v风轮中心处的暴风风速，取v=60m/s ss 由matlab计算得，F二1905000N as ！ 丹麦RIS公式 式中：P风轮单位扫掠面积上的平均风压，通常取P=300N/m2 A风轮的扫掠面积，单位m2 s 由matlab计算得，F二2356200N as 2）欧美国家塔架静态强度设计的一般载荷条件欧美国家塔架静态强度设计的一般载荷条件 风载条件：风速65m/s，（2s,平均）风轮停转，叶片顺浆，风向沿机舱横向作用在塔架上。 正常运行工况条件的地震载荷：考虑额定风速时产生的风轮最大轴向力，同

22、时根据均匀建筑物由地震产生的水平载荷因子，将其产生的惯性力附加在风轮轴向推力上。 【 机组的最大运行载荷：一般为额定风速下正常运行载荷的2倍。 3 3）确定塔架设计载荷的要求）确定塔架设计载荷的要求 设计载荷需要正确的分析塔架承受的载荷及其作用，大致可分为三种类型设计载荷： 最大极限载荷：指塔架可能承受的最大载荷。 疲劳载荷：指塔架结构能够承受交变载荷次数的能力。 共振激励载荷：塔架结构系统的共振响应。 7、联轴器载荷计算、联轴器载荷计算 联轴器所承受的载荷以转矩为主，所以可以利用公式来计算。 低速轴联轴器载荷：=2250630Nm 髙速轴联轴器载荷：18454Nm 阶段性总结阶段性总结 气动

23、特性 攻角 5 升力系数 】 阻力系数 发电机 额定功率 3402kW 1650r/min 额定转速 齿轮箱 齿轮箱传动比 120 低速轴转速 min 髙速轴转速 1650r/min 齿轮箱功率 3244kW 联轴器 低速轴功率 3400kW 低速轴载荷 2250630Nm 髙速轴功率 3500kW 髙速轴载荷 18454Nm 变流器 变流器功率 1500kW 偏航电机 偏航电机功率 4kW & 叶片载荷 作用在叶片上的离心力 00N 作用在叶片上的风压力 755580N 作用在叶片上的气动力矩 & 144250Nm 作用在叶片上的陀螺力矩 00Nm 风轮载荷 作用在整个风轮上

24、的轴向推力 718370N 】 作用在整个风轮上的转矩 4613600Nm 主轴 低速轴转矩 m 低速轴直径 m 髙速轴转矩 髙速轴直径 塔架载荷 作用在塔架上的载荷 前苏联法捷耶夫公式 1108400N 荷兰ECN的公式 1473800N 德国DFVLR公式 1905000N 丹麦RIS公式 2356200N 三、塔架根部截面应力计算三、塔架根部截面应力计算 塔架通常采用等强度变截面的设计，危险截面一般位于塔架的根部。塔筒根部的结构强度分析是确定塔架整体结构尺寸的基本设计依据。 下面考虑塔架高度折减系数的强度计算： 塔架根部的截面应力可以表示为： 其中：W2塔架根部抗弯截面模数，单位是cm2

25、 A2塔根部截面积，单位cm2 G2塔架本身所受重力，单位为N 变截面塔架的长度折减系数，这里取 式中：卩与塔架截面变化有关的折算长度修正系数，可以根据之比确定 J塔架顶部截面惯性矩，单位cm4 min J塔架根部截面惯性矩，单位cm4 max H塔架髙度，单位为cm Y2塔架根部截面的惯性半径，单位是cm F取1905000Nas h=1/20*R= H=90m=9000cm 塔架上部直径取D=40m，下部直径取D=50m，壁厚取45/20mm 其中D=50m，a=()/=代入数值，计算的W2= A= 2 G+G2为塔顶质量与塔架本身重量之和，取400tJ= min J= max 因此，=，

26、所以卩取值为 即， 塔架的根部截面应力为：40Mpa 四、课程设计结论四、课程设计结论 所设计的风电机组的技术参数见下表： 机组 风力机等级 IECIIIA 额定功率 3MW 切入风速 ? 3m/s 切出风速 25m/s 额定风速 13m/s 设计寿命 20年 设计标准 IEC61400 运行温度 -1545 ! 叶轮 叶轮直径 92m 叶片长度 叶片数量 3 额定转速 min 齿轮箱 结构形式 2级行星+1级平行轴 传动比 120 发电机 型式 双馈异步感应电机 额定输出电压 690V 额定功率 3402kW 频率 50Hz 功率因数 容性感性 额定转速/范围 1650/(10002000)

27、rpm 变浆系统 驱动控制 电机 偏航系统 型式 主动式对风式 控制方式 电机+减速机 控制系统 控制方式 远程监控 塔筒 型式 钢制锥形塔筒 塔筒髙 90m 五、参考文献五、参考文献 【1】姚兴佳、田德，宋俊、莴晓明，风力发电机组设计与制造，华北电力大学，校内试用教材 【2】徐大平、柳亦兵、吕跃刚，风力发电原理，机械工业出版社 【3】朱永强、张旭，风电场电气系统，机械工业出版社 【4】敬照亮，MATLAB教程与应用清华大学出版社北京交通大学出版社 六、附录六、附录 叶片设计和各部分载荷计算的MATLAB程序： D二ceil(sqrt(8*1000000/pi/(133)； n=14; lamda=5