汽车动力性计算matlab程序Word格式.docx
《汽车动力性计算matlab程序Word格式.docx》由会员分享，可在线阅读，更多相关《汽车动力性计算matlab程序Word格式.docx（9页珍藏版）》请在冰豆网上搜索。
发动机转速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）;
ylabel（AX
（2）,'
发动机功率Pe（Kw）'
发动机外特性曲线'
%%各挡位速度曲线
%计算各挡位车速
fori=1:
length（ig）;
ua（i,:
）=0.377*r*n/ig（i）/i0;
%计算各档位最高车速
uamax=ua（:
length（ua（1,:
）））;
figure（3）
plot（n,ua（i,:
）,'
转速n（r/min）'
各挡位车速（km/h）'
各挡位车速-转速表'
1挡车速'
2挡车速'
3挡车速'
4挡车速'
5挡车速'
%%驱动力和行驶阻力平衡图
%计算滚动阻力系数
f（i,:
）=f0+f1*（ua（i,:
）/100）+f4*（ua（i,:
）/100）.A4;
end%计算滚动阻力
Ff（i,:
）=c*M*g*f（i,:
%计算空气阻力
Fw（i,:
）=CD*A*（ua（i,:
）.A2）/21.15;
%计算行驶阻力
F（i,:
）=Ff（i,:
）+Fw（i,:
%计算汽车驱动力
Ft（i,:
）=Ttq*ig（i）*i0*eta/r;
figure（4）
plot（ua（i,:
）,Ft（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,:
）.A3/（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）
holdonfori=1:
ifi==1
）,I（i,:
else
最大爬坡度（%）'
I捋高速档'
爬坡度曲线图'
%%加速度曲线图
deta=1+1/M*4*Iw/rA2+1/M*If*ig.A2*i0A2*eta/rA2;
a（i,:
）=（Ft（i,:
）-Ff（i,:
）-Fw（i,:
））./deta（i）/M;
ifi==5
forj=1:
length（n）
irawAo
aw）占
e-seendendendendfigured）
ho-donfornl」engfh（ig=
p_of（ua（i..）£
D（i..）<
=
e-se
p_of（ua（L）a（L）-3
ho-don
x_abe_（->
B（km、h））
y-abe-（-nuBMa（m、s>
2））
I档高速档'
加速度曲线图'
axis（[012001.5]）
%%加速度倒数曲线
b（i,j）=1./a（i,j）;
figure（9）
length（ig）
）,b（i,:
车速（km/h）'
各档加速度倒数1/a'
各档加速度倒数1/a曲线'
各档加速度倒数曲线图'
axis（[0120010]）ad1=b（1,:
ad2=ua（1,:
（length（ig）-1）;
ifua（i+1,j）>
=ua（i,length（n））
flag（i）=j;
break;
ad1=[ad1b（i+1,j:
length（n））];
ad2=[ad2ua（i+1,j:
figure（10）
plot（ad2,ad1,'
加速度倒数1/a'
legend（'
加速度倒数1/a曲线'
title（'
加速度倒数曲线图'
axis（[0120010]）
%%加速时间曲线
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];
endend
figure（11）
plot（totalat,ad2,'
时间（s）'
加速时间'
加速时间曲线图'
axis（[01000120]）
WelcomeTo
Download!
!
欢迎您的下载，资料仅供参考！
升级会员 