最新汽车理论13和27matlab编程答案.docx
《最新汽车理论13和27matlab编程答案.docx》由会员分享，可在线阅读，更多相关《最新汽车理论13和27matlab编程答案.docx（14页珍藏版）》请在冰豆网上搜索。
最新汽车理论13和27matlab编程答案
孙野20081268
1.3
（1）绘制汽车驱动力与行驶阻力平衡图
选用5挡变速器进行整车性能计算
发动机转速与汽车行驶速度之间的关系：
发动机使用外特性的Tq-n曲线拟合公式：
汽车驱动力：
行驶阻力：
用matlab编写程序如下：
m1=2000;
m2=1800;
mz=3880;
g=9.81;
r=0.367;
CdA=2.77;
f=0.013;
nT=0.85;
ig=[5.562.7691.6441.000.793];
i0=5.83;
If=0.218;
Iw1=1.798;
Iw2=3.598;
%输入已知参数
Iw=2*Iw1+4*Iw2;
fori=1:
69
n（i）=（i+11）*50;
Ttq（i）=-19.313+295.27*（n（i）/1000）-165.44*（n（i）/1000）^2+40.874*（n（i）/1000）^3-3.8445*（n（i）/1000）^4;
end
%从600~4000rpm油拟合公式计算发动机转距
forj=1:
5
fori=1:
69
Ft（i,j）=Ttq（i）*ig（j）*i0*nT/r;
ua（i,j）=0.377*r*n（i）/（ig（j）*i0）;
Fz（i,j）=CdA*ua（i,j）^2/21.15+mz*g*f;
end
end
%计算各档位下的驱动力和行驶阻力
plot（ua,Ft,ua,Fz）
title（'汽车驱动力与行驶阻力平衡图'）;
xlabel（'ua（km/h）'）;
ylabel（'Ft（N）'）;
%作汽车驱动力与行驶阻力平衡图
得到汽车驱动力与行驶阻力平衡图如下：
（2）
①计算汽车最高车速
Matlab程序如下：
fork=1:
175
n1（k）=3300+k*0.1;
Ttq（k）=-19.313+295.27*（n1（k）/1000）-165.44*（n1（k）/1000）^2+40.874*（n1（k）/1000）^33.8445*（n1（k）/1000）^4;
Ft（k）=Ttq（k）*ig（5）*i0*nT/r;
ua（k）=0.377*r*n1（k）/（ig（5）*i0）;
Fz（k）=CdA*ua（k）^2/21.15+mz*g*f;
E（k）=abs（（Ft（k）-Fz（k）））;
end
%从3300rpm开始以0.1rpm为步长搜索直到找到
fork=1:
175
%驱动力与行驶阻力之差绝对值最小点（Ⅴ档处）
if（E（k）==min（E））
%通过其他检验程序可估计出该点大致位置（3310rpm附近）
disp（'汽车最高车速='）;
disp（ua（k））;
disp（'km/h'）;
end
%输出结果
end
汽车最高车速=
98.7627
km/h
计算得汽车最高车速为98.7627km/h。
②计算最大爬坡度
Matlab程序如下：
forp=1:
150
%从2000rpm开始以0.5为步长搜索找驱动力
n2（p）=2000+p*0.5;
%与行驶阻力之差最大点（Ⅰ档处）
Ttq（p）=-19.313+295.27*（n2（p）/1000）-165.44*（n2（p）/1000）^2+40.874*（n2（p）/1000）^3-3.8445*（n2（p）/1000）^4;
Ft（p）=Ttq（p）*ig
（1）*i0*nT/r;
ua（p）=0.377*r*n2（p）/（ig
（1）*i0）;
Fz（p）=CdA*ua（p）^2/21.15+mz*g*f;
af（p）=asin（（Ft（p）-Fz（p））/（mz*g））;
end
%通过其他检验程序可估计出该点大致位置（2000rpm附近）
forp=1:
150
if（af（p）==max（af））
i=tan（af（p））;
disp（'汽车最大爬坡度='）;
disp（i）;
end
end
汽车最大爬坡度=
0.3518
计算得最大爬坡度为0.3518
（3）
①绘制汽车行驶加速度倒数曲线
Matlab程序如下：
fori=1:
69
%开始计算绘制加速度倒数曲线所需数据
n（i）=（i+11）*50;
Ttq（i）=-19.313+295.27*（n（i）/1000）-165.44*（n（i）/1000）^2+40.874*（n（i）/1000）^3-3.8445*（n（i）/1000）^4;
end
forj=1:
5
fori=1:
69
deta=1+Iw/（mz*r^2）+If*ig（j）^2*i0^2*nT/（mz*r^2）;
ua（i,j）=0.377*r*n（i）/（ig（j）*i0）;
a（i,j）=（Ttq（i）*ig（j）*i0*nT/r-CdA*ua（i,j）^2/21.15-mz*g*f）/（deta*mz）;
if（a（i,j）<=0
a（i,j）=a（i-1,j）;
end
if（a（i,j）>0.05）
%剔除加速度值小于0.05的数据
b1（i,j）=a（i,j）;
u1（i,j）=ua（i,j）;
else
b1（i,j）=a（i-1,j）;
u1（i,j）=ua（i-1,j）;
end
b（i,j）=1/b1（i,j）;
end
end
x1=u1（:
1）;y1=b（:
1）;
x2=u1（:
2）;y2=b（:
2）;
x3=u1（:
3）;y3=b（:
3）;
x4=u1（:
4）;y4=b（:
4）;
x5=u1（:
5）;y5=b（:
5）;
plot（x1,y1,x2,y2,x3,y3,x4,y4,x5,y5）;
title（'加速度倒数曲线'）;
axis（[0120030]）;
xlabel（'ua（km/h）'）;
ylabel（'1/aj'）;
%作加速度倒数曲线
绘制汽车行驶加速度倒数曲线如下图所示：
②计算汽车用2挡起步加速行驶至70km/h的加速时间
Matlab程序如下：
fori=1:
69
%从作出的加速度倒数曲线知,各档之间无交点
A=ua（i,3）-ua（69,2）;
%因此换档时刻按低档最高转速对应的车速计
if（A<1&A>0）
j=i;
end
B=ua（i,4）-ua（69,3）;
if（B<2&B>0）
k=i;
end
if（ua（i,4）<=70）
m=i;
end
end
%以上是计算三个换档时刻
t=ua（1,2）*b（1,2）;
%开始计算加速时间,这里是以二档起步得到的原始时间
forp1=2:
69
t1（p1）=（ua（p1,2）-ua（p1-1,2））*（b（p1,2）+b（p1-1,2））*0.5;
t=t+t1（p1）;
end
%为减小计算误差,采用梯形算法累加,以下均是
forp2=j:
69
t2（p2）=（ua（p2,3）-ua（p2-1,3））*（b（p2,3）+b（p2-1,3））*0.5;
t=t+t2（p2）;
end
forp3=k:
m
t3（p3）=（ua（p3,4）-ua（p3-1,4））*（b（p3,4）+b（p3-1,4））*0.5;
t=t+t3（p3）;
end
t=t+（ua（j,3）-ua（69,2））*b（69,2）+（ua（k,4）-ua（69,3））*b（69,3）+（70-ua（m,4））*b（m,4）;
%加上了换档处未计算的那部分时间
tz=t/3.6;
%换算单位
disp（'加速时间='）;
disp（tz）;
disp（'s'）;
%输出计算结果
t=
100.5903
加速时间=
27.9417
s
计算得加速时间为17.9417s。
2.7
（1）汽车功率平衡图
发动机在各挡下的功率Pe：
汽车的阻力功率对车速的关系：
Matlab程序如下：
m1=2000;
m2=1800;
mz=3880;
g=9.81;
r=0.367;
nT=0.85;
f=0.013;
CdA=2.77;
i0=5.83;
If=0.218;
Iw1=1.798;
Iw2=3.598;
ig=[5.562.7691.6441.000.793];
n1=[8151207161420122603300634033804];
Iw=2*Iw1+4*Iw2;
nd=400;Qid=0.299;%输入已知参数
forj=1:
5
fori=1:
69
n（i）=（i+11）*50;
Ttq（i）=-19.313+295.27*（n（i）/1000）-165.44*（n（i）/1000）^2+40.874*（n（i）/1000）^3-3.8445*（n（i）/1000）^4;
Pe（i）=Ttq（i）*n（i）/9550;
ua（i,j）=0.377*r*n（i）/（ig（j）*i0）;
Pz（i,j）=（mz*f*g*ua（i,j）/3600+CdA*ua（i,j）^3/76140）/nT;
end
end%计算各档阻力功率
plot（ua,Pe,ua,Pz）;
xlabel（'ua（km/h）'）;
ylabel（'Pe（kW）'）;
title（'汽车功率平衡图'）;
text（12,65,'Pe1'）;
text（26,65,'Pe2'）;
text（46,65,'Pe3'）;
text（80,65,'Pe4'）;
text（115,65,'Pe5'）;
text（98,78,'（Pf+Pw）/ηT'）;%作功率平衡图
绘制汽车功率平衡图如下：
（2）最高档与次高档的等速百公里油耗曲线
先确定最高档次和次高档次的发动机转速的范围，然后利用，求出对应档位的车速。
由于汽车是等速行驶，因此发动机发出的功率应该比汽车受到的阻力功率折合到曲轴上的功率相等，即。
然后根据不同的和n，用题中给出的拟合公式
求出对应工况的燃油消耗率。
先利用表中的数据，使用插值法，求出每个n值对应的拟合系数：
。
在这里为了保证曲线的光滑性，使用了三次样条插值。
利用求得的各个车速对应下的功率求出对应的耗油量燃油消耗率b。
利用公式，即可求出对应的车速的百公里油耗（L/100km）。
Matlab程序如下：
forj=1:
5
fori=1:
8
Td（i）=-19.313+295.27*（n（i）/1000）-165.44*（n（i）/1000）^2+40.874*（n（i）/1000）^3-3.8445*（n（i）/1000）^4;
Pd（i）=n1（i）*Td（i）/9549;
u（i,j）=0.377*n1（i）*r/（ig（j）*i0）;
end
end%计算在给定转速下的发动机功率
b
（1）=0.17768*Pd
（1）^4-5.8629*Pd
（1）^3+72.379*Pd
（1）^2-416.46*Pd
（1）+1326.8;
b
（2）=0.043072*Pd
（2）^4-2.0553*Pd
（2）^3+36.657*Pd
（2）^2-303.98*Pd
（2）+1354.7;
b（3）=0.0068164*Pd（3）^4-0.51184*Pd（3）^3+14.524*Pd（3）^2-189.75*Pd（3）+1284.4;
b（4）=0.0018555*Pd（4）^4-0.18517*Pd（4）^3+7.0035*Pd（4）^2-121.59*Pd（4）+1122.9;
b（5）=0.
升级会员 