2. Department of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China
电控单体泵是一种能够满足当前柴油机排放法规和经济性要求的新型燃油喷射系统,可实现较高的喷油压力及良好的燃料雾化[1]。高速电磁阀 (简称电磁阀) 是电控单体泵的核心组件之一,它的响应速度决定了喷油压力的建立与喷射后油压卸载速度等特性,从而会影响到喷油系统的喷油定时、循环喷油量等关键特性[2]。电磁阀较大的响应延迟会引起喷油定时失准和循环喷油量的精度变差,从而导致柴油机排放超标及油耗增大。为进一步提高喷油控制的精确性,需要对影响电磁阀响应的关键特性参数进行优化设计,以减小电磁阀的响应延迟时间。
目前,国内在电磁阀的铁芯材质、驱动电路设计等方面进行了较多的实验研究和优化设计。张奇等对电控柴油机的电磁阀驱动电路进行了分析,通过有限元软件Ansys对电磁阀进行有限元建模,模拟了电磁阀关闭动态响应过程,并对驱动电路进行优化设计,降低了电磁阀闭合响应时间[3]。张廷羽等通过Ansys分析电磁阀的电磁部分,利用AMESim建立了电磁阀整体仿真模型,对影响电磁阀的各个因素进行了计算和分析,并提出了适合电磁阀铁心材质、线圈等优化设计的方案[4]。李铁栓等采用模拟退火算法,结合电磁阀Ansys有限元仿真模型,通过多目标优化平台modeFrontier对高压共轨电磁阀的开启、关闭延迟时间进行优化设计,降低了电磁阀的开启和关闭延迟时间[5]。
以上所述对电磁阀的优化设计皆把电磁阀作为独立系统进行优化设计,本文把电磁阀作为电控单体泵系统中的一部分进行电磁阀延迟响应时间优化,并分析优化后对电控单体泵系统喷射特性的影响。本文在AMESim环境下建立电控单体泵仿真模型,并结合试验数据对模型的准确性进行验证。应用实验设计方法,对电控单体泵电磁阀部分各特性参与电磁阀响应延迟时间的相关性进行深入分析,得到电磁阀响应的关键影响参数。进而利用多目标优化平台modeFRONTIER与AMESim进行联合仿真,对电磁阀响应时间的关键影响参数进行优化,以获得电磁阀与电信号的开启和关闭延迟最小响应时间。
1 电磁阀结构组成及工作原理 1.1 电磁阀工作原理电控单体泵及电磁阀的结构如图 1所示,主要包括电磁阀控制部分和柱塞加压部分:柱塞加压部分包括柱塞、柱塞套和柱塞弹簧;电磁阀控制部分主要包括电磁铁、衔铁、控制阀杆、衔铁复位弹簧、出油堵头等零部件。通电后,电磁铁吸合衔铁,拉动控制阀杆,关闭密封锥面,切断燃油回路,从而在泵腔内建立起燃油喷射所需的高压;断电后,复位弹簧迫使衔铁推动控制阀杆复位,开启密封锥面,卸载高压燃油,停止燃油喷射。该方式实现了对燃油喷射过程的数字控制,改变了传统喷油泵复杂的机械控制方式,通过调节控制阀杆的闭合时间和闭合时刻,可实现对循环喷油量和喷油定时的灵活控制。
1.2 电磁阀响应延迟特性在电磁阀关闭和开启阶段,由于电流作用时间、磁滞现象、阀杆动作延迟、阀杆运动过程等因素导致衔铁相对于控制电流 (信号) 存在滞后所产生的延迟时间。在电磁阀开启阶段由于较高压力的燃油迅速泄压,在燃油流经处会产生局部压降现象,使得燃油会对锥阀产生与开启方向相反的液动力,从而导致开启延迟时间增加。
根据衔铁受电磁力作用过程,可以将电磁阀的运动过程分为3个阶段 (如图 2所示):
1) 关闭延迟阶段A:线圈通电起始时间到衔铁完全吸合时间即电磁阀到达最大升程时间。
2) 关闭阶段C。衔铁被吸合后,与电磁铁间隙较小,电磁力一直大于弹簧力,电磁阀保持完全关闭状态。
3) 开启延迟阶段B。电磁阀断电时间到衔铁复位即电磁阀完全打开时间。
2 仿真模型电控单体泵系统是集电场、磁场、机械运动和流场于一体的复杂系统,各物理场通过各自的控制方程及相互作用变量耦合在一起,每个场的参数都会直接或间接的影响最终的喷射特性。单纯的实验研究难以实现对系统喷射特性的深入分析,需采用数值仿真的方法,以揭示影响喷射特性的本质规律。因此,本文在AMESim环境中建立了数值模型,如图 3所示。
为验证模型的准确性,图 4中 (a) 是凸轮转速为900 r/min,循环喷油量为150 mm3时,在相同控制电流下,喷油规律的实验和仿真结果对比曲线。
图 4(b)是在不同凸轮转速下,满足实验和仿真相同循环喷油量 (某种机型外特性下循环喷油量) 的前提下,泵端压力和嘴端压力的实验和仿真结果。可见仿真模型能够准确的预测泵端压力、嘴端压力和喷油规律,由图中可知系统在喷油过程中的燃油喷射持续期和喷油提前角是一致的,喷油时序也是一致的,而且在任何转速下泵端压力和嘴端压力都能得到很好的预测,最大偏差为7%,因此该模型能够准确预测系统的喷射特性参数[6]。
3 电磁阀响应延迟时间相关性预测电磁阀部分的特性参数会影响电磁阀的响应延迟特性,为筛选关键因素,在保持控制策略及电磁阀固有物性参数如升程、质量和材质不变的前提下,选取了如表 1所示的影响参数对电磁阀响应延迟经行了相关性预测分析,其中表 1中参数基准值是实际测量值,取值范围是在不影响电控单体泵及电磁阀等各部分正常工作所允许的波动范围。
参数 | 基准值 | 取值范围 |
弹簧预紧力/N | 70 | 50~90 |
阻尼/(N·(m·s-1)-1) | 50 | 0~120 |
锥阀半锥角/(°) | 68 | 56~80 |
衔铁残余气隙/mm | 0.12 | 0.08~0.15 |
阀杆直径/mm | 5.5 | 5~6 |
锥阀直径/mm | 7.6 | 7.5~8.2 |
阀杆孔径/mm | 7 | 6.3~7.3 |
电磁阀容积腔/mm3 | 0.41 | 0.3~0.5 |
弹簧刚度/(N·m-1) | 10 000 | 8 000~12 000 |
相关性分析是考察两个变量之间线性关系的一种分析方法,即,当一个变量发生变化时,另一变量随之如何变化,此时需要通过计算相关系数来做深入的定量考察。相关性分析需要从两方面来考察两变量间的线性关系,一是相关的强度,二是相关的方向。
设两变量有如下两组观测值:
X:x1,x2,…,xn;Y:y1,y2,…,yn,则
$ r = \frac{{\sum\limits_i^n {\left( {{X_i}-\bar X} \right)\left( {{Y_i}-\bar Y} \right)} }}{{\sqrt {\sum\limits_i^n {{{\left( {{X_i}-\bar X} \right)}^2}} } \sqrt {\sum\limits_i^n {{{\left( {{Y_i} - \bar Y} \right)}^2}} } }} $ | (1) |
式中:r即为变量X与Y的相关系数,r的取值范围为-1≤r≤1。相关系数的绝对值大小 (即|r|) 表示两变量直间的线性相关强度,相关系数r的正负表示相关的方向,r为正数,表示两变量间为正相关;r为负数,表示两变量间为负相关。
典型工况 (凸轮转速为900 r/min,喷油脉宽为8℃aA凸轮轴转角) 下,应用试验设计的思想,在MODDE环境下采取面中央合成设计 (central composite face,CCF) 的响应面 (response surface modeling,RSM) 方法,将各影响因素作为自变量,将电磁阀响应延迟时间作为响应变量,对各特性参数与电磁阀开启和关闭延迟时间的相关性进行了分析。
如图 5所示,特性参数对电磁阀关闭延迟时间都是正相关的,说明关闭延迟时间随特性参数增大而增加;弹簧预紧力、残余气隙为、阻尼、锥阀直径相关系数分别为0.55、0.50、0.42、0.33,表明这些参数会对关闭延迟时间影响较大。而对于开启延迟时间,阻尼、阀杆直径、弹簧刚度、电磁阀容积腔是正相关的,说明开启延迟时间随这些特性参数增大而增加;其余是负相关,说明开启延迟时间随这些特性参数增大而减小。其中弹簧预紧力、锥阀半锥角、阻尼、阀杆直径相关系数分别为-0.67、-0.5、0.38、0.33,表明该些参数会对开启延迟时间影响较大;电磁阀容积腔、弹簧刚度与关闭延迟时间的相关系数分别为0.039、0.036,对开启延迟时间相关系数分别为0.018、0.014,表明它们对电磁阀响应影响微弱。其中阻尼是由于燃油的粘滞性、流动性及机械偶件之间相互摩擦所导致,即使对电磁阀响应延迟影响较大,但在优化设计中无法将其设计到最小。综上所述,本文选取优化的关键特性参数:衔铁残余气隙、弹簧预紧力、锥阀半锥角、阀杆直径及锥阀直径。
4 电磁阀响应延迟多目标优化 4.1 NSGA-Ⅱ遗传算法NSGA-Ⅱ[5, 7-11]是一种精英策略非支配排序遗传算法,在NSGA的基础上加上了精英策略、密度值估计策略和快速非支配排序策略。其优点是:计算复杂性从O(mN3) 降至O(Mn2),其中,m表示目标函数的数目、N表示种群中个体的数目;具有最优保留机制并且不需要确定一个共享参数;提高了算法的收敛性。NSGA-Ⅱ算法得到的非劣解分布均匀,收敛性好,已经成为进化多目标优化领域的基准算法之一。NSGA-Ⅱ算法流程如图 6[12]所示。
4.2 多目标优化数学模型本文优化目标为开启延迟时间fopen(x) 和关闭延迟时间fclose(x) 两个函数值最小时对应的x值,优化的关键特性参数:锥阀直径D1、阀杆直径D2、弹簧预紧力Y、锥阀半锥角J、衔铁残余气隙R,则电磁阀响应延迟的优化目标函数用数学描述为
Minimize fopen(x)
Minimize fclose(x)
根据设计中对各参数的选取范围,确定了各关键特性参数的约束条件:7.5≤D1≤8.2;5≤D2≤6;50≤Y≤90;56≤J≤80;0.08≤R≤0.15,fi(x) 为第i(i=1,2,…,n) 个目标函数,其中x为关键特性参数,x=(D1,D2,Y,J,R)T。
最终优化目标的解:取得出解集中离原点距离最近的点为最优集,即
$ d = \min \sqrt {{f_{{\rm{open}}}}{{\left( {{x_0}} \right)}^2} + {f_{{\rm{close}}}}{{\left( {{x_0}} \right)}^2}} $ | (2) |
式中:x0=(D10,D20,Y0,J0,R0)T为多目标优化结果。
4.3 多目标优化模型图 7所示为多目标优化平台modeFRONTIER与AMESim联合仿真模型,首先实验设计随机产生输入变量的相互组合,即赋予遗传算法初始值,经过AMESim计算得出输出变量的值,优化算法再根据目标函数的要求进行选择优化参数的变异和交换等操作,从而达到优化目的。本研究共有5个输入变量,电控单体泵在典型工况下进行仿真,将开启、关闭延迟时间作为输出变量,并将开启、关闭延迟时间最小作为最终优化目标。
4.4 多目标优化结果分析图 8(a)所示为电磁阀响应延迟时间多目标优化解集,从中选取符合电磁阀开启、关闭延迟时间都最小的Pareto最优解集前沿,如图 8(b)所示,在Pareto前沿上选取满足优化目标的点,点A为优化前电磁阀响应延迟时间点,点C为关闭延迟最小点,点D为开启延迟最小点,但这两个点的参数配置均只能实现对单个目标的优化,为了同时实现对电磁阀开启和关闭延迟时间优化的目标,应用式 (2) 计算出d最小值的x0点,即点B作为最终优化结果点。
从优化结果中可看出,表 2给出了优化前后设计参数对比值,表 3所示为优化前后电磁阀关闭延迟时间和开启延迟时间。优化后,关闭延迟时间减小了6%,开启延迟时间减小了17.7%,循环喷油量减小了1.50 mm3,喷油压力峰值增加0.63 MPa。说明经过电磁阀延迟响应多目标优化,随着电磁阀响应的加快,喷油压力有所增加,对循环喷油量的控制也更加精确,可进一步改善电控单体泵的喷射特性。
参数 | 优化前 | 优化后 |
锥阀直径/mm | 7.60 | 8.20 |
阀杆直径/mm | 5.50 | 5.55 |
衔铁残余气隙/mm | 0.12 | 0.08 |
锥阀半锥角/(°) | 68.00 | 80.00 |
弹簧预紧力/N | 70.00 | 88.00 |
参数 | 关闭延迟时间/ms | 开启延迟时间/ms | 循环喷油量/mm3 | 喷油压力峰值/MPa |
优化前 | 0.499 | 0.650 | 104.43 | 113.50 |
优化后 | 0.469 | 0.535 | 102.93 | 114.13 |
优化率/% | 6 | 17.7 | — | — |
图 9给出了电磁阀优化前后电磁阀升程的对比,从图 8(b)中可以看出,优化后的电磁阀开启时间提前了0.621℃aA凸轮轴转角即0.115 ms,电磁阀关闭时间早于优化前0.162℃aA凸轮轴转角即0.03 ms,响应速度快于优化前。
优化前后电控单体泵的喷油压力及喷油速率对比如图 10和图 11所示,由于优化后,电磁阀关闭延迟时间减小,关闭时间提前,从而导致喷油时刻提前,优化后的喷油压力及喷油速率在同一凸轮转角时高于优化前。当喷油压力和喷油速率到达最大值后,由于电磁阀开启延迟时间减小,开启时间提前,燃油泄压速度加快,从而喷油压力和喷油速率迅速降低,更有利于循环喷油量的精确控制。
5 结论
1) 应用实验设计方法,通过对电磁阀延迟响应因素的相关性分析,得出电磁阀延迟响应关键因素:衔铁残余气隙、弹簧预紧力、锥阀半锥角、锥阀直径、阀杆直径及其阻尼。
2) 采用多目标优化算法NSGA-Ⅱ对电控单体泵电磁阀响应特性进行了优化,优化结果表明:关闭延迟时间减小了6%,关闭时间提前,喷油压力及喷油速率在同一凸轮转角时高于优化前。
3) 经过优化后,开启延迟时间减小了17.7%。开启时间提前,燃油泄压速度加快,喷油压力和喷油速率迅速降低,更有利于循环喷油量的精确控制。
[1] | MULEMANE A, HAN Joongsub, LU Paihsiu, et al. Modeling dynamic behavior of diesel fuel injection systems[C]//SAE Paper, 2004-01-0536, 2004. |
[2] | FAN L Y, MA X Z, TIAN B Q, et al. Quantitative analysis on cycle fuel injection quantity fluctuation of diesel engine electronic in-line pump system[C]//SAE Paper, 2010-01-0875, 2010. |
[3] |
张奇, 张科勋, 李建秋, 等. 电控柴油机电磁阀驱动电路优化设计[J].
内燃机工程, 2005, 26(2): 1–4.
ZHANG Qi, ZHANG Kexun, LI Jianqiu, et al. Optimization design of solenoid drive circuit for electronically controlled diesel engine[J]. Chinese internal combustion engine engineering, 2005, 26(2): 1–4. |
[4] |
张廷羽, 张国贤. 高速开关电磁阀的性能分析及优化研究[J].
机床与液压, 2006(9): 139–142.
ZHANG Tingyu, ZHANG Guoxian. Performance analysis and optimization of high speed solenoid valve[J]. Machine tool & hydraulics, 2006(9): 139–142. |
[5] |
范立云, 田丙奇, 马修真, 等. 电控单体泵喷射特性关键影响因素研究[J].
农业机械学报, 2011(9): 20–26.
FAN Liyun, TIAN Bingqi, MA Xiuzhen, et al. Key influence factors investigation on the electronic unit pump injection characteristics[J]. Transactions of the Chinese society of agricultural machinery, 2011(9): 20–26. |
[6] | DEB K, PBATAP A, AGARWAL S. A fast and elitist multiobjective genetic algorithm:NSGA-Ⅱ[J]. IEEE transactions on evolutionary computation, 2002, 6(2): 182–197. DOI:10.1109/4235.996017 |
[7] | JIANG Wenying, LIN Yan, CHEN Ming, et al. An ant colony optimization-genetic algorithm approach for ship pipe route design[J]. International shipbuilding progress, 2014(3/4): . |
[8] | GIAGKIOZIS I, FLEMING P J. Methods for multi-objective optimization:an analysis[J]. Information sciences, 2015: . |
[9] | SUI Haiteng, NIU Wentie, NIU Yaxiao, et al. Pipe-assembly approach for ships using modified NSGA-Ⅱ algorithm[J]. Computer aided drafting, design and manufacturing, 2016(2): 34–42. |
[10] | LIU H, GU F. A improved NSGA-Ⅱ algorithm based on sub-regional search[C]//Evolutionary computation 2011, "CEC"11. New Orleans, USA, 2011. |
[11] | DEB K, AGRAWAL S, PRATAP A, et al. A fast elitist non-dominated sorting genetic algorithm for multi-objective optimization:NSGA-Ⅱ[C]//Parallel problem solving from nature. Berlin, Germany:Springer-Verlag, 2000:849-858. |
[12] |
杜海峰. 基于DE和NSGA-Ⅱ的进化多目标优化算法及应用[D]. 西安: 西安电子科技大学, 2012.
DU Haifeng. Evolutionary multi-objective optimization algorithm and its application based on DE and NSGA-Ⅱ[D]. Xi'an:Xidian University, 2012. |
[13] | DHANALAKSHMI S, KANNAN S, MAHADEVAN K, et al. Application of modified NSGA-Ⅱ algorithm to combined economic and emission dispatch problem[J]. International journal of electrical power and energy systems, 2011, 33(4): 992–1002. DOI:10.1016/j.ijepes.2011.01.014 |
[14] | LI An'da, HE Zhen, ZHANG Yang. Bi-objective variable selection for key quality characteristics selection based on a modified NSGA-Ⅱ and the ideal point method[J]. Computers in industry, 2016, 82: 95–103. DOI:10.1016/j.compind.2016.05.008 |
[15] | NAJMEH A, SEYED M M, RAJA A R G, et al. Application of the NSGA-Ⅱ algorithm to a multi-period inventory-redundancy allocation problem in a series-parallel system[J]. Reliability engineering and system safety, 2016, 160: 1–10. |
[16] | VO-DUY T, DUONG-GIA D, HO-HUU V, et al. Multi-objective optimization of laminated composite beam structures using NSGA-Ⅱ algorithm[J]. Composite structures, 2017, 168: 498–509. DOI:10.1016/j.compstruct.2017.02.038 |
[17] | LIANG Xiaozhen, YAN Guangrong, CHU Hongzen. A research on optimization of cutting parameters based on NSGA-Ⅱ[J]. Advanced materials research, 2014, 936: 1687–1693. DOI:10.4028/www.scientific.net/AMR.936 |
[18] | LI Hongtao, NIU Wentie, FU Shengli, et al. Multiobjective optimization of steering mechanism for rotary steering system using modified NSGA-Ⅱ and fuzzy set theory[J]. Mathematical problems in engineering, 2015(1): 1–13. |