扭矩外特性计算最高车速

clear all
clc

%------------------------
%   发动机外特性
%------------------------
n=[1500:500:5500];
T=[78.59 83.04 85.01 86.63 87.09 85.87 84.67 82.50 80.54];
dt=polyfit(n,T,2);
n1=800:100:5500;
t=polyval(dt,n1);
figure(1)
plot(n1,t,n,T,'o'),grid on
title('发动机外特性-扭矩')
xlabel('转速(rpm)')
ylabel('扭矩(N·m)')


%------------------------
%	汽车驱动力计算
%------------------------
m=input('整车质量(kg),m=');
G=m*9.8;
i=[3.416 1.894 1.28 0.914 0.757];%速比
k=1:5;%5个前进档
r=0.272;%车轮滚动半径
i0=4.388;%主减速比
eff=0.9;%传动机械效率

%计算每个档位的最低行驶速度，即怠速下的车速
vmin=0.377*r*800./(i(k)*i0);

%计算每个档位的最高行驶速度（非最高车速），即额定转速下的车速
vmax=0.377*r*5500./(i(k)*i0);

%依次计算5个档位的驱动力
for k=1:5;
    v=vmin(k):vmax(k);
    n=i(k)*i0.*v./r/0.377;
    t=polyval(dt,n);
    Fd=t.*i(k)*i0*eff/r;
    figure(2)
    plot(v,Fd)
    hold on,grid on
end
title('汽车车速-驱动力曲线')
xlabel('车速(km/h)')
ylabel('驱动力(N)')


%------------------------
%   汽车行驶阻力计算
%------------------------
%行驶阻力包括①滚动阻力Ff、②空气阻力Fw、③坡道阻力Fi和④加速阻力Fj
%此处不考虑上坡和加速
f0=0.009;
f1=0.001;
f4=0.0003;
v=vmin(k):vmax(k);
f=f0+(v/100)*f1+(v/100).^4*f4;%滚阻系数
Ff=G*f;
Cd=0.3;%风阻系数
A=2.3;%迎风面积(m^2)
Fw=v.^2*A*Cd/21.15;%风阻
F=Ff+Fw;
plot(v,F,'k','LineWidth',1.5);
title('汽车车速-驱动力&阻力曲线')
ylabel('驱动力&阻力(N)')
hold on,grid on

%------------------------
%	用图解法求解最高车速
%------------------------
for v=150:0.1:160;k=5;
    n=i(k)*i0.*v./r/0.377;
    t=polyval(dt,n);
    Fd=t.*i(k)*i0*eff/r;
    f=f0+(v/100)*f1+(v/100).^4*f4;
    Ff=G*f;
    Fw=v.^2*A*Cd/21.15;
    F=Ff+Fw;
    if abs(Fd-F)<0.2%该值与速度步长有关，共同调整逼近精度
        v_max=v;
        break
    end
end
disp '=====汽车动力性计算结果====='
fprintf('最高车速=%.2f km/h \n',v_max)