风力发电机设计与制造课程设计

1、一总体参数设计总体参数是设计风力发电机组总体结构和功能的基本参数，主要包括额定功率、发电机额定转速、风轮转速、设计寿命等。1. 额定功率、设计寿命根据设计任务书选定额定功率Pr =3.5MW；一般风力机组设计寿命至少为20 年，这里选 20 年设计寿命。2. 切出风速、切入风速、额定风速切入风速 取 V in = 3m/s切出风速取 V out = 25m/s额定风速 V r = 12m/s （对于一般变桨距风力发电机组（选3.5MW）的额定风速与平均风速之比为 1.70 左右 ,V r =1.70Vave=1.70 ×7.0 12m/s）3. 重要几何尺寸(1) 风轮直径和扫掠面积

2、由风力发电机组输出功率得叶片直径：D8Pr8 350000031.225 123104mVrC p 1 2 30.45 0.95 0.96 0.95其中：Pr风力发电机组额定输出功率，取3.5MW ；错误 !未找到引用源。 空气密度（一般取标准大气状态），取1.225kg/m3；V r额定风速，取12m/s；D风轮直径；1 传动系统效率，取0.95；2 发电机效率，取0.96；错误 !未找到引用源。3 变流器效率，取0.95；Cp额定功率下风能利用系数,取 0.45。由直径计算可得扫掠面积：AD 210428482m2错误 !未找到引用源。错误 ! 未找到引用源。44综上可得风轮直径D=104

3、m，扫掠面积 A=8482 m2精选文库4. 功率曲线自然界风速的变化是随机的, 符合马尔可夫过程的特征, 下一时刻的风速和上一时刻的结果没什么可预测的规律。由于风速的这种特性, 可以把风力发电机组的功率随风速的变化用如下的模型来表示：P ( t )P sta t ( t )P ( t )P(t ) 在真实湍流风作用下每一时刻产生的功率, 它由 t 时刻的 V(t)决定 ;Pstat (t ) 在给定时间段内V(t)的平均值所对应的功率；P (t) 表示 t 时刻由于风湍流引起的功率波动。对功率曲线的绘制 , 主要在于对风速模型的处理。 若假定上式表示的风模型中Pstat(t)的始终为零 ,

4、即视风速为不随时间变化的稳定值 , 在切入风速到切出风速的范围内按照设定的风速步长, 得到对应风速下的最佳叶尖速比和功率系数 ,带入式：P132C P V r D12381 传动系统效率，取0.95；2 发电机效率，取0.96；错误 !未找到引用源。3 变流器效率，取0.95；错误 !未找到引用源。 空气密度（一般取标准大气状态），取1.225kg/m3；V r额定风速，取12m/s；D风轮直径；Cp额定功率下风能利用系数,取 0.45。2精选文库由以上公式，使用excel 计算出不同风速对应的功率值,见表 1表 1 风速功率关系风速（ m/s）34567891011功率（ w）5474412

5、97632534444379526954521038109147809020275582698680风速（ m/s）121314151617181920功率（ w）350000035000003500000350000035000003500000350000035000003500000风速（ m/s）2122232425功率（ w）35000003500000350000035000003500000将得到的数据对绘制成静态风功率曲线，如图一图 1Pv 静态功率曲线3精选文库5. 风轮额定转速三叶片风力发电机组的风轮叶尖速比0 一般在 6 至 8 之间，不同攻角下的风能利用系数随叶尖速比的

6、变化曲线即CP0 曲线如图 。图C P0曲线由 Cp0 曲线可得出0 =7.5，则风轮额定转速可由下式计算得到：Rnr30Vrnr600Vr60 7.5 12 16.5r / minD1046. 叶片数现代风力发电机的实度比较小，一般需要1-3 个叶片。选择风轮叶片数时考虑风电机组性能和载荷、风轮和传动成本、风力机气动噪声及景观影响等因素。3 叶片较 1、2 叶片风轮有如下优点：平衡简单、动态载荷小。基本消除了系统的周期载荷，输出较稳定转矩；能提供较佳的效率；更加美观；噪声较小；4精选文库轮毂较简单等。综上所述，叶片数选择3。7、功率控制方式、制动系统形式功率控制方式选择主动变桨距控制；制动系

