汽车理论1.3和2.7matlab编程答案.doc

《汽车理论1.3和2.7matlab编程答案.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、孙野200812681。3(1)绘制汽车驱动力与行驶阻力平衡图选用5挡变速器进行整车性能计算发动机转速与汽车行驶速度之间的关系:发动机使用外特性的Tq—n曲线拟合公式:汽车驱动力：行驶阻力：用matlab编写程序如下:m1=2000；m2=1800；mz=3880；g=9。81;r=0。367;CdA=2.77;f=0.013;nT=0.85;ig=［5。562.7691。6441。000.793］;i0=5.83；If=0。218；Iw1=1。798；Iw2=3。598;％输入已知参数Iw=2*Iw1+4*Iw2;fori=1：69n(i）=(

2、i+11）＊50；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%从600～4000rpm油拟合公式计算发动机转距forj=1：5fori=1：69Ft(i，j)=Ttq（i)*ig（j)＊i0＊nT/r;ua(i,j)=0.377*r＊n（i)/(ig（j）*i0）；Fz（i，j）=CdA*ua(i,j)^2/21.15+mz*g＊f;endend％计算各档位下的驱动力和行驶阻力plot（ua

3、,Ft，ua,Fz）title（'汽车驱动力与行驶阻力平衡图')；xlabel(’ua（km/h）’）；ylabel('Ft（N）'）;%作汽车驱动力与行驶阻力平衡图本文为互联网收集，请勿用作商业用途得到汽车驱动力与行驶阻力平衡图如下：（2)①计算汽车最高车速Matlab程序如下：fork=1：175n1（k)=3300+k*0.1；Ttq(k）=-19.313+295。27＊(n1（k）/1000)—165.44＊(n1(k)/1000)^2+40.874*（n1（k)/1000）^33。8445＊（n1(k)/1000）^4；Ft(k）=Tt

4、q(k）*ig(5）＊i0＊nT/r;ua(k)=0.377＊r＊n1（k)/(ig（5）＊i0);Fz（k）=CdA*ua(k)^2/21.15+mz*g＊f;E（k）=abs(（Ft(k）—Fz（k)));end％从3300rpm开始以0。1rpm为步长搜索直到找到fork=1:175％驱动力与行驶阻力之差绝对值最小点(Ⅴ档处)if(E(k）==min（E）)％通过其他检验程序可估计出该点大致位置（3310rpm附近)disp（'汽车最高车速='）;disp(ua(k））;disp（'km/h’)；end％输出结果end汽车最高车速=98.7

5、627km/h计算得汽车最高车速为98.7627km/h。②计算最大爬坡度Matlab程序如下:forp=1:150%从2000rpm开始以0。5为步长搜索找驱动力n2（p）=2000+p＊0。5;％与行驶阻力之差最大点（Ⅰ档处)Ttq(p)=-19。313+295。27＊(n2(p）/1000)-165.44＊（n2(p)/1000）^2+40.874*(n2（p)/1000）^3—3.8445＊(n2(p)/1000）^4;Ft（p）=Ttq(p)*ig（1)*i0*nT/r;ua(p)=0.377*r＊n2(p）/（ig（1）*i0）;Fz

6、（p）=CdA*ua(p）^2/21。15+mz*g*f;af(p）=asin((Ft（p）—Fz(p））/（mz*g））；end％通过其他检验程序可估计出该点大致位置（2000rpm附近）forp=1:150if(af（p）==max(af))i=tan（af(p）)；disp(’汽车最大爬坡度=')；disp(i)；endend汽车最大爬坡度=0.3518计算得最大爬坡度为0。3518（3）①绘制汽车行驶加速度倒数曲线Matlab程序如下：fori=1:69%开始计算绘制加速度倒数曲线所需数据n（i）=(i+11）*50;Ttq(i）=-19

7、。313+295。27＊(n（i)/1000）-165.44*（n(i)/1000)^2+40.874*（n（i）/1000）^3—3。8445*(n（i)/1000)^4;endforj=1：5fori=1:69deta=1+Iw/（mz*r^2）+If＊ig(j)^2＊i0^2＊nT/(mz*r^2)；ua（i，j）=0。377＊r＊n(i）/(ig（j）*i0)；a（i，j）=(Ttq（i)＊ig（j）＊i0*nT/r—CdA＊ua（i,j）^2/21。15-mz*g＊f）/（deta*mz);if(a(i，j)<=0a(i，j)=a（i-

8、1,j）；endif(a（i,j）>0。05）％剔除加速度值小于0.05的数据b1（i，j)=a（i,j）；u1（i，j）=ua（i,j