CN101706366A - 一种发动机万有特性曲线绘制方法 - Google Patents

一种发动机万有特性曲线绘制方法 Download PDF

Info

Publication number
CN101706366A
CN101706366A CN200910185315A CN200910185315A CN101706366A CN 101706366 A CN101706366 A CN 101706366A CN 200910185315 A CN200910185315 A CN 200910185315A CN 200910185315 A CN200910185315 A CN 200910185315A CN 101706366 A CN101706366 A CN 101706366A
Authority
CN
China
Prior art keywords
matrix
interpolation
rotating speed
fuel consumption
torque
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN200910185315A
Other languages
English (en)
Inventor
李红蕾
宋丰源
刘红领
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chery Automobile Co Ltd
Original Assignee
SAIC Chery Automobile Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by SAIC Chery Automobile Co Ltd filed Critical SAIC Chery Automobile Co Ltd
Priority to CN200910185315A priority Critical patent/CN101706366A/zh
Publication of CN101706366A publication Critical patent/CN101706366A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

本发明公开了一种发动机万有特性曲线绘制方法,本方法利用计算机使用MATLAB软件,对经过试验获得的发动机特性参数进行数据处理,从而绘制发动机特性曲线,这些发动机特性参数包括转速、油门开度量、转矩和燃油消耗量。本发明的技术方案中,利用计算机对发动机台架试验获得的数据绘制成发动机万有特性曲线图,克服了用传统的作图法来绘制发动机万有特性图工作量大,并且由于人为因素,精度较低的缺陷。

Description

