发动机及驱动力曲线绘制.ppt

1、发动机特性曲线matlab绘图详解,杜常清,绘图内容,负荷特性,速度特性,耗油率MAP,转矩MAP,油耗万有特性,效率万有特性,各挡最大驱动力曲线，行驶阻力曲线,一、绘图前的数据准备工作,Clc 清屏 eng_throttle=10:10:100; 设置油门范围及间隔 eng_spd=1250:250:4250;% r/min-1 设置发动机转速范围及间隔，单位：r/min-1,一、绘图前的数据准备工作,输入发动机油耗数据，单位：g/kW*hr,eng_ge_map = 555.29815.8849.33850851852853854855856857858859 328.76376.5238

2、7.86419.53438.6512.3641.73784.29785786787788789 295.34293.35296.48299.25307.71323.72348.02374.9390.05417.27432.7478.45504.14 294.89287.62284.22277.66280.21280.8286.6293.91300.25307.75316.2329.65334.55 292.89285.9283.23279.02275.03274.01278.23277.49273.71273.25279.61283.71288.06 286.47277.96293.98284

3、.19279.67265.29264.44256.23262.42257.08261.55265.85267.26 307.32273.87266.34258.63243.83245.6260.72264.96268.48264.2268.57269.08267.58 304.73296.42286.82285.85275.84280.79260.77264.73263.99261.26265.05267.89269.91 293.03290.28284.22286.97272.78278.02289.94291.42292.14294.79304.01304.79300.54 288.722

4、85.54281.51279.46271.34277.85286.97289.77290.18291.45303.34303.06299.7,一、绘图前的数据准备工作,发动机扭矩数据， 单位：N.m,eng_trq_map=23.813.511.56.7000000000 5642.13730.526.921.115.812.58.37000 80.77469.263.557.150.544.538.834.530.627.623.821 98.796.595.893.288.384.279.674.468.562.557.553.349.9 99.5100.2102.3106.2104.81

5、04.3102.997.795.493.390.38784.5 100.5102.3104.2108.5109.6110.2110.6112.1108.8108.2106.2103.9102.4 101.5104.6107.2110.6114.2116.5117.2118.1119.1119.8117.6116.9115.7 103.2106.4109.2111.5115.9118.1119.5121.2123.4123.1121.5119.6118.4 107.8109.9112.2114.1118.5120.9122.1124.5126.3126.3124.5122.6120.9 109.

6、6112.5113.6116.2119.4121.6123.9125.6127.8127.9125.4123.6121.8;,eng_eff_map=0.147340.100290.0963330.0645000000000 0.248870.21730.210950.195020.186540.159710.12750.104320.076280.068153000 0.277030.278910.275970.273410.265890.252740.23510.218240.209760.196080.189090.171010.16229 0.278450.284470.287870.

7、294670.291990.291380.285480.278380.27250.265860.258750.24820.24456 0.278550.286180.288880.293230.297490.30260.294070.294850.298920.299430.292620.288390.28403 0.278410.288850.288310.29490.300.308410.30940.319320.311780.318260.312820.307760.30614 0.276230.28900.289350.295350.308550.311140.313820.30879

8、0.314750.319680.304640.304070.30577 0.275490.289020.288260.294230.298610.295390.313760.309060.309930.313170.308690.305420.30313 0.274210.288860.287870.294110.295940.294290.287890.283760.289060.287550.269130.268440.27224 0.273380.286540.287640.292770.293150.2910.285110.282360.281960.280730.269720.269

9、970.273;,效率数据,一、绘图前的数据准备工作,二、速度特性曲线的绘制,二、速度特性曲线的绘制,a=size(eng_trq_map,1); figure(color,white) hold on; grid on; for i=1:a plot(eng_spd,eng_trq_map(i,:); end axis(1250,4250,0,140); xlabel(转速n/r.min-1); ylabel(转矩T/N.m); title(速度特性);,三、负荷特性曲线的绘制,三、负荷特性曲线的绘制,figure(color,white) hold on; grid on; b=size(

10、eng_ge_map,2); for i=1:b plot(eng_throttle,eng_ge_map(:,i) end ylabel(耗油率/(g/kW.h); xlabel(节气门开度/%); title(负荷特性);,四、转矩、油耗、效率map二维插值,X,Y=meshgrid(eng_spd,eng_throttle); xi=linspace(1250,4250,31); yi=linspace(10,100,46); X1,Y1=meshgrid(xi,yi); z1=interp2(X,Y,eng_trq_map,X1,Y1,*cubic); z2=interp2(X,Y,e

11、ng_ge_map,X1,Y1,*cubic); z3=interp2(X,Y,eng_eff_map,X1,Y1,*cubic);,五、转矩MAP的绘制,五、转矩MAP的绘制,figure(color,white); title(插值前); surf(X,Y,eng_trq_map) xlabel(转速n/r.min-1); ylabel(节气门开度/%); zlabel(转矩T/N.m); axis(1250,4250,0,100,0,140);,五、转矩MAP的绘制,figure(color,white) title(插值后); surf(X1,Y1,z1); xlabel(转速n/r.

12、min-1); ylabel(节气门开度/%); zlabel(转矩T/N.m); axis(1250,4250,0,100,0,140);,六、油耗MAP的绘制,六、油耗MAP的绘制,figure(color,white) title(插值前); surf(X,Y,eng_ge_map); xlabel(转速n/r.min-1); ylabel(节气门开度/%); zlabel(耗油率/(g/kW.h); axis(1250,4250,0,100,0,1000);,六、油耗MAP的绘制,figure(color,white) title(插值后); surf(X1,Y1,z2); xlabe

13、l(转速n/r.min-1); ylabel(节气门开度/%); zlabel(耗油率/(g/kW.h); axis(1250,4250,0,100,0,1000);,七、油耗率万有特性的绘制,d=size(z1,1); for i=1:d eng_pwr_map(i,:) = z1(i,:).*xi/9550; %(Kw) end figure(color,white) hold on; C,h = contour(X1,z1,z2,20); clabel(C,h); contour(X1,z1,eng_pwr_map,g-); xlabel(n(r/min); ylabel(Ttq(N.m

14、); axis(1250,4250,20,140); grid on; plot(eng_spd,eng_trq_map(end,:),r,LineWidth,2);,八、效率MAP的绘制,figure(color,white) hold on; C,h = contour(X1,z1,z3,20); clabel(C,h); contour(X1,z1,eng_pwr_map,g-); xlabel(n(r/min); ylabel(Ttq(N.m); axis(1250,4250,20,140); grid on; plot(eng_spd,eng_trq_map(end,:),r,LineWidth,2);,九、画各挡最大驱动力曲线，行驶阻力曲线,九、画各挡最大驱动力曲线，行驶阻力曲线,ig=3.31 2.05 1.4 1 0.84; io=3.5; Rr=0.25; nT=0.9; %传动效率 Ga=12000; Cd=0.3; A=2; f=0.013;,数据准备,九、画各挡最大驱动力曲线，行驶阻力曲线,vehspd=(eng_spd*(1./ig)./io*0.377.*Rr; %(Km/h) max_trq