| 专利名称： | 一种电控柴油发动机单点工况优化方法 | ||
| 专利名称(英文)： | An electrically controlled diesel engine single point mode optimization method for | ||
| 专利号： | CN201410795751.7 | 申请时间： | 20141218 |
| 公开号： | CN104573195A | 公开时间： | 20150429 |
| 申请人： | 东风康明斯发动机有限公司 | ||
| 申请地址： | 441004 湖北省襄樊市高新技术产业开发区 | ||
| 发明人： | 舒咏强; 刘利 | ||
| 分类号： | G06F17/50 | 主分类号： | G06F17/50 |
| 代理机构： | 武汉开元知识产权代理有限公司 42104 | 代理人： | 胡红林 |
| 摘要： | 本发明公开了一种电控柴油发动机单点工况优化方法，该方法包括：通过Cobra生成关键因子设计试验方案，再以该方案采集的台架试验数据作为输入建立基于多元分析的回归模型，在指定参数范围条件下对单点工况进行优化，从而得到最优参数组合及对应的优化结果。设计试验方案及优化过程均通过自编程序实现，实现了单点工况的自动优化过程，在电控柴油发动机开发过程中，极大的减小开发时间及资源成本的同时，方便快捷的实现了兼顾经济性、排放及其他机械限值的参数组合优化。 | ||
| 摘要(英文)： | The invention discloses an electrically-controlled diesel engine single point mode optimization method, the method comprises : generating Cobra key factor design of the test programme, and the programme collecting bench test data as input a regression model based on multivariate analysis, under the condition of the specified parameter range is optimized to the single-point condition, so as to obtain the best parameter combination and corresponding and optimize the result. Design test programme and the optimization process is realized through the self-programming, the single operating condition of the automatic optimization process, in an electronically-controlled diesel engine in the development process, greatly reducing development time and resources at the same time cost, convenient and rapid of the realization of the balance between economy, emissions and other mechanical limit value of the parameter combination optimization. | ||
1.一种电控柴油发动机单点工况优化方法，其特征在于：通过Cobra 生成关键因子设计试验方案，再以该方案采集的台架试验数据作为输入建 立基于多元分析的回归模型，在指定参数范围条件下对单点工况进行优化， 从而得到最优参数组合及对应的优化结果。
2.根据权利要求1所述的电控柴油发动机单点工况优化方法，其特征 在于：所述关键因子设计试验方案由Cobra生成，包括如下步骤： 1)定义所需的关键因子，包括提前角SOI、轨压PRS、预喷油量PIF、 预喷提前角PIT、后喷油量POF、后喷提前角POT、EGR率或EGR阀开度、 可变截面开度VGP，关键因子个数为ivs-2，ivs为回归模型自变量ivi的个 数，试验样本数N，N大于m，其中m为建立多元回归模型最小数据点数； 2)输入外特性扭矩曲线，其中转速按升序排列； 3)通过指定转速范围及负荷范围定义工况区； 4)指定工况区四顶点的各关键因子取值及幅动范围； 5)通过Cobra识别各关键因子在工况区内幅动范围，在此区域内随机 分布试验点，试验点各关键因子在区域取值范围内非重复性的随机分布。
3.根据权利要求1所述的电控柴油发动机单点工况优化方法，其特征 在于：通过Cobra生成的关键因子设计试验方案，在台架进行数据采集后， 作为建立回归模型源数据，回归模型为基于多元分析的多元回归模型，包 括以下步骤： 1)从源数据中定义回归模型的自变量ivi，包括转速speed、扭矩torque 以及指定的关键因子，自变量个数为ivs；从源数据中定义回归模型的因变 量dv，包括喷油量total_fueling、油耗fuel_rate及其他性能排放参数，因变 量个数为dvs； 2)建立回归模型，其模型表达式为： y(dv1,dv2,dv3,…dvdvs)＝f(iv1,iv2,iv3,…ivivs) 其中自变量有各自取值范围 ivimin<ivi<ivimax 3)N组试验数据为回归模型观测矩阵，模型表示为： y(dv1i,dv2i,dv3i,…dvdvsi)＝f(iv1i,iv2i,iv3i,…ivivsi)，其中0<i<N+1。
4.根据权利要求3所述的电控柴油发动机单点工况优化方法，其特征 在于：建立的多元回归模型，针对此区域内需优化的单点工况，在指定约 束条件下进行寻优求解，指定约束条件包括极大值、极小值及区间值三种 类型，所述用作约束条件的参数包括比油耗BSFC、烟度smoke、氮氧化物 NOx、颗粒fsDPM、爆压PCP、增压器转速TCS、涡前排温TIT、压气机 出口温度COT及其他限值参数。
5.根据权利要求4所述的电控柴油发动机单点工况优化方法，其特征 在于：颗粒fsDPM由AVL415烟度转化而得。
6.根据权利要求4或5所述的电控柴油发动机单点工况优化方法，其 特征在于：结合约束条件后建立的多元回归模型表达式为： 所建立的多元回归模型，通过回归整理得到因变量与自变量关系式， 指向优化的单点工况(kps,kpt)后，在此点进行展开得到优化目标函数为：
根据目标函数计算各因变量dvj，在满足各自变量约束同时使BSFC最 小，由此得到的iv3～ivivs即为此工况点输出的最优参数组合，计算而得的 各因变量dvj即为此工况预测的最优结果。
1.一种电控柴油发动机单点工况优化方法，其特征在于：通过Cobra 生成关键因子设计试验方案，再以该方案采集的台架试验数据作为输入建 立基于多元分析的回归模型，在指定参数范围条件下对单点工况进行优化， 从而得到最优参数组合及对应的优化结果。
2.根据权利要求1所述的电控柴油发动机单点工况优化方法，其特征 在于：所述关键因子设计试验方案由Cobra生成，包括如下步骤： 1)定义所需的关键因子，包括提前角SOI、轨压PRS、预喷油量PIF、 预喷提前角PIT、后喷油量POF、后喷提前角POT、EGR率或EGR阀开度、 可变截面开度VGP，关键因子个数为ivs-2，ivs为回归模型自变量ivi的个 数，试验样本数N，N大于m，其中m为建立多元回归模型最小数据点数； 2)输入外特性扭矩曲线，其中转速按升序排列； 3)通过指定转速范围及负荷范围定义工况区； 4)指定工况区四顶点的各关键因子取值及幅动范围； 5)通过Cobra识别各关键因子在工况区内幅动范围，在此区域内随机 分布试验点，试验点各关键因子在区域取值范围内非重复性的随机分布。
3.根据权利要求1所述的电控柴油发动机单点工况优化方法，其特征 在于：通过Cobra生成的关键因子设计试验方案，在台架进行数据采集后， 作为建立回归模型源数据，回归模型为基于多元分析的多元回归模型，包 括以下步骤： 1)从源数据中定义回归模型的自变量ivi，包括转速speed、扭矩torque 以及指定的关键因子，自变量个数为ivs；从源数据中定义回归模型的因变 量dv，包括喷油量total_fueling、油耗fuel_rate及其他性能排放参数，因变 量个数为dvs； 2)建立回归模型，其模型表达式为： y(dv1,dv2,dv3,…dvdvs)＝f(iv1,iv2,iv3,…ivivs) 其中自变量有各自取值范围 ivimin<ivi<ivimax 3)N组试验数据为回归模型观测矩阵，模型表示为： y(dv1i,dv2i,dv3i,…dvdvsi)＝f(iv1i,iv2i,iv3i,…ivivsi)，其中0<i<N+1。
4.根据权利要求3所述的电控柴油发动机单点工况优化方法，其特征 在于：建立的多元回归模型，针对此区域内需优化的单点工况，在指定约 束条件下进行寻优求解，指定约束条件包括极大值、极小值及区间值三种 类型，所述用作约束条件的参数包括比油耗BSFC、烟度smoke、氮氧化物 NOx、颗粒fsDPM、爆压PCP、增压器转速TCS、涡前排温TIT、压气机 出口温度COT及其他限值参数。
5.根据权利要求4所述的电控柴油发动机单点工况优化方法，其特征 在于：颗粒fsDPM由AVL415烟度转化而得。
6.根据权利要求4或5所述的电控柴油发动机单点工况优化方法，其 特征在于：结合约束条件后建立的多元回归模型表达式为： 所建立的多元回归模型，通过回归整理得到因变量与自变量关系式， 指向优化的单点工况(kps,kpt)后，在此点进行展开得到优化目标函数为：
根据目标函数计算各因变量dvj，在满足各自变量约束同时使BSFC最 小，由此得到的iv3～ivivs即为此工况点输出的最优参数组合，计算而得的 各因变量dvj即为此工况预测的最优结果。
技术领域
本发明涉及柴油机的优化方法，具体地说是一种电控柴油发动机单点 工况优化方法。
背景技术
随着社会的不断发展和进步，“节能环保”概念日益重要，在越来越多 的场合和技术领域中都在强调，发动机的开发过程亦不例外。排放法规的 升级，对电控柴油机燃油经济性优化提出严峻挑战，如何在满足相关排放 法规及可靠性指标等要求的同时，提升燃油经济性是电控柴油机行业一直 需要倾注大量资源和时间的一件大事。
电控柴油机有众多指标要求，包括动力性、经济性、可靠性及排放， 使得发动机开发过程需关注多输入多输出，其工作量之大及寻优困难一直 是开发者头疼的关键环节。传统的开发方法，主要依靠人为经验，对电控 柴油发动机燃烧控制参数进行组合。电控柴油发动机燃烧控制参数包括提 前角、轨压、预喷提前角、预喷油量、后喷提前角、后喷油量、EGR率或 EGR阀开度及VG开度等，因参数较多，依靠人工经验进行组合亦较为困 难，需耗费大量时间及台架资源，且不易找出最佳参数组合方案。
发明内容
本发明目的在于克服上述现有技术的不足而提供一种电控柴油发动机 单点工况优化方法，该方法根据输入的外特性曲线、转速扭矩定义的工况 区以及指定的燃烧控制参数关键因子自动生成设计试验样本的软件Cobra， 于台架采集设计试验样本数据后在Cobra中建立基于多元分析的回归模型， 结合输入的自变量、因变量的约束条件，由软件Cobra自动输出最佳参数 组合方案的单点工况优化方法。
实现本发明方法采用的技术方案是一种电控柴油发动机单点工况优化 方法，该方法包括：通过Cobra生成关键因子设计试验方案，再以该方案 采集的台架试验数据作为输入建立基于多元分析的回归模型，在指定参数 范围条件下对单点工况进行优化，从而得到最优参数组合及对应的优化结 果。
在上述技术方案中，所述关键因子设计试验方案由Cobra生成，包括 如下步骤：
1)定义所需的关键因子，包括提前角SOI、轨压PRS、预喷油量PIF、 预喷提前角PIT、后喷油量POF、后喷提前角POT、EGR率或EGR阀开度、 可变截面开度VGP，关键因子个数为ivs-2，ivs为回归模型自变量ivi的个 数，试验样本数N，N大于m，其中m为建立多元回归模型最小数据点数；
2)输入外特性扭矩曲线，其中转速按升序排列；
3)通过指定转速范围及负荷范围定义工况区；
4)指定工况区四顶点的各关键因子取值及幅动范围；
5)通过Cobra识别各关键因子在工况区内幅动范围，在此区域内随机 分布试验点，试验点各关键因子在区域取值范围内非重复性的随机分布。
在上述技术方案中，通过Cobra生成的关键因子设计试验方案，在台 架进行数据采集后，作为建立回归模型源数据，回归模型为基于多元分析 的多元回归模型，包括以下步骤：
1)从源数据中定义回归模型的自变量ivi，包括转速speed、扭矩torque 以及指定的关键因子，自变量个数为ivs；从源数据中定义回归模型的因变 量dv，包括喷油量total_fueling、油耗fuel_rate及其他性能排放参数，因变 量个数为dvs；
2)建立回归模型，其模型表达式为
y(dv1,dv2,dv3,…dvdvs)＝f(iv1,iv2,iv3,…ivivs)
其中自变量有各自取值范围
ivimin<ivi<ivimax
3)N组试验数据为回归模型观测矩阵，模型表示为
y(dv1i,dv2i,dv3i,…dvdvsi)＝f(iv1i,iv2i,iv3i,…ivivsi)，其中0<i<N+1。
在上述技术方案中，建立的多元回归模型，针对此区域内需优化的单 点工况，在指定约束条件下进行寻优求解，指定约束条件包括极大值、极 小值及区间值三种类型，所述用作约束条件的参数包括比油耗BSFC、烟度 smoke、氮氧化物NOx、颗粒fsDPM、爆压PCP、增压器转速TCS、涡前 排温TIT、压气机出口温度COT及其他限值参数。
在上述技术方案中，颗粒fsDPM由AVL415烟度转化而得。
进一步地，结合约束条件后建立的多元回归模型表达式为：
所建立的多元回归模型，通过回归整理得到因变量与自变量关系式， 指向优化的单点工况(kps,kpt)后，在此点进行展开得到优化目标函数为
根据目标函数计算各因变量dvj，在满足各自变量约束同时使BSFC最 小，由此得到的iv3～ivivs即为此工况点输出的最优参数组合，计算而得的 各因变量dvj即为此工况预测的最优结果。
本发明是根据输入的外特性曲线，以及根据转速宽度及负荷宽度定义 工况区，结合指定的燃烧控制参数关键因子及试验样本数自动生成设计试 验样本，由设计试验样本于台架采集的数据为源数据建立基于多元分析的 回归模型，得到比油耗、机械限值、排放及其他因变量与转速、扭矩及燃 烧控制参数关键因子等自变量的关系式，具体到需优化的单点工况生成目 标函数，结合施加的约束条件进行寻优求解，设计试验方案及优化过程均 通过自编程序实现，实现了单点工况的自动优化过程，在电控柴油发动机 开发过程中，极大的减小开发时间及资源成本的同时，方便快捷的实现了 兼顾经济性、排放及其他机械限值的参数组合优化。
附图说明
图1为本发明电控柴油单点工况优化方法示意图。
图2设计试验方案生成流程图。
图3源数据定义示意图。
图4多元回归模型建立流程图。
图5目标函数寻优求解流程图。
具体实施方式
下面结合附图和具体实施例对本发明技术方案作进一步的说明。
本发明电控柴油机单点工况方法，如图1所示，包括：
通过Cobra设计试验模块生成指定燃烧控制参数关键因子的设计试验 方案，再由该方案于台架采集而得试验数据作为多元回归模型的源数据， 结合指定工况点的约束条件，得到指定工况的燃烧控制参数组合；本发明 设计试验方案针对一个区域，该区域由高低转速、高低负荷4点定义，指 定工况点需在此区域内部；指定工况约束条件包括机械限值约束、比油耗 极值约束以及其他性能排放参数约束。
下面对通过Cobra设计试验模块生成指定燃烧控制参数的关键因子设 计试验方案的方法进行详细说明，如图2所示，包括以下具体步骤：
S101、定义所需的关键因子，即发动机燃烧控制参数，包括提前角SOI、 轨压PRS、预喷油量PIF、预喷提前角PIT、后喷油量POF、后喷提前角POT、 EGR率或EGR阀开度、可变截面开度VGP。关键因子个数为ivs-2，ivs为 回归模型自变量ivi的个数，试验样本数N，N大于m，其中m为建立多 元回归模型最小数据点数。
S102、输入外特性扭矩曲线，其中转速按升序排列。
S103、设置试验样本数N，
S104、指定转速范围speedA～speedB及负荷范围loadC～loadB。
S105、结合输入的外特性曲线生成工况区，于步骤S104范围内随机生 成N组试验工况点(spdi,tqi)。
S106、根据工况区四顶点各关键因子中值及偏差，基于参数iv中值宽 度及高度计算名义中值meaniv。
S107、基于参数iv偏差宽度及高度计算名义偏差devsiv，参数iv幅动 范围为(meaniv－devsiv，meaniv+devsiv)。
S108、于步骤S107范围内生成连续N组值valuei。
S109、输出关键因子设计试验方案样本(spdi,tqi,valuei)。
多元回归模型建立流程包括源数据定义及多元数据回归，其中源数据 定义如图3所示：
根据设计试验方案于台架采集而得的试验数据，作为建立基于多元分 析的多元回归模型的源数据，源数据包括转速speed、扭矩torque、比油耗 BSFC、机械限值(如爆压PCP、增压器转速TCS、涡前排温TIT、压气机 出口温度COT)、喷油量total_fueling以及各燃烧控制参数关键因子。源数 据被定义为自变量及因变量两大类，其中转速speed、扭矩torque以及各燃 烧控制参数关键因子为自变量，其余为因变量。
如图4所示，基于输入的源数据及对源数据因变量及自变量定义，生 成回归模型，具体流程如下：
S201、读入源数据，按列存储为userdata[i][j]，其中userdata[0][j]为参 数名，userdata数组行向量为源数据工况点数pamN。
S202、显示参数名userdata[0][j]，将参数分类定义为因变量及自变量。 自变量包括转速speed、扭矩torque、提前角SOI、轨压PRS、预喷油量PIF、 预喷提前角PIT、后喷油量POF、后喷提前角POT、EGR率或EGR阀开度、 可变截面开度VGP等，因变量包括比油耗BSFC、机械限值(如爆压PCP、 增压器转速TCS、涡前排温TIT、压气机出口温度COT)、喷油量 total_fueling、排放及烟度等。自变量个数为ivs，因变量个数为dvs。
S203、观测模型初始化y(dv1,dv2,dv3,…dvdvs)＝f(iv1,iv2,iv3,…ivivs)。
S204、自变量取值范围识别ivimin<ivi<ivimax，ivimin及ivimax分别 存入数组ivmin[i]及ivmax[j]。
S205、对观测模型进行展开，建立因变量与自变量关系式。
S206、源数据代入，生成观测矩阵
dvji＝f(ivji,ivji,ivji,…ivji)其中，0<j<ivs,0<i<N+1。
S207、对观测矩阵进行线性转化
S208、对线性观测矩阵进行回归分析，最小二乘估计求解各系数
S209、输出多元回归模型，其表达式为：
所建立的多元回归模型，通过回归整理得到因变量与自变量关系式， 指向优化的单点工况(kps,kpt)后，在此点进行展开得到优化目标函数为：
当j＝dvs时，dvj＝BSFC，且BSFC取极小值为必选约束，亦即保证在 其他各约束条件满足时比油耗BSFC最小，由此得到的iv3～ivivs即为此工 况点输出的最优参数组合，计算而得的各因变量dvj即为此工况预测的最优 结果。
目标函数寻优求解流程如图5所示。
对优化目标函数在指定单点工况(kps,kpt)展开，其中iv1＝kps，iv2 ＝kpt，iv3～ivivs为其余各自变量。单点工况以转速/扭矩数据存储于txt格 式文本，需在指定工况空间区内。
定义约束条件，包括比油耗BSFC取极小值，烟度smoke、氮氧化物 NOx、颗粒fsDPM、爆压PCP、增压器转速TCS、涡前排温TIT、压气机 出口温度COT及其他限值参数小于给定限值limit[i][j]，根据实际需要进行 拟定。其中，本实施例中，颗粒fsDPM由AVL415烟度转化而得，比油耗 BSFC取极小值为所建立的多元回归模型寻优求解的必选约束，目的在于使 优化工况点在满足各种强制约束条件时获得最佳燃油经济性。
对自变量ivivs随机取值进行组合，全部组合数为ivn，对ivn组组合代 入目标函数计算，若各限值参数均在限值limit[i][j]范围内，记录此时计算 而出的BSFC及dvj，记录于数组optr[ivn][j+ivs],其中第一列为BSFC。对 optr数组第一列求最小值，并同步输出对应的自变量及因变量值，此时iv3～ ivivs即为各约束条件满足且BSFC最优时对应的关键因子组合，对应的各 因变量dvj即为此工况的预测值。
