汽车理论动力性仿真-MATLAB

这篇具有很好参考价值的文章主要介绍了汽车理论动力性仿真-MATLAB。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

引言

汽车动力性是指汽车在良好路面上直线行驶时,由汽车受到的纵向外力决定的、所能达到的平均行驶速度。汽车动力性评价指标主要有汽车的最高车速、加速时间和最大爬坡度。

一、汽车动力性仿真参数

1.汽车基本参数:

汽车总质量m/kg 滾动阻力系数 f f f 空气阻力系数 C D C_D CD 迎风正面面积A/ m 2 m^2 m2 滚动半径r/ m m m 旋转质量换算系数 δ δ δ 传动效率 η t η_t ηt 主减速器传动比 i 0 i_0 i0 变速器各挡传动比 i g i_g ig
936 0.012 0.3 1.75 0.272 1.03+0.04* i g i_g ig 0.9 4.388 [3.416,1.894,1.280,1.000,0.757]

2.发动机转速与转矩的数据:

最低转速:800 r/min;最高转速:6800r/min。
已知转速与转矩数据如下表:

转速/(r/min) 1000 1500 2000 2500 3000 3500 4000 4500 5000
转矩/N • m 78.6 83.0 85.0 86. 6 87. 1 85.9 84.7 82.5 80. 5

利用上表的转速与转矩数据,在MATLAB中编写发动机转矩与转速关系曲线拟合程序如下:

n=[1000,1500,2000,2500,3000,3500,4000,4500,5000];
T=[78.6,83.0,85.0,86.6,87.1,85.9,84.7,82.5,80.5];
cftool;

程序运行后
汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
在1中利用“X 数据”和“Y 数据”的下拉菜单读取转速数据n和转矩数据T;选择多项式函数,在2中选择拟合次数为3,拟合图窗口显示出拟合曲线图,在3中显示拟合公式的参数p1/p2/p3,根据该参数得到发动机转矩与转速的关系:

Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07

注:如已知该公式,在动力性仿真中直接使用,不必重新拟合。

二、驱动力-行驶阻力平衡图

1.驱动力-行驶阻力平衡图

1.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
f=0.012;%滚动阻力系数
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
%汽车的基本参数
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
ua=0:5:220;
Ff=m*g*f;%滚动阻力
Fw=CD*A*ua.^2/21.25;%空气阻力
Fx=Ff+Fw;%行驶阻力
plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fx);
text(35,4100,'Ft1');
text(55,2400,'Ft2');
text(100,1600,'Ft3');
text(135,1300,'Ft4');
text(200,950,'Ft5');
text(175,1200,'Ft+Fw');
xlabel('车速/(km/h)');
ylabel('汽车驱动力/N,行驶阻力/N');
%驱动力-行驶阻力平衡图
Fw5=CD*A*ua5.^2/21.25;
Fx2=Ff+Fw5;
k=find(abs(Ft5-Fx2)<1.5);
uamax=ua5(k);
fprintf('最高车速Vmax = %.2f km/h\n',uamax);
%根据驱动力-行驶阻力平衡图可知驱动力曲线与行驶阻力曲线有交点,交点所对应的车速就是最高车速。

1.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最高车速 V m a x = 173.49 k m / h V_{max} = 173.49 km/h Vmax=173.49km/h

2.各档加速度曲线图

2.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
f=0.012;%滚动阻力系数
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
Fw1=CD*A*ua1.^2/21.25;
Fw2=CD*A*ua2.^2/21.25;
Fw3=CD*A*ua3.^2/21.25;
Fw4=CD*A*ua4.^2/21.25;
Fw5=CD*A*ua5.^2/21.25;
%各档位的迎风阻力
Ff=m*g*f;
dt1=1.03+0.04*ig(1);
dt2=1.03+0.04*ig(2);
dt3=1.03+0.04*ig(3);
dt4=1.03+0.04*ig(4);
dt5=1.03+0.04*ig(5);
%各档位的旋转质量换算系数
a1=(Ft1-Ff-Fw1)/dt1/m;
a2=(Ft2-Ff-Fw2)/dt2/m;
a3=(Ft3-Ff-Fw3)/dt3/m;
a4=(Ft4-Ff-Fw4)/dt4/m;
a5=(Ft5-Ff-Fw5)/dt5/m;
%各档位的加速度
plot(ua1,a1,ua2,a2,ua3,a3,ua4,a4,ua5,a5);
text(35,3.7,'一档');
text(55,2.2,'二档');
text(85,1.4,'三档');
text(130,0.8,'四档');
text(165,0.3,'五档');
ylim([0 4]);
xlabel('车速/(km/h)');
ylabel('加速度/(m/s^2)');
%加速度曲线图
k=max(a1);
amax=k;
fprintf('最大加速度amax = %.2f m/s^2\n',amax);
%由加速度曲线图可知一档时有最大加速度

