《风力发电机组设计与制造》课程设计任务书Word格式.docx
- 文档编号:21046655
- 上传时间:2023-01-27
- 格式:DOCX
- 页数:23
- 大小:178.12KB
《风力发电机组设计与制造》课程设计任务书Word格式.docx
《《风力发电机组设计与制造》课程设计任务书Word格式.docx》由会员分享,可在线阅读,更多相关《《风力发电机组设计与制造》课程设计任务书Word格式.docx(23页珍藏版)》请在冰豆网上搜索。
天
2
风电机组气动特性初步计算
2天
3
机组及部件载荷计算
|
4
齿轮箱、发电机、变流器技术参数
塔架根部截面应力计算
*
1天
5
报告撰写
6
课程设计答辩
五、设计(实验)成果要求
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取
发电机效率η2取
变流器效率η3取
7、风能利用系数Cp
额定功率下风能利用系数Cp取
8、风轮直径
由风力发电机组输出功率得叶片直径:
D=
=
=92m
其中:
Pr--风力发电机组额定输出功率,取3000kW;
--空气密度(一般取标准大气状态),取m3;
Vr--额定风速,取13m/s;
D--风轮直径;
η1--传动系统效率,取;
η2--发电机效率,取;
η3--变流器效率,取;
Cp--额定功率下风能利用系数,取。
9、功率Pw曲线,风能利用系数CP和推力系数CT曲线的确定
风电机组的设计需要给出功率输出特性曲线。
如果确定了风电机组的切入风速,额定风速和切出风速,可以利用下面的方程计算输出功率,查阅相关资料,得到其风能利用系数随着风速的变化值,得到下表。
(其中所用公式为:
Cp=4a(1-a)2Ct=4a(1-a)
`
D=
)
利用上述的公式,由excel的函数运算,计算得到随着风速的变化,Cp、和Ct的值得变化
风速v[m/s]
功率Pw[kW]
风能利用系数Cp
轴向诱导因子a
推力系数CT
叶尖速比
¥
22.
16.
>
13.
11.
】
7
8
—
9
10
、
11
12
}
5.
13
14
?
15
16
17
18
》
19
20
【
21
22
23
24
.
25
~
10、风轮转速
风力发电机组产生的气动噪声正比于叶尖速度的5次方。
通常将路基风力发电机组的叶尖速度限制在65m/s左右。
这里取叶尖速V=65m/s。
]
由v=60得到n=60v/,将V=65m/s带入得,r/min。
11、叶尖速比λ
由公式,带入额定转速n=min可以得到
参考CP-λ曲线,当Cp=时对应的λ的值约为5,因此所选数值合理。
12、功率控制方式
功率控制方式选取当今应用较为广泛的变浆距控制。
这种机组当风速过高时,通过减小叶片翼型上合成气流与翼型几何弦的夹角(攻角),改变风力发电机组获得的空气动力转矩,能使功率输出保持稳定。
同时,风机在起动过程也需要通过变浆距来获得足够的起动转矩。
采用变浆距技术的风力发电机组还可以使叶片和整机的受力状况大为改善,这对于大型风力发电机组十分的有利。
13、传动系统
[
传动系统选定由以下四部分组成:
主轴系统、增速传动系统(齿轮箱)、轴系的支撑与连接(如轴承和联轴器)、制动装置。
14、制动系统
制动系统的形式以空气动力制动和机械制动为主。
15、电气系统
风电机组的电气系统包括以下四个主要部分:
变浆距驱动,偏航驱动,发电机部分,控制和监控部分
16、塔架高度
由于风速与距地面高度有关,增加塔架高度可使风轮获取更多的风能,但制造更高的塔架也需要更多的材料,使其造价相应增加。
大型机组的塔架高度H可以按下式初步确定:
。
式中D为风轮直径,单位为m。
塔架的高度参数的选择与地形和地貌有关,如陆地与海上风电机组会有所不同。
陆地的地表性对粗糙且风速随高度变化缓慢,可能需要较高的塔架;
而海平面较光滑,风速沿着高度方向变化梯度大。
一般而论,即使相同容量的机组,在近海安装的机组塔架高度相对较低。
考虑地貌因素的塔架最低设计高度一般可用下式估算:
H=h+C+R
式中:
h——机组附近的障碍物高度,单位m;
C——障碍物最高点到风轮扫掠面最低点的距离(最小取),单位m。
综上考虑,塔架高度H取90m
17、设计寿命
按照惯例,设计寿命取20年。
18、设计标准
设计标准为IEC61400-1
19、风力机等级
70m年平均风速为s。
去IECⅢA,参考下表:
风电机组等级
第Ⅰ级
第Ⅱ级
第Ⅲ级
第Ⅳ级
参考风速Vr=5Va(m/s)
50
30
年平均风速Va(m/s)
·
湍流强度等级
A
ε15
a
B
b
阶段性总结
叶片数
切入风速
3m/s
功率控制
变浆距控制
切出风速
25m/s
制动系统形式
空气动力制动和机械制动
额定风速
13m/s
塔架高度
90m
额定功率
3MW
设计寿命
20年
风能利用系数
设计标准
IEC61400
风轮转速
min
风力机等级
IECⅢA
风轮直径
92m
/
二、关键部件气动载荷的计算
包括设定几种风轮的Cp曲线和Ct曲线,计算几种关键零部件的载荷(叶片载荷、风轮载荷、主轴载荷、连轴器载荷和塔架载荷等);
最后提交有关的分析报告。
根据所提供的63418翼型的升力系数、阻力系数数据,得到攻角与升力系数、阻力系数的关系曲线。
由图可以观察出,攻角约等于13°
时升力系数最大,大约为左右。
攻角约等于9°
时,升力系数约为最大升力系数的倍,为,此时的阻力系数为。
当攻角为5°
时升阻比最大,此时升力系数为。
在叶片处选择10段截面,分别在10%R,20%R,30%R……90%R,100%R。
由以下公式:
a、给出r,攻角i已知,CL=已知;
b、由式
求得
角;
c、由式
得到k;
d、{
e、由式
得到h;
f、由式
得到I;
g、
,
;
h、将a--g编程进行运算得到数据--表2。
(程序见附录)
通过excel计算得到以上各个参数的数值:
位置
半径r
&
λ
Ψ
k
h
φ
β
C
10%
20%
30%
40%
50%
60%
70%
35
80%
40
90%
45
$
100%
1、主要部件功率计算
1)发电机
发电机类型:
双馈异步变速恒频式发电机
额定功率:
P发电机=Pr/η3*=3000/*=3402kW(预留10%的裕度)
极对数:
由公式60f=np可以得到发电机的额定转速为1650m/min(预留10%的裕度)
2)齿轮箱
齿轮箱类型:
增速,2级行星齿轮+1级平行齿轮
齿轮箱效率
传动比:
120
低速轴转速:
r/min;
高速轴转速:
1650r/min
齿轮箱功率:
PGB=3581kw
3)联轴器
低速轴联轴器动率:
3400kW
高速轴联轴器动率:
3500kW
4)变流器
变流器的功率通常为风电机组的额定功率的1/2~1/3,考虑到风电机组的可靠性,通常为风电机组额定功率的1/2,即变流器的功率为1500kW。
5)偏航电机
偏航电机的功率由最大偏航扭矩确定。
在此先初步选定4台4kW的偏航电机。
2、主要部件的载荷计算
1)《
2)叶片的载荷计算
叶片绕风轮轴旋转时,有离心力作用在叶片上。
离心力的方向是自旋转中心沿半径向外。
作用在叶片上的离心力Fc可按照下式求解:
ρy——叶片的密度,单位是kg/m3,取m3
Ar——叶素处的叶片截面面积,单位是m2
Ω——风轮的角速度,单位是rad/s
r0——叶片起始处的旋转半径
R——叶片结束处的旋转半径
其中:
Ω=rad/s
由matlab计算的:
Fc=00N。
3)作用在叶片上的风压力
风压力是作用在叶片上沿风速方向的气动力。
风轮转动时的风压力:
设FV的作用点几粒风轮轴的距离为rm,则:
由matlab计算得:
Fv=755580N,rm=。
4)作用在叶片上的气动力矩Mb
Mb是使风轮转动的力矩,由matlab计算得:
Mb=144250N·
m。
5)作用在叶片上的陀螺力矩
Mk是风轮对风调向时产生的惯性力矩。
当风向改变时,风轮除以角速度Ω绕OX轴转动外,还要以角速度ω绕OZ轴转动。
整个叶片的转动惯量为:
由matlab计算得,J=0kg·
m2
对风调向时,叶片可看成以角速度ω绕OZ轴转动的牵连运动、以Ω绕OX轴转动的相对运动的复合运动,这时要产生柯式加速度ak:
=2w*v
这里是惯性半径,β是Vr与ω的夹角。
当叶片旋转到铅垂位置时,β=0°
,柯式加速度达到最大值:
由于柯式加速度,就产生了柯式角加速度
由动量矩定理知,叶片受到的惯性力矩Mk的作用,这个力矩称为陀螺力矩:
由matlab计算得,Mk=00N·
m
3、风轮载荷计算
轴向诱导因子,周向诱导因子
作用在整个风轮上的轴向推力
作用在整个风轮上的转矩由matlab计算得,T=718370N,M=4613600N·
4、(
5、主轴载荷计算
低速轴角速度为:
s
高速轴角速度为:
低速轴功率为:
高速轴功率为:
低速轴转矩为:
高速轴转矩为:
假设所选择的实心钢轴的推荐最大应力为55MPa,则:
低速轴的直径:
高速轴的直径:
6、塔架载荷计算
除了风载荷外,风电机组的几乎所有的载荷也都将传给塔架,按风电机组的载荷源分类又要有:
风轮等构件承受的空气动力载荷
重力和惯性载荷:
由重力、振荡、旋转以及地震等引起的静态和动态载荷
操作载荷:
在机组运行和控制过程中产生的载荷。
如功率变化、偏航、变浆以及制动过程产生的载荷等
其他载荷:
诸如尾迹载荷、冲击载荷、覆冰载荷等、
下面只讨论与塔架结构强度计算有关的两种载荷:
即由风轮作用的最大气动推力以及塔架本身所承受最大风压产生的载荷。
1)暴风工况的风轮气动推力计算
暴风工况所对应的最大风速vs与当地的年平均风速vave有关。
由于当地的年平均风速为s,所以vs取60m/s。
前苏联的法捷耶夫公式:
Fas=
As——叶片的投影面积,单位为m2
Vs——风轮中心处的暴风风速,单位为m/s
B——风轮叶片数
@
叶片的投影面积
其中σ为风轮的实度,风轮实度与叶尖速比有关,当λ=5时,可以近似的认为σ=
由matlab计算得,Fas=1108400N。
荷兰ECN的公式
Ct——推力系数,取Ct=
q——动态风压,单位是N/m2
ψ——动态系数,取ψ=
S——安全系数,取S=
其中的q是随着高度变化,这里q的值为1390N/m2
由matlab计算得,Fas=1473800N
德国DFVLR公式
vs——风轮中心处的暴风风速,取vs=60m/s
由matlab计算得,Fas=1905000N
!
丹麦RIS公式
Pl——风轮单位扫掠面积上的平均风压,通常取Pl=300N/m2
As——风轮的扫掠面积,单位m2
由matlab计算得,Fas=2356200N
2)欧美国家塔架静态强度设计的一般载荷条件
风载条件:
风速65m/s,(2s,平均)风轮停转,叶片顺浆,风向沿机舱横向作用在塔架上。
正常运行工况条件的地震载荷:
考虑额定风速时产生的风轮最大轴向力,同时根据均匀建筑物由地震产生的水平载荷因子,将其产生的惯性力附加在风轮轴向推力上。
【
机组的最大运行载荷:
一般为额定风速下正常运行载荷的2倍。
3)确定塔架设计载荷的要求
设计载荷需要正确的分析塔架承受的载荷及其作用,大致可分为三种类型设计载荷:
最大极限载荷:
指塔架可能承受的最大载荷。
疲劳载荷:
指塔架结构能够承受交变载荷次数的能力。
共振激励载荷:
塔架结构系统的共振响应。
7、联轴器载荷计算
联轴器所承受的载荷以转矩为主,所以可以利用公式来计算。
低速轴联轴器载荷:
=2250630N·
高速轴联轴器载荷:
18454N·
气动特性
攻角
5°
升力系数
阻力系数
发电机
3402kW
额定转速
1650r/min
齿轮箱
齿轮箱传动比
低速轴转速
高速轴转速
齿轮箱功率
3244kW
联轴器
低速轴功率
低速轴载荷
2250630N·
高速轴功率
高速轴载荷
变流器
变流器功率
1500kW
偏航电机
偏航电机功率
4kW
叶片载荷
作用在叶片上的离心力
00N
作用在叶片上的风压力
755580N
作用在叶片上的气动力矩
144250N·
作用在叶片上的陀螺力矩
00N·
风轮载荷
718370N
作用在整个风轮上的转矩
4613600N·
主轴
低速轴转矩
低速轴直径
高速轴转矩
高速轴直径
塔架载荷
作用在塔架上的载荷
前苏联法捷耶夫公式
1108400N
荷兰ECN的公式
1473800N
德国DFVLR公式
1905000N
丹麦RIS公式
2356200N
三、塔架根部截面应力计算
塔架通常采用等强度变截面的设计,危险截面一般位于塔架的根部。
塔筒根部的结构强度分析是确定塔架整体结构尺寸的基本设计依据。
下面考虑塔架高度折减系数的强度计算:
塔架根部的截面应力可以表示为:
W2——塔架根部抗弯截面模数,单位是cm2
A2——塔根部截面积,单位cm2
G2——塔架本身所受重力,单位为N
ψ2——变截面塔架的长度折减系数,这里取
μ——与塔架截面变化有关的折算长度修正系数,可以根据之比确定
Jmin——塔架顶部截面惯性矩,单位cm4
Jmax——塔架根部截面惯性矩,单位cm4
H——塔架高度,单位为cm
γ2——塔架根部截面的惯性半径,单位是cm
Fas取1905000N
h1=1/20*R=
H=90m=9000cm
塔架上部直径取D1=40m,下部直径取D2=50m,壁厚取45/20mm
其中D=50m,α=()/=
代入数值,计算的W2=
A2=
G1+G2为塔顶质量与塔架本身重量之和,取400t
Jmin=
Jmax=
因此,=,所以μ取值为
即,
塔架的根部截面应力为:
40Mpa
四、课程设计结论
所设计的风电机组的技术参数见下表:
:
机组
运行温度
-15°
~45°
叶轮
叶轮直径
叶片长度
叶片数量
结构形式
2级行星+1级平行轴
传动比
型式
双馈异步感应电机
额定输出电压
690V
频率
50Hz
功率因数
容性~感性
额定转速/范围
1650/(1000~2000)rpm
变浆系统
驱动控制
电机
偏航系统
主动式对风式
控制方式
电机+减速机
控制系统
远程监控
塔筒
钢制锥形塔筒
塔筒高
五、参考文献
【1】姚兴佳、田德,宋俊、莴晓明,《风力发电机组设计与制造》,华北电力大学,校内试用教材
【2】徐大平、柳亦兵、吕跃刚,《风力发电原理》,机械工业出版社
【3】朱永强、张旭,《风电场电气系统》,机械工业出版社
【4】敬照亮,《MATLAB教程与应用》,清华大学出版社·
北京交通大学出版社
六、附录
叶片设计和各部分载荷计算的MATLAB程序:
D=ceil(sqrt(8**1000000/pi/(13^3)));
n=14;
lamda=5;
a=[:
:
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/(h+1);
R=D/2;
r=a*R;
omega=lamda*13/R;
yl=r.^2.*C;
Fc=*1000*omega^2*trapz(r,yl)
cl=;
cd=;
y2=(1+(1./tan(I)).^2).*(cl*cos(I)+cd*sin(I)).*C;
Fv=1/2**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**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;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 风力发电机组设计与制造 风力 发电 机组 设计 制造 课程设计 任务书
![提示](https://static.bdocx.com/images/bang_tan.gif)