二级斜齿圆柱齿轮减速器优化设计_斜齿轮二级减速器设计
二级斜齿圆柱齿轮减速器优化设计由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“斜齿轮二级减速器设计”。
二级斜齿圆柱齿轮减速机
优化设计
1.题目
二级斜齿圆柱齿轮减速机。高速轴输入功率R=6.2kW,高速轴转速n1=1450r/min,总传动比iΣ=31.5,齿轮的齿宽系数Φa=0.4;齿轮材料和热处理;大齿轮45号钢正火硬度为187~207HBS,小齿轮45号钢调质硬度为228~255HBS。总工作时间不小于10年。要求按照总中心距最小确定总体方案中的主要参数。
2.已知条件
已知高速轴输入功率R=6.2kW,高速轴转速n1=1450r/min,总传动比iΣ=31.5,齿轮的齿宽系数Φa=0.4。
3.建立优化模型
3.1问题分析及设计变量的确定
由已知条件求在满足使用要求的情况下,使减速机的总中心距最小,二级减速机的总中心距为:
aa1a2其中mn1z11i1mn2z31i22cos
mn1、mn2zz分别为高速级和低速级齿轮副的模数,1、3分
ii别为高速级和低速级小齿轮齿数,1、2分别为高速级和低速级传动比,为齿轮副螺旋角。所以与总中心距a相关的独立参数为:mn1、mn2、z1、z3、i1(i231.5)、。则设计变量可取为:
i1x=[mn1 mn2 z1 z3 i1 ]T=[x1 x2 x3 x4 x5 x6]T 3.2目标函数为
fxx1x31x5x2x4131.5x52cosx6
为了减速机能平稳运转,所以必须满足以下条件:
2mn15、3.5mn26、14z12216z322、5.8i17、815
3.3约束条件的建立3.3.1线性不等式约束条件
g1x2x10 g2xx150 g3x3.5x20 g4xx260 g5x14x30 g6xx3220 g7x16x40 g8xx4220 g9x5.8x50 g10xx570 g11x8x60 g12xx6150
3.3.2非线性不等式约束条件
1)齿轮的接触应力不得大于许用应力值,得
H1925i113K1T1925b1i1H1 H2 H2即 i213K2T2b2i2 2 H1mn31z13i189252K1T12cos30
H2m28925KT33n232222zicos302)齿轮的弯曲应力不得大于许用弯曲应力值,得
F11.5K1T1F1bd1mn1Y2Y1F2Y2
F2F1即
F1Y13K1T11i1mn31z12cos20 F2Y23K1T1和 1i1mn31z12cos201i2mn32z32cos20 F3Y33K2T2F4Y43K2T21i2mn32z32cos20其中齿形系数的计算如下:
Y10.1690.006666z10.0000854z122Y20.1690.006666z20.0000854z2Y30.1690.006666z30.0000854z23
2Y40.1690.006666z40.0000854z43)高速级齿轮和低速级齿轮不得发生干涉,得:
mn2z31i22cosEmn1mn1z1i10
E为低速轴轴线与高速级大齿轮齿顶圆之间的距离,单位为mm。
大齿轮45号钢正火硬度为187~207HBS,查表得计算可得小齿轮45号钢调质硬度为228~255HBS H1H2518.75N/mm2,F1F3153.5N/mm2 2F2F4141.6N/mm高速轴输入功率R=6.2kW,高速轴转速n1=1450r/min,计算可得
T141690N/mm、T240440i1N/mm、查表可得:
K10.248、K21.204、Y10.248、Y20.302、Y30.256、Y40.302、E=50mm整理可得:
3g13xcos3x63.079106x13x3x50233g14xx5cos3x61.701104x2x403g15xcos2x69.939105x13x31x50232g16xx5cos2x61.70610431.5x5x2x402g17xx52x50cosx6x1x2x51x2x4x531.502g18xcos2x61.116104x13x31x50232g19xx5cos2x61.17110431.5x5x2x40
4.优化方法的选择
由于该问题有6个设计变量,19个约束条件的优化设计问题,采用传统的优化设计方法比较繁琐,比较复杂,所以选用Matlab优化工具箱中的fmincon函数来求解此非线性优化问题,避免了较为繁重的计算过程。5.数学模型的求解
5.1.1将已知及数据代入上式,该优化设计的数学优化模型表示为:
minfxx1x31x5x2x4131.5x52cosx6
Subject to: 4 g1x2x10g2xx150g3x3.5x20g4xx260g5x14x30g6xx3220g7x16x40g8xx4220g9x5.8x50g10xx570g11x8x60g12xx61503g13xcos3x63.079106x13x3x50233g14xx5cos3x61.701104x2x403g15xcos2x69.939105x13x31x50232g16xx5cos2x61.70610431.5x5x2x402g17xx52x50cosx6x1x2x51x2x4x531.502g18xcos2x61.116104x13x31x50232g19xx5cos2x61.17110431.5x5x2x40
5.1.2运用Matlab优化工具箱对数学模型进行程序求解
首先在Matlab优化工具箱中编写目标函数的M文件 myfun.m,返回x处的函数值f:
function f = myfun(x)f=(x(1)*x(3)*(1+x(5))+x(2)*x(4)*(1+31.5/x(5)))/(2*cos(x(6)))由于约束条件中有非线性约束,故需要编写一个描述非线性约束条件的M文件myobj.m:
function[c,ceq]=myobj(x)c=[2-x(1);x(1)-5;3.5-x(2);x(2)-6;14-x(3);x(3)-22;16-x(4);x(4)-22;5.8-x(5);x(5)-7;8-x(6);x(6)-15;(cos(x(6)))^3-3.079*10^(-6)*x(1)^3*x(3)^3*x(5);x(5)^2*(cos(x(6)))^3-1.701*10^(-4)*x(2)^3*x(4)^3;(cos(x(6)))^2-9.939*10^(-5)*x(1)^3*x(3)^3*(1+x(5));x(5)^2*(cos(x(6)))^2-1.706*10^(-4)*(31.5+x(5))*x(2)^3*x(4)^2;x(5)*(2*(x(1)+50)*(cos(x(6)))^2+x(1)*x(2)*x(5))-x(2)*x(4)*(x(5)+31.5);(cos(x(6)))^2-1.116*10^(-4)*x(1)^3*x(3)^2*(1+x(5));x(5)^2*(cos(x(6)))^2-1.171*10^(-4)*(31.5+x(5))*x(2)^3*x(4)^2];ceq=[];最后在command window里输入: x0=[3;4;20;20;6;10];%给定初始值
[x,fval,exitflag,output]=fmincon(@myfun,x0,[],[],[],[],[],[],@myobj)%调用优化过程
5.1.3最优解以及结果分析 运行结果如下图所示:由图可知,优化后的结果为:
x=[2.0471 3.6059 18.5067 16 5.8 8] f(x)= 3.1742e+03 圆整得:
x=[2 3.5 19 16 5.8 8] f(x)= 3.1750e+03 6.结果对比分析
如按常规设计,即x=[3 5 19 17 6.3 11],则中心距为470,相比中心距减少了32.5%。说明优化结果是成功的。7.总结体会
做优化难点在各个约束的确定,特别是非线性约束,需要对减速机的知识有较深的认识,另外Matlab软件的运用也相当重要,使用过程当中难免出现问题,要能够解决过程中出现的问题。
在将约束写进Matlab command window时要细心,不然会出现很多问题,得不到正确的结果。8.参考文献
[1]濮良贵,纪名刚.机械设计.8版.北京:高等教育出版社,2006.[2] 孙桓,陈作模,葛文杰.机械原理.7版.北京:高等教育出版社,2006.[3] 胡新华.单级圆柱齿轮减速器的优化设计【J】组合机床与自动化加工技术,2006(7)