汽车理论课程设计报告.docx

1、汽车理论课程设计报告汽车理论设计报告汽车动力性、经济性的计算机模拟 学号20152805037姓名张少波学院交通学院班级车辆本1501班2017 年 12 月 18 日1汽车驱动力图 11.1汽车驱动力图简介 11.2汽车驱动力图 22汽车驱动力-行驶阻力平衡图 22.1汽车驱动力-行驶阻力平衡图简介 22.1.1滚动阻力 32.1.2空气阻力 32.2汽车驱动力-行驶阻力平衡图 33汽车爬坡度图 33.1汽车爬坡度图简介 33.2汽车爬坡度图 44汽车加速度图 54.1汽车加速度图简介 54.2汽车加速度图 55汽车动力特性图 65.1汽车动力特性图简介 65.2汽车动力特性图 66汽车功率

2、平衡图 66.1汽车功率平衡图简介 66.2汽车功率平衡图 77汽车百公里油耗图 77.1汽车百公里油耗图简介 77.2汽车百公里油耗图 8 8 8汽车动力性、经济性的计算机模拟张少波（交通学院，车辆工程，1501班，20152805037）摘要：通过MATLAB计算机软件进行汽车动力性、经济性的计算机模拟，模拟得出汽车的各种特性曲线。包括：驱动力图，驱动力-行驶阻力平衡图，爬坡度图，加速度图，动力特性图，功率平衡图，百公里油耗图。同时对汽车特性曲线的计算公式加以统计汇总学习。关键词：驱动力图，驱动力-行驶阻力平衡图，爬坡度图，加速度图，动力特性图，功率平衡图，百公里油耗图，汽车特性曲线的计算

3、公式1 汽车驱动力图1.1 汽车驱动力图简介驱动力是由发动机的转矩经传动系统传至驱动轮上得到的。汽车发动机产生的转矩，经传动系传至驱动轮上。此时作用于驱动轮上的转矩Tt产生一对地面的圆周力F0,地面对驱动轮的反作用力Ft（方向与F0相反）既是驱动汽车的外力，此外力称为汽车的驱动力。其数值为 （1.1）式中，Tt为作用于驱动轮上的转矩；r为车轮半径。作用于驱动轮上的转矩Tt是由发动机产生的转矩经传动系传至车轮上的。若令Ttq表示发动机转矩，ig表示变速器的传动比，i0表示主减速器的传动比，表示传动系的机械效率，则有 （1.2）对于装有分动器、轮边减速器、液力传动等装置的汽车，上式应计入相应的传动

4、比和机械效率。因此驱动力为 （1.3）汽油发动机使用外特性的Tq-n曲线的拟合公式为 （1.4）一般用根据发动机外特性确定的驱动力与车速之间的函数关系曲线Ft-ua来全面表示汽车的驱动力，称为汽车的驱动力图。设计中的汽车有了发动机的外特性曲线、传动系的传动比、传动效率、车轮半径等参数后，即可用式（1.1）求出各个档位的Ft值，再根据发动机转速与汽车行驶速度之间转换关系求出ua，即可求得各个档位的Ft于ua曲线。发动机转速与汽车行驶速度之间的关系式为 （1.5）式中，ua为汽车行驶速度（km/h）;n为发动机转速（r/min）;r为车轮半径（m）;ig为变速器传动比；i0为主减速器传动比。另：实

5、际行驶中，发动机常在节气门部分开启下工作，相应的驱动力要比它小。1.2汽车驱动力图其车辆基本参数见附录1，其Matlab程序见附录2。附图：2 汽车驱动力-行驶阻力平衡图2.1 汽车驱动力-行驶阻力平衡图简介汽车在水平道路上等速行驶时，必须克服来自地面的滚动阻力和来自空气的空气阻力。滚动阻力以符号Ff表示，空气阻力以符号Fw表示。当汽车在坡道上上坡行驶时，还必须克服重力沿坡道的分力，称为坡道阻力，以符号Fi表示。汽车加速行驶时还需要克服加速阻力，以符号Fj表示。因此，汽车行驶的总阻力为 （2.1） 在汽车驱动力图上把汽车行驶中经常遇到的滚动阻力和空气阻力一并画出，做出汽车驱动力-行驶阻力平衡图

6、。2.1.1 滚动阻力滚动阻力等于滚动阻力系数与车轮负载之乘积。即 （2.2）其中滚动阻力系数由试验确定。滚动阻力系数与路面的种类、行驶车速以及轮胎的构造、材料、气压等有关。这里选取滚动阻力系数为良好的沥青或混凝土路面滚动阻力系数。（详见附录1）。2.1.2 空气阻力汽车直线行驶时受到的空气作用力在行驶方向上的分力称为空气阻力。此处讨论无风条件下汽车的运动，空气阻力计算公式为 （2.3）式中，A为迎风面积，即汽车行驶方向的投影面积（m2）；CD为空气阻力系数；ua2为车速，km/h。（详见附录1）。2.2 汽车驱动力-行驶阻力平衡图其车辆基本参数见附录1，其Matlab程序见附录3。3 汽车爬

