《汽车理论》第五版1.3题Matlabcode

1、精选优质文档-倾情为你奉上运行结果>> Untitled3Result1 = 99.0589km/h-最高车速Result2 = 8.6862km/h-最大爬坡度时的车速A = 0.3386rad-最大坡角I = 0.3522-最大爬坡度Cf1 = 1.2041-前驱附着率（>1，舍弃）Cf2 = 0.4978-后驱克服该坡度的附着率time70 = 26.2088s-2档起步加速至70km/h的时间>>图1 汽车驱动力-行驶阻力平衡图图2 汽车行驶加速度曲线图3 汽车行驶加速度倒数曲线MATLAB-Code%1.3.1绘制汽车驱动力与行驶阻力平衡图专心-专注-专

2、业m=3880;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=5.56,2.769,1.644,1.00,0.793;L=3.2;a=1.947;hg=0.9;g=9.8;n=1:1:4000;%为了得到完整的图线，避免矢量长度不符，故从转速1开始Tq=-19.313+295.27*(n./1000)-165.44*(n./1000).2+40.874*(n./1000).3-. 3.8445*(n/1000).4;%发动机转矩拟合公式for k=1:5%求出各档位时发动机600-4000转对应

3、的驱动力、车速和空气阻力 for t=1:4000 Ft(k,t)=Tq(t)*ig(k)*i0*nT/r;%求各档的驱动力 ua(k,t)=0.377*n(t)*r/ig(k)/i0;%求各档的车速范围 Fw(k,t)=CdA*(ua(k,t).2)/21.15;%求不同车速对应的空气阻力 if t<600 Ft(k,t)=NaN;%把1-600转的辅助数据点设置为非数 end endendFf=m*g*f;%滚动阻力figure(1)%图1plot(ua(1,:),Ft(1,:),ua(2,:),Ft(2,:),ua(3,:),Ft(3,:),ua(4,:),Ft(4,:),. ua

4、(5,:),Ft(5,:),ua(5,:),Ff,ua(5,:),Ff+Fw(5,:);%输出汽车驱动力-行驶阻力平衡图legend('F_t1','F_t2','F_t3','F_t4','F_t5','F_w','F_f+F_w')%图例title('汽车驱动力-行驶阻力平衡图');%图表标题xlabel('u_a/km*h-1')%横坐标ylabel('F/N')%纵坐标%1.3.2.1求最高车速%利用图形数据游标，可以得出最高车

5、速在96-100km/h之间,空气阻力在1725-1823N之间，从工作区查得转速在3239-3380之间min=;%存储驱动力与行驶阻力最小差值，以找出最高车速点for k=1:160%find the fastest speed Ft5(k)=Tq(3200+k)*ig(5)*i0*nT/r;%5档的驱动力 Fx(k)=CdA*(0.377*(3200+k)*r/ig(5)/i0).2)/21.15+m*g*f;%5档时的行驶阻力 if Ft5(k)>=Fx(k)%驱动力小于行驶阻力的数据无效 min1=Ft5(k)-Fx(k);%当前车速下驱动力与行驶阻力差值 end if min

6、>min1 min=min1;%找到最小差值 result1=k,Ft5(k),Fx(k);%存储最高车速对应的转速（3200+k）、驱动力和行驶阻力 end endtext(0.377*(3200+result1(1)*r/ig(5)/i0,result1(2),'M')%标出最高车速点Result1=0.377*(3200+result1(1)*r/ig(5)/i0%得出最高车速值%1.3.2.2求最大爬坡度%由图得，1档时有最大驱动力，在车速为7-10km/h之间，Ft>13100N,查得转速在1850-2280之间max=0;%存储驱动力与行驶阻力最大差值，

7、以找出最大爬坡度点for k=1:430%find the fastest speed Ft1(k)=Tq(1850+k)*ig(1)*i0*nT/r;%1档时的驱动力 Fx(k)=CdA*(0.377*(1850+k)*r/ig(1)/i0).2)/21.15+m*g*f;%1档的行驶阻力 max1=Ft1(k)-Fx(k);%当前车速下的驱动力与行驶阻力差值 if max<max1 max=max1;%找出最大差值 result2=k,Ft1(k),Fx(k);%存储最大爬坡度对应的转速（1850+k）、驱动力和行驶阻力 end endtext(0.377*(1850+result2

8、(1)*r/ig(1)/i0,result2(2),'I')%标出最大爬坡度点Result2=0.377*(1850+result2(1)*r/ig(1)/i0%得出爬最大坡时的车速值A=asin(max/m/g)%能爬上的最大坡角I=tan(A)%最大爬坡度I%1.3.2.3求克服该坡度相应的附着率q=I+(1/cos(A)*(1/g)*(max-m*g*I*cos(A)/m);%等效坡度Cf2=q/(a/L+hg/L*q);%后轮附着率Cf1=q/(L-a)/L-hg/L*q);%前轮附着率Cf1Cf2%1.3.3.1绘制汽车行驶加速度倒数曲线for k=1:5 de(k)

9、=1+1/m*(Iw1+Iw2)2/r+1/m*If*ig(k)2*i02*nT/r2;%汽车旋转质量换算系数 for t=1:4000 acc(k,t)=1/(de(k)*m)*(Ft(k,t)-Ff-Fw(k,t);%各档的加速度 if t<600 acc(k,t)=NaN;%把1-600转的辅助数据点设置为非数endif k=5&t>3310 acc(k,t)=NaN;%最高车速对应的转速为3310，之后就没有加速能力了 end endendfigure(2)%图2plot(ua(1,:),acc(1,:),ua(2,:),acc(2,:),ua(3,:),acc(3

10、,:),ua(4,:),acc(4,:),. ua(5,:),acc(5,:);%输出加速度图形legend('1','2','3','4','5')%图例title('汽车行驶加速度曲线');xlabel('u_a/km*h-1')ylabel('a/m*s-2')for k=1:5 for t=1:4000 if t<600 acc1(k,t)=NaN;%前600个辅助点设为非数 else acc1(k,t)=1/acc(k,t);%各档的加速度的倒数 en

11、dif k=5&t>3310 acc1(k,t)=NaN;%最高车速对应的转速为3310，之后就没有加速能力了 end endend figure(3)%图3plot(ua(1,:),acc1(1,:),ua(2,:),acc1(2,:),ua(3,:),acc1(3,:),. ua(4,:),acc1(4,:),ua(5,:),acc1(5,:);%输出加速度倒数曲线图legend('1','2','3','4','5')%图例title('汽车行驶加速度倒数曲线');axis(0

12、100 0 10);%设置图形界限，最高车速之后就没有加速能力了xlabel('u_a')ylabel('1/a')%1.3.3.2积分法求2档起步加速至70公里的时间time=0;%二档起步加速至70公里用时；z=1;%循环标记for k=2:4%由图中看出2档起步，4档即可加速至70公里 for t=600:3999 if k=2%二档加速至2档最高车速然后换挡 time1=(ua(k,t+1)-ua(k,t)/3.6*(acc1(k,t+1)+acc1(k,t)/2;%时间积分微元 time(z+1)=time(z)+time1;%时间累加 z=z+1; end if k=3&ua(k-1,4000)<ua(k,t)%3档加速至3档最高车速然后换挡，从3档对应的2档最高车速开始 time1=(ua(k,t+1)-ua(k,t)/3.6*(acc1(k,t+1)+acc1(k,t)/2; time(z+1)=time(z)+time1; z=z+1;