汽车理论2.7

汽车理论2.7

%汽车参数设置
m1=2000;    %装载质量
m2=1800;    %整车整备质量
m3=3880;    %总质量
r=0.367;    %车轮半径
hT=0.85;    %传动系机械效率
f=0.013;    %滚动阻力系数
CdA=2.77;   %空气阻力系数*迎风面积
i0=5.81;    %主减速器传动比
If=0.218;   %飞轮转动惯量
Iw1=1.798;  %二前轮转动惯量
Iw2=3.598;  %四后轮转动惯量
Ig=[5.522;2.769;1.644;1.00;0.793];  %变速器传动比
L=3.2;      %轴距
a=1.947;    %质心至前轴距离（满载）
hg=0.9;     %质心高（满载）
g=9.8;      %重力加速度
pg=7.0;     %汽油的pg

%%绘制汽车功率平衡图
figure(1)
%计算各挡功率
for i=1:5
    ig=Ig(i);
    n=600:0.1:4000;
    ua=0.377*n*r/(i0*ig);
    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*i0*hT/r;
    Pe=Ft.*ua/1000;
    plot(ua,Pe,'LineWidth',3)
    hold on
end
%计算阻力功率
for i=1:5
    ig=Ig(i);
    n=600:0.1:4000;
    ua=0.377*n*r/(i0*ig);
    Pw=CdA*ua.^3/21.15;
    Pf=m3*g*f*ua;
    P=(Pw+Pf)/hT/1000;   
    plot(ua,P,'c','LineWidth',3) 
    hold on
    legend('一档','二档','三挡','四挡','五挡','阻力功率')
end
axis([0 140 0 250]);
xlabel('ua/(Km/h)');
ylabel('Pe/Kw');
title('汽车功率平衡图','FontName','Microsoft YaHei UI');

%%绘制百公里油耗曲线
figure(2)
B0=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7];
B1=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291];
B2=[72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113];
B3=[-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215];
B4=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568];
N=[815 1207 1614 2012 2603 3006 3403 3804];

for j=4:5
    ig=Ig(j);
    i=1:8;
    n=N(i);
    ua=0.337*r*n/(ig*i0);
    Pe=(m3.*g.*f.*ua./3600+CdA.*ua.^3./76140)./hT;
    b=B0(i)+B1(i).*Pe+B2(i).*Pe.^2+B3(i).*Pe.^3+B4(i).*Pe.^4;
    Qs=Pe.*b./(1.02.*ua.*pg);
    p = polyfit(ua,Qs,3);
    Qs_=polyval(p,ua);
    plot(ua,Qs_,'LineWidth',3)
    hold on
    legend('四挡','五挡')
end

axis([0 120 0 35]);
xlabel('ua/(Km/h)');
ylabel('Qs/(L/100km)');
title('百公里油耗曲线','FontName','Microsoft YaHei UI');

运行代码
百公里油耗曲线

汽车功率平衡图