CN104535325A - 一种电控柴油机参数敏感性分析方法 - Google Patents
一种电控柴油机参数敏感性分析方法 Download PDFInfo
- Publication number
- CN104535325A CN104535325A CN201410794541.6A CN201410794541A CN104535325A CN 104535325 A CN104535325 A CN 104535325A CN 201410794541 A CN201410794541 A CN 201410794541A CN 104535325 A CN104535325 A CN 104535325A
- Authority
- CN
- China
- Prior art keywords
- parameter
- operating point
- electronic controlled
- controlled diesel
- data
- 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.)
- Granted
Links
Landscapes
- Combined Controls Of Internal Combustion Engines (AREA)
- Testing Of Engines (AREA)
Abstract
本发明涉及一种电控柴油机参数敏感性分析方法,该方法包括:通过PreDoE计算出电控柴油机运行工况区内关键工况点;再以指定发动机参数对象及因子水平得到参数组合方案,对所述关键工况点指定所述参数组合方案中的参数组合,并结合台架试验对各参数组合进行数据采集,根据采集数据分析各参数敏感性,对于敏感参数,保留至后一阶段的开发方案中,对于非敏感参数可直接剔除。本发明方法适用于各类电控柴油机参数敏感性分析,在发动机性能开发前拟定参数组合方案,使后续发动机性能开发更具针对性,并节约开发时间及成本。
Description
技术领域
本发明涉及电控柴油机参数的分析,具体地指一种电控柴油机参数敏感性分析方法,用于电控柴油发动机的燃烧控制参数敏感性分析。
背景技术
在电控发动机开发过程中,如何通过众多发动机参数组合及优化来实现排放、燃油经济性、动力性及可靠性目标一直是核心课题,此部分工作既是电控发动机开发关键环节,亦是难点环节。在常规开发流程中,通常使用经验判断或大量DoE试验,对工程技术背景及试验成本要求较高,且不易找出最优方案。其缘由在于电控柴油发动机燃烧控制参数包括预喷提前角、预喷油量、后喷提前角、后喷油量及EGR率等参数,各参数取值范围宽窄不一,参数组合数量庞大,若无先期一定敏感性分析判断,开发环节将耗费大量时间,且DoE试验不充足时易遗漏最优方案。
电控柴油机开发环节前进行必要的参数敏感性分析至关重要,可先期确立参数组合方案,缩小参数范围,并基于各参数对发动机燃烧影响敏感性不同,有针对性进行参数调整,从而实现排放、燃油经济性、动力性及可靠性目标。因此,有必要研究一种电控柴油机参数敏感性分析方法,以助于电控柴油机的开发。
发明内容
本发明目的采用的技术方案是一种电控柴油机参数敏感性分析方法,该方法根据输入扭矩曲线、附件功曲线及车辆运行循环工况等数据自动识别和计算关键工况点,在此基础上制定参数因子水平及组合方案,进而输出参数敏感性分析用设计试验方案,最后通过处理设计试验采集数据得出参数敏感性。
实现本发明目的采用的技术方案是一种电控柴油机参数敏感性分析方法,该方法包括:
通过PreDoE计算出电控柴油机运行工况区内关键工况点;
再以指定发动机参数对象及因子水平输出参数组合方案;
对所述关键工况点指定所述参数组合方案中的参数组合,并结合台架试验对各参数组合进行数据采集;
根据采集数据分析各参数敏感性,对于敏感参数,保留至后一阶段的开发方案中;对于非敏感参数可直接剔除。
在上述技术方案中,所述指定的发动机参数对象为电控柴油机燃烧控制参数,包括预喷提前角PIT、预喷油量PIF、后喷提前角POT、后喷油量POF以及EGR率。
进一步地,所述6个电控柴油机燃烧控制参数根据需要自行指定,若参数被指定,则同步对参数因子水平及值进行赋值,继而作为关键工况点燃烧控制参数。
进一步地,各参数的因子水平乘积为kpn,对各关键工况点指定发动机参数时,该工况点对应的试验方案总数为kpn,参数组合方式为完全正交分布。
进一步地,所述的电控柴油机运行工况区包括排放控制区及整车运行工况区。
在上述技术方案中,所述通过PreDoE计算电控柴油机运行工况区内关键工况点包括以下步骤:
1)根据外特性扭矩曲线及附件功曲线插值计算ESC法规中A、B、C转速n及对应的扭矩Tq,所得工况点以二维数组esc[i][j]存储;
2)根据车辆运行循环工况数据,以及划分的转速及负载区识别权重较高的运行工况点,所得工况点以二维数组cycle[i][j]存储;
3)综合二维数组esc[i][j]、cycle[i][j]及指定输入的工况点ikp[i][j]输出关键工况点,所得关键工况点以二维数组kp[i][j]存储。
本发明具有以下优点:
(1)通过PreDoE软件自动计算得到发动机在特定车辆运行循环工况下的运行权重高的工况点;
(2)通过PreDoE软件可自定义燃烧控制参数,方便快捷的定义各工况的参数组合方案,由此迅速输出用于参数敏感性分析用的设计试验方案,有针对性的提高了电控柴油机性能开发的效率。
附图说明
图1为本发明电控柴油机参数敏感性分析方法的流程示意图。
图2为PreDoE计算ESC工况点的流程图。
图3为运行权重高的工况点计算流程图。
图4为参数组合方案输出流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明。
如图1所示,本发明电控柴油机参数敏感性分析方法,包括以下步骤:
S100、通过PreDoE计算出电控柴油机运行工况区内关键工况点,其中,本发明电控柴油机运行工况区包括排放控制区及整车运行工况区。关键工况点的计算过程具体包括以下步骤:
S101、根据外特性扭矩曲线及附件功曲线插值计算ESC法规中A、B、C转速n及对应的扭矩Tq,所得工况点以二维数组esc[i][j]存储。
S102、根据车辆运行循环工况数据,以及划分的转速及负载区识别权重较高的运行工况点,所得工况点以二维数组cycle[i][j]存储。
S103、综合二维数组esc[i][j]、cycle[i][j]及指定输入的工况点ikp[i][j](即用户手动输入工况点,是PreDoE中的自定义工况点,)输出关键工况点,所得关键工况点以二维数组kp[i][j]存储。
其中,上述步骤S101的操作流程如图2所示,包括以下具体步骤:
S1011、在PreDoE中输入外特性扭矩数据,数据按转速从小到大排列,需包含低怠速转速、拐点转速及高怠速转速;
S1012、读入附件功数据,数据按转速从小到大排列,数据按转速从小到大排列,需包含低怠速转速、拐点转速及高怠速转速;
S1013、根据Pe=Tq·n/9550计算外特性功率,除去附件功后得到外特性净功率,并输出最大净功率Pnet;
S1014、计算Pnhi=0.7*Pnet及Pnho=0.5*Pnet,在外特性净功率曲线上反插值得出nhi及nho,其中nhi取大,nho取小;
S1015、根据如下公式计算ESC对应的A、B及C转速,
esca=nho+0.25*(nhi-nho)
escb=nho+0.50*(nhi-nho)
escc=nho+0.75*(nhi-nho);
S1016、在外特性扭矩曲线上插值得到A、B及C转速对应扭矩TQA、TQB及TQC,继而A转速100%负荷点为(A,TQA),B转速100%负荷点为(B,TQB),C转速100%负荷点为(C,TQC),由此可得到各转速对应75%、50%及25%负荷工况点。
S102、根据车辆运行循环工况数据,以及划分的转速及负载区识别权重较高的运行工况点,所得工况点以二维数组cycle[i][j]存储。对于给定车辆运行循环工况数据,将数据点分布于由500~3000rpm及外特性最大扭矩Tqmax组成的外特性扩展区域之上,并将外特性扩展区域划分为若干单区,统计分布于单区之中数据点数,数据点数达到一定值时,此单区被识别为运行权重较高工况区,继而输出权重较高的运行工况点。具体操作如图3所示,包括以下步骤:
S1021、外特性数据传递,以及输入整车运行工况数据;
S1022、读入车辆运行循环工况中转速及扭矩,存入二维数组DCPT[i][j]中,数据点总数为k;
S1023、根据外特性传递数据划分单区,以转速200rpm,扭矩100N.m定义外特性扩展曲线包围下的区域,单区编号为fn[i],单区总数为fnk;
单区详细定义如下:
4-1)从外特性中识别最大扭矩分别为Tqmax,外特性扩展区域由左转速600rpm、右转速3000rpm及最大扭矩Tqmax包围区域组成;
4-2)以转速200rpm,扭矩100N.m定义单区,单区编号为fn[i],区域纵向单元个数为fnk1=12,横向单元个数为fnk2=[Tqmax/100],单区总数fnk=12*fnk2,其中[]为取整。fn[i]其上下左右四个顶点分别为:
左上顶点:转速500+200*(i-j*fnk1)-200扭矩0+j*100
右上顶点:转速500+200*(i-j*fnk1)扭矩0+j*100
左下顶点:转速500+200*(i-j*fnk1)-200扭矩0+j*100-100
右下顶点:转速500+200*(i-j*fnk1)扭矩0+j*100-100
单区中心:转速500+200*(i-j*fnk1)-100扭矩0+j*100-50
其中j=fnk2;
S1024、对数组DCPT[i][j]中转速、扭矩值逐一判断是否处于单区fn[i],统计单区数据组数fnp[i];
S1025、若fnp[i]>fnk/k,则此单区识别为关键单区,单区中心点输出为fnc[i][j];若fnp[i]>fnk/k,则此单区识别为非关键单区,单区中心点输出为fnc[i][j]=[00];
S1026、将统计完之后的中心点排除[00]后输出至二维数组cycle[i][j]。
S200、以指定发动机参数对象及因子水平得出参数组合方案,本发明指定的发动机参数对象为电控柴油机燃烧控制参数,包括预喷提前角PIT、预喷油量PIF、后喷提前角POT、后喷油量POF以及EGR率,上述6个发动机参数可根据需要自行指定,若参数被指定,需同步对参数因子水平(参数因子水平范围为0~3)及值进行赋值,继而作为关键工况点燃烧控制参数。本实施例以指定上述6个发动机参数为例,具体过程如图4所示,包括以下具体步骤:
S201、读入各参数,包括预喷油量因子水平piff、预喷提前角因子水平pitf、后喷油量因子水平poff、后喷提前角因子水平potf、EGR率因子水平egrf。
S202、计算总参数组合kpn=piff*pitf*poff*potf*egr;各参数的因子水平乘积为kpn,对各关键工况点指定发动机参数时,该工况点对应的试验方案总数为kpn,参数组合方式为完全正交分布。
S203、参数组合矩阵pa[m][n]定义,其中m=1~kpn;n=0~5。
S204、以完全正交方式逐一输出参数组合至excel。
S300、对所述关键工况点指定所述参数组合方案中的参数组合,并结合台架试验对各参数组合进行数据采集;根据采集数据分析各参数敏感性,即根据参数组合方案采集而得数据,针对某关键工况点,观察指定因子水平的各参数对排放、燃油经济性等量影响,若影响有限,则判定为非敏感参数,若影响显著则判定为敏感参数。
对于敏感参数,需保留至后一阶段的开发方案中,并逐步细化;非敏感参数可直接剔除,减小后续开发时间及资源消耗。
Claims (8)
1.一种电控柴油机参数敏感性分析方法,其特征在于,包括:
通过PreDoE计算出电控柴油机运行工况区内关键工况点;
再以指定发动机参数对象及因子水平得到参数组合方案,
对所述关键工况点指定所述参数组合方案中的参数组合,并结合台架试验对各参数组合进行数据采集,根据采集数据分析各参数敏感性,对于敏感参数,保留至后一阶段的开发方案中,对于非敏感参数可直接剔除。
2.根据权利要求1所述电控柴油机参数敏感性分析方法,其特征在于所述指定的发动机参数对象为电控柴油机燃烧控制参数,包括预喷提前角PIT、预喷油量PIF、后喷提前角POT、后喷油量POF以及EGR率。
3.根据权利要求2所述电控柴油机参数敏感性分析方法,其特征在于:所述6个电控柴油机燃烧控制参数根据需要自行指定,若参数被指定,则同步对参数因子水平及值进行赋值,继而作为关键工况点燃烧控制参数。
4.根据权利要求3所述电控柴油机参数敏感性分析方法,其特征在于:各参数的因子水平乘积为kpn,对各关键工况点指定发动机参数时,该工况点对应的试验方案总数为kpn,参数组合方式为完全正交分布。
5.根据权利要求4所述电控柴油机参数敏感性分析方法,其特征在于:所述的电控柴油机运行工况区包括排放控制区及整车运行工况区。
6.根据权利要求1~5任一项所述电控柴油机参数敏感性分析方法,其特征在于所述通过PreDoE计算电控柴油机运行工况区内关键工况点包括以下步骤:
1)根据外特性扭矩曲线及附件功曲线插值计算ESC法规中A、B、C转速n及对应的扭矩Tq,所得工况点以二维数组esc[i][j]存储;
2)根据车辆运行循环工况数据,以及划分的转速及负载区识别权重较高的运行工况点,所得工况点以二维数组cycle[i][j]存储;
3)综合二维数组esc[i][j]、cycle[i][j]及指定输入的工况点ikp[i][j]输出关键工况点,所得关键工况点以二维数组kp[i][j]存储。
7.根据权利要求6所述电控柴油机参数敏感性分析方法,其特征在于所述步骤1)包括以下步骤:
1)在PreDoE中输入外特性扭矩数据,数据按转速从小到大排列;
2)读入附件功数据,数据按转速从小到大排列,数据按转速从小到大排列;
3)根据Pe=Tq·n/9550计算外特性功率,除去附件功后得到外特性净功率,并输出最大净功率Pnet;
4)计算Pnhi=0.7*Pnet及Pnho=0.5*Pnet,在外特性净功率曲线上反插值得出nhi及nho,其中nhi取大,nho取小;
5)根据如下公式计算ESC对应的A、B及C转速,
esca=nho+0.25*(nhi-nho)
escb=nho+0.50*(nhi-nho)
escc=nho+0.75*(nhi-nho);
6)在外特性扭矩曲线上插值得到A、B及C转速对应扭矩TQA、TQB及TQC。
8.根据权利要求6所述电控柴油机参数敏感性分析方法,其特征在于所述步骤2)包括:
对于给定车辆运行循环工况数据,将数据点分布于由500~3000rpm及外特性最大扭矩Tqmax组成的外特性扩展区域之上,并将外特性扩展区域划分为若干单区,统计分布于单区之中数据点数,数据点数达到一定值时,此单区被识别为运行权重较高工况区,继而输出权重较高的运行工况点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410794541.6A CN104535325B (zh) | 2014-12-18 | 2014-12-18 | 一种电控柴油机参数敏感性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410794541.6A CN104535325B (zh) | 2014-12-18 | 2014-12-18 | 一种电控柴油机参数敏感性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104535325A true CN104535325A (zh) | 2015-04-22 |
CN104535325B CN104535325B (zh) | 2017-10-20 |
Family
ID=52850899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410794541.6A Active CN104535325B (zh) | 2014-12-18 | 2014-12-18 | 一种电控柴油机参数敏感性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104535325B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111046114A (zh) * | 2019-12-24 | 2020-04-21 | 天津雷沃发动机有限公司 | 一种基于Excel的路谱采集分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0319224D0 (en) * | 2003-08-15 | 2003-09-17 | Omitec Group Ltd | Testing diesel engines |
JP2006244042A (ja) * | 2005-03-02 | 2006-09-14 | Denso Corp | エンジン制御パラメータの適合方法及び適合装置 |
CN103615299A (zh) * | 2013-11-21 | 2014-03-05 | 同济大学 | 设计柴油机后处理系统的方法及系统 |
CN103748340A (zh) * | 2011-08-19 | 2014-04-23 | 特斯托股份公司 | 用于确定内燃机中NOx和/或CO2的特定和/或绝对排放值的方法及测量装置 |
-
2014
- 2014-12-18 CN CN201410794541.6A patent/CN104535325B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0319224D0 (en) * | 2003-08-15 | 2003-09-17 | Omitec Group Ltd | Testing diesel engines |
JP2006244042A (ja) * | 2005-03-02 | 2006-09-14 | Denso Corp | エンジン制御パラメータの適合方法及び適合装置 |
CN103748340A (zh) * | 2011-08-19 | 2014-04-23 | 特斯托股份公司 | 用于确定内燃机中NOx和/或CO2的特定和/或绝对排放值的方法及测量装置 |
CN103615299A (zh) * | 2013-11-21 | 2014-03-05 | 同济大学 | 设计柴油机后处理系统的方法及系统 |
Non-Patent Citations (2)
Title |
---|
安士杰,等: "基于DOE与遗传算法相结合的共轨柴油机喷射参数优化", 《武汉理工大学学报(交通科学与工程版)》 * |
李铁栓,等: "柴油机ESC排放自动测试系统开发", 《车辆与动力技术》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111046114A (zh) * | 2019-12-24 | 2020-04-21 | 天津雷沃发动机有限公司 | 一种基于Excel的路谱采集分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104535325B (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104850629A (zh) | 一种基于改进k-means算法的海量智能用电数据分析方法 | |
CN101047316A (zh) | 确定与电网相关联的参数值的系统、方法和制品 | |
TW200925015A (en) | Mixed type vehicle power system and method of forming multidimentional data of fuel consumption | |
JP2012523620A5 (zh) | ||
CN101047315A (zh) | 控制发电系统运行的系统、方法和制品 | |
CN105958603A (zh) | 一种电池的充电方法和装置 | |
CN106786977B (zh) | 一种电动汽车充电站的充电调度方法 | |
CN103034744A (zh) | 一种低速柴油机气缸盖温度场的模拟方法 | |
CN104763585B (zh) | 基于分布式数据收集的风电机组动态重构方法 | |
CN104391748A (zh) | 一种mapreduce计算过程优化方法 | |
CN106374513A (zh) | 一种基于主从博弈的多微网联络线功率优化方法 | |
CN102541622B (zh) | 一种负载相关的虚拟机放置方法 | |
CN105046593A (zh) | 一种符合低碳能源政策的智能配用电评价方法 | |
Sun et al. | Multi-objective component sizing for a battery-supercapacitor power supply considering the use of a power converter | |
CN113852135A (zh) | 一种虚拟电厂能源调度方法、装置、存储介质及平台 | |
CN104535325A (zh) | 一种电控柴油机参数敏感性分析方法 | |
CN111489009A (zh) | 一种电动汽车充电站运行方式的优化计算方法及装置 | |
CN108959793B (zh) | 一种伞齿轮双频分段循环感应加热数值模拟方法 | |
CN102831230B (zh) | 一种混合动力总成的部件选型方法及系统 | |
CN111539138B (zh) | 基于阶跃函数的结构动力学峰值时域响应灵敏度求解方法 | |
CN104573195B (zh) | 一种电控柴油发动机单点工况优化方法 | |
CN103744296B (zh) | 一种机电控一体化优化设计方法 | |
CN115642655B (zh) | 一种水电厂多机组负荷分配优化方法及系统 | |
CN106096180B (zh) | 一种差异化配置末级叶片模块的低压缸改造方法 | |
CN110569590B (zh) | 基于整车转鼓试验的phev部件工况数据统计与工况构建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |