汽車理論習題Maab程序

1、1.3確定一輕型貨車的動力性能（貨車可裝用4擋或5擋變速器，任選其中的一種進行整車性能計算）：1）繪制汽車驅(qū)動力與行駛阻力平衡圖。2）求汽車最高車速，最大爬坡度及克服該坡度時相應的附著率。3） 繪制汽車行駛加速度倒數(shù)曲線，用圖解積分法求汽車用2檔起步加速行 駛至70km/h的車速一時間曲線，或者用計算機求汽車用 2檔起步加速行駛至 70km/h的加速時間。輕型貨車的有關數(shù)據(jù)：汽油發(fā)動機使用外特性的Tq-n曲線的擬合公式為Tq19.313 295.27（）165.44（）2 40.874（）3 3.8445（）41000 1000 1000 1000式中，Tq為發(fā)動機轉矩（N?m）n為發(fā)動機轉速

2、（r/min ）。發(fā)動機的最低轉速nmin=600r/min,最高轉速nmax=4000r/min。裝載質(zhì)量2000kg整車整備質(zhì)量1800kg總質(zhì)量3880kg車輪半徑0.367m傳動系機械效率n=0.85f=0.013空氣阻力系數(shù)X迎風面積CA=2.77m2主減速器傳動比i0=5.83飛輪轉動慣量lf=0.218kg?m一前輪轉動慣量Iw1=1.798kg?m四后輪轉動慣量Iw2=3.598kg?m變速器傳動比ig（數(shù)據(jù)如下表）滾動阻力系數(shù)I檔U檔川檔W檔V檔四檔變速器6.093.091.711.00-五檔變速器5.562.7691.6441.000.793軸距L=3.2m質(zhì)心至前軸距離（

3、滿載）a=1.974m質(zhì)心高（滿載）hg=0.9m解：Matlab程序:（1）求汽車驅(qū)動力與行駛阻力平衡圖和汽車最高車速程序:n=600:10:4000;Tq=-19.313+295.27*( n/1000)-165.44*(門/1000)八2+40.874*(門/1000)八3-3.8445*( n/1000).A4；m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=5.56 2.769 1.644 1.00 0.793;nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218

4、;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;ua=0:5:120;Ff=G*f;Fw=CDA*ua.A2/21.15;Fz=Ff+Fw;plot(ua1,Ft1,

5、ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);title(驅(qū)動力-行駛阻力平衡圖);xlabel(ua(km/s);ylabel(Ft(N);gtext(Ft1),gtext(Ft2),gtext(Ft3),gtext(Ft4),gtext(Ft5),gtext(Ff+Fw);zoom on;x,y=ginput(1);zoom off;disp(汽車最高車速=);disp(x);disp(km/h);汽車最高車速 =99.3006 km/h(2) 求汽車最大爬坡度程序：n=600:10:4000;Tq=-19.313+295.27*( n/1000)-165.

6、44*(門/1000)八2+40.874*(門/1000)八3-3.8445*( n/1000).A4；m=3880;g=9.8; nmi n=600; nm ax=4000;G=m*g;ig=5.56 2.769 1.644 1.00 0.793; nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;lf=0.218;lw1= 1.798;lw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw仁CDA*ua2/21.15;Fz1=Ff+Fw1;Fi1

7、=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G);disp(汽車最大爬坡度=);disp(imax);disp(%);汽車最大爬坡度 =35.2197%(3) 求最大爬坡度相應的附著率和求汽車行駛加速度倒數(shù)曲線程序：clear n=600:10:4000;Tq=-19.313+295.27*( n/1000)-165.44*(門/1000)八2+40.874*(門/1000)八3-3.8445*( n/1000).A4；m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=5.56 2.769 1.644 1.00 0.79

8、3;nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*

9、r*n/ig(5)/i0;Fw1=CDA*ua1.A2/21.15;Fw2=CDA*ua2.A2/21.15;Fw3=CDA*ua3.A2/21.15;Fw4=CDA*ua4.A2/21.15;Fw5=CDA*ua5.A2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*rA2)+(If*ig(1)A2*i0A2*nT)/(m*rA2);deta2=1+(Iw1+Iw2)/(m*rA2)+(If*ig(2)A2*i0A2*nT)/(m*rA2);deta3=1+(Iw1+Iw2)/(m*rA2)+(If*ig(3)A2*i0A2*nT)/(m*rA2);deta4=1+(I

10、w1+Iw2)/(m*rA2)+(If*ig(4)A2*i0A2*nT)/(m*rA2);deta5=1+(Iw1+Iw2)/(m*rA2)+(If*ig(5)A2*i0A2*nT)/(m*rA2);a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2

