汽車理論編程題

1、第一章1）m1=2000 %裝載質量 m2=1800 %整車整備質量 m=3880 ； %總質量 r0=0.367 ； %車輪半徑 gt=0.85 ； % 傳動系機械效率 f=0.013 ； % 波動阻力系數(shù) CDA=2.77 ； % 空氣阻力系數(shù)×迎風面積 i0=5.83 ； %主減速器傳動比 If=0.218 ； %飛輪轉功慣量 Iw1=1.798 ； % 二前輪轉動慣量 Iw2=3.598 ； % 四后輪轉功慣量 Ig5=5.56 2.769 1.644 1.00 0.793；%變速器傳動比%繪制功率平衡圖，以五檔為例for i=1:3401; n(i)=i+599; Ttq(

2、i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).2+40.874*(n(i)/1000).3-3.8445*(n(i)/1000).4;%轉矩擬合曲線end for i=1:3401; for j=1:5; Ft(i,j)=Ttq(i)*i0*Ig5(j)*gt/r0;%計算驅動力 ua(i,j)=0.377*r0*n(i)./(Ig5(j)*i0);%計算車速 F(i,j)=f*m*9.8+CDA*ua(i,j).2/21.15;%計算行駛阻力 end end plot(ua,Ft,ua,F)xlabel('ua/(km/h)'

3、;);%對 汽車驅動力-行駛阻力平衡圖 進行標注 ylabel('F/N');title('汽車驅動力-行駛阻力平衡圖');gtext('Ft1')gtext('Ft2')gtext('Ft3')gtext('Ft4')gtext('Ft5')gtext('Ff+Fw')gridaxis(90 100 1500 2000)ginput(2)由作圖法可求出交點 ans = 1.0e+003 * 0.0991 1.78100.0991 1.7791用ginput求最大速度

4、為99.1km/h2）用符號函數(shù)求最大速度syms n ss='(-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/1000)4)*5.83*0.793*0.85/0.367=0.013*3880*9.8+2.77*(0.377*0.367*n/(5.83*0.793)2/21.15'a=solve(s,n);b=double(a)ua=0.377*0.367*b/(5.83*0.793)b = 1.0e+003 * 0.2597 3.3108 3.5307 - 2.7101i 3.5307

