学位论文-—风力发电机组设计与制造课程设计.doc

课程设计(综合实验)报告
( 2012 – 2013 年度第二学期)

名 称： 风力发电机组设计与制造
题 目： 风电机组总体技术设计
院 系： 可 再生 能 源 学 院
班 级： 风 能 1002 班
学 号： 1101540209
学生姓名： 刘 凤 魁
指导教师： 田 德、王 永
设计周数： 19周-20周

成 绩：
日 期： 2013年 7 月10 日

目 录
课程设计任务书 1
第一章 风电机组整体参数设计 5
1.1 额定功率 5
1.2 设计寿命 5
1.3 切出风速、切入风速、额定风速 5
1.4 发电机额定转速和转速范围 5
1.5 重要几何尺寸 6
1.5.1 风轮直径和扫掠面积 6
1.5.2 轮毂高度 7
1.6 叶片数 7
1.7 风轮转速 8
1.8 功率曲线、Cp曲线、Ct曲线、攻角ɑ 8
1.9 功率控制方式 14
1.10 制动系统形式 14
1.11 风力机等级 14
第二章 风电机组气动特性初步计算 14
2.1 叶片的设计理论 14
2.1.1 动量理论 14
2.1.2 叶素理论 15
2.2 叶片设计 16
第三章 机组及部件载荷计算 18
3.1 叶片载荷计算 18
3.1.1 作用在叶片上的离心力FC 18
3.1.2 作用在叶片上的风压力FV 19
3.1.3 作用在叶片上的气动力矩Mb 19
3.1.4 作用在叶片上的陀螺力矩MK 19
3.2 风轮轮毂的载荷计算 20
3.3 主轴载荷计算 20
3.5 塔架载荷计算 21
3.5.1 暴风工况的风轮气动推力计算 21
3.5.2 欧美国家塔架静态强度设计的一般载荷条件 22
3.5.3 确定塔架设计载荷的要求 23
3.6 联轴器载荷计算 23
第四章 主要部件技术参数 23
4.1 齿轮箱 23
4.2 发电机 24
4.3 变流器 24
4.4 联轴器 25
4.5 偏航执行机构 25
4.6 变桨距执行机构 26
4.7 塔架根部截面应力计算 26
风力机总体布局图： 28
参考文献： 28

风力发电机组设计与制造课程设计
课程设计任务书
一、 设计内容
风电机组总体技术设计
二、 目的与任务
主要目的：
1. 以大型水平轴风力机为研究对象，掌握系统的总体设计方法；
2. 熟悉相关的工程设计软件；
3. 掌握科研报告的撰写方法。
主要任务：
每位同学独立完成风电机组总体技术设计，包括：
1. 确定风电机组的总体技术参数；
2. 关键零部件（齿轮箱、发电机和变流器）技术参数；
3. 计算关键零部件（叶片、风轮、主轴、连轴器和塔架等）载荷和技术参数；
4. 完成叶片设计任务；
5. 确定塔架的设计方案。
6. 每人撰写一份课程设计报告。
三、 主要内容
每人选择功率范围在1.5MW至6MW之间的风电机组进行设计。
1）原始参数：风力机的安装场地50米高度年平均风速为7.0m/s，60米高度年平均风速为7.3m/s，70米高度年平均风速为7.6 m/s，当地历史最大风速为49m/s，用户希望安装1.5 MW至6MW之间的风力机。采用63418翼型，63418翼型的升力系数、阻力系数数据如表1所示。空气密度设定为1.225kg/m3。
2）设计内容
（1）本参数包括叶片数、风轮直径、额定风速、切入风速、切出风速、功率控制方式、传动系统、电气系统、制动系统形式和塔架高度等，根据标准确定风力机等级；
（2）关键部件气动载荷的计算。设定几种风轮的Cp曲线和Ct曲线，计算几种关键零部件的载荷（叶片载荷、风轮载荷、主轴载荷、连轴器载荷和塔架载荷等）；根据载荷和功率确定所选定机型主要部件的技术参数（齿轮箱、发电机、变流器、连轴器、偏航和变桨距电机等）和型式。以上内容建议用计算机编程实现，确定整机和各部件（系统）的主要技术参数。
（3）塔架根部截面应力计算。计算暴风工况下风轮的气动推力，参考风电机组的整体设计参数，计算塔架根部截面的应力。最后提交有关的分析计算报告。

