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

课程设计(综合实验)报告(2012–2013年度第二学期)名称：题目：院系：班级：学号：学生姓名：指导教师：设计周数：成绩：日期：2013年月日１《风力发电机组设计与制造》课程设计任务书一、设计内容风电机组总体技术设计二、目的与任务主要目的：1.以大型水平轴风力机为研究对象，掌握系统的总体设计方法；2.熟悉相关的工程设计软件；3.掌握科研报告的撰写方法。主要任务：1.确定风电机组的总体技术参数；2.关键零部件（齿轮箱、发电机和变流器）技术参数；3.计算关键零部件（叶片、风轮、主轴、连轴器和塔架等）载荷和技术参数；4.完成叶片设计任务；5.确定塔架的#设计#。6.撰写一份课程设计报告。三、主要内容选择功率范围在1.5MW至6MW之间的风电机组进行设计。1）原始参数：风力机的安装场地50米高度年平均风速为7.0m/s，60米高度年平均风速为7.3m/s，70米高度年平均风速为7.6m/s，当地历史最大风速为49m/s，用户希望安装1.5MW至6MW之间的风力机。采用63418翼型，63418翼型的升力系数、阻力系数数据如表1所示。空气密度设定为1.225kg/m3。2）设计内容（）确定整机设计的技术参数。设定几种风力机的曲线和曲线，风力机基本参数包括叶1CpCt片数、风轮直径、额定风速、切入风速、切出风速、功率控制方式、传动系统、电气系统、制动系统形式和塔架高度等，根据确定风力机等级；（）关键部件气动载荷的计算。设定几种风轮的曲线和曲线，计算几种关键零部件的载2CpCt荷（叶片载荷、风轮载荷、主轴载荷、连轴器载荷和塔架载荷等）；根据载荷和功率确定所选定机型主要部件的技术参数（齿轮箱、发电机、变流器、连轴器、偏航和变桨距电机等）和型式。２以上内容建议用计算机编程实现，确定整机和各部件（系统）的主要技术参数。（3）塔架根部截面应力计算。计算暴风工况下风轮的气动推力，参考风电机组的整体设计参数，计算塔架根部截面的应力。最后提交有关的分析计算报告。四、进度计划序号设计(实验)内容完成时间备注1风电机组整体参数设计2.5天2风电机组气动特性初步计算2天3机组及部件载荷计算2天4齿轮箱、发电机、变流器技术参数1.5天4塔架根部截面应力计算1天5报告撰写1.5天6课程设计答辩1.5天五、设计（实验）成果要求1．提供设计的风电机组的性能计算结果；2．绘制整机总体布局工程图。六、考核方式提交一份课程设计报告；准备课程设计PPT，答辩。风电机组总体技术设计一、确定风电机组的总体技术参数３1、叶片数B=3风轮的叶片数取决于风轮尖速比，一般来说，要得到很大的输出扭矩就需要较大的叶片实度；现代风力发电机组实度较少，一般只需要1~3个叶片。叶片数多的风力机在低尖速比运行时有较高的风能利用系数，既有较大的转矩，而且启动风速低，适用于提水。而叶片数少的风力机则在高尖速比运行时有较高的风能利用系数，但启动风速高，因此适用于风力发电。3叶片风轮有以下优点：平衡简单，动载荷小；3叶片使风力发电机组系统运行平稳，基本消除了系统的周期载荷，输出较稳定的转矩；3叶片风轮通常能提供较佳的效率；3叶片风轮从审美的角度来说更令人满意；3叶片风轮的受力平衡好，轮毂可以简单些；噪声相对较小。综上所述，叶片数B选取3。2、机组额定功率及种类选取3MW级别的双馈变速恒频式风力发电机组3、切入风速取切入风速Vin=3m/s。4、切出风速取切出风速Vout=25m/s。5、额定风速取额定风速Vr=13m/s。6、各部分效率传动系统效率η1取0.95发电机效率η2取0.96变流器效率η取30.967、风能利用系数Cp额定功率下风能利用系数Cp取0.388、风轮直径由风力发电机组输出功率得叶片直径：8PD=r==92mV3Crp123其中：风力发电机组额定输出功率，取；Pr--3000kW--空气密度（一般取标准大气状态），取1.225kg/m3；额定风速，取；Vr--13m/sD--风轮直径；η1--传动系统效率，取0.95；η2--发电机效率，取0.96；η变流器效率，取；3--0.96额定功率下风能利用系数取。Cp--,0.389、功率Pw曲线，风能利用系数CP和推力系数CT曲线的确定风电机组的设计需要给出功率输出特性曲线。如果确定了风电机组的切入风速，额定风速和切出风速，可以利用下面的方程计算输出功率，查阅相关资料，得到其风能利用系数随着风速的变化值，４得到下表。（其中所用为：Cp=4a(1-a)2Ct=4a(1-a)8PD=rV3Crp123)利用上述的公式，由excel的函数运算，计算得到随着风速的变化，Cp、和Ct的值得变化风速风能利用系数轴向诱导因功率Pw[kW]推力系数CT叶尖速比v[m/s]Cp子a30.000000.000000.000000.0000022.479840764131.250000.476930.174940.5773316.859880575300.000000.558150.243340.7365113.487904466506.250000.545060.229120.7065111.239920387750.000000.508520.197270.633419.63421746981031.250000.468420.169710.563638.42994028691350.000000.430670.148440.505637.493280254101706.250000.396810.131510.456856.743952228112100.000000.366930.117510.414806.130865662122531.250000.340670.106520.380685.61996019133000.000000.317560.097330.351425.18765556143000.000000.254260.074130.274534.817108735153000.000000.206720.058270.219494.495968152163000.000000.170330.046610.177734.214970143173000.000000.142010.038260.147203.967030723183000.000000.119630.031850.123333.746640127193000.000000.101720.026820.104413.549448541203000.000000.087210.022820.089193.371976114213000.000000.075340.019590.076823.211405823223000.000000.065520.016710.065733.065432831233000.000000.057340.014630.057652.932153143243000.000000.050470.012870.050832.809980095253000.000000.044650.011390.045042.697580891５10、风轮转速风力发电机组产生的气动噪声正比于叶尖速度的5次方。通常将路基风力发电机组的叶尖速度限制在65m/s左右。这里取叶尖速V=65m/s。由v=3.14nD/60得到n=60v/3.14D，将V=65m/s带入得，13.5r/min。11、叶尖速比λ由公式，带入额定转速n=13.5r/min可以得到参考CP-λ曲线，当Cp=0.38时对应的λ的值约为5，因此所选数值合理。12、功率控制方式功率控制方式选取当今应用较为广泛的变浆距控制。这种机组当风速过高时，通过减小叶片翼型上合成气流与翼型几何弦的夹角（攻角），改变风力发电机组获得的空气动力转矩，能使功率输出保持稳定。同时，风机在起动过程也需要通过变浆距来获得足够的起动转矩。采用变浆距技术的风力发电机组还可以使叶片和整机的受力状况大为改善，这对于大型风力发电机组十分的有利。６13、传动系统传动系统选定由以下四部分组成：主轴系统、增速传动系统（齿轮箱）、轴系的支撑与连接（如轴承和联轴器）、制动装置。14、制动系统制动系统的形式以空气动力制动和机械制动为主。15、电气系统风电机组的电气系统包括以下四个主要部分：变浆距驱动，偏航驱动，发电机部分，控制和监控部分16、塔架高度由于风速与距地面高度有关，增加塔架高度可使风轮获取更多的风能，但制造更高的塔架也需要更多的材料，使其造价相应增加。大型机组的塔架高度H可以按下式初步确定：式中D为风轮直径，单位为m。塔架的高度参数的选择与地形和地貌有关，如陆地与海上风电机组会有所不同。陆地的地表性对粗糙且风速随高度变化缓慢，可能需要较高的塔架；而海平面较光滑，风速沿着高度方向变化梯度大。一般而论，即使相同容量的机组，在近海安装的机组塔架高度相对较低。考虑地貌因素的塔架最低设计高度一般可用下式估算：H=h+C+R式中：h——机组附近的障碍物高度，单位m；C——障碍物最高点到风轮扫掠面最低点的距离（最小取1.5-2.0m），单位m。综上考虑，塔架高度H取90m17、设计寿命按照惯例，设计寿命取20年。18、设计标准设计标准为IEC61400-119、风力机等级70m年平均风速为7.6m/s。去IECⅢA，参考下表：风电机组等级第Ⅰ级第Ⅱ级第Ⅲ级第Ⅳ级参考风速Vr=5Va(m/s)5042.537.530年平均风速Va(m/s)108.57.56ε150.180.180.180.18A湍流强度a2222等级ε150.160.160.160.16Bb3333７阶段性总结叶片数3叶尖速比5切入风速3m/s功率控制变浆距控制切出风速25m/s制动系统形式空气动力制动和机械制动额定风速13m/s塔架高度90m额定功率3MW设计寿命20年风能利用系数0.32设计标准IEC61400风轮转速13.5r/min风力机等级IECⅢA风轮直径92m二、关键部件气动载荷的计算包括设定几种风轮的Cp曲线和Ct曲线，计算几种关键零部件的载荷（叶片载荷、风轮载荷、主轴载荷、连轴器载荷和塔架载荷等）；根据载荷和功率确定所选定机型主要部件的技术参数（齿轮箱、发电机、变流器、连轴器、偏航和变桨距电机等）和型式。以上内容建议用计算机编程实现，确定整机和各部件（系统）的主要技术参数。最后提交有关的分析报告。根据所提供的63418翼型的升力系数、阻力系数数据，得到攻角与升力系数、阻力系数的关系曲线。由图可以观察出，攻角约等于13°时升力系数最大，大约为1.400左右。攻角约等于9°时，升力系数约为最大升力系数的0.9倍，为1.244，此时的阻力系数为0.015。当攻角为5°时升阻比最大，此时升力系数为0.904。在叶片处选择10段截面，分别在10%R，20%R，30%R……90%R，100%R。由以下公式：16R/rN9r42()200R9８、给出，攻角已知，已知；ariCL=0.9041b、由式arctg求得角；33c、由式k(21)cos得到k；1k2d、由式h1得到h；21he、由式Iarcctg()得到I；1kf、Ii，C[8r(h1)cosI]/[BC(h1)]；Lg、将a--g编程进行运算得到数据--表2。（程序见附录）通过excel计算得到以上各个参数的数值：位置半径rλΨkhφβC10%50.51.201750.403312.085510.738100.5811012.0590420%1011.309000.366031.366030.523600.3666012.4157430%151.51.374800.351091.178840.392000.2350010.5443340%2021.416250.344211.104710.309100.152108.7838150%252.51.443960.340601.068380.253670.096677.4143260%3031.463550.338511.048040.214500.057506.3713570%353.51.478030.337191.035540.185530.028535.5665380%4041.489140.336321.027340.163320.006324.9327490%454.51.497910.335711.021670.14578-0.011224.42339100%5051.505000.335271.017600.13160-0.025404.006411、主要部件功率计算1)发电机发电机类型：双馈异步变速恒频式发电机额定功率：P=P/η*1.1=3000/0.97*1.1=3402kW（预留10%的裕度）发电机r3极对数：2由公式60f=np可以得到发电机的额定转速为1650m/min（预留10%的裕度）2)齿轮箱齿轮箱类型：增速，2级行星齿轮+1级平行齿轮齿轮箱效率0.95传动比：120低速轴转速：13.5r/min；高速轴转速：1650r/min齿轮箱功率：PGB=3581kw９3)联轴器低速轴联轴器动率：3400kW高速轴联轴器动率：3500kW4)变流器变流器的功率通常为风电机组的额定功率的1/2~1/3，考虑到风电机组的可靠性，通常为风电机组额定功率的1/2，即变流器的功率为1500kW。5)偏航电机偏航电机的功率由最大偏航扭矩确定。在此先初步选定4台4kW的偏航电机。2、主要部件的载荷计算1)叶片的载荷计算叶片绕风轮轴旋转时，有离心力作用在叶片上。离心力的方向是自旋转中心沿半径向外。作用在叶片上的离心力可按照下式求解：Fcρ——叶片的密度，单位是3，取3ykg/m1.8kg/m——叶素处的叶片截面面积，单位是2ArmΩ——风轮的角速度，单位是rad/s——叶片起始处的旋转半径r0R——叶片结束处的旋转半径其中：Ω=rad/s由计算的：。matlabFc=1026500000N2)作用在叶片上的风压力风压力是作用在叶片上沿风速方向的气动力。风轮转动时的风压力：设的作用点几粒风轮轴的距离为，则：FVrm由matlab计算得：Fv=755580N，rm=36.9690m。作用在叶片上的气动力矩3)Mb是使风轮转动的力矩，由计算得·。Mbmatlab:Mb=144250Nm4)作用在叶片上的陀螺力矩是风轮对风调向时产生的惯性力矩。当风向改变时，风轮除以角速度Ω绕轴转动外，还要MkOX以角速度ω绕OZ轴转动。整个叶片的转动惯量为：由matlab计算得，J=607400000kg·m2对风调向时，叶片可看成以角速度ω绕OZ轴转动的牵连运动、以Ω绕OX轴转动的相对运动的１０复合运动，这时要产生柯式加速度：ak=2w*v这里是惯性半径，β是与ω的夹角。当叶片旋转到铅垂位置时，β°，柯式加速度达到最Vr=0大值：由于柯式加速度，就产生了柯式角加速度由动量矩定理知，叶片受到的惯性力矩的作用，这个力矩称为陀螺力矩：Mk由计算得，·matlabMk=2053000000Nm3、风轮载荷计算轴向诱导因子，周向诱导因子作用在整个风轮上的轴向推力作用在整个风轮上的转矩由matlab计算得，T=718370N，M=4613600N·m。4、主轴载荷计算低速轴角速度为：1.466rad/s高速轴角速度为：172.788rad/s低速轴功率为：3299.353kW高速轴功率为：3188.437kW低速轴转矩为：2250581.865N·m高速轴转矩为：18452.884M·m假设所选择的实心钢轴的推荐最大应力为55MPa，则：低速轴的直径：0.5929m高速轴的直径：0.1196m5、塔架载荷计算除了风载荷外，风电机组的几乎所有的载荷也都将传给塔架，按风电机组的载荷源分类又要有：风轮等构件承受的空气动力载荷重力和惯性载荷：由重力、振荡、旋转以及地震等引起的静态和动态载荷操作载荷：在机组运行和控制过程中产生的载荷。如功率变化、偏航、变浆以及制动过程产生的载荷等其他载荷：诸如尾迹载荷、冲击载荷、覆冰载荷等、下面只讨论与塔架结构强度计算有关的两种载荷：即由风轮作用的最大气动推力以及塔架本身所承受最大风压产生的载荷。1）暴风工况的风轮气动推力计算暴风工况所对应的最大风速与当地的年平均风速有关。由于当地的年平均风速为，vsvave7.6m/s１１所以取。vs60m/s前苏联的法捷耶夫公式：2Fas=0.784AbvsB式中：——叶片的投影面积，单位为2Asm——风轮中心处的暴风风速，单位为Vsm/sB——风轮叶片数叶片的投影面积其中σ为风轮的实度，风轮实度与叶尖速比有关，当λ=5时，可以近似的认为σ=0.05由计算得，。matlabFas=1108400N荷兰ECN的公式式中：——推力系数，取CtCt=1.5q——动态风压，单位是N/m2ψ——动态系数，取ψ=1.2S——安全系数，取S=1.5其中的q是随着高度变化，这里q的值为1390N/m2由计算得，matlabFas=1473800N德国DFVLR公式式中：——推力系数，取CtCt=2.2——风轮中心处的暴风风速，取vsvs=60m/s由计算得，matlabFas=1905000N丹麦RIS公式式中：——风轮单位扫掠面积上的平均风压，通常取2PlPl=300N/m——风轮的扫掠面积，单位2Asm由计算得，matlabFas=2356200N2)欧美国家塔架静态强度设计的一般载荷条件➢风载条件：风速65m/s，（2s，平均）风轮停转，叶片顺浆，风向沿机舱横向作用在塔架上。➢正常运行工况条件的地震载荷：考虑额定风速时产生的风轮最大轴向力，同时根据均匀建筑物由地震产生的水平载荷因子，将其产生的惯性力附加在风轮轴向推力上。➢机组的最大运行载荷：一般为额定风速下正常运行载荷的2倍。3)确定塔架设计载荷的要求设计载荷需要正确的分析塔架承受的载荷及其作用，大致可分为三种类型设计载荷：最大极限载荷：指塔架可能承受的最大载荷。疲劳载荷：指塔架结构能够承受交变载荷次数的能力。共振激励载荷：塔架结构系统的共振响应。１２6、联轴器载荷计算联轴器所承受的载荷以转矩为主，所以可以利用公式来计算。低速轴联轴器载荷：=2250630N·m高速轴联轴器载荷：18454N·m阶段性总结攻角5°气动特性升力系数0.904阻力系数0.007额定功率3402kW发电机额定转速1650r/min齿轮箱传动比120低速轴转速13.5r/min齿轮箱高速轴转速1650r/min齿轮箱功率3244kW低速轴功率3400kW低速轴载荷2250630N·m联轴器高速轴功率3500kW高速轴载荷18454N·m变流器变流器功率1500kW偏航电机偏航电机功率4kW作用在叶片上的离心力1026500000N作用在叶片上的风压力755580N叶片载荷作用在叶片上的气动力矩144250N·m作用在叶片上的陀螺力矩2053000000N·m作用在整个风轮上的轴向推力718370N风轮载荷作用在整个风轮上的转矩4613600N·m低速轴转矩2250581.865N·m低速轴直径18452.884M·m主轴高速轴转矩0.5929m高速轴直径0.1196m前苏联法捷耶夫公式1108400N作用在荷兰ECN的公式1473800N塔架载荷塔架上德国DFVLR公式1905000N的载荷丹麦RIS公式2356200N三、塔架根部截面应力计算１３塔架通常采用等强度变截面的设计，危险截面一般位于塔架的根部。塔筒根部的结构强度分析是确定塔架整体结构尺寸的基本设计依据。下面考虑塔架高度折减系数的强度计算：塔架根部的截面应力可以表示为：其中：——塔架根部抗弯截面模数，单位是2W2cm——塔根部截面积，单位2A2cm——塔架本身所受重力，单位为G2Nψ——变截面塔架的长度折减系数，这里取20.60式中：μ——与塔架截面变化有关的折算长度修正系数，可以根据之比确定——塔架顶部截面惯性矩，单位4Jmincm——塔架根部截面惯性矩，单位4JmaxcmH——塔架高度，单位为cmγ——塔架根部截面的惯性半径，单位是2cm取Fas1905000Nh1=1/20*R=2.3mH=90m=9000cm塔架上部直径取，下部直径取，壁厚取D1=40mD2=50m45/20mm其中D=50m，α=（2.3-0.045）/2.3=0.98代入数值，计算的3W2=850.01m2A2=65.05m为塔顶质量与塔架本身重量之和，取G1+G2400t4Jmin=8608.35m4Jmax=21420.12m因此，=0.4096，所以μ取值为1.24即，塔架的根部截面应力为：40Mpa四、课程设计结论所设计的风电机组的技术参数见下表：风力机等级IECⅢA额定功率3MW机组切入风速3m/s切出风速25m/s１４额定风速13m/s设计寿命20年设计标准IEC61400运行温度-15°~45°叶轮直径92m叶片长度47.5m叶轮叶片数量3额定转速13.5r/min结构形式2级行星+1级平行轴齿轮箱传动比120型式双馈异步感应电机额定输出电压690V额定功率3402kW发电机频率50Hz功率因数容性0.95~感性0.9额定转速/范围1650/(1000~2000)rpm变浆系统驱动控制电机型式主动式对风式偏航系统控制方式电机+减速机控制系统控制方式远程监控型式钢制锥形塔筒塔筒塔筒高90m五、参考文献【1】姚兴佳、田德，宋俊、莴晓明，《风力发电机组设计与制造》，华北电力大学，校内试用教材【2】徐大平、柳亦兵、吕跃刚，《风力发电原理》，机械工业出版社【3】朱永强、张旭，《风电场电气系统》，机械工业出版社【4】敬照亮，《MATLAB教程与应用》，清华大学出版社·北京交通大学出版社六、附录叶片设计和各部分载荷计算的MATLAB程序：>>D=ceil(sqrt(8*3.0*1000000/0.95/0.97/0.97/0.32/pi/1.225/(13^3)));>>n=14;>>lamda=5;>>a=[0.1:0.001:1];>>b=lamda*a;>>psi=1/3*atan(b)+pi/3;１５>>k=sqrt(b.^2+1).*cos(psi);>>h=sqrt((1-k.^2)./b.^2+1);>>I=atan((1+k)./(1+h)./b);>>theta=I-5*pi/180;>>C=8*pi*a.*(D/2).*(h-1).*cos(I)/3/0.904/(h+1);>>R=D/2;>>r=a*R;>>omega=lamda*13/R;>>yl=r.^2.*C;>>Fc=1.8*1000*omega^2*trapz(r,yl)>>cl=0.904;>>cd=0.007;>>y2=(1+(1./tan(I)).^2).*(cl*cos(I)+cd*sin(I)).*C;>>Fv=1/2*1.225*13^2*trapz(r,y2);>>y3=y2.*r;>>rm=trapz(r,y3)./trapz(r,y2)>>y21=(1+(1./tan(I)).^2).*(cl*sin(I)-cd*cos(I)).*C;>>Mb=1/2*1.225*13^2*trapz(r,y21);>>J=Fc/omega^2;>>epsilon=2*omega*lamda*13/R;>>Mk=J*epsilon>>r1=[0,r];>>a1=[0,a];>>b1=lamda*a1;>>psil=1/3*atan(b1)+pi/3;>>k1=sqrt(b1.^2+1).*cos(psil);>>y4=(1-k1.^2).*r1;>>T=pi*1.225*13^2*trapz(r1,y4);>>h1=sqrt((1-k1.^2)./b1.^2+1);>>y5=(h-1).*(k+1).*r.^3;>>M=pi*1.225*omega*13*trapz(r,y5)>>Fas1=pi*D^2/4*0.05/3*0.784*60^2*3>>Fas2=1.5*1390*1.2*0.05*pi*D^2/4*1.5>>Fas3=2.2*0.5*1.225*60^2*0.05*pi*D^2/4>>Fas4=300*pi*D^2/4七、相关图片1.偏航系统结构平面图１６2.机舱结构图１７