汽車理論maab作業(yè)

1、確定一輕型貨車的動力性能1) 繪制汽車驅動力與行駛阻力平衡圖；2) 求汽車最高車速與最大爬坡度；3) 繪制汽車行駛加速度倒數(shù)曲線；用計算機求汽車用H檔起步加速行駛至70km/h所需的加速時間。已知數(shù)據(jù)略。(參見汽車理論習題第一章第3題)解題程序如下：用Mat lab語言(1)繪制汽車驅動力與行駛阻力平衡圖m 仁 2000;m2二1800;mz二3880;g二9. 81;r二0.367;CdA二2. 77; f二0.013;nT二0.85;ig二5. 56 2. 769 1. 644 1.00 0. 793;i0二5. 83;If二0.218;Iwl=l. 798;Iw2二3. 598;Iw=2

2、*lwl+4*lw2;for i=l:69n (i)二(i+ll)*50;Ttq(i)=-19. 313+295. 27*(n(i)/1000) -165. 44*(n(i)/1000)A2+40. 874*(n(i)/1000)A3-3. 8445*(n(i)/1000)M;endfor j=l:5for i二1:69Ft(i, j)=Ttq(i)*ig(j)*iO* nT/r;ua(i, j)=0. 377*r* n( i)/(ig(j)*iO);Fz(i, j)=CdA*ua(i, j)A2/21. 15+mz*g*f;endendplot (ua, Ft, ua, Ff, ua, Ff

3、+Fw)titleC汽車驅動力與行駛阻力平衡圖，)；xlabel ( ua (km/h);ylabelC Ft (N);gtext( Ft 1)gtext( Ft2)gtext ( Ft3)gtext( Ft4)gtext( Ft5)gtext ( Ff+Fw)(2)求最大速度和最大爬坡度for k=l:175nl (k)=3300+k*0. 1;Ttq(k)=-19. 313+295. 27*(nl (k)/1000) -165. 44*(nl (k)/1000)A2+40.874* (n 1 (k)/1000) A33. 8445* (n 1 (k)/1000) A4;Ft(k)=Ttq(

4、k)*ig (5)*i0* nT/r;ua(k)=0. 377*r* nl(k)/(ig (5)*i0);Fz(k)=CdA*ua(k)A2/21. 15+mz*g*f;E(k)=abs(Ft(kFz(k);endfor k=l:175if (E(k) =mi n(E)disp(，汽車最高車速9);disp(ua (k);disp( km/h);endendfor p=l:150n2 (p)=2000+p*0. 5;Ttq (p) =-19. 313+295. 27*(n2 (p)/1000) -165. 44* (n2 (p)/1000) A2+40. 874* (n2 (p)/1000)

5、A3-3.8445*( n2(p)/1000)A4;Ft (p)=Ttq(p)*ig(1)*i0* nT/r;ua(p)=0. 377*r* n2 (p)/ (ig(l)*i0);Fz(p)=CdA*ua(p)A2/21. 15+mz*g*f;af (p)=asi n( (Ft(p)-Fz(p)/(mz*g);endfor p=l:150if(af(p)=max(af)i=ta n(af(p);dispC汽車最大爬坡度二)；disp(i);endend汽車最高車速二99. 0679km/h 汽車最大爬坡度二0. 3518 計算2檔起步加速到70km/h所需時間for i=l:69n(i) =

6、(i+ll)*50;-3.8445*(n(i)/100Ttq(i) =-19. 313+295. 27*(n(i)/1000) -165. 44*(n(i)/1000)A2+40. 874*(n(i)/1000)A30) A4;end for j=l：5for i二1:69deta=l+Iw/(mz*r*2) +If*ig(j)A2*iOA2* nT/(mz*A2) ; ua(i, j)=0. 377*r*n(i)/(ig(j)*iO);a(i, j) = (Ttq(i)*ig(j)*iO*nT/r -CdA*ua(i, j)A2/21. 15-mz*g*f)/(deta*mz);辻(a(i,

7、 j)0. 05)bl(i, j)=a(i, j);ul (i, j)=ua(i, j);elsebl(i, j)=a(i-l, j);ul (i, j)=ua(i -1, j);endb(i, j)=l/bl(i, j);endendxl=ul(:, 1) ;yl=b(:, 1);x2=ul(:, 2) ;y2=b(:, 2);x3=ul(:, 3) ;y3=b(:, 3);x4=ul (:, 4) ;y4=b(:,4);x5=ul(:,5);y5=b(:,5); plot (xl, yl, x2, y2, x3, y3, x4, y4, x5, y5); t itleC 加速度倒數(shù)時間曲線

8、，); axis(0 120 0 30);xlabel ( ua (km/h);ylabel ( l/aj,);gtext (1/al)gtext (J l/a2)gtext (T/a3)gtext C l/a4)gtext C l/a5)for i=l:69A=ua(i, 3)-ua(69, 2);if (A0)J=i；endB=ua(i, 4)-ua(69, 3);if(B0)k二 i;endif (ua(i,4)=70) m=i;endend t=ua(l, 2)*b(l, 2):for pl=2:69tl (pl) = (ua(pl, 2) 一ua(plT, 2)*(b(pl, 2)+

9、b(pl -1, 2)*0. 5; t=t+tl(pl);endfor p2=j:69t2 (p2) = (ua (p2, 3) -ua (p2l, 3) *(b (p2, 3) +b (p2 -1, 3) *0. 5;t=t+t2(p2);endfor p3=k:mt3 (p3) = (ua (p3, 4) -ua (p3l, 4) * (b (p3,4) +b (p3 -1, 4)*0. 5;t二t+t3(p3);endt=t+(ua(j, 3) -ua (69, 2) *b (69, 2) + (ua (k, 4) -ua (69, 3) *b (69, 3)+ (70-ua (m, 4

10、) *b (m, 4);tz=t/3. 6;dispC加速時間二)；disp(tz);disp( s)；加速時間=29. 0585s二、計算與繪制題1中貨車的1)汽車功率平衡圖;2）最高檔與次高檔的等速百公里油耗曲線。己知數(shù)據(jù)略。（參見汽車理論習題第二章第7題）解題程序如下：用Matlab語言ml二2000; m2二1800; mz3880 g9.81r二0.367; CdA二2. 77; f二0.013; nT=0. 85;ig=5. 56 2. 769 1. 644 1.00 0. 793:i0=5. 83; If二0.218; Iwl=l. 798; Iw2=3. 598:nl二815

11、1207 1614 2012 2603 3006 3403 3804;Iw=2*Iwl+4*Iw2;nd=400; Qid二0. 299;for j=l：5for i=l：69n(i) = (i+U)*50;Ttq(i)=-19. 313+295. 27*(n(i)/1000) -165.44*(n(o /1000) A 2+40. 874*(77 (Q /1000) A 3-3. 8445*(n(i)/1000)A4;Pe(i)=n(i)*Ttq (i)/9549;ua(i, j)=0. 377*r*n(i)/(ig(j)*i0);Pz(i, j) = (mz*g*f*ua(i, j)/36

12、00. +CdA*ua(i, j)A3/76140. )/nT;endendplot (ua, Pe, ua, Pz);titleC汽車功率平衡圖)，)；xlabel ( ua (km/h);ylabel ( Pe, Pz (kw);gtext (T)gtext(II) gtext (5 IlT ) gtext(IV) gtext(V) gtext (J P 阻)for j 二 1:5 for i=l ：8 Td(i)=-19. 313+295. 27*(nl (i)/1000. 0) -165. 44*(nl (i)/1000. 0) A2+40. 874*(nl (i)/1000. 0)A

13、3-3. 8445*(nl (i)/1000. 0)A4;Pd(i)=nl (i)*Td(i)/9549;u(i, j) =0. 377*nl(i)*r/(ig(j)*i0);endendb(l)=0. 17768*Pd(l)A4 -5. 8629*Pd(l)A3+72. 379*Pd(l)A2 -416. 46*Pd (1)+1326. 8; b(2)=0. 043072*Pd (2) A4 -2. 0553*Pd(2)A3+36. 657*Pd(2)A2 -303. 98水Pd (2)+1354. 7; b(3)=0. 0068164*Pd(3) A4 -0. 5U84*Pd(3) A3+

14、14. 524*Pd(3)A2 -189. 75*Pd (3)+1284. 4; b(4)=0. 0018555*Pd(4)A4 -0. 18517*Pd(4)A3+7. 0035*Pd(4)A2 -121. 59*Pd(4)+1122. 9;b(5)=0. 00068906*Pd(5)A4 -0. 091077*Pd(5)A3+4. 4763*Pd(5)A2 -98. 893*Pd(5)+1141. 0; b(6)=0. 00035032*Pd(6)A4 -0. 05138*Pd(6)A3+2. 8593*Pd(6)A2 -73. 714水Pd(6)+1051. 2; b(7)=0. 000

15、28230*Pd(7)A4 -0. 047449*Pd(7)A3+2. 9788*Pd(7)A2 -84. 478*Pd(7)+1233. 9; b(8)二-0. 000038568*Pd(8)A40. 00075215*Pd(8)A3+0. 71113*Pd(8)A245. 291*Pd(8) +1129. 7;ul二u(:, 1);u2二u (:,2);u3二u(:, 3)，;u4二u(:,4);u5二u(:,5);Bl二polyfit (ul, b, 3);B2二polyfit (u2, b, 3);B3=polyf it (u3, b, 3);B4=polyfit(u4, b, 3);

16、B5二polyfit (u5, b, 3);for q二1:69bh (q, 1) =polyval (Bl, ua (q, 1);bh (q, 2)=polyval (B2, ua (q, 2);bh (q, 3) =polyval (B3, ua (q, 3);bh(q,4)=polyval(B4, ua(q, 4);bh (q, 5) =polyval (Bo, ua (q, 5);endfor i=l：5for q二1:69Q(q, i) =Pz(q, i)*bh(q, i)/(1. 02*ua(q, i)*7. 05);endend plot (ua (:, 4), Q (:, 4)

17、, ua (:, 5), Q (:, 5); t四檔五檔等速百公里油耗圖，);xlabel (ua(km/h);ylabel ( Qs (L/lOOkm);三、改變1.3題中輕型貨車的主減速器傳動比，做出 i。為5. 17、5.43、5. 83.6. 17.6. 33時的燃油經濟性一加速時間曲線，討論不同值對汽車性能的影響。Mat lab 程序：ml二2000 m2二1800 m二3880 r0=0. 367 gt=0. 85 f=0. 013 CDA=2. 77 i0=5. 83 If=0. 218Iwl=l. 798Iw2=3. 598Ig5=5. 56 2. 769 1. 644 1.0

18、0 0. 793:Ig0=5. 17 5.43 5. 83 6. 17 6. 33;B二1326.8 -416.46 72. 379 -5. 8629 0. 17768;1354. 7 -303. 98 36. 657-2.0353 0. 043072;1284. 4 -189. 75 14. 524n -1101 n ccqoiq-0. 01184 0. 0068164;1122.9-121. 59 7. 00351141.0 -98. 893 4. 4763- 1817 0- 0018oo- 1051. 2 -73. 714 2. 85930. 091077 0. 00068906;123

19、3.9-84. 478 2. 9788-q. 05138 0. 00035032;-0.047449 0.00028230;1129.7 -45. 291 0.71113 -0. 00075215 -0.000038568;n二815 1207 1614 2012 2603 3006 3403 3804;for i=l：5for k二1:8 ua(i,k)=0. 377*0. 367*n(k)/(IgO(i)*Ig5(5);-3. 8445.*(n(i)/Ttq(i)=-19. 313+295. 27. *(n(i)/1000) 一 165. 44. *(n(i)/1000) 92+40. 8

20、74. *(n(i)/1000)93 1000).4;F5(i, k)=0. 013*3880*9. 8+2. 77. *ua(i, k)A2/21. 15;Pe (i, k) =F5(i, k)*ua(i,k)/(3600*0. 85);b5 (i, k)二B (k, 1)+B (k, 2)*Pe (i , k)+B (k, 3)*Pe (i, k) A2+B(k, 4) *Pe (i, k) A3+B (k, 5)*Pe (i, k)A4; endendudl二25; si二50;Fa5=0. 013*3880*9. 8+2. 77. *ual. A2/21. 15;Pe5二Fd5. *u

21、al/(3600*0. 85);dl二polyfit (Pe(l, ：), b5(l, ：), 3);bal=polyval(dl, Pe5);d2=polyfit(Pe(2, :), b5(2, :),3);ba2=polyval(d2, Pe5);d3=polyfit(Pe(3, :), b5(3, :), 3);ba3=polyval(d3, Pe5);d4=polyfit(Pe(4, :), b5(4, :), 3);ba4=polyval(d4, Pe5); d5=polyfit(Pe(5, :), b5(5, :), 3);ba5=polyval(d5, Pe5) ; ba=bal

22、 ba2 ba3 ba4 ba5;Qal=Pe5. *ba*50/(ua1*102*7)ua2=25:40;Q2二1+(2*Iwl+4*Iw2)/(m*rOA2)+If*Ig5(2). A2*IgO(2)A2*gt/(m*rOA2)Fb5=0. 013*3880*9. 8+2. 77. *ua2. A2/21. 15+Q2*m*0. 25;Pb5=Fb5. *ua2/(3600*0. 85);dbl=polyfit (Pe (1, :), b5 (1, :),3);bbl=polyval(dbl, Pb5);db2二polyfit (Pe(2, :), b5(2, :),3);bb2二poly

23、val(db2, Pb5);db3=polyfit(Pe(3, :), b5(3, :),3);bb3=polyval(db3, Pb5);db4二polyfit (Pe (4, :), b5 (4, ：), 3);bb4=polyval(db4, Pb5); db5=polyfit(Pe(5, :),b5(5, :), 3);bb5=polyval(db5, Pb5);bb=bblbb2bb3bb4bb5;Pb二Pb5Pb5Pb5Pb5Pb5;Qb二Pb. *bb/(367. 1*7);for i=l：5for j二 1:15 qb (i, j) =Qb (i, j) +Qb (i, j+1

24、);endendQb2二sum(qb)ud3二40; s2=250;Fc5=0. 013*3880*9. 8+2. 77. *ua3. A2/21. 15;Pc5=Fc5. *ua3/(3600*0. 85);del二polyfit (Pe (1, ：), b5 (1, ：), 3);bcl=polyval(del, Pc5); dc2=polyfit(Pe(2, :), b5(2, :), 3);bc2=polyval(dc2, Pc5);dc3二polyfit (Pe(3, :), b5(3, :),3);bc3=polyval(dc3, Pc5);dc4=polyfit(Pe(4, :)

25、, b5(4,:),3); bc4=polyval(dc4,Pc5);dc5二polyfit(Pe(5,:), b5(5, :), 3); bc5=polyval(dc5, Pc5);bc=bcl bc2 bc3 bc4 bc5;Qc3二Pc5. *bc*250/(ua3*102*7)ua4=40:50;Q4=l + (2*Iwl+4*Iw2) / (m*rOA2) +If*Ig5 (4). A2*IgO (4) A2*gt/(m*rOA2)Fd5二0. 013*3880*9. 8+2. 77. *ua442/21. 15+Q4*m*0. 2; Pd5二Fd5. *ua4/(3600*0. 8

26、5); ddl=polyfit (Pe (1, :), b5 (1, :),3);bdl=polyval(ddl, Pd5);dd2二polyfit (Pe(2, :), b5(2, :),3);bd2=polyval(dd2, Pd5);dd3=polyfit (Pe(3, :), b5(3, :),3);bd3=polyval(dd3, Pd5);dd4=polyfit(Pe(4, :), b5(4, :), 3);bd4=polyval(dd4, Pd5);dd5=polyfit(Pe(5, :), b5(5, :), 3);bd5=polyval(dd5, Pd5);bd二bdlbd2

27、bd3bd4 bd5;Pd二Pd5Pd5Pd5Pd5Pd5;Qd二Pd. *bd/(367. 1*7):for i=l：5for j二 1:10 qd(i, j)=Qd(i, j)+Qd(i, j+1);end end Qd4二sum(qd) ua5二50; s2=250;Ff5=0. 013*3880*9. 8+2. 77. *ua5. A2/21. 15;Pf5=Ff5. *ua5/(3600*0. 85); dfl二polyfit (Pe(1, :), b5(1, :), 3); bfl=polyval (dfl, Pf5); df2二polyfit (Pe(2, :), b5(2, :

28、),3);bf2=polyval (df2, Pf5); df3二polyfit (Pe(3, :), b5(3, :),3); bf3=polyval (df3, Pf5); df4二polyfit (Pe(4, :), b5(4, :), 3);bf4=polyval(df4, Pf5); df5=polyfit(Pe(5, :),b5(5, :), 3); bf5=polyval(df5, Pf5) bf=bfl bf2 bf3 bf4 bf5;Qf5二Pf5. *bf*250/(ua3*102*7) Qi二0.299;Qg二(50-25)/3. 6/0. 36*0. 299;Qg6=Q

29、g Qg Qg Qg QgQ二QalQb2Qc3Qd4Qf5Qg6;Qz=sum(Q)/1075*100 for k二1:5 for i二1:3401;for j二1:5; n(i)二i+599;ua (i, j)=0. 377*r0*n(i) ./(Ig5(j)*IgO (k);Q (j)二 1+ (2*lwl+4*lw2)/(m*rOA2) +lf*lg5 (j). A2*lg0 (k)A2*gt/ (m* (rO)A2);endenduamax=max(ua);uamin=min(ua);ua2=uamin(2):uamax(2);n=ua2*lg0(k)*lg5 (2)/(0. 377

30、*0. 367);Ttq=-19. 313+295. 27*(n/1000) -165. 44*(n/1000).人 2+40. 874*(n/1000)93Ft2=Ttq*lgO(k)*lg5(2)*0. 85/0.367;F2二0. 013*3880*9. 8+2. 77*ua2. A2/21. 15;a2= (Ft2-F2). / (Q (2) *m);t2=trapz (ua2, 1. /a2)/3. 6;ua3=uamax(2):uamax(3);n3=ua3*lg0(k)*lg5(3)/(0. 377*0. 367);Ttq3=-19. 313+295. 27*(n3/1000)

31、-165.44*(n3/1000). A2+40. 874*(n3/1000). A3A4;-3. 8445* (n/1000)44;-3. 8445*(n3/1000).Ft3=Ttq3*lgO (k)*lg5(3)*0.85/0.367;F3=0. 013*3880*9. 8+2. 77*ua3. A2/21. 15;a3=(Ft3-F3)./(Q(3)*m);t3=trapz(ua3, 1. /a3)/3. 6;ua4=uamax(3):70;n4=ua4*lg0(k)*lg5(4)/(0. 377*0. 367);Ttq4二-19. 313+295. 27*(n4/1000) -165

32、. 44*(n4/1000). A2+40. 874*(n4/1000). A3-3. 8445*(n4/1000).A4;Ft4=Ttq4*lgO(k)*lg5(4)*0. 85/0.367;F4二0. 013*3880*9. 8+2. 77*ua4. A2/21. 15;a4= (Ft4-F4). / (Q (4) *m);t4=trapz (ua4, 1. /a4)/3. 6;t (k)二t2+t3+t4;endf二polyfit(Qz, t, 2);QQ二Qz(l):0. 001:Qz(5);b=polyval(f, QQ)plot (QQ, b, Qz, t, *)gtext ( 5

33、. 17)gtext ( 5. 43)gtext ( 5. 83)gtext ( 6. 17)gtext ( 6. 33)titlef燃油經濟性-動力性曲線JxlabelC 燃油經濟性 /L( 100km) -11，)ylabel ( t/s)四、一中型貨車裝有前后制動器分開的雙管路制動系，其有關參 數(shù)如下：載荷質量(kg)質心高hg/m軸距L/m質心至前軸距離a/m制動力分配系數(shù)B空載40800. 8453. 9502. 1000. 38、卄 滿載92901. 1703. 9502.9500. 381) 計算并繪制利用附著系數(shù)曲線和制動效率曲線2) 求行駛車速Ua二30km/h,在二0. 8

34、0路面上車輪不抱死的制動距離。計算時取制動系反應時間2二0. 02s,制動減速度上升時間2二0. 02so3)求制動系前部管路損壞時汽車的制動距離s,制動系后部管路損壞時汽車的制動距離S oMat lab 程序：m 仁 4080;hgl=0. 845;al=2. 100;m2=9290;hg2=l. 17;a2=2. 95;beta=0. 38;L=3. 95;z=0:0. 05:1gf 仁 beta. *z*L. / (L-al+z*hgl);gf2=beta. *z*L. /(La2+z*hg2);grl=(l-beta). *z*L. /(alz*hgl);gr2=(lbeta). *z

35、*L. /(a2z*hg2);薩z;for i二1:21if (z(i)0. 15);g3(i)=z(i)+0. 08;endif (z(i)二0. 3);g3(i)=0. 38+(z(i)-0. 3)/0. 74;endendzl=0. 15:0.01:0.3;g4=zl-0. 08;plot (z, gfl, z, gf2, z, grl, z, gr2, z, g, z, g3, xk axis(0 1 0 1. 2), zl, g4,，x)title(利用附著系數(shù)與制動強度的關系曲線)xlabel (，制動強度z/g，)ylabel (利用附著系數(shù)g ) gtext (，空車前軸)gt

36、ext(空車后軸)gtext (，滿載前軸) gtext (滿載后軸)gtext (，ECE 法規(guī) )C=0:0. 05:1; Erl=(al/L)./(1-beta)+C*hgl/L)*100;Ef=(L-a2)/L. /(beta-C*hg2/L)*100: Er=(a2/L). /(1-beta)+C*hg2/L)*100; plot(C,Er,C,Ef,C,Erl) axis(0 1 0 100) title(前后附著效率曲線)xlabel ( *附著系數(shù)Cf ) ylabel (，制動效率() ) gtext (滿載) gtext ( Ef ) gtext ( Er ) gtext

37、(空載)gtext ( Er)Cl二0. 8El=(akl/L)./(1 -beta)+Cl*hgl/L); E2=(am2/L)/(1 -beta)+Cl*hg2/L);al=El*Cl*9. 8; a2=E2*Cl*9. 8;ud二30;i21二0. 02;i22二0. 02;sl=(i21+i22/2)*ua/3. 6+uaA2/(25. 92*akl); s2=(i21+i22/2)*ua/3. 6+uaA2/(25. 92灼m2);dispC滿載時不抱死的制動距離=，)disp(s2)dispC空載時不抱死的制動距離二，)disp(sl)滿載時不抱死的制動距離=5. 3319空載時不

38、抱死的制動距離=6. 8406beta3=lbeta4=0Ekr=(al/L)/(lbeta4)+Cl*hgl/L);Ekf=(Lal)/(beta3*L -Cl*hgl);Emf=(L-a2)/L./(beta3 -Cl*hg2/L);Emr=(a2/L). /(I beta4)+Cl*hg2/L); akr=0. 8*9. 8*Ekr;akf=0. 8*9. 8*Ekf;amrO. 8*9. 8*Emr:amf二0. 8*9. 8*Emf;skr=(i21+i22/2)*ua/3. 6+uaA2/(25. 92*akr);skf=(i21+i22/2)*ua/3. 6+uaA2/(25.

39、92*akf);smf=(i21+i22/2)*ua/3. 6+u2/(25. 92*amf);smr=(i21+i22/2) *ua/3. 6+uaA2/ (25. 92*amr);dispC空車后管路失效時制動距離Jdisp(skf)dispC空車前管路失效時制動距離，)disp(skr)dispC滿載后管路失效時制動距離，)disp(smf)dispC滿載前管路失效時制動距離)disp(smr)運行結果為：空車后管路失效時制動距離8. 0879空車前管路失效時制動距離10. 0061滿載后管路失效時制動距離13. 5986滿載前管路失效時制動距離7. 5854五、二自由度轎車模型的有關參

40、數(shù)如下：總質量m=1818. 2kg繞Oz軸轉動慣量Iz 3885kg m2軸距L=3. 048m質心至前軸距離a二 1. 463m質心至后軸距離b二 1. 585m前輪總側偏剛度ki=-62618N/rad后輪總側偏剛度2二kiiniQRN/eri轉向系總傳動比1=20試求：2)穩(wěn)態(tài)橫擺角速度增益曲線二茁、車速u=22. 35m/s時的轉向靈敏度二。SW3) 靜態(tài)儲備系數(shù)S. M.,側向加速度為0. 4g時的前、后輪側偏角絕對值之差12與轉彎半徑的比值R/Ro (R0=15m) o4) 車速u二30. 56m/s時，瞬態(tài)響應的橫擺角速度波動的固有(圓)頻率。、阻尼比、反應時間與峰值反應時間M

41、at lab 程序：m二 1818. 2; lz二3885; L=3. 048 ;a=l. 463;b=l. 585;kl = -62618;k2二-110185;i二20;g=9.8;R0=15;ul=30. 56;K=m/LA2* (a/k2 -b/kl);uch二 1/KU1/2);u二0:0. 05:25;s二u/L./(l+K*u. A2);disp (，穩(wěn)定因素K);disp(K);dispC特征車速uch=，);disp(uch);plot(u,s);xlabel ( ua/ (m/s);ylabelC穩(wěn)態(tài)橫擺增益，)；titlef汽車的穩(wěn)態(tài)橫態(tài)擺角速度增益曲線，)；dispCu

42、a二22.35m/s時，轉向靈敏度為)；disp (s(448);SM=k2/(kl+k2)-a/L;A二K*0. 4*g*L:Q二L/RO;R=L/(Q-A);r二R/RO;disp(靜態(tài)儲備系數(shù)S.M二)；disp(SM);dispC前后輪側偏角絕對值Z差(dl-a2)3 );disp (A);disp (，轉彎半徑的比值J;disp(r);wO二L/ul*(kl*k2*(l+K*uU2)/(m*lz)A(l/2);n= (m*(aA2*kl+bA2*k2) -Iz*(kl+k2)/(2*L*(m*Iz*kl*k2*(1+K*ul*2)A(1/2);t=atan(l 41 人 2)人,(l

43、/2)/( -m*ul*a*w0/(L*k2) -n)/(w0*(l -nA2)A(l/2);e=ata n (1-門人 2)人(2)/n)/(w0*(l -n 人 2)人(1/2)+匕dispC橫擺角速度波動時的固有頻率為，)；disp (wO);dispC阻尼比為);disp( n);disp(反應時間為，)；disp(t);disp(，達到第一峰值的時間為，);disp(e);運算結果：穩(wěn)定因素K=0. 0024特征車速uch=20. 6053ua=22. 35m/s時，轉向靈敏度為3. 3690靜態(tài)儲備系數(shù)S.M.=0. 1576前后輪側偏角絕對值之差(al-a2)=0. 0281轉彎

44、半徑的比值二1. 1608橫擺角速度波動時的固有頻率為5. 5758阻尼比為0. 5892反應時間為0. 1811達到第一峰值的時間為0. 3899六、車身-車輪雙質量系統(tǒng)參數(shù)：f。1.5Hz,0. 25,9,10o“人體-座椅”系統(tǒng)參數(shù)：fs 3Hz, s0. 25。車速u 20m/s,路 面不平度系數(shù)Gqn。 2. 56 10 m3,參考空間頻率 n- =0. Im% 計算時頻率步長f 0. 2Hz,計算頻率點數(shù)N 180o1)計算并畫出幅頻特性乙/q、Zz/乙、q/z2和均方根值譜.Gz, r、Gr、Gaf譜圖。進一步計算0. 125 0. 5 參數(shù)：2、)改變“人體-座椅”系統(tǒng)參數(shù):f

45、sl56Hz,s 分析aw、Law值隨fs、s的變化。3)分別改變車身-車輪雙質量系統(tǒng)f0253Hz, 0. 1250. 5,4. 51 &5 20。繪制Z八三個響應量均方根值隨以上四個系統(tǒng)參數(shù)變化的曲線。Mat lab 程序：f0=l. 5;g=0. 25:r=9;u=10;fs=3;gs=0. 25;v二20;Gqn0二2. 56e -008;n0=0. 1;df二0. 2;N二180;f二0:0.2:36;ff0=f/f0;d=(l -ff0/2).*(l+r -l/u*ff0. a2 -1) .A2+4*gA2*ffO. a2. *r -(l/u+l)*ff0. A2.A2; zlq=

46、r*sqrt (1 -ffO. A2). A2+4*gA2*ffO.A2). /d);z2z 仁 sqrt(1+(2*g*ff0). A2). /(1 -ff0. A2). A2+(2*g*ff0). A2);ffs=f/fs;pz2=sqrt(l+(2*gs*ffs). A2). /(1 -ffs. A2). A2+(2*gs*ffs). A2)loglog(f, zlq, f, z2zl,., f, pz2);axis(0 100 0 10)set(gca, xtick, 0 0. 1 1 10 100)set(gca, ytick, 0 0. 1 1 101)xlabel (激振頻率 f

47、/HZjw=2*pi*f;Gzl=w. A2. *zlq. * (2*3. 14. /w). *sqrt (Gq n0*n 0 人 2 憐)；loglog(f, Gzl)hold onGz2=w. A2. *zlq. *z2zl. *(2*3. 14. /w). *sqrt (Gq n0*n 0 丿J 2;title (，車輪部分加速度均方根值,)xlabel (，激振頻率 f/Hzf)ylabel ( Gzl )loglog(f, Gz2)titlef車身部分加速度均方根值,)xlabel (激振頻率 f/Hz?)ylabel ( Gz2)Gza=w. A2. *zlq. *z2zl. *p

48、z2. * (2*3. 14. /w). *sqrt (Gq n0*n 0 人 20) ; loglog (f, Gza)title (，人體加速度均方根值,)xlabel (激振頻率 f/Hz)ylabelC i Gza )format longf0=1.5;g=0.25;r=9;u=10;fs=3;gs=0.25;v=20;Gq | 0二2. 56*10 人(-8) ;n0二0. 1; df二0. 2;N二 180;f二0. 000001:0.2:36;ffO=f/fO;ffs=f/fs;w=2*pi*f;d二(1-ffO. a2) . * (1+r -ffO. a2/u) -1). A2