7、坡度图3.1 汽车爬坡度图简介汽车的上坡能力是用满载（或某一载质量）时汽车在良好路面上的最大爬坡度imax表示的。另外imax代表了汽车的极限爬坡能力，它应比实际行驶中遇到的道路最大坡度超出很多，这是因为应考虑到在实际坡道行驶时，在坡道上停车后顺利起步加速、克服松软坡道路面的大阻力、克服坡道上崎岖不平路面的局部大阻力等要求的缘故。也可以理解汽车的爬坡能力，指汽车在良好路面上克服Ff+Fw后的余力全部用来（即等速）克服坡度阻力时能爬上的坡度。因此 Fi=Ft-（Ff+Fw） （3.1）一般汽车最大爬坡度达30%左右，因此利用汽车行驶方程式确定1档及低档爬坡能力时，应采用Gsin作为坡道阻力，即上

8、式应为 （3.2）即 （3.3）根据驱动力-行驶阻力图求出汽车能爬上之坡道角，相应地根据tan=i可求出坡度值。3.2 汽车爬坡度图其车辆基本参数见附录1，其Matlab程序见附录4。4 汽车加速度图4.1 汽车加速度图简介汽车的加速能力可用它在水平良好路面上行驶时能产生的加速度来评价。而汽车加速行驶时，需克服的其质量加速运动时的惯性力，就是加速阻力Fj。汽车的质量分为平移质量和旋转质量两部分。为了简易计算旋转质量，常用系数作为计算计入旋转质量惯性力偶矩后的汽车旋转质量换算系数。主要与飞轮的转动惯量、车轮的转动惯量以及传动系的传动比有关。其公式如下： （4.1） 式中，Iw为车轮的转动惯量（k

9、gm2）；If为飞轮的转动惯量（kgm2）；i0为主传动比；ig为变速器的速比。汽车加速时其驱动力需全力加速，即坡道阻力为0 。 4.2 汽车加速度图其车辆基本参数见附录1，其Matlab程序见附录5。5 汽车动力特性图5.1 汽车动力特性图简介用动力特性图来分析汽车动力性。将汽车行驶方程两边除以汽车重力并整理如下： （5.1） (5.2)令 为汽车的动力因数并以符号D表示，则 （5.3） 汽车在各档下的动力因数与车速的关系曲线称为动力特性图。5.2 汽车动力特性图其车辆基本参数见附录1，其Matlab程序见附录6。6 汽车功率平衡图6.1 汽车功率平衡图简介在汽车行驶的每一瞬间，发动机发出的

10、功率始终等于机械传动损失功率与全部运动阻力多消耗的功率。汽车运动阻力所消耗的功率有滚动阻力功率Pf、空气阻力功率Pw、坡度阻力功率Pi及加速阻力功率Pj。将汽车行驶方程式两边乘以行驶车速ua，并经单位换算整理出汽车功率平衡方程式（式中功率单位为kW），如下 （6.1）将汽车经常遇到的阻力功率 （Pf+Pw）。发动机功率对车速的曲线绘制。得到汽车功率平衡图。6.2 汽车功率平衡图其车辆基本参数见附录1，其Matlab程序见附录7。7 汽车百公里油耗图7.1 汽车百公里油耗图简介根据参考文献1，等速行驶工况燃油消耗量的计算。在万有特性图上有等燃油消耗率曲线。根据曲线可以确定发动机在一定转速n,发出

11、一定功率Pe时的燃油消耗率b。 （7.1） （7.2） （7.3）7.2 汽车百公里油耗图其车辆基本参数见附录1，其Matlab程序见附录8。1 余志生.汽车理论M.北京：机械工业出版社，2009.3.附录1车辆基本参数为发动机最低转速 600r/min ， 最高转速 4000r/min ， 装载质量 2000kg ，整车整备质量 1800kg ,车轮半径 0.367m ，传动系效率 0.85 ,滚动阻力系数 0.013 ，空气阻力系数与迎风面积之积为 2.77m2 ， 主减速器传动比5.83 ,飞轮转动惯量 0.218kg m2 ,前轮总转动惯量为 1.798kg m2 ,后轮总的转动惯量为