四、 进度计划
序号
设计(实验)内容
完成时间
备注
1
风电机组整体参数设计
2.5天

2
风电机组气动特性初步计算
2天

3
机组及部件载荷计算
2天

4
齿轮箱、发电机、变流器技术参数
1.5天

4
塔架根部截面应力计算
1天

5
报告撰写
1.5天

6
课程设计答辩
1.5天

五、 设计（实验）成果要求
1． 提供设计的风电机组的性能计算结果；
2． 绘制整机总体布局工程图。
六、 考核方式
每人提交一份课程设计报告；准备课程设计PPT，答辩。
2013年 6 月 28 日

表1 63418翼型的升力系数和阻力系数
攻角
升力系数
阻力系数
扭矩系数
-11.04
-0.77
0.031
0
-9
-0.745
0.012
0
-8
-0.626
0.01
0
-7
-0.508
0.01
0
-6
-0.39
0.009
0
-5
-0.263
0.008
0
-4
-0.138
0.007
0
-3
-0.03
0.006
0
-2
0.079
0.006
0
-1
0.175
0.006
0
0
0.313
0.006
0
1
0.446
0.006
0
2
0.551
0.006
0
3
0.666
0.007
0
4
0.777
0.007
0
5
0.904
0.007
0
6
0.984
0.008
0
7
1.074
0.01
0
8
1.17
0.013
0
9
1.244
0.015
0
10
1.307
0.018
0
12
1.395
0.025
0
14
1.395
0.032
0
16
1.37
0.039
0
20
1
0.185
0
25
1.04
0.315
0
30
1.02
0.442
0
40
0.95
0.68
0
45
0.9
0.78
0
50
0.82
0.891
0
55.54
0.697
1.022
0
60.67
0.626
1.095
0
66.93
0.495
1.144
0
71.07
0.401
1.179
0
75.5
0.303
1.235
0
80.75
0.199
1.268
0

攻角
升力系数
阻力系数
扭矩系数
-180
-0.475
0.028
0
-179.9
-0.475
0.028
0
-177.79
0.083
0.039
0
-172.18
0.701
0.07
0
-168.55
0.917
0.172
0
-165.1
0.661
0.268
0
-159.59
0.477
0.344
0
-154.91
0.612
0.426
0
-150.61
0.697
0.536
0
-145.49
0.757
0.681
0
-141.27
0.716
0.772
0
-134.54
0.643
0.829
0
-127.58
0.572
0.886
0
-124.17
0.537
0.938
0
-120.56
0.466
1
0
-115.24
0.361
1.069
0
-109.23
0.287
1.108
0
-104.33
0.179
1.145
0
-100.65
0.094
1.167
0
-96.72
0.023
1.183
0
-88.29
-0.099
1.18
0
-82.73
-0.192
1.166
0
-78.79
-0.261
1.148
0
-74.68
-0.332
1.112
0
-68.72
-0.437
1.044
0
-62.9
-0.538
1.003
0
-58.21
-0.597
0.965
0
-53.67
-0.656
0.902
0
-50.77
-0.673
0.825
0
-44.08
-0.68
0.702
0
-38.09
-0.698
0.602
0
-32.99
-0.655
0.508
0
-29.64
-0.535
0.429
0
-24.44
-0.386
0.28
0
-19
-0.354
0.191
0
-12.78
-0.463
0.092
0
-11.58
-0.8
0.036
0

攻角
升力系数
阻力系数
扭矩系数
86.9
0.099
1.272
0
91.87
-0.014
1.265
0
95.79
-0.112
1.26
0
99.84
-0.214
1.254
0
104.93
-0.317
1.205
0
110.44
-0.427
1.127
0
115.85
-0.539
1.094
0
124.05
-0.629
0.978
0
129.73
-0.695
0.872
0
136.34
-0.796
0.805
0
140.53
-0.77
0.732
0
144.46
-0.729
0.633
0
145.07
-0.721
0.616
0
150.54
-0.674
0.466
0
155.32
-0.597
0.308
0
160.62
-0.843
0.271
0
163.66
-1.124
0.183
0
168.46
-0.88
0.073
0
175.01
-0.475
0.028
0
179.9
0.083
0.039
0
180
0.083
0.039
0

