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

機械優(yōu)化設(shè)計課程設(shè)計題目：齒輪減速器最優(yōu)化設(shè)計班級：機械 班成員 2013年6月19日一 設(shè)計題目：二級斜齒圓柱減速器的最優(yōu)化設(shè)計二 設(shè)計要求：要求減速器有最小的體積和最緊湊的結(jié)構(gòu)三 原始數(shù)據(jù)：沖頭壓力（KN）生產(chǎn)效率Piece/min沖頭行程（mm）原動機傳動方式1104090三相交流異步電動機二級斜齒圓柱減速器四 設(shè)計內(nèi)容1. 設(shè)計方案的擬定及說明2. 電動機的選擇及參數(shù)計算3. 帶輪的初選與計算4. 計算圓柱斜齒輪的輸入轉(zhuǎn)矩、傳動比、轉(zhuǎn)速，然后建立數(shù)學(xué)模型編寫matlab語言程序，運行程序包括geardesign. m 齒輪系統(tǒng)設(shè)計主程序 Gearobjfun. m目標(biāo)函數(shù)子程序 Gearconstr. m 約束條件子程序 Gearparameter.m許用應(yīng)力計算子程序5. 輸出結(jié)果1. 該減速器為二級斜齒圓柱減速器，低速級采用二級斜齒圓柱齒輪傳動，選擇三相交流異步電動機，v帶傳動2. 確定電動機的容量:選擇電動機的容量應(yīng)保證電動機的額定功率大于等于工作機所需要的功率電動機參數(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 其中齒輪傳動1=0.96滾動軸承2=0.98齒輪傳動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三，確定電動機的轉(zhuǎn)速已知壓片機的轉(zhuǎn)速40piece/min帶傳動的傳動比i1=24二級齒輪減速器的傳動比i2=840，所以電動機的轉(zhuǎn)動范圍n=i1i2n=6406400r/min 可行方案如下確定電動機的轉(zhuǎn)速 具體數(shù)據(jù)如下型號功率kw滿載轉(zhuǎn)速r/min額定轉(zhuǎn)矩最大轉(zhuǎn)矩Y100L2-4314202.22.2計算減速器輸入轉(zhuǎn)矩T1，輸入轉(zhuǎn)速n，總傳動比ii=1420/40=35.5 初值i0=3計算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）傳動比效率輸入輸出電動機軸2.5126.43142030.96輸入軸2.933.3,76.12473.3311.83輸入轉(zhuǎn)矩T1=76.12N.m，輸入轉(zhuǎn)速n1=473.33r/min,總傳動比11.83建立二級減速器的數(shù)學(xué)模型，并用MATLAB語言編制求解程序，如下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);%非線性約束問題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)過適當(dāng)?shù)膱A整和標(biāo)準(zhǔn)化修正后的齒輪傳動系統(tǒng)設(shè)計參數(shù)見表如下：參數(shù) 齒輪齒輪1齒