《汽車理論》第五版1.3題Matlabcode

1、精選優(yōu)質文檔-傾情為你奉上運行結果>> 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 汽車行駛加速度倒數(shù)曲線MATLAB-Code%1.3.1繪制汽車驅動力與行駛阻力平衡圖專心-專注-專

2、業(yè)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;%發(fā)動機轉矩擬合公式for k=1:5%求出各檔位時發(fā)動機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轉的輔助數(shù)據(jù)點設置為非數(shù) 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求最高車速%利用圖形數(shù)據(jù)游標，可以得出最高車

5、速在96-100km/h之間,空氣阻力在1725-1823N之間，從工作區(qū)查得轉速在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)%驅動力小于行駛阻力的數(shù)據(jù)無效 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之間，F(xiàn)t>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繪制汽車行駛加速度倒數(shù)曲線for k=1:5 de(k)

9、=1+1/m*(Iw1+Iw2)2/r+1/m*If*ig(k)2*i02*nT/r2;%汽車旋轉質量換算系數(shù) 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轉的輔助數(shù)據(jù)點設置為非數(shù)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個輔助點設為非數(shù) else acc1(k,t)=1/acc(k,t);%各檔的加速度的倒數(shù) 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,:);%輸出加速度倒數(shù)曲線圖legend('1','2','3','4','5')%圖例title('汽車行駛加速度倒數(shù)曲線');axis(0

12、100 0 10);%設置圖形界限，最高車速之后就沒有加速能力了xlabel('u_a')ylabel('1/a')%1.3.3.2積分法求2檔起步加速至70公里的時間time=0;%二檔起步加速至70公里用時；z=1;%循環(huán)標記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;