第一章 风电机组整体参数设计
1.1 额定功率
根据《设计任务书》选定额定功率为3MW。
1.2 设计寿命
风力发电机组安全等级I～III的设计寿命至少为20年。这里选20年设计寿命。
1.3 切出风速、切入风速、额定风速
切入风速 取 vin=3m/s；
切出风速 取 vout=25m/s；
额定风速 取 vr=13m/s；
对于一般变桨距风力发电机组（选3MW）的额定风速羽平均风速之比为1.70左右。
1.4 发电机额定转速和转速范围
应用高速发电机的风力发电机组多为异步发电机。根据电机学理论，当异步发电机接入频率恒定的电网上时，由定子三相绕组中的电流产生的旋转磁场的同步速取决于电网的额定频率和发电机绕组的极对数，三者关系为，
其中：
n1—同步转速，单位为r/min；
f1—电网频率，单位为Hz；
P—发电机绕组的极对数。
异步发电机中旋转场和转子之间的相对转速为，相对转速与同步速的比值称为异步发电机的转差率，用s表示，即。
当异步发电机的转子在风力机的拖动下，以高于同步速旋转时，发电机运行在发电状态，发电机中的电磁转矩为制动转矩，阻碍发电机旋转，此时发电机需要从外部吸收无功电流建立磁场（如由电容提供无功电流），而将从风力机中获得的机械能转化为电能提供给电网。喜事发电机的转差率为负值，一般其绝对值在2%~5%之间，并网运行的较大容量的异步发电机的转子转速一般在（1~1.05）n1之间。
应用双馈式发电机，转子转速一般在（130%）n1之间变化。
1.5 重要几何尺寸
1.5.1 风轮直径和扫掠面积
由风力发电机组输出功率得叶片直径：

其中：
Pr—风力发电机组额定输出功率，取3000kW；
—空气密度（一般取标准大气状态），取1.225kg/m3；
Vr—额定风速，取13m/s；
D—风轮直径，当额定风速为13m/s时，风轮直径为97m；
—传动系统效率，取0.96；
—发电系统效率，取0.94（发电机效率0.97，变流器效率0.97）；
Cp--额定功率下风能利用系数，取0.33。
由直径计算可得扫掠面积：
由大型风力发电机组风轮扫掠面积与额定功率关系的经验曲线图可知：每平方米风轮扫掠面积产生的额定功率为405W/m2，
风轮扫掠面积与额定功率关系的经验曲线图
综上可得风轮直径D=97m，扫掠面积S=7390m2。
1.5.2 轮毂高度
轮毂高度是从地面到风轮扫掠面中心的高度，用Zhub表示，单位为m。
Zhub=Zt+H，
其中：
H—塔架高度，单位为m；
Zt—塔顶平面到风轮扫掠面中心的高度，单位为m。
对于风轮直径25m以上的机组，其轮毂中心高度与风轮直径的比基本为1：1，
取Zhub=97m。
1.6 叶片数
现代风力发电机的实度比较小，一般需要1-3个叶片。选择风轮叶片数时考虑风电机组性能和载荷、风轮和传动成本、风力机气动噪声及景观影响等因素。
3叶片较1、2叶片风轮有如下优点：
1. 平衡简单、动态载荷小。基本消除了系统的周期载荷，输出较稳定转矩；
2. 能提供较佳的效率；
3. 更加美观；
4. 噪声较小；
5. 轮毂较简单等。
综上所述，叶片数选择B=3。
1.7 风轮转速
风轮的输出功率与其转矩和转速有关。由风轮转矩M=9550PS/n可知，当额定功率PS(kW)已选定时，风轮转速n（r/min）与转矩（N·m）成反比。为降低转矩M，应提高风轮转速。但转速过高对风轮不利，所以应选择一个合适的风轮转速。通常将陆上风力发电机组的叶尖速限制在65m/s左右，三叶片风力发电机组的风轮叶尖速比一般在6.0至8.0之间。此处取6.5。

所以，风轮额定转速为17r/min。

1.8 功率曲线、Cp曲线、Ct曲线、攻角ɑ

其中：
к为威布尔分布形状参数。根据大量的观测，中国地区k值通常在1.0-2.6之间。这里取k=2.5。
风能利用系数Cp计算公式如下：

推力系数CT计算公式如下：

