车辆动力性经济性 matlab计算

一、发动机特性数据输入

t_p_pa_b=[650 170.2 1258.2 232.2
569 149.0 1101.4 232.7
488 127.7 944.6 236.0
406 106.3 785.9 244.1
325 85.1 629.1 257.5
244 63.9 472.3 276.8
163 42.7 315.5 317.4
706 170.0 1366.6 227.1
618 148.8 1196.2 226.7
530 127.6 1025.9 228.8
441 106.2 853.6 234.7
353 85.0 683.3 245.9
265 63.8 512.9 262.8
177 42.6 342.6 297.3
775 170.4 1500.1 223.6
678 149.1 1312.4 225.3
581 127.8 1124.6 229.3
484 106.4 936.8 234.7
388 85.3 751.0 243.4
291 64.0 563.3 258.0
194 42.7 375.5 290.9
830 165.1 1606.6 214.3
726 144.4 1405.3 216.4
623 123.9 1205.9 217.5
519 103.3 1004.6 224.2
415 82.6 803.3 230.5
311 61.9 602.0 241.2
208 41.4 402.6 273.1
850 160.2 1645.3 210.8
744 140.2 1440.1 211.9
638 120.3 1234.9 213.1
531 100.1 1027.8 217.1
425 80.1 822.6 223.4
319 60.1 617.5 236.1
213 40.1 412.3 262.6
850 151.3 1645.3 206.5
744 132.4 1440.1 208.5
638 113.6 1234.9 210.9
531 94.5 1027.8 214.0
425 75.7 822.6 221.9
319 56.8 617.5 231.6
213 37.9 412.3 257.9
850 142.4 1645.3 203.9
744 124.6 1440.1 205.6
638 106.9 1234.9 209.9
531 89.0 1027.8 212.9
425 71.2 822.6 218.8
319 53.4 617.5 229.0
213 35.7 412.3 254.0
850 133.5 1645.3 203.3
744 116.9 1440.1 203.8
638 100.2 1234.9 209.5
531 83.4 1027.8 210.7
425 66.8 822.6 216.2
319 50.1 617.5 226.7
213 33.5 412.3 250.2
850 124.6 1645.3 200.6
744 109.1 1440.1 202.6
638 93.5 1234.9 205.2
531 77.8 1027.8 209.6
425 62.3 822.6 212.3
319 46.8 617.5 220.1
213 31.2 412.3 238.1
850 115.7 1645.3 202.5
744 101.3 1440.1 204.5
638 86.8 1234.9 205.0
531 72.3 1027.8 207.6
425 57.9 822.6 206.2
319 43.4 617.5 216.2
213 29.0 412.3 235.4
850 106.8 1645.3 207.9
744 93.5 1440.1 206.9
638 80.2 1234.9 206.8
531 66.7 1027.8 207.7
425 53.4 822.6 208.7
319 40.1 617.5 215.4
213 26.8 412.3 232.0
785 90.4 1519.5 216.3
689 79.4 1333.7 212.6
589 67.8 1140.1 207.6
491 56.6 950.4 209.5
393 45.3 760.7 213.0
294 33.9 569.1 219.0
196 22.6 379.4 241.2
660 62.2 1277.5 232.4
578 54.5 1118.8 220.8
495 46.6 958.1 217.01
413 38.9 799.4 215.4
330 31.1 638.8 223.08
248 23.4 480.0 233.43
165 15.5 319.4 262.49];
n_temp=[2500
        2500
        2500
        2500
        2500
        2500
        2500
        2300
        2300
        2300
        2300
        2300
        2300
        2300
        2100
        2100
        2100
        2100
        2100
        2100
        2100
        1900
        1900
        1900
        1900
        1900
        1900
        1900
        1800
        1800
        1800
        1800
        1800
        1800
        1800
        1700
        1700
        1700
        1700
        1700
        1700
        1700
        1600
        1600
        1600
        1600
        1600
        1600
        1600
        1500
        1500
        1500
        1500
        1500
        1500
        1500
        1400
        1400
        1400
        1400
        1400
        1400
        1400
        1300
        1300
        1300
        1300
        1300
        1300
        1300
        1200
        1200
        1200
        1200
        1200
        1200
        1200
        1100
        1100
        1100
        1100
        1100
        1100
        1100
         900
         900
         900
         900
         900
         900
         900];
n_t_p_a_b=[n_temp,t_p_pa_b];
T=n_t_p_a_b(:,2);
n=n_t_p_a_b(:,1);
b=n_t_p_a_b(:,5);
%T=T;%!!!!!!!!!台架特性！！！！！！！！！
 Ti=linspace(min(min(T)),max(max(T)),30);%%%%%%%%%%
 ni=linspace(min(min(n)),max(max(n)),30);%%%%%%%%%%
[Ti,ni]=meshgrid(Ti,ni);
bi=griddata(T,n,b,Ti,ni);
subplot(2,2,1)
mesh(ni,Ti,bi);
hold on;
plot3(n,T,b,'ko');
xlabel('n(rpm)');
ylabel('T(Nm)');
zlabel('b(g/kW.h)');
title('万有特性曲线');
%q=interp2(Ti,ni,bi,150,900);
%q1=interp2(Ti,ni,bi,150,900,'cubic');
%q2=interp2(Ti,ni,bi,150,900,'nearest');
v=[200,202,205,210,215,220,225,235,245,265,285];
subplot(2,2,2)
[C,h]=contour(ni,Ti,bi,v,'r*');
clabel(C,h);
xlabel('n(rpm)');
ylabel('T(Nm)');
title('万有特性曲线');
%b_full=max(b')';
%[T_full,ji]=max(T');
n0_T_b=n_t_p_a_b(1:7:90,:);
n0=n0_T_b(:,1);
T_full=n0_T_b(:,2);%!!!!!!!!!台架特性！！！！！！！！！
b_full=n0_T_b(:,5);
Pe=T_full.*n0/9549;
subplot(2,2,3)
plot(n0,[T_full,b_full,Pe]);
xlabel('n(rpm)');
text(1500,450,'T(Nm)');
text(1500,300,'b(g/kW.h)');
text(1500,100,'P(kw)');
title('发动机外特性曲线');

二、车辆动力性经济性计算

ig=[10.338,6.212,4.070,2.56,1.62,1,0.847];%%%%DC7J120T
nn=0.95*0.96;
i0=6.83;
CD=0.8;%%%%%卡车取0.8 客车0.65
G=[60000,9455];%%%额定7850
A=1.92*3.09;
f=0.011;
r=0.502;
qq=1.027+0.000331*i0^2.*ig.^2;
for i=1:length(ig)
    V(:,i)=0.377*r*n0/ig(i)/i0;
    Ft(:,i)=T_full*ig(i)*i0*nn/r*0.95;%*0.95使用特性！！！！！！！！！！！！
    Fw(:,i)=CD*A*V(:,i).^2/21.15;
end
h1=figure;
subplot(2,2,1)
plot(V,n0);
xlabel('V(km/h)');
ylabel('n(rpm)');
title('发动机转速车速曲线');
grid on
G=G*9.8;
D=(Ft-Fw)/G(1);
Ff=G(1)*f;

for i=1:length(ig)
    j(:,i)=9.8*(D(:,i)-f)/qq(i);
end
I=(D.*sqrt(1+f^2-D.^2)-f)./(1-D.^2)*100;
j(find(j<0))=0;
j_1=1./j;
%%%%%%%%%%%%%%%%求加速时间%%%%%%%%%%%%%%%%%！！！！！！！！！！！！！！！
VII=min(V(:,2))+1:1:max(V(:,length(ig)))*0.8;%%%%%2档起步%%%%%%%%%%
VVI=min(V(:,length(ig)-1))+1:1:max(V(:,length(ig)))*0.8;%%%%%直接档%%%%%
for i=1:length(VII)
    TII(i)=j_1(1,2)*V(1,2)/3.6+ t_v_fun(V,j_1,2,min(V(:,2)),VII(i));%%%%%2档起步
end
for i=1:length(VVI)
    TVI(i)= t_v_fun(V,j_1,length(ig)-1,min(V(:,length(ig)-1)),VVI(i));%直接档
end
%%%%%%%%%%%%%%%%求加速时间%%%%%%%%%%%%%%%%%！！！！！！！！！！！！！！！
subplot(2,2,2)
plot(V,Ft,V,Fw+Ff);
xlabel('V(km/h)');
ylabel('Ft,Fw+Ff(N)');
title('驱动力阻力平衡图');
grid on
subplot(2,2,3)
plot(V,D);
xlabel('V(km/h)');
ylabel('D');
title('动力因数图');
grid on
subplot(2,2,4)
plot(V,I);
xlabel('V(km/h)');
ylabel('I（%）');
title('爬坡度图');
grid on
h2=figure;
subplot(2,2,1)
plot(V,j);
xlabel('V(km/h)');
ylabel('j（m/s^2）');
title('加速度图II档起步与直接档');
grid on
subplot(2,2,2)
plot(VII,TII,VVI,TVI);
xlabel('V(km/h)');
ylabel('t（s）');
title('加速时间图');
grid on

Pz=(Fw+Ff).*V/3.6/nn/1000/0.95;%/0.95使用特性变换成万有发动机台架特性！！！！！
%%%%%%%%%%%%%%%%求燃油消耗率%%%%%%%%%%%%%%%%%！！！！！！！！！！！！！！！
n_6=[n0,n0,n0,n0,n0,n0,n0];%%%%%%%7个档
Tz=Pz./n_6*9549;
bii=interp2(Ti,ni,bi,Tz,n_6);
Qii=Pz.*bii./V/1.02/8.33;%7.94-8.13密度8.5
subplot(2,2,3)
plot(V,bii);
xlabel('V(km/h)');
ylabel('b（g/kw.h）');
title('燃油消耗率图');
grid on
subplot(2,2,4)
plot(V,Qii);
xlabel('V(km/h)');
ylabel('Q（L/100km）');
title('燃油消耗量图');
grid on



function t=t_v_fun(V,j_1,xNo,xVmin,xVmax)
[num_V,num_ig]=size(V);
if(xNo<1|xNo>num_ig)
    error('ig_error')
end
if(xVmin<min(V(:,xNo))|xVmin>max(V(:,xNo)))
    error('min_chesu_error')
end
if(xVmax>max(V(:,num_ig))|xVmax<min(V(:,xNo)))
    error('max_chesu_error')
end
t=0;
V_min=xVmin;
i_xNo=xNo;
while(1)
    if(xVmax<=max(V(:,i_xNo)))
        V_max=xVmax;
        if(V_max-V_min>=0.1)
            xi=V_min:0.1:V_max;
            yi=interp1(V(:,i_xNo),j_1(:,i_xNo),xi);
            t=t+trapz(xi,yi);
        end
        t=t/3.6;
    return
    else
        V_max=max(V(:,i_xNo));
        if(V_max-V_min>=0.1)
            xi=V_min:0.1:V_max;
            yi=interp1(V(:,i_xNo),j_1(:,i_xNo),xi);
            t=t+trapz(xi,yi);
        end
        i_xNo=i_xNo+1;
        V_min=V_max;
    end
end