11、,ua3,ad3,ua4,ad4,ua5,ad5);axis(0 99 0 10);title（汽車的加速度倒數(shù)曲線）;xlabel(ua(km/h);ylabel(1/a);gtext(1/a1);gtext(1/a2);gtext(1/a3);gtext(1/a4);gtext(1/a5);a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp(假設后輪驅(qū)動，最大爬坡度相應的附著率=);disp(C);假設后輪驅(qū)動，最大爬坡度相應的附著率0.4219(4) clear nT=0.85;r=0.367;f=

12、0.013;CDA=2.77;i0=5.83;lf=0.218;lw1=1.798;lw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=5.56 2.769 1.644 1.00 0.793;nmin=600;nm ax=4000;u1=0.377*r* nmin ./ig/i0;u2=0.377*r* nm ax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(lw1+lw2)/(m*rA2)+(lf*(ig(i)F2*iOA2* nT)/(m*A2);endua=6:O.O1:99;N=length(ua);

13、n=O;Tq=O;Ft=O;inv_a=O*ua;delta=O*ua;Ff=G*f;Fw=CDA*ua.A2/21.15;for i=1:Nk=i;if ua(i)=u2(2)n=ua(i)*(ig(2)*iO/r)/O.377;Tq=-19.313+295.27*( n/1000)-165.44*( n/1000F2+40.874*( n/1000F3-3.8445*( n/1000)A4;Ft=Tq*ig(2)*iO*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i);delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)=u2(3)n=

14、ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*( n/1000)-165.44*( n/1000F2+40.874*( n/1000F3-3.8445*( n/1000)A4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i);delta(i)=0.01*inv_a(i)/3.6;elseif ua(i) ginput ans =25.822370.073725.746770.0737所以汽車2檔原地起步換擋加速行駛至70km/h的加速時間約為25.8s汽車衛(wèi)輕慮地起步換苕抑連時同曲緘10 2030

15、4350時間i (s)60703000範307OS05O4O302O1012.7已知貨車裝用汽油發(fā)動機的負荷特性與萬有特性。負荷特性曲線的擬合公式為：b Bo BiPe B2Pe2 民巳3 B4FV其中，b為燃油消耗率g/（kW?h）Pe為發(fā)動機凈功率（kW）；擬合式中的系數(shù)隨轉速n變化。怠速油耗Qid 0.299mL/s （怠速轉速400r/min ）。計算與繪制題1.3中貨車的1）汽車功率平衡圖。2）最高檔與次高檔的等速百公里油耗曲線。或利用計算機求貨車按JB3352-83規(guī)定的六工況循環(huán)行駛的百公里油耗。計算中確定燃油消耗值b時, 若發(fā)動機轉速與負荷特性中給定的轉速不相等， 可由相鄰轉速

16、的兩根曲線用插值 法求得。解：Matlab程序：(1) 汽車功率平衡圖程序：clearn=600:10:4000;Tq=-19.313+295.27*( n/1000)-165.44*(門/1000)八2+40.874*(門/1000)八3-3.8 445* (n/1000).A4;m=3880;g=9.8;G=m*g;ig=5.56 2.769 1.644 1.00 0.793;nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;ua1=0.377*r*n/ig(

17、1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Pe1=Tq.*ig(1)*i0.*ua1./(3600*r);Pe2=Tq.*ig(2)*i0.*ua2./(3600*r);Pe3=Tq.*ig(3)*i0.*ua3./(3600*r);Pe4=Tq.*ig(4)*i0.*ua4./(3600*r);Pe5=Tq.*ig(5)*i0.*ua5./(3600*r);ua=0:0.35:119;Ff=G*f;Fw=CDA*ua.A2/21.15;Pf

18、=Ff*ua/3600;Pw=Fw.*ua/3600;Pe0=(Pf+Pw)./nT;Pe=max(Pe1);plot(ua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua5,Pe5,ua,Pe0,ua,Pe);axis(0 119 0 100);title（汽車功率平衡圖）;xlabel(ua(km/h);ylabel(Pe(kw);gtext(1),gtext(2),gtext(3),gtext(4),gtext(5),gtext(Pf+Pw)/et),gtext(Pe);(2)最高檔與次高檔的等速百公里油耗曲線程序:clearn=600:1:4000;m=3880;g=9

19、.8;G=m*g;ig=5.56 2.769 1.644 1.00 0.793;nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;n0=815 1207 1614 2012 2603 3006 3403 3804;B00=1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7;B10=-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291

20、;B20=72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113;B30=-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215;B40=0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568;B0=spline(n0,B00,n);B1=spline(n0,B10,n);B2=spline(n0,B20,n);B3=spline(n0

21、,B30,n);B4=spline(n0,B40,n);Ff=G*f;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fz4=Ff+CDA*(ua4.A2)/21.15;Fz5=Ff+CDA*(ua5.A2)/21.15;Pe4=Fz4.*ua4./(nT*3.6*1000);Pe5=Fz5.*ua5./(nT*3.6*1000);for i=1:1:3401b4(i)=B0(i)+B1(i)*Pe4(i)+B2(i)*Pe4(i)A2+B3(i)*Pe4(i)A3+B4(i)*Pe4(i)A4;b5(i)=B0(i)+B1(i)*Pe5(i)+B2

22、(i)*Pe5(i)A2+B3(i)*Pe5(i)A3+B4(i)*Pe5(i)A4;endpg=7.0;Q4=Pe4.*b4./(1.02.*ua4.*pg);Q5=Pe5.*b5./(1.02.*ua5.*pg);plot(ua4,Q4,ua5,Q5);axis(0 100 10 30);title(最高檔與次高檔等速百公里油耗曲線)；xlabel(ua(km/h);ylabel(百公里油耗(L/100km);gtext(4),gtext(5);最高杪與次高秒尊速百公里油託曲線3C262624222018161412W1020 3D 405060708090 1DQua(km/h)3.1改

23、變1.3題中輕型貨車的主減速器傳動比，做出i。為5.17、5.43 5.83 6.17、6.33時的燃油經(jīng)濟性一加速時間曲線，討論不同io值對汽車性能的影響。解：Matlab程序：主程序：i0=5.17,5.43,5.83,6.17,6.33;輸入主傳動比的數(shù)據(jù)for i=1:1:5y(i)=jiasushijian(i0(i); %求加速時間endy;for i=1:1:5b(i)=youhao(iO(i); %求對應iO的六工況百公里油耗endb;plot(b,y,+r)hold onb1=linspace(b(1),b(5),1OO);y1=spline(b,y,b1); % 三次樣條插

24、值plot(b1,y1); % 繪制燃油經(jīng)濟性 -加速時間曲線title(燃油經(jīng)濟性一加速時間曲線);xlabel(百公里油耗(L/1OOkm);ylabel(加速時間s);gtext(iO=5.17),gtext(iO=5.43),gtext(iO=5.83),gtext(iO=6.17),gtext(iO=6.33);子程序：(1) function y=jiasushijian(i0) %求加速時間的處理函數(shù) n1=linspace(0,5000); %先求各個檔位的驅(qū)動力nmax=4000;nmin=600;r=0.367;yita=0.85;CDA=2.77;f=0.013;G=(3

25、880)*9.8;ig=6.09,3.09,1.71,1.00;%i0=5.83%i為檔數(shù)%計算各個檔位的最大速度與最for i=1:1:4uamax(i)=chesu(nmax,r,ig(i),i0);小速度uamin(i)=chesu(nmin,r,ig(i),i0);ua(i,:)=linspace(uamin(i),uamax(i),100);%計算各個檔位的轉速范圍%求出各檔位的轉矩范圍%求出驅(qū)動力%求出滾動阻力和空氣阻n(i,:)=zhuansu(ua(i,:),r,ig(i),i0);Ttq(i,:)=zhuanju(n(i,:);Ft(i,:)=qudongli(Ttq(i,:

26、),ig(i),i0,yita,r);F(i,:)=f*G+CDA*(ua(i,:)A2)/21.15;力的和delta(i,:)=1+(1.798+3.598+0.218*(ig(i)A2)*(i0A2)*yita)/(3880W2); %求轉動質(zhì)量換算系數(shù)a(i,:)=1./(delta(i,:).*3880./(Ft(i,:)-F(i,:); %求出加速度F2(i,:)=Ft(i,:)-F(i,:);end %下面分各個檔位進行積分，求出加速時間temp1(1,:)=ua(2,:)/3.6;temp1(2,:)=1./a(2,:);n1=1;for j1=1:1:100if ua(3,j

27、1)max(ua(2,:)&ua(3,j1)max(ua(3,:)&ua(4,j1)=70;temp3(1,n2)=ua(4,j1)/3.6;temp3(2,n2)=1./a(4,j1);endendy=temp1(1,1)*temp1(2,1)+qiuji(temp1(1,:),temp1(2,:)+qiuji(temp2(1,:),temp2(2,:)+qiuji (temp3(1,:),temp3(2,:);end(2) function ua=chesu(n,r,ig,i0); %由轉速計算車速ua=0.377*r.*n/(ig*i0);(3) function n=zhuansu(u

28、a,r,ig,i0); %求轉速n=ig*i0.*ua./(0.377*r);end(4) function y=zhuanju(n);%求轉矩函數(shù)y=-19.313+295.27.*( n./1000)-165.44.*(門./1000)八2+40.874.*(門./1000)八3-3.8445*( n./1 000)八4;(5) function y=qudongli(Ttq,ig,i0,yita,r); %求驅(qū)動力函數(shù) y=(ig*i0*yita.*Ttq)/r;end(6) function p=qiuji(x0,y0)%求積分函數(shù)n0=size(x0);n=n0(2);x=linsp

29、ace(x0(1),x0(n),200) ;y=spline(x0,y0,x);%插值% figure;plot(x,y);p=trapz(x,y) ;end求不同iO下的六工況油耗function b=youhao(i0);global f G CDA yita m r If Iw1 Iw2 pg BO B1 B2 B3 B4 n %聲明全局變量ig=6.09,3.09,1.71,1.00;r=0.367;yita=0.85;CDA=2.77;f=0.013;%i0=5.83;G=(3880)*9.8;If=0.218;Iw1=1.798;Iw2=3.598;m=3880; %汽車的基本參數(shù)

30、設定n0=815 1207 1614 2012 2603 3006 3403 3804;B00=1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7;B10=-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291;B20=72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113;B30=-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449

31、-0.00075215;B40=0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230-0.000038568;n=600:1:4000;B0=spline(n0,B00,n);B1=spline(n0,B10,n);B2=spline(n0,B20,n);%使用三次樣條插值，保證曲線的光滑連續(xù)B3=spline(n0,B30,n);B4=spline(n0,B40,n);ua4=0.377*r.*n./(i0*ig(4);%求出發(fā)動機轉速范圍內(nèi)對應的 III 、IV 檔車速F4=f*G+CDA*(ua4.A2

32、)/21.15;%求出滾動阻力和空氣阻力的和%求出阻力功率for i=1:1:3401%用擬合公式求出各個燃油消耗率P_fw4=F4.*ua4./(yita*3.6*1000);b4(i)=B0(i)+B1(i)*P_fw4(i)+B2(i)*(P_fw4(i)A2+B3(i)*(P_fw4(i)F3+B4(i)*(P_fw4(i)F4;endpg=7.06;%汽油的重度取 7.06N/Lua4_m=25,40,50;%勻速階段的車速s_m=50,250,250;%每段勻速走過的距離b4_m=spline(ua4,b4,ua4_m);%插值得出對應速度的燃油消耗率%車速對應的阻力F4_m=f*

33、G+CDA*(ua4_m.A2)/21.15;P_fw4_m=F4_m.*ua4_m./(yita*3.6*1000); %發(fā)動機功率Q4_m=P_fw4_m.*b4_m.*s_m./(102.*ua4_m.*pg) ;Q4_a1=jiasu(40,25,ig(4),0.25,ua4,i0);Q4_a2=jiasu(50,40,ig(4),0.2,ua4,i0);Qid=0.299;tid=19.3;s=1075;Q_i=Qid*tid;%求出減速階段的燃油消耗量Q4all=(sum(Q4_m)+Q4_a1+Q4_a2+Q*100/s; %IV 檔六工況百公里燃油消耗量b=Q4all;(8)加

34、速階段處理函數(shù)function q=jiasu(umax,umin,ig,a,ua0,i0);global f G CDA yita m r If Iw1 Iw2 pg B0 B1 B2 B3 B4 n; %i0 ;ua1=umin:1:umax;%把速度范圍以 1km/h 為間隔進行劃分delta=1+(Iw1+Iw2)/(m*rA2)+(If*igA2*i0A2*yita)/(m*rA2);PO=(G*f.*uaO./36OO+CDA*uaO.A3/7614O+(delta*m.*uaO/36OO)*a)/yita;P=(G*f.*ua1/3600+CDA.*ua1.A3/76140+(d

35、elta*m.*ua1/3600)*a)/yita;dt=1/(3.6*a) ;%速度每增加 1km/h 所需要的時間for i=1:1:3401%重新利用擬合公式求出b與ua的關系bO(i)=BO(i)+B1(i)*PO(i)+B2(i)*(PO(i)A2+B3(i)*(PO(i)A3+B4(i)*(PO(i)A4;endb1=interp1(uaO,bO,ua1);%插值出各個速度節(jié)點的燃油消耗率Qt=P.*b1./(367.1.*pg);%求出各個速度節(jié)點的燃油消耗率i1=size(Qt);i=i1(2);Qt1= Qt(2:i-1);q=(Qt(1)+Qt(i)*dt./2+sum(Q

36、t1)*dt;%求該加速階段的燃油消耗量燃油經(jīng)濟性一扣速時間曲線2827.827.5227 22726.826.526.42613.313 413.513.613.713.81.91414.114 214.3百公里油耗(LnOOkm)4.3 一中型貨車裝有前后制動器分開的雙管路制動系，其有關參數(shù)如下:載荷質(zhì)量（kg）質(zhì)心高hg/m軸距L/m質(zhì)心至前軸距離a/m制動力分配系數(shù)B空載40800.8453.9502.1000.38滿載92901.1703.9502.9500.381）計算并繪制利用附著系數(shù)曲線和制動效率曲線2）求行駛車速Ua= 30km/h，在 =0.80路面上車輪不抱死的制動距離。

37、計算時取制動系反應時間2 = 0.02s，制動減速度上升時間2 = 0.02s。3）求制動系前部管路損壞時汽車的制動距離s,制動系后部管路損壞時汽車的制動距離 s 。解：Matlab程序：(1) 求利用附著系數(shù)曲線和制動效率曲線程序：cleark=4080;hgk=0.845;Lk=3.950;ak=2.10;betak=0.38;bk=Lk-ak;空載時的參數(shù) mm=9290;hgm=1.170;Lm=3.950;am=2.950;betam=0.38;bm=Lm-am;滿載時的參數(shù)z=0:0.01:1.0;figure(1);fai=z;fai_fk=betak*z*Lk./(bk+z*h

38、gk);% 空載時前軸的 f fai_fm=betam*z*Lm./(bm+z*hgm);% 滿載時前軸的 f fai_rk=(1-betak)*z*Lk./(ak-z*hgk);% 空載時后軸的 r fai_rm=(1-betam)*z*Lm./(am-z*hgm);%滿載時后軸的 r plot(z,fai_fk,b-,z,fai_fm,r,z,fai_rk,b-,z,fai_rm,r,z,fai,k);title（利用附著系數(shù)與制動強度的關系曲線）;xlabel(制動強度(z/g);ylabel(利用附著系數(shù))；gtext(|r(空載),gtext(|r(滿載),gtext(=z),gte

39、xt(旅空載),gtext(旅滿載); figure(2);Efk=z./fai_fk*100;%空載時前軸的制動效率Efm=z./fai_fm*100;Erk=z./fai_rk*100;Erm=z./fai_rm*100;plot(fai_fk,Efk,b,fai_fm,Efm,r,fai_rk,Erk,b,fai_rm,Erm,r);axis(0 1 0 100);title(前.后制動效率曲線);xlabel(附著系數(shù)妒)；ylabel(制動效率%);gtext(Ef),gtext(Er),gtext(Er),gtext滿載),gtext(空 載);利用附著系數(shù)與制動強度的關系曲線on

40、ooOQ.10.20.30.40,50.6制動強度(z/ g)070,30.912 8 6/.21.8E1 10.0.4 2 QO.后制動效率曲線0.1020.304050.S0.7080.91附看系數(shù)屮9080403020(2)問和(3)問程序:clear mk=4080;hgk=0.845;Lk=3.950;ak=2.10;betak=0.38;bk=Lk-ak;空 載 時 的參數(shù)mm=9290;hgm=1.170;Lm=3.950;am=2.950;betam=0.38;bm=Lm-am;%滿 載時的參數(shù)z=0:0.01:1;fai_fk=betak*z*Lk./(bk+z*hgk);%

41、 空載時前軸的 f fai_fm=betam*z*Lm./(bm+z*hgm);% 滿載時前軸的 f fai_rk=(1-betak)*z*Lk./(ak-z*hgk);% 空載時后軸的 r fai_rm=(1-betam)*z*Lm./(am-z*hgm);%滿載時后軸的 rEfk=z./fai_fk*100;%空載時前軸的制動效率Efm=z./fai_fm*100;Erk=z./fai_rk*100;Erm=z./fai_rm*100;t1=0.02;t2=0.02;ua0=30;fai=0.80;g=9.8;ak1=Erk(81)*g*fai/100;am1=Erm(81)*g*fai/

42、100;Sk1=(t1+t2/2)*ua0/3.6+ua0A2/(25.92*ak1);%制動距離Sm仁(t1+t2/2)*ua0/3.6+ua0A2/(25.92*am1);disp(空載時，汽車制動距離Sk1=);disp(Sk1);disp(滿載時，汽車制動距離Sm1=);disp(Sm1);ak2=fai*g*ak/(Lk+fai*hgk);am2=fai*g*am/(Lm+fai*hgm);ak3=fai*g*bk/(Lk-fai*hgk);am3=fai*g*bm/(Lk-fai*hgm);Sk2=(t1+t2/2)*ua0/3.6+ua0A2/(25.92*ak2);%制 動距離

43、Sm2=(t1+t2/2)*ua0/3.6+ua0A2/(25.92*am2);Sk3=(t1+t2/2)*ua0/3.6+ua0A2/(25.92*ak3);Sm3=(t1+t2/2)*ua0/3.6+ua0A2/(25.92*am3);disp(空載時，前制動器損壞，汽車制動距離Sk2=);disp(Sk2);disp(滿載時，前制動器損壞，汽車制動距離Sm2=);disp(Sm2);Sk3=);Sm3=);disp(空載時，后制動器損壞，汽車制動距離disp(Sk3);disp(滿載時，后制動器損壞，汽車制動距離disp(Sm3);空載時，汽車制動距離 Sk1=7.8668滿載時，汽車制

44、動距離 Sm1=5.6354空載時，前制動器損壞，汽車制動距離 Sk2=10.0061滿載時，前制動器損壞，汽車制動距離 Sm2=7.5854空載時，后制動器損壞，汽車制動距離 Sk3=8.0879滿載時，后制動器損壞，汽車制動距離 Sm3=13.59865.11二自由度轎車模型的有關參數(shù)如下:總質(zhì)量繞Oz軸轉動慣量軸距質(zhì)心至前軸距離質(zhì)心至后軸距離前輪總側偏剛度后輪總側偏剛度轉向系總傳動比試求：m=1818.2kg2I z3885 kg mL=3.048ma=1.463mb=1.585m ki=-62618N/rad k2=-110185N/rad i=201) 穩(wěn)定性因數(shù)K、特征車速Uch。

45、2)穩(wěn)態(tài)橫擺角速度增益曲線Ua、車速u=22.35m/s時的轉向靈敏度r。sw3) 靜態(tài)儲備系數(shù) S.M.,側向加速度為0.4g時的前、后輪側偏角絕對值之差與轉彎半徑的比值 R/R 0(R0=15m)。0、阻尼4) 車速 u=30.56m/s 時，瞬態(tài)響應的橫擺角速度波動的固有(圓)頻率比 、反應時間 與峰值反應時間解：Matlab程序：m=1818.2;Iz=3885;L=3.048;a=1.463;b=1.585;k1=-62618;k2=-110185;i=20;g=9.8;R0=15;u1=30.56;K=m*(a/k2-b/k1)/LA2;Uch=(1/K)/2);% 特征車速dis

46、p(穩(wěn)定性因數(shù)(sA2/mA2)K=);disp(K);disp(特征車速(m/s)Uch=);disp(Uch);u=0:0.05:30;S=u./(L*(1+K*u.A2);% 穩(wěn)態(tài)橫擺角速度增益plot(u,S);title(汽車穩(wěn)態(tài)橫擺角速度增益曲線);xlabel(車速 u(m/s);ylabel(穩(wěn)態(tài)橫擺角速度增益);disp（u=22.35m/s 時，轉向靈敏度為 ）;disp(S(448);SM=k2/(k1+k2)-a/L;ay=0.4*g;A=K*ay*L;B=L/R0;R=L/(B-A);C=R/RO;%轉彎半徑比disp(靜態(tài)儲備系數(shù)S.M.=);disp(SM);di

47、sp(側向加速度為0.4g時前、后輪側偏角絕對值之差(rad) a1-a2=); disp(A);disp(側向加速度為0.4g時轉彎半徑比值R/RO=);disp(C);W0=L/u1*(k1*k2/(m*lz)*(1+K*uM2)F(1/2);% 固有(圓)頻率D=(-m*(k1*aA2+k2*bA2)-Iz*(k1+k2)/(2*L*(m*Iz*k1*k2*(1+K*u1A2)A(1/2);%阻 尼比峰值反應時間);E=ata n( (1-DA2F(1/2)/D)/(W0*(1-DA2F(1/2)+t;% disp(車速u=30.56m/s時的瞬態(tài)響應參數(shù)分別為: disp(橫擺角速度波

48、動的固有(圓)頻率(rad)為); disp(W0);disp(阻尼比為);disp(D);disp(反應時間(s)為)；disp(t);disp(峰值反應時間(s)為)； disp(E);穩(wěn)定性因數(shù)(sA2/mA2)K=0.0024特征車速 (m/s)Uch=20.6053u=22.35m/s 時，轉向靈敏度為3.3690靜態(tài)儲備系數(shù) S.M.=0.1576(rad) a1-a2=側向加速度為 0.4g 時前、后輪側偏角絕對值之差0.0281側向加速度為 0.4g 時轉彎半徑比值 R/R0=1.1608車速 u=30.56m/s 時的瞬態(tài)響應參數(shù)分別為：橫擺角速度波動的固有(圓)頻率(rad

49、)為5.5758阻尼比為0.5892反應時間（s）為0.1811峰值反應時間（S）為0.3899汽車穩(wěn)態(tài)橫擺角速度堵益曲錢車 iS d(in/s)6.5車身-車輪雙質(zhì)量系統(tǒng)參數(shù)：f0 1.5Hz, 0.25,9,10?！叭梭w-座椅”系統(tǒng)參數(shù):3Hz, s 0.25。車速 u20m/ s，路面不平度系數(shù)3）分別改變車身-車輪雙質(zhì)量系統(tǒng)參數(shù):f00.25 3Hz, 0.125 0.5，Gq n02.56 10 8m3，參考空間頻率 n0=0.1m-1。計算時頻率步長f 0.2 Hz，計算頻率點數(shù) N 180。1）計算并畫出幅頻特性z,/q、.Gz、. Gr譜圖。進一步計算2）改變“人體-座椅”系統(tǒng)

50、參數(shù)：fsZ2/Z1、q/z2和均方根值譜 vGz f 、q、勺、Z2、a、aw、Law 值1.5 -6Hz, s0.125 0.5。分析 aw、Law值隨fs、 s的變化。4.5 18,520 。繪制z2fdFd/G 三個響應量均方根值隨以上四個系統(tǒng)參數(shù)變化的曲線。解：Matlab程序(1)問yps=0.25;%阻尼比Zgama=9;%剛度比丫mu=10;%質(zhì)量比卩fs=3;ypss=0.25;g=9.8;a0=10A(-6);f0=1.5;ua=20;Gq n0=2.56*10(8); n0=0.1;detaf=0.2;N=180;f=detaf*0:N;lamta=f/f0;lamtas=f/fs;Wf=0*f;deta=(1-lamta.A2).*(1+gama-1/mu*lamta.A2)-1).A2+4*ypsA2*lamta.A2.*(gama-(1/mu+1)*lamta.A2).A2;z1_q=gama*sqrt(1-lamta.A2).A2+4*ypsA2*lamta.A2)./deta);z2_z 仁 sqrt(1+4*ypsA2*lamta.A2)./(1-lam