汽车动力性计算matlab程序Word文档下载推荐.docx
《汽车动力性计算matlab程序Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《汽车动力性计算matlab程序Word文档下载推荐.docx（10页珍藏版）》请在冰点文库上搜索。
发动机转速n（r/min）'
ylabel（'
发动机转矩Ttq（N*m）'
title（'
发动机转矩曲线'
legend（'
测试曲线'
'
拟合曲线'
%%所给车型动力总成相关参数
ig=[3.072.161.481.00.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）
fori=length（n）;
Pe=Ttq.*n/9550;
end
[AX,H1,H2]=plotyy（n,Ttq,n,Pe）;
xlabel（'
ylabel（AX
（2）,'
发动机功率Pe（Kw）'
发动机外特性曲线'
%%各挡位速度曲线
%计算各挡位车速
fori=1:
length（ig）;
ua（i,:
）=0.377*r*n/ig（i）/i0;
end
%计算各档位最高车速
uamax=ua（:
length（ua（1,:
）））;
figure（3）
plot（n,ua（i,:
）,'
转速n（r/min）'
各挡位车速（km/h）'
title（'
各挡位车速-转速表'
1挡车速'
2挡车速'
3挡车速'
4挡车速'
5挡车速'
%%驱动力和行驶阻力平衡图
%计算滚动阻力系数
f（i,:
）=f0+f1*（ua（i,:
）/100）+f4*（ua（i,:
）/100）.^4;
%计算滚动阻力
Ff（i,:
）=c*M*g*f（i,:
%计算空气阻力
Fw（i,:
）=CD*A*（ua（i,:
）.^2）/21.15;
%计算行驶阻力
F（i,:
）=Ff（i,:
）+Fw（i,:
%计算汽车驱动力
Ft（i,:
）=Ttq*ig（i）*i0*eta/r;
figure（4）
holdon
plot（ua（i,:
）,Ft（i,:
plot（ua（i,:
）,F（i,:
r'
）,Ff（i,:
b'
车速（km/h）'
驱动力、行驶阻力（N）'
驱动力Ft'
行驶阻力Ff+Fw'
滚动阻力Ff'
驱动力-行驶阻力平衡图'
%%汽车功率平衡图
%计算各档位功率
P（i,:
）=Ft（i,:
）.*ua（i,:
）/（3600*eta）;
%计算风阻阻力功率
Pw（i,:
）=CD*A*ua（i,:
）.^3/（76140*eta）;
%计算滚动阻力功率
Pf（i,:
）=M*g*f（i,:
%计算总阻力功率
Pz（i,:
）=Pw（i,:
）+Pf（i,:
figure（5）
）,P（i,:
）,Pz（i,:
发动机功率、阻力功率（kW）'
发动机功率P'
阻力功率Pz'
Location'
NorthWest'
功率平衡图'
%%动力特性图（动力因数图）
D（i,:
）=（Ft（i,:
）-Fw（i,:
））/M/g;
figure（6）
）,D（i,:
）,f（i,:
动力因数D'
滚动阻力系数f'
动力特性图'
%%爬坡度曲线图
I（i,:
）=（tan（asin（（Ft（i,:
）-（Ff（i,:
）））/（M*g））））*100;
figure（7）
ifi==1
）,I（i,:
else
最大爬坡度（%）'
Ⅰ挡'
高速档'
爬坡度曲线图'
%%加速度曲线图
deta=1+1/M*4*Iw/r^2+1/M*If*ig.^2*i0^2*eta/r^2;
a（i,:
）=（Ft（i,:
）-Ff（i,:
）-Fw（i,:
））./deta（i）/M;
ifi==5
forj=1:
length（n）
ifa（i,j）<
0
a（i,j）=0;
else
end
end
figure（8）
plot（ua（i,:
）,a（i,:
gridon
加速度a（m/s^2）'
Ⅰ档'
加速度曲线图'
axis（[012001.5]）
%%加速度倒数曲线
forj=1:
length（n）
b（i,j）=1./a（i,j）;
figure（9）
length（ig）
plot（ua（i,:
）,b（i,:
各档加速度倒数1/a'
各档加速度倒数1/a曲线'
各档加速度倒数曲线图'
axis（[0120010]）
ad1=b（1,:
ad2=ua（1,:
（length（ig）-1）;
length（n）
ifua（i+1,j）>
=ua（i,length（n））
flag（i）=j;
break;
end
ad1=[ad1b（i+1,j:
length（n））];
ad2=[ad2ua（i+1,j:
figure（10）
plot（ad2,ad1,'
加速度倒数1/a'
加速度倒数1/a曲线'
加速度倒数曲线图'
%%加速时间曲线
k=length（n）;
t（i,1）=0;
forj=2:
k
t（i,j）=abs（ua（i,j）-ua（i,j-1））*（b（i,j）+b（i,j-1））/2;
at（i,j）=sum（t（i,1:
j））/3.6;
totalat=at（1,:
forj=flag（i）:
totalat=[totalattotalat（length（totalat））+t（i+1,j）/3.6];
figure（11）
plot（totalat,ad2,'
时间（s）'
加速时间'
加速时间曲线图'
axis（[01000120]）
升级会员 