动力计算

本设计的动力性计算的四种方案都采用同一个MATLAB程序，代入不同的发动机和变速器参数后得出不同结果。以下给出动力性计算的基本程序：

%---------基本常数------%

m=18000; g=9.8; G=m*g; r=0.522; yt=0.85; f=0.013; CA=4.9; f=0.018; ig1=x; ig2=x; ig3=x; ig4=x; ig5=x;ig6=x;

%------------------------%

n=n1:1:nmax;

T=xxxxxxxxxxx;

%驱动力%

Ft1=T*ig1*io*yt/r;

Ft2=T*ig2*io*yt/r;

Ft3=T*ig3*io*yt/r;

Ft4=T*ig4*io*yt/r;

io=4.14;

Ft5=T*ig5*io*yt/r;

Ft6=T*ig6*io*yt/r;

%车速u与转速n关系式%

u1=0.377*r*n./ig1/io;

u2=0.377*r*n./ig2/io;

u3=0.377*r*n./ig3/io;

u4=0.377*r*n./ig4/io;

u5=0.377*r*n./ig4/io;

u6=0.377*r*n./ig4/io;

%旋转质量换算系数%

b1=1+0.04 *ig1.^2;

b2=1+0.04 *ig2.^2;

b3=1+0.04 *ig3.^2;

b4=1+0.04 *ig4.^2;

b5=1+0.04 *ig5.^2;

b6=1+0.04 *ig6.^2;

%常遇阻力：滚阻和风阻%

U=0:130;

Ff=G*f;

Ffw=Ff+CA/21.15*U.^2;

%加速度倒数

A1=delta1*m./(Ft1-Ffw1);

A2=delta2*m./(Ft2-Ffw2);

A3=delta3*m./(Ft3-Ffw3);

A4=delta4*m./(Ft4-Ffw4);

%----汽车驱动力-行驶阻力平衡图---%

plot(u1,Ft1,u2,Ft2,u3,Ft3,u4,Ft4,u5,Ft5,u6,Ft6,U,Ffw)

%---汽车行驶加速度倒数曲线图-----%

figure; %新建一张空白图%

plot(u1,A1,u2,A2,u3,A3,u4,A4,u5,A5,u6,A6)

%---------求最高车速Uamax--------%

Uamax=0.377*r*nmax/ig6/io;

%---------求最大爬坡度------------%

syms n %定义n为常量%

T=xxxxxxxxxx;

Ft1=T*ig1*io*yt/r;

u1=0.377*r*n/ig1/io;

Fw=CA*u1^2/21.15;

D=(Ft1-Fw)/G; %I挡动力因素

dD=diff(D,n); %将表达式D对n进行一阶求导

nm=solve(dD,n); %求D最大时n的值

nm=nm(1); %取出实根

Dmax=subs(D,n,nm); %将实根代入D，求得Dmax

a=(Dmax-f*(1+f^2-Dmax^2)^0.5)/(1+f^2);

imax=tan(a) %求最大坡度

%---------求加速时间----------%

syms ig u %定义ig、u为常量

n=ig*io*u/0.377/r;

T=xxxxxxxxxxxxxx;

b=1+0.04+0.04*ig.^2; %旋转质量换算系数

Ft=T*ig*io*yt/r; Ff=G*f;

Ffw=Ff+CA/21.15*u^2;

A=b*m/(Ft-Ffw); %加速度倒数

A1=subs(A,ig,ig1); %将表达式A中的ig用ig2代入

A2=subs(A,ig,ig2);

A3=subs(A,ig,ig3);

A4=subs(A,ig,ig4);

A5=subs(A,ig,ig5);

A6=subs(A,ig,ig6);

u1max=0.377*r*nmax/ig1/io/3.6; %单位 m/s）

u2max=0.377*r*nmax/ig2/io/3.6;

u3max=0.377*r*4000/ig3/io/3.6;

u4max=0.377*r*4000/ig4/io/3.6;

u5max=0.377*r*4000/ig5/io/3.6;

uo=100/3.6;

t1=vpa(int(A1,0,u1max) %Ⅰ挡加速时间

t2=vpa(int(A2,u1max,u2max) %Ⅱ挡加速时间

t3=vpa(int(A3,u2max,u3max) %Ⅲ挡加速时间

t4=vpa(int(A3,u3max,u4max) %Ⅳ挡加速时间

t5=vpa(int(A3,u4max,u5max) %Ⅴ挡加速时间

t6=vpa(int(A3,u5max,u0) %Ⅵ挡加速时间

t=t1+t2+t3+t4+t5+t6 %加速总时间