汽车理论课后作业matlab编程详解带注释范本模板Word文档下载推荐.docx
《汽车理论课后作业matlab编程详解带注释范本模板Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《汽车理论课后作业matlab编程详解带注释范本模板Word文档下载推荐.docx(15页珍藏版)》请在冰豆网上搜索。
![汽车理论课后作业matlab编程详解带注释范本模板Word文档下载推荐.docx](https://file1.bdocx.com/fileroot1/2023-2/3/e59d18fd-0f62-47df-affd-79a28621c702/e59d18fd-0f62-47df-affd-79a28621c7021.gif)
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0。
377*r*n/ig(5)/i0;
%计算各档位驱动力。
Tq=-19。
313+295.27*(n/1000)-165.44*(n/1000)。
^2+40。
874*(n/1000)。
^3—3.8445*(n/1000).^4;
Ft1=Tq*ig
(1)*i0*y/r;
Ft2=Tq*ig
(2)*i0*y/r;
Ft3=Tq*ig(3)*i0*y/r;
Ft4=Tq*ig(4)*i0*y/r;
Ft5=Tq*ig(5)*i0*y/r;
%计算行驶阻力。
Fz1=m*g*f+2.77*ua1.^2/21。
15;
Fz2=m*g*f+2.77*ua2。
^2/21。
Fz3=m*g*f+2.77*ua3.^2/21。
Fz4=m*g*f+2。
77*ua4.^2/21。
Fz5=m*g*f+2。
77*ua5.^2/21.15;
%驱动力—行驶阻力平衡图.plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua1,Fz1,’k'
ua2,Fz2,’k’,ua3,Fz3,'
k'
,ua4,Fz4,'
,ua5,Fz5,'
);
title(’驱动力-行驶阻力平衡图'
xlabel(’ua(km/s)’);
ylabel('
Ft(N)’);
legend(’Ft1’,’Ft2'
,'
Ft3'
'
Ft4’,’Ft5’,’Ff+Fw'
)
%求最高车速。
[x,y]=ginput
(1);
disp(‘汽车最高车速=’);
disp(x);
disp(’km/h'
(2)
%汽车最大爬坡度
clear
clc
%汽车相关参数.
m=3880;
g=9。
8;
nmin=600;
G=m*g;
562.7691。
6441.000.793];
85;
r=0.367;
77;
L=3.2;
947;
If=0.218;
Iw1=1.798;
Iw2=3。
598;
Tq=—19。
313+295.27*(n/1000)—165.44*(n/1000)。
^2+40.874*(n/1000).^3-3。
8445*(n/1000)。
^4;
Ft1=Tq*ig
(1)*i0*y/r;
%一档时爬坡能力最强。
377*r*n/ig
(1)/i0;
Ff=G*f;
%滚动阻力。
Fw1=CdA*ua1。
15;
%空气阻力。
Fz1=Ff+Fw1;
Fi1=Ft1—Fz1;
Fmax=max(Fi1/G);
imax=100*tan(asin(Fmax));
disp('
汽车最大爬坡度='
disp(imax),disp(’%'
);
(3)
%求汽车加速度倒数曲线和最大爬坡度显得附着率。
clear
m=3880;
%汽车相关参数。
8;
G=m*g;
562.7691.6441.000。
793];
y=0.85;
013;
i0=5.83;
L=3.2;
hg=0.9;
If=0。
218;
Iw1=1。
798;
%发动机转速范围。
313+295。
27*(n/1000)-165。
44*(n/1000)。
^2+40.874*(n/1000)。
^3-3。
8445*(n/1000)。
^4;
Ft1=Tq*ig
(1)*i0*y/r;
%各档位下的驱动力。
Ft3=Tq*ig(3)*i0*y/r;
Ft4=Tq*ig(4)*i0*y/r;
Ft5=Tq*ig(5)*i0*y/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;
^2/21.15;
%空气阻力。
Fw2=CdA*ua2.^2/21。
Fw3=CdA*ua3。
Fw4=CdA*ua4。
Fw5=CdA*ua5.^2/21。
Ff=G*f;
%滚动阻力.
q1=1+(Iw1+Iw2)/(m*r^2)+(If*ig
(1)^2*i0^2*y)/(m*r^2);
%求各档位下的汽车旋转质量换算系数.
q2=1+(Iw1+Iw2)/(m*r^2)+(If*ig
(2)^2*i0^2*y)/(m*r^2);
q3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*y)/(m*r^2);
q4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*y)/(m*r^2);
q5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*y)/(m*r^2);
a1=(Ft1-Ff—Fw1)/(q1*m);
%求出汽车加速度。
a2=(Ft2—Ff—Fw2)/(q2*m);
a3=(Ft3—Ff—Fw3)/(q3*m);
a4=(Ft4—Ff-Fw4)/(q4*m);
a5=(Ft5-Ff-Fw5)/(q5*m);
ad1=1./a1;
;
%求加速度倒数。
ad2=1。
/a2;
ad3=1./a3;
ad4=1。
/a4;
ad5=1./a5;
plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);
%绘制加速度倒数曲线。
axis([0100010]);
title(’汽车的加速度倒数曲线’);
xlabel('
ua(km/h)'
1/a’);
legend(’ad1’,'
ad2’,'
ad3'
’ad4'
ad5’)
%求最大爬坡度时的附着率.
Fz=Ft1-Ff—Fw1;
aa=asin(max(Fz)/G);
%最大爬坡角度。
C=tan(aa)/(a/L+hg*tan(aa)/L);
%附着率。
假设后轮驱动,最大爬坡度相应的附着率=);
disp(C);
(4)%求汽车2档原地起步换挡加速行驶至70km/h的加速时间
nT=0。
CDA=2.77;
83;
598;
2;
a=1.947;
ig=[5.562。
000。
793];
u1=0.377*r*nmin./ig/i0;
u2=0。
377*r*nmax./ig/i0;
deta=0*ig;
fori=1:
5
deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);
end
ua=6:
0。
01:
99;
N=length(ua);
n=0;
Tq=0;
Ft=0;
inv_a=0*ua;
delta=0*ua;
Fw=CDA*ua。
N
k=i;
ifua(i)〈=u2
(2)
n=ua(i)*(ig
(2)*i0/r)/0.377;
Tq=—19。
27*(n/1000)-165。
44*(n/1000)^2+40.874*(n/1000)^3—3。
8445*(n/1000)^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;
elseifua(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/1000)^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;
elseifua(i)<
=u2(4)
n=ua(i)*(ig(4)*i0/r)/0。
377;
Tq=-19。
27*(n/1000)—165。
44*(n/1000)^2+40。
874*(n/1000)^3—3.8445*(n/1000)^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;
else
n=ua(i)*(ig(5)*i0/r)/0.377;
313+295.27*(n/1000)—165.44*(n/1000)^2+40。
874*(n/1000)^3—3。
8445*(n/1000)^4;
Ft=Tq*ig(5)*i0*nT/r;
inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));
01*inv_a(i)/3.6;
end
a=delta(1:
k);
t(i)=sum(a);
plot(t,ua);
axis([0800100]);
title('
汽车2档原地起步换挡加速时间曲线'
时间t(s)'
速度ua(km/h)’);
2。
7matlab程序:
%绘制汽车功率平衡图。
%汽车参数。
000.793];
367;
f=0.013;
hg=0.9;
798;
G=m*g;
n=600:
10:
Tq=-19.313+295。
44*(n/1000).^2+40.874*(n/1000).^3-3。
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;
P1=Tq。
*ig
(1)*i0。
*ua1。
/(3600*r);
%各档位下发动机的输出功率。
P2=Tq.*ig
(2)*i0.*ua2。
/(3600*r);
P3=Tq.*ig(3)*i0。
*ua3./(3600*r);
P4=Tq.*ig(4)*i0.*ua4./(3600*r);
P5=Tq.*ig(5)*i0.*ua5。
ua=0:
0.2:
120;
Ff=G*f;
Fw=CdA*ua.^2/21.15;
Pf=Ff*ua/3600;
%克服滚动阻力的功率
Pw=Fw.*ua/3600;
%克服空气阻力的功率。
P6=(Pf+Pw)./y;
plot(ua1,P1,ua2,P2,ua3,P3,ua4,P4,ua5,P5,ua,P6);
%绘制汽车功率平衡图。
axis([01200100]);
title(’汽车功率平衡图’);
ylabel(’Pe(kw)’);
legend(’Pe1'
Pe2'
’Pe3'
Pe4’,'
Pe5’,'
P(Ff+Fw)'
%最高档与次高档等速百公里油耗曲线图
1:
%货车相关参数.
ig=[5.562.7691.6441.000。
793];
CdA=2.77;
i0=5.83;
a=1.947;
n0=[8151207161420122603300634033804];
B00=[1326.81354。
71284。
41122.91141.01051。
21233.91129.7];
B10=[-416.46—303。
98-189。
75—121.59—98。
893—73。
714-84。
478-45.291];
B20=[72。
37936。
65714。
5247.00354。
47632.85932。
97880。
71113];
B30=[—5.8629-2.0553-0。
51184-0。
18517—0.091077—0.05138—0。
047449-0.00075215];
B40=[0。
177680。
0430720。
00681640.00185550。
000689060.000350320.00028230-0。
000038568];
B0=spline(n0,B00,n);
%插值拟合系数曲线。
B1=spline(n0,B10,n);
B2=spline(n0,B20,n);
B3=spline(n0,B30,n);
B4=spline(n0,B40,n);
377*r*n/ig(4)/i0;
%4、5档车速。
Fz4=Ff+CdA*(ua4。
^2)/21.15;
%4、5档空气阻力。
Fz5=Ff+CdA*(ua5.^2)/21.15;
Pe4=Fz4。
*ua4。
/(y*3.6*1000);
%4、5档克服行驶阻力所需功率.
Pe5=Fz5。
*ua5。
/(y*3。
6*1000);
fori=1:
3401%计算燃油消耗率。
b4(i)=B0(i)+B1(i)*Pe4(i)+B2(i)*Pe4(i).^2+B3(i)*Pe4(i).^3+B4(i)*Pe4(i).^4;
b5(i)=B0(i)+B1(i)*Pe5(i)+B2(i)*Pe5(i)。
^2+B3(i)*Pe5(i)。
^3+B4(i)*Pe5(i).^4;
pg=7.0;
Q4=Pe4。
*b4./(1。
02.*ua4.*pg);
%计算百公里油耗。
Q5=Pe5.*b5./(1。
02。
*ua5.*pg);
plot(ua4,Q4,ua5,Q5);
%绘制曲线。
axis([01001030]);
title(’最高档与次高档等速百公里油耗曲线'
xlabel(’ua(km/h)’);
百公里油耗(L/100km)'
legend(’4'
5’)