matlab与优化设计.docx

机械优化设计课程设计题目：齿轮减速器最优化设计班级：机械班成员 2013年6月19日设计题目：二级斜齿圆柱减速器的最优化设计设计要求：要求减速器有最小的体积和最紧凑的结构原始数据：冲头压力（KN）生产效率Piece/min冲头行程（mm）原动机传动方式1104090三相交流异步电动机二级斜齿圆柱减速器设计内容设计方案的拟定及说明电动机的选择及参数计算带轮的初选与计算计算圆柱斜齿轮的输入转矩、传动比、转速，然后建立数学模型编写matlab语言程序，运行程序包括geardesign. m 齿轮系统设计主程序Gearobjfun. m目标函数子程序Gearconstr. m 约束条件子程序Gearparameter.m许用应力计算子程序输出结果该减速器为二级斜齿圆柱减速器，低速级采用二级斜齿圆柱齿轮传动，选择三相交流异步电动机，v带传动确定电动机的容量:选择电动机的容量应保证电动机的额定功率大于等于工作机所需要的功率电动机参数 t=60/40=1.5s v=s/t=6.0*10^-3m/sP=0.5FV=0.5*110*60*0.001=3.3kwη=η1* η32*η23*η4*η5 其中齿轮传动η1=0.96滚动轴承η2=0.98齿轮传动η3=0.97联轴器η4=0.99卷筒η5=1.0 η=0.96*0.98^3*0.97^2*0.99*1.00=0.84Pd=3.3\0.84=2.93kw三，确定电动机的转速已知压片机的转速40piece/min带传动的传动比i1=2~4二级齿轮减速器的传动比i2=8~40，所以电动机的转动范围n=i1i2n=640~6400r/min 可行方案如下确定电动机的转速 具体数据如下型号功率kw满载转速r/min额定转矩最大转矩Y100L2-4314202.22.2计算减速器输入转矩T1，输入转速n，总传动比iiα=1420/40=35.5 初值i0=3计算i=iα/i0=35.5/3=10.52转速n=nm/i0=1420/3=473.33r/min输入转矩Td=9550*3.93/1420=26.43N.m T=Td*i0*η=26.43*3*0.96=72.12N.m轴名效率P（KW）转矩（N.M）转速（n.min）传动比效率输入输出电动机轴2.512696输入轴2.933.3,76.12473.3311.83输入转矩T1=76.12N.m，输入转速n1=473.33r/min,总传动比11.83建立二级减速器的数学模型，并用MATLAB语言编制求解程序，如下程序如下：(1)%geardesign.mglobal rat;globalkk;global tt1;global n0;global t;globalhbs;globalsigmahp;globalsigmafp;kk=1.2;tt1=7.612e4;rat01=4;rat=11.83;n0=473.33;t=24000;hbs=[280,240,280,240];[sigmahp,sigmafp]=gearparameter(rat01);x0=[2,3,30,25,4,10,15,];options=optimset(Display,iter,TolFun,1e-8);lbnd=[2,2,20,20,3,8,8];ubnd=[4,4,40,40,5,16,16];[x,f,exitflag,output]=fmincon(gearobjfun,x0,[],[],[],[],lbnd,ubnd,gearconstr,options)(2)%gearobjfun.mfunction f=gearobjfun(x)global rat;f=(x(1)*x(3)*(1+x(5)))/cos(x(6)*pi/180)+(x(2)*x(4)*(1+rat/x(5)))/cos(x(7)*pi/180);(3)%约束条件子程序gearconstr.mfunction [c,ceq]=gearconstr(x)global rat;global tt1;globalkk;globalsigmahp;globalsigmafp;beta(1)=x(6)*pi/180;beta(2)=x(7)*pi/180;epa(1)=(1.88-3.2*(x(5)+1)/x(3)/x(5))*cos(beta(1));epa(2)=(1.88-3.2*(rat+x(5))/x(4)/rat)*cos(beta(2));fori=1:2yep(i)=1.053/epa(i);zep(i)=sqrt(yep(i));zh(i)=2.9027-0.3702*exp(beta(i));endz