2.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最大加速度 a m a x = 3.83 m / s 2 a_{max}= 3.83 m/s^2 amax=3.83m/s2

3.各档爬坡度曲线图

3.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
f=0.012;%滚动阻力系数
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
Fw1=CD*A*ua1.^2/21.25;
Fw2=CD*A*ua2.^2/21.25;
Fw3=CD*A*ua3.^2/21.25;
Fw4=CD*A*ua4.^2/21.25;
Fw5=CD*A*ua5.^2/21.25;
%各档位的迎风阻力
Ff=m*g*f;
i1=asin((Ft1-Ff-Fw1)/m/g);
i2=asin((Ft2-Ff-Fw2)/m/g);
i3=asin((Ft3-Ff-Fw3)/m/g);
i4=asin((Ft4-Ff-Fw4)/m/g);
i5=asin((Ft5-Ff-Fw5)/m/g);
%各档位的爬坡度
plot(ua1,i1,ua2,i2,ua3,i3,ua4,i4,ua5,i5);
text(35,0.45,'一档');
text(55,0.25,'二档');
text(90,0.16,'三档');
text(130,0.09,'四档');
text(165,0.03,'五档');
ylim([0 0.5]);
xlabel('车速/(km/h)');
ylabel('爬坡度');
%爬坡度曲线图
k=max(i1);
imax=k;
anglemax=imax*45;
fprintf('最大爬坡度imax = %.2f\n',imax);
fprintf('最大坡度角anglemax = %.2f°\n',anglemax);
%由爬坡度曲线图可知一档时有最大爬坡度和最大坡度角

3.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最大爬坡度 i m a x = 0.47 i_{max} = 0.47 imax=0.47
最大坡度角 a n g l e m a x = 21.33 ° angle_{max}= 21.33° anglemax=21.33°

三、动力特性图

1.动力特性图

1.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
ua=0:5:220;
Fw1=CD*A*ua1.^2/21.25;
Fw2=CD*A*ua2.^2/21.25;
Fw3=CD*A*ua3.^2/21.25;
Fw4=CD*A*ua4.^2/21.25;
Fw5=CD*A*ua5.^2/21.25;
%各档位的迎风阻力
D1=(Ft1-Fw1)/m/g;
D2=(Ft2-Fw2)/m/g;
D3=(Ft3-Fw3)/m/g;
D4=(Ft4-Fw4)/m/g;
D5=(Ft5-Fw5)/m/g;
%各档位的动力因数
c=1.2;%良好沥青路面
f0=0.008;%范围:0.0072~0.0120以上
f1=0.00025;%范围:0.00025~0.00280
f4=0.002;%范围:0.00065~0.002以上
f=c*(f0+f1*(ua/100)+f4*(ua/100).^4);%良好道路上轿车轮胎的滚动阻力系数
plot(ua1,D1,ua2,D2,ua3,D3,ua4,D4,ua5,D5,ua,f);
text(35,0.45,'一档');
text(68,0.25,'二档');
text(104,0.15,'三档');
text(144,0.09,'四档');
text(180,0.02,'五档');
text(60,0.03,'f');
ylim([0 0.5]);
xlabel('速度/(km/h)');
ylabel('动力因数');
%动力特性图
f5=c*(f0+f1*(ua5/100)+f4*(ua5/100).^4);
k=find(abs(D5-f5)<0.0002);
uamax=ua5(k);
fprintf('最高车速Vmax = %.2f km/h\n',uamax);
%由动力特性图可知五档时有最高车速值

1.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最高车速 V m a x = 160.53 k m / h V_{max} = 160.53 km/h Vmax=160.53km/h

2.各档加速度曲线图