表格 1
风速v（m/s）
额定功率Pr（kw）
输出功率Pw
（kw）
Cp
a
Ct
3
3000
0.0000
0.0000
0.0000
0.0000
4
3000
114.9352
0.4399
0.1920
0.6205
5
3000
234.3764
0.4593
0.2991
0.8386
6
3000
395.6851
0.4487
0.2667
0.7823
7
3000
602.8693
0.4305
0.2202
0.6868
8
3000
859.6122
0.4113
0.1858
0.6052
9
3000
1169.3409
0.3929
0.1610
0.5404
10
3000
1535.2735
0.3761
0.1419
0.4869
11
3000
1960.4531
0.3608
0.1267
0.4425
12
3000
2447.7736
0.3470
0.1144
0.4051
13
3000
3000.0000
0.3345
0.1042
0.3734
14
3000
3000.0000
0.2678
0.0789
0.2907
15
3000
3000.0000
0.2177
0.0617
0.2317
16
3000
3000.0000
0.1794
0.0496
0.1886
17
3000
3000.0000
0.1496
0.0406
0.1558
18
3000
3000.0000
0.1260
0.0337
0.1303
19
3000
3000.0000
0.1071
0.0284
0.1102
20
3000
3000.0000
0.0919
0.0241
0.0941
21
3000
3000.0000
0.0793
0.0207
0.0810
22
3000
3000.0000
0.0690
0.0179
0.0703
23
3000
3000.0000
0.0604
0.0156
0.0614
24
3000
3000.0000
0.0532
0.0137
0.0539
25
3000
3000.0000
0.0470
0.0120
0.0476

表格 2
攻角
升力系数
阻力系数
升阻比
攻角
升力系数
阻力系数
升阻比
-180
-0.475
0.028
-16.964
0
0.313
0.006
52.167
-179.9
-0.475
0.028
-16.964
1
0.446
0.006
74.333
-177.79
0.083
0.039
2.128
2
0.551
0.006
91.833
-172.18
0.701
0.07
10.014
3
0.666
0.007
95.143
-168.55
0.917
0.172
5.331
4
0.777
0.007
111.000
-165.1
0.661
0.268
2.466
5
0.904
0.007
129.143
-159.59
0.477
0.344
1.387
6
0.984
0.008
123.000
-154.91
0.612
0.426
1.437
7
1.074
0.01
107.400
-150.61
0.697
0.536
1.300
8
1.17
0.013
90.000
-145.49
0.757
0.681
1.112
9
1.244
0.015
82.933
-141.27
0.716
0.772
0.927
10
1.307
0.018
72.611
-134.54
0.643
0.829
0.776
12
1.395
0.025
55.800
-127.58
0.572
0.886
0.646
14
1.395
0.032
43.594
-124.17
0.537
0.938
0.572
16
1.37
0.039
35.128
-120.56
0.466
1
0.466
20
1
0.185
5.405
-115.24
0.361
1.069
0.338
25
1.04
0.315
3.302
-109.23
0.287
1.108
0.259
30
1.02
0.442
2.308
-104.33
0.179
1.145
0.156
40
0.95
0.68
1.397
-100.65
0.094
1.167
0.081
45
0.9
0.78
1.154
-96.72
0.023
1.183
0.019
50
0.82
0.891
0.920
-88.29
-0.099
1.18
-0.084
55.54
0.697
1.022
0.682
-82.73
-0.192
1.166
-0.165
60.67
0.626
1.095
0.572
-78.79
-0.261
1.148
-0.227
66.93
0.495
1.144
0.433
-74.68
-0.332
1.112
-0.299
71.07
0.401
1.179
0.340
-68.72
-0.437
1.044
-0.419
75.5
0.303
1.235
0.245
-62.9
-0.538
1.003
-0.536
80.75
0.199
1.268
0.157
-58.21
-0.597
0.965
-0.619
86.9
0.099
1.272
0.078
-53.67
-0.656
0.902
-0.727
91.87
-0.014
1.265
-0.011
-50.77
-0.673
0.825
-0.816
95.79
-0.112
1.26
-0.089
-44.08
-0.68
0.702
-0.969
99.84
-0.214
1.254
-0.171
-38.09
-0.698
0.602
-1.159
104.93
-0.317
1.205
-0.263
-32.99
-0.655
0.508
-1.289
110.44
-0.427
1.127
-0.379
-29.64
-0.535
0.429
-1.247
115.85
-0.539
1.094
-0.493
-24.44
-0.386
0.28
-1.379
124.05
-0.629
0.978
-0.643
-19
-0.354
0.191
-1.853
129.73
-0.695
0.872
-0.797
-12.78
-0.463
0.092
-5.033
136.34
-0.796
0.805
-0.989
-11.58
-0.8
0.036
-22.222
140.53
-0.77
0.732
-1.052
-11.04
-0.77
0.031
-24.839
144.46
-0.729
0.633
-1.152
-9
-0.745
0.012
-62.083
145.07
-0.721
0.616
-1.170
-8
-0.626
0.01
-62.600
150.54
-0.674
0.466
-1.446
-7
-0.508
0.01
-50.800
155.32
-0.597
0.308
-1.938
-6
-0.39
0.009
-43.333
160.62
-0.843
0.271
-3.111
-5
-0.263
0.008
-32.875
163.66
-1.124
0.183
-6.142
-4
-0.138
0.007
-19.714
168.46
-0.88
0.073
-12.055
-3
-0.03
0.006
-5.000
175.01
-0.475
0.028
-16.964
-2
0.079
0.006
13.167
179.9
0.083
0.039
2.128
-1
0.175
0.006
29.167
180
0.083
0.039
2.128

