汽车运用工程汽车动力性的应用

1、汽车运用工程大作业 1 汽车动力性计算（加速性能、最高车速、动力因素加速时间和加速距离）已知：发动机外特性、装载质量、整备质量、总质量、车轮半径、传动系效率、滚动阻力系数、空气阻力系数×迎风面积、主减速 器速比飞轮转动惯量、两个前轮转动惯量、四个后轮转动惯量、变速器速比、轴距、质心至前轴距离、质心高度。例如，已知数据如下表:项目单位数值装载质量kg1995整备质量kg1700总质量kg3695车轮半径m0.37传动系效率 0.87+0.024=0.894滚动阻力系数 0.012-0.0024=0.0096空气阻力系数×迎风面积m22.60主减速器速比&#

2、160;5.83飞轮转动惯量kgm20.22两前轮转动惯量kgm21.80四后轮转动惯量kgm23.60轴距m3.40质心至前轴距离m2.0质心高度m0.9变速器速比变速器I档II档III档IV档V档四档6.123.111.691.00 五档3.552.771.651.000.83发动机外特性 Tq=a+a1(n/c)+a2(n/c)2+a2(n/c)3+a3(n/c)4参数aa1a2a2a3c数值-19.31296-16540.9-3.8510001)  绘制汽车驱动力与行驶阻力平衡图；2)  求汽车的最高车速、最大爬坡度；3)  绘制加速度倒数曲线4

3、)  用图解法或编程绘制汽车动力因数特性曲线5)  图解手工或编程绘制II档起步，加速至70km/h的车速时间曲线以及距离时间曲线。解:轻型货车的有关数据：汽油发动机使用外特性的Tq-n曲线的拟合公式为Tq=a+a1(n/c)+a2(n/c)*(n/c)+a3(n/c)*(n/c)*(n/c)+a4(n/c)*(n/c)*(n/c)*(n/c)式中，Tq为发动机转矩（Nm）;n为发动机转速（r/min）。（ km/h）其中：为车轮滚动半径，m; 汽车的牵引力： （ N ）汽车的空气阻力：（ N ）其中：为空气阻力系数，A为汽车迎风面积，m2。 汽车的滚动阻力：（ N ） 其

4、中：mg 为满载或空载汽车总重(N)，为滚动阻尼系数 汽车的行驶阻力之和：（ N ）Matlab程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序：五档车程序如下 n=600:10:4000;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;nmax=4000;G=m*g;ig=3.55 2.77 1.65 1.00 0.83;nT=0.8924;r=0.37;f=0.0096;CDA=2.60;i0=5.83;L=3.40;a=2.0;hg=0.9;If

5、=0.22;Iw1=1.80;Iw2=3.60;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;ua=0:5:120;Ff=G*f;Fw=CDA*ua.2/21.15;Fz=Ff+Fw;plot(ua1,Ft

6、1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw');zoom on;x,y=ginput(1);zoom off;disp('汽车最高车速=

7、9;);disp(x);disp('km/h');汽车最高车速=124km/h四档车程序如下：q=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;nmax=4000;G=m*g;ig=6.12 3.11 1.69 1.00 ;nT=0.894;r=0.37;f=0.0096;CDA=2.6;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.80;Iw2=3.60; Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*i

8、g(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua=0:5:120;Ff=G*f;Fw=CDA*ua.2/21.15;Fz=Ff+Fw;plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua,Fz);title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel(&

9、#39;Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ff+Fw');zoom on;x,y=ginput(1);zoom off;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;nmax=4000;G=m*g;ig=6.12 3.11 1.69 1.00 ;nT=0.894;r=0.37;

10、f=0.0096;CDA=2.6;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.80;Iw2=3.60; Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua=0:5:120;Ff=G*f;Fw=CDA*ua.2/21.15;Fz=Ff+Fw;plot(ua1

11、,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua,Fz);title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ff+Fw');zoom on;x,y=ginput(1);zoom off;disp('汽车最高车速=');disp(x);disp('km/h

12、');汽车最高车速=116km/h(2)其中：为第档对应不同转速（或车速）下的爬坡度求汽车最大爬坡度程序：五档车：n=600:10:4000;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;nmax=4000;G=m*g;ig=3.55 2.77 1.65 1.00 0.83;nT=0.87;r=0.37;f=0.0096;CDA=2.60;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.8;Iw2=3.6;Ft1=Tq*

13、ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G);disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度= 23.0844%四档车： n=600:10:4000;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=

