汽車理論編程作業(yè)1.2.3章

1、精選優(yōu)質(zhì)文檔-傾情為你奉上汽車理論作業(yè)：第一章1.3選四檔變速器，用MATLAB語言驅(qū)動力與行駛阻力平衡圖m1=2000;m2=1800;mz=3880;g=9.81;f=0.013;n=600:1:4000;nt=0.85;r=0.367;i0=5.83;ig1=6.09;ig2=3.09;ig3=1.71;ig4=1.00;ua1=0.377*r.*n/(ig1.*i0);ua2=0.377*r.*n/(ig2.*i0);ua3=0.377*r.*n/(ig3.*i0);ua4=0.377*r.*n/(ig4.*i0);tq=-19.313+295.27*(n/1000)-165.44*(

2、n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;ft1=tq.*ig1.*i0.*nt/r;ft2=tq.*ig2.*i0.*nt/r;ft3=tq.*ig3.*i0.*nt/r;ft4=tq.*ig4.*i0.*nt/r;fz1=mz.*g.*f+2.77*ua1.2/21.15; %fz=Ff+Fwfz2=mz.*g.*f+2.77*ua2.2/21.15;fz3=mz.*g.*f+2.77*ua3.2/21.15;fz4=mz.*g.*f+2.77*ua4.2/21.15;plot(ua1,ft1,ua2,ft2,ua3,ft3,ua4,ft4

3、,ua1,fz1,ua2,fz2,ua3,fz3,ua4,fz4);title(驅(qū)動力與行駛阻力平衡圖);Xlabel(ua(km/h);Ylabel(ft(N);最高車速m1=2000;m2=1800;mz=3880;g=9.81;f=0.013;k=600:1:4000;n(k)=k;nt=0.85;r=0.367;i0=5.83;ig1=6.09;ig2=3.09;ig3=1.71;ig4=1.00;ua1(k)=0.377*r.*n(k)/(ig1.*i0);ua2(k)=0.377*r.*n(k)/(ig2.*i0);ua3(k)=0.377*r.*n(k)/(ig3.*i0);ua

4、4(k)=0.377*r.*n(k)/(ig4.*i0);tq(k)=-19.313+295.27*(n(k)/1000)-165.44*(n(k)/1000).2+40.874*(n(k)/1000).3-3.8445*(n(k)/1000).4;ft1(k)=tq(k).*ig1.*i0.*nt/r;ft2(k)=tq(k).*ig2.*i0.*nt/r;ft3(k)=tq(k).*ig3.*i0.*nt/r;ft4(k)=tq(k).*ig4.*i0.*nt/r;fz1(k)=mz.*g.*f+2.77*ua1(k).2/21.15; %fz=Ff+Fwfz2(k)=mz.*g.*f+2

5、.77*ua2(k).2/21.15;fz3(k)=mz.*g.*f+2.77*ua3(k).2/21.15;fz4(k)=mz.*g.*f+2.77*ua4(k).2/21.15;ua4=max(ua4(k)disp(汽車的最高車速=);disp(ua4);disp(km/h);得到汽車的最高車速= 94.9290(km/h)最大坡度和相應的附著率k=600:1:4000;f=0.013;l=3.2;a=1.947;n(k)=knt=0.85;r=0.367;f=0.013;i0=5.83;m=3880;g=9.81;G=m.*g;CDA=2.77;If=0.218;Iw1=1.798;Iw

6、2=3.598;ig1=6.09;ua1(k)=0.377*r.*n(k)/(ig1.*i0);tq(k)=-19.313+295.27*(n(k)/1000)-165.44*(n(k)/1000).2+40.874*(n(k)/1000).3-3.8445*(n(k)/1000).4;b1=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig1.2.*i0.2.*nt)/(m.*r.2);ft1(k)=tq(k).*ig1.*i0.*nt/r;fw1(k)=CDA.*(ua1(k).2)/21.15;D1(k)=(ft1(k)-fw1(k)./G;i(k)=D1(k)-fi=max(i

7、(k);q=i.*l./a;disp(汽車的最大爬坡度=);disp(i);disp(最大爬坡度相應的附著率=);disp(q);得計算結(jié)果：汽車的最大爬坡度= 0.3648最大爬坡度相應的附著率= 0.5995加速度倒數(shù)曲線n=600:1:4000;nt=0.85;r=0.367;f=0.013;i0=5.83;m=3880;g=9.81;G=m.*g;CDA=2.77;If=0.218;Iw1=1.798;Iw2=3.598;ig1=6.09;ig2=3.09;ig3=1.71;ig4=1.00;ua1=0.377*r.*n/(ig1.*i0);ua2=0.377*r.*n/(ig2.*i

8、0);ua3=0.377*r.*n/(ig3.*i0);ua4=0.377*r.*n/(ig4.*i0);tq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;b1=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig1.2.*i0.2.*nt)/(m.*r.2);b2=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig2.2.*i0.2.*nt)/(m.*r.2);b3=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig3.2.*i0.2.*nt)/(m.*r.2

9、);b4=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig4.2.*i0.2.*nt)/(m.*r.2);ft1=tq.*ig1.*i0.*nt/r;ft2=tq.*ig2.*i0.*nt/r;ft3=tq.*ig3.*i0.*nt/r;ft4=tq.*ig4.*i0.*nt/r;fw1=CDA.*(ua1.2)/21.15;fw2=CDA.*(ua2.2)/21.15;fw3=CDA.*(ua3.2)/21.15;fw4=CDA.*(ua4.2)/21.15;D1=(ft1-fw1)./G;D2=(ft2-fw2)./G;D3=(ft3-fw3)./G;D4=(ft4-fw4)./

10、G;a1=g.*(D1-f)/b1;a2=g.*(D2-f)/b2;a3=g.*(D3-f)/b3;a4=g.*(D4-f)/b4;plot(ua1,1./a1,ua2,1./a2,ua3,1./a3,ua4,1./a4);title(加速度倒數(shù)曲線圖);Xlabel(ua(km/h);Ylabel(1/a(m/s2);加速時間k=600:1:4000;n(k)=k;nt=0.85;r=0.367;f=0.013;i0=5.83;m=3880;g=9.81;G=m.*g;CDA=2.77;If=0.218;Iw1=1.798;Iw2=3.598;ig2=3.09;ig3=1.71;ig4=1.

11、00;ua2(k)=0.377*r.*n(k)/(ig2.*i0);ua3(k)=0.377*r.*n(k)/(ig3.*i0);ua4(k)=0.377*r.*n(k)/(ig4.*i0);tq(k)=-19.313+295.27*(n(k)/1000)-165.44*(n(k)/1000).2+40.874*(n(k)/1000).3-3.8445*(n(k)/1000).4;b2=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig2.2.*i0.2.*nt)/(m.*r.2);b3=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig3.2.*i0.2.*nt)/(m.*r

12、.2);b4=1+(Iw1+Iw2)./(m.*r.2)+(If.*ig4.2.*i0.2.*nt)/(m.*r.2);ft2(k)=tq(k).*ig2.*i0.*nt/r;ft3(k)=tq(k).*ig3.*i0.*nt/r;ft4(k)=tq(k).*ig4.*i0.*nt/r;fw2(k)=CDA.*(ua2(k).2)/21.15;fw3(k)=CDA.*(ua3(k).2)/21.15;fw4(k)=CDA.*(ua4(k).2)/21.15;D2(k)=(ft2(k)-fw2(k)./G;D3(k)=(ft3(k)-fw3(k)./G;D4(k)=(ft4(k)-fw4(k).

13、/G;a2(k)=g.*(D2(k)-f)/b2;a3(k)=g.*(D3(k)-f)/b3;a4(k)=g.*(D4(k)-f)/b4;t2=trapz(ua2(k),a2(k);t3=trapz(ua3(k),a3(k);t4=trapz(ua4(k),a4(k);tz=(t2+t3+t4)/3.6;disp(加速時間為);disp(tz);得加速時間為：27.4156 s第二章2.7題的程序:m=3880;g=9.8;r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;ig=6.09 3.09 1.7

14、1 1;nd=400;Qid=0.299;nn=815 1207 1614 2012 2603 3006 3403 3804;for i=1:69 n(i)=(i+11)*50; Tq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).2+40.874*(n(i)/1000).3-3.8445*(n(i)/1000).4;endfor j=1:4 for i=1:69 ua(i,j)=0.377.*r.*n(i)/(ig(j).*i0); Pe(i,j)=Tq(i).*ig(j).*i0.*nt.*ua(i,j)/(3600.*r); Pfw(i

15、,j)=(m.*g.*f.*ua(i,j)/3600+CdA.*ua(i,j).3/76140)/nt; endendplot(ua,Pe,ua,Pfw)title(功率平衡圖);xlabel(ua(km/h); ylabel(Pe,Pfw(kw);功率平衡圖如圖： 最高檔與次高檔的等速百公里油耗曲線m=3880;g=9.8;r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;ig=6.09 3.09 1.71 1;nd=400;Qid=0.299;nn=815 1207 1614 2012 2603 30

16、06 3403 3804;for j=1:4 for i=1:8 Tq1(i)=-19.313+295.27*(nn(i)/1000)-165.44*(nn(i)/1000).2+40.874*(nn(i)/1000).3-3.8445*(nn(i)/1000).4; ua1(i,j)=0.377*nn(i).*r./(ig(j).*i0); Pe1(i,j)=Tq(i).*ig(j).*i0.*ua1(i,j)/(3600.*r); Pfw(i,j)=(m.*g.*f.*ua1(i,j)/3600+CdA.*ua1(i,j).3/76140)/nt; endendb(1)=0.17768*P

17、e1(1)4-5.8629*Pe1(1)3+72.379*Pe1(1)2-416.46*Pe1(1)+1326.8; b(2)=0.*Pe1(2)4-2.0553*Pe1(2)3+36.657*Pe1(2)2-303.98*Pe1(2)+1354.7; b(3)=0.*Pe1(3)4-0.51184*Pe1(3)3+14.524*Pe1(3)2-189.75*Pe1(3)+1284.4; b(4)=0.*Pe1(4)4-0.18517*Pe1(4)3+7.0035*Pe1(4)2-121.59*Pe1(4)+1122.9; b(5)=0.*Pe1(5)4-0.*Pe1(5)3+4.4763*P

18、e1(5)2-98.893*Pe1(5)+1141.0; b(6)=0.*Pe1(6)4-0.05138*Pe1(6)3+2.8593*Pe1(6)2-73.714*Pe1(6)+1051.2; b(7)=0.*Pe1(7)4-0.*Pe1(7)3+2.9788*Pe1(7)2-84.478*Pe1(7)+1233.9; b(8)=-0.*Pe1(8)4-0.*Pe1(8)3+0.71113*Pe1(8)2-45.291*Pe1(8)+1129.7;v3=ua1(:,3); v4=ua1(:,4);for s=1:8 Qs3(s)=Pfw(s,3)/0.85.*b(s)/(1.02*v3(s)

19、*7);%次高擋，%取密度和g的乘積為7.0kg*m/(L*s2) Qs4(s)=Pfw(s,4)/0.85.*b(s)/(1.02*v4(s)*7);%最高擋endfigureplot(v3,Qs3,-,markersize,14);hold onplot(v4,Qs4,-,markersize,14);hold offtitle(等速百公里油耗曲線);xlabel(ua(km/h); ylabel(Qs(L/100km);四檔三檔第三章m=3880;g=9.8;r=0.367;nt=0.85;f=0.013;CdA=2.77;If=0.218;Iw1=1.798;Iw2=3.598;ig=

20、6.09 3.09 1.71 1;L=3.2;a=1.947;hg=0.9;i0=5.17 5.43 5.83 6.17 6.33;%選取最高擋，轉(zhuǎn)速n=3403r/min,Pe=60kw，對應的油耗比較經(jīng)濟Pe=60;B=0.*Pe4-0.*Pe3+2.9788*Pe2-84.478*Pe+1233.9; for s=1:5ua0(s)=0.377*r*3403/(ig(4)*i0(s); %設各個i0值對應的速度在該情況下為ua0Qs(s)=Pe*B/(1.02*ua0(s)*7.0); %求出各i0值的等速百公里油耗endfor i=1:69 n(i)=(i+11)*50;Tq(i)=-

21、19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).2+40.874*(n(i)/1000).3-3.8445*(n(i)/1000).4;endfor s=1:5 %求各個i0值對應的加速至最高速時間，從1擋起步 for j=1:4 %求加速度倒數(shù)曲線 for i=1:69 deta(j)=1+(Iw1+Iw2) /(m.*r.2)+If.*ig(j).2.*i0(s).2.*nt/(m.*r.2); ua(i,j)=0.377.*r.*n(i)/(ig(j).*i0(s); Ft(i,j)=Tq(i).*ig(j).*i0(s).*nt/r; Fw(i

22、,j)=CdA.*ua(i,j).2/21.15; a(i,j)=(Ft(i,j)-Fw(i,j)-m.*g.*f)/(deta(j)*m); b(i,j)=1/a(i,j); %b(i,j)為加速度倒數(shù)函數(shù) endendk=-1;i=1;while(k0)k=ua(i,2)-ua(69,1);i=i+1;endw=i-1; %w為從1擋換至2擋時的轉(zhuǎn)速的自變量 k=-1;i=1; while (k0)k=ua(i,3)-ua(69,2);i=i+1;endx=i-1; %x為換至3擋時的轉(zhuǎn)速的自變量 k=-1;i=1;while (k0)k=ua(i,4)-ua(69,3);i=i+1;en