您好,欢迎访问三七文档
当前位置:首页 > 机械/制造/汽车 > 汽车理论 > 《汽车理论》第五版1.3题Matlabcode
运行结果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;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;%发动机转矩拟合公式fork=1:5%求出各档位时发动机600-4000转对应的驱动力、车速和空气阻力fort=1:4000Ft(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;%求不同车速对应的空气阻力ift600Ft(k,t)=NaN;%把1-600转的辅助数据点设置为非数endendendFf=m*g*f;%滚动阻力figure(1)%图1plot(ua(1,:),Ft(1,:),ua(2,:),Ft(2,:),ua(3,:),Ft(3,:),ua(4,:),Ft(4,:),...ua(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求最高车速%利用图形数据游标,可以得出最高车速在96-100km/h之间,空气阻力在1725-1823N之间,从工作区查得转速在3239-3380之间min=100000;%存储驱动力与行驶阻力最小差值,以找出最高车速点fork=1:160%findthefastestspeedFt5(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档时的行驶阻力ifFt5(k)=Fx(k)%驱动力小于行驶阻力的数据无效min1=Ft5(k)-Fx(k);%当前车速下驱动力与行驶阻力差值endifminmin1min=min1;%找到最小差值result1=[k,Ft5(k),Fx(k)];%存储最高车速对应的转速(3200+k)、驱动力和行驶阻力endendtext(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之间,Ft13100N,查得转速在1850-2280之间max=0;%存储驱动力与行驶阻力最大差值,以找出最大爬坡度点fork=1:430%findthefastestspeedFt1(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);%当前车速下的驱动力与行驶阻力差值ifmaxmax1max=max1;%找出最大差值result2=[k,Ft1(k),Fx(k)];%存储最大爬坡度对应的转速(1850+k)、驱动力和行驶阻力endendtext(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)%最大爬坡度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绘制汽车行驶加速度倒数曲线fork=1:5de(k)=1+1/m*(Iw1+Iw2)^2/r+1/m*If*ig(k)^2*i0^2*nT/r^2;%汽车旋转质量换算系数fort=1:4000acc(k,t)=1/(de(k)*m)*(Ft(k,t)-Ff-Fw(k,t));%各档的加速度ift600acc(k,t)=NaN;%把1-600转的辅助数据点设置为非数endifk==5&t3310acc(k,t)=NaN;%最高车速对应的转速为3310,之后就没有加速能力了endendendfigure(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}')fork=1:5fort=1:4000ift600acc1(k,t)=NaN;%前600个辅助点设为非数elseacc1(k,t)=1/acc(k,t);%各档的加速度的倒数endifk==5&t3310acc1(k,t)=NaN;%最高车速对应的转速为3310,之后就没有加速能力了endendendfigure(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([0100010]);%设置图形界限,最高车速之后就没有加速能力了xlabel('u_{a}')ylabel('1/a')%1.3.3.2积分法求2档起步加速至70公里的时间time=[0];%二档起步加速至70公里用时;z=1;%循环标记fork=2:4%由图中看出2档起步,4档即可加速至70公里fort=600:3999ifk==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;endifk==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;endifk==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;endendendtime70=time(z-1)
本文标题:《汽车理论》第五版1.3题Matlabcode
链接地址:https://www.777doc.com/doc-5910310 .html