电动汽车动力性MATLAB仿真

电动汽车动力性仿真

i=2;%速比范围
r=0.3;%轮胎半径
fai=1.5;%路面附着系数
Tmax=210;%单电机最大转矩
yita=0.98;%传动系统效率
m=1500;%车重KG
Tn=100;
f=0.012;%滚动阻力系数
Cd=0.3;%风阻系数
A=2.1;%迎风面积
Pn=30;%电机额定功率
g=9.8;%重力加速度
%nn=Pn*9550/Tn;
Ff=m*g*f%滚动阻力
for j=1:1901
    v(j)=0.1*j-0.1;
    n(j)=i*v(j)/r/0.377;
    if n(j)<nn
        Ft(j)=Tmax*i*yita/r;
    else 
        Ft(j)=(Tn*9550/n(j))*i*yita/r;
    end
    Fw(j)=Cd*A*(v(j).^2)/21.5;
    F(j)=Fw(j)+Ff;
    if abs(Ft(j)-F(j))<0.5
        vmax=v(j);
       end
    a(j)=(Ft(j)-F(j))/m;
end
for k=1:1901
    va(1)=0;
    s(1)=0;
    na(k)=i*va(k)/r/0.377;
    if na(k)<nn
        Fta(k)=Tmax*i*yita/r;
    else 
        Fta(k)=(Tn*9550/na(k))*i*yita/r;
    end
    Fwa(k)=Cd*A*(va(k).^2)/21.5;
    Fa(k)=Fwa(k)+Ff;
    acc(k)=(Fta(k)-Fa(k))/yita/m;
    va(k+1)=va(k)+acc(k)*0.1*3.6;
    s(k+1)=s(k)+va(k)*0.1;
        if abs(va(k)-100)<0.5
        ta=(k-1)*0.1;
       
         end
   
    end

figure(1)
plot(v,Ft,v,F)%绘制驱动力和行驶阻力的曲线

figure(2)
t=0:1901;
plot(t*0.1,va)%绘制速度时间曲线

fprintf('vmax=%.2fkm/h\n',vmax)  %写出最高车速
fprintf('ta=%.2fs\n',ta)%写出百公里加速时间

另外，也可以用SIMULINK模型搭建出来