基于MATLAB编程绘制动力特性图分析汽车动力性

10.16638/ki.1671-7988.2017.13.020

基于MATLAB编程绘制动力特性图分析汽车动力性

尹佣博

（武汉理工大学汽车工程学院，湖北武汉430070 ）

摘要：汽车动力性主要由最高车速、加速时间、最大爬坡度三个指标来评价。汽车动力性试验主要包括最高车速、起步连续换挡加速与超车加速和汽车最大爬坡度三大内容。上述试验项目需要特定的试验路段，且耗时较长。基于MA TLAB编程可以绘出发动机外特性曲线图和动力特性图，从动力特性图的角度研究汽车的动力性，不仅可以准确获得最高车速、加速时间和最大爬坡度等指标，还可以通过图像观察相关变化趋势，为进一步研究汽车的动力性提供了参考。

关键词：汽车动力性；MATLAB编程；动力特性图

中图分类号：U461.2 文献标识码：A 文章编号：1671-7988 (2017)13-65-04

Drawing dynamic characteristic diagram based on MATLAB programming to

analyze automobile power

Yin Yongbo

( Wuhan University of Technology. School of Automotive Engineering, Hubei Wuhan 430070 )

Abstract:Vehicle dynamics mainly by the highest speed, acceleration time, the maximum climbing three indicators to evaluate. Vehicle power test includes the maximum speed, starting continuous shift acceleration and overtaking acceleration and the maximum car climbing three content. The pilot project requires a specific test section, and take a long time. Based on the MA TLAB programming, we can plot the engine characteristic curve and the dynamic characteristic diagram. From the view of the dynamic characteristic diagram, we can study the dynamic performance of the vehicle, not only can get the maximum vehicle speed, acceleration time and maximum climbing index, but also observe the correlation And provides a reference for further research on the dynamics of the vehicle.

Keywords: Vehicle dynamics; MATLAB programming; Dynamic characteristics map

CLC NO.: U461.2 Document Code: A Article ID: 1671-7988 (2017)13-65-04

前言

汽车是高效率的交通工具，在汽车的众多性能中，动力性是汽车最重要的性能。传统动力性试验不仅试验项目多、试验周期长并且要特定的试验路段，基于MATLAB编程绘制动力特性图来分析汽车动力性，不用进行动力学试验就可以获得动力性能试验的相关结果，计算结果可靠，被广泛采用。同时为相关工作人员进一步分析汽车的动力性提供了一些参考。

1 汽车的动力性评价指标

汽车动力性评价指标有最高车速、加速时间和最大爬坡度等，相应的试验内容有最高车速、汽车起步连续换挡加速与超车加速以及汽车最大爬坡度等。货车经常在不同地区的各种道路上行驶，必须具有足够的爬坡能力，一般最大爬坡

作者简介：尹佣博，就读于武汉理工大学汽车工程学院。

尹佣博：基于MATLAB编程绘制动力特性图分析汽车动力性66 2017年第13期

度在30%即16.7°左右。根据汽车行驶时在纵向受到的驱动力和行驶阻力的平衡关系可建立汽车行驶方程式，整理行驶方程式得到动力因素，绘制动力特性图来得到最高车速、加速时间和最大爬坡度来分析汽车的动力性[1]。

2 汽车的最高车速、加速时间、最大爬坡度求解

2.1 汽车行驶方程式及动力因素

根据汽车纵向受力平衡关系建立汽车行驶方程式：

行驶方程式具体化：

式中：F t为驱动力（N）；T ta为发动机输出转矩（N·m）；