7、统形式为第一制动采用气动刹车，第二制动采用高速轴机械刹车。8 、风力机等级由 IEC 标准，如 表 2，选择风力机等级为 IECIIIA 。WTGS 等级IIIIIISvref (m / s)5042.537.5AI ref0.16设计值由设计者选定BI ref0.14CI ref0.12表 2风机等级规范表注：表中数据为轮毂高度处值，其中：A 表示较高湍流特性级；参考风速Vref 为 10min 平均风速； B 表示中等湍流特性级；I 15 风速为 15m/s 时的湍流强度特性值。 C 表示较低湍流特性级；除表基本参数外，在风力发电机组设计中，还需要某些更重要的参数来规定外部条件。对风力发电

8、机组 IA IIIC 级，统称为风力发电机组的标准等级。阶段性总结表总体参数设计值总体参数设计值叶片数B=3风轮直径D=104m额定输出功率P=3.5MW轮毂高度Zhub=106.25m设计寿命20 年风能利用系数Cp=0.45切入风速V in=3m/s叶尖速比7.5切出风速V out功率控制方式主动变桨距控制=25m/s额定风速V r=12m/s制动形式气动刹车、机械刹车风轮额定转速nr=16.5r/min传动系统高传动比齿轮箱传动风力机等级IECIIIA电气系统双馈发电机 +变流器5精选文库二叶片设计1. 叶片材料选择叶片选用 T-700 碳纤维，相比玻璃纤维，叶片密度较小，发电效率更高，

9、密度为 1800 kg / m3 。2. 计算各剖面的叶尖速比将叶片分为 20 个叶素，每个叶素间隔 0.05R，其中 5%半径处叶片是筒状， 10%-60%半径处采用钝后缘叶片， 65%-100%半径处 采用通用风电机组叶片翼型。叶片内圈采用钝后缘翼型，外圈采用 63415 翼型。 根据下式求各叶素的叶尖速比 。r0 R叶素位置和叶尖速比数值见下表 2：表 2不同叶素位置的叶尖速比叶素位置 /%5101520253035404550叶尖速比0.3750.7501.1251.5001.8752.2502.6253.0003.3753.750叶素位置 /%55606570758085909510

10、0叶尖速比4.1254.5004.8755.2505.6256.0006.3756.7507.12575003. 根据翼型确定叶片最佳攻角，升力系数 Cl，Cd风力机翼型为 NACA63-415 ，图 3图 3 NACA63-415 翼型图计算雷诺数 Re在 20 ，压强为标准大气压101.325kPa 时，空气的动力粘度17.910 66VL iR e精选文库1.22512526.28 10717 .9106根据所得雷诺数查得Cl/Cd 、Cl/alpha ，见图 4图 4Cl-C d 曲线和 Cl-alpha 曲线图 5 C l/C d alpha 图从图中可以得出翼型取得最大升阻比时，最

11、佳攻角5.25，此时升力系数 Cl=0.9461，Cd=0.00791，最佳升阻比 C l / Cd118.7,本次设计选取最佳攻角5.25 错误 ! 未找到引用源。，则升力系数和阻力系数分别为Cl， d。叶片每个截面的升力系数相同，为=0.9461 C =0.00791Cl=0.9461。4. 叶片弦长计算步骤7精选文库通过下面的计算，可以得到沿叶片各径向位置r 上的弦长 C 和叶素桨距角，即可完成叶片的初步设计，但要是想对叶片进一步优化，还需对翼型、叶根、叶尖风进行气动优化设计和工艺优化设计，在本次设计报告中，只对叶片作了初步设计。如下：（1）求利用公式1 arctan33（2）求轴向干扰

