汽车理论编程题

1、第一章1）m1=2000 %装载质量 m2=1800 %整车整备质量 m=3880 ； %总质量 r0=0.367 ； %车轮半径 gt=0.85 ； % 传动系机械效率 f=0.013 ； % 波动阻力系数 CDA=2.77 ； % 空气阻力系数×迎风面积 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）用符号函数求最大速度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);%汽车旋转质量换算系数 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('加速度倒数时间曲线'

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(j)=1+(2*Iw1+4*Iw2)/(m*r02)+If*Ig5(j).2

8、*i02*gt/(m*(r0)2);%汽车旋转质量换算系数 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 ; %波动阻力系数 CDA=2.77 ; % 空气阻力系数×迎风面积 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 ; %波动阻力系数 CDA=2.77 ; % 空气阻力系数×迎风面积 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 %六工况循环行驶的百公里油耗%第一工况 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(j)=1+(2*Iw1+4*Iw2)/(m*r02)+If*Ig5(j).2*Ig0(k)2*gt/(m*(r0)2);%汽车旋转质量换算系数 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('燃油经济性-动力性曲线')xlabel('燃油经济性/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制动法规 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('利用附着系数与制动强度的关系曲线')xlabel('制动强度z/g')ylabel('利用附着系数g')gtext('空车前轴')gtext('空车后轴')gtext(&

36、#39;满载前轴')gtext('满载后轴')gtext('ECE法规')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('附着系数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;%求稳定性因素和特征车速K=m*(a/k2-b/k1)/L2;Uch=sqrt(1/K);disp('稳定性因素 K=');disp(K);disp('特征车速Uch=

41、9;);disp(Uch);运行结果为：稳定性因素 K= 0.0024特征车速Uch= 20.6053%绘制稳态横摆角速度增益曲线u=0:0.1:100;wr=(u/L)./(1+K*u.2);%转向不足时稳态横摆角速度增益K1=0;wr1=(u/L)./(1+K1*u.2);%中性转向时稳态横摆角速度增益plot(u,wr,u,wr1)axis(0 100 0 20);title('稳态横摆角速度增益曲线') 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%计算静态储备系数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('静态储备系数S.M.=')disp(SM)disp('前后轮侧偏角绝对值之差（

43、a1-a2）=')disp(f)disp('转弯半径之比R/R0')disp(b)运行结果为：静态储备系数S.M.= 0.1576前后轮侧偏角绝对值之差（a1-a2）= 0.0281转弯半径之比R/R0 1.1608%车速Ua=30.56m/s,瞬态响应的横摆角速度波动的固有(圆)频率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,瞬态响应的横摆角速度波动的固有(圆)频率w0')disp(w0)disp('阻尼比g')disp(g)disp('反应时间I')disp(I)disp('峰值反应时间 t')disp(t)运行结果为：车速Ua=30.56m/s,瞬态响应的横摆角速度波动的固有(圆)频率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: