汽车动力性计算matlab程序Word文档下载推荐.docx
- 文档编号:320943
- 上传时间:2023-04-28
- 格式:DOCX
- 页数:10
- 大小:15.45KB
汽车动力性计算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])
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 汽车 动力 计算 matlab 程序