12、因子 k 利用公式k21 cos（3）求切向干扰因子 h 利用公式12hk12(4）求入流角利用公式1harc cot()（5）求叶素桨距角（6）计算叶片弦长C8 r (h1) cosCBCl (h1)8精选文库叶片气动特性通过excel 计算，得到叶片各个截面气动特性参数，如表 3：表 3叶片气动特性参数位置 (%)半径 r(m)叶尖速比khC(修正)52.60.3751.1670.4202.6190.8080.7167.116105.20.751.2620.3801.5880.6180.5278.522157.81.1251.3290.3611.2990.4840.3937.9472010

13、.41.51.3750.3511.1790.3920.3006.98525131.8751.4070.3461.1180.3270.2356.0863015.62.251.4310.3421.0840.2790.1875.3353518.22.6251.4490.3401.0620.2430.1514.7214020.831.4640.3391.0480.2150.1234.2214523.43.3751.4750.3371.0380.1920.1003.80950263.751.4840.3371.0310.1740.0823.4665528.64.1251.4920.3361.0260.1

14、590.0673.1776031.24.51.4980.3361.0220.1460.0542.9306533.84.8751.5030.3351.0190.1350.0432.7187036.45.251.5080.3351.0160.1250.0342.53475395.6251.5120.3351.0140.1170.0262.3738041.661.5160.3351.0120.1100.0182.2308544.26.3751.5190.3351.0110.1040.0122.1049046.86.751.5220.3341.0100.0980.0061.9919549.47.125

15、1.5240.3341.0090.0930.0011.889100527.51.5270.3341.0080.088-0.0031.7975. 叶片根部载荷计算与材料选择叶片根部处理方式：距叶根0 5m 处制作成直径为3m 的圆柱结构处理，且根部采用金属法兰连接。见图69精选文库图 6 金属法兰连接表 4 增强材料力学性能根据表 4 材料选择为 T700 碳纤维，抗拉强度为4.9Gpa取0.3 4.9 109 pa 1.47 109 pa32(T r m ) 2M 232( 66386034 .912 ) 230186002343941 .07 md min(1)1 .47 10(1 0.96

16、 )取 d1 c11 7.11 3.555m dmin22所以风轮根部直径选择3.6m10精选文库三 . 确定主要部件1发电机发电机类型：双馈异步变速恒频发电机；额定功率： 3.5MW ；额定转速： 1500r/min；发电机极对数为2，发电机主轴转矩 T发电机主轴 为：T发电机主轴 9550Pr955035002.443 104 N mn231500 0.960.95选择刚轴推荐最大扭剪应力：fs55MPa则发电机的主轴直径D 发电机为：D发电机22T发电机主轴232 2.4431040.1313m73f s5.510取发电机主轴直径D 为 0.15m.2变流器变流器功率通常为风力发电机组的

17、 1/21/3，为保证机组可靠性，通常为额定功率的 1/2, 所以变流器功率为 1500kW。3齿轮箱方式：行星齿轮传动两级NGW ；低速轴转速： nl=16.5r/min11精选文库高速轴转速： nh =1500r/min传动比： i = 90齿轮箱效率：330.95 0.9831齿轮箱功率：PGBPr35000006210.950.9633.9 10 W10.9534联轴器低速轴联轴器功率PmPr35000003.9 106W2223 310.960.953 0.95错误 !未找到引用源。高速轴联轴器功率：Pr35000006Pt0.963.8410 W230.955主轴低速轴角速度为：2

18、 n1 216.5m1.727rad / s6060高速轴角速度为：2 nh 21500t157rad / s6060低速轴功率为：PmPr35000003.9 106W320.960.953210.9523高速轴功率为：Pr3500000106WPt0.963.84230.9512精选文库低速轴转矩为：Pm3900000T m2258251 .3N mm1.727高速轴转矩为：Pt3840000T t24458.6 N mt157低速轴直径：D L22 T3fm232 2258251 .30 .59 ms55000000高速轴直径：D H22T t23 2 2.44510 40.1314 m

19、35.57fs10综上可得，低速轴直径取0.7m，高速轴直径取 0.18m。6偏航系统类型：主动偏航，并选用强制外置6 电机偏航；偏航范围：800o 800o偏航角速度： 0.6o / s偏航轴承： 4 点接触球轴承；偏航驱动： 6 个 3kW 偏航电机；偏航制动：液压控制摩擦制动；大齿轮齿数： 135；小齿轮齿数： 16减速箱传动比： i=140结构简图见图 713精选文库图 7 偏航机构结构简图7变桨系统根据调整桨距角调整风能利用系数，桨距角与风能利用系数曲线见图8图 8 风能利用系数与桨距角关系曲线变桨类型： 3 叶片独立变桨控制，采用电动驱动装置。供电方式：超级电容供电发电机转速：永磁

20、电动机1500r/min，电机功率： 310kw变桨范围： 0 90 o （主要变桨范围 030 o ）变桨速度： 7 o /s14精选文库传动方式：齿形带传动变桨机构结构图见图990 度限位开关0度接近开关变桨电机 1变桨柜 1滑环旋转编码器电磁刹车动力电源线连接器90 度限位开关0度接近开关变桨电机 2变桨柜 1旋转编码器电磁刹车动力电源线90 度限位开关0度接近开关变桨电机 3旋转编码器变桨柜 1电磁刹车动力电源线图 9 变桨机构结四风机载荷计算1叶片载荷计算（1） 作用在叶片上的离心力Fc叶片绕风轮旋转时，有离心力作用在叶片上。方向是自旋转中心沿半径想外。在半径r处，从叶片上取长为dr

21、 的一个叶素，该叶素上的离心力为dFc ，则叶片上的离心力为：错误 ! 未找到引用源。其中： 风轮角速度；r0 叶片起始处旋转半径，约为R 的 1/20,即为 1.75m；y 叶片的密度，叶片内部为空心梁结构， 查得碳纤维密度为，1800kg/m3取 900kg/m3；Ar 叶素处的叶片截面积；翼型为NACA63-415, 取 Ar 0.1c22 nr216.5601.727rad / s60用 matlab 计算得：Fc4139100N15精选文库（2） 叶轮转动时的风压力Fv：风压力是作用在叶片上沿风速方向的气动力。12R(12)(CL cos CD sin ) cdrFVvr0cot2其

22、中：c 叶素的弦长，单位为m； 来流角；设 Fv 作用点距叶轮轴的距离为rm ，则有：Rcot2)(CL cosC D sin)c rdr(1r mr 0R(1cot2C D sin)cdrr)(CL cos0用 matlab 计算得：Fv494560Nrm34.912m（3）作用在叶片上的气动力矩MbMb 是使风轮转动的力矩，可由下式求出12 R2M b2v r 0(1cot)(CL cos CD sin ) c rdrM b195330 N m（4）作用在叶片上的陀螺力矩M kM k 是风轮对风调向时产生的惯性力矩。当风向改变时， 风轮除以角速度绕水平主轴转动外，还以角速度w 垂直于水平主

23、轴并通过塔筒的轴线转动。整个叶片的转动惯量为：JR2dryAr rr 0其中：y 叶片的密度，为900kg/m3；Ar 叶素处的叶片截面积，根据翼型图得Ar0.1c2 ，单位为 m2；设科氏加速度为：16精选文库ak20x v sinrR其中：x 惯性半径，值为J / m ；r Vr 与w 的夹角， Vr =0x v ；R由动量矩定理知，叶片受到惯性力矩M k 的作用，这个力矩称为陀螺力矩。用 matlab 计算得：J = 31994000 m4M kJ ak191680000 N m2风轮载荷计算（1）轴向诱导因子 错误 !未找到引用源。 ，周向诱导因子 错误 !未找到引用源。则作用在风轮上

24、的轴向推力T 可表示为：T42 Ra)rdr2R2v 0 a(1v 0 (1k )rdr用 matlab 计算得：T663860Nm（2）作用在整个风轮上的转矩M 可表示为：M 4vR3 drvR(h - 1（)k1）3 dr0b(1 a)0rr用 matlab 计算得：M3018600N m阶段性总结表叶片上离心力 Fc4139100NFv 作用点距轴心距离34.912m叶片转动风压力 Fv494560N17精选文库气动力矩 Mb195330m4转动惯量 J31994000m4惯性力矩 Mk191680000N m风轮轴向推力 T663860 N m风轮上转矩 M3018600N m五塔架设

25、计1塔架高度塔架高度参数的选择与地形和地貌有关，陆地和海上风机组会有所不同，陆地地表相对粗糙，风速随高度变化缓慢，需要高的塔架，相反海平面较光滑，风速沿高度变化梯度大。此外，随着风轮直径的增大，塔架高度逐渐减小，对于风轮直径25m 以上的机组，轮毂中心高与风轮直径为1:1，所以选取塔架高度为90m。轮毂高度是从地面到风轮扫掠面中心的高度，用Zhub 表示zhubztH2.259092.25mzt 塔顶平面到风轮扫掠面中心的高度；H 塔架高度。2. 塔架的基本结构形式塔架形式主要分为钢筋混凝土、桁架、钢筒三种结构，目前大型风机组主要采用钢筒结构，钢筒结构中又分为直通和锥筒，在满足相同的要求的前提

26、下，锥筒可以减少材料消耗，降低成本，所以本设计选取锥形钢筒结构。同时采用柔性塔，可以降低塔架质量和成本，柔塔即塔架固有频率在风轮旋转与叶片通过频率之间。18精选文库3塔架的结构尺寸、重量和材料连接方式：法兰连接；材料： Q460c，屈服强度是 460MPa；机组质量： m1 =220t；塔筒质量： m2 = 200t；塔筒壁厚：从根部到顶部选用壁厚由30-20mm 过渡；塔筒根部外径： D1 = 6000mm；塔筒顶部外径： d1 = 4000mm；塔筒根部内径： D2 = 5940mm;塔筒顶部内径： d2= 3960mm；参照上海电气风电设备有限公司研发的3.6MW的风机，额定风速 12m

27、/s，风轮直径 116m，机舱重 150t ，塔筒分为三段，上塔筒重约48.7t ，中塔筒重约 90.4t ，下塔筒重约 86.3t ，叶片 54t ，轮毂 38t ，机组重量 467.475t 。塔筒根部外径 5056mm，顶部外径 3815mm。，4. 塔架载荷分析作用在塔架上的载荷有以下几类：（1） 风轮等构件承受的空气动力载荷；（2） 重力和惯性载荷：由重力、振动、旋转以及地震等引起的静态和动态载荷；（3） 操作载荷：在机组运行和控制过程中产生的载荷。如功率变化、偏航、变桨以及制动过程产生的载荷等；（4） 其它载荷：诸如尾迹载荷、冲击载荷、覆冰载荷等；（5） 下面只讨论与塔架结构强度计

28、算有关的两种载荷，即由风轮作用的最大气动推力以及塔架本身所承受最大风压产生的载荷；5. 塔架强度校核塔架根部为最危险截面，接下来将对根部进行受力分析和强度校核。为了确保在暴风作用条件下塔架不倾倒，校核强度时均要按照暴风工况考虑风轮的气动推力和塔架的风压力。暴风工况下最大风速 V s 和年平均风速 V ave 有关。如下表 5：表 5 年平均风速和最大风速关系V ave<=5.5m/s5.5m/s< Vave<7m/sV ave>=7m/sV s =40m/sV s=50m/sV s=60m/s本次设计中， V ave=7.0m/s>7m/s，所以选取 V s=60

29、m/s。19精选文库（1） 塔架的受力分析如下图10：图 10塔架受力塔架根部截面应力可表示为：F as (h1H ) F tsH2(G1 G2)W 22 A2其中：Fas 塔架风压力，单位 N；Fts 塔架气动推力 ,单位 N ；W2 塔架根部抗弯截面系数 ,单位为 m3；A2 塔架根部截面积，单位为m2 ；G2 塔架自重，单位为 N ；G1 机舱总成质量，单位为N；2 变截面塔架的长度折减系数，可根据2H 从图 11得到220精选文库图 112- 2关系曲线 与塔架截面变化有关的折算长度修正系数，可根据Jmin 错误 ! 未找到引用源。 之比J max由表 6 的选择参考设计值， J mi

30、n 为塔架顶部截面惯性矩， 单位为 m4， Jmax 为塔架根部截面惯性矩，单位为 m4；J min0.10.20.30.40.50.60.70.80.91.0J max1.651.451.331.241.181.141.101.061.031.00表 6-Jmin关系表J max2 塔架根部的截面惯性半径，单位为m，Jmax ；A2H 是塔架高度，即H=90m， h1 是轮毂高度，由前面知道为h1=2.25m；21精选文库（2）塔架风压力 Fas ，有以下几种计算方式: 前苏联的法捷耶夫公式:Fas0.784Ab vs2B其中：A b 叶片的投影面积，单位m2。 AbA，其中为风轮实度，风轮

31、实度与叶尖速比B有关，0 =7.5 错误 ! 未找到引用源。时，近似认为=0.05；V s 风轮中心处的暴风风速,单位为 m/s； 荷兰 ECN 的公式FasCt qAb BS其中：Ct 推力系数，取 Ct1.5；q 动态风压 ,单位为 N / m ， q 随高度变化，风轮中心高度Zhub=92.25m 处对应的q=1410N/m2。错误 !未找到引用源。 动态系数，取 =1.2；S安全系数，取 S 1.5 ； 丹麦 RIS 公式FasPl A其中：Pl 风轮单位扫掠面积上的平均风压，通常取 Pl =300N/m2； A 风轮的扫掠面积。单位为 m2；（3）塔架气动推力FtsF1v2 Ats2

32、s t其中： 空气密度， 1.225kg/m3；At 气动推力作用于塔筒的面积，单位为m2 ；22精选文库 取 0.7；V s 风轮中心处的暴风风速,单位为 m/s；（4）计算过程气动推力作用于塔筒的面积AtD1d1H6490450m 222叶片的投影面积AbA10420.05141.51m2B4334634抗弯截面系数W2D1 1D215.940.8353232D16塔架根部截面积A2D12D 22625.9420.5623m2444d24443.964塔架顶部截面惯性矩J min114d1640.495m64d1444644塔架根部截面惯性矩J maxD11D215.942.506m464

33、64D16J0.424min0.197查表 6得=1.45，Jmax2.5062.114 ,则H1.45 90,材料为低合A20.5603261.7322.114金钢，根据图5 得变截面塔架的长度折减系数2 =0.75机组的总重量为G1=2200000N ，塔架的总重量为G2=2000000N苏联的法捷耶夫公式：Fas10.784 Ab vs2 B0.784141.51602 31198193 N荷兰 ECN 的公式：FasCqA B S1.51410141.5131.2Ntb1.5 16161857丹麦 RIS 公式： Fasp1 Ad3005222547168N选取丹麦 RIS 公式得出的

34、 Fas = 2547.168KN塔架气动推力 : Fts1vs2 At11.2256020.7 396611226 N22代入计算得：23精选文库Fas (1 H )FtsH(G1 G2)2547168 92.25 611226 90( 22000002000000)h22W 22 A20.8350.75324MPa0.5623强度校核：s460MPa324MPa所以塔架最危险截面，即根部满足要求。六风电机组布局见附录整机工程图七设计总结1总体参数风机总体设计参数总结，如表7表 7总体参数表叶片数B=3风轮直径额定输出功率P=3.5MW扫掠面积设计寿命20 年风能利用系数切入风速V in=3

35、m/s叶尖速比切出风速V out功率控制方式=25m/s额定风速V r=12m/s制动形式风轮额定转速nr=16.5r/min传动系统D=104mA=8482m2Cp=0.450 7.5主动变桨距控制气动刹车、机械刹车高传动比齿轮箱传动24精选文库风力机等级IECIIIA电气系统双馈发电机 +变流器叶片数B=3风轮直径D=70m2功率、气动特性和载荷计算总结风机的功率、叶片气动特性和载荷总结，如表8表 8 风机的功率、叶片气动特性和载荷攻角5.25°气动特性升力系数0.9461阻力系数0.00791发电机额定功率3.5MW发电机发电机额定转速1500r/min发电机轴直径0.18m形

36、式行星齿轮传动两级 NGW齿轮箱传动比90齿轮箱低速轴转速16.5r/min高速轴转速1500r/min齿轮箱功率3.9MW联轴器低速轴联轴器功率3.9MW高速轴联轴器功率3.84MW变流器变流器功率1.5MW变桨电机数量 变桨电机功率3 10KW偏航电机数量 错误 ! 未找到引用源。6 3KW偏航电机功率离心力4139.1KN叶片载荷风压力494.56N气动力矩195330N·m风压力作用点距叶轮轴距离34.912m风轮载荷轴向推力663860N转矩3018600N· m主轴载荷低速轴转矩2258251.3N·m25精选文库高速轴转矩24458.6N·m低速轴直径0.7m高速轴直径0.18m作用在塔架上载荷Fas=2547168N塔架载荷作用在在塔架上风压力Fts=611226N塔架根部截面应力325MPa3其余部分设计总结整机部分设计总结，如表9图 9 整机部分设计形式机组运行环境温度-30°45°生存环境温度-45°45