汽車主減速器傳動比優(yōu)化

汽車主減速器傳動比優(yōu)化汽車主減速器傳動比優(yōu)化某輕型貨車，車有關(guān)數(shù)據(jù)如下：發(fā)動機的最低轉(zhuǎn)速nmin＝600r/min，最高轉(zhuǎn)速nmax=4000r/min?？傎|(zhì)量3650kg；車輪半徑0.367m；傳動系機械效率nT＝0.85；滾動阻力系數(shù)f＝0.013；空氣阻力系數(shù)*迎風(fēng)面積CDA＝2.77m2；飛輪轉(zhuǎn)動慣量If＝0.218kg.m2；二前輪轉(zhuǎn)動慣量Iw1＝1.798kg.m2;軸距L=3.2m;質(zhì)心至前軸距離（滿載）a=1.947m；質(zhì)心高（滿載）hg=0.9m；四后輪轉(zhuǎn)動慣量Iw2＝3.598kg.m2；T=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;怠速油耗Qid＝0.299mL/s（怠速轉(zhuǎn)速400r/min）。負荷特性曲線擬合公式為：b=B0+B1*Pe+B2*Pe.^2+B3*Pe.^3+B4*Pe.^4.N/(r/min)B0B1B2B3B48151326.8-416.4672.379-5.86290.1776812071354.7-303.9836.657-2.05530.04307216141284.4-189.751452.4-0.511840.006816420121122.9-121.597.0035-0.185170.001855526031141.0-98.8934.4763-0.0910770.0006890630061051.2-73.7142.8593-0.051380.0003503234031233.9-84.4782.9788-0.0474490.0002823038041129.7-45.2910.71113-0.00075215-0.0000385681檔2檔3檔4檔5檔變速器傳動比5.064.0163.091.711計算1計算加速時間Matlab程序：n=600:1:4000;m=3650;g=9.8;G=m*g;r=0.367;i0=5.089;nT=0.85;f=0.013;CdA=2.77;Iw1=1.798;Iw2=3.598;Iw=Iw1+Iw2;If=0.218;a=1.947;hg=0.9;L=3.2;Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;%以2擋加速時間計算u(1)=0.377*r*600/i0/4.016;dt=0.01;w=1;t(1)=0;ig=4.016;while(u(w)<99.08)n(w)=u(w)*ig*i0/0.377/r;Ttq(w)=-19.313+295.27*n(w)/1000-165.44*(n(w)/1000)^2+40.874*(n(w)/1000)^3-3.8445*(n(w)/1000)^4;s=[50,250,250];Q1=s.*Pa.*b./(102*ua1*pg);disp(Q1);Q1=Q1(1)+Q1(2)+Q1(3);%勻加速行駛的油耗%第一段ua=25:1:40;a=0.25;q=1+Iw/(m*r^2)+If*ig^2*i0^2*nT/(m*r^2);Pa=(G*f*ua/3600+CdA*ua.^3/76140+q*m*ua.*a/3600)/nT;n=ua*ig*i0/0.377/r;B0=spline(n1,b0,n);B1=spline(n1,b1,n);B2=spline(n1,b2,n);B3=spline(n1,b3,n);B4=spline(n1,b4,n);b=B0+B1.*Pa+B2.*Pa.^2+B3.*Pa.^3+B4.*Pa.^4;b=spline(Pa,b,Pa);Qs=Pa.*b/(367.1*pg);t=1/(3.6*a);Q21=0;forn=1:15Q21=Q21+(Qs(n)+Qs(n+1))*t/2;enddisp(Q21);%第二段ua=40:1:50;a=0.20;q=1+Iw/(m*r^2)+If*ig^2*i0^2*nT/(m*r^2);Pa=(G*f*ua/3600+CdA*ua.^3/76140+q*m*ua.*a/3600)/nT;n=ua*ig*i0/0.377/r;B0=spline(n1,b0,n);B1=spline(n1,b1,n);B2=spline(n1,b2,n);B3=spline(n1,b3,n);B4=spline(n1,b4,n);b=B0+B1.*Pa+B2.*Pa.^2+B3.*Pa.^3+B4.*Pa.^4;b=spline(Pa,b,Pa);Qs=Pa.*b/(367.1*pg);t=1/(3.6*a);Q22=0;forn=1:10Q22=Q22+(Qs(n)+Qs(n+1))*t/2;enddisp(Q22);Q2=Q21+Q22;%勻減速行駛的油耗Qid=0.299;t3=19.3;Q3=Qid*t3;disp(Q3);%總油耗s=1075;Q=(Q1+Q2+Q3)*100/1075;disp(Q);io5.0895.445.856.016.17Q(L)13.000613.408613.739913.853513.9682繪制燃油經(jīng)濟性-加速時間曲線Q=[13.000613.408613.739913.853513.9682];t=[24.055324.239624.055323.871023.6866];plot(Q,t,'+r');holdonq1=linspace(Q(1),Q(5),100);t1=spline(Q,t,q1);%三次樣條插值plot(q1,t1);%繪制燃油經(jīng)濟性-加速時間曲線title('燃油經(jīng)濟性—加速時間曲線');xlabel('百公里油耗/(L/100km)');ylabel('加速時間/s');gtext('i0=5.089'),gtext('i0=5.44'),gtext('i0=5.85'),gtext('i0=6.01'),gtext('i0=6.17');計算最高車速n=600:1:4000;r=0.367;i0=5.089;nT=0.85;f=0.013;m=3650;g=9.8;G=m*g;CdA=2.77;a=1.947;hg=0.9;L=3.2;If=0.218;Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;%驅(qū)動力行駛阻力平衡圖forig=[5.06,4.016,3.09,1.71,1]Ua=0.377*r*n/ig/i0;Ft=Ttq*ig*i0*nT/r;plot(Ua,Ft);holdon;endFf=G*f;ua=0:0.1:max(Ua);Fw=CdA*ua.^2/21.15;plot(ua,(Ff+Fw));title('驅(qū)動力-行駛阻力平衡圖');xlabel('Ua/(km/h)');ylabel('Ft/N');i05.0895.445.856.016.17umax102.4102.1100.299.297.95.繪制i0-umax曲線（求最高車速約束條件下臨界i0值）i0=[5.0895.445.856.016.17];umax=[102.4102.1100.299.297.9];plot(i0,umax,'+r');holdoni01=linspace(i0(1),i0(5));umax1=spline(i0,umax,i01);plot(i01,umax1);title('i0-umax曲線');xlabel('i0');ylabel('最高車速km/h');[x,y]=ginput(1);disp('臨界i0=');disp(x);%根據(jù)最高車速約束條件，求臨界i0值gtext('i0=5.089'),gtext('i0=5.44'),gtext('i0=5.85'),gtext('i0=6.01'),gtext('i0=6.17');臨界i0=5.8855最高車速大于100km/h.所以io取5.0895.445.85繪制i0-t圖（求加速時間約束下臨界i0值）i0=[5.0895.445.856.016.17];t=[24.055324.239624.055323.871025.6866];plot(i0,t,'+r');holdon；i01=linspace(i0(1),i0(5));umax1=spline(i0,t,i01);plot(i01,umax1);title('i0-t曲線');xlabel('i0');ylabel('加速時間s');zoomon;[x,y]=ginput(1);zoomoff;disp('臨界i0=');disp(x);gtext('i0=5.089'),gtext('i0=5.44'),gtext('i0=5.85'),gtext('i0=6.01'),gtext('i0=6.17');加速時間不超過25s所以io可以 取5.0895.445.856.017最大爬坡度計算n=600:1:4000;r=0.367;i0=5.089;eff=0.85;f=0.013;m=3650;g=9.8;G=m*g;CdA=2.77;a=1.947;hg=0.9;L=3.2;If=0.218;Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;ig=[5.064.0163.091.711];Ft1=Ttq*ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw1=CdA*ua1.^2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;imax=100*tan(asin(max(Fi1/G)));disp('汽車最大爬坡度=');disp(imax);disp('%');i05.0895.445.856.016.17i(%)28.973131.274134.013335.101736.2011繪制i0-i曲線（求最大爬坡度約束下i0)i0=[5.0895.445.856.016.17];i=[28.973131.274134.013335.101736.2011];plot(i0,i,'+r');holdoni01=linspace(i0(1),i0(5),100);i1=spline(i0,i,i01);plot(i01,i1);title('i0-imax曲線');xlabel('i0');ylabel('imax(%)');[x,y]=ginput(1);disp('臨界i0=');disp(x);gtext('i0=5.089'),gtext('i0=5.44'),gtext('i0=5.85'),gtext('i0=6.01'),gtext('i0=6.17');最大的爬坡度大于32%，因此io大于5.55,io可取5.856.016.17。8最高檔動力因數(shù)校核一最大轉(zhuǎn)矩計算n=600:1:4000;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;plot(n,Tq);holdon;title('n-Tq曲線');xlabel('轉(zhuǎn)速nr/min');ylabel('轉(zhuǎn)矩Tqn*m');[max_Tq,nn]=max(Tq);disp('最大值:');disp(max_Tq);disp('對應(yīng)的n的值:');disp(n(nn));最大轉(zhuǎn)矩174.9667最大轉(zhuǎn)矩的轉(zhuǎn)速2041二校核最高檔的動力因數(shù)Tmax=174.9667;n=2041;CdA=2.77;nT=0.85;ig=1;r=0.367;m=3650;g=9.8;G=m*g;i0=[5.445.85];F=Tmax*i0*ig*nT/r;uat=0.377