14、600;nmax=4000;G=m*g;ig=6.12 3.11 1.69 1.00 ;nT=0.87;r=0.37;f=0.0096;CDA=2.60;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.8;Iw2=3.6;Ft1=Tq*ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G);disp('汽车最大爬坡度=');disp(imax);disp(

15、'%');汽车最大爬坡度= 43.0871%(3)求汽车行驶加速度倒数曲线程序：五档车：n=600:10:4000;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;nmax=4000;G=m*g;ig=3.55 2.771 1.65 1.00 0.83 ;nT=0.87;r=0.37;f=0.0096;CDA=2.60;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.8;Iw2=3.6;Ft1=Tq*ig(1)*i

16、0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fw1=CDA*ua1.2/21.15;Fw2=CDA*ua2.2/21.15;Fw3=CDA*ua3.2/21.15;Fw4=CDA*ua4.2/21.15;Fw5=CDA*ua5.2/

17、21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r2)+(If*ig(1)2*i02*nT)/(m*r2);deta2=1+(Iw1+Iw2)/(m*r2)+(If*ig(2)2*i02*nT)/(m*r2);deta3=1+(Iw1+Iw2)/(m*r2)+(If*ig(3)2*i02*nT)/(m*r2);deta4=1+(Iw1+Iw2)/(m*r2)+(If*ig(4)2*i02*nT)/(m*r2);deta5=1+(Iw1+Iw2)/(m*r2)+(If*ig(5)2*i02*nT)/(m*r2);a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1.

18、/a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis(0 99 0 10);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/

19、a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5');a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);四档车：n=600:10:4000;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;nmax=4000;G=m*g;ig=6.12 3.11

20、1.69 1.00 nT=0.87;r=0.37;f=0.0096;CDA=2.60;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.8;Iw2=3.6;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;Fw1=CDA*ua1.2/21.15;Fw2=CDA*ua

21、2.2/21.15;Fw3=CDA*ua3.2/21.15;Fw4=CDA*ua4.2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r2)+(If*ig(1)2*i02*nT)/(m*r2);deta2=1+(Iw1+Iw2)/(m*r2)+(If*ig(2)2*i02*nT)/(m*r2);deta3=1+(Iw1+Iw2)/(m*r2)+(If*ig(3)2*i02*nT)/(m*r2);deta4=1+(Iw1+Iw2)/(m*r2)+(If*ig(4)2*i02*nT)/(m*r2);a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=

22、(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4);axis(0 99 0 10);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3'

23、);gtext('1/a4');a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);4 图解法或编程绘制汽车动力因数特性曲线 Ft=Ff+Fi+Fw+Fj令为汽车的动力因数并以符号D表示（1）动力因素的计算同理，在发动机外特性曲线上选取8个点（与上面的相同），计算对应的驱动力与空气阻力，然后根据公式5算出各档对应的动力因素点。数据见表2-3所示。（2）滚动阻力系数的确定f=f0(1+ua2/19440)取f0=0.014据公式6算出对应车速下的滚动阻力系数，数据见表2.2所示。（3） 基本数据表四-

24、1四-2四-3四-4五-1五-2五-3五-4五-5a-19.3-19.3-19.3-19.3-19.3-19.3-19.3-19.3-19.3a1296296296296296296296296296a2-165-165-165-165-165-165-165-165-165a340.940.940.940.940.940.940.940.940.9a4-3.85-3.85-3.85-3.85-3.85-3.85-3.85-3.85-3.85c100010001000100010001000100010001000i05.835.835.835.835.835.835.835.835.83ig

25、6.123.111.6913.552.771.6510.830.8940.8940.8940.8940.8940.8940.8940.8940.894r0.370.370.370.370.370.370.370.370.37Cd*A2.62.62.62.62.62.62.62.62.6f0.00960.00960.00960.00960.00960.00960.00960.00960.0096m369536953695369536953695369536953695If0.220.220.220.220.220.220.220.220.225.45.45.45.45.45.45.45.45.4

26、v18406580154060801001.5061.1381.0481.0231.1771.1121.0461.0231.019n4604.14549.34591.13346.12225.54630.894137.723346.613468.99Ttq108.11114.12109.35174.49178.6104.29149.46173.34171.75Ff354.72354.72354.72354.72354.72354.72354.72354.72354.72Fw39.829196.69519.38786.727.659196.690442.553786.76121229.3144Ft

27、9320.14999.42603.22457.98931.24069.363473.862441.7602008.0715D0.25110.1290.05630.04520.240.1040.0820.0440.02105 图解或编程绘制II档起步，加速至70km/h的车速时间曲线以及距离时间曲线。编程如下：五档车编程如下：>>clearnT=0.894;r=0.37;f=0.0096;CDA=2.60;i0=5.83;If=0.22;Iw1=1.8;Iw2=3.4;L=3.4;a=2.0;hg=0.9;m=3695;g=10;G=m*g; ig=3.55 2.77 1.65 1.