5、+ 2.7101iua = 1.0e+002 * 0.0777 0.9908 1.0566 - 0.8111i 1.0566 + 0.8111i所以最大速度為90.08km/h%計算最大爬坡度D=(Ft-CDA*ua.2/21.15)/(m*9.8);C=(D-f*sqrt(1-D.2+f2)/(1+f2)a=asin(C);i=tan(a);imax=max(i)運行結果為 imax =所以最大爬坡度為0.353%計算附著率L=3.2hg=0.9a=1.947imax=0.353Ca2=imax/(a/L+imax*hg/L)運行結果為 Ca2 = 0.49878571113172for j

6、=1:5for i=1:3401Q(j)=1+(2*Iw1+4*Iw2)/(m*r02)+If*Ig5(j).2*i02*gt/(m*(r0)2);%汽車旋轉質量換算系數(shù) a(i,j)=(Ft(i,j)-F(i,j)/(Q(j)*m); if (a(i,j)<0.005) a(i,j)=a(i,j-1); end endendplot(ua(:,1),1./a(:,1),ua(:,2),1./a(:,2),ua(:,3),1./a(:,3),ua(:,4),1./a(:,4),ua(:,5),1./a(:,5)axis(0 98 0 10)title('加速度倒數(shù)時間曲線'

7、;)xlabel('ua(km/h)')ylabel('1/a(s2/m2)')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext('1/a5')3）%計算從二檔起步加速到70km/h需時間for i=1:3401; for j=1:5; n(i)=i+599; ua(i,j)=0.377*r0*n(i)./(Ig5(j)*i0);%計算車速Q(mào)(j)=1+(2*Iw1+4*Iw2)/(m*r02)+If*Ig5(j).2

8、*i02*gt/(m*(r0)2);%汽車旋轉質量換算系數(shù) end enduamax=max(ua)uamin=min(ua)ua2=uamin(2):uamax(2);n=ua2*5.83*Ig5(2)/(0.377*0.367);Ttq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;%轉矩擬合曲線Ft2=Ttq*5.83*Ig5(2)*0.85/0.367;F2=0.013*3880*9.8+2.77*ua2.2/21.15;a2=(Ft2-F2)./(Q(2)*m);t2=tra

9、pz(ua2,1./a2)/3.6ua3=uamax(2):uamax(3);n3=ua3*5.83*Ig5(3)/(0.377*0.367);Ttq3=-19.313+295.27*(n3/1000)-165.44*(n3/1000).2+40.874*(n3/1000).3-3.8445*(n3/1000).4;%轉矩擬合曲線Ft3=Ttq3*5.83*Ig5(3)*0.85/0.367;F3=0.013*3880*9.8+2.77*ua3.2/21.15;a3=(Ft3-F3)./(Q(3)*m);t3=trapz(ua3,1./a3)/3.6ua4=uamax(3):70n4=ua4*

10、5.83*Ig5(4)/(0.377*0.367);Ttq4=-19.313+295.27*(n4/1000)-165.44*(n4/1000).2+40.874*(n4/1000).3-3.8445*(n4/1000).4;%轉矩擬合曲線Ft4=Ttq4*5.83*Ig5(4)*0.85/0.367;F4=0.013*3880*9.8+2.77*ua4.2/21.15;a4=(Ft4-F4)./(Q(4)*m);t4=trapz(ua4,1./a4)/3.6t=t2+t3+t4%加速到70km/h的時間運行結果為：t2 =6.2033t3 = 8.5885t4 = 11.1905t = 25

11、.9823所以從二檔起步加速到70km/h需25.9823s。第二章m1=2000 ; %裝載質量 m2=1800 ; %整車整備質量 m=3880 ; %總質量 r0=0.367 ;%車輪半徑 gt=0.85 ; %傳動系機械效率 f=0.013 ; %波動阻力系數(shù) CDA=2.77 ; % 空氣阻力系數(shù)×迎風面積 i0=5.83 ;%主減速器傳動比 If=0.218 ;%飛輪轉功慣量 Iw1=1.798 ;% 二前輪轉動慣量 Iw2=3.598 ;% 四后輪轉功慣量 Ig5=5.56 2.769 1.644 1.00 0.793;%變速器傳動比%繪制功率平衡圖，以五檔為例for

12、i=1:3401; n(i)=i+599; Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).2+40.874*(n(i)/1000).3-3.8445*(n(i)/1000).4;%轉矩擬合曲線end for i=1:3401; for j=1:5; Ft(i,j)=Ttq(i)*i0*Ig5(j)*gt/r0;%計算驅動力 ua(i,j)=0.377*r0*n(i)./(Ig5(j)*i0);%計算車速 F(i,j)=f*m*9.8+CDA*ua(i,j).2/21.15;%計算行駛阻力 Pe(i,j)=Ft(i,j)*ua(i,j)

13、/(3600*gt); P(i,j)=F(i,j)*ua(i,j)/(3600*gt); end endplot(ua,Pe,ua,P)title('汽車功率平衡圖')ylabel('Pe/kw')xlabel('ua/(km/h)')gtext('I')gtext('II')gtext('III')gtext('IV')gtext('V')gtext('P阻')format long B=1326.8 -416.46 72.379 -5.8629 0

14、.17768; 1354.7 -303.98 36.657 -2.0553 0.043072; 1284.4 -189.75 14.524 -0.51184 0.0068164; 1122.9 -121.59 7.0035 -0.18517 0.0018555; 1141.0 -98.893 4.4763 -0.091077 0.00068906; 1051.2 -73.714 2.8593 -0.05138 0.00035032; 1233.9 -84.478 2.9788 -0.047449 0.00028230; 1129.7 -45.291 0.71113 -0.00075215 -0

15、.000038568;n=815 1207 1614 2012 2603 3006 3403 3804;ua4=0.377*0.367*n./(5.83*Ig5(4)%計算車速ua5=0.377*0.367*n./(5.83*Ig5(5)%計算車速Ttq=-19.313+295.27.*(n/1000)-165.44.*(n/1000).2+40.874.*(n/1000).3-3.8445.*(n/1000).4;%轉矩擬合曲線F4=0.013*3880*9.8+2.77.*ua4.2/21.15;%計算行駛阻力F5=0.013*3880*9.8+2.77.*ua5.2/21.15;Pe4=

16、F4.*ua4/(3600*0.85) %阻力功率Pe5=F5.*ua5/(3600*0.85)for i=1:8 b4(i)=B(i,1)+B(i,2)*Pe4(i)+B(i,3)*Pe4(i)2+B(i,4)*Pe4(i)3+B(i,5)*Pe4(i)4; b5(i)=B(i,1)+B(i,2)*Pe5(i)+B(i,3)*Pe5(i)2+B(i,4)*Pe5(i)3+B(i,5)*Pe5(i)4;endQs4=b4.*Pe4./(ua4*7*1.02)Qs5=b5.*Pe5./(ua5*7*1.02)Ua5=ua5(find(ua5<110)QS5=Qs5(find(Qs5>

17、7)f4=polyfit(ua4,Qs4,3);f5=polyfit(Ua5,QS5,3);for i=1:3301 n(i)=i+599; UA4(i)=0.377*0.367*n(i)/(5.83*Ig5(4); UA5(i)=0.377*0.367*n(i)/(5.83*Ig5(5);endQ4=polyval(f4,UA4);Q5=polyval(f5,UA5);plot(UA4,Q4,UA5,Q5)gtext('四檔')gtext('五檔')title('四檔和五檔等速百公里油耗曲線')xlabel('ua/(km/h)'

18、;)ylabel('Q/(L/100km)')第三章m1=2000 ; %裝載質量 m2=1800 ; %整車整備質量 m=3880 ; %總質量 r0=0.367 ;%車輪半徑 gt=0.85 ; %傳動系機械效率 f=0.013 ; %波動阻力系數(shù) CDA=2.77 ; % 空氣阻力系數(shù)×迎風面積 i0=5.83 ;%主減速器傳動比 If=0.218 ;%飛輪轉功慣量 Iw1=1.798 ;% 二前輪轉動慣量 Iw2=3.598 ;% 四后輪轉功慣量 Ig5=5.56 2.769 1.644 1.00 0.793;%變速器傳動比Ig0=5.17 5.43 5.83

19、 6.17 6.33;%主減速器傳動比 B=1326.8 -416.46 72.379 -5.8629 0.17768; 1354.7 -303.98 36.657 -2.0553 0.043072; 1284.4 -189.75 14.524 -0.51184 0.0068164; 1122.9 -121.59 7.0035 -0.18517 0.0018555; 1141.0 -98.893 4.4763 -0.091077 0.00068906; 1051.2 -73.714 2.8593 -0.05138 0.00035032; 1233.9 -84.478 2.9788 -0.047

20、449 0.00028230; 1129.7 -45.291 0.71113 -0.00075215 -0.000038568;n=815 1207 1614 2012 2603 3006 3403 3804;for i=1:5 for k=1:8ua(i,k)=0.377*0.367*n(k)/(Ig0(i)*Ig5(5);%計算車速Ttq(i)=-19.313+295.27.*(n(i)/1000)-165.44.*(n(i)/1000).2+40.874.*(n(i)/1000).3-3.8445.*(n(i)/1000).4;%轉矩擬合曲線F5(i,k)=0.013*3880*9.8+

21、2.77.*ua(i,k)2/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)2+B(k,4)*Pe(i,k)3+B(k,5)*Pe(i,k)4; endend %六工況循環(huán)行駛的百公里油耗%第一工況 ua1=25km/h s1=50mua1=25; s1=50;Fa5=0.013*3880*9.8+2.77.*ua1.2/21.15;Pe5=Fa5.*ua1/(3600*0.85);d1=polyfit(Pe(1,:),b5(1,:),3);ba1=polyval(d

22、1,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=ba1 ba2 ba3 ba4 ba5;Qa1=Pe5.*ba*50/(ua1*102*7)%第二工況 勻加速度為0.25m/s.2ua2=25:40;Q2=1+(2*Iw1+4

23、*Iw2)/(m*r02)+If*Ig5(2).2*Ig0(2)2*gt/(m*r02)Fb5=0.013*3880*9.8+2.77.*ua2.2/21.15+Q2*m*0.25;Pb5=Fb5.*ua2/(3600*0.85);db1=polyfit(Pe(1,:),b5(1,:),3);bb1=polyval(db1,Pb5);db2=polyfit(Pe(2,:),b5(2,:),3);bb2=polyval(db2,Pb5);db3=polyfit(Pe(3,:),b5(3,:),3);bb3=polyval(db3,Pb5);db4=polyfit(Pe(4,:),b5(4,:),

24、3);bb4=polyval(db4,Pb5);db5=polyfit(Pe(5,:),b5(5,:),3);bb5=polyval(db5,Pb5);bb=bb1 bb2 bb3 bb4 bb5;Pb=Pb5 Pb5 Pb5 Pb5 Pb5;Qb=Pb.*bb/(367.1*7);for i=1:5 for j=1:15 qb(i,j)=Qb(i,j)+Qb(i,j+1); end end Qb2=sum(qb')%第三工況 勻速40km/hua3=40; s2=250;Fc5=0.013*3880*9.8+2.77.*ua3.2/21.15;Pc5=Fc5.*ua3/(3600*0

25、.85);dc1=polyfit(Pe(1,:),b5(1,:),3);bc1=polyval(dc1,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,:),b5(4,:),3);bc4=polyval(dc4,Pc5);dc5=polyfit(Pe(5,:),b5(5,:),3);bc5=polyval(dc5,Pc5) ;bc=bc1 bc2 bc3 bc4 bc5;Qc3=Pc5.*bc

26、*250/(ua3*102*7)%第四工況 勻加速度0.2m/s.2ua4=40:50;Q4=1+(2*Iw1+4*Iw2)/(m*r02)+If*Ig5(4).2*Ig0(4)2*gt/(m*r02)Fd5=0.013*3880*9.8+2.77.*ua4.2/21.15+Q4*m*0.2;Pd5=Fd5.*ua4/(3600*0.85);dd1=polyfit(Pe(1,:),b5(1,:),3);bd1=polyval(dd1,Pd5);dd2=polyfit(Pe(2,:),b5(2,:),3);bd2=polyval(dd2,Pd5);dd3=polyfit(Pe(3,:),b5(3

27、,:),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=bd1 bd2 bd3 bd4 bd5;Pd=Pd5 Pd5 Pd5 Pd5 Pd5;Qd=Pd.*bd/(367.1*7);for i=1:5 for j=1:10 qd(i,j)=Qd(i,j)+Qd(i,j+1); end end Qd4=sum(qd')%第五工況 ua=50km/h s=250ua5=50; s

28、2=250;Ff5=0.013*3880*9.8+2.77.*ua5.2/21.15;Pf5=Ff5.*ua5/(3600*0.85);df1=polyfit(Pe(1,:),b5(1,:),3);bf1=polyval(df1,Pf5);df2=polyfit(Pe(2,:),b5(2,:),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

29、,:),3);bf5=polyval(df5,Pf5) ;bf=bf1 bf2 bf3 bf4 bf5;Qf5=Pf5.*bf*250/(ua3*102*7)%第六工況 減速 減速度為0.36m/s.2Qi=0.299;Qg=(50-25)/3.6/0.36*0.299;Qg6=Qg Qg Qg Qg QgQ=Qa1 Qb2 Qc3 Qd4 Qf5 Qg6;Qz=sum(Q)/1075*100%計算從二檔起步加速到70km/h需時間for k=1:5for i=1:3401; for j=1:5; n(i)=i+599; ua(i,j)=0.377*r0*n(i)./(Ig5(j)*Ig0(k

30、);%計算車速Q(mào)(j)=1+(2*Iw1+4*Iw2)/(m*r02)+If*Ig5(j).2*Ig0(k)2*gt/(m*(r0)2);%汽車旋轉質量換算系數(shù) endenduamax=max(ua);uamin=min(ua);ua2=uamin(2):uamax(2);n=ua2*Ig0(k)*Ig5(2)/(0.377*0.367);Ttq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;%轉矩擬合曲線Ft2=Ttq*Ig0(k)*Ig5(2)*0.85/0.367;F2=0.01

31、3*3880*9.8+2.77*ua2.2/21.15;a2=(Ft2-F2)./(Q(2)*m);t2=trapz(ua2,1./a2)/3.6;ua3=uamax(2):uamax(3);n3=ua3*Ig0(k)*Ig5(3)/(0.377*0.367);Ttq3=-19.313+295.27*(n3/1000)-165.44*(n3/1000).2+40.874*(n3/1000).3-3.8445*(n3/1000).4;%轉矩擬合曲線Ft3=Ttq3*Ig0(k)*Ig5(3)*0.85/0.367;F3=0.013*3880*9.8+2.77*ua3.2/21.15;a3=(Ft

32、3-F3)./(Q(3)*m);t3=trapz(ua3,1./a3)/3.6;ua4=uamax(3):70;n4=ua4*Ig0(k)*Ig5(4)/(0.377*0.367);Ttq4=-19.313+295.27*(n4/1000)-165.44*(n4/1000).2+40.874*(n4/1000).3-3.8445*(n4/1000).4;%轉矩擬合曲線Ft4=Ttq4*Ig0(k)*Ig5(4)*0.85/0.367;F4=0.013*3880*9.8+2.77*ua4.2/21.15;a4=(Ft4-F4)./(Q(4)*m);t4=trapz(ua4,1./a4)/3.6;

33、t(k)=t2+t3+t4;end f=polyfit(Qz,t,2);QQ=Qz(1):0.001:Qz(5);b=polyval(f,QQ)plot(QQ,b,Qz,t,'*') gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')title('燃油經(jīng)濟性-動力性曲線')xlabel('燃油經(jīng)濟性/L(100km)-1')ylabel('t/s')第四章m1=4080;h

34、g1=0.845;a1=2.100; %空載時m2=9290;hg2=1.17;a2=2.95; %滿載beta=0.38;L=3.95;z=0:0.05:1gf1=beta.*z*L./(L-a1+z*hg1);gf2=beta.*z*L./(L-a2+z*hg2);gr1=(1-beta).*z*L./(a1-z*hg1);gr2=(1-beta).*z*L./(a2-z*hg2);g=z;for i=1:21%ECE制動法規(guī) if (z(i)<0.3&z(i)>0.15); g3(i)=z(i)+0.08; end if(z(i)>=0.3); g3(i)=0.

35、38+(z(i)-0.3)/0.74; end end z1=0.15:0.01:0.3; g4=z1-0.08; plot(z,gf1,'-.',z,gf2,z,gr1,'-.',z,gr2,z,g,z,g3,'xk',z1,g4,'x')axis(0 1 0 1.2)title('利用附著系數(shù)與制動強度的關系曲線')xlabel('制動強度z/g')ylabel('利用附著系數(shù)g')gtext('空車前軸')gtext('空車后軸')gtext(&

36、#39;滿載前軸')gtext('滿載后軸')gtext('ECE法規(guī)')C=0:0.05:1;Er1=(a1/L)./(1-beta)+C*hg1/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,Er1)axis(0 1 0 100)title('前后附著效率曲線')xlabel('附著系數(shù)C')ylabel('制動效率（%）')gtext('滿載')gtext

37、('Ef')gtext('Er')gtext('空載')gtext('Er')當C=0.8時 前軸先抱死%計算附著率C=0.8時，ua=30km/h時制動距離C1=0.8E1=(ak1/L)./(1-beta)+C1*hg1/L);E2=(am2/L)/(1-beta)+C1*hg2/L);a1=E1*C1*9.8;a2=E2*C1*9.8;ua=30;i21=0.02;i22=0.02;s1=(i21+i22/2)*ua/3.6+ua2/(25.92*ak1);s2=(i21+i22/2)*ua/3.6+ua2/(25.92*a

38、m2);disp('滿載時不抱死的制動距離=')disp(s2)disp('空載時不抱死的制動距離=')disp(s1)滿載時不抱死的制動距離= 5.3319空載時不抱死的制動距離=6.8406%前后管路損壞時汽車的制動距離beta3=1 %后管路損壞beta4=0 %前管路損壞Ekr=(a1/L)/(1-beta4)+C1*hg1/L);Ekf=(L-a1)/(beta3*L-C1*hg1);Emf=(L-a2)/L./(beta3-C1*hg2/L);Emr=(a2/L)./(1-beta4)+C1*hg2/L);akr=0.8*9.8*Ekr;akf=0.

39、8*9.8*Ekf;amr=0.8*9.8*Emr;amf=0.8*9.8*Emf;skr=(i21+i22/2)*ua/3.6+ua2/(25.92*akr);skf=(i21+i22/2)*ua/3.6+ua2/(25.92*akf);smf=(i21+i22/2)*ua/3.6+ua2/(25.92*amf);smr=(i21+i22/2)*ua/3.6+ua2/(25.92*amr);disp('空車后管路失效時制動距離')disp(skf)disp('空車前管路失效時制動距離')disp(skr)disp('滿載后管路失效時制動距離')

40、disp(smf)disp('滿載前管路失效時制動距離')disp(smr)運行結果為：空車后管路失效時制動距離 8.0879空車前管路失效時制動距離 10.0061滿載后管路失效時制動距離 13.5986滿載前管路失效時制動距離 7.5854第五章、m=1818.2;Iz=3885;L=3.048;a=1.463;b=1.585;k1=-62618;k2=-110185;i=20;%求穩(wěn)定性因素和特征車速K=m*(a/k2-b/k1)/L2;Uch=sqrt(1/K);disp('穩(wěn)定性因素 K=');disp(K);disp('特征車速Uch=

41、9;);disp(Uch);運行結果為：穩(wěn)定性因素 K= 0.0024特征車速Uch= 20.6053%繪制穩(wěn)態(tài)橫擺角速度增益曲線u=0:0.1:100;wr=(u/L)./(1+K*u.2);%轉向不足時穩(wěn)態(tài)橫擺角速度增益K1=0;wr1=(u/L)./(1+K1*u.2);%中性轉向時穩(wěn)態(tài)橫擺角速度增益plot(u,wr,u,wr1)axis(0 100 0 20);title('穩(wěn)態(tài)橫擺角速度增益曲線') xlabel('u(km/h)')ylabel('橫擺角速度增益')gtext('K=0')gtext('K=0

42、.0024')%計算當速度為22.35m/s時的轉向靈敏度ua=22.35*3.6;Wr=(ua/L)/(1+K*ua2);disp('速度為22.35m/s時的轉向靈敏度Wr=')disp(Wr)運行結果為：速度為22.35m/s時的轉向靈敏度Wr=1.6247%計算靜態(tài)儲備系數(shù)S.M.及前后輪側偏角絕對值之差（a1-a2），和轉彎半徑之比R/R0ay=0.4*9.8;R=15;SM=k2/(k1+k2)-a/L;f=K*ay*L;b=L/（L-R*f）;disp('靜態(tài)儲備系數(shù)S.M.=')disp(SM)disp('前后輪側偏角絕對值之差（

43、a1-a2）=')disp(f)disp('轉彎半徑之比R/R0')disp(b)運行結果為：靜態(tài)儲備系數(shù)S.M.= 0.1576前后輪側偏角絕對值之差（a1-a2）= 0.0281轉彎半徑之比R/R0 1.1608%車速Ua=30.56m/s,瞬態(tài)響應的橫擺角速度波動的固有(圓)頻率w0、阻尼比g、反應時間I 與峰值反應時間 tUa=30.56*3.6;w0=L/Ua*sqrt(k1*k2/(m*Iz)*(1+K*Ua2);d=(-m)*(a2*k1+b2*k2)-Iz*(k1+k2);c=2*L*sqrt(m*Iz*k1*k2*(1+K*Ua2);g=d/c;I=a

44、tan(sqrt(1-g2)/(-m*Ua*a*w0/(L*k2)-g)/(w0*sqrt(1-g2);t=atan(sqrt(1-g2)/g)/(w0*(1-g2)+I;disp('車速Ua=30.56m/s,瞬態(tài)響應的橫擺角速度波動的固有(圓)頻率w0')disp(w0)disp('阻尼比g')disp(g)disp('反應時間I')disp(I)disp('峰值反應時間 t')disp(t)運行結果為：車速Ua=30.56m/s,瞬態(tài)響應的橫擺角速度波動的固有(圓)頻率w0 4.7035阻尼比g 0.1621反應時間I 0.

45、0529峰值反應時間 t0.3603第六章 f0=1.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=(1-ff0.2).*(1+r-1/u*ff0.2-1).2+4*g2*ff0.2.*r-(1/u+1)*ff0.2.2;z1q=r*sqrt(1-ff0.2).2+4*g2*ff0.2)./d);z2z1=sqrt(1+(2*g*ff0).2)./(1-ff0.2).2+(2*g*ff0).2);ffs=f/fs;pz2=sqrt(1+(2*gs*ffs)

46、.2)./(1-ffs.2).2+(2*gs*ffs).2)loglog(f,z1q,'*',f,z2z1,'.',f,pz2);axis(0 100 0 10)set(gca,'xtick',0 0.1 1 10 100)set(gca,'ytick',0 0.1 1 10)xlabel('激振頻率f/HZ') w=2*pi*f;Gz1=w.2.*z1q.*(2*3.14./w).*sqrt(Gqn0*n02*v);loglog(f,Gz1)hold onGz2=w.2.*z1q.*z2z1.*(2*3.14./

47、w).*sqrt(Gqn0*n02*v);title('車輪部分加速度均方根值')xlabel('激振頻率f/Hz')ylabel('|Gz1|')loglog(f,Gz2)title('車身部分加速度均方根值')xlabel('激振頻率f/Hz')ylabel('|Gz2|')Gza=w.2.*z1q.*z2z1.*pz2.*(2*3.14./w).*sqrt(Gqn0*n02*v);loglog(f,Gza)title('人體加速度均方根值')xlabel('激振頻率f/

48、Hz')ylabel('|Gza|')format longf0=1.5;g=0.25;r=9;u=10;fs=3;gs=0.25;v=20;Gqn0=2.56*10(-8);n0=0.1;df=0.2;N=180;f=0.000001:0.2:36;ff0=f/f0;ffs=f/fs;w=2*pi*f;d=(1-ff0.2).*(1+r-ff0.2/u)-1).2+4.*g2.*ff0.2.*(r-(1/u+1).*ff0.2).2;z1q=r*sqrt(1-ff0.2).2+4*g2*ff0.2)./d);z2z1=sqrt(1+(2*g*ff0).2)./(1-f

49、f0.2).2+(2*g*ff0).2);ffs=f/fs;pz2=sqrt(1+(2*gs*ffs).2)./(1-ffs.2).2+(2*gs*ffs).2);ez=sqrt(sum(16.*(3.14)4.*f.2.*Gqn0.*n02.*v*df);disp('路面不平度加速度均方根值=')disp(ez)ez1=sqrt(sum(z1q.2.*16.*pi4.*Gqn0.*v.*f.2*df*n02);disp('車輪加速度均方根值=')disp(ez1)ez2=sqrt(sum(z2z1.2.*z1q.2.*16.*pi4.*Gqn0.*v.*f.2

50、*n02*df);disp('車身加速度均方根值=')disp(ez2)Gza=w.2.*z1q.*z2z1.*pz2.*(2*3.14./w).*sqrt(Gqn0*n02*v);eza=sqrt(sum(Gza.2*df)disp('人體加速度均方根值=')disp(eza)Gza=w.2.*z1q.*z2z1.*pz2.*(2*3.14./w).*sqrt(Gqn0*n02*v);for i=1:180 if(f(i)<2&f(i)>0.5) f1(i)=f(i);Gza(i)=w(i).2.*z1q(i).*z2z1(i).*pz2(

51、i).*(2*3.14./w(i).*sqrt(Gqn0*n02*v); a1(i)=0.5*Gza(i).2*0.2; end if(f(i)>2&f(i)<4) a2(i)=f(i)/4*Gza(i).2; f2(i)=f(i); endif(f(i)>4&f(i)<12.5); a3(i)=Gza(i).2; f3(i)=f(i); endif(f(i)>12.5&f(i)<80); a4(i)=12.5/f(i)*Gza(i).2; f4(i)=f(i);endend aw1=trapz(f1,a1);aw2=trapz(f2

52、,a2); aw3=trapz(f3,a3); aw4=trapz(f4,a4); aw=aw1+aw2+aw3+aw4disp('人體加權加速度均方根值=')disp(aw)Law=20*log10(aw*106)disp('加權振級=')disp(Law)路面不平度加速度均方根值= 0.35045611862153車輪加速度均方根值= 0.23907567188580車身加速度均方根值=eza =人體加速度均方根值=aw = 1.758474239477950e-004人體加權加速度均方根值= 1.758474239477950e-004Law = 44.90272021126767加權振級= 44.902720211267672）>> fs1=linspace(1.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;Gqn0=2.56*10(-8);n0=0.1;df=0.2;N=180;f=0.000001: