车辆爬坡能力计算仿真计算

ig=[3.416 1.894 1.28 0.914 0.757];
r=0.272;i0=4.388;eta=0.9;
G=970*9.8;L=2.335;a=1.12;hg=0.71;b=L-a;
k=1:5;
ngk=[800 800 800 800 800 800];
ngm=[5500 5500 5500 5500 5500];
ugk=0.377.*r.*ngk(k)./(ig(k).*i0);
ukm=0.377.*r.*ngm(k)./(ig(k).*i0);

f0=0.009;f1=0.002;f4=0.0003;
Cd=0.3;A=2.3;
for k=1:5
if k==1
u=ugk(k):ukm(k);
n=ig(k)*i0.*u./r./0.377;
t=57.793+1.8363.*(n./100)-2.8832.*(n./1000).^2-8.7594e-11.*n.^3;
Ft=t*ig(k)*i0/r;
f=f0+f1.*(u./100)+f4.*(u./100).^4;
Ff=G*f;
Fw=Cd.*A.*u.^2./21.15;
alpha=asin((Ft-Ff-Fw)/G);
I=tan(alpha);
I1max=max(I);
alphamax=max(alpha)*180/pi;
fprintf('1档的最大爬坡度为I1max=%3.3f\n',I1max)
fprintf('1档的最大爬坡角为alphamax=%3.2f度\n',alphamax)
plot(u,I);
hold on
grid on

elseif k==2
u=ugk(k):ukm(k);
n=ig(k)*i0.*u./r./0.377;
t=57.793+1.8363.*(n./100)-2.8832.*(n./1000).^2-8.7594e-11.*n.^3;
Ft=t*ig(k)*i0/r;
f=f0+f1.*(u./100)+f4.*(u./100).^4;
Ff=G*f;
Fw=Cd.*A.*u.^2./21.15;
alpha=asin((Ft-Ff-Fw)/G);
I=tan(alpha);
I2max=max(I);
plot(u,I);
hold on
grid on

elseif k==3
u=ugk(k):ukm(k);
n=ig(k)*i0.*u./r./0.377;
t=57.793+1.8363.*(n./100)-2.8832.*(n./1000).^2-8.7594e-11.*n.^3;
Ft=t*ig(k)*i0/r;
f=f0+f1.*(u./100)+f4.*(u./100).^4;
Ff=G*f;
Fw=Cd.*A.*u.^2./21.15;
alpha=asin((Ft-Ff-Fw)/G);
I=tan(alpha);
plot(u,I);
hold on
grid on

elseif k==4
u=ugk(k):ukm(k);
n=ig(k)*i0.*u./r./0.377;
t=57.793+1.8363.*(n./100)-2.8832.*(n./1000).^2-8.7594e-11.*n.^3;
Ft=t*ig(k)*i0/r;
f=f0+f1.*(u./100)+f4.*(u./100).^4;
Ff=G*f;
Fw=Cd.*A.*u.^2./21.15;
alpha=asin((Ft-Ff-Fw)/G);
I=tan(alpha);
plot(u,I);
hold on
grid on

elseif k==5
u=ugk(k):ukm(k);
n=ig(k)*i0.*u./r./0.377;
t=57.793+1.8363.*(n./100)-2.8832.*(n./1000).^2-8.7594e-11.*n.^3;
Ft=t*ig(k)*i0/r;
f=f0+f1.*(u./100)+f4.*(u./100).^4;
Ff=G*f;
Fw=Cd.*A.*u.^2./21.15;
alpha=asin((Ft-Ff-Fw)/G);
I=tan(alpha);
plot(u,I);
title('汽车爬坡度图')
xlabel('车速/(km/h)');ylabel('爬坡度/(°)')
hold on
grid on
end
end

q=I1max;
Cfai1=q*L/(b-hg*q);
fai=0.7;
fprintf('前驱汽车匀速上坡时要求的附着系数为Cfai1=%3.2f\n',Cfai1)
if Cfai1>fai
q=b/L/(1/fai+hg/L);
alpha=atan(q)*180/pi;
fprintf('汽车能通过的实际最大坡度为i=%3.2f\n',i)
fprintf('汽车能通过的实际最大坡度角为alpha=%3.2f度\n',alpha)
else
fprintf('1档的最大坡度为I1max=%3.2f\n',I1max)
fprintf('1档的最大坡度角为alphamax=%3.2f度\n',alphamax)
end