matlab與優(yōu)化設(shè)計(jì).docx

機(jī)械優(yōu)化設(shè)計(jì)課程設(shè)計(jì)題目：齒輪減速器最優(yōu)化設(shè)計(jì)班級(jí)：機(jī)械 班成員 2013年6月19日一 設(shè)計(jì)題目：二級(jí)斜齒圓柱減速器的最優(yōu)化設(shè)計(jì)二 設(shè)計(jì)要求：要求減速器有最小的體積和最緊湊的結(jié)構(gòu)三 原始數(shù)據(jù)：沖頭壓力（KN）生產(chǎn)效率Piece/min沖頭行程（mm）原動(dòng)機(jī)傳動(dòng)方式1104090三相交流異步電動(dòng)機(jī)二級(jí)斜齒圓柱減速器四 設(shè)計(jì)內(nèi)容1. 設(shè)計(jì)方案的擬定及說(shuō)明2. 電動(dòng)機(jī)的選擇及參數(shù)計(jì)算3. 帶輪的初選與計(jì)算4. 計(jì)算圓柱斜齒輪的輸入轉(zhuǎn)矩、傳動(dòng)比、轉(zhuǎn)速，然后建立數(shù)學(xué)模型編寫matlab語(yǔ)言程序，運(yùn)行程序包括geardesign. m 齒輪系統(tǒng)設(shè)計(jì)主程序 Gearobjfun. m目標(biāo)函數(shù)子程序 Gearconstr. m 約束條件子程序 Gearparameter.m許用應(yīng)力計(jì)算子程序5. 輸出結(jié)果1. 該減速器為二級(jí)斜齒圓柱減速器，低速級(jí)采用二級(jí)斜齒圓柱齒輪傳動(dòng)，選擇三相交流異步電動(dòng)機(jī)，v帶傳動(dòng)2. 確定電動(dòng)機(jī)的容量:選擇電動(dòng)機(jī)的容量應(yīng)保證電動(dòng)機(jī)的額定功率大于等于工作機(jī)所需要的功率電動(dòng)機(jī)參數(shù) t=60/40=1.5s v=s/t=6.0*10-3m/s P=0.5FV=0.5*110*60*0.001=3.3kw1. =1* 32*23*4*5 其中齒輪傳動(dòng)1=0.96滾動(dòng)軸承2=0.98齒輪傳動(dòng)3=0.97聯(lián)軸器4=0.99卷筒5=1.0 =0.96*0.983*0.972*0.99*1.00=0.84 Pd=3.30.84=2.93kw三，確定電動(dòng)機(jī)的轉(zhuǎn)速已知壓片機(jī)的轉(zhuǎn)速40piece/min帶傳動(dòng)的傳動(dòng)比i1=24二級(jí)齒輪減速器的傳動(dòng)比i2=840，所以電動(dòng)機(jī)的轉(zhuǎn)動(dòng)范圍n=i1i2n=6406400r/min 可行方案如下確定電動(dòng)機(jī)的轉(zhuǎn)速 具體數(shù)據(jù)如下型號(hào)功率kw滿載轉(zhuǎn)速r/min額定轉(zhuǎn)矩最大轉(zhuǎn)矩Y100L2-4314202.22.2計(jì)算減速器輸入轉(zhuǎn)矩T1，輸入轉(zhuǎn)速n，總傳動(dòng)比ii=1420/40=35.5 初值i0=3計(jì)算i=i/i0=35.5/3=10.52轉(zhuǎn)速n=nm/i0=1420/3=473.33r/min輸入轉(zhuǎn)矩Td=9550*3.93/1420=26.43N.m T=Td*i0*=26.43*3*0.96=72.12N.m軸名效率P（KW）轉(zhuǎn)矩（N.M）轉(zhuǎn)速（n.min）傳動(dòng)比效率輸入輸出電動(dòng)機(jī)軸2.5126.43142030.96輸入軸2.933.3,76.12473.3311.83輸入轉(zhuǎn)矩T1=76.12N.m，輸入轉(zhuǎn)速n1=473.33r/min,總傳動(dòng)比11.83建立二級(jí)減速器的數(shù)學(xué)模型，并用MATLAB語(yǔ)言編制求解程序，如下3. 程序如下：(1)%geardesign.m4. global rat;5. global kk;6. global tt1;7. global n0;8. global t;9. global hbs;10. global sigmahp;11. global sigmafp;12. kk=1.2;13. tt1=7.612e4;14. rat01=4;15. rat=11.83;16. n0=473.33;17. t=24000;18. hbs=280,240,280,240;19. sigmahp,sigmafp=gearparameter(rat01);20. x0=2,3,30,25,4,10,15,;21. options=optimset(Display,iter,TolFun,1e-8);22. lbnd=2,2,20,20,3,8,8;23. ubnd=4,4,40,40,5,16,16;24. 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;global kk;global sigmahp;global sigmafp;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);for i=1:2 yep(i)=1.053/epa(i); zep(i)=sqrt(yep(i); zh(i)=2.9027-0.3702*exp(beta(i);endzv(1)=x(3)/(cos(beta(1)3;zv(2)=x(5)*zv(1);zv(3)=x(4)/(cos(beta(2)3;zv(4)=rat/x(5)*zv(3);for i=1:4 yfa(i)=3.8193-0.3628*log(zv(i);%齒形系數(shù) ysa(i)=1.1233+0.146*log(zv(i);%應(yīng)力校正系數(shù)endfor i=1:2 epsbeta(i)=0.318*x(i+2)*tan(beta(i); ybeta(i)=1.0-epsbeta(i)*x(i+5)/120; if(ybeta(i)0.75)ybeta(i)=0.75; endendc(1)=zh(1)*zep(1)*(cos(beta(1)2*sqrt(1+1/x(5)/(x(1)*x(3)3)*sqrt(kk*tt1)-3.726e-3*sigmahp(2);c(2)=zh(2)*zep(2)*(cos(beta(2)2*sqrt(x(5)*(1+x(5)/rat)/(x(2)*x(4)3)*sqrt(kk*tt1)-3.802e-3*sigmahp(4);c(3)=(cos(beta(1)2/x(3)2/x(1)3*yfa(1)*ysa(1)*yep(1)*ybeta(1)*kk*tt1-0.5*sigmafp(1);c(4)=(cos(beta(1)2/x(3)2/x(1)3*yfa(2)*ysa(2)*yep(1)*ybeta(1)*kk*tt1-0.5*sigmafp(2);c(5)=x(5)*(cos(beta(2)2/x(4)2/x(2)3*yfa(3)*ysa(3)*yep(2)*ybeta(2)*kk*tt1-0.5208*sigmafp(3);c(6)=x(5)*(cos(beta(2)2/x(4)2/x(2)3*yfa(4)*ysa(4)*yep(2)*ybeta(2)*kk*tt1-0.5208*sigmafp(4);%非線性約束問(wèn)題d1=x(1)*x(3)/cos(beta(1);d2=x(1)*x(3)*x(5)/cos(beta(1);d3=x(2)*x(4)/cos(beta(2);d4=x(2)*x(4)*rat*x(5)/cos(beta(2);a1=0.5*x(1)*x(3)*(1+x(5)/cos(beta(1);a2=0.5*x(2)*x(4)*(1+rat/x(5)/cos(beta(2);c(7)=d4-d2-20;c(8)=40+0.5*d2-a2;c(9)=d2-350;c(10)=d4-350;c(11)=50-d1;c(12)=abs(2*x(1)-round(2*x(1)-0.001;c(13)=abs(2*x(2)-round(2*x(2)-0.001;ceq=;(4) %gearparameter.mfunction sigmahp,sigmafp=gearparameter(rat01)global n0;global rat;global kk;global t;global hbs;for i=1:4 sigmah(i)=495+0.889*(hbs(i)-150); sigmaf(i)=170+0.3*(hbs(i)-150);endnr(1)=n0;nr(2)=nr(1)/rat01;nr(3)=nr(1)/rat01;nr(4)=n0/rat;for i=1:4 nn(i)=60*nr(i)*t; if(nn(i)1e+5) zn(i)=1.6; elseif(nn(i)1e+7) zn(i)=2.75-0.1006*log(nn(i); else zn(i)=1.0; end if(nn(i)1e+4) yn(i)=2.5; elseif(nn(i) 輸出最優(yōu)解如下：x = 2.5225 1.2726 40.4671 19.2726 2.2726 7.3134 7.2726f = 490.2449經(jīng)過(guò)適當(dāng)?shù)膱A整和標(biāo)準(zhǔn)化修正后的齒輪傳動(dòng)系統(tǒng)設(shè)計(jì)參數(shù)見表如下：參數(shù) 齒輪齒輪1齒