一种发动机万有特性曲线绘制方法
技术领域
本发明涉及对汽车发动机进行试验的数据进行数据处理领域,特别涉及一种利用计算机对汽车发动机通过台架试验或者其它试验而获得的发动机特性参数如转速、油门开度量、转矩和燃油消耗量等进行数据处理,最后得出能直观反映出汽车发动机特性的扭矩特性图、燃油特性图、万有特性图、负荷特性图和MAP图等汽车发动机的万有特性曲线。
背景技术
发动机的万有特性曲线就是在一张图上用多参数的特性曲线较全面的表示发动机的性能。通常是以转速n为横坐标,以转矩T或平均有效压力为纵坐标,在图上画出许多等油耗率ge曲线和等功率曲线。这些曲线图主要有:扭矩特性图、燃油特性图、万有特性图、负荷特性图和MAP图。通过万有特性可以方便的查出发动机各种工况下的性能。
在汽车设计与开发工作中,常需要根据发动机台架试验得到的发动机万有特性图与汽车功率平衡图,对汽车燃油经济性进行估算。传统的发动机万有特性图绘制方法主要是将各负荷特性曲线都绘制在同一张纸上,形成曲线簇,然后从曲线簇上把等油耗点逐一投影到万有特性图上,并圆滑地连接成等油耗线。用传统的作图法来绘制发动机万有特性图工作量大,并且由于人为因素,精度较低;随着科学技术的发展,越来越多的人采用MATLAB处理发动机万有特性试验数据,主要采用两种方法实行曲线拟合:一种采用直接曲面拟合法,一种是多项式与插值法相结合的数据拟合法;这两种方法对试验数据的形式要求较高,而且拟合回归方程公式复杂,容易出现图像失真的情况。
发明内容
本发明利用计算机使用MATLAB软件,编写M文件,对发动机试验数据进行大量的矩阵运算和一维插值,方便快捷的绘制发动机万有特性图,使绘制过程流程化,易于操作,并且提高了绘图的精度。
本发明解决其技术问题所采用的技术方案是:一种发动机万有特性曲线绘制方法,本方法利用计算机使用MATLAB软件,对经过试验获得的发动机特性参数进行数据处理,从而绘制发动机特性曲线,所述的发动机特性参数包括转速、油门开度量、转矩和燃油消耗量;包括以下步骤:
A、建立转速的列向量,同一转速下的油门开度量、转矩和燃油消耗量建立行向量;
B、以油门开度量为插值向量,对转矩和燃油消耗量进行插值,得到插值后的转矩矩阵和插值后的燃油消耗量矩阵,然后将转速列向量和油门开度量扩展成扩展的转速矩阵和扩展的油门开度量矩阵;计算机利用扩展的转速矩阵、扩展的油门开度量矩阵、插值后的转矩矩阵和插值后的燃油消耗量矩阵绘制扭矩特性图、燃油特性图和万有特性图;
C、将插值后的燃油消耗量矩阵换算为相应的燃油消耗率矩阵;
D、求出插值后的转矩矩阵的最大值,在所述的最大值范围内,建立插值向量,对同一转速下的燃油消耗率矩阵进行插值,获得第一次燃油消耗率插值矩阵;
E、对所述的转速列向量进行扩充形成扩充的转速向量,以扩充后的转速列向量为插值向量,对第一次燃油消耗率插值矩阵进行插值获得第二次燃油消耗率插值矩阵,并将第二次燃油消耗率插值矩阵换算成燃油消耗量;
F、将所述的插值向量和扩充后的转速列向量扩展成与燃油消耗量和第二次燃油消耗率插值矩阵相同大小的转矩矩阵和转速矩阵,计算机利用转速矩阵、转矩矩阵、燃油消耗量和第二次燃油消耗率插值矩阵绘制负荷特性图和MAP图。
进一步的,上述的一种发动机万有特性曲线绘制方法中:所述的步骤B中的插值过程如下:
将油门开度的最大值和最小值之间设置增量值,油门开度的最大值和最小值之间形成插值向量,以同一转速下的油门开度量作为自变量,转矩和燃油消耗量分别作为因变量,采用分段三次多项式法进行插值,插值结果为插值转矩和插值燃油消耗量。
更进一步的,上述的一种发动机万有特性曲线绘制方法中:所述的油门开度的最大值、最小值和增量值分别为100、0、10。
本发明的技术方案中,利用计算机对发动机台架试验获得的数据进行数据处理,进而绘制成发动机万有特性曲线图,克服了用传统的作图法来绘制发动机万有特性图工作量大,并且由于人为因素,精度较低的缺陷。
以下将结合附图和实施例,对本发明进行较为详细的说明。
附图说明
图1为实施例1绘制发动机万有特性曲线图流程图。
具体实施方式
如图1所示:一种发动机万有特性曲线绘制方法,本方法利用计算机使用MATLAB软件,其中,MATLAB是Mathwork公司于1984年推出的一套高性能的数值计算和可视化软件,它集数值分析、矩阵运算、信号处理和图形显示于一体,可方便地应用于数学计算、数据分析和可视化、科学和工程绘图等方面。本发明运用MATLAB强大的数值处理和绘图功能,提出了一种新的绘制发动机万有特性图的方法。
本方法利用计算机使用MATLAB软件对经过试验获得的发动机特性参数进行运算,从而绘制发动机特性曲线,所述的发动机特性参数包括转速n0、油门开度量k0、转矩T0和燃油消耗量Q0。包括以下步骤:
A、建立转速的列向量n0,同一转速下的油门开度量k0、转矩T0和燃油消耗量Q0建立行向量。本实施例中,将转速列向量n0接升序排列。
B、以油门开度量k0为插值向量,对转矩T0和燃油消耗量Q0进行插值,得到插值后的转矩矩阵T1和插值后的燃油消耗量矩阵Q1,然后将转速列向量n0和油门开度量k0扩展成扩展的转速矩阵n1和扩展的油门开度量矩阵k1;计算机利用扩展的转速矩阵n1、扩展的油门开度量矩阵k1、插值后的转矩矩阵T1和插值后的燃油消耗量矩阵Q1绘制扭矩特性图、燃油特性图和万有特性图.本实施例中,将油门开度的最大值k_max比如10和最小值k_min比如0之间设置增量值dk为1,油门开度的最大值k_max和最小值k_min之间形成插值向量中k=(0、10、20、30、40、50、60、70、80、90、100),当然,为了提高数据处理精度这里k可以有更多阶,比如k可以有100以内的数据,当然为了计算方便可以是其它值只要k_max与k_min之间.以同一转速下的油门开度量作为自变量,转矩T0和燃油消耗量Q0分别作为因变量,采用分段三次多项式法进行插值,插值结果为插值转矩T1和插值燃油消耗量Q1.
C、将插值后的燃油消耗量矩阵Q1换算为相应的燃油消耗率矩阵ge。
D、求出插值后的转矩矩阵T1的最大值T_max,在所述的最大值T_max范围内,建立插值向量Ti,对同一转速下的燃油消耗率矩阵ge进行插值,获得第一次燃油消耗率插值矩阵z。
E、对所述的转速列向量进行扩充形成扩充的转速向量ni,以扩充后的转速列向量ni为插值向量,对第一次燃油消耗率插值矩阵z进行插值获得第二次燃油消耗率插值矩阵Z,并将第二次燃油消耗率插值矩阵Z换算成燃油消耗量Q。
F、将所述的插值向量Ti和扩充后的转速列向量ni扩展成与燃油消耗量Q和第二次燃油消耗率插值矩阵Z相同大小的转矩矩阵Y和转速矩阵X,计算机利用转速矩阵X、转矩矩阵Y、燃油消耗量Q和第二次燃油消耗率插值矩阵Z绘制负荷特性图和MAP图。
下面以某发动机为例,该发动机台架试验的原始数据如表1所示,进行数据处理。
(1)输入发动机的试验数据,并将燃油消耗率换算成燃油消耗量,处理后得到的数据原表如表1所示;建立列向量n,向量的值为[n1;n2;n3;…ni];将k0、T0、Q0转换成行向量,其对应形式如表2所示。
(2)建立插值点向量k,向量的值为[k1,k2,k3,…kj],对同一转速下的T0和Q0分别进行插值,得到矩阵T1和Q1,其形式如表3所示。
(3)将n和k扩展成大小为i×j的矩阵n1和k1,其对应形式如表4所示;对进行绘图,得到发动机扭矩特性曲线图、发动机油耗特性曲线图、发动机万有特性曲面图。
(4)将Q1换算成燃油消耗率ge,求出矩阵T1的最大转矩T_max,并且选择合适的增量值ΔT,建立插值点向量Ti,向量的值为[t1,t2,t3,…ts],对同一转速下的ge进行插值,得到矩阵z,其形式如表5中z所示。
(5)选择合适的增量值Δn,建立向量ni,其值为[N1;N2;N3;…Nr]。
(6)将发动机外特性数据n和T1_max,用ni插值,得到插值结果Ti_max。
(7)依次取z中的每一列,用ni插值,并将超出Ti_max的Ti的元素对应的结果设为NaN,得到矩阵Z,其形式如表5中Z所示。
(8)换算出燃油消耗量Q,并将ni和Ti扩展成大小为r×s的矩阵X、Y,矩阵X、Y、Z、Q的对应形式如表6所示。
(9)对矩阵绘图得到发动机负荷特性图和发动机等油耗线图。
表1:数据源表:
Figure G2009101853157D0000061
表2:数据向量表
Figure G2009101853157D0000062
表3:矩阵T1和Q1
Figure G2009101853157D0000063
表4:矩阵n1和k1
Figure G2009101853157D0000064
表5:矩阵z和Z
Figure G2009101853157D0000065
表6:矩阵X、Y、Q
Figure G2009101853157D0000071

Claims (4)

1.一种发动机万有特性曲线绘制方法,本方法利用计算机使用MATLAB软件,对经过试验获得的发动机特性参数进行数据处理,从而绘制发动机特性曲线,所述的发动机特性参数包括转速(n0)、油门开度量(k0)、转矩(T0)和燃油消耗量(Q0);其特征在于:包括以下步骤:
A、建立转速的列向量(n0),同一转速下的油门开度量(k0)、转矩(T0)和燃油消耗量(Q0)建立行向量;
B、以油门开度量(k0)为插值向量,对转矩(T0)和燃油消耗量(Q0)进行插值,得到插值后的转矩矩阵(T1)和插值后的燃油消耗量矩阵(Q1),然后将转速列向量(n0)和油门开度量(k0)扩展成扩展的转速矩阵(n1)和扩展的油门开度量矩阵(k1);计算机利用扩展的转速矩阵(n1)、扩展的油门开度量矩阵(k1)、插值后的转矩矩阵(T1)和插值后的燃油消耗量矩阵(Q1)绘制扭矩特性图、燃油特性图和万有特性图;
C、将插值后的燃油消耗量矩阵(Q1)换算为相应的燃油消耗率矩阵(ge);
D、求出插值后的转矩矩阵(T1)的最大值(T_max),在所述的最大值(T_max)范围内,建立插值向量(Ti),对同一转速下的燃油消耗率矩阵(ge)进行插值,获得第一次燃油消耗率插值矩阵(z);
E、对所述的转速列向量进行扩充形成扩充的转速向量(ni),以扩充后的转速列向量(ni)为插值向量,对第一次燃油消耗率插值矩阵(z)进行插值获得第二次燃油消耗率插值矩阵(Z),并将第二次燃油消耗率插值矩阵(Z)换算成燃油消耗量(Q);
F、将所述的插值向量(Ti)和扩充后的转速列向量(ni)扩展成与燃油消耗量(Q)和第二次燃油消耗率插值矩阵(Z)相同大小的转矩矩阵(Y)和转速矩阵(X),计算机利用转速矩阵(X)、转矩矩阵(Y)、燃油消耗量(Q)和第二次燃油消耗率插值矩阵(Z)绘制负荷特性图和MAP图。
2.根据权利要求1所述的一种发动机万有特性曲线绘制方法,其特征在于:所述的步骤A中,转速列向量(n0)是按升序排列的。
3.根据权利要求1所述的一种发动机万有特性曲线绘制方法,其特征在于:所述的步骤B中的插值过程如下:
将油门开度的最大值(k_max)和最小值(k_min)之间设置增量值(dk),油门开度的最大值(k_max)和最小值(k_min)之间形成插值向量(k),以同一转速下的油门开度量作为自变量,转矩(T0)和燃油消耗量(Q0)分别作为因变量,采用分段三次多项式法进行插值,插值结果为插值转矩(T1)和插值燃油消耗量(Q1)。
4.根据权利要求3所述的一种发动机万有特性曲线绘制方法,其特征在于:所述的油门开度的最大值(k_max)、最小值(k_min)和增量值(dk)分别为100、0、10。
CN200910185315A 2009-11-06 2009-11-06 一种发动机万有特性曲线绘制方法 Pending CN101706366A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910185315A CN101706366A (zh) 2009-11-06 2009-11-06 一种发动机万有特性曲线绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910185315A CN101706366A (zh) 2009-11-06 2009-11-06 一种发动机万有特性曲线绘制方法

Publications (1)

Publication Number Publication Date
CN101706366A true CN101706366A (zh) 2010-05-12

Family

ID=42376603

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910185315A Pending CN101706366A (zh) 2009-11-06 2009-11-06 一种发动机万有特性曲线绘制方法

Country Status (1)

Country Link
CN (1) CN101706366A (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622722A (zh) * 2012-01-10 2012-08-01 中联重科股份有限公司 发动机特性曲线转化方法及装置
CN103277201A (zh) * 2013-04-24 2013-09-04 三一重机有限公司 电控发动机及其节能方法、节能装置和工程机械
CN104102770A (zh) * 2014-06-28 2014-10-15 中国北方发动机研究所(天津) 一种柴油机实时模型的压气机流量特性拟合方法
CN104196995A (zh) * 2014-08-22 2014-12-10 奇瑞汽车股份有限公司 汽车ecu换档提示方法
CN104454206A (zh) * 2013-09-18 2015-03-25 北汽福田汽车股份有限公司 汽车的控制方法、控制系统及具有该控制系统的汽车
CN105043774A (zh) * 2015-07-14 2015-11-11 中国运载火箭技术研究院 一种采用插值计算舵机功耗的方法
CN105160072A (zh) * 2015-08-06 2015-12-16 中国北方发动机研究所(天津) 柴油机万有特性等油耗曲线绘制用数据的处理方法
CN105486512A (zh) * 2015-11-24 2016-04-13 奇瑞汽车股份有限公司 基于整车环境下的发动机外特性测试方法
US9726279B2 (en) 2015-12-03 2017-08-08 Allison Transmission, Inc. System and method to control the operation of a transmission using engine patterns
US9890851B2 (en) 2015-12-03 2018-02-13 Allison Transmission, Inc. System and method to control the operation of a transmission using engine fuel consumption data
CN108759952A (zh) * 2018-03-09 2018-11-06 程奇 一种分段式基于汽车发动机转速和扭矩的油耗估计方法
CN109115502A (zh) * 2018-08-22 2019-01-01 安徽江淮汽车集团股份有限公司 一种发动机万有特性检测方法及模块
CN110015154A (zh) * 2017-10-25 2019-07-16 上海汽车集团股份有限公司 一种混合动力汽车的怠速充电方法及装置
CN110926822A (zh) * 2019-12-04 2020-03-27 昆山三一动力有限公司 发动机扭矩模型标定方法及发动机扭矩模型标定系统
CN110987450A (zh) * 2019-12-16 2020-04-10 洲际联合超伦科技(北京)有限公司 小型活塞发动机的功率测算方法
CN112417671A (zh) * 2020-11-18 2021-02-26 奇瑞汽车股份有限公司 获得发动机万有特性曲线数据的方法、装置及存储介质
CN112709641A (zh) * 2020-12-15 2021-04-27 中国航发控制系统研究所 二维流量表的反插方法
CN112907694A (zh) * 2019-12-04 2021-06-04 中车时代电动汽车股份有限公司 用于驱动电机系统的效率试验数据处理方法及装置

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622722A (zh) * 2012-01-10 2012-08-01 中联重科股份有限公司 发动机特性曲线转化方法及装置
CN103277201A (zh) * 2013-04-24 2013-09-04 三一重机有限公司 电控发动机及其节能方法、节能装置和工程机械
CN103277201B (zh) * 2013-04-24 2016-05-11 三一重机有限公司 电控发动机及其节能方法、节能装置和工程机械
CN104454206A (zh) * 2013-09-18 2015-03-25 北汽福田汽车股份有限公司 汽车的控制方法、控制系统及具有该控制系统的汽车
CN104102770A (zh) * 2014-06-28 2014-10-15 中国北方发动机研究所(天津) 一种柴油机实时模型的压气机流量特性拟合方法
CN104196995A (zh) * 2014-08-22 2014-12-10 奇瑞汽车股份有限公司 汽车ecu换档提示方法
CN105043774B (zh) * 2015-07-14 2017-12-22 中国运载火箭技术研究院 一种采用插值计算舵机功耗的方法
CN105043774A (zh) * 2015-07-14 2015-11-11 中国运载火箭技术研究院 一种采用插值计算舵机功耗的方法
CN105160072A (zh) * 2015-08-06 2015-12-16 中国北方发动机研究所(天津) 柴油机万有特性等油耗曲线绘制用数据的处理方法
CN105160072B (zh) * 2015-08-06 2018-05-18 中国北方发动机研究所(天津) 柴油机万有特性等油耗曲线绘制用数据的处理方法
CN105486512B (zh) * 2015-11-24 2018-03-16 奇瑞汽车股份有限公司 基于整车环境下的发动机外特性测试方法
CN105486512A (zh) * 2015-11-24 2016-04-13 奇瑞汽车股份有限公司 基于整车环境下的发动机外特性测试方法
US10174832B2 (en) 2015-12-03 2019-01-08 Allison Transmission, Inc. System and method to control the operation of a transmission using engine fuel consumption data
US9726279B2 (en) 2015-12-03 2017-08-08 Allison Transmission, Inc. System and method to control the operation of a transmission using engine patterns
US9890851B2 (en) 2015-12-03 2018-02-13 Allison Transmission, Inc. System and method to control the operation of a transmission using engine fuel consumption data
CN110015154A (zh) * 2017-10-25 2019-07-16 上海汽车集团股份有限公司 一种混合动力汽车的怠速充电方法及装置
CN108759952A (zh) * 2018-03-09 2018-11-06 程奇 一种分段式基于汽车发动机转速和扭矩的油耗估计方法
CN109115502A (zh) * 2018-08-22 2019-01-01 安徽江淮汽车集团股份有限公司 一种发动机万有特性检测方法及模块
CN110926822A (zh) * 2019-12-04 2020-03-27 昆山三一动力有限公司 发动机扭矩模型标定方法及发动机扭矩模型标定系统
CN112907694A (zh) * 2019-12-04 2021-06-04 中车时代电动汽车股份有限公司 用于驱动电机系统的效率试验数据处理方法及装置
CN110926822B (zh) * 2019-12-04 2021-08-06 湖南道依茨动力有限公司 发动机扭矩模型标定方法及发动机扭矩模型标定系统
CN110987450A (zh) * 2019-12-16 2020-04-10 洲际联合超伦科技(北京)有限公司 小型活塞发动机的功率测算方法
CN110987450B (zh) * 2019-12-16 2021-12-07 洲际联合超伦科技(北京)有限公司 小型活塞发动机的功率测算方法
CN112417671A (zh) * 2020-11-18 2021-02-26 奇瑞汽车股份有限公司 获得发动机万有特性曲线数据的方法、装置及存储介质
CN112709641A (zh) * 2020-12-15 2021-04-27 中国航发控制系统研究所 二维流量表的反插方法
CN112709641B (zh) * 2020-12-15 2021-12-07 中国航发控制系统研究所 基于二维流量表确定油针位置的方法

Similar Documents

Publication Publication Date Title
CN101706366A (zh) 一种发动机万有特性曲线绘制方法
CN105844061B (zh) 生成换挡map的方法和装置
CN109409016A (zh) 一种航空发动机压气机非定常流动可视化方法
CN103644036A (zh) 发动机高原动力性能控制方法
Llamas et al. Parameterizing compact and extensible compressor models using orthogonal distance minimization
CN109100148B (zh) 柴油-天然气双燃料商用车燃油经济性计算方法及装置
CN108730044B (zh) 一种发动机的可变气门正时控制方法及装置
CN110765675A (zh) 一种汽车排气系统噪声预测方法
CN109684704B (zh) 一种基于速度密度模型的发动机进气流量在线标定方法
CN113673089B (zh) 一种发动机性能确定方法、装置以及电子设备
CN115288849B (zh) 一种往复式发动机燃烧室磨损状态评价方法及装置
CN113962026B (zh) 航空燃气涡轮过渡态性能相似方法及装置
CN105301652A (zh) 一种三维地震数据固定轴向切片动态判断体绘制方法
CN107719187A (zh) 一种电动汽车电驱动系统扭矩控制提高扭矩精度的方法
CN111159909A (zh) 一种压气机特性数据处理及特性云图绘制方法
CN110941877A (zh) 一种针对纯电动汽车电机工作区域的分析方法
Radu et al. Alternative solutions for supercharging with aggregates of turbocharger type
CN105300702B (zh) 车用柴油机发动机因子的确定方法
CN117780522B (zh) 一种氢气发动机控制方法、装置、车辆及存储介质
CN216899638U (zh) 一种信号采集装置
Wu et al. Research of universal characteristic curve using MATLAB
CN118153183A (zh) 基于响应面分析的航空发动机管路的动刚度设计方法
Tian et al. Application of neural-fuzzy system in twin-turbo hydraulic torque converter's performance testing
CN116738713A (zh) 增压器诊断方法、装置、电子设备及存储介质
Dai et al. Optimization Based on Computer Fluid Dynamics Finite Element Method Automotive Cooling Fan Structure

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20100512