汽车动力性计算matlab程序
汽车理论课程设计：基于Matlab的汽车动力性的仿真
2009 届海南大学机电工程学院汽车工程系汽车理论课程设计题目：汽车动力性的仿真学院：机电工程学院专业：09级交通运输姓名：黄生锐学号：20090504指导教师：编号名称件数页数编号名称件数页数1 课程设计论文 1 3Matlab编程源程序 12 设计任务书 12012年6月20日成绩汽车理论课程设计任务书姓名黄生锐学号20090504 专业09交通运输课程设计题目汽车动力性的仿真内容摘要：本设计的任务是对一台Passat 1.8T手动标准型汽车的动力性能进行仿真。
采用MATLAB编程仿真其性能，其优点是：一是能过降低实际成本，提高效率；二是获得较好的参数模拟，对汽车动力性能提供理论依据。
主要任务：根据该车的外形、轮距、轴距、最小离地间隙、最小转弯半径、车辆重量、满载重量以及最高车速等参数，结合自己选择的适合于该车的发动机型号求出发动机的最大功率、最大扭矩、排量等重要的参数。
并结合整车的基本参数，选择适当的主减速比。
依据GB、所求参数，结合汽车设计、汽车理论、机械设计等相关知识，计算出变速器参数，进行设计。
论证设计的合理性。
设计要求：1、动力性分析：1)绘制汽车驱动力与行驶阻力平衡图；2)求汽车的最高车速、最大爬坡度；3)用图解法或编程绘制汽车动力特性曲线4)汽车加速时间曲线。
2、燃油经济性分析：1) 汽车功率平衡图；完成内容：1．Matlab编程汽车驱动力与行驶阻力平衡图2．编程绘制汽车动力特性曲线图3．编程汽车加速时间曲线图4．课程设计论文1份汽车动力性仿真摘要本文是对Passat 1.8T 手动标准型汽车的动力性能采用matlab 编制程序，对汽车动力性进行计算。
从而对汽车各个参数做出准确的仿真研究，为研究汽车动力性提供理论依据，本文主要进行的汽车动力性仿真有：最高车速、加速时间和最大爬坡度。
及相关汽车燃油性经济。
关键词：汽车；动力性；试验仿真；matlab1. Passat 1.8T 手动标准型汽车参数功率Pe （kw ）转速n （r/min ）15 1000 36 1750 50 2200 66 2850 80 3300 90 4000 110 5100 1055500各档传动比主减速器传动比第1档 3.665 4.778第2档 1.999 第3档 1.407 第4档 1 第5档 0.472 车轮半径0.316（m ）传动机械效率0.91 假设在良好沥青或水泥路面上行驶，滚动阻力系数 0.014 整车质量1522kgC D A2.4m22. 最高车速汽车的最高车速是指汽车标准满载状态，在水平良好的路面(混凝土或沥青路面)上所能达到的最高行驶速度。
汽车动力性matlab仿真源程序
汽车动力性matlab仿真源程序clcn=[1500:500:5500];%转速范围T=[78.59 83.04 85.01 86.63 87.09 85.87 84.67 82.50 80.54];%对应各转矩dt=polyfit(n,T,3);%对发动机输出转矩特性进行多项式拟合，阶数取4n1=1000:100:5500;%t=polyval(dt,n1);figure(1)title('发动机外特性')plot(n1,t,n,T,'o'),grid on%图示发动机输出转矩特性%汽车驱动力计算G=input('整车重力/N，G=');%输入970*9.8ig=[3.416 1.894 1.28 0.914 0.757];%变速器速比k=1:5;%5个前进档r=0.272;i0=4.388;eta=0.9;ngk=[800 800 800 800 800];ngm=[5500 5500 5500 5500 5500];ugk=0.377.*r.*ngk(k)./(ig(k).*i0);%计算每一档发动机800rpm 时的最低行驶速度ugm=0.377.*r.*ngm(k)./(ig(k).*i0);%计算每一档发动机5400rpm最高行驶速度for k=1:5%依次计算5个档的驱动力u=ugk(k):ugm(k);n=ig(k)*i0.*u./r/0.377;t=54.8179+2.2441.*(n./100)-4.8003.*(n./1000).^2+2.815e-10.*n.^3Ft=t.*ig(k).*i0*eta/r;figure(2)plot(u,Ft)hold on,grid on %保证K的每次循环的图形都保留显示end%行驶阻力计算f0=0.009;f1=0.002;f4=0.0003;%三者都是轿车滚动阻力系数% disp'空气阻力系数Cd=0.3--0.41,迎风面积A=1.7--2.1'Cd=input('空气阻力系数Cd=');%输入0.3A=input('迎风面积/m2,A=');%输入2.3u=0:10:180;f=f0+f1.*(u./100)+f4.*(u./100).^4;Ff=G*f;%计算滚动阻力Fw=Cd*A.*u.^2./21.15;%计算空气阻力F=Ff+Fw;%滚动阻力、空气阻力之和title('驱动力-阻力图（五档速比为3.416 1.894 1.28 0.914 0.757）')plot(u,F,'mo-');grid on%图解法求最高车速for u=50:180;k=5;n=ig(k)*i0.*u./r/0.377;t=54.8179+2.2441.*(n./100)-4.8003.*(n./1000).^2+2.815e-10.*n.^3;Ft=t.*ig(k).*i0*eta/r;f=f0+f1.*(u./100)+f4.*(u./100).^4;Ff=G*f;Fw=Cd*A.*u.^2./21.15;F=Ff+Fw;if abs(Ft-F)<1;umax=u;breakendenddisp('== == =汽车动力性能仿真计算结果== == =')disp('驱动力-行驶阻力平衡图及最高车速')fprintf('汽车理论最高车速（驱动力与行驶阻力曲线交点）Vmax=%3.3f km/h\n',umax)（注：文档可能无法思考全面，请浏览后下载，供参考。
汽车整体动力性MATLAB仿真计算
车各种性能中最基本、最重要的性能。
动力性评价指标
汽车动力性评价指标有最高车速、加速 时间、最大爬坡度等，与之对应的试验内容 有最高车速的测试、汽车起步连续换挡加速 时间与超车加速时间的测试和汽车最大爬坡 度的测试。另外，按照我国标准，动力性评 价试验均在满载情况下进行。
M文件
发动机输出功率计算，具体程序如下：
for k=1:5 u=ugk(k):ukm(k); n=ig(k)*i0*u/r/0.377; pe=polyval(dp,n); figure(4) plot(u,pe) hold on,grid on%保证每次循环的图形都保留显示
end
M文件
发动机阻力功率计算，具体程序如下：
M文件
图解法求解最高车速，具体程序如下：
for u=50:0.1:180; k=5; n=ig(k)*i0*u/r/0.377; t=polyval(dt,n); Ft=t*ig(k)*i0*eta/r; f=f0+f1*(u/100)+f4*(u/100).^4; Ff=G*f; %计算滚动阻力 Fw=Cd*A*u.^2./21.15; %计算空气阻力 F=Ff+Fw; if abs(Ft-F)<1; %当驱动力与行驶阻力差值小于1N时，近似认为为理论最高车速 umax=u;break end
仍然取先前国产某小型轿车进行仿真计算，整车综合性 能参数如前所述。发动机输出功率特性可由台架试验获得， 如下表2所示：
表2 发动机输出功率特性
转矩/(r/min) 1500 2000 2500 3000 3500 4000 4500 5000 5500 输出功率/kW 10.6 13.4 17.6 21.8 25.2 27.5 31.1 35.4 31.6
毕业设计(论文)-基于Matlab的汽车动力性与经济性分析计算
提供全套毕业论文图纸，欢迎咨询编号毕业设计（论文）题目基于Matlab的汽车动力性与经济性分析计算二级学院专业车辆工程班级学生姓名学号指导教师职称时间目录摘要 (I)Abstract (II)1绪论 (1)1.1选题背景 (1)1.2国内外汽车动力性经济性仿真研究发展过程与现状 (2)1.3课题研究主要内容与意义 (3)2汽车动力性经济性计算中发动机模拟的理论基础 (5)2.1发动机数学模型的建立 (5)2.1.1发动机外特性 (5)2.1.2发动机万有特性 (8)2.2本章小结 (10)3汽车动力性模拟计算 (12)3.1汽车最高车速的计算 (13)3.2加速时间的计算 (16)3.2.1原地起步加速时间 (17)3.2.2超车加速时间 (19)3.3最大爬坡度的计算 (19)3.4各档动力因数的计算 (21)3.5小结 (23)4汽车燃油经济性模拟计算 (24)4.1汽车燃油经济性的评价指标 (24)4.2不同行驶工况的汽车燃油经济性计算 (25)4.2.1等速工况 (25)4.2.2加速工况 (25)4.2.3减速工况 (27)4.2.4怠速工况 (27)4.2.5多工况循环百公里油耗 (27)4.3等速行驶工况的汽车燃油经济性计算 (27)4.4影响汽车燃油经济性的因素 (31)4.5小结 (38)5动力性和燃油经济性的参数敏感性分析 (39)5.1汽车动力性的参数敏感度分析 (39)5.2汽车燃油经济性的灵敏度分析 (40)5.3本章小结 (42)6结论 (44)致谢 (46)参考文献 (47)附录 (49)附录一 (49)附录二 (51)附录三 (62)附录四 (68)文献综述 (85)摘要本文首先利用了最小二乘法和回归分析法对所给试验数据拟合，得到了一个函数化的发动机外特性模型和万有特性模型，为后文汽车动力与燃油经济性的仿真奠定了基础。
然后，我们建立了汽车的动力性数学模型，详细分析了汽车动力性的几个评价指标，给出各个指标的计算方法及计算公式，并绘制出汽车的动力特性图，在完成动力性仿真分析之后，我们建立了汽车燃油经济性的数学模型，分析不同工况下汽车的燃油经济性并重点分析等速百公里的燃油消耗量，得到不同档位下不同车速的等速百公里燃油消耗量，并从汽车的使用方面和汽车结构方面详细的分析了影响汽车燃油经济性的因素。
汽车性能仿真计算实验实验报告
汽车性能仿真计算实验实验报告实验⼀汽车动⼒性仿真计算实验⽬的1.掌握汽车动⼒性评价指标和评价⽅法2.学会使⽤matlab 对汽车动⼒性指标进⾏计算实验内容1.学习汽车动⼒性理论2.编写计算程序3.绘制汽车动⼒性图形实验设备硬件环境：汽车虚拟仿真实验室软件环境：matlab2016a 及以上版本实验步骤1.学习汽车动⼒性理论2.编写计算程序3.绘制汽车动⼒性图形实验报告1. 运⽤matlab 解决《汽车理论》第⼀章习题1.31）绘制汽车驱动⼒与⾏驶阻⼒平衡图汽车驱动⼒Ft=ri i T to g tq η⾏驶阻⼒F f +F w ＋F i +F j ＝G ?f +2D 21.12A C a u +G ?i+dt dum δ发动机转速与汽车⾏驶速度之间的关系式为：0g i nr 0.377ua i ?= 由本题的已知条件,即可求得汽车驱动⼒和⾏驶阻⼒与车速的关系,编程即可得到汽车驱动⼒与⾏驶阻⼒平衡图。
2）求汽车最⾼车速，最⼤爬坡度及克服该坡度时相应的附着率①由1）得驱动⼒与⾏驶阻⼒平衡图，汽车的最⾼车速出现在5档时汽车的驱动⼒曲线与⾏驶阻⼒曲线的交点处，Ua max ＝99.08m/s 2。
②汽车的爬坡能⼒，指汽车在良好路⾯上克服w f F F +后的余⼒全部⽤来（等速）克服坡度阻⼒时能爬上的坡度，此时0=dt du，因此有()w f t i F F F F +-=，可得到汽车爬坡度与车速的关系式：()+-=G F F F i w f t arcsin tan ；⽽汽车最⼤爬坡度为Ⅰ档时的最⼤爬坡度。
利⽤MATLAB 计算可得，352.0max =i 。
③如是前轮驱动，1?C ＝qb hg q L L -；相应的附着率1?C 为1.20，不合理，舍去。
如是后轮驱动，2?C ＝qa hg q L L+；相应的附着率2?C 为0.50。
3）绘制汽车⾏驶加速度倒数曲线，求加速时间利⽤MATLAB 画出汽车的⾏驶加速度图和汽车的加速度倒数曲线图：忽略原地起步时的离合器打滑过程，假设在初时刻时，汽车已具有Ⅱ档的最低车速。
基于MATLAB的汽车动力性及燃油经济性的计算机仿真
基于MATLAB的汽车动力性及燃油经济性的计算机仿真作者：张静来源：《科技视界》2012年第23期【摘要】建立汽车动力性和燃油经济性仿真计算的数学模型，利用MATLAB计算机仿真程序，为车辆的前期设计提供理论计算依据。
【关键词】动力性；经济性；MATLAB仿真０引言动力性和燃油经济性是汽车性能的重要指标，石油价格的飞速上涨，对汽车性能有了更高的要求。
动力性和燃油经济性的计算机仿真能准确、快速、有效的预测性能指标。
节省实车试验中不必要的巨额浪费及实车道路试验中驾驶员、道路环境、气候等因素对汽车使用性能测定的影响，在新车设计中迅速且经济地选择最佳方案。
1发动机数学模型发动机数学模型是整车动力性和燃油经济性仿真计算的重要依据，包括外特性数学模型和万有特性数学模型。
本文以发动机台架实验数据为依据，采用插值法描述发动机万有特性；采用最小二乘法曲线拟合描述发动机外特性。
1.1发动机的外特性在进行动力性估算时，一般仍沿用稳态工况时发动机台架试验所得到的使用外特性中的功率与转矩曲线。
稳定工况时发动机转矩曲线基本呈抛物线形状，并且为转速的一元函数，所以采用最小二乘法曲线拟合法描述。
Ｔ■＝ａ■＋ａ■ｎ＋ａ■ｎ■＋…＋ａ■ｎ■式中n为发动机转速（r/min）；T■为稳定工况下发动机转矩（N·m）；系数a■,a■,a■, …,a■可由最小二乘法来确定；拟合阶数k随特性曲线而异，一般在2、3、4、5中选取。
1.2发动机万有特性发动机的万有特性是个二元函数，燃油消耗率b是发动机转速n和功率p的函数，国内外多采用试验数据的矩阵描述方法，需要时插值提取；国内也有采用曲面拟合法的，但要小心其系数矩阵可能出现病态。
本文采用插值法比较迅速，精度的高低取决于数据点的疏密程度，可真实的反映万有特性的局部特点。
二元插值公式：ｂ＝ｂ（ｎ，ｐ）＝■■■■?鄢ｂ（ｎ■，ｐ■）?鄢■■其中n■，p■，b(n■，p■)为给定的万有特性上的节点。
汽车动力性计算matlab程序
%% 汽车动力性计算（自己编的动力性计算程序，供大家计算动力性时参考，具体参数大家根据所给程序对应输入，并对坐标轴数值按需要进行修改）clc; clear;close all;%%根据所给发动机数据拟合外特性曲线（发动机数据按照你所得到的数据进行输入）n_test=[500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 2200];T_test=[975 1108 1298 1496 1546 1620 1670 1785 1974 1974 1970 1889 1829 1748 1669 1700 1524 1105];figure(1)plot(n_test,T_test,'g');hold ongrid on%p=polyfit(n_test,T_test,7);p=polyfit(n_test,T_test,2);n=[450:1:2200];Ttq=polyval(p,n);plot(n,Ttq,'k');xlabel('发动机转速n(r/min)');ylabel('发动机转矩Ttq(N*m)');title('发动机转矩曲线');legend('测试曲线','拟合曲线');%%所给车型动力总成相关参数ig=[3.07 2.16 1.48 1.0 0.82];i0=4.0; eta=0.78; r=0.57; M=25000; g=9.8; c=1.5; f0=0.01; f1=0.0002; f4=0.0005; CD=1; A=8;Iw=3.6;If=0.04;%% 发动机外特性曲线图figure(2)hold ongrid onfor i=length(n);Pe=Ttq.*n/9550;end[AX,H1,H2]=plotyy(n,Ttq,n,Pe);xlabel('发动机转速n(r/min)');ylabel('发动机转矩Ttq(N*m)');ylabel(AX(2),'发动机功率Pe(Kw)');title('发动机外特性曲线');%% 各挡位速度曲线%计算各挡位车速for i=1:length(ig);ua(i,:)=0.377*r*n/ig(i)/i0;end%计算各档位最高车速uamax=ua(:,length(ua(1,:)));figure(3)hold onfor i=1:length(ig);plot(n,ua(i,:),'k');endhold ongrid onxlabel('转速n(r/min)');ylabel('各挡位车速(km/h)');title('各挡位车速-转速表');legend('1挡车速','2挡车速','3挡车速','4挡车速','5挡车速'); %% 驱动力和行驶阻力平衡图%计算滚动阻力系数for i=1:length(ig);f(i,:)=f0+f1*(ua(i,:)/100)+f4*(ua(i,:)/100).^4;end%计算滚动阻力for i=1:length(ig);Ff(i,:)=c*M*g*f(i,:);end%计算空气阻力for i=1:length(ig);Fw(i,:)=CD*A*(ua(i,:).^2)/21.15; end%计算行驶阻力for i=1:length(ig);F(i,:)=Ff(i,:)+Fw(i,:);end%计算汽车驱动力for i=1:length(ig);Ft(i,:)=Ttq*ig(i)*i0*eta/r;endfigure(4)hold onfor i=1:length(ig);plot(ua(i,:), Ft(i,:),'k');plot(ua(i,:), F(i,:),'r');plot(ua(i,:), Ff(i,:),'b');endhold ongrid onxlabel('车速(km/h)');ylabel('驱动力、行驶阻力(N)');legend('驱动力Ft','行驶阻力Ff+Fw','滚动阻力Ff'); title('驱动力-行驶阻力平衡图');%% 汽车功率平衡图%计算各档位功率for i=1:length(ig);P(i,:)=Ft(i,:).*ua(i,:)/(3600*eta);end%计算风阻阻力功率for i=1:length(ig);Pw(i,:)=CD*A*ua(i,:).^3/(76140*eta);end%计算滚动阻力功率for i=1:length(ig);Pf(i,:)=M*g*f(i,:).*ua(i,:)/(3600*eta);end%计算总阻力功率for i=1:length(ig);Pz(i,:)=Pw(i,:)+Pf(i,:);endfigure(5)hold onfor i=1:length(ig);plot(ua(i,:), P(i,:),'k');plot(ua(i,:), Pz(i,:),'r');endhold ongrid onxlabel('车速(km/h)');ylabel('发动机功率、阻力功率(kW)');legend('发动机功率P','阻力功率Pz','Location','NorthWest'); title('功率平衡图');%% 动力特性图（动力因数图）for i=1:length(ig);D(i,:)= (Ft(i,:)- Fw(i,:))/M/g;endfigure(6)hold onfor i=1:length(ig);plot(ua(i,:), D(i,:),'k');plot(ua(i,:), f(i,:),'r');endhold ongrid onxlabel('车速(km/h)');ylabel('动力因数D');legend('动力因数D','滚动阻力系数f');title('动力特性图');%% 爬坡度曲线图for i=1:length(ig);I(i,:)= (tan(asin((Ft(i,:)-(Ff(i,:)+Fw(i,:)))/(M*g))))*100; endfigure(7)hold onfor i=1:length(ig);if i==1plot(ua(i,:),I(i,:),'r');elseplot(ua(i,:),I(i,:),'k');endendhold ongrid onxlabel('车速(km/h)');ylabel('最大爬坡度（%）');legend('Ⅰ挡','高速档');title('爬坡度曲线图');%% 加速度曲线图deta=1+1/M*4*Iw/r^2+1/M*If*ig.^2*i0^2*eta/r^2; for i=1:length(ig);a(i,:)=(Ft(i,:)-Ff(i,:)-Fw(i,:))./deta(i)/M;if i==5for j=1:length(n)if a(i,j)<0a(i,j)=0;elseendendendendfigure(8)hold onfor i=1:length(ig);if i==1plot(ua(i,:),a(i,:),'r'); elseplot(ua(i,:),a(i,:),'k');endendhold ongrid onxlabel('车速(km/h)');ylabel('加速度a（m/s^2）');legend('Ⅰ档','高速档');title('加速度曲线图');axis([0 120 0 1.5])%% 加速度倒数曲线for i=1:length(ig);for j=1:length(n)b(i,j)=1./a(i,j);endendfigure(9)hold onfor i=1:length(ig)plot(ua(i,:),b(i,:),'k');endhold ongrid onxlabel('车速(km/h)');ylabel('各档加速度倒数1/a');legend('各档加速度倒数1/a曲线','Location','NorthWest'); title('各档加速度倒数曲线图');axis([0 120 0 10])ad1=b(1,:);ad2=ua(1,:);for i=1:(length(ig)-1);for j=1:length(n)if ua(i+1,j)>=ua(i,length(n)) flag(i)=j;break;endendad1=[ad1 b(i+1,j:length(n))]; ad2=[ad2 ua(i+1,j:length(n))];endfigure(10)hold onplot(ad2,ad1,'k');hold ongrid onxlabel('车速(km/h)');ylabel('加速度倒数1/a');legend('加速度倒数1/a曲线','Location','NorthWest');title('加速度倒数曲线图');axis([0 120 0 10])%% 加速时间曲线k=length(n);for i=1:length(ig);t(i,1)=0;for j=2:kt(i,j)=abs(ua(i,j)-ua(i,j-1))*(b(i,j)+b(i,j-1))/2;endendfor i=1:length(ig);for j=1:kat(i,j)=sum(t(i,1:j))/3.6;endendtotalat=at(1,:);for i=1:(length(ig)-1);for j=flag(i):ktotalat=[totalat totalat(length(totalat))+t(i+1,j)/3.6];endendfigure(11)hold onplot(totalat,ad2,'k');hold ongrid onxlabel('时间（s）');ylabel('车速(km/h)');legend('加速时间','Location','NorthWest'); title('加速时间曲线图');axis([0 100 0 120])Welcome To Download !!!欢迎您的下载，资料仅供参考！。
汽车理论课后习题MATLAB编程-武汉理工版
汽车理论课后习题MATLAB编程1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h的加速时间。
解：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];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;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,Ft1,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('汽车最高车速=');disp(x);disp('km/h');汽车最高车速=99.3006km/h(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];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;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('%');汽车最大爬坡度=35.2197%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];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;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/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r^2)+(If*ig(1)^2*i0^2*nT)/(m*r^2);deta2=1+(Iw1+Iw2)/(m*r^2)+(If*ig(2)^2*i0^2*nT)/(m*r^2);deta3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*nT)/(m*r^2);deta4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*nT)/(m*r^2);deta5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*nT)/(m*r^2);a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./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/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);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);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:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^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/100 0)^4;Ft=Tq*ig(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/100 0)^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;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^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;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%% 汽车动力性计算（自己编的动力性计算程序，供大家计算动力性时参考，具体参数大家根据所给程序对应输入，并对坐标轴数值按需要进行修改）clc; clear;close all;%%根据所给发动机数据拟合外特性曲线（发动机数据按照你所得到的数据进行输入）n_test=[500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 2200];T_test=[975 1108 1298 1496 1546 1620 1670 1785 1974 1974 1970 1889 1829 1748 1669 1700 1524 1105];figure(1)plot(n_test,T_test,'g');hold ongrid on%p=polyfit(n_test,T_test,7);p=polyfit(n_test,T_test,2);n=[450:1:2200];Ttq=polyval(p,n);plot(n,Ttq,'k');xlabel('发动机转速n(r/min)');ylabel('发动机转矩Ttq(N*m)');title('发动机转矩曲线');legend('测试曲线','拟合曲线');%%所给车型动力总成相关参数ig=[3.07 2.16 1.48 1.0 0.82];i0=4.0; eta=0.78; r=0.57; M=25000; g=9.8; c=1.5; f0=0.01; f1=0.0002;f4=0.0005; CD=1;A=8;Iw=3.6;If=0.04;%% 发动机外特性曲线图figure(2)hold ongrid onfor i=length(n);Pe=Ttq.*n/9550;end[AX,H1,H2]=plotyy(n,Ttq,n,Pe);xlabel('发动机转速n(r/min)');ylabel('发动机转矩Ttq(N*m)');ylabel(AX(2),'发动机功率Pe(Kw)');title('发动机外特性曲线');%% 各挡位速度曲线%计算各挡位车速for i=1:length(ig);ua(i,:)=0.377*r*n/ig(i)/i0;end%计算各档位最高车速uamax=ua(:,length(ua(1,:)));figure(3)hold onfor i=1:length(ig);plot(n,ua(i,:),'k');endhold ongrid onxlabel('转速n(r/min)');ylabel('各挡位车速(km/h)');title('各挡位车速-转速表');legend('1挡车速','2挡车速','3挡车速','4挡车速','5挡车速'); %% 驱动力和行驶阻力平衡图%计算滚动阻力系数for i=1:length(ig);f(i,:)=f0+f1*(ua(i,:)/100)+f4*(ua(i,:)/100).^4;end%计算滚动阻力for i=1:length(ig);Ff(i,:)=c*M*g*f(i,:);end%计算空气阻力for i=1:length(ig);Fw(i,:)=CD*A*(ua(i,:).^2)/21.15;end%计算行驶阻力for i=1:length(ig);F(i,:)=Ff(i,:)+Fw(i,:);end%计算汽车驱动力for i=1:length(ig);Ft(i,:)=Ttq*ig(i)*i0*eta/r;endfigure(4)hold onfor i=1:length(ig);plot(ua(i,:), Ft(i,:),'k');plot(ua(i,:), F(i,:),'r');plot(ua(i,:), Ff(i,:),'b');endhold ongrid onxlabel('车速(km/h)');ylabel('驱动力、行驶阻力(N)');legend('驱动力Ft','行驶阻力Ff+Fw','滚动阻力Ff');title('驱动力-行驶阻力平衡图');%% 汽车功率平衡图%计算各档位功率for i=1:length(ig);P(i,:)=Ft(i,:).*ua(i,:)/(3600*eta);end%计算风阻阻力功率for i=1:length(ig);Pw(i,:)=CD*A*ua(i,:).^3/(76140*eta);end%计算滚动阻力功率for i=1:length(ig);Pf(i,:)=M*g*f(i,:).*ua(i,:)/(3600*eta);end%计算总阻力功率for i=1:length(ig);Pz(i,:)=Pw(i,:)+Pf(i,:);endfigure(5)hold onfor i=1:length(ig);plot(ua(i,:), P(i,:),'k');plot(ua(i,:), Pz(i,:),'r');endhold ongrid onxlabel('车速(km/h)');ylabel('发动机功率、阻力功率(kW)');legend('发动机功率P','阻力功率Pz','Location','NorthWest'); title('功率平衡图');%% 动力特性图（动力因数图）for i=1:length(ig);D(i,:)= (Ft(i,:)- Fw(i,:))/M/g;endfigure(6)hold onfor i=1:length(ig);plot(ua(i,:), D(i,:),'k');plot(ua(i,:), f(i,:),'r');endhold ongrid onxlabel('车速(km/h)');ylabel('动力因数D');legend('动力因数D','滚动阻力系数f');title('动力特性图');%% 爬坡度曲线图for i=1:length(ig);I(i,:)= (tan(asin((Ft(i,:)-(Ff(i,:)+Fw(i,:)))/(M*g))))*100; endfigure(7)hold onfor i=1:length(ig);if i==1plot(ua(i,:),I(i,:),'r');elseplot(ua(i,:),I(i,:),'k');endendhold ongrid onxlabel('车速(km/h)');ylabel('最大爬坡度（%）');legend('Ⅰ挡','高速档');title('爬坡度曲线图');%% 加速度曲线图deta=1+1/M*4*Iw/r^2+1/M*If*ig.^2*i0^2*eta/r^2;for i=1:length(ig);a(i,:)=(Ft(i,:)-Ff(i,:)-Fw(i,:))./deta(i)/M;if i==5for j=1:length(n)if a(i,j)<0a(i,j)=0;elseendendendendfigure(8)hold onfor i=1:length(ig);if i==1plot(ua(i,:),a(i,:),'r');elseplot(ua(i,:),a(i,:),'k');endendhold ongrid onxlabel('车速(km/h)');ylabel('加速度a（m/s^2）');legend('Ⅰ档','高速档');title('加速度曲线图');axis([0 120 0 1.5])%% 加速度倒数曲线for i=1:length(ig);for j=1:length(n)b(i,j)=1./a(i,j);endendfigure(9)hold onfor i=1:length(ig)plot(ua(i,:),b(i,:),'k');endhold ongrid onxlabel('车速(km/h)');ylabel('各档加速度倒数1/a');legend('各档加速度倒数1/a曲线','Location','NorthWest'); title('各档加速度倒数曲线图');axis([0 120 0 10])ad1=b(1,:);ad2=ua(1,:);for i=1:(length(ig)-1);for j=1:length(n)if ua(i+1,j)>=ua(i,length(n))flag(i)=j;break;endendad1=[ad1 b(i+1,j:length(n))];ad2=[ad2 ua(i+1,j:length(n))];endfigure(10)hold onplot(ad2,ad1,'k');hold ongrid onxlabel('车速(km/h)');ylabel('加速度倒数1/a');legend('加速度倒数1/a曲线','Location','NorthWest');title('加速度倒数曲线图');axis([0 120 0 10])%% 加速时间曲线k=length(n);for i=1:length(ig);t(i,1)=0;for j=2:kt(i,j)=abs(ua(i,j)-ua(i,j-1))*(b(i,j)+b(i,j-1))/2;endendfor i=1:length(ig);for j=1:kat(i,j)=sum(t(i,1:j))/3.6;endendtotalat=at(1,:);for i=1:(length(ig)-1);for j=flag(i):ktotalat=[totalat totalat(length(totalat))+t(i+1,j)/3.6]; endendfigure(11)hold onplot(totalat,ad2,'k');hold ongrid onxlabel('时间（s）');ylabel('车速(km/h)');legend('加速时间','Location','NorthWest'); title('加速时间曲线图');axis([0 100 0 120])。