49、+4. *gA2. *ffO. A2. * (r -(1/u+l). *ffO. A2). A2;zlq=r*sqrt(1 -ff0. A2). A2+4*gA2*ffO. A2). /d);z2zl=sqrt(l+(2*g*ff0). A2). / (1 -ffO. A2). A2+(2*g*ff0). A2);ffs=f/fs;pz2=sqrt(l+(2*gs*ffs). A2). /(1 -ffs. A2). A2+(2*gs*ffs). A2);ez二sqrt (sum(16. *(3. 14) A4. *f. A2. *Gqn0. *n0A2. *v*df);dispf路面不平度加速

50、度均方根值三)disp(ez)ezl=sqrt (sum(zlq. A2. *16. *piA4. *Gqn0. *v. *f. A2*df*n0A2);dispf車輪加速度均方根值9)disp(ezl) ez2=sqrt (sum(z2zl. A2. *zlq. A2. *16. *piA4. *Gqn0. *v. *f. A2*nOA2*df) ; disp(車身力口速度 均方根值三)disp(ez2) Gza=w. A2. *zlq. *z2zl. *pz2. *(2*3. 14. /w). *sqrt (GqnO*nOA2*v) ; eza=sqrt (sum(Gza. A2*df)

51、dispe人體加速度均方根值二)disp(eza)Gza=w. A2. *zlq. *z2zl. *pz2. *(2*3. 14. /w). *sqrt(Gqn0*n0A2*v);for i=l：180if(f (i)0. 5)fl(i)=f(i);Gza(i)=w(i). A2. *zlq(i). *z2zl (i). *pz2(i). *(2*3. 14. /w(i). *sqrt (GqnO*nOA2*v);al(i)=0. 5*Gza(i).A2*0. 2;endif(f(i)2&f(i)4&f(i)12. 5&f(i) fsl=linspace(l. 5, 6, 100);gs2=linspace(0.125, 0. 5, 10);for j二1:100f0=1.5;g=0.25;r=9;u=10;fs=3;gs=0.25;v=20;Gq | 0二2. 56*10 人(-8) ;n0二0. 1; df二0. 2;N二 180;f=0. 000001:0.2:36;ff0二f/f0;ffsl=f/fsl(j);w=2*pi*f;d=(l-ffO. A2).*(l+r -ffO.*2/u) -l).A2+4. *gA2. *ff0?2.*(r -(1/u+l). *ffO. A2). A2;zlq=