Matlab程序模拟汽车理论中的数据特性图1.doc

车辆理论作业第三组 第一大题确定一辆轻型货车的动力性能 1、绘制发动机的使用外特性曲线 n=600:1:4000 N=n/1000 T=-19.313+295.27*n/1000-165.44*N.*N+40.874*N.*N.*N-3.8445*N.*N.*N.*N Pe=T.*n/9550 %发动机的净功率 figure(1) plot(n,T) xlabel(转速n) ylabel(转矩T) figure(2) plot(n,Pe) xlabel(转速n) ylabel(功率Pe) 2、绘制驱动力---行驶阻力平衡图 i0=6.17; %减速器的传动比 ig(1)=6.09 % 一档的传动比 ig(2)=3.09 ig(3)=1.71 ig(4)=1.00 nT=0.85 r=0.367 f=0.013 m1=3880 m2=9590 CdA=2.77 Ff1=m1*9.8*f Ff2=m2*9.8*f ft1=T*ig(1)*i0*nT/r %一档所产生的驱动力 ft2=T*ig(2)*i0*nT/r ft3=T*ig(3)*i0*nT/r ft4=T*ig(4)*i0*nT/r Ft=[ft1;ft2;ft3;ft4]; ua1=0.377*r*n/(ig(1)*i0) %一档的速度 ua2=0.377*r*n/(ig(2)*i0) ua3=0.377*r*n/(ig(3)*i0) ua4=0.377*r*n/(ig(4)*i0) ua=[ua1;ua2;ua3;ua4]; fw1=CdA*ua1.^2/21.15+Ff1 %一档的行驶阻力 fw2=CdA*ua2.^2/21.15+Ff1 fw3=CdA*ua3.^2/21.15+Ff1 fw4=CdA*ua4.^2/21.15+Ff1 Fz1=[fw1;fw2;fw3;fw4]; figure(3) plot(ua,Ft) %绘制各档的驱动力曲线 hold on plot(ua,Fz1) %绘制行驶阻力曲线 hold on plot(ua4,Ff2) %绘制滚动阻力曲线 xlabel(ua/(km.h)) ylabel(F/N) title(汽车驱动力-行驶阻力平衡图) text(20,12000,Ft1); text(30,6000,Ft2); text(40,4000,Ft3); text(80,3000,Ft4); text(90,2200,（空载）Ff+fw); text(105,1450,满载Ff) 一档的驱动力很大，可用于爬坡且坡度较大；二档用于直接原地起步加速 3、绘制动力特性图 D=(Ft-Fw)/(m1*9.8) %求动力因数 f1=0.0076+0.000056*ua %一定车速条件下的滚动阻力 figure(4) plot(ua,D) hold on plot(ua,f1) text(20,0.3,I) text(20,0.154,II) text(40,0.1,III) text(80,0.05,IV) text(110,0.025,f) title(汽车动力特性图) xlabel(ua/(km.h)) ylabel(D) If=0.218 Iw1=1.798 Iw2=3.598 die1=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(1).^2)/(m1*r^2);%一档的旋转质量换算系数 die2=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(2).^2)/(m1*r^2); die3=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(3).^2)/(m1*r^2); die4=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(4).^2)/(m1*r^2); die=[die1;die2;die3;die4]; a1=9.8*D(:,1)/die1-9.8*f/die1; %求一档的加速度 a2=9.8*D(:,2)/die2-9.8*f/die2; a3=9.8*D(:,3)/die3-9.8*f/die3; a4=9.8*D(:,4)/die4-9.8*f/die4; aI=1./a1 %求一档的加速度倒数 aII=1./a2 aIII=1./a3 aIV=1./a4 aa=[aI,aII,aIII,aIV]; figure(5) %绘制各档的加速度倒数曲线 plot(ua1,aI) hold on plot(ua2,aII) hold on plo