n=linspace(600,4000,100); %均分计算指令 ,600 最低转速 ,4000 最高转速 ,均分为 100 等分
r=0.367;
i0=5.83;
nt=0.85;
G=3880*9.8;
f=0.013;
CDA=2.77;
If=0.218;
Iw1=1.798;
Iw2=3.598;
m=3880;
L=3.2;
a=1.947;
hg=0.9;
ig=[6.09,3.09,1.71,1.00]; %输入已知参数
ua1=0.377*r*n/i0/ig(1);
ua2=0.377*r*n/i0/ig(2);
ua3=0.377*r*n/i0/ig(3);
ua4=0.377*r*n/i0/ig(4); %各转速各挡位下的速度
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; % 从
600~4000rpm 油拟合公式计算发动机转距
Ft1=Tq*i0*ig(1)*nt/r;
Ft2=Tq*i0*ig(2)*nt/r;
Ft3=Tq*i0*ig(3)*nt/r;
Ft4=Tq*i0*ig(4)*nt/r; %从 600~4000rpm 各挡位的驱动力
Ff=G*f;
ua=linspace(0,200,100);
Fw=CDA*ua.*ua/21.15; %空气阻力
plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua,Ff+Fw); %画出各挡位的 Ua-Ft,及 Ua-Ff+Ft
xlabel('ua/ km/h');
ylabel('F/N'); %标注横纵轴
title(' 汽车驱动力 -行驶阻力平衡图 ');%标注图形题目
gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ff+Fw'); %给每根线条添加符号
legend('Ft1','Ft2','Ft3','Ft4','Ff+Fw'); %标注图例
umax=max(ua4);
disp('汽车最高车速 =');
disp(umax);disp('km/h');
imax=tan(asin(max((Ft1-(Ff+Fw))/G))); %最大爬坡度的公式
disp('汽车最大爬坡度 =');
disp(imax);%输出最高车速 ,与最大爬坡度的结果
评论0
最新资源