i g、i0为变速器传动比、主减速器传动比；ηt为传动系的机械效率；r为车轮滚动半径（m）；F f为滚动阻力（N）；G为汽车重力（N）；f为滚动阻力系数；F w为空气阻力（N）；C D为空气阻力系数；A为迎风面积(m2）；u a为车速（km/h）；F i 为坡度阻力（N）；α为坡道角；F j为加速阻力（N）；δ为汽车旋转质量换算系数；m为质量（kg）；i为道路坡度；D为动力因素；ψ为道路阻力系数g为重力加速度，取9.8m/s2。

2.2 汽车发动机特性参数及货车的特性参数

所研究的汽油发动机使用外特性T ta-n曲线的拟合公式为：

所选择的发动机最低转速600r/min，最高转速4000r/ min。

所研究货车的动力参数如：装载质量2000kg，整车整备质量1800kg，总质量3800kg，车轮半径0.367m，传动系机械效率0.85，滚动阻力系数0.013，空气阻力×迎风面积2.77m2，主减速器传动比5.83，飞轮转动惯量0.218Kg.m2。变速器的传动比如下：一档5.56，二档2.769，三档1.644，四档1，五档0.793。所研究的货车的滚动阻力系数随车速的变化关系为：

2.3 汽车的最高车速

最高车速时，加速阻力和坡道阻力为0。在动力特性图上滚动阻力曲线和动力因素曲线的交点是最高车速。每一档变速器的传动比确定之后可以根据转速范围得到每一档的车速范围：一档变速器传动比5.56，车速范围[2.561 17.074]；二档变速器传动比2.769，车速范围[5.142 34.283]；三档变速器传动比1.644，车速范围[8.661 57.74]；四档变速器传动比1，车速范围[14.239 94.929]；五档变速器传动比0.793，车速范围[17.956 119.709]。用MA TLAB编程，得到动力特性图：clear；

figure(1)；

m=3880; %汽车总质量

r=0.367; %车轮半径

Nt=0.85; %机械传动效率

f=0.013; %滚动阻力系数

CdA=2.77; %空气阻力系数与迎风面积的乘积

i0=5.83; %主减速器传动比

If=0.218; %飞轮转动惯量

Iw=5.396; %车轮转动惯量之和

ig1=5.56;%一档驱动力

Ua1=linspace(2.561,17.074,10000);%一档车速

n1=i0*ig1*Ua1/0.377/0.367;%一档转速

Ttq1=-19.313+295.27*(n1/1000)-165.44*(n1/1000).^2+40. 874*(n1/1000).^3-3.8445*(n1/1000).^4;%一档转矩

Ft1=Ttq1*ig1*i0*Nt/r; %一档驱动力

Fw1=CdA*Ua1.^2/21.15;%一档空气阻力

D1=(Ft1-Fw1)/(m*9.8);%一档动力因素

ig2=2.769; %二档驱动力

Ua2=linspace(5.142,34.283,10000); %二档车速

n2=i0*ig2*Ua2/0.377/0.367; %二档转速

Ttq2=-19.313+295.27*(n2/1000)-165.44*(n2/1000).^2+40. 874*(n2/1000).^3-3.8445*(n2/1000).^4; %二档转矩

Ft2=Ttq2*ig2*i0*Nt/r; %二档驱动力

Fw2=CdA*Ua2.^2/21.15;%二档空气阻力

D2=(Ft2-Fw2)/(m*9.8);%二档动力因素

ig3=1.644; %三档驱动力

Ua3=linspace(8.661,57.74,10000);%三档车速

n3=i0*ig3*Ua3/0.377/0.367; %三档转速

Ttq3=-19.313+295.27*(n3/1000)-165.44*(n3/1000).^2+40. 874*(n3/1000).^3-3.8445*(n3/1000).^4; %三档转矩

Ft3=Ttq3*ig3*i0*Nt/r; %三档驱动力

Fw3=CdA*Ua3.^2/21.15;%三档空气阻力

D3=(Ft3-Fw3)/(m*9.8);%三档动力因素

ig4=1.0; %四档驱动力

Ua4=linspace(14.239,94.929,10000); %四档车速

n4=i0*ig4*Ua4/0.377/0.367; %四档转速

Ttq4=-19.313+295.27*(n4/1000)-165.44*(n4/1000).^2+40. 874*(n4/1000).^3-3.8445*(n4/1000).^4; %四档转矩

Ft4=Ttq4*ig4*i0*Nt/r; %四档驱动力

Fw4=CdA*Ua4.^2/21.15;%四档空气阻力

D4=(Ft4-Fw4)/(m*9.8);%四档动力因素

ig5=0.793; %五档驱动力

Ua5=linspace(17.956,119.709,10000); %五档车速

n5=i0*ig5*Ua5/0.377/0.367;%五档转速

Ttq5=-19.313+295.27*(n5/1000)-165.44*(n5/1000).^2+40.