结合图表可以看出，攻角约等于5度时，升阻比最大为129.143，此时升力系数Cl=0.904，阻力系数Cd=0.007，本设计中攻角取9度，升力系数Cl=1.244，阻力系数Cd=0.015，升阻比为82.933。
1.9 功率控制方式
采用主动变桨距控制。
1.10 制动系统形式
第一制动采用气动刹车，第二制动采用高速轴机械刹车。
1.11 风力机等级
由IEC标准，如下表，选择风力机等级为IECIIA。
WTGS等级
I
II
III
S
vref[m/s]
50
42.5
37.5
由设计者规定各参数
A Iref（—）
0.16
B Iref（—）
0.14
C Iref（—）
0.12

第二章 风电机组气动特性初步计算
2.1 叶片的设计理论
2.1.1 动量理论
如图4所示，由动量理论对风轮进行计算：
在风轮r—r+dr的环域内应用动量定理，则作用在环域上的
轴向推力为： （2-1）
轴向推力又可表示为： （2-2）
利用伯努利方程式，
二式相减并带入（2-2）可得： （2-3）
联立（2-1）、（2-3）可得： （2-4）
引入轴向诱导因子，则 （2-5）
（2-6）气流流过风轮后，还产生与风轮转动方向相反的角速度，因而作用在环域上的转矩为： （2-7）
引入切向诱导因子，则转矩表示为：
（2-8）
从而可知，作用在环域上的气流提供的动量为：
（2-9）
2.1.2 叶素理论

如上图所示，设，则 （2-10）
因此，升力和阻力可表示为：假设，则：
那么，作用在风轮r—r+dr的环域上的推理和转矩分别为：
（2-11）
（2-12）
2.2 叶片设计
1) 给出r，攻角i已知，CL=0.904已知；
2) ；
3) 由式，得到角；
4) 由式，得到k；
5) 由式，得到h；
6) 由式，得到；
7) ；
8) ；
将1）~8）在excel中编辑公式计算，得出相应的数据。
其中：
B—叶片数（此处为3）；
k—轴向干扰因子；
h—切向干扰因子；
C—叶片弦长；
i—攻角；
—入流角；
—桨矩角。
表格 3
叶素位置/%
5
10
15
20
25
30
35
40
45
50

6.5
6.5
6.5
6.5
6.5
6.5
6.5
6.5
6.5
6.5
叶尖速比
0.325
0.65
0.975
1.3
1.625
1.95
2.275
2.6
2.925
3.25
叶素位置/%
55
60
65
70
75
80
85
90
95
100

6.5
6.5
6.5
6.5
6.5
6.5
6.5
6.5
6.5
6.5
叶尖速比
3.575
3.9
4.225
4.55
4.875
5.2
5.525
5.85
6.175
6.5

表格 4
r
位置/%

尖速比

k
h

