2024年4月1日发(作者:宣如)
第29卷第2期(下)
2013年2月
赤峰学院学报(自然科学版)
Journal of Chifeng University(Natural Science Edition)
VoI_29No.2
Feb.2013
风力发电机齿轮箱建模与优化设计
乔印虎,张春燕,陈杰平,缑瑞宾
(安徽科技学院
摘
机电与车辆工程学院,安徽凤阳233100)
要:阐述了风力发电齿轮箱的构造以及对齿轮箱进行三维建模和优化的过程.其中主要设计过程包括齿轮箱的结构
和工作过程的分析;齿轮箱内各齿轮参数、基本尺寸和强度等的计算.尤其是对齿轮箱内部结构的三维实体建模和齿轮箱优
化的方法与过程.最终使齿轮箱在满足各种强度的条件下使其体积最小,质量最轻以达到提高质量降低生产成本的目的.
关键词:齿轮箱;优化;最小;提高质量
中图分类号:TM315 文献标识码:A 文章编号:1673—260X(2013)02—0011—04
能源是人类社会存在与发展的物质基础,人们在物质
生活和精神生活不断提高的同时,能源危机正向人类袭来.
而风能作为一种清洁能源已越来越受到全世界各国人民的
欢迎和重视IlJ.
来得到想要的结果,这种方法可以省去大量的人力和时间
而且准确.目前优化设计软件已得到了广泛的应用[51.
2齿轮箱计算与三维建模
风力发电机齿轮箱结构简图如图1(a)所示,从图中我们
风力发电机组中的齿轮箱是一个重要的机械部件,其
可以看到,该风力机齿轮箱为一级行星轮加两级平行轴的
混合结构.低速级采用了行星架浮动式均载机构.采用该机
主要功用是将风轮在风力作用下所产生的动力传递给发电
机并使其得到相应的转速.目前我国使用的国内外风电齿轮
箱,主要配套的有GE、维德等公司齿轮箱,以及在此基础上
自行设计的国产风电齿轮箱.而在国外尤其是一些发达国家
在这方面已经取得了很大的成就,通过对齿轮箱的优化设
构可以省去行星轮的支撑,简化了结构,减小了齿轮箱尺
寸,便于在风力机上应用.
啮合鞭率(H
蒋期 第1缀 第1I鳐l箱IlI蹑
1 52 32 196 39 530 25
2 104.64 392.78 1060 50
计不仅能满足其设计要求而且在此基础上使齿轮箱的体积
更小、质量最轻.从而大大降低了设计成本12].
1 齿轮箱工作原理与优化设计理论
旋转j瞎啦(卜I萄
倍频 nz6
'
2
2 9
5 9
31 z5.8 n2 z7
9 8
19 6
25.3
50 5
1.1齿轮箱工作原理
(8)
风力发电机齿轮箱结构简图如图1(a)【3]所示:工作原
理:在叶轮转动下首先带动行星轮系的行星架转动,通过行
星轮与太阳轮的啮合使太阳轮速度提高,然后再通过二级
定轴轮系将速度提高到预定要求,最终将输出转速传递给
发电机 .
1.2齿轮箱优化方法
图1 (a)风力发电机齿轮箱结构简图(b)啮合和旋转频率
为了提高承载能力,齿轮、轴均采用合金钢制造_夕 齿轮
推荐采20CrMnMo、15CrNi6等材料.齿轮箱内用作主传动的
齿轮精度:外齿轮不低于5级(GB/T10095),内齿轮不低于6
级(GB/TlO095).通常最终的热处理方法是渗碳淬火,齿表面
硬度达到HRC60+/-2,具有良好的抗磨损接触强度,轮齿心
通常优化方法有两种,一种是解析法:就是通过人为的
部则具有相对较低的硬度和较好的韧性,能提高抗弯曲强
度,通常齿轮的最终加工是采用磨齿工艺[61.
2.1基本尺寸的计算
对模型进行分析建模列出符合要求的参数方程及一系列约
束方程应用不同的优化方法对其进行分析计算最终得到优
化结果.如:惩罚函数法、基因遗传算法等;另一种就是优化
方法与有限元分析和计算机辅助设计的组合,即通过某种
如图10o)所示:根据该齿轮箱的啮合频率fm和旋转频
率‘的已知数据:
对于行星轮系tm=Zr(n ±n ̄)/60 fr=n/60
计算机软件对已绘制好的实体零件可以直接对其分析计算
基金项目:安徽省教育厅自然科学基金项目(KJ2010B051)
行星轮系实体传动简图如图5所示:
要求以重量最轻、体积最小为目标,对其进行优化设计.
目标函数和设计变量的确定
该行星轮系的重量可取太阳轮和C个行星轮重量之和
来代替,因此目标函数可简化为
f(x)=0.19635m ̄z.2b[4(t*一2)2c]
式中z,一中心轮1的齿数;m一模数,单位为(mm);b—齿
宽,单位为(mm);c一行星轮的个数; 一轮系的传动比.
影响目标函数的独立参数z 、b、m、c应列为设计变量,
且Ⅱ
x=Ix1 X2 X3 x4]T_[z b m c】T
在通常情况下,行星轮个数可以根据机构类型事先选
定,这样,设计变量为
X=Ix1 X2 xdT=[z b m】T
(2)约束条件的建立
1)保证小齿z。不根切,得gl(x)=17一X ≤O
2)限制齿宽最小值,得g2(x)=lO—x ≤0
3)限制模数最小值,得g3(x)=2一X ≤0
4)限制齿宽系数b/m的范围:5≤b/m ̄<17,得
g4(x)=5x3一x2≤0 g5(x)=x2—17x3 ̄0
5)满足接触强度要求,得
既(x)=750937.3/(x1x2、/ )一【 ]H
式中f盯】H__一许用接触应力.
6)满足弯曲强度要求,得
g7(x)=1482000yryJ(x1x2x32)一[盯】F
式中Yv,Y广一齿轮的齿形系数和应力校正系数;
【cr] —许用弯曲应力.
由上面目标函数和约束条件可以用惩罚函数法来对其
优化则方程可写成如下所示
minffx)=O.19635mZza ̄b[4+(p ̄-2)2c]
s.t.gj(x) ̄O(j=l,2,…,6)求该方程的最优解
解构造内点惩罚函数
(x,r)=0.19635m22 [4+( 一2)z—r∑In[一舀(x)]对于任意给
J=1
定的惩罚因子(r(r>0),函数 (x,r)为凸函数.
用解析法求函数 (x,r)的极小值,即令V p(x,r)=0得方程
普=28XIXzX32- _0 (1)
善=14x1 ̄32一 一 1一 =0 (2)
普 dX 一 一X3-z ̄ 1X,一3X' l一 / Xa -X9=0 (3)
联立以上方程组求得:
XI(r)=10±、/lO0+r
,、
158±N/24964+16r
x2(r J=———————— ——————一
(r): 盈
,
当x1(r)=10一vT667时不满足约束条件gllo(x)=17一x ≤0
应舍去
当 (r): 时不满足约束条件 (x):2-x3≤0
应舍去
当X2(r)= 坠 时不满足约束条件勖(X)
=
5x 一X240应舍去
则剩下的便为该点的无约束极值点x (r)、X2*(r)、X3*(r)
经过7次迭代得到最优解x (r)=【21.25 52.94 5.1】T
将最优方案圆整得x (r)=[22 54 8]T
目标函数值f(x (r))=2.03×10
初始点原设计方案x (r)=【21 60 8]T
目标函数值f(x (r))=2.34×10
从结果可以看出,目标函数值下降了15.2%且各项约束
都得到满足.
经过7次迭代得到最优解x (r)=[21.25 52.94 5.1r
将最优方案圆整得x (r)=[22 54 8丁r
目标函数值f(x (r))=2.03×10
初始点原设计方案x (r)=[21 60 8/r
目标函数值f(x (r))=2.34×10
从结果可以看出,目标函数值下降了15.2%且各项约束
都得到满足.
3.2对定轴轮系进行优化
定轴轮系实体图如图6所示,它是两对齿轮啮合所组
成二级平行定轴轮系传动构成的增速装置,从左到右大齿
轮为z ,绿色小齿轮为z ,大齿轮为z。,黄色齿轮为z,.
已知风力机组给定传递的功率p、总传动比i和输出的
转速n.要求在满足强度的条件下,使其体积最小,以达到使
结构紧凑、质量最小的目的.
上面提到,设计时要使体积最小,这就是本优化问题追
求的目标函数.它可以归结为使其总的中心距A为最小,写
成
A 【mnt 1+il) mr 8( i2)卜 mi“
保证总中心距为最小时应满足的条件时本优化设计问
一
13—
题的约束条件,它们是:齿面的接触强度和齿根的弯曲强度
以及中间轴上的小齿轮Z 不与高速轴发生干涉.
图6定轴轮系结构图
(1)齿回援触强度计算给出
[o"a
万]2 ̄ am.13z63il一。s p≥O和
一
cos 0
式中 l1卜一许用接触应力;crH=799.33MP,T 一低速轴的
扭矩;Tr_974000 ̄p
=974000 =643735kg.mm,rr厂中间
nn l,斗
轴的扭矩;T2=974000
P-=974000 =190493kg
.
mm,K 、
n
K厂一载荷系数;取K1=2.76、K2=2.07,中(厂齿宽系数.取 d=0.8
(2)齿根弯曲强度计算给出
高速级大、小齿轮的齿根弯曲强度条件为
4 (1+i )m -3z62一c。s p≥0和
(1舢 cos 13 ̄>o
低速级大、小齿轮的齿根弯曲强度条件为
(1+i2) 3zR2_c0s p≥0
式中 w】一是齿轮的许用弯曲应力;【g ̄=484.72MP,由
机械设计手册P
得:Y6=1.93、Y5=2.27、Y8=1.85、Y7=2.09
(3)根据不干涉条件
_(m 器+s)
式中s—低速轴的轴线和中间轴上的小齿轮齿顶间的
距离.可取¥=5mm.则得mn2z8(1+i2)一2cos13(5+m.2)一in 1z6il >t0.
由上面目标函数和约束条件可以用惩罚函数法来对其
优化刚方程可写成如下所示!
一
14一
minfx) 者[mnlz ( i )栅n2zs( i2)
s.t.&(x)≥0 0=1,2,…,6)
经7次迭代得
x= l X2 x3 X4 X5 X6 x7J
=
[10。46’4.23 70.5 0.26 5.37 60]"r
将最优解圆整得 =【l1 5 70 0-3 6 60]T,目标函数
值f(x3=554.76
原设计方案x =[10 8 67 0.29 6 54]r,目标函数
值f(x3:571.35
目标函数值下降3%且各项约束都得到满足.
3-3整体方案论证
通过以上对齿轮箱的优化可以得出,优化过后的齿轮
箱既能满足齿轮箱的传动要求又能满足齿轮间强度要求.应
选择优化过后的齿轮箱参数来对齿轮箱进行设计计算.即对
于行星轮系按x (r)=[22 54 8r中的参数来设计;对于定
轴轮系按x =[11 5 70 0.3 6 60] 中的参数来设计.
4总结
本文主要针对风力发电机齿轮箱的基本尺寸计算、强
度计算、三位实体建模以及对齿轮箱的优化设计过程,齿轮
箱尺寸和强度计算并绘出齿轮箱的三维实体模型,通过解
析法对其结构进行优化.
参考文献:
[1]汤克平.风电增速箱结构设计叙谈卟机械传动,2004,05.
[2]关立山.世界风力发电现状及展望U].全球科技经济展
望.2004.
[3]唐新安,谢志明,王哲,等.风力机齿轮箱故障诊断.噪声
与振动制,2007,01(2):120—124.
[4]施鹏飞.从世界发展趋势展望我国风力发电前景 中国
电力.2003.
[5]张国瑞,张展.行星传动技术IM].上海交通大学出版社,
1989.
[6]王承勋,张源.风力发电[M].中国电力出版社,2003.
[7]毛谦德,李振清.袖珍机械设计师手册[M].北京:机械工业
出版社.1996.
[8]饶镇纲.行星传动机构设计【M】.国防工业出版社,1994.
2024年4月1日发(作者:宣如)
第29卷第2期(下)
2013年2月
赤峰学院学报(自然科学版)
Journal of Chifeng University(Natural Science Edition)
VoI_29No.2
Feb.2013
风力发电机齿轮箱建模与优化设计
乔印虎,张春燕,陈杰平,缑瑞宾
(安徽科技学院
摘
机电与车辆工程学院,安徽凤阳233100)
要:阐述了风力发电齿轮箱的构造以及对齿轮箱进行三维建模和优化的过程.其中主要设计过程包括齿轮箱的结构
和工作过程的分析;齿轮箱内各齿轮参数、基本尺寸和强度等的计算.尤其是对齿轮箱内部结构的三维实体建模和齿轮箱优
化的方法与过程.最终使齿轮箱在满足各种强度的条件下使其体积最小,质量最轻以达到提高质量降低生产成本的目的.
关键词:齿轮箱;优化;最小;提高质量
中图分类号:TM315 文献标识码:A 文章编号:1673—260X(2013)02—0011—04
能源是人类社会存在与发展的物质基础,人们在物质
生活和精神生活不断提高的同时,能源危机正向人类袭来.
而风能作为一种清洁能源已越来越受到全世界各国人民的
欢迎和重视IlJ.
来得到想要的结果,这种方法可以省去大量的人力和时间
而且准确.目前优化设计软件已得到了广泛的应用[51.
2齿轮箱计算与三维建模
风力发电机齿轮箱结构简图如图1(a)所示,从图中我们
风力发电机组中的齿轮箱是一个重要的机械部件,其
可以看到,该风力机齿轮箱为一级行星轮加两级平行轴的
混合结构.低速级采用了行星架浮动式均载机构.采用该机
主要功用是将风轮在风力作用下所产生的动力传递给发电
机并使其得到相应的转速.目前我国使用的国内外风电齿轮
箱,主要配套的有GE、维德等公司齿轮箱,以及在此基础上
自行设计的国产风电齿轮箱.而在国外尤其是一些发达国家
在这方面已经取得了很大的成就,通过对齿轮箱的优化设
构可以省去行星轮的支撑,简化了结构,减小了齿轮箱尺
寸,便于在风力机上应用.
啮合鞭率(H
蒋期 第1缀 第1I鳐l箱IlI蹑
1 52 32 196 39 530 25
2 104.64 392.78 1060 50
计不仅能满足其设计要求而且在此基础上使齿轮箱的体积
更小、质量最轻.从而大大降低了设计成本12].
1 齿轮箱工作原理与优化设计理论
旋转j瞎啦(卜I萄
倍频 nz6
'
2
2 9
5 9
31 z5.8 n2 z7
9 8
19 6
25.3
50 5
1.1齿轮箱工作原理
(8)
风力发电机齿轮箱结构简图如图1(a)【3]所示:工作原
理:在叶轮转动下首先带动行星轮系的行星架转动,通过行
星轮与太阳轮的啮合使太阳轮速度提高,然后再通过二级
定轴轮系将速度提高到预定要求,最终将输出转速传递给
发电机 .
1.2齿轮箱优化方法
图1 (a)风力发电机齿轮箱结构简图(b)啮合和旋转频率
为了提高承载能力,齿轮、轴均采用合金钢制造_夕 齿轮
推荐采20CrMnMo、15CrNi6等材料.齿轮箱内用作主传动的
齿轮精度:外齿轮不低于5级(GB/T10095),内齿轮不低于6
级(GB/TlO095).通常最终的热处理方法是渗碳淬火,齿表面
硬度达到HRC60+/-2,具有良好的抗磨损接触强度,轮齿心
通常优化方法有两种,一种是解析法:就是通过人为的
部则具有相对较低的硬度和较好的韧性,能提高抗弯曲强
度,通常齿轮的最终加工是采用磨齿工艺[61.
2.1基本尺寸的计算
对模型进行分析建模列出符合要求的参数方程及一系列约
束方程应用不同的优化方法对其进行分析计算最终得到优
化结果.如:惩罚函数法、基因遗传算法等;另一种就是优化
方法与有限元分析和计算机辅助设计的组合,即通过某种
如图10o)所示:根据该齿轮箱的啮合频率fm和旋转频
率‘的已知数据:
对于行星轮系tm=Zr(n ±n ̄)/60 fr=n/60
计算机软件对已绘制好的实体零件可以直接对其分析计算
基金项目:安徽省教育厅自然科学基金项目(KJ2010B051)
行星轮系实体传动简图如图5所示:
要求以重量最轻、体积最小为目标,对其进行优化设计.
目标函数和设计变量的确定
该行星轮系的重量可取太阳轮和C个行星轮重量之和
来代替,因此目标函数可简化为
f(x)=0.19635m ̄z.2b[4(t*一2)2c]
式中z,一中心轮1的齿数;m一模数,单位为(mm);b—齿
宽,单位为(mm);c一行星轮的个数; 一轮系的传动比.
影响目标函数的独立参数z 、b、m、c应列为设计变量,
且Ⅱ
x=Ix1 X2 X3 x4]T_[z b m c】T
在通常情况下,行星轮个数可以根据机构类型事先选
定,这样,设计变量为
X=Ix1 X2 xdT=[z b m】T
(2)约束条件的建立
1)保证小齿z。不根切,得gl(x)=17一X ≤O
2)限制齿宽最小值,得g2(x)=lO—x ≤0
3)限制模数最小值,得g3(x)=2一X ≤0
4)限制齿宽系数b/m的范围:5≤b/m ̄<17,得
g4(x)=5x3一x2≤0 g5(x)=x2—17x3 ̄0
5)满足接触强度要求,得
既(x)=750937.3/(x1x2、/ )一【 ]H
式中f盯】H__一许用接触应力.
6)满足弯曲强度要求,得
g7(x)=1482000yryJ(x1x2x32)一[盯】F
式中Yv,Y广一齿轮的齿形系数和应力校正系数;
【cr] —许用弯曲应力.
由上面目标函数和约束条件可以用惩罚函数法来对其
优化则方程可写成如下所示
minffx)=O.19635mZza ̄b[4+(p ̄-2)2c]
s.t.gj(x) ̄O(j=l,2,…,6)求该方程的最优解
解构造内点惩罚函数
(x,r)=0.19635m22 [4+( 一2)z—r∑In[一舀(x)]对于任意给
J=1
定的惩罚因子(r(r>0),函数 (x,r)为凸函数.
用解析法求函数 (x,r)的极小值,即令V p(x,r)=0得方程
普=28XIXzX32- _0 (1)
善=14x1 ̄32一 一 1一 =0 (2)
普 dX 一 一X3-z ̄ 1X,一3X' l一 / Xa -X9=0 (3)
联立以上方程组求得:
XI(r)=10±、/lO0+r
,、
158±N/24964+16r
x2(r J=———————— ——————一
(r): 盈
,
当x1(r)=10一vT667时不满足约束条件gllo(x)=17一x ≤0
应舍去
当 (r): 时不满足约束条件 (x):2-x3≤0
应舍去
当X2(r)= 坠 时不满足约束条件勖(X)
=
5x 一X240应舍去
则剩下的便为该点的无约束极值点x (r)、X2*(r)、X3*(r)
经过7次迭代得到最优解x (r)=【21.25 52.94 5.1】T
将最优方案圆整得x (r)=[22 54 8]T
目标函数值f(x (r))=2.03×10
初始点原设计方案x (r)=【21 60 8]T
目标函数值f(x (r))=2.34×10
从结果可以看出,目标函数值下降了15.2%且各项约束
都得到满足.
经过7次迭代得到最优解x (r)=[21.25 52.94 5.1r
将最优方案圆整得x (r)=[22 54 8丁r
目标函数值f(x (r))=2.03×10
初始点原设计方案x (r)=[21 60 8/r
目标函数值f(x (r))=2.34×10
从结果可以看出,目标函数值下降了15.2%且各项约束
都得到满足.
3.2对定轴轮系进行优化
定轴轮系实体图如图6所示,它是两对齿轮啮合所组
成二级平行定轴轮系传动构成的增速装置,从左到右大齿
轮为z ,绿色小齿轮为z ,大齿轮为z。,黄色齿轮为z,.
已知风力机组给定传递的功率p、总传动比i和输出的
转速n.要求在满足强度的条件下,使其体积最小,以达到使
结构紧凑、质量最小的目的.
上面提到,设计时要使体积最小,这就是本优化问题追
求的目标函数.它可以归结为使其总的中心距A为最小,写
成
A 【mnt 1+il) mr 8( i2)卜 mi“
保证总中心距为最小时应满足的条件时本优化设计问
一
13—
题的约束条件,它们是:齿面的接触强度和齿根的弯曲强度
以及中间轴上的小齿轮Z 不与高速轴发生干涉.
图6定轴轮系结构图
(1)齿回援触强度计算给出
[o"a
万]2 ̄ am.13z63il一。s p≥O和
一
cos 0
式中 l1卜一许用接触应力;crH=799.33MP,T 一低速轴的
扭矩;Tr_974000 ̄p
=974000 =643735kg.mm,rr厂中间
nn l,斗
轴的扭矩;T2=974000
P-=974000 =190493kg
.
mm,K 、
n
K厂一载荷系数;取K1=2.76、K2=2.07,中(厂齿宽系数.取 d=0.8
(2)齿根弯曲强度计算给出
高速级大、小齿轮的齿根弯曲强度条件为
4 (1+i )m -3z62一c。s p≥0和
(1舢 cos 13 ̄>o
低速级大、小齿轮的齿根弯曲强度条件为
(1+i2) 3zR2_c0s p≥0
式中 w】一是齿轮的许用弯曲应力;【g ̄=484.72MP,由
机械设计手册P
得:Y6=1.93、Y5=2.27、Y8=1.85、Y7=2.09
(3)根据不干涉条件
_(m 器+s)
式中s—低速轴的轴线和中间轴上的小齿轮齿顶间的
距离.可取¥=5mm.则得mn2z8(1+i2)一2cos13(5+m.2)一in 1z6il >t0.
由上面目标函数和约束条件可以用惩罚函数法来对其
优化刚方程可写成如下所示!
一
14一
minfx) 者[mnlz ( i )栅n2zs( i2)
s.t.&(x)≥0 0=1,2,…,6)
经7次迭代得
x= l X2 x3 X4 X5 X6 x7J
=
[10。46’4.23 70.5 0.26 5.37 60]"r
将最优解圆整得 =【l1 5 70 0-3 6 60]T,目标函数
值f(x3=554.76
原设计方案x =[10 8 67 0.29 6 54]r,目标函数
值f(x3:571.35
目标函数值下降3%且各项约束都得到满足.
3-3整体方案论证
通过以上对齿轮箱的优化可以得出,优化过后的齿轮
箱既能满足齿轮箱的传动要求又能满足齿轮间强度要求.应
选择优化过后的齿轮箱参数来对齿轮箱进行设计计算.即对
于行星轮系按x (r)=[22 54 8r中的参数来设计;对于定
轴轮系按x =[11 5 70 0.3 6 60] 中的参数来设计.
4总结
本文主要针对风力发电机齿轮箱的基本尺寸计算、强
度计算、三位实体建模以及对齿轮箱的优化设计过程,齿轮
箱尺寸和强度计算并绘出齿轮箱的三维实体模型,通过解
析法对其结构进行优化.
参考文献:
[1]汤克平.风电增速箱结构设计叙谈卟机械传动,2004,05.
[2]关立山.世界风力发电现状及展望U].全球科技经济展
望.2004.
[3]唐新安,谢志明,王哲,等.风力机齿轮箱故障诊断.噪声
与振动制,2007,01(2):120—124.
[4]施鹏飞.从世界发展趋势展望我国风力发电前景 中国
电力.2003.
[5]张国瑞,张展.行星传动技术IM].上海交通大学出版社,
1989.
[6]王承勋,张源.风力发电[M].中国电力出版社,2003.
[7]毛谦德,李振清.袖珍机械设计师手册[M].北京:机械工业
出版社.1996.
[8]饶镇纲.行星传动机构设计【M】.国防工业出版社,1994.