28、00 0.83;nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5 deta(i)=1+(Iw1+Iw2)/(m*r2)+(If*(ig(i)2*i02*nT)/(m*r2);endua=6:0.01:99;N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.2/21.15;for i=1:N k=i; if ua(i)<=u2(2) n=ua(i)*(ig(2)*i0/r)/0.377; T

29、q=-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/1000)4; Ft=Tq*ig(2)*i0*nT/r; inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i); delta(i)=0.01*inv_a(i)/3.6; elseif ua(i)<=u2(3) n=ua(i)*(ig(3)*i0/r)/0.377; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/1000)4; Ft=Tq*ig

30、(3)*i0*nT/r; inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i); delta(i)=0.01*inv_a(i)/3.6; elseif ua(i)<=u2(4) n=ua(i)*(ig(4)*i0/r)/0.377; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/1000)4; Ft=Tq*ig(4)*i0*nT/r; inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i); delta(i)=0.01*inv_a(i)/3.6; else n=ua(

31、i)*(ig(5)*i0/r)/0.377; Tq=-19.3+296*(n/1000)-165*(n/1000)2+40.9*(n/1000)3-3.85*(n/1000)4; Ft=Tq*ig(5)*i0*nT/r; inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i); delta(i)=0.01*inv_a(i)/3.6; end a=delta(1:k); t(i)=sum(a);endplot(t,ua);axis(0 80 0 100);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel

32、('速度ua（km/h）');>> ginputans =2 汽车燃料经济性计算（等速、多工况）例如，已知汽车结构参数同题1，负 荷特性拟合公式为N/(rmin-1)B0B1B2B3B48201427-41672-6.00.1812001355-30437-2.050.04316001284-19015-0.510.00720001123-121.67.0-0.1850.001926001141-994.5-0.0910.000730001051-742.86-0.0510.0003534001233-853.0-0.04750.000338801130-450.7

33、1-0.0008-0.00004   怠速油耗为0.3ml/s（400r/min），求1) 汽车功率平衡图；2) 最高档与次高档的等速百公里油耗曲线。3) 按6工况循环行驶的百公里油耗。计算中确定燃油消耗率值b，若发动机转速与负荷特性中给定的转速不相等，可用相邻转速的两根曲线用插值法求得。解：1.五档车汽车功率平衡图MATLAP程序图：n=600:10:4000;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;nmax=4000;G=m*g;i

34、g=3.55 2.77 1.65 1.00 0.83nT=0.87;r=0.37;f=0.0096;CDA=2.60;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.8;Iw2=3.6;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Pe1=Tq.*ig(1)*i0.*ua1./(3600*r);Pe2=Tq.*ig(2)*i0.*ua2./(3600*r);Pe3=Tq.*ig(

35、3)*i0.*ua3./(3600*r);Pe4=Tq.*ig(4)*i0.*ua4./(3600*r);Pe5=Tq.*ig(5)*i0.*ua5./(3600*r);ua=0:0.35:119;Ff=G*f;Fw=CDA*ua.2/21.15;Pf=Ff*ua/3600;Pw=Fw.*ua/3600;Pe0=(Pf+Pw)./nT;Pe=max(Pe1);plot(ua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua5,Pe5,ua,Pe0,ua,Pe);axis(0 119 0 100);title('汽车功率平衡图');xlabel('ua(km

36、/h)');ylabel('Pe(kw)');gtext('1'),gtext('2'),gtext('3'),gtext('4'),gtext('5'),gtext('(Pf+Pw)/et'),gtext('Pe');四档车汽车功率平衡图程序：clearn=600:10:4000;Tq=-19.3+296*(n/1000)-165*(n/1000).2+40.9*(n/1000).3-3.85*(n/1000).4;m=3695;g=10;nmin=600;

37、nmax=4000;G=m*g; ig=6.12 3.11 1.69 1.00 nT=0.87;r=0.37;f=0.0096;CDA=2.60;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.8;Iw2=3.6;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;Pe1=Tq.*ig(1)*i0.*ua1./(3600*r);Pe2=Tq.*ig(2)*i0.*ua2./(3600*r);Pe3=Tq.*ig(3)*i0.*ua

38、3./(3600*r);Pe4=Tq.*ig(4)*i0.*ua4./(3600*r);ua=0:0.35:119;Ff=G*f;Fw=CDA*ua.2/21.15;Pf=Ff*ua/3600;Pw=Fw.*ua/3600;Pe0=(Pf+Pw)./nT;Pe=max(Pe1);plot(ua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua,Pe0,ua,Pe);axis(0 119 0 100);title('汽车功率平衡图');xlabel('ua(km/h)');ylabel('Pe(kw)');gtext('1&#

39、39;),gtext('2'),gtext('3'),gtext('4'),gtext('(Pf+Pw)/et'),gtext('Pe');（2）计算燃油消耗率按下式计算行驶阻尼功率 ： （kw） 按下式计算对应的发动机转速： （r/min） 其中：直接档（或超速档）传动比根据万有特性，通过插值求出对应于上面计算结果、的有效燃油消耗率（g/kw·h） 计算直接档等速百公里油耗 （L/100km） 其中： 为燃料的重度最高档与次高档的等速百公里油耗曲线程序：五档车程序设计 clear n=600:1:400

