《汽车理论matlab作业.docx》由会员分享，可在线阅读，更多相关《汽车理论matlab作业.docx（56页珍藏版）》请在三一文库上搜索。
1、确定一轻型货车的动力性能。1)绘制汽车驱动力与行驶阻力平衡图;2)求汽车最高车速与最大爬坡度；7/543)驶至70km/h所需绘制汽车行驶加速度倒数曲线；用计算机求汽车用n档起步加速行的加速时间。巳知数据略。（参见汽车理论习题第一章第3题） 解题程序如下:用Matlab语言 绘制汽车驱动力与行驶阻力平衡图ml=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 = 1.798; Iw2=3.598;Iw=2
2、*Iwl+4*Iw2;for i=l:69 n(i)=(i+l 1)*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)A4;end for j=l:5for i=l:69Ft(i,j)=Ttq(i)*ig(j)*iO*nT/r;ua(i, j)=0.377*r*n(i) / (ig(j) *i0);Fz(i,j)=CdA*ua(i,j)A2/21.15+mz*g*f;endendplot(ua,Ft,ua,Ff,ua,Ff+Fw)title(汽
3、车驱动力与行驶阻力平衡图)；xlabelfuaflan/h),)；ylabelfFt(N),)；gtextCFtl)gtext(*Ft2)gtextCFt3)gtext(,Ft4,)gtext(Ft5)gtext(Ff+Fw)(2)求最大速度和最大爬坡度for k=l:175 nl(k)=3300+k*0.1;Ttq(k)=-19.313+295.27*(nl(k)/1000)-165.44*(nl(k)/1000)A2 +40.874*(nl(k)/1000)A33.8445*(nl(k)/1000)A4;Ft(k)=Ttq(k)*ig(5)*iO*nT/r;ua(k)=O. 377*r*n
4、 l(k) / (ig(5)*i0);Fz(k)=CdA*ua(k)A2/21.15+mz*g*f;E(k)=abs(Ft(k)-Fz(k);endfor k=l:175if(E(k)=min(E)dispC汽车最高车速=);disp(ua(k);dispfkm/h1)；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)A3-3.8445*(n2(p)/ 1000)A4;Ft(p)=Ttq(p)*ig(l)*iO*nT/r;
5、ua(p)=0.377*r*n2 (p) /(ig(l)*i0);Fz(p)=CdA*ua(p)A2/ 21.15+mz*g*f;af(p)=asin(Ft(p)-Fz(p)/(niz*g);endfor p=l: 150if(af(p)=max(af)i=tan(af(p);dispC汽车最大爬坡度=)；disp(i);endend汽车最高车速=99.0679km/h汽车最大爬坡度=0.3518(3)计算2档起步加速到70km/h所需时Rfor i=l:69n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)A
6、2+40.874* (n(i)/ 1000)A3-3.8445*(n(i)/ 1000)A4;end for j=l:5for i=l:69deta=l+Iw/ (mz*rA2)+If*ig(j)A2*iOA2*nT/(niz*rA2);ua(i,j)=0.377*r*n(i)/(ig(j)*i0);a(i,j)=(Ttq(i)*ig(j)*iO*nT / r-CdA*ua(iJ) A2/ 21.15-niz*g*f)/ (deta*mz);end elsebl(i,j)=a(i-l,j);ul(i,j)=ua(i-l,j);endb(i,j)=l/bl(i,j);endendxl=ul(:,
7、l);yl=b(:,l);x2=ul(:,2);y2=b(:,2);x3=ul(:,3);y3=b(:,3);x4=ul(:,4);y4=b(:,4)；x5=ul(:,5);y5=b(:,5);plotfxl ,y 1 ,x2,y2,x3,y3,x4,y4,x5,y5);titled加速度倒数时间曲线)；axis(O 120 0 30);xlabel(ua(lm/h);ylabelfl/aj)；gtextfl/al)gtext(,l/a2)gtext(*l/a3)gtext(l/a4)gtextfl/aS*)32100102030405060708090u a (km/h)10987加速度倒数
8、时间曲线6 5 4 (ewes)gLfor i=l:69A=ua(i,3)-ua(69,2);if (A0) J=l;endB=ua(i,4)-ua(69,3);if(B0)k=i;endif(ua(i.4)=70)m=i;endt=ua(l,2)*b(l,2);for pl=2:6911 (p l)=(ua(p 1,2)-ua(p 1-1,2)*(b(pl ,2)+b(p 1-1,2)*0.5； t=t+tl(pl);endfor p2=j:69t2(p2)=(ua(p2,3)-ua(p2-l,3)*(b(p2,3)+b(p2-l,3)*0.5;t=t+t2(p2);endfor p3=k:
9、mt3(p3)=(ua(p3,4)-ua(p3-l,4)*(b(p3,4)+b(p3-l,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)*b(m,4)；tz=t/3.6;dispf加速时间=，);disp(tz);disp(s);加速时间=29.0585s二、计算与绘制题1中货车的1)汽车功率平衡图；2)最高档与次高档的等速百公里油耗曲线。巳知数据略。(参见汽车理论习题第二章第7题)解题程序如下:用Matlab语言ml=2000; m2=1800; mz=3880;
10、 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 = 1.798; Iw2=3.598;nl=815 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+11)*50;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)A2+40.874*(n(i)/ 1000)
11、A3-3.8445*(n(i)/ 1000)A4;Pe(i)=n(i)*Ttq(i)/9549;ua(i,j)=0.377*r*n(i) / (ig(j)*iO);Pz(iJ)=(niz*g*f*ua(iJ)/ 3600.+CdA*ua(i,j)A3/76140.)/nT;end 9/54endplot(ua,Pe,ua,Pz);titlef汽车功率平衡图力; xlabelfuatlmi/h),)； ylabel(Pe,Pz(kw),)；gtext(T)gtext(*ir)gtextCUP)gtext(IV)gtext(V)gtextCP 阻)10 / 54for j=l:5for i=l:8
12、Td(i)=-19.313+295.27*(nl(i)/1000.0)-165.44*(nl(i)/1000.0)A2+40 .874*(nl(i)/1000.0)A3-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*P d(l)+1326.8;b(2)=0.043072*Pd(2)A4-2.0553*Pd(2)A3+36.657*Pd(2)A2-303.98
13、*Pd +1354.7;b(3)=0.0068164*Pd(3)A4-0.51184*Pd(3)A3+14.524*Pd(3)A2-189.7 5*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)+l 122.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.7 14
14、*Pd(6)+1051.2;b(7)=0.00028230*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(:,l),;u2=u(:,2);u3=u(:,3);u4=u(:,4);u5=u(:,5);B l=polyfit(u l,b,3);B2=polyfit(u2,b,3);B3=polyfit(u3,b,3);B4=polyfit(u4,b,3);B
15、5=polyfit(u5,b,3);for q=l:69bh(q,l)=polyval(B 1 ,ua(q,l);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(B5,ua(q,5);end for i=l:5for q=l:69Q(q,i)=Pz(q,i)*bh(q,i)/(1.02*ua(q,i)*7.05);endendplot(ua(:,4),Q(:,4),ua(:,5),Q(:,5);title(四档五档等速百公里油耗图上xlab e
16、lfuaflmi/h),)；ylabel(Qs(L/ 100km),)；40100四档和五档等速百公里油耗曲线1520406080100120ua/(km/h)5 0 5 03 3 2 2Eaoo 3)0三、改变L3题中轻型货车的主减速器传动比，做出为5.17. 5.43、 583、617、6.33时的燃油经济性一加速时间曲线，讨论不同。值对汽车性能的影响。Matlab 程序：ml=2000;m2=1800;m=3880;r0=0.367;gt=0.85;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iwl = 1.798;Iw2=3.598;Ig5=5.56 2.769
17、 1.644 1.00 0.793;Ig0=5.17 5.43 5.83 6.17 6.33;B=1326.8 -416.46 72.3791354.7 -303.98 36.6571284.4 -189.75 14.5241122.9 -121.59 7.00351141.0 -98.893 4.47631051.2 -73.714 2.8593-5.8629 0.17768;-2.0553 0.043072;-0.51184 0.0068164;-0.18517 0.0018555;-0.091077 0.00068906;-0.05138 0.00035032;1233.9 -84.47
18、8 2.9788 -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=l: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).A2+40.874.*(n(i)/1000). A3-3.8445.*(n(i)/1000). A4;F5(i,k)=0.013*388
19、0*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,l)+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;endendual=25; sl=50;Fa5=0.013*3880*9.8+2.77.*ual.A2/21.15;Pe5=Fa5.*ual/(3600*0.85);dl=polyfit(Pe(l,:),b5(l,:),3);ba l=polyval(d 1 ,Pe5);d2=polyfit(Pe(2,:
20、),b5(2,:),3);ba2=polyval(d2 ,Pe5);15 / 54d3=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 ba2 ba3 ba4 ba5;Qa l=Pe5.*ba*50/(ua 1 * 102*7)ua2=25:40;Q2=l+(2*Iwl+4*Iw2)/(m*r0A2)+If*Ig5(2).A2*Ig0(
21、2)A2*gt/(m*r0A2)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(l,:),b5(l,:),3);bb l=polyval(db l,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,:),3);bb4=polyval(db4 ,Pb5);
22、db5=polyfit(Pe(5,:),b5(5,:),3);bb5=polyval(db5,Pb5);16 / 54bb=bblbb2bb3bb4bb5;Pb=Pb5Pb5Pb5Pb5Pb5;Qb=Pb.*bb/(367.1*7);for i=l:5for j=l: 15qb(i,j)=Qb(i,j)+Qb(i,j+1);endendQb2=sum(qb)ua3=40; s2=250;Fc5=0.013*3880*9.8+2.77.*ua3.A2/21.15;Pc5=Fc5.*ua3/(3600*0.85);dcl=polyfit(Pe(l,:),b5(l,:),3);be 1 =polyv
23、al(dc 1, Pc5);17 / 54dd5=polyfit(Pe(5,:),b5(5,:),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=bcl bc2 bc3 bc4 bc5;Qc3=Pc5.*bc*250/(ua3* 102*7)ua4=40:50;Q4=l+(2*Iwl+
24、4*Iw2)/(m*r0A2)+If*Ig5(4).A2*Ig0(4)A2*gt/(m*r0A2)Fd5=0.013*3880*9.8+2.77.*ua4.A2/21.15+Q4*m*0.2;Pd5=Fd5.*ua4/(3600*0.85);dd 1 =polyfit(Pe( 1,:) ,b5( 1,:) ,3);bd l=polyval(dd 1 ,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
25、(4,:),b5(4,:),3);bd4=polyval(dd4 ,Pd5);bd5=polyval(dd5,Pd5);bd=bdlbd2bd3bd4bd5;Pd=Pd5Pd5Pd5Pd5Pd5;Qd=Pd.*bd/(367.1*7);for i=l:5for j=l: 10qd(i,j)=Qd(i,j)+Qd(i,j+1);endendQd4=sum(qd)ua5=50; s2=250;Ff5=0.013*3880*9.8+2.77.*ua5.A2/21.15;Pf5=Ff5.*ua5/(3600*0.85);19 / 54dfl=polyfit(Pe(l,:),b5(l,:),3);bfl
26、=polyval(dfl ,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,:),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/
27、0.36*0.299;Qg6=Qg Qg Qg Qg QgQ=QalQb2Qc3Qd4Qf5Qg6;Qz=sum(Q)/1075*100for k=l:5 for i=l:3401;for j=l:5;n(i)=i+599;ua(i,j)=0.377*r0*n(i)./(Ig5(j)*IgO(k);Q(j)=l + (2*Iwl+4*Iw2)/(m*r0A2)+If*Ig5(j).A2*Ig0(k)A2*gt/(m*(r0)A2)；endenduamax=max(ua);uamin=min(ua);ua2=uamin(2):uamax(2);n=ua2*Ig0(k)*Ig5(2)/(0.377
28、*0.367);Ttq=-19.313+295.27*(n/1000)-165.44*(n/1000).A2+40.874*(n/l000).A3-3.8445*(n/1000). A4;Ft2=Ttq*Ig0(k)*Ig5(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 /2)/ 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
29、.44*(n3/1000).A2+40.874*( n3/1000). A3-3.8445*(n3/1000). A4;Ft3=Ttq3*Ig0(k)*Ig5(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*Ig0(k)*Ig5(4)/(0.377*0.367);Ttq4=-19.313+295.27*(n4/1000)-165.44*(n4/1000).A2+40.874*( n4/1000).
30、 A3-3.8445*(n4/1000). A4;Ft4=Ttq4*IgO(k)*Ig5(4)*O.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 -/h4)/ 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.17)gtext(543)gtext(583)gtext(6.17)gtext(,6.33,) title(燃油经济性-动
31、力性曲线)xlabel(燃油经济性/L( 1 OOkm)-1)ylabelft/s)四、一中型货车装有前后制动器分开的双管路制动系，其有关参数如下:载荷质量(kg)质心商 hg/m轴距L/m质心至前轴距离制动力分配系数0a/m空栽40800.8453.9502.1000.38满载92901.1703.9502.9500.381)计算并绘制利用附着系数曲线和制动效率曲线2)求行驶车速Ua=30km/h,在夕= 0.80路面上车轮不抱死的制动距离。计算时取制动系反应时间G=0.02s,制动减速度上升时间r；= 0.02so3)求制动系前部管路损坏时汽车的制动距离s,制动系后部管路损坏时汽车的制动距
32、离Matlab 程序：ml=4080;hgl=0.845;al=2.100;m2=9290;hg2= 1.17;a2=2.95;beta=0.38;L=3.95;z=0:0.05:1gf l=beta.*z*L./(L-al+z*hgl);gf2=beta.*z*L./(L-a2+z*hg2);gr 1=( 1 -beta) .*z*L. / (al -z*hgl);gr2=(l-beta).*z*L./(a2-z*hg2);g=z；for i=l: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.7
33、4;endendzl=O. 15:0.01:0.3;g4=zl-0.08;plot(z,gfl,-.,z,gf2,z,grl,-.,z,gr2,z,g,z,g3,xk,zl,g4,x,)axis(0 1 0 1.2)title(利用附着系数与制动强度的关系曲线)xlabel(，制动强度z/)ylabelf利用附着点数g)gtext(空车前轴)gtext(空车后轴)gtext(满载前轴)gtext(满载后轴)gtext(ECE 法规)利用附着系数与制动强度的关系曲线1 0.8 0.6 0.4 0.200.10.20.30.40.50.60.70.80.91制动强度z/g6熬瞩*超田=长C=0:0
34、.05:1;Erl=(al/L)./(l-beta)+C*hgl/L)*100;Ef= (L-a2)/L./(beta-C*hg2/L)*1OO;Er= (a2/L)./(l-beta)+C*hg2/L)*1OO;plot(C,Er,C,Ef,C,Erl)axis(O 1 0 100)title(前后附着效率曲线，)xlabelf附着系数C)ylabel(制动效率(%)1)gtext(满载)gtext(Ef)gtext(Er)gtext(空载)gtext(Er)56 / 5400.10.20.30.40.50.60.70.80.91附着系数C前后附着效率曲线 100908070605040配O
35、Cl=0.8El=(akl/L)./(l-beta)+Cl*hgl/L);E2=(am2 / L) / (1 -b eta)+C1 *hg2 / L);al=El*Cl*9.8;a2=E2*Cl*9.8;ua=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*am2);disp(，满载时不抱死的制动距离刁disp(s2)dispf空载时不抱死的制动距离=)disp(sl)满载时不抱死的制动距离=5.3319空轴时不抱死的制动距离=6.8406beta3=lb
36、eta4=0Elr=(al/L)/(l-beta4)+Cl*hgl/L);Ekf= (L-a 1) / (beta3*L-C l*hgl);Emf=(L-a2)/L./(beta3-C 1 *hg2 / L);Emr=(a2/L)./(l-beta4)+C 1 *hg2 / L);akr=0.8*9.8*Ekr;akf=0.8*9.8*Ekf;amr=0.8*9.8*Emr;anif=0.8*9.8*Emf;skr=(i21+i22/2)*ua/3.6+uaA2/(25.92*akr);sld= (i2 l+i22/2)*ua/3.6+uaA2/(2 5.92*akf);smf=(i21+i2
37、2/2)*ua/3.6+uaA2/(25.92*amf);smr=(i21+i22/2)*ua/3.6+uaA2/(25.92*amr);disp(空车后管路失效时制动距离)disp(sld)disp(空车前管路失效时制动距离)disp(skr)disp(满载后管路失效时制动距离)disp(smf)disp(满载前管路失效时制动距离)disp(smr)运行结果为：空车后管路失效时制动距离8.0879空车前管路失效时制动距离10.0061满载后管路失效时制动距离13.5986满轴前管路失效时制动距离7.5854五、二自由度轿车模型的有关参数如下:总质绕Oz轴转动惯轴距质心至前轴距离质心至后轴距离
38、前轮总侧偏刚度后轮总侧偏刚度转向系总传动比试求：m= 1818.2kg/. =3885Ag /L=3.048ma= 1.463mb= 1.585mki=-62618N/rad k2=-110185N/radi=201)稳定性因数K、特征车速Ueh。2）稳态横摆角速度增益曲线外飞、车速S35m/s时的转向灵敏度冬。3）静态储备系数S.M.,侧向加速度为0.4g时的前、后轮侧偏角绝对值之差%-%与转弯半径的比值R/Ro（Ro=15m）e4）车速u=30.56m/s时,瞬态响应的横摆角速度波动的固有（圆）频率必、阻尼比八反应时间，与峰值反应时间Matlab 程序：m=1818.2;Iz=3885;L
39、=3.048;a=1.463;b=1.585;kl=-62618;k2=-110185;i=20;g=9.8;R0=15;ul=30.56;K=m/LA2*(a/k2-b/kl);uch=l/KA(l/2);u=0:0.05:25;s=u/L./(l+K*u.A2);disp(稳定因素K=);disp(K);dispf特征车速uch=);disp(uch);plot(u,s);xlabelfua/ (m/s),);ylabel(稳态横摆增益，);title(汽车的稳态横态摆角速度增益曲线)；disp(a=22.35m/s时，转向灵敏度为)；disp(s(448);SM=k2/(kl+k2)-a
40、/L;A=K*0.4*g*L;Q=L/RO;R=L/(Q-A);r=R/RO;disp(静态储备系数S.M.=);disp(SM);disp(，前后轮侧偏角绝对值之差(al-a2)刁；disp (A);disp(单弯半径的比值引；disp(r);wO=L/ul*(kl*k2*(l+K*ulA2)/(m*Iz)A(l/2);n=(-m*(aA2*kl+bA2*k2)-Iz*(kl+k2)/ (2*L*(m*Iz*kl*k2*(l+K*ulA2)A(l/2)；t=atan(l-nA2)A(l/2)/(-m*ul*a*wO/(L*k2)-n)/(wO*(l-nA2)A(l/2);e=atan(l-n
41、A2)A(l/2)/n)/(wO*(l-nA2)A(l/2)+t;disp(，横摆角速度波动时的固有频率为)；disp(wO);disp（阻尼比为）;disp(n);disp（反应时间为）;disp(t);dispf达到第一峰值的时间为）;disp(e);3.5汽车的稳态横态摆角速度增益曲线32.521.510.5520251015ua/(m/s)运算结果：稳定因素K= 0.0024特征车速uch= 20.6053 ua=22.35m/s时，转向灵敏度为 3.3690静态储备系数S.M尸 0.1576前后轮侧偏角绝对值之差（al-a2）=0.0281转弯半径的比值= 1.1608横摆角速度波动
42、时的固有频率为5.5758阻尼比为0.5892反应时间为0.1811达到第一峰值的时间为 0.3899六、车身-车轮双质系统参数：/）=1.5比4=0.25,7 = 9, = 10。“人体-座椅”系统参数： f,. = 3小，=0.25 o车速=20m/s ,路面不平度 系数G/%）=2.56xl0-s,参考空间频率no=O.lm-】。计算时频率步长9=0.2Hz, 计算频率点数N = 180。1）计算并画出帽频特性|./“、|马/。|、和均方根值谱向、Gnf）、，GQ ）谱图。进一步计算?、%、4、七、4卬值2）改变“人体-座椅”系统参数：，=1.56=0.1250.5。分析、心值随人、乙的
43、变化。3）分别改变车身-车轮双质，系统参数：人=0.25 3%，S = 0.125 0.5 , y = 4.5 18, = 5 。用、三个响应均方根值随以上四个系统参数变化的曲线。Matlab 程序：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=(l-ff0.A2).*(l+r-l/u*ff0.A2-l).A2+4*gA2*ff0.A2.*r-(l/u+l)*ff0.A2.A2;zlq=r*sqrt(l-ffO.A2).A2+4*gA2*ffO.A2)./ d);z2zl=sqrt(l+(2*g*ff0).八 2)./(l-ff0.八 2).八 2+(2*g*ff0).八 2);ffs=f/fs;pz2=sqrt(l+(2*gs*ffs).A2)./(l-ffs.A2).A2+(2*gs*ffs).A2)Ioglog(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 10)xlabel(，激振频率 f/HZ)|10w=2*pi*f;Gzl=w.A2.*