matlab与优化设计

《matlab与优化设计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、机械优化设计课程设计题目：齿轮减速器最优化设计班级：机械班成员2013年6月19日一．设计题目：二级斜齿圆柱减速器的最优化设计二．设计要求：要求减速器有最小的体积和最紧凑的结构三．原始数据：冲头压力（KN）生产效率Piece/min冲头行程（mm）原动机传动方式1104090三相交流异步电动机二级斜齿圆柱减速器四．设计内容1.设计方案的拟定及说明2.电动机的选择及参数计算3.带轮的初选与计算4.计算圆柱斜齿轮的输入转矩、传动比、转速，然后建立数学模型编写matlab语言程序，运行程序包括gearde

2、sign.m齿轮系统设计主程序Gearobjfun.m目标函数子程序Gearconstr.m约束条件子程序Gearparameter.m许用应力计算子程序1.输出结果1.该减速器为二级斜齿圆柱减速器，低速级采用二级斜齿圆柱齿轮传动，选择三相交流异步电动机，v带传动2.确定电动机的容量:选择电动机的容量应保证电动机的额定功率大于等于工作机所需要的功率电动机参数t=60/40=1.5sv=s/t=6.0*10^-3m/sP=0.5FV=0.5*110*60*0.001=3.3kw1.η=η1*η32*η

3、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.84=2.93kw三，确定电动机的转速已知压片机的转速40piece/min带传动的传动比i1=2~4二级齿轮减速器的传动比i2=8~40，所以电动机的转动范围n=i1i2n=640~6400r/min可行方案如下确定电动机的转速具体数据如下型号功率kw满载转速r/min额定转矩最大转矩Y10

4、0L2-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.mT=Td*i0*η=26.43*3*0.96=72.12N.m轴名效率P（KW）转矩（N.M）转速（n.min）传动比效率输入输出电动机轴2.5126.43142030.96输入轴2.933.3,76.12473.3311.

5、83输入转矩T1=76.12N.m，输入转速n1=473.33r/min,总传动比11.83建立二级减速器的数学模型，并用MATLAB语言编制求解程序，如下1.程序如下：(1)%geardesign.m2.globalrat;3.globalkk;4.globaltt1;5.globaln0;6.globalt;7.globalhbs;8.globalsigmahp;9.globalsigmafp;10.kk=1.2;11.tt1=7.612e4;12.rat01=4;13.rat=11.83;14

6、.n0=473.33;15.t=24000;16.hbs=[280,240,280,240];17.[sigmahp,sigmafp]=gearparameter(rat01);18.x0=[2,3,30,25,4,10,15,];19.options=optimset('Display','iter','TolFun',1e-8);20.lbnd=[2,2,20,20,3,8,8];21.ubnd=[4,4,40,40,5,16,16];22.[x,f,exitflag,output]=fminc

7、on('gearobjfun',x0,[],[],[],[],lbnd,ubnd,'gearconstr',options)(2)%gearobjfun.mfunctionf=gearobjfun(x)globalrat;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)globalr

8、at;globaltt1;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