汽車?yán)碚?-4MATLAB編程

PAGEPAGE4汽車?yán)碚?-4MATLAB編程figure(2);plot(u1,h1,u2,h2,u3,h3,u4,h4,'LineWidth',2);title('加速度倒數(shù)-速度曲線圖');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')%由加速度倒數(shù)-速度曲線圖可知u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=70;x1=[];x2=[];x3=[];x4=[];y=3401;fori=1:3401;ifu3(i)<=u3min;x1=[i];endendq1=max(x1);ua3=u3(q1:y);a3=h3(q1:y);fori=1:3401;ifu4(i)<=u4min;x2=[i];elseifu4(i)<=u4max;x3=[i];endendq2=max(x2);q3=max(x3);ua4=u4(q2:q3);a4=h4(q2:q3);s1=trapz(h2,u2);%二擋運(yùn)行時(shí)間s2=trapz(ua3,a3);s3=trapz(ua4,a4);s=[s1s2s3];disp('積分得')t=sum(s)*1000/3600%總時(shí)間《第二章》Pe1=Ft1.*u1./3600;%計(jì)算各檔對(duì)應(yīng)轉(zhuǎn)速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;P1=F1.*u1./(3600*x);%計(jì)算各檔對(duì)應(yīng)的各個(gè)車速下的行駛功率P2=F2.*u2./(3600*x);P3=F3.*u3./(3600*x);P4=F4.*u4./(3600*x);figure(3);plot(u1,Pe1,'-r',u2,Pe2,'-m',u3,Pe3,'-k',u4,Pe4,'-b',u1,P1,'k',u2,P2,'k',u3,P3,'k',u4,P4,'k','linewidth',2);gtext('Pe1')gtext('Pe2')gtext('Pe3')gtext('Pe4')xlabel('u/(km/h)');ylabel('P/kW');title('汽車功率平衡圖');n=[8151207161420122603300634033804];Ua=[];Ua=0.377*r*n./(io*ig(4))ft=[];ft=m*g*f+(2.77/21.15)*Ua.^2;%計(jì)算各轉(zhuǎn)速對(duì)應(yīng)的各個(gè)車速下的行駛阻力Pe(1)=ft(1).*Ua(1)./(3600*x);%計(jì)算各轉(zhuǎn)速對(duì)應(yīng)的各個(gè)車速下的行駛阻力功率Pe(2)=ft(2).*Ua(2)./(3600*x);Pe(3)=ft(3).*Ua(3)./(3600*x);Pe(4)=ft(4).*Ua(4)./(3600*x);Pe(5)=ft(5).*Ua(5)./(3600*x);Pe(6)=ft(6).*Ua(6)./(3600*x);Pe(7)=ft(7).*Ua(7)./(3600*x);Pe(8)=ft(8).*Ua(8)./(3600*x)B0=[1326.81354.71284.41122.91141.01051.21233.91129.7];B1=[-416.46-303.98-189.75-121.59-98.893-73.714-84.478-45.291];B2=[72.73936.65714.5257.00354.47632.85932.97880.7113];B3=[-5.8629-2.0533-0.51184-0.18517-0.091077-0.05138-0.047449-0.00075215];B4=[0.177680.0430720.00681640.00185550.000689060.000350320.00028230-0.000038568];b1=(B0(1))+(B1(1)*Pe(1))+(B2(1)*Pe(1)^2)+(B3(1)*Pe(1)^3)+(B4(1)*Pe(1)^4);b2=(B0(2))+(B1(2)*Pe(2))+(B2(2)*Pe(2)^2)+(B3(2)*Pe(2)^3)+(B4(2)*Pe(2)^4);b3=(B0(3))+(B1(3)*Pe(3))+(B2(3)*Pe(3)^2)+(B3(3)*Pe(3)^3)+(B4(3)*Pe(3)^4);b4=(B0(4))+(B1(4)*Pe(4))+(B2(4)*Pe(4)^2)+(B3(4)*Pe(4)^3)+(B4(4)*Pe(4)^4);b5=(B0(5))+(B1(5)*Pe(5))+(B2(5)*Pe(5)^2)+(B3(5)*Pe(5)^3)+(B4(5)*Pe(5)^4);b6=(B0(6))+(B1(6)*Pe(6))+(B2(6)*Pe(6)^2)+(B3(6)*Pe(6)^3)+(B4(6)*Pe(6)^4);b7=(B0(7))+(B1(7)*Pe(7))+(B2(7)*Pe(7)^2)+(B3(7)*Pe(7)^3)+(B4(7)*Pe(7)^4);b8=(B0(8))+(B1(8)*Pe(8))+(B2(8)*Pe(8)^2)+(B3(8)*Pe(8)^3)+(B4(8)*Pe(8)^4);p=0.7;Qs=[];Qs(1)=(Pe(1)*b1)/(1.02*Ua(1).*p*g);Qs(2)=(Pe(2)*b2)/(1.02*Ua(2).*p*g);Qs(3)=(Pe(3)*b3)/(1.02*Ua(3).*p*g);Qs(4)=(Pe(4)*b4)/(1.02*Ua(4).*p*g);Qs(5)=(Pe(5)*b5)/(1.02*Ua(5).*p*g);Qs(6)=(Pe(6)*b6)/(1.02*Ua(6).*p*g);Qs(7)=(Pe(7)*b7)/(1.02*Ua(7).*p*g);Qs(8)=(Pe(8)*b8)/(1.02*Ua(8).*p*g);M=polyfit(Ua,Qs,2);UA=0.377*r*600/(ig(4)*io):1:0.377*r*4000/(ig(4)*io);%UA表示車速Q(mào)S=polyval(M,UA);%QS表示油耗figure(4);plot(UA,QS,'linewidth',2);title('最高檔等速百公里油耗曲線');xlabel('Ua/(km/h)');ylabel('Qs/L');《第三章》io=[5.175.435.836.176.33];Va=0.377*r*n(7)./(io.*ig(4));Ps=46.9366;disp('假設(shè)以最高檔，較高轉(zhuǎn)速（n取3401），最經(jīng)濟(jì)負(fù)荷（即90%負(fù)荷大約46.9366Kw）行駛時(shí)油耗')qs=[];qs(1)=(Ps*b7)/(1.02*Va(1).*p*g);qs(2)=(Ps*b7)/(1.02*Va(2).*p*g);qs(3)=(Ps*b7)/(1.02*Va(3).*p*g);qs(4)=(Ps*b7)/(1.02*Va(4).*p*g);qs(5)=(Ps*b7)/(1.02*Va(5).*p*g);st=[17.581316.212114.512613.377512.9185];%加速時(shí)間：（這里以最高檔〈四檔〉、速度由0加速到94.93Km/h的時(shí)間）因與題1.3第三問(wèn)求法相同，這里不在累述，可直接有計(jì)算機(jī)求得：figure(5);plot(qs,st,'+','linewidth',2)holdonplot(qs,st,'linewidth',2);gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')title('燃油積極性-加速時(shí)間曲線');xlabel('燃油經(jīng)濟(jì)性(qs/L)');ylabel('動(dòng)力性--原地起步加速時(shí)間(st/s)');《第四章》4-3載荷質(zhì)量m/kg質(zhì)心高h(yuǎn)g/m軸距L/m質(zhì)心至前軸距離a/m制動(dòng)力分配系數(shù)β空載40800.8453.9502.1000.38滿載92901.1703.9502.9500.381）前軸利用附著系數(shù)為：后軸利用附著系數(shù)為：空載時(shí)：=所以空載時(shí)后輪總是先抱死。滿載時(shí)：=時(shí)：前輪先抱死時(shí)：后輪先抱死利用MATLAB作圖得到下圖：2）由MATLAB計(jì)算知當(dāng)φ=0.8時(shí)，(空載)制動(dòng)效率為0.6720所以其最大制動(dòng)減速度為代入公式：=6.57m計(jì)算可得：滿載時(shí)制動(dòng)效率為0.87因此其最大動(dòng)減速度制動(dòng)距離=5.34mA.若制動(dòng)系前部管路損壞后軸利用附著系數(shù)后軸制動(dòng)效率代入數(shù)據(jù)得：空載時(shí)：=0.45滿載時(shí)：=0.60a)空載時(shí)其最大動(dòng)減速度代入公式：=10.09mb)滿載時(shí)其最大動(dòng)減速度代入公式：=7.63mB．若制動(dòng)系后部管路損壞前軸利用附著系數(shù)前軸制動(dòng)效率代入數(shù)據(jù)空載時(shí)：=0.57滿載時(shí)：=0.33a)空載時(shí)其最大動(dòng)減速度代入公式：=8.02mb)滿載時(shí)其最大動(dòng)減速度代入公式：=13.67m4-51）同步附著系數(shù)2）因所以前輪先抱死==0.9513）最大制動(dòng)減速度：=4）a)1失效2失效b)1失效2失效c)1失效2失效5）a)1失效后軸利用附著系數(shù)后軸制動(dòng)效率0.46最大動(dòng)減速度2