《《汽车理论》第五版13题Matlabcode》由会员分享，可在线阅读，更多相关《《汽车理论》第五版13题Matlabcode（9页珍藏版）》请在金锄头文库上搜索。
1、.word可编辑.运行结果 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绘制汽车驱动力与行驶阻力平衡图. 专业.专注 .m=3880;r=0.367;nT=0.85
2、;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转对应的驱动力、车速和空气阻力 for t=1:40
3、00 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=Fx(k)%驱动力小于行驶阻力的数据无效 min1=Ft5(k)-Fx(k);%当前车速下驱动力与行驶阻力差值 end if minmin1 min=min1;%找到最小差值 result1=k,Ft5(k),Fx(k);%存储最高车速对应的转速（3200+k）、驱动力和行驶阻力 end endtext(0.377*(3200+result
4、1(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之间，Ft13100N,查得转速在1850-2280之间max=0;%存储驱动力与行驶阻力最大差值，以找出最大爬坡度点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).
5、2)/21.15+m*g*f;%1档的行驶阻力 max1=Ft1(k)-Fx(k);%当前车速下的驱动力与行驶阻力差值 if maxmax1 max=max1;%找出最大差值 result2=k,Ft1(k),Fx(k);%存储最大爬坡度对应的转速（1850+k）、驱动力和行驶阻力 end endtext(0.377*(1850+result2(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)%最大爬坡
6、度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)=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 t3310 acc(k,t)=N
7、aN;%最高车速对应的转速为3310，之后就没有加速能力了 end endendfigure(2)%图2plot(ua(1,:),acc(1,:),ua(2,:),acc(2,:),ua(3,:),acc(3,:),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 t3310 acc1(k,t)=NaN;%最高车速对应的转速为3310，之后就没有加速能力了 end e
8、ndend 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 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档
9、起步，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; end if k=4&ua(k-1,4000)ua(k,t)&ua(k,t)70%4档加速至70公里，从4档对应的3档最高车速开始 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 endend time70=time(z-1)