c
2.425
5
6.5
0.325
1.152
0.428
2.956
0.838
39.038
5.399
4.85
10
6.5
0.65
1.239
0.388
1.735
0.663
29.006
6.914
7.275
15
6.5
0.975
1.305
0.367
1.382
0.532
21.497
6.768
9.7
20
6.5
1.3
1.352
0.356
1.232
0.437
16.051
6.149
12.125
25
6.5
1.625
1.387
0.349
1.154
0.368
12.096
5.444
14.55
30
6.5
1.95
1.413
0.345
1.110
0.316
9.115
4.853
16.975
35
6.5
2.275
1.433
0.342
1.082
0.276
6.822
4.330
19.4
40
6.5
2.6
1.448
0.340
1.063
0.245
5.045
3.869
21.825
45
6.5
2.925
1.461
0.339
1.050
0.220
3.611
3.497
24.25
50
6.5
3.25
1.471
0.338
1.041
0.199
2.408
3.214
26.675
55
6.5
3.575
1.480
0.337
1.034
0.182
1.433
2.952
29.1
60
6.5
3.9
1.487
0.336
1.029
0.167
0.573
2.761
31.525
65
6.5
4.225
1.493
0.336
1.025
0.155
-0.115
2.588
33.95
70
6.5
4.55
1.499
0.336
1.021
0.144
-0.745
2.350
36.375
75
6.5
4.875
1.503
0.335
1.019
0.135
-1.261
2.283
38.8
80
6.5
5.2
1.507
0.335
1.016
0.127
-1.720
2.056
41.225
85
6.5
5.525
1.511
0.335
1.014
0.119
-2.178
1.915
43.65
90
6.5
5.85
1.514
0.335
1.013
0.113
-2.522
1.885
46.075
95
6.5
6.175
1.517
0.335
1.012
0.107
-2.866
1.839
48.5
100
6.5
6.5
1.520
0.334
1.010
0.102
-3.153
1.616

第三章 机组及部件载荷计算
3.1 叶片载荷计算
3.1.1 作用在叶片上的离心力FC
叶片绕风轮轴旋转是，有离心力作用在叶片上。离心力的方向是自旋转中心沿半径向外。则叶片的离心力为：
其中：
—叶片的密度，单位为kg/m3 ，本设计为530 kg/m3；
—叶素处的叶片截面积，单位为m2；
—风轮角速度，单位为rad/s；
—叶片起始处的旋转半径；
—叶片结束处的旋转半径。
利用gambit软件粗略计算出，用excel表格计算出各截面的面积，然后再利用matlab软件拟合出半径与叶素截面面积的函数，计算积分，求得 39786kN。

r
2.425
4.85
7.275
9.7
12.125
14.55
16.975
19.4
21.825
24.25
c
5.399
6.914
6.768
6.149
5.444
4.853
4.330
3.869
3.497
3.214
Ar
2.43
3.98
3.81
3.15
2.47
1.96
1.56
1.25
1.02
0.86
r
26.675
29.1
31.525
33.95
36.375
38.8
41.225
43.65
46.075
48.5
c
2.952
2.761
2.588
2.350
2.283
2.056
1.915
1.885
1.839
1.616
Ar
0.73
0.63
0.56
0.46
0.43
0.35
0.31
0.30
0.28
0.22

3.1.2 作用在叶片上的风压力FV
风压力是作用在叶片上沿风速方向的气动力。风轮静止和转动时，风压力大小不相等。
风轮静止时风压力为：
设FV的作用点距风轮轴距离为，则有：

风轮转动时的风压力为：

用Matlab计算得出FV1=289.3962N，rm1=14.01m；FV2=1969.4952kN，rm2=33.95m。
3.1.3 作用在叶片上的气动力矩Mb
Mb是使风轮转动的力矩，由下式求出：

用matlab计算得出：Mb =3568.7628kN·m
3.1.4 作用在叶片上的陀螺力矩MK
Mk是风轮对风调向时产生的惯性力矩。当风向改变时，风轮除以角速度绕Ox轴转动外，还要以角速度绕Oz轴转动。整个叶片的转动惯量为
对风调向时产生的最大科氏加速度为：由于科氏加速度，产生了柯氏角加速度，由动量矩定理知，叶片收到惯性力矩Mk的作用，这个力矩称为陀螺力矩，Mk=。
用matlab计算得出J=512.47Mkg*m2；，Mk=23.3686MN·m
3.2 风轮轮毂的载荷计算
作用在整个风轮上的轴向推力可表示为。
作用在整个风轮上的转矩可表示为：。
用用matlab计算得出T =627.4529kN；M