《汽车理论》习题5-6章答案余志生4版-6页文档资料
汽车理论习题答案
汽车理论课后习题答案(完整版)从第一章到第六章（附图文）第一章1.1、试说明轮胎滚动阻力的定义、产生机理和作用形式?答：1）定义：汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动阻力。
2）产生机理：由于轮胎内部摩擦产生弹性轮胎在硬支撑路面上行驶时加载变形曲线和卸载变形曲线不重合会有能量损失，即弹性物质的迟滞损失。
这种迟滞损失表现为一种阻力偶。
当车轮不滚动时，地面对车轮的法向反作用力的分布是前后对称的；当车轮滚动时，由于弹性迟滞现象，处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力，这样，地面法向反作用力的分布前后不对称，而使他们的合力Fa 相对于法线前移一个距离a, 它随弹性迟滞损失的增大而变大。
即滚动时有滚动阻力偶矩 a F T z f = 阻碍车轮滚动。
3）作用形式：滚动阻力fw F f = rT F ff = （f 为滚动阻力系数）1.2、滚动阻力系数与哪些因素有关？提示：滚动阻力系数与路面种类、行驶车速以及轮胎的构造、材料、气压等有关。
1.3、解答：1）（取四档为例） 由u F n u n Tq Tq F t t →⇒⎪⎭⎪⎬⎫→→→ 即 ri i T F To g q t η=432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-=og i i rnu 377.0=行驶阻力为w f F F +： 215.21a D w f U A C Gf F F +=+ 2131.0312.494aU +=由计算机作图有※本题也可采用描点法做图：由发动机转速在m in /600n min r =，m in /4000n max r =，取六个点分别代入公式：……………………………… 2）⑴最高车速：有w f t F F F +=⇒2131.0312.494a t U F += 分别代入a U 和t F 公式：2)09.6*83.53697.0*377.0(131.0312.494367.085.0*83.5*9.6*n T q += 把q T 的拟和公式也代入可得： n>4000而4000m ax =n r/min∴93.9483.5*0.14000*367.0*377.0max ==U Km/h ⑵最大爬坡度：挂Ⅰ档时速度慢，Fw 可忽略： ⇒)(max w f t i F F F F +-=⇒GfF Git -=max⇒013.08.9*388014400max max -=-=f G F i t =0.366（3）克服该坡度时相应的附着率 zxF F =ϕ 忽略空气阻力和滚动阻力得：6.0947.12.3*366.0/=====a il l a i F Fi zϕ 3）①绘制汽车行驶加速倒数曲线（已装货）：40.0626)(1f D g du dt a -==δ （G FwFt D -=为动力因素） Ⅱ时，22022111r i i I m r I m Tg f wηδ++=∑2222367.085.0*83.5*09.3*218.038001367.0598.3798.1380011+++= =1.128ri i T F To g q t η=432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-=215.21a D w U A C F =由以上关系可由计算机作出图为：②用计算机求汽车用Ⅳ档起步加速至70km/h 的加速时间。
汽车理论课后题答案清华大学余志生
第一章1.1、试说明轮胎滚动阻力的定义、产生机理和作用形式?答：1）定义：汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动阻力。
2）产生机理：由于轮胎内部摩擦产生弹性轮胎在硬支撑路面上行驶时加载变形曲线和卸载变形曲线不重合会有能量损失，即弹性物质的迟滞损失。
这种迟滞损失表现为一种阻力偶。
当车轮不滚动时，地面对车轮的法向反作用力的分布是前后对称的；当车轮滚动时，由于弹性迟滞现象，处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力，这样，地面法向反作用力的分布前后不对称，而使他们的合力Fa 相对于法线前移一个距离a, 它随弹性迟滞损失的增大而变大。
即滚动时有滚动阻力偶矩 a F T z f = 阻碍车轮滚动。
3）作用形式：滚动阻力 fw F f = r T F ff =（f 为滚动阻力系数）1.2、滚动阻力系数与哪些因素有关？提示：滚动阻力系数与路面种类、行驶车速以及轮胎的构造、材料、气压等有关。
1.3、确定一轻型货车的动力性能（货车可装用4档或5档变速器，任选其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车的最高车速、最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速倒数曲线，用图解积分法求汽车有Ⅱ档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用Ⅱ档起步加速至70km/h 的加速时间。
轻型货车的有关数据：汽油发动机使用外特性的Tq —n 曲线的拟合公式为432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-= 式中, Tq为发功机转矩(N ·m)；n 为发动机转速(r ／min)。
发动机的最低转速n min =600r/min ，最高转速n max =4000 r ／min装载质量 2000kg 整车整备质量 1800kg 总质量 3880 kg 车轮半径 0.367 m传动系机械效率 ηт＝0.85 波动阻力系数 f ＝0.013 空气阻力系数×迎风面积 C D A ＝2.772m 主减速器传动比 i0＝5.83飞轮转功惯量 I f ＝0.218kg ·2m 二前轮转动惯量 I w1＝1.798kg ·2m四后轮转功惯量 I w2＝3.598kg ·2m变速器传动比 i g (数据如下表)轴距 L ＝3.2m质心至前铀距离(满载) α＝1.947m 质心高(满载) h g ＝0.9m解答：1）（取四档为例）由uF n u n Tq Tq F t t →⇒⎪⎭⎪⎬⎫→→→ 即ri i T F To g q t η=432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-= og i i rn u 377.0=行驶阻力为w fF F +：215.21a D w f U A C Gf F F +=+ 2131.0312.494aU +=由计算机作图有※本题也可采用描点法做图：由发动机转速在m in /600n min r =，m in /4000n max r =，取六个点分别代入公式：……………………………… 2）⑴最高车速：有w f tF F F +=⇒2131.0312.494a t U F += 分别代入a U 和t F 公式：2)09.6*83.53697.0*377.0(131.0312.494367.085.0*83.5*9.6*n T q +=把q T 的拟和公式也代入可得： n>4000而4000m ax =n r/min∴93.9483.5*0.14000*367.0*377.0max ==U Km/h⑵最大爬坡度：挂Ⅰ档时速度慢，Fw 可忽略：⇒)(max w f t i F F F F +-=⇒GfF Gi t -=max⇒013.08.9*388014400max max-=-=f G F i t =0.366（3）克服该坡度时相应的附着率 zxF F =ϕ 忽略空气阻力和滚动阻力得：6.0947.12.3*366.0/=====a il l a i F Fi z ϕ 3）①绘制汽车行驶加速倒数曲线（已装货）：40.0626)(1f D g du dt a -==δ（GFwFt D -=为动力因素）Ⅱ时，22022111r i i I m r ImTg f wηδ++=∑2222367.085.0*83.5*09.3*218.038001367.0598.3798.1380011+++= =1.128ri i T F To g q t η=432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-=215.21a D wU A C F = 由以上关系可由计算机作出图为：②用计算机求汽车用Ⅳ档起步加速至70km/h 的加速时间。
汽车理论课后习题答案(余志生版)()之欧阳计创编
第一章时间：2021.02.11 创作：欧阳计1.1、试说明轮胎滚动阻力的定义、产生机理和作用形式? 答：1）定义：汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动阻力。
2）产生机理：由于轮胎内部摩擦产生弹性轮胎在硬支撑路面上行驶时加载变形曲线和卸载变形曲线不重合会有能量损失，即弹性物质的迟滞损失。
这种迟滞损失表现为一种阻力偶。
当车轮不滚动时，地面对车轮的法向反作用力的分布是前后对称的；当车轮滚动时，由于弹性迟滞现象，处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力，这样，地面法向反作用力的分布前后不对称，而使他们的合力Fa相对于法线前移一个距离a, 它随弹性迟滞损失的增大而变大。
即滚动时有滚动阻力偶矩阻碍车轮滚动。
3）作用形式：滚动阻力（f为滚动阻力系数）1.2、滚动阻力系数与哪些因素有关？提示：滚动阻力系数与路面种类、行驶车速以及轮胎的构造、材料、气压等有关。
1.3、解答：1）（取四档为例）由即行驶阻力为：由计算机作图有※本题也可采用描点法做图：由发动机转速在，，取六个点分别代入公式：………………………………2）⑴最高车速：有分别代入和公式：把的拟和公式也代入可得：n>4000而r/min∴ Km/h⑵最大爬坡度：挂Ⅰ档时速度慢，可忽略：=0.366（3）克服该坡度时相应的附着率忽略空气阻力和滚动阻力得：3）①绘制汽车行驶加速倒数曲线（已装货）：40.0626（为动力因素）Ⅱ时，1.128由以上关系可由计算机作出图为：②用计算机求汽车用Ⅳ档起步加速至70km/h的加速时间。
（注：载货时汽车用Ⅱ档起步加速不能至70km/h）由运动学可知：即加速时间可用计算机进行积分计算求出，且曲线下两速度间的面积就是通过此速度去件的加速时间。
经计算的时间为： 146.0535s1.4、空车、满载时汽车动力性有无变化？为什么？答：汽车的动力性指汽车在良好路面上直线行驶时，由纵向外力决定的所能达到的平均行驶速度。
《汽车理论》清华大学余志生版--期末考试复习资料
汽车理论习题集必考试题一、单项选择题（在每小题列出的四个备选项中,只有一项是最符合题目要求的，请将其代码写在该小题后的括号内）1、评价汽车动力性的指标是（A ）A．汽车的最高车速、加速时间和汽车能爬上的最大坡度B．汽车的最高车速、加速时间和传动系最大传动比C．汽车的最高车速、加速时间和传动系最小传动比D．汽车的最高车速、加速时间和最大驱动力2、汽车行驶速度( B ）A．与发动机转速、车轮半径和传动系传动比成正比B．与发动机转速和车轮半径成正比，与传动系传动比成反比C．与发动机转速和传动系传动比成正比，与车轮半径成反比D．与发动机转速成正比，与车轮半径和传动系传动比成反比3、汽车在水平路面上加速行驶时，其行驶阻力包括（B ）。
A。
滚动阻力、空气阻力、坡度阻力B．滚动阻力、空气阻力、加速阻力C．空气阻力、坡度阻力、加速阻力D。
滚动阻力、空气阻力、坡度阻力、加速阻力4、汽车等速上坡行驶时，其行驶阻力包括（ A ）。
A. 滚动阻力、空气阻力、坡度阻力B．滚动阻力、空气阻力、加速阻力C．空气阻力、坡度阻力、加速阻力 D. 滚动阻力、空气阻力、坡度阻力、加速阻力5、汽车加速上坡行驶时，其行驶阻力包括（ D ）.A。
滚动阻力、空气阻力、坡度阻力B．滚动阻力、空气阻力、加速阻力C．空气阻力、坡度阻力、加速阻力 D. 滚动阻力、空气阻力、坡度阻力、加速阻力6、汽车行驶时的空气阻力包括（D ）。
A．摩擦阻力和形状阻力 B. 摩擦阻力和干扰阻力C．形状阻力和干扰阻力D。
摩擦阻力和压力阻力7、汽车行驶时的空气阻力（B ）。
A. 与车速成正比B。
与车速的平方成正比C. 与车速的3次方成正比D. 与车速的4次方成正比8、汽车行驶时的空气阻力(C )。
A. 与迎风面积和车速成正比B。
与迎风面积的平方和车速成正比C. 与迎风面积和车速的平方成正比D。
与迎风面积的平方和车速的平方成正比9、同一辆汽车，其行驶车速提高1倍,空气阻力增大(C ）。
汽车理论习题答案及MATLAB编程
汽车理论课后习题答案（部分题目）教材：汽车理论（第5版）清华大学余志生主编机械工业出版社出版本文档包含习题答案：第一章：1.1 1.2 1.7 1.8第二章：2.2 2.3 2.4第三章第四章：4.1 4.2 4.3 4.5第五章：5.1 5.8 5.10 5.11 5.12 5.14 5.17第六章：6.1~6.5其中第三章 4.3 5.11 6.2 6.4 6.5所要求的MATLAB编程见本作者另一百度文库文档。
获取方法：点击左上角本作者的昵称，进入本作者的文库，《汽车理论课后习题答案—附MATLAB编程》第三章：解：（完整解题含MATLAB编程见另一文件）根据“货车最高车速为100km/h左右”，为满足五种主减速器与两种变速器匹配后车速均能达到100km/h，参照表中发动机转速，可选3006r/min、3403r/min、3804r/min三种。
结合万有特性图，考虑到动力性与燃油经济性的统一，选择转速3403r/min进行计算。
由于本人能力有限，暂难以完成EPA循环工况下的燃油消耗量计算，故以“货车采用最高档、90％负荷率工况下”的等速百公里油耗代替表示燃油经济性。
取汽油密度为0.7kg/L。
燃油经济性以等速百公里油耗Qs(L/100km)表示；动力性以2挡起步加速行驶至96.6km/h车速所用时间t(s)表示。
下面运用MATLAB编程实现计算和绘图，编程（m文件）列于最后，运行结果（4挡）如下：t=48.4235t=46.3511t=41.5557t=33.3796t=30.556490％负荷Pe=55.5647Qs4=25.103126.365628.307829.958730.7356图1加速度倒数-速度曲线图2燃油经济性-加速时间曲线由此可得出：主减速器传动比越大，动力性越好，燃油经济性越差。
与4挡变速器相比，采用5挡变速器时动力性和燃油经济性有所提高和改善。
（完整解题含MATLAB 编程见另一文件）图1图25.8某种汽车的质心位置、轴距和前后轮胎的型号已定。
汽车理论课后习题答案及MATLAB编程
汽车理论课后习题答案—附MATLAB编程教材：汽车理论（第5版）清华大学余志生主编机械工业出版社出版本文档包含第三章 4.3 5.11 6.2 6.4 6.5所要求的MATLAB编程其它习题答案：第一章：1.1 1.2 1.7 1.8第二章：2.2 2.3 2.4第三章第四章：4.1 4.2 4.3 4.5第五章：5.1 5.8 5.10 5.11 5.12 5.14 5.17第六章：6.1~6.5见本作者另一百度文库文档。
获取方法：点击左上角本作者的昵称，进入本作者的文库，《汽车理论课后习题答案（含MATLAB编程）》第三章：m=3880;%总质量g=9.8;%重力加速度r=0.367;%车轮半径eta_t=0.85;%传动系机械效率f=0.013;%滚动阻力系数CDA=2.77;%空气阻力系数*迎风面积i0=[5.175.435.836.176.33];%主减速器传动比If=0.218;%飞轮转动惯量Iw1=1.798;Iw2=3.598;%两前轮/四后轮转动惯量Iw=Iw1+Iw2;ig4=[6.09 3.091.711];%4挡变速器传动比ig5=[5.56 2.7691.64410.793];%5挡变速器传动比n=600:1:4000;%发动机转速Tq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;%4挡发动机for i=1:5Ft1=Tq*ig4(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig4(2)*i0(i)/r;Ft3=Tq*ig4(3)*i0(i)/r;Ft4=Tq*ig4(4)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig4(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig4(2));u3=0.377*r*n/(i0(i)*ig4(3));u4=0.377*r*n/(i0(i)*ig4(4));F1=m*g*f+CDA*u1.^2/21.15;%各档对应转速下的行驶阻力F2=m*g*f+CDA*u2.^2/21.15;F3=m*g*f+CDA*u3.^2/21.15;F4=m*g*f+CDA*u4.^2/21.15;delta1=1+(Iw+If.*ig4(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig4(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig4(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig4(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;figure(1)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=96.6;x1=[];x2=[];x3=[];x4=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];elseif u4(k)<=u4maxx3=[k];endendq2=max(x2);q3=max(x3);ua4=u4(q2:q3);a4=b4(q2:q3);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s=[s2s3s4];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs4(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs4=');disp(Qs4);tt=[48.423546.351141.555733.379630.5564];figure(2)plot(Qs4,tt,'*')hold onxi=25:0.001:31;pp=interp1(Qs4,tt,xi,'cubic');plot(xi,pp)title('燃油经济性-加速时间曲线');xlabel('燃油经济性--等速百公里油耗(L/100km)');ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')%5挡发动机for i=1:5Ft1=Tq*ig5(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig5(2)*i0(i)/r;Ft3=Tq*ig5(3)*i0(i)/r;Ft4=Tq*ig5(4)*i0(i)/r;Ft5=Tq*ig5(5)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig5(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig5(2));u3=0.377*r*n/(i0(i)*ig5(3));u4=0.377*r*n/(i0(i)*ig5(4));u5=0.377*r*n/(i0(i)*ig5(5));F1=m*g*f+CDA*u1.^2/21.15;%各档对应转速下的行驶阻力F2=m*g*f+CDA*u2.^2/21.15;F3=m*g*f+CDA*u3.^2/21.15;F4=m*g*f+CDA*u4.^2/21.15;F5=m*g*f+CDA*u5.^2/21.15;delta1=1+(Iw+If.*ig5(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig5(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig5(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig5(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta5=1+(Iw+If.*ig5(5).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);a5=(Ft5-F5)/(delta5*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;b5=1./a5;figure(3)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4,u5,b5);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext('1/a5')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=max(u4);u5min=u4max;u5max=96.6;x1=[];x2=[];x3=[];x4=[];x5=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];endendq2=max(x2);ua4=u4(q2:y);a4=b4(q2:y);for l=1:3401if u5(l)<=u5minx3=[l];elseif u5(l)<=u5maxx4=[l];endendq2=max(x2);q3=max(x3);q4=max(x4);ua5=u5(q3:q4);a5=b5(q3:q4);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s5=trapz(ua5,a5);%五挡运行时间s=[s2s3s4s5];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;Pe5=Ft5.*u5./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs5(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs5=');disp(Qs5);tt=[93.973058.283444.678445.366744.9793];figure(2)plot(Qs5,tt,'*')hold onxi=30:0.001:100;pp=interp1(Qs5,tt,xi,'cubic');plot(xi,pp)%title('燃油经济性-加速时间曲线');%xlabel('燃油经济性--等速百公里油耗(L/100km)');%ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')运行结果：t=48.4235t=46.3511t=41.5557t=33.3796t=30.556490％负荷Pe=55.5647Qs4=25.103126.365628.307829.958730.7356图1加速度倒数-速度曲线图2燃油经济性-加速时间曲线4.3%空载(no load)—1；满载(full load)—2m1=4080;m2=9290;hg1=0.845;hg2=1.170;L=3.950;a1=2.100;a2=2.950;b1=L-a1;b2=L-a2;beta=0.38;%利用附着系数与制动强度的关系曲线z=0:0.01:1;phi_f1=L*beta.*z./(b1+z*hg1);%前轮利用附着系数(空载)phi_r1=L*(1-beta).*z./(a1-z*hg1);%空载时后轮利用附着系数(空载)phi_f2=L*beta.*z./(b2+z*hg2);%前轮利用附着系数(满载)phi_r2=L*(1-beta).*z./(a2-z*hg2);%后轮利用附着系数(满载)phi=z;%最理想情况figure(1);plot(z,phi_f1,'g--',z,phi_f2,'g-',z,phi_r1,'b--',z,phi_r2,'b-',z,phi,'k-'); axis([0,1,0,1]);grid on;box off;axis square;title('利用附着系数与制动强度的关系曲线');xlabel('制动强度z/g');ylabel('利用附着系数{\phi}');%ECE法规要求界限hold on;z1=0.2:0.01:0.8;z2=0.15:0.01:0.3;z3=0.3:0.01:1;phi1=(z1+0.07)/0.85;phi2=z2-0.08;phi3=z2+0.08;phi4=(z3-0.02)/0.74;plot(z1,phi1,'r-.',z2,phi2,'r-.',z2,phi3,'r-.',z3,phi4,'r-.'); legend('{\phi}_{f}(空载)','{\phi}_{f}(满载)','{\phi}_{r}(空载)',...'{\phi}_{r}(满载)','{\phi}=z','ECE法规要求界限');legend('Location','northwest');%制动效率曲线Ef1=z./phi_f1*100;Er1=z./phi_r1*100;Ef2=z./phi_f2*100;Er2=z./phi_r2*100;figure(2);plot(phi_f1,Ef1,'k',phi_r1,Er1,'k',phi_f2,Ef2,'b',phi_r2,Er2,'b'); axis([0,1,0,100]);grid on;box off;axis square;title('制动效率曲线');xlabel('附着系数{\phi}');ylabel('制动效率(%)');text(0.25,90,'E_f');text(0.6,88,'满载');text(0.8,90,'E_r');text(0.7,65,'空载');text(0.72,72,'E_r');程序运行结果如下：图1图2 5.11L=3.048;K=0.0024;i=20;ua=0:0.01:120;f=1/L.*(1./(3.6./ua+K.*ua./3.6));syms uF(u)=(1/L.*(1./(1./u+K.*u)))./i;plot(ua,f,'b-');xlabel('u_a(km/h)');ylabel('\omega_r/\delta)_s');F(22.35)运行结果：ans=0.1667386094155036.2f=0.1:0.1:100;Gq1=5.12*10^(-5)./f.^2;Gq2=2.02*10^(-3).*f./f;Gq3=7.98*10^(-2).*f.^2;figure(1);loglog(f,Gq1,'b-');box off;axis square;title('位移功率谱密度'); xlabel('f/Hz');ylabel('G_q(f)/(m^2s)');figure(2);loglog(f,Gq2,'b-');box off;axis square;title('速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''(f)/(m^2/s)');figure(3);loglog(f,Gq3,'b-');box off;axis square;title('加速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''_''(f)/(m^2/s^3)');运行结果如下：6.4f1=0.1:0.01:1;f2=1:0.01:10;G1=0.010106.*f1.^2;G2=0.010106.*f2./f2;loglog(f1,G1,'b-');hold onloglog(f2,G2,'b-');title('车身加速度的功率谱密度'); xlabel('f/Hz');ylabel('G_z_''_''(f)/m^-^1'); axis([0.1,10,0.0001,0.1]);运行结果如下：6.51.计算说明1）①幅频特性()122222114z q λζλγ⎡⎤-+⎢=⎢⎥∆⎣⎦，其中0ωλω=，()2222222111141λλγζλγλμμ⎡⎤⎡⎤⎛⎫⎛⎫∆=-+--+-⎢⎥⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎣⎦⎣⎦；()12222222211414z z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦；()1222222221414s s s s s p z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦，其中s s ωλω=；②均方根谱()()221~j 2z q H f f ωππ=；()()222~j 2z qH f f ωππ== ；()()2221~1j 2p qz z H f f z qωππ== ；③其他值()0.5360q q G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536110z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536220z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.5360a a G f df σ⎡⎤=⎢⎥⎣⎦⎰()()()()()()13622w 012222412.5362202412.5d 12.5 0.5d d 1d d 4a a a a a a W f G f f f G f f G f f G f f G f f f ⎡⎤=⎢⎥⎣⎦⎡⎤⎛⎫⎛⎫=+++⎢⎥ ⎪ ⎪⎝⎭⎢⎥⎝⎭⎣⎦⎰⎰⎰⎰⎰其中，()()()()()0.5 0.52 2441 412.512.512.536f f f W f f f f <<⎧⎪⎪<<⎪=⎨<<⎪⎪<<⎪⎩()aw w 020lg L a =2）见1）中计算式3）()()()22004q q q G f G G n n uωπ== ()()()1122363622d ~00d d dd f q q f q f H f G f f G f f q σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122d 12f q f γλπ⎡⎤=⎢⎥∆⎣⎦ ()()()1122363622d //~00d d dd F Gq q F G q F H f G f f G f f Gq σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122222d 1412F f Gq g λζλμπγ⎡⎤⎛⎫-+⎢⎥⎪+⎢⎥⎝⎭=⎢⎥∆⎢⎥⎢⎥⎣⎦其余见1）中计算式2.程序清单1）f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi*fs;lambda=omega./omega0;lambda_s=omega./omega_s;%计算并绘制幅频特性delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;figure(1)loglog(f,z1_q,'b-',0:0.1:100,(0:0.1:100)./(0:0.1:100),'g-');axis([0.1,100,0.1,10]);grid on ;title('z_1~q 的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_1}{q}|$$','Interpreter','latex');gtext('0:1');figure(2)loglog(f,z2_z1,'b-',0:0.1:1,(0:0.1:1)./(0:0.1:1),'g-',1:0.1:100,1./...(1:0.1:100),'g-');axis([0.1,100,0.01,10]);grid on;title('z_2~z_1的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_2}{z_1}|$$','Interpreter','latex'); gtext('0:1');gtext('-1:1');figure(3)loglog(f,p_z2,'b-');axis([0.1,100,0.01,10]);grid on;title('p~z_2的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{p}{z_2}|$$','Interpreter','latex');%计算并绘制均方根谱sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q; figure(4)loglog(f,sqrt_Gz1,'b-');grid on;title('车轮位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''1}(f)}$','Interpreter','latex'); figure(5)loglog(f,sqrt_Gz2,'b-');grid on;title('车身位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''2}(f)}$','Interpreter','latex'); figure(6)loglog(f,sqrt_Ga,'b-');grid on;title('传至人体的位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_a(f)}$','Interpreter','latex');%计算其它值sigma_q=sqrt(trapz(f,sqrt_Gq.^2));%路面不平度加速度均方根值sigma_z1=sqrt(trapz(f,sqrt_Gz1.^2));%车轮加速度均方根值sigma_z2=sqrt(trapz(f,sqrt_Gz2.^2));%车身加速度均方根值sigma_a=sqrt(trapz(f,sqrt_Ga.^2));%传至人体的加速度均方根值for i=1:Nif f(i)<=2W(i)=0.5;elseif f(i)<=4W(i)=f(i)/4;elseif f(i)<=12.5W(i)=1;elseW(i)=12.5/f(i);endendaw=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));%加权加速度均方根值a0=10^(-6);Law=20*log10(aw/a0);%加权振级format shortdisp('路面不平度加速度均方根值=');disp(sigma_q);disp('车轮加速度均方根值=');disp(sigma_z1);disp('车身加速度均方根值=');disp(sigma_z2);disp('传至人体的加速度均方根值=');disp(sigma_a);disp('加权加速度均方根值=');disp(aw);disp('加权振级=');disp(Law);2）%随fs变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=1.5:0.01:6;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;lambda=omega./omega0;for i=1:length(fs)omega_s=2*pi.*fs(i);lambda_s=omega./omega_s;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(fs,aw);grid ontitle('a_w随f_s的变化曲线');xlabel('f_s/Hz');ylabel('a_w/m·s^-^2');figure(2)plot(fs,Law);grid ontitle('L_aw随f_s的变化曲线');xlabel('f_s/Hz');ylabel('L_a_w/dB');%随zeta_s变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_ss=0.125:0.001:0.5;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi.*fs;lambda=omega./omega0;lambda_s=omega./omega_s;for i=1:length(zeta_ss)zeta_s=zeta_ss(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(zeta_ss,aw);grid ontitle('a_w随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('a_w/m·s^-^2');figure(2)plot(zeta_ss,Law);grid ontitle('L_aw随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('L_a_w/dB');3）%随f0变化f0=0.25:0.01:3;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;lambda_s=omega./omega_s;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(f0)omega0=2*pi*f0(i);lambda=omega./omega0;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(f0==1.5);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(f0,sm_z2,f0,sm_fd,f0,sm_FdG);axis([0.25,3,-25,10]);grid on;title('各响应量均方根值随f_0变化的曲线');xlabel('f_0/Hz');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随zeta变化f0=1.5;zeta0=0.125:0.001:0.5;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(zeta0)zeta=zeta0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(zeta0==0.25);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(zeta0,sm_z2,zeta0,sm_fd,zeta0,sm_FdG);axis([0.125,0.5,-4,4]);grid on;title('各响应量均方根值随\zeta变化的曲线');xlabel('\zeta');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');%随gamma变化f0=1.5;zeta=0.25;gamma0=4.5:0.1:18;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(gamma0)gamma=gamma0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(gamma0==9);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(gamma0,sm_z2,gamma0,sm_fd,gamma0,sm_FdG);axis([4.5,18,-6,6]);grid on;title('各响应量均方根值随\gamma变化的曲线');xlabel('\gamma');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随mu变化f0=1.5;zeta=0.25;gamma=9;mu0=5:0.1:20;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(mu0)mu=mu0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(mu0==10);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(mu0,sm_z2,mu0,sm_fd,mu0,sm_FdG);axis([5,20,-2,2]);grid on;title('各响应量均方根值随\mu变化的曲线');xlabel('\mu');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB');legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');3.结果分析1）路面不平度加速度均方根值=0.3523车轮加速度均方根值=0.2391车身加速度均方根值=0.0168传至人体的加速度均方根值=0.0161加权加速度均方根值=0.0100加权振级=80.0287分析：根据课本中表6-2，a w=0.0100<0.315，L aw=80.0287<110，故乘客没有不舒适。
汽车理论课后习题答案及MATLAB编程
汽车理论课后习题答案—附MATLAB编程教材：汽车理论（第5版）清华大学余志生主编机械工业出版社出版本文档包含第三章 4.3 5.11 6.2 6.4 6.5所要求的MATLAB编程其它习题答案：第一章：1.1 1.2 1.7 1.8第二章：2.2 2.3 2.4第三章第四章：4.1 4.2 4.3 4.5第五章：5.1 5.8 5.10 5.11 5.12 5.14 5.17第六章：6.1~6.5见本作者另一百度文库文档。
获取方法：点击左上角本作者的昵称，进入本作者的文库，《汽车理论课后习题答案（含MATLAB编程）》第三章：m=3880;%总质量g=9.8;%重力加速度r=0.367;%车轮半径eta_t=0.85;%传动系机械效率f=0.013;%滚动阻力系数CDA=2.77;%空气阻力系数*迎风面积i0=[5.175.435.836.176.33];%主减速器传动比If=0.218;%飞轮转动惯量Iw1=1.798;Iw2=3.598;%两前轮/四后轮转动惯量Iw=Iw1+Iw2;ig4=[6.09 3.091.711];%4挡变速器传动比ig5=[5.56 2.7691.64410.793];%5挡变速器传动比n=600:1:4000;%发动机转速Tq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;%4挡发动机for i=1:5Ft1=Tq*ig4(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig4(2)*i0(i)/r;Ft3=Tq*ig4(3)*i0(i)/r;Ft4=Tq*ig4(4)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig4(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig4(2));u3=0.377*r*n/(i0(i)*ig4(3));u4=0.377*r*n/(i0(i)*ig4(4));F1=m*g*f+CDA*u1.^2/21.15;%各档对应转速下的行驶阻力F2=m*g*f+CDA*u2.^2/21.15;F3=m*g*f+CDA*u3.^2/21.15;F4=m*g*f+CDA*u4.^2/21.15;delta1=1+(Iw+If.*ig4(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig4(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig4(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig4(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;figure(1)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=96.6;x1=[];x2=[];x3=[];x4=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];elseif u4(k)<=u4maxx3=[k];endendq2=max(x2);q3=max(x3);ua4=u4(q2:q3);a4=b4(q2:q3);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s=[s2s3s4];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs4(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs4=');disp(Qs4);tt=[48.423546.351141.555733.379630.5564];figure(2)plot(Qs4,tt,'*')hold onxi=25:0.001:31;pp=interp1(Qs4,tt,xi,'cubic');plot(xi,pp)title('燃油经济性-加速时间曲线');xlabel('燃油经济性--等速百公里油耗(L/100km)');ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')%5挡发动机for i=1:5Ft1=Tq*ig5(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig5(2)*i0(i)/r;Ft3=Tq*ig5(3)*i0(i)/r;Ft4=Tq*ig5(4)*i0(i)/r;Ft5=Tq*ig5(5)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig5(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig5(2));u3=0.377*r*n/(i0(i)*ig5(3));u4=0.377*r*n/(i0(i)*ig5(4));u5=0.377*r*n/(i0(i)*ig5(5));F1=m*g*f+CDA*u1.^2/21.15;%各档对应转速下的行驶阻力F2=m*g*f+CDA*u2.^2/21.15;F3=m*g*f+CDA*u3.^2/21.15;F4=m*g*f+CDA*u4.^2/21.15;F5=m*g*f+CDA*u5.^2/21.15;delta1=1+(Iw+If.*ig5(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig5(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig5(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig5(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta5=1+(Iw+If.*ig5(5).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);a5=(Ft5-F5)/(delta5*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;b5=1./a5;figure(3)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4,u5,b5);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext('1/a5')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=max(u4);u5min=u4max;u5max=96.6;x1=[];x2=[];x3=[];x4=[];x5=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];endendq2=max(x2);ua4=u4(q2:y);a4=b4(q2:y);for l=1:3401if u5(l)<=u5minx3=[l];elseif u5(l)<=u5maxx4=[l];endendq2=max(x2);q3=max(x3);q4=max(x4);ua5=u5(q3:q4);a5=b5(q3:q4);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s5=trapz(ua5,a5);%五挡运行时间s=[s2s3s4s5];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;Pe5=Ft5.*u5./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs5(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs5=');disp(Qs5);tt=[93.973058.283444.678445.366744.9793];figure(2)plot(Qs5,tt,'*')hold onxi=30:0.001:100;pp=interp1(Qs5,tt,xi,'cubic');plot(xi,pp)%title('燃油经济性-加速时间曲线');%xlabel('燃油经济性--等速百公里油耗(L/100km)');%ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')运行结果：t=48.4235t=46.3511t=41.5557t=33.3796t=30.556490％负荷Pe=55.5647Qs4=25.103126.365628.307829.958730.7356图1加速度倒数-速度曲线图2燃油经济性-加速时间曲线4.3%空载(no load)—1；满载(full load)—2m1=4080;m2=9290;hg1=0.845;hg2=1.170;L=3.950;a1=2.100;a2=2.950;b1=L-a1;b2=L-a2;beta=0.38;%利用附着系数与制动强度的关系曲线z=0:0.01:1;phi_f1=L*beta.*z./(b1+z*hg1);%前轮利用附着系数(空载)phi_r1=L*(1-beta).*z./(a1-z*hg1);%空载时后轮利用附着系数(空载)phi_f2=L*beta.*z./(b2+z*hg2);%前轮利用附着系数(满载)phi_r2=L*(1-beta).*z./(a2-z*hg2);%后轮利用附着系数(满载)phi=z;%最理想情况figure(1);plot(z,phi_f1,'g--',z,phi_f2,'g-',z,phi_r1,'b--',z,phi_r2,'b-',z,phi,'k-'); axis([0,1,0,1]);grid on;box off;axis square;title('利用附着系数与制动强度的关系曲线');xlabel('制动强度z/g');ylabel('利用附着系数{\phi}');%ECE法规要求界限hold on;z1=0.2:0.01:0.8;z2=0.15:0.01:0.3;z3=0.3:0.01:1;phi1=(z1+0.07)/0.85;phi2=z2-0.08;phi3=z2+0.08;phi4=(z3-0.02)/0.74;plot(z1,phi1,'r-.',z2,phi2,'r-.',z2,phi3,'r-.',z3,phi4,'r-.'); legend('{\phi}_{f}(空载)','{\phi}_{f}(满载)','{\phi}_{r}(空载)',...'{\phi}_{r}(满载)','{\phi}=z','ECE法规要求界限');legend('Location','northwest');%制动效率曲线Ef1=z./phi_f1*100;Er1=z./phi_r1*100;Ef2=z./phi_f2*100;Er2=z./phi_r2*100;figure(2);plot(phi_f1,Ef1,'k',phi_r1,Er1,'k',phi_f2,Ef2,'b',phi_r2,Er2,'b'); axis([0,1,0,100]);grid on;box off;axis square;title('制动效率曲线');xlabel('附着系数{\phi}');ylabel('制动效率(%)');text(0.25,90,'E_f');text(0.6,88,'满载');text(0.8,90,'E_r');text(0.7,65,'空载');text(0.72,72,'E_r');程序运行结果如下：图1图2 5.11L=3.048;K=0.0024;i=20;ua=0:0.01:120;f=1/L.*(1./(3.6./ua+K.*ua./3.6));syms uF(u)=(1/L.*(1./(1./u+K.*u)))./i;plot(ua,f,'b-');xlabel('u_a(km/h)');ylabel('\omega_r/\delta)_s');F(22.35)运行结果：ans=0.1667386094155036.2f=0.1:0.1:100;Gq1=5.12*10^(-5)./f.^2;Gq2=2.02*10^(-3).*f./f;Gq3=7.98*10^(-2).*f.^2;figure(1);loglog(f,Gq1,'b-');box off;axis square;title('位移功率谱密度'); xlabel('f/Hz');ylabel('G_q(f)/(m^2s)');figure(2);loglog(f,Gq2,'b-');box off;axis square;title('速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''(f)/(m^2/s)');figure(3);loglog(f,Gq3,'b-');box off;axis square;title('加速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''_''(f)/(m^2/s^3)');运行结果如下：6.4f1=0.1:0.01:1;f2=1:0.01:10;G1=0.010106.*f1.^2;G2=0.010106.*f2./f2;loglog(f1,G1,'b-');hold onloglog(f2,G2,'b-');title('车身加速度的功率谱密度'); xlabel('f/Hz');ylabel('G_z_''_''(f)/m^-^1'); axis([0.1,10,0.0001,0.1]);运行结果如下：6.51.计算说明1）①幅频特性()122222114z q λζλγ⎡⎤-+⎢=⎢⎥∆⎣⎦，其中0ωλω=，()2222222111141λλγζλγλμμ⎡⎤⎡⎤⎛⎫⎛⎫∆=-+--+-⎢⎥⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎣⎦⎣⎦；()12222222211414z z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦；()1222222221414s s s s s p z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦，其中s s ωλω=；②均方根谱()()221~j 2z q H f f ωππ=；()()222~j 2z qH f f ωππ== ；()()2221~1j 2p qz z H f f z qωππ== ；③其他值()0.5360q q G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536110z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536220z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.5360a a G f df σ⎡⎤=⎢⎥⎣⎦⎰()()()()()()13622w 012222412.5362202412.5d 12.5 0.5d d 1d d 4a a a a a a W f G f f f G f f G f f G f f G f f f ⎡⎤=⎢⎥⎣⎦⎡⎤⎛⎫⎛⎫=+++⎢⎥ ⎪ ⎪⎝⎭⎢⎥⎝⎭⎣⎦⎰⎰⎰⎰⎰其中，()()()()()0.5 0.52 2441 412.512.512.536f f f W f f f f <<⎧⎪⎪<<⎪=⎨<<⎪⎪<<⎪⎩()aw w 020lg L a =2）见1）中计算式3）()()()22004q q q G f G G n n uωπ== ()()()1122363622d ~00d d dd f q q f q f H f G f f G f f q σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122d 12f q f γλπ⎡⎤=⎢⎥∆⎣⎦ ()()()1122363622d //~00d d dd F Gq q F G q F H f G f f G f f Gq σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122222d 1412F f Gq g λζλμπγ⎡⎤⎛⎫-+⎢⎥⎪+⎢⎥⎝⎭=⎢⎥∆⎢⎥⎢⎥⎣⎦其余见1）中计算式2.程序清单1）f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi*fs;lambda=omega./omega0;lambda_s=omega./omega_s;%计算并绘制幅频特性delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;figure(1)loglog(f,z1_q,'b-',0:0.1:100,(0:0.1:100)./(0:0.1:100),'g-');axis([0.1,100,0.1,10]);grid on ;title('z_1~q 的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_1}{q}|$$','Interpreter','latex');gtext('0:1');figure(2)loglog(f,z2_z1,'b-',0:0.1:1,(0:0.1:1)./(0:0.1:1),'g-',1:0.1:100,1./...(1:0.1:100),'g-');axis([0.1,100,0.01,10]);grid on;title('z_2~z_1的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_2}{z_1}|$$','Interpreter','latex'); gtext('0:1');gtext('-1:1');figure(3)loglog(f,p_z2,'b-');axis([0.1,100,0.01,10]);grid on;title('p~z_2的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{p}{z_2}|$$','Interpreter','latex');%计算并绘制均方根谱sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q; figure(4)loglog(f,sqrt_Gz1,'b-');grid on;title('车轮位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''1}(f)}$','Interpreter','latex'); figure(5)loglog(f,sqrt_Gz2,'b-');grid on;title('车身位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''2}(f)}$','Interpreter','latex'); figure(6)loglog(f,sqrt_Ga,'b-');grid on;title('传至人体的位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_a(f)}$','Interpreter','latex');%计算其它值sigma_q=sqrt(trapz(f,sqrt_Gq.^2));%路面不平度加速度均方根值sigma_z1=sqrt(trapz(f,sqrt_Gz1.^2));%车轮加速度均方根值sigma_z2=sqrt(trapz(f,sqrt_Gz2.^2));%车身加速度均方根值sigma_a=sqrt(trapz(f,sqrt_Ga.^2));%传至人体的加速度均方根值for i=1:Nif f(i)<=2W(i)=0.5;elseif f(i)<=4W(i)=f(i)/4;elseif f(i)<=12.5W(i)=1;elseW(i)=12.5/f(i);endendaw=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));%加权加速度均方根值a0=10^(-6);Law=20*log10(aw/a0);%加权振级format shortdisp('路面不平度加速度均方根值=');disp(sigma_q);disp('车轮加速度均方根值=');disp(sigma_z1);disp('车身加速度均方根值=');disp(sigma_z2);disp('传至人体的加速度均方根值=');disp(sigma_a);disp('加权加速度均方根值=');disp(aw);disp('加权振级=');disp(Law);2）%随fs变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=1.5:0.01:6;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;lambda=omega./omega0;for i=1:length(fs)omega_s=2*pi.*fs(i);lambda_s=omega./omega_s;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(fs,aw);grid ontitle('a_w随f_s的变化曲线');xlabel('f_s/Hz');ylabel('a_w/m·s^-^2');figure(2)plot(fs,Law);grid ontitle('L_aw随f_s的变化曲线');xlabel('f_s/Hz');ylabel('L_a_w/dB');%随zeta_s变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_ss=0.125:0.001:0.5;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi.*fs;lambda=omega./omega0;lambda_s=omega./omega_s;for i=1:length(zeta_ss)zeta_s=zeta_ss(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(zeta_ss,aw);grid ontitle('a_w随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('a_w/m·s^-^2');figure(2)plot(zeta_ss,Law);grid ontitle('L_aw随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('L_a_w/dB');3）%随f0变化f0=0.25:0.01:3;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;lambda_s=omega./omega_s;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(f0)omega0=2*pi*f0(i);lambda=omega./omega0;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(f0==1.5);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(f0,sm_z2,f0,sm_fd,f0,sm_FdG);axis([0.25,3,-25,10]);grid on;title('各响应量均方根值随f_0变化的曲线');xlabel('f_0/Hz');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随zeta变化f0=1.5;zeta0=0.125:0.001:0.5;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(zeta0)zeta=zeta0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(zeta0==0.25);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(zeta0,sm_z2,zeta0,sm_fd,zeta0,sm_FdG);axis([0.125,0.5,-4,4]);grid on;title('各响应量均方根值随\zeta变化的曲线');xlabel('\zeta');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');%随gamma变化f0=1.5;zeta=0.25;gamma0=4.5:0.1:18;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(gamma0)gamma=gamma0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(gamma0==9);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(gamma0,sm_z2,gamma0,sm_fd,gamma0,sm_FdG);axis([4.5,18,-6,6]);grid on;title('各响应量均方根值随\gamma变化的曲线');xlabel('\gamma');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随mu变化f0=1.5;zeta=0.25;gamma=9;mu0=5:0.1:20;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(mu0)mu=mu0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(mu0==10);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(mu0,sm_z2,mu0,sm_fd,mu0,sm_FdG);axis([5,20,-2,2]);grid on;title('各响应量均方根值随\mu变化的曲线');xlabel('\mu');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB');legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');3.结果分析1）路面不平度加速度均方根值=0.3523车轮加速度均方根值=0.2391车身加速度均方根值=0.0168传至人体的加速度均方根值=0.0161加权加速度均方根值=0.0100加权振级=80.0287分析：根据课本中表6-2，a w=0.0100<0.315，L aw=80.0287<110，故乘客没有不舒适。
最新汽车理论-余志生-课后习题答案详解(全)
第一章1.1、试说明轮胎滚动阻力的定义、产生机理和作用形式?答：1）定义：汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动阻力。
2）产生机理：由于轮胎内部摩擦产生弹性轮胎在硬支撑路面上行驶时加载变形曲线和卸载变形曲线不重合会有能量损失，即弹性物质的迟滞损失。
这种迟滞损失表现为一种阻力偶。
当车轮不滚动时，地面对车轮的法向反作用力的分布是前后对称的；当车轮滚动时，由于弹性迟滞现象，处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力，这样，地面法向反作用力的分布前后不对称，而使他们的合力Fa 相对于法线前移一个距离a, 它随弹性迟滞损失的增大而变大。
即滚动时有滚动阻力偶矩a F T z f = 阻碍车轮滚动。
3）作用形式：滚动阻力 fw F f = rT F f f =（f 为滚动阻力系数）1.2、滚动阻力系数与哪些因素有关？提示：滚动阻力系数与路面种类、行驶车速以及轮胎的构造、材料、气压等有关。
1.3、确定一轻型货车的动力性能（货车可装用4档或5档变速器，任选其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车的最高车速、最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速倒数曲线，用图解积分法求汽车有Ⅱ档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用Ⅱ档起步加速至70km/h 的加速时间。
轻型货车的有关数据：汽油发动机使用外特性的Tq —n 曲线的拟合公式为432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-= 式中, Tq 为发功机转矩(N ·m)；n 为发动机转速(r ／min)。
发动机的最低转速nmin=600r/min ，最高转速nmax=4000 r ／min 装载质量 2000kg 整车整备质量 1800kg 总质量 3880 kg 车轮半径 0.367 m传动系机械效率 ηт＝0.85 波动阻力系数 f ＝0.013 空气阻力系数×迎风面积 CDA ＝2.772m主减速器传动比 i 0＝5.83 飞轮转功惯量 If ＝0.218kg ·2m 二前轮转动惯量 Iw1＝1.798kg ·2m四后轮转功惯量 Iw2＝3.598kg ·2m变速器传动比 ig (数据如下表)轴距 L ＝3.2m 质心至前铀距离(满载) α＝1.947m质心高(满载) h g ＝0.9m 解答：1）（取四档为例） 由uF n u n Tq Tq F t t →⇒⎪⎭⎪⎬⎫→→→即ri i T F To g q t η=432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-=og i i rn u 377.0=行驶阻力为w fF F +：215.21a D w f U A C Gf F F +=+2131.0312.494aU +=由计算机作图有※本题也可采用描点法做图： 由发动机转速在m in/600n minr =，m in/4000n maxr =，取六个点分别代入公式：……………………………… 2）⑴最高车速： 有w f t F F F +=⇒2131.0312.494a t U F +=分别代入a U 和t F 公式：2)09.6*83.53697.0*377.0(131.0312.494367.085.0*83.5*9.6*n T q +=把q T 的拟和公式也代入可得： n>4000 而4000m ax=n r/min∴93.9483.5*0.14000*367.0*377.0max ==U Km/h⑵最大爬坡度：挂Ⅰ档时速度慢，Fw 可忽略：⇒)(max w f t i F F F F +-=⇒GfF Gi t -=max⇒013.08.9*388014400max max -=-=f G F i t =0.366（3）克服该坡度时相应的附着率zxF F =ϕ忽略空气阻力和滚动阻力得：6.0947.12.3*366.0/=====a il l a i F Fi z ϕ 3）①绘制汽车行驶加速倒数曲线（已装货）：40.0626)(1f D g du dt a -==δ（GFwFt D -=为动力因素）Ⅱ时，22022111r i i I m r ImT g f wηδ++=∑2222367.085.0*83.5*09.3*218.038001367.0598.3798.1380011+++==1.128ri i T F To g q t η=432)1000(8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-=215.21a D w U A C F =由以上关系可由计算机作出图为：②用计算机求汽车用Ⅳ档起步加速至70km/h 的加速时间。
汽车理论课后习题答案部分
第一章 汽车的动力性1.1试说明轮胎滚动阻力的定义，产生机理和作用形式。
答：车轮滚动时,由于车轮的弹性变形、路面变形和车辙摩擦等原因所产生的阻碍汽车行驶的力称为轮胎滚动阻力。
产生机理和作用形式：(1)弹性轮胎在硬路面上滚动时，轮胎的变形是主要的，由于轮胎有内部摩擦，产生弹性迟滞损失，使轮胎变形时对它做的功不能全部回收。
由于弹性迟滞，地面对车轮的法向作用力并不是前后对称的，这样形成的合力z F 并不沿车轮中心（向车轮前进方向偏移a ）。
如果将法向反作用力平移至与通过车轮中心的垂线重合，则有一附加的滚动阻力偶矩f z T F a =⋅。
为克服该滚动阻力偶矩，需要在车轮中心加一推力P F 与地面切向反作用力构成一力偶矩。
（2）轮胎在松软路面上滚动时，由于车轮使地面变形下陷，在车轮前方实际形成了具有一定坡度的斜面，对车轮前进产生阻力。
（3）轮胎在松软地面滚动时，轮辙摩擦会引起附加阻力。
（4）车轮行驶在不平路面上时，引起车身振荡、减振器压缩和伸长时做功，也是滚动阻力的作用形式。
1.2滚动阻力系数与哪些因素有关？答：滚动阻力系数与路面的种类、行驶车速以及轮胎的构造、材料和气压有关。
这些因素对滚动阻力系数的具体影响参考课本P9。
第三章 汽车动力装置参数的选定3.1改变1.3题中轻型货车的主减速器传动比，做出0i 为5.17、5.43、5.83、6.17、6.33时的燃油经济性—加速时间曲线，讨论不同0i 值对汽车性能的影响。
解：加速时间的结算思路与方法：在算加速时间的时候，关键是要知道在加速的过程中，汽车的行驶加速度dudt随着车速的变化。
由汽车行驶方程式：0221.15tq g TD a T i i C A duGf Gi u m rdtηδ=+++，可以的到： 021[()]21.15tq g TD a T i i C A du Gf u dt m r ηδ=-+（0i F =） 由于对于不同的变速器档位，车速a u 与发动机转速n 的对应关系不同，所以要针对不同的变速器档位，求出加速度a 随着车速a u 变化的关系。
汽车理论第五版课后习题答案
第一章汽车的动力性试说明轮胎滚动阻力的定义，产生机理和作用形式。
定义:汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动祖力。
2)产生机理:由于轮胎内部摩擦产生弹性轮胎在硬支律路面上行驶时加载变形曲线和卸载变形曲线不重合会有能全损失，即弹性物质的迟滞损失。
这种迟滞损失表现为一种阻力偶。
当车轮不滚动时，地面对丰轮的法向反作用力的分布是前后对称的;当车轮滚动时，由于弹性阻滞现象，处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力，这样，地面法向反作用力的分布前后不对称，而使他们的合力F Z相对于法线前移一个距离a，它随弹性迟滞损失的增大变大。
即滚动时有滚动阻力偶T f = F Z.a阻碍车轮滚动。
3]作用形式: T f = Wf，T f = T f/r滚动阻力系数与哪些因素有关答：滚动阻力系数与路面的种类、行驶车速以及轮胎的构造、材料和气压有关。
=+由计算机作图有：空车、满载时汽车动力性有无变化为什么答：动力性会发生变化。
因为满载时汽车的质量会增大，重心的位置也会发生改变。
质量增大，滚动阻力、坡度阻力和加速阻力都会增大，加速时间会增加，最高车速降低。
重心位置的改变会影响车轮附着率，从而影响最大爬坡度。
如何选择汽车发动机功率答：发动机功率的选择常先从保证汽车预期的最高车速来初步确定。
若给出了期望的最高车速，选择的发动机功率应大体等于，但不小于以最高车速行驶时的行驶阻力功率之和。
发动机的最大功率但也不宜过大，否则会因发动机负荷率偏低影响汽车的燃油经济性。
在实际工作中，还利用现有汽车统计数据初步估计汽车比功率来确定发动机应有功率。
不少国家还对车辆应有的最小比功率作出规定，以保证路上行驶车辆的动力性不低于一定水平，防止某些性能差的车辆阻碍车流。
超车时该不该换入低一挡的排挡答：超车时排挡的选择，应该使车辆在最短的时间内加速到较高的车速，所以是否应该换入低一挡的排挡应该由汽车的加速度倒数曲线决定。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第 五 章 汽车的操纵稳定性5．1一轿车(每个)前轮胎的侧偏刚度为-50176N ／rad 、外倾刚度为-7665N ／rad 。
若轿车向左转弯，将使两前轮均产生正的外倾角，其大小为40。
设侧偏刚度与外倾刚度均不受左、右轮载荷转移的影响．试求由外倾角引起的前轮侧偏角。
答： 由题意：F Y =k α+k γγ=0故由外倾角引起的前轮侧偏角：α=- k γγ/k=-7665⨯4/-50176=0.61105．2 6450轻型客车在试验中发现过多转向和中性转向现象，工程师们在前悬架上加装前横向稳定杆以提高前悬架的侧倾角刚度，结果汽车的转向特性变为不足转向。
试分析其理论根据(要求有必要的公式和曲线)。
答： 稳定性系数：⎪⎪⎭⎫ ⎝⎛-=122k b k a L mK1k 、2k 变化，原来K ≤0,现在K>0,即变为不足转向。
5．3汽车的稳态响应有哪几种类型?表征稳态响应的具体参数有哪些?它们彼此之间的关系如何(要求有必要的公式和曲线)？答： 汽车稳态响应有三种类型 ：中性转向、不足转向、过多转向。
几个表征稳态转向的参数：1．前后轮侧偏角绝对值之差(α1-α2);2. 转向半径的比R/R 0;3．静态储备系数S.M.彼此之间的关系见参考书公式（5-13）（5-16）（5-17）。
5．4举出三种表示汽车稳态转向特性的方法，并说明汽车重心前后位置和内、外轮负荷转移如何影响稳态转向特性？答：方法：1.α1-α2 >0时为不足转向，α1-α2 =0时 为中性转向，α1-α2 <0时为过多转向；2. R/R0>1时为不足转向，R/R0=1时为中性转向，R/R0<1时为过多转向；3 .S.M.>0时为不足转向，S.M.=0时为中性转向，S.M.<0时为过多转向。
汽车重心前后位置和内、外轮负荷转移使得 汽车质心至前后轴距离a 、b 发生变化，K 也发生变化。
5．5汽车转弯时车轮行驶阻力是否与直线行驶时一样？答：否，因转弯时车轮受到的侧偏力，轮胎产生侧偏现象，行驶阻力不一样。
5．6主销内倾角和后倾角的功能有何不同?答：主销外倾角可以产生回正力矩，保证汽车直线行驶；主销内倾角除产生回正力矩外，还有使得转向轻便的功能。
5．7横向稳定杆起什么作用?为什么有的车装在前恳架，有的装在后悬架，有的前后都装？答： 横向稳定杆用以提高悬架的侧倾角刚度。
装在前悬架是使汽车稳定性因数K 变大，装在后悬架使K 变小，前后悬架都装则使前后悬架侧倾角刚度同时增大。
5．8某种汽车的质心位置、轴距和前后轮胎的型号已定。
按照二自由度操纵稳定性模型，其稳态转向特性为过多转向，试找出五种改善其转向特性的方法。
答： 即要K 变大，可在前悬架上装横向稳定杆，后悬架不变；前悬架不变，减小后悬架侧倾角刚度；同时在前后悬架装横向稳定杆，但保证a/k2-b/k1变大；同时减小前后悬架侧倾角刚度，但保证a/k2-b/k1变大；增大汽车质量。
5．9汽车空载和满载是否具有相同的操纵稳定性?答： 否，m 不同，空载时的汽车m 小于满载时的 m ，即满载时的K 更大，操纵稳定性更好。
5．10试用有关计算公式说明汽车质心位置对主要描述和评价汽车操纵稳定性、稳态响应指标的影响。
答：稳定性系数公式：⎪⎪⎭⎫ ⎝⎛-=122k b k a L mK汽车质心位置变化，则a 、b 变化，即K 也随之改变。
5．11二自由度轿车模型的有关参数如下：总质量 m ＝1818．2kg绕z o 轴转动惯量 23885z I kg m =⋅轴距 L=3.048m质心至前轴距离 a=1.463m质心至后轴距离 b=1.585m前轮总侧偏刚度 162618/k N rad =后轮总侧偏刚度 2110185k =- /N rad转向系总传动比 i=20试求：1)稳定性因数K 、特征车速ch u 。
2)稳态横摆角速度增益曲线r ωδ)s ---a u 车速u=22.35m/s 时的转向灵敏度r sw ωδ。
3)静态储备系数S.M.，侧向加速度为0．4g 时的前、后轮侧偏角绝对值之差12a a -与转弯半径的比值R/R 0(R 0=15m)。
4)车速u=30.56m/s,瞬态响应的横摆角速度波动的固有(圆)频率0ω、阻尼比ς、反应时间τ与峰值反应时间ε。
提示： 1) 稳定性系数：⎪⎪⎭⎫ ⎝⎛-=122k b k a L mK特征车速Ku ch 1= 2) 转向灵敏度21KuL u s r +=⎪⎭⎫δω 3) ()211αα-=La K y ⇒ 21αα-， 4) 固有圆频率 m c '=0ω 阻尼比m h'=02ωξ 反应时间ωτΦ-=峰值反应时间 ωξωωεΦ-=0arctg5．12稳态响应中横摆角速度增益达到最大值时的车速称为特征车速ch u 。
证明： 特征车速ch uL 中性转向汽车横摆角速度增益的一半。
答： 转向灵敏度21KuL u s r +=⎪⎭⎫δω 特征车速K u ch 1=⇒ s r ⎪⎭⎫δωL u =,中性转向时s r ⎪⎭⎫δωLu = 得证。
5．13测定汽车稳态转向特性常用两种方法，一为固定转向盘转角法，并以0/R R —y a 曲线来表示汽车的转向特性(见第五章第三节二)；另一为固定圆周法。
试验时在场地上画一圆，驾驶员以低速沿圆周行驶，记录转向盘转角0sw δ，，然后驾驶员控制转向盘使汽车始终在圆周上以低速持续加速行驶。
随着车速的提高，转向盘转角sw δ（一般）将随之加大。
记录下sw δ角，并以0sw y sw a δδ-曲线来评价汽车的转向特性，试证201sw sw Ku δδ=+，说明如何根据0sw y sw a δδ-曲线来判断汽车的转向特性。
证明：转向半径()()222111Ku L Ku Ku L u u u R r +=+=⎪⎭⎫ ⎝⎛+==δδω0R 5．14习题图4是滑柱连杆式独立悬架(常称为Mc Pherson strut suspnsion)示意图。
试证：1)R.C.为侧倾中心。
2)悬架的侧倾角刚度为2()s mp K r k nϕ=，式中s k 为一个弹簧的(线)刚度。
提示： 1)画出地面对于车厢的瞬时转动中心，即为侧倾中心R.C.2)证明参考书P135-1365．15试求计算稳态响应质心侧偏角增益)s βδ的公式，并求题5．11中轿车在u=31．3m/s(70 mile ／h)、0.4y a g =时的质心侧偏角。
计算u=31．3m ／s 时的瞬态响应峰值反应时间ε和轿车的汽车因数T.B.值。
提示：将方程组（5-9）两式联立，v =0, r ω =0,消去r ω⇒)sβδ 5．16为什么有些小轿车后轮也设计有安装前束角和外倾角?答：因为轿车后轮安装前束角和外倾角是为提高操纵稳定性。
5．18转向盘力特性与哪些因素有关，试分析之。
答： 转向盘力特性决定于以下因素：转向器传动比及其变化规律、转向器效率、动力转向器的转向盘操作力特性、转向杆系传动比、转向杆系效率、由悬架导向杆系决定的主销位置、轮胎上的载荷、轮胎气压、轮胎力学特性、地面附着条件、转向盘转动惯性、转向柱摩擦阻力及汽车整体动力学特性。
5．19地面作用于轮胎的切向反作用力是如何控制转向特性的?提示：参考书P152-155。
第六章 汽车的平顺性6.l 、设通过座椅支承面传至人体垂直加速度的谱密度为一白噪声，Ga ( f )=0.132m -⋅s 。
求在0.5～80HZ 频率范围加权加速度均方根值a w 和加权振级Law，并由表6-2查出相应人的主观感觉。
答：21 805.02])()([dffGfWaaw⎰⋅=查173P图知：人的主观感觉为极不舒适。
6.2、设车速u=20m／s ，路面不平度系38q10*56.2)(G mn-=，参考空间频率no=0.1-1m。
画出路面垂直位移、速度和加速度)(Gqf、)(Gqf&、)(Gqf&&的谱图。
画图时要求用双对数坐标，选好坐标刻度值，并注明单位。
解：22822q20*1.0*10*56.2)()(GffunnGfq-==画出图形为：6.3、设车身－车轮二自由度汽车模型，其车身部分固有频率fo ＝2Hz。
它行驶在波长λ=5m 的水泥接缝路上，求引起车身部分共振时的车速un(km/h)。
该汽车车轮部分的固有频率 f t=10Hz，在砂石路上常用车速为30km/h。
问由于车轮部分共振时，车轮对路面作用的动载所形成的搓板路的波长λ=?答：①当激振力等于车辆固有频率时，发生共振，所以发生共振时的车速为：②搓板路的波长 ： m 65106.3/30==λ 6.4、设车身单质量系统的幅频 |z ／q |用双对数坐标表示时如习题图6所示。
路面输入谱与题6.2相同。
求车身加速度的谱密度)(f G z&&，画出其谱图，并计算0.1～10Hz 频率范围车身加速度的均方根值z &&σ。
答：① )()2()(02n G qz f f G q z ⋅⋅=π&&6.5、上机计算作业(报告应包括：题目、计算说明、程序清单、结果分析)。
车身－车轮双质最系统参数：f o=1.5Hz 、ζ=0.25、γ=9、μ=10。
“人体—座椅”系统参数：f s=3Hz 、ζs ＝0.25。
车速u =20m/s ，路面不平度系数G q (no)=2.56×-810，参考空间频率no=0.1-1m 。
计算时频率步长△f ＝0.2Hz ，计算频率点数N ＝180。
1）计算并画出幅频特性|z1／q |、|z2／z1|、|p ／z2|和均方根值谱)(G zf &&、)(G 2z f &&、)(G a f 谱图。
进—步计算q&&σ、1z &&σ、2z &&σ、a σ、w a 、aw L 值。
2)改变“人体—座椅”系统参数：f s=1.5～6Hz 、ζs ＝0.125～0.5。
分析w a 、aw L 值随f s 、ζs 的变化。
3)分别改变车身—车轮双质量系统参数：f o ＝0.25～3Hz 、ζ=0.125～0.5、γ=4.5～18、μ=5～20。
绘制z &&σ、σfd 、σFd/G 三个响应量均方根值随以上四个系统参数变化的曲线。
提示：本题可简单利用matlab 软件求出各数值，并作出相应的图。
6.6、设前、后车轮两个输入的双轴汽车模型行驶在随机输入的路面上，其质量分配系数ε＝1，前、后车身局部系统的固有频率均为f o ＝2Hz ，轴距L ＝2.5m 。
问引起车身俯仰角共振时的车速u a=? 相应随机路面输入的λ=?答：。