汽车驱动力及电机效率绘制

一、汽车驱动力绘制

例1：已知发动机转矩与转速的关系为

汽车变速器五档传动比分别为3.769，2.095，1.433，1.079，0.851，主减速器传动比为4.188，传动系统效率为0.92，车轮滚动半径为0.29，绘制汽车驱动力曲线。

n=800:10:6800;%定义转速范围
Tq=65.07+0.01756*n-(4.225e-6)*n.^2+(2.714e-10)*n.^3;%计算发动机转矩
r=0.29;%汽车变量赋值
nt=0.92;
i0=4.188;
ig=[3.769 2.095 1.433 1.079 0.851];%变速器传动比赋值
Ft1=Tq*ig(1)*i0*nt/r;%计算一档驱动力
Ft2=Tq*ig(2)*i0*nt/r;%计算二档驱动力
Ft3=Tq*ig(3)*i0*nt/r;%计算三档驱动力
Ft4=Tq*ig(4)*i0*nt/r;%计算四档驱动力
Ft5=Tq*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;%计算五档速度
plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5);%绘制
xlabel('车速/(km/h)')
ylabel('汽车驱动力/N')
text(35,4100,'一档')
text(55,2400,'二档')
text(100,1600,'三档')
text(140,1200,'四档')
text(200,1000,'五档')
print(gcf,'-r600','-djpeg','图1-29.jpg')

二、电机效率MAP图绘制

例2：某电机转速、转矩和效率数据见表1-4，绘制电机二维MAP图和三维MAP图。

表1-4 电机转速、转矩和效率数据

转速/(r/min)	300	300	700	700	700	900	900	900	1500	1500
转矩/(N·m)	138	50	328	138	50	328	138	50	138	50
效率	74.8	84.4	74.9	92.2	92.5	79.4	92.7	92.6	92.2	92.1
转速/(r/min)	3000	3000	4000	4000	5000	5000	6000	6000	7000	7000
转矩/(N·m)	138	50	50	30	50	30	50	30	40	30
效率	95.2	94.1	93.9	91.5	89.8	88.6	88.2	87.9	84.6	84.1

n=[300 300 700 700 700 900 900 900 1500 1500 3000 3000 4000 4000 5000 5000 6000 6000 7000 7000];%转速赋值
T=[138 50 328 138 50 328 138 50 138 50 138 50 50 30 50 30 50 30 40 30];%转矩赋值
nt=[0.748 0.844 0.749 0.922 0.925 0.794 0.927 0.926 0.922 0.921 0.952 0.941 0.939 0.915 0.898 0.886 0.882 0.879 0.846 0.841];%效率赋值
figure(1)%设置图形窗口1
n1=linspace(300,7000,1000);%生成等间隔横坐标
T1=linspace(30,350,1000);%生成等间隔纵坐标
[n1,T1]=meshgrid(n1,T1);%生成二维网格矩阵
nt1=griddata(n,T,nt,n1,T1);%插入分散数据拟合
pcolor(n1,T1,nt1)%绘制伪彩色图
shading interp%设置颜色着色属性
hold on%保存图形
colorbar%生成颜色栏
[c,h]=contourf(n1,T1,nt1,7);%绘制等高线
xlabel('转速/(r/min)')
ylabel('转矩/N·m')
print(gcf,'-r600','-djpeg','图1-30(a).jpg')
%%
figure(2)%设置图形窗口2
n1=linspace(300,7000,50);
T1=linspace(30,350,50);
[n1,T1]=meshgrid(n1,T1);
nt1=griddata(n,T,nt,n1,T1);
mesh(n1,T1,nt1);%绘制三维图
xlabel('转速/(r/min)')
ylabel('转矩/N·m')
zlabel('效率')
print(gcf,'-r600','-djpeg','图1-30(b).jpg')