2.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
ua=0:5:220;
Fw1=CD*A*ua1.^2/21.25;
Fw2=CD*A*ua2.^2/21.25;
Fw3=CD*A*ua3.^2/21.25;
Fw4=CD*A*ua4.^2/21.25;
Fw5=CD*A*ua5.^2/21.25;
%各档位的迎风阻力
D1=(Ft1-Fw1)/m/g;
D2=(Ft2-Fw2)/m/g;
D3=(Ft3-Fw3)/m/g;
D4=(Ft4-Fw4)/m/g;
D5=(Ft5-Fw5)/m/g;
%各档位的动力因数
c=1.2;%良好沥青路面
f0=0.008;%范围:0.0072~0.0120以上
f1=0.00025;%范围:0.00025~0.00280
f4=0.002;%范围:0.00065~0.002以上
f=c*(f0+f1*(ua/100)+f4*(ua/100).^4);%良好道路上轿车轮胎的滚动阻力系数
f21=c*(f0+f1*(ua1/100)+f4*(ua1/100).^4);
f22=c*(f0+f1*(ua2/100)+f4*(ua2/100).^4);
f23=c*(f0+f1*(ua3/100)+f4*(ua3/100).^4);
f24=c*(f0+f1*(ua4/100)+f4*(ua4/100).^4);
f25=c*(f0+f1*(ua5/100)+f4*(ua5/100).^4);
%各档位的滚动阻力系数
dt1=1.03+0.04*ig(1);
dt2=1.03+0.04*ig(2);
dt3=1.03+0.04*ig(3);
dt4=1.03+0.04*ig(4);
dt5=1.03+0.04*ig(5);
%各档位的旋转质量换算系数
a1=g/dt1*(D1-f21);
a2=g/dt2*(D2-f22);
a3=g/dt3*(D3-f23);
a4=g/dt4*(D4-f24);
a5=g/dt5*(D5-f25);
%各档位的加速度
plot(ua1,a1,ua2,a2,ua3,a3,ua4,a4,ua5,a5);
text(40,3.5,'一档');
text(70,2,'二档');
text(108,1.2,'三档');
text(140,0.6,'四档');
text(160,0.15,'五档');
ylim([0 4]);
xlabel('速度/(km/h)');
ylabel('加速度/(m/s^2)');
%加速度曲线图
k=max(D1-f21);
amax=g/dt1*k;
fprintf('最大加速度amax = %.2f m/s^2\n',amax);
%由加速度曲线图可知一档时有最大加速度

2.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最大加速度 a m a x = 3.85 m / s 2 a_{max} = 3.85 m/s^2 amax=3.85m/s2

3.各档爬坡度曲线图

3.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
f=0.012;%滚动阻力系数
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
ua=0:5:220;
Fw1=CD*A*ua1.^2/21.25;
Fw2=CD*A*ua2.^2/21.25;
Fw3=CD*A*ua3.^2/21.25;
Fw4=CD*A*ua4.^2/21.25;
Fw5=CD*A*ua5.^2/21.25;
%各档位的迎风阻力
D1=(Ft1-Fw1)/m/g;
D2=(Ft2-Fw2)/m/g;
D3=(Ft3-Fw3)/m/g;
D4=(Ft4-Fw4)/m/g;
D5=(Ft5-Fw5)/m/g;
%各档位的动力因数
i1=D1-f;
i2=D2-f;
i3=D3-f;
i4=D4-f;
i5=D5-f;
%各档位的爬坡度
plot(ua1,i1,ua2,i2,ua3,i3,ua4,i4,ua5,i5);
text(38,0.42,'一档');
text(70,0.23,'二档');
text(106,0.14,'三档');
text(136,0.08,'四档');
text(165,0.03,'五档');
ylim([0 0.5]);
xlabel('速度/(km/h)');
ylabel('爬坡度');
%爬坡度曲线图
k=max(D1);
amax=asin((k-f*sqrt(1-k.^2+f.^2))/(1+f.^2));
imax=tan(amax);
anglemax=imax*45;
fprintf('最大爬坡度imax = %.2f\n',imax);
fprintf('最大坡度角anglemax = %.2f°\n',anglemax);
%由爬坡度曲线图可知一档时有最大爬坡度和最大坡度角

