最后优化改了单级减速器.docVIP

2 优化设计实例 一级圆柱齿轮减速器如图1 所示。已知输入扭矩为T=2056 n.m,输入转速为1000r/m in,传动比i＝5，大齿轮材质45号钢正火，HB1=162~217；小齿轮材质45号钢调质，HB2=217~255，要求在保证齿轮减速器承载能力的条 件下，重量按体积计算最轻。 2.1 建立数学模型 根据设计要求，取设计变量x=[x1,x2,x3,x4,x5,x6] =[b,z1,m ,l,d1,d2]。 建立目标函数f(x)： f(x)=0.785539815(4.75x1x2^2x3^2+85x1x2x3^2 -85x1x3^2+0.92x1x6^2-x1x5^2 +0.8x1x2x3x6-1.6x1x3x6+x4x5^2 +x4x6^2+28x5^2+32x6^2) 建立约束条件: g1(x)=-x2+17≤0 齿数约束 g2(x)=-x1+16x3≤0 间隔限制约束 g3(x)=x1-35x3≤0 间隔限制约束 g4(x)=0.2-x3≤0 模数约束 g5(x)=x2x3-30≤0 间隔限制约束 g6(x)=-x6+13≤0 直径约束 g7(x)=-x5+10≤0 直径约束 g8(x)=x5-15≤0 直径约束 g9(x)=x6-20≤0 直径约束 g10(x)=x1-x4+0.5x6+40≤0 轴规格约束 g11(x)= 383255/x2x3(x1)^(1/2)-8550≤0 接触应力约束 g12(x)=53456/x1x2x3^2(0.169+0.006666x2-0.000854x2^2)-2610≤0 弯曲应力约束 g13(x)= 53456/x1x2x3^2(0.2824+02 -0.0000015765x2)^2-2130≤0 弯曲应力约束 g14(x)=0.009286x4^3/x2x3x5^4-0.003x4≤0 挠度约束 g15(x)=[(21991x4/x2x3)^2+(0.58*20560)^2]^(1/2)/x5^3-55≤0 轴的许用弯曲应力约束 g16(x)=[(21991x4/x2x3)^2+(0.58*5*20560)^2]^(1/2)/x5^3-55≤0 轴的许用弯曲应力约束 2.2 MATLAB 编程计算 fm incon 用到的两个函数定义如下： 目标函数： Function f=G earO bjFunct(x) f=0.785539815(4.75*x(1)*x(2)^2*x(3)^2+85*x(1)*x(2)x(3)^2-85*x(1)*x(3)^2+0.92*x(1)*x(6)^2-x(1)*x(5)^2+ 0.8*x(1)*x(2)*x(3)*x(6)-1.6*x(1)*x(3)*x(6) +x(4)*x(5)^2+x(4)*x(6)^2+28*x(5)^2+32*x(6)^2) 非线性不等式约束函数： function [C,Ceq]=G earN onLinConstr(x) C(1)=-x(2)+17; C(2)=-x(1)+16*x(3); C(3)=x(1)-35*x(3); C(4)=0.2-x(3); C(5)=x(2)*x(3)-30; C(6)=-x(6)+13; C(7)=-x(5)+10; C(8)=x(5)-15; C(9)=x(6)-20; C(10)=x(1)-x(4)+0.5*x(6)+40; C(11)=383255/(x(2)*x(3)*sqrt(x(1)))-8550; C(12)=53456/(x(1)*x(2)*x(3)^2*(0.169+0.006666*x(2)-0.000854*x(2)^2))-2610; C(13)= 53456/(x(1)*x(2)*x(3)^2*(0.2824+0x(2 )-0.0000015765*x(2))^2)-2130; C(14)=0.009286*x(4)^3/(x(2)*x(3)*x(5)^4)-0.003*x(4); C(15)=((21991*x(4)/(x(2)*x(3)))^2+(0.58*20560)^2)^(1/2)/x(5)^3-55; C(16)=((21991*x(4)/(x(2)*x(3)))^2+(0.58*5*20560)^2)^(1/2)/x(5)^3-55; Ceq=[]; 将上面的两个函数分别保存为：G earO bjFunct.m 和 G earN onLinConstr.m 文件。 最后在M A TLA B 命令行窗口调用上面的两个函数： ％ 程序如下： x0=[3.2,17,0.2,15,10,13]; LBnd= [3.2,17,0.2,15,10,13]; UBnd= [7,30,0.3