12、 3.598kg m2,轴距为 3.2m 质心至前轴距离 1.947m ，质心高 0.9m ， 四挡变速器的变速比分别为 6.09 ， 3.09 ， 1.71 ， 1 。附录2close all;n=linspace(600,4000,100);%设定转速范围(linspace(x1,x2,N)在x1x2范围内选取N个点，均分。)ua1=0.377*0.367*n/(5.83*6.09);%公式计算各档车速范围ua2=0.377*0.367*n/(5.83*3.09);ua3=0.377*0.367*n/(5.83*1.71);ua4=0.377*0.367*n/(5.83*1);Ttq=(-

13、19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4);%?Ttq-n曲线拟合公式x1=6.09;x2=3.09;x3=1.71;x4=1;Ft1=x1*Ttq*5.83*0.85/0.367;%公式计算各档对应转速下的驱动力Ft2=x2*Ttq*5.83*0.85/0.367;Ft3=x3*Ttq*5.83*0.85/0.367; Ft4=x4*Ttq*5.83*0.85/0.367;plot(ua1,Ft1,g,ua2,Ft2,r,ua3,Ft3,y,ua4,Ft4,b);%汽车驱动力-行驶阻

14、力平衡图xlabel(ua/(km/h);%汽车驱动力-行驶阻力平衡图ylabel(F/N);legend(1,2,3,4);grid on;title(汽车驱动力图);附录3close all;n=linspace(600,4000);%设定转速范围ua1=0.377*0.367*n/(5.83*6.09);%公式计算各档车速范围ua2=0.377*0.367*n/(5.83*3.09);ua3=0.377*0.367*n/(5.83*1.71);ua4=0.377*0.367*n/(5.83*1);Ttq=(-19.13+259.27*(n/1000)-165.44*(n/1000).2+

15、40.874*(n/1000).3-3.8445*(n/1000).4);%?Ttq-n曲线拟合公式x1=6.09;x2=3.09;x3=1.71;x4=1;Ft1=x1*Ttq*5.83*0.85/0.367;%公式计算各档对应转速下的驱动力Ft2=x2*Ttq*5.83*0.85/0.367;Ft3=x3*Ttq*5.83*0.85/0.367;Ft4=x4*Ttq*5.83*0.85/0.367;F1=3800*9.8*0.013+(2.77/21.15)*ua1.2;%公式计算各档对应的各个车速下的行驶阻力F2=3800*9.8*0.013+(2.77/21.15)*ua2.2;F3=

16、3800*9.8*0.013+(2.77/21.15)*ua3.2;F4=3800*9.8*0.013+(2.77/21.15)*ua4.2;plot(ua1,Ft1,g,ua2,Ft2,b,ua3,Ft3,y,ua4,Ft4,m,ua1,F1,r,ua2,F2,r,ua3,F3,r,ua4,F4,r);%汽车驱动力-行驶阻力平衡图?xlabel(ua/(km/h);%汽车驱动力-行驶阻力平衡图ylabel(F/N);legend(1,2,3,4,行驶阻力);grid on;title(汽车驱动力-行驶阻力平衡图);附录4close all;n=600:10:4000;Tq=-19.13+25

17、9.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;m=3800;g=9.8;nmin=600;nmax=4000;G=m*g;ig=6.09 3.09 1.71 1 ;nT=0.85;%传动系效率r=0.367;%车轮半径f=0.013;%滚动阻力系数CDA=2.77;%空气阻力系数与迎风面积之积i0=5.83;%主减速器传动比L=3.2;%轴距a=1.947;%质心至前轴距离hg=0.9;%质心高If=0.218;%飞轮转动惯量Iw1=1.798;%前轮总转动惯量Iw2=3.598;%后轮总转动惯量Ft1=

18、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*ua2.2/21.15;Fw3=CDA*ua3.2/21.15;Fw4=CDA*ua4.2/21.15;Ff=G*f;i1=asin(Ft1-Ff-Fw1)/G);i2=asin(Ft2-Ff-Fw2

19、)/G);i3=asin(Ft3-Ff-Fw3)/G);i4=asin(Ft4-Ff-Fw4)/G);plot(ua1,i1,ua2,i2,ua3,i3,ua4,i4);title(汽车的爬坡度图);xlabel(ua/(km*h-1);ylabel(i/%);legend(1,2,3,4);grid on;附录5close all;n=600:10:4000;Tq=-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;

20、ig=6.09 3.09 1.71 1.00;nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;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*ua

21、1.2/21.15;Fw2=CDA*ua2.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

22、-Ff-Fw1)/(deta1*m);a2=(Ft2-Ff-Fw2)/(deta2*m);a3=(Ft3-Ff-Fw3)/(deta3*m);a4=(Ft4-Ff-Fw4)/(deta4*m);plot(ua1,a1,ua2,a2,ua3,a3,ua4,a4);title(汽车的行驶加速度曲线);xlabel(ua/(km*h-1);ylabel(a/(m*s-2);legend(1,2,3,4);grid on;附录6close all;n=linspace(600,4000);%设定转速范围 ua1=0.377*0.367*n/(5.83*6.09);%公式计算各档车速范围 ua2=0.3

23、77*0.367*n/(5.83*3.09); ua3=0.377*0.367*n/(5.83*1.71); ua4=0.377*0.367*n/(5.83*1.00); Ttq=(-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4);%Ttq-n 曲线拟合公式 m=3800;g=9.8;G=m*g;x1=6.69;x2=3.09; x3=1.71;x4=1; Ft1=x1*Ttq*5.83*0.85/0.367;%公式计算各档对应转速下的驱动力Ft2=x2*Ttq*5.83*0.85/0.3

24、67; Ft3=x3*Ttq*5.83*0.85/0.367;Ft4=x4*Ttq*5.83*0.85/0.367; F1=(2.77/21.15)*ua1.2;%公式计算各档对应的各个车速下的空气阻力F2=(2.77/21.15)*ua2.2;F3=(2.77/21.15)*ua3.2;F4=(2.77/21.15)*ua4.2;D1=(Ft1-F1)/G;D2=(Ft2-F2)/G;D3=(Ft2-F2)/G;D4=(Ft3-F3)/G;plot(ua1,D1,ua2,D2,ua3,D3,ua4,D4);title(汽车的动力特性图);xlabel(ua/(km*h-1);ylabel(D

25、);legend(1,2,3,4)grid on;附录7close all;n=linspace(600,4000);%设定转速范围?ua1=0.377*0.367*n/(5.83*6.09);%公式计算各档车速范围?ua2=0.377*0.367*n/(5.83*3.09);ua3=0.377*0.367*n/(5.83*1.71);ua4=0.377*0.367*n/(5.83*1.00);Ttq=(-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4);%?Ttq-n?曲线拟合公式?x1=

26、6.09;x2=3.09;x3=1.71;x4=1;Ft1=x1*Ttq*5.83*0.85/0.367;%公式计算各档对应转速下的驱动力?Ft2=x2*Ttq*5.83*0.85/0.367;Ft3=x3*Ttq*5.83*0.85/0.367;Ft4=x4*Ttq*5.83*0.85/0.367;F1=3800*9.8*0.013+(2.77/21.15)*ua1.2;%公式计算各档对应的各个车速下的行驶阻力?F2=3800*9.8*0.013+(2.77/21.15)*ua2.2;F3=3800*9.8*0.013+(2.77/21.15)*ua3.2;F4=3800*9.8*0.013

27、+(2.77/21.15)*ua4.2;Pe1=Ft1.*ua1./(3600*0.85);%计算各档对应转速下的功率?Pe2=Ft2.*ua2./(3600*0.85);Pe3=Ft3.*ua3./(3600*0.85);Pe4=Ft4.*ua4./(3600*0.85);P1=F1.*ua1./(3600*0.85);%计算各档对应的各个车速下的行驶阻力功率?P2=F2.*ua2./(3600*0.85);P3=F3.*ua3./(3600*0.85);P4=F4.*ua4./(3600*0.85);plot(ua1,Pe1,b,ua2,Pe2,y,ua3,Pe3,m,ua4,Pe4,y,

28、ua1,P1,r,ua2,P2,r,ua3,P3,r,ua4,P4,r);%汽车功率平衡图?xlabel(ua/(km/h);%汽车功率平衡图ylabel(P/kW);grid on;legend(1,2,3,4);title(汽车功率平衡图);附录8n=600:1:4000;r=0.367;i0=5.83;eff=0.85;f=0.013;CdA=2.77;m=3880;g=9.8;G=m*g;Ttq=-19.13+259.27*n/1000-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;n1=815,1207,1614,2012

29、,2603,3006,3403,3494;b0=1300,1354.56,1280.4,1125.9,1146.0,1053.2,1240.9,1130.7;b1=-418.46,-300.98,-190.75,-120,-98.5,-74.8,-85.56,-45.23;b2=72.4,36.8,15.21,7.0,4.56,2.92,3,0.765;b3=-6,-2.042,-0.6114,-0.2,-0.1,-0.052,-0.05,-0.0008;b4=0.2,0.043,0.0068,0.0019,0.0007,0.000352,0.0003,-0.00004;B0=spline(n1,b0,n); B1=spline(n1,b1,n); B2=spline(n1,b2,n); B3=spline(n1,b3,n); B4=spline(n1,b4,n); deng=8; figure;for ig=6.09,3.09,1.71,1.00 Ua=0.377*r*n/ig/i0; Pf=G*f*Ua/3600; Pw=CdA*Ua.3/76140; P=(Pf+Pw)/eff; b=B0+B1.*P+B2.*P.2+B3.*P.3+B4.*P.4; Qs=P.*b./(1.02*Ua*deng); plot(Ua,Qs); hold