3.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最大爬坡度 i m a x = 0.52 i_{max} = 0.52 imax=0.52
最大坡度角 a n g l e m a x = 23.18 ° angle_{max}= 23.18° anglemax=23.18°

四、功率平衡图

1.功率平衡图

1.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
f=0.012;%滚动阻力系数
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
Pe1=Ft1.*ua1/3600;
Pe2=Ft2.*ua2/3600;
Pe3=Ft3.*ua3/3600;
Pe4=Ft4.*ua4/3600;
Pe5=Ft5.*ua5/3600;
%各档位的发动机功率
ua=0:5:220;
Pf=m*g*f*ua/3600;
Pw=CD*A*ua.^3/76140;
Pz=(Pf+Pw)/nt;
%阻力功率
plot(ua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua5,Pe5,ua,Pz);
text(20,42,'一档');
text(50,42,'二档');
text(80,42,'三档');
text(110,42,'四档');
text(185,42,'五档');
text(120,15,'(Pf+Pw)/ηt');
ylim([0 60]);
xlabel('速度/(km/h)');
ylabel('功率/(kw)');
%功率平衡图
Pf5=m*g*f*ua5/3600;
Pw5=CD*A*ua5.^3/76140;
Pz5=(Pf5+Pw5)/nt;
k=find((Pe5-Pz5)>0&(Pe5-Pz5)<0.1);
uamax=ua5(k);
fprintf('最高车速Vmax = %.2f km/h\n',uamax);
%由功率平衡图可知发动机功率曲线与阻力功率曲线相交处对应的车速为最高车速
%由功率平衡图可知发动机功率曲线与阻力功率曲线相交处对应的车速为最高车速

1.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最高车速 V m a x = 164.23 k m / h V_{max} = 164.23 km/h Vmax=164.23km/h

2.各档加速度曲线图

2.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
f=0.012;%滚动阻力系数
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
ua=0:5:220;
Pe1=Ft1.*ua1/3600;
Pe2=Ft2.*ua2/3600;
Pe3=Ft3.*ua3/3600;
Pe4=Ft4.*ua4/3600;
Pe5=Ft5.*ua5/3600;
%各档位的发动机功率
Pf1=m*g*f*ua1/3600;
Pf2=m*g*f*ua2/3600;
Pf3=m*g*f*ua3/3600;
Pf4=m*g*f*ua4/3600;
Pf5=m*g*f*ua5/3600;
%各档位的滚动阻力功率
Pw1=CD*A*ua1.^3/76140;
Pw2=CD*A*ua2.^3/76140;
Pw3=CD*A*ua3.^3/76140;
Pw4=CD*A*ua4.^3/76140;
Pw5=CD*A*ua5.^3/76140;
%各档位的迎风阻力功率
dt1=1.03+0.04*ig(1);
dt2=1.03+0.04*ig(2);
dt3=1.03+0.04*ig(3);
dt4=1.03+0.04*ig(4);
dt5=1.03+0.04*ig(5);
%各档位的旋转质量换算系数
a1=(Pe1-(Pf1+Pw1)/nt)*3600/m/dt1/ua1;
a2=(Pe2-(Pf2+Pw2)/nt)*3600/m/dt2/ua2;
a3=(Pe3-(Pf3+Pw3)/nt)*3600/m/dt3/ua3;
a4=(Pe4-(Pf4+Pw4)/nt)*3600/m/dt4/ua4;
a5=(Pe5-(Pf5+Pw5)/nt)*3600/m/dt5/ua5;
%各档位的加速度
plot(ua1,a1,'.',ua2,a2,'.',ua3,a3,'.',ua4,a4,'.',ua5,a5,'.');
text(50,3.5,'一档');
text(88,1.9,'二档');
text(128,1.1,'三档');
text(164,0.65,'四档');
text(215,0.1,'五档');
ylim([0 4]);
xlabel('速度/(km/h)');
ylabel('加速度/(m/s^2)');
%加速度曲线图
amax=a1;
fprintf('最大加速度amax = %.2f\n',amax);
%由加速度曲线图可知一档时有最大加速度

2.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最大加速度 a m a x = 3.45 m / s 2 a_{max} = 3.45m/s^2 amax=3.45m/s2

3.各档爬坡度曲线图

3.1 MATLAB程序

