汽车理论编程作业123章.doc

汽车理论作业：第一章1.3选四档变速器，用MATLAB语言驱动力与行驶阻力平衡图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*(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,ua1,fz1,ua2,fz2,ua3,fz3,ua4,fz4);title(驱动力与行驶阻力平衡图);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);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;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.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;Iw2=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(k);q=i.*l./a;disp(汽车的最大爬坡度=);disp(i);disp(最大爬坡度相应的附着率=);disp(q);得计算结果：汽车的最大爬坡度= 0.3648最大爬坡度相应的附着率= 0.5995加速度倒数曲线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.*i0);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);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)./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(加速度倒数曲线图);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.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.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)./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.71 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,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 3006 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*Pe1(1)4-5.8629*Pe1(1)3+72.379*Pe1(1)2-416.46*Pe1(1)+1326.8; b(2)=0.043072*Pe1(2)4-2.0553*Pe1(2)3+36.657*Pe1(2)2-303.98*Pe1(2)+1354.7; b(3)=0.0068164*Pe1(3)4-0.51184*Pe1(3)3+14.524*Pe1(3)2-189.75*Pe1(3)+1284.4; b(4)=0.0018555*Pe1(4)4-0.18517*Pe1(4)3+7.0035*Pe1(4)2-121.59*Pe1(4)+1122.9; b(5)=0.00068906*Pe1(5)4-0.091077*Pe1(5)3+4.4763*Pe1(5)2-98.893*Pe1(5)+1141.0; b(6)=0.00035032*Pe1(6)4-0.05138*Pe1(6)3+2.8593*Pe1(6)2-73.714*Pe1(6)+1051.2; b(7)=0.00028230*Pe1(7)4-0.047449*Pe1(7)3+2.9788*Pe1(7)2-84.478*Pe1(7)+1233.9; b(8)=-0.000038568*Pe1(8)4-0.00075215*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)*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=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;%选取最高挡，转速n=3403r/min,Pe=60kw，对应的油耗比较经济Pe=60;B=0.00028230*Pe4-0.047449*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)=-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 %求加速度倒数曲线 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,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)为加速度倒数函数 endendk=-1;i=1;while(k0)k=ua(i,2)-ua(69,1);i=i+1;endw=i-1; %w为从1挡换至2挡时的转速的自变量 k=-1;i=1; while (k0)k=ua(i,3)-ua(69,2);i=i+1;endx=i-1; %x为换至3挡时的转速的自变量 k=-1;i=1;while (k0)k=ua(i,4)-ua(