40、0;m=3695;g=10;G=m*g;ig=3.55 2.77 1.65 1.00 0.83;nT=0.894;r=0.37;f=0.0096;CDA=2.6;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.22;Iw1=1.8;Iw2=3.6;n0=820 1200 1600 2000 2600 3000 3400 3880;B00=1427 1355 1284 1123 1141 1051 1233 1130;B10=-416 -304 -190 -121.6 -99 -74 -85 -45;B20=72 37 15 7 4.5 2.86 3.0 0.71;B30=-6.

41、0 -2.05 -0.51 -0.185 -0.091 -0.051 -0.0475 -0.0008;B40=0.18 0.043 0.007 0.0019 0.0007 0.00035 0.0003 -0.00004;B0=spline(n0,B00,n);B1=spline(n0,B10,n);B2=spline(n0,B20,n); B3=spline(n0,B30,n);B4=spline(n0,B40,n);Ff=G*f;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fz4=Ff+CDA*(ua4.2)/21.15;Fz5=Ff+CDA*

42、(ua5.2)/21.15;Pe4=Fz4.*ua4./(nT*3.6*1000);Pe5=Fz5.*ua5./(nT*3.6*1000);for i=1:1:3401 b4(i)=B0(i)+B1(i)*Pe4(i)+B2(i)*Pe4(i).2+B3(i)*Pe4(i).3+B4(i)*Pe4(i).4;b5(i)=B0(i)+B1(i)*Pe5(i)+B2(i)*Pe5(i).2+B3(i)*Pe5(i).3+B4(i)*Pe5(i).4;endpg=7.0;Q4=Pe4.*b4./(1.02.*ua4.*pg); Q5=Pe5.*b5./(1.02.*ua5.*pg);plot(ua4

43、,Q4,ua5,Q5);axis(0 100 10 30);title('最高档与次高档等速百公里油耗曲线');xlabel('ua(km/h)');ylabel('百公里油耗（L/100km）');gtext('4'),gtext('5'); 四档车程序设计如下 clear n=600:1:4000;m=3695;g=10;G=m*g;ig=6.12 3.11 1.69 1.00 ;nT=0.894;r=0.37;f=0.0096;CDA=2.6;i0=5.83;L=3.4;a=2.0;hg=0.9;If=0.2

44、2;Iw1=1.8;Iw2=3.6;n0=820 1200 1600 2000 2600 3000 3400 3880;B00=1427 1355 1284 1123 1141 1051 1233 1130;B10=-416 -304 -190 -121.6 -99 -74 -85 -45;B20=72 37 15 7 4.5 2.86 3.0 0.71;B30=-6.0 -2.05 -0.51 -0.185 -0.091 -0.051 -0.0475 -0.0008;B40=0.18 0.043 0.007 0.0019 0.0007 0.00035 0.0003 -0.00004;B0=s

45、pline(n0,B00,n);B1=spline(n0,B10,n);B2=spline(n0,B20,n); B3=spline(n0,B30,n);B4=spline(n0,B40,n);Ff=G*f;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;Fz4=Ff+CDA*(ua4.2)/21.15;Fz3=Ff+CDA*(ua3.2)/21.15;Pe4=Fz4.*ua4./(nT*3.6*1000);Pe3=Fz3.*ua3./(nT*3.6*1000);for i=1:1:3401 b4(i)=B0(i)+B1(i)*Pe4(i)+B2(i

46、)*Pe4(i).2+B3(i)*Pe4(i).3+B4(i)*Pe4(i).4;b3(i)=B0(i)+B1(i)*Pe3(i)+B2(i)*Pe3(i).2+B3(i)*Pe3(i).3+B4(i)*Pe3(i).4;endpg=7.0;Q3=Pe3.*b3./(1.02.*ua3.*pg); Q4=Pe4.*b4./(1.02.*ua4.*pg);plot(ua3,Q3,ua4,Q4);axis(0 100 10 30);title('最高档与次高档等速百公里油耗曲线');xlabel('ua(km/h)');ylabel('百公里油耗（L/100km）');gtext('3'),gtext('4'); 怠速工况的计算燃油消耗量： （ml） 其中：为怠速燃油消耗率为怠速工况时间（由标准工况决定）行驶距离 （m） 匀速工况的计算 燃油消耗量发动机功率发动机转速 其中：车轮滚动半径由、根据万有特性插值得到燃油消耗率按下式计算燃油消耗量其中：t为匀速工况时间行驶距离 其中