n=800:10:6800;
Ttq=(2.714e-10)*n.^3-(4.255e-06)*n.^2+0.01756*n+65.07;
m=936;%质量
r=0.272;%滚动半径
nt=0.9;%传动效率
CD=0.3;%空气阻力系数
A=1.75;%迎风面积
f=0.012;%滚动阻力系数
i0=4.388;%主减速器传动比
ig=[3.416,1.894,1.28,1.00,0.757];%变速器各档传动比
g=9.8;%重力加速度
Ft1=Ttq*ig(1)*i0*nt/r;
Ft2=Ttq*ig(2)*i0*nt/r;
Ft3=Ttq*ig(3)*i0*nt/r;
Ft4=Ttq*ig(4)*i0*nt/r;
Ft5=Ttq*ig(5)*i0*nt/r;
%各档位的驱动力
ua1=0.377*r*n/ig(1)/i0;
ua2=0.377*r*n/ig(2)/i0;
ua3=0.377*r*n/ig(3)/i0;
ua4=0.377*r*n/ig(4)/i0;
ua5=0.377*r*n/ig(5)/i0;
%各档位的速度与转速之间的关系
ua=0:5:220;
Pe1=Ft1.*ua1/3600;
Pe2=Ft2.*ua2/3600;
Pe3=Ft3.*ua3/3600;
Pe4=Ft4.*ua4/3600;
Pe5=Ft5.*ua5/3600;
%各档位的发动机功率
Pf1=m*g*f*ua1/3600;
Pf2=m*g*f*ua2/3600;
Pf3=m*g*f*ua3/3600;
Pf4=m*g*f*ua4/3600;
Pf5=m*g*f*ua5/3600;
%各档位的滚动阻力功率
Pw1=CD*A*ua1.^3/76140;
Pw2=CD*A*ua2.^3/76140;
Pw3=CD*A*ua3.^3/76140;
Pw4=CD*A*ua4.^3/76140;
Pw5=CD*A*ua5.^3/76140;
%各档位的迎风阻力功率
i1=(Pe1-(Pf1+Pw1)/nt)*3600/m/g/ua1;
i2=(Pe2-(Pf2+Pw2)/nt)*3600/m/g/ua2;
i3=(Pe3-(Pf3+Pw3)/nt)*3600/m/g/ua3;
i4=(Pe4-(Pf4+Pw4)/nt)*3600/m/g/ua4;
i5=(Pe5-(Pf5+Pw5)/nt)*3600/m/g/ua5;
%各档的爬坡度
plot(ua1,i1,'.',ua2,i2,'.',ua3,i3,'.',ua4,i4,'.',ua5,i5,'.');
text(50,0.41,'一档');
text(87,0.22,'二档');
text(128,0.13,'三档');
text(165,0.07,'四档');
text(215,0.01,'五档');
ylim([0 0.45]);
xlabel('速度/(km/h)');
ylabel('爬坡度');
%爬坡度曲线图
imax=i1;
anglemax=imax*45;
fprintf('最大爬坡度imax = %.2f\n',imax);
fprintf('最大坡度角anglemax = %.2f°\n',anglemax);
%由爬坡度曲线图可知一档时有最大爬坡度和最大坡度角

3.2 运行结果

汽车matlab,汽车理论动力性MATLAB仿真,matlab,汽车
最大爬坡度 i m a x = 0.41 i_{max} = 0.41 imax=0.41
最大坡度角 a n g l e m a x = 18.50 ° angle_{max} = 18.50° anglemax=18.50°

总结

由仿真结果可知,一般情况下,最高车速是在五档;最大加速度(最短加速时间)是在一档;最大爬坡度是在一档。文章来源地址https://www.toymoban.com/news/detail-529420.html

到了这里,关于汽车理论动力性仿真-MATLAB的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 建模分析 | 平面2R机器人(二连杆)运动学与动力学建模(附Matlab仿真)

    🔥附C++/Python/Matlab全套代码🔥课程设计、毕业设计、创新竞赛必备!详细介绍全局规划(图搜索、采样法、智能算法等);局部规划(DWA、APF等);曲线优化(贝塞尔曲线、B样条曲线等)。 🚀详情:图解自动驾驶中的运动规划(Motion Planning),附几十种规

    2024年02月05日
    浏览(59)
  • 基于Matlab的汽车安全应用轨道融合仿真(附源码)

    目录 一、目的 二、轨道到轨道架构 三、定义方案 4.1 在仿真开始时分析跟踪

    2024年02月06日
    浏览(42)
  • 基于证据理论的多源数据融合仿真实验matlab代码

    matlab源码链接见文章末尾 D-S证据理论是一种不确定性推理方法,所处理的数据信息大多是具有不确定性的,该方法能够摆脱对先验概率的依赖,把难以下手分析和处理的完整问题分解成很多易于处理的子问题,利用D-S证据理论将子问题的分析结论综合起来,采用数学推理推理

    2024年02月03日
    浏览(29)
  • 3D-MIMO信道模型的理论概述和matlab仿真

    目录 一、理论基础 二、核心程序 三、测试结果        信道模型的建模,其在不同场景中所对应的参数和分别都是不同的,因此,通过修改信道模型在角度域和延迟域[40]的特定的参数和分布,就可以获得不同的信道模型。        3D-MIMO(Three-Dimensional Multiple-Input Multiple-O

    2024年02月07日
    浏览(38)
  • 基于条纹投影的物体三维形貌测量理论与仿真实验-含Matlab代码

    基于光栅投影的光学三维面形测量具有非接触、精度高、速度快、低成本、全场测量和易于实现等特点,在逆向工程、工业检测、质量控制、虚拟现实、医学诊断、文物修复和物体识别等领域应用广泛。相移轮廓术(PSM)和傅里叶变换轮廓术(FTP)是两种传统方法。相移法计算简单

    2024年02月03日
    浏览(141)
  • 基于Matlab的Robotics Toolbox工具箱的机器人仿真函数介绍(空间位姿表示与动力学)

    随着我们了解到机器人如何建立运动学模型和动力学模型之后,我们可以使用Matlab中的仿真工具箱内来对模型的准确性进行验证,并且可以通过内置的函数进行简单的轨迹规划和可视化观察,本节涉及到的工具箱是MATLAB自带的Robotics Toolbox工具箱。 平移与旋转函数: SE2(x,y,t

    2024年02月07日
    浏览(49)
  • MATLAB仿真Gough-Stewart并联机器人斯图尔特6自由度并联机器人逆运动学仿真 动力学控制pid控制

    MATLAB仿真Gough-Stewart并联机器人斯图尔特6自由度并联机器人逆运动学仿真 动力学控制pid控制 1.搭建了六自由度Stewart并联机器人simulink simscape仿真模型 2.建立了逆向运动学仿真 输入位置和姿态求解各个杆长 3.运用pid控制器进行动力学跟踪控制 使用MATLAB进行了Gough-Stewart并联机器

    2024年01月16日
    浏览(55)
  • 调频连续波(FMCW)波形设计、真实道路场景仿真及汽车自适应巡航控制信号处理(Matlab代码实现)

            目录 💥1 概述 📚2 运行结果 🎉3 参考文献 👨‍💻4 Matlab代码 以下是关于调频连续波(FMCW)波形设计、真实道路场景仿真以及汽车自适应巡航控制信号处理的概述,以及Matlab代码实现的示例: 调频连续波(FMCW)波形设计:FMCW波形是一种特殊的雷达波形,通过改

    2024年02月16日
    浏览(38)
  • 自相关算法,协方差算法,后向加窗算法,前向加窗算法以及QR分解法的理论介绍与matlab仿真分析

    目录 1.自相关算法 2.协方差算法 3.后向加窗算法 4.前向加窗算法 5.QR分解法        自相关算法是一种在信号处理中用来描述信号特性的算法,它主要用于估计一个信号的功率谱。对于一个离散信号x[n],其自相关函数定义为: Rxx[n] = E[x[n+m]*x[m]]       其中E[]表示期望。可以看

    2024年04月09日
    浏览(48)
  • MATLAB - 四旋翼飞行器动力学方程

      本例演示了如何使用 Symbolic Math Toolbox™(符号数学工具箱)推导四旋翼飞行器的连续时间非线性模型。具体来说,本例讨论了 getQuadrotorDynamicsAndJacobian 脚本,该脚本可生成四旋翼状态函数及其雅各布函数。这些函数将在使用非线性模型预测控制(模型预测控制工具箱)控制

    2024年01月22日
    浏览(64)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包