CN107269408B - 柴油机优化燃烧控制器及对仿真模型控制方法 - Google Patents
柴油机优化燃烧控制器及对仿真模型控制方法 Download PDFInfo
- Publication number
- CN107269408B CN107269408B CN201710338848.9A CN201710338848A CN107269408B CN 107269408 B CN107269408 B CN 107269408B CN 201710338848 A CN201710338848 A CN 201710338848A CN 107269408 B CN107269408 B CN 107269408B
- Authority
- CN
- China
- Prior art keywords
- output
- diesel engine
- time
- input
- formula
- 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.)
- Active
Links
- 238000002485 combustion reaction Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000004088 simulation Methods 0.000 title description 18
- 239000000446 fuel Substances 0.000 claims abstract description 69
- 238000002347 injection Methods 0.000 claims abstract description 45
- 239000007924 injection Substances 0.000 claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 230000005284 excitation Effects 0.000 claims abstract description 10
- 238000005457 optimization Methods 0.000 claims abstract description 10
- 239000000243 solution Substances 0.000 claims abstract description 5
- 230000008859 change Effects 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000003912 environmental pollution Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000010200 validation analysis Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013100 final test Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D37/00—Non-electrical conjoint control of two or more functions of engines, not otherwise provided for
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02B—INTERNAL-COMBUSTION PISTON ENGINES; COMBUSTION ENGINES IN GENERAL
- F02B37/00—Engines characterised by provision of pumps driven at least for part of the time by exhaust
- F02B37/12—Control of the pumps
- F02B37/24—Control of the pumps by using pumps or turbines with adjustable guide vanes
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/0025—Controlling engines characterised by use of non-liquid fuels, pluralities of fuels, or non-fuel substances added to the combustible mixtures
- F02D41/0047—Controlling exhaust gas recirculation [EGR]
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D41/1406—Introducing closed-loop corrections characterised by the control or regulation method with use of a optimisation method, e.g. iteration
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
- F02D41/1444—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases
- F02D41/146—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an NOx content or concentration
- F02D41/1461—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an NOx content or concentration of the exhaust gases emitted by the engine
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/30—Controlling fuel injection
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1412—Introducing closed-loop corrections characterised by the control or regulation method using a predictive controller
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1433—Introducing closed-loop corrections characterised by the control or regulation method using a model or simulation of the system
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/02—Input parameters for engine control the parameters being related to the engine
- F02D2200/06—Fuel or fuel supply system parameters
- F02D2200/0625—Fuel consumption, e.g. measured in fuel liters per 100 kms or miles per gallon
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/02—Input parameters for engine control the parameters being related to the engine
- F02D2200/10—Parameters related to the engine output, e.g. engine torque or engine speed
- F02D2200/1002—Output torque
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
一种柴油机优化燃烧控制器,属于柴油机技术领域。本发明的目的通过对柴油机的喷油量、喷油时刻、EGR阀的开度及VGT截面积进行控制,从而使得柴油机的NOx排放满足国IV排放标准,同时使其具有良好的动力性及燃油经济性能的柴油机优化燃烧控制器。本发明步骤是:激励数据的选取、通过输入输出Hankel矩阵构造***增量形式的预测输出方程、在得到增量型的预测矩阵、数据驱动预测控制器设计。本发明使用***的输入输出数据离线辨识得到***面向控制的模型,这大大地减少了传统的机理建模的时间;同时由于辨识得到的面向控制的模型具有线性的结构,这有利于MPC控制器的设计及其优化问题的在线求解。
Description
技术领域
本发明属于柴油机技术领域。
背景技术
随着汽车工业的迅速发展,世界范围内汽车数量也在急剧增加,2015年全球汽车保有量将超过11亿辆,而2050年将预计增长到35亿辆。如此庞大数量的汽车每年将消耗10亿多吨燃油,已超过世界石油年产量的1/3。以目前的开采速度计算,全球石油资源将会在50年后枯竭,并且美国能源部研究预测,2020年以后全球石油需求与常规石油供给之间将出现净缺口。可见能源危机迫在眉睫。伴随燃油的使用,汽车排放的尾气造成了严重的环境污染,并且危害着人类的健康。我国大部分地区出现的雾霾问题已引起了人们的高度重视[2]。能源短缺、环境污染和气候变暖是当前汽车和能源产业共同面临的巨大挑战。针对汽车的生产和销售,我国出台了最新的国Ⅴ排放标准,将对汽车的排放进行了更为严格的限制。对柴油机进行优化燃烧控制,能有效的减低柴油机的排放,同时能够有效提高柴油机的经济性能与动力性能,这对于缓解我国的环境污染问题及资源枯竭问题具有重要的意义。
发明内容
本发明的目的通过对柴油机的喷油量、喷油时刻、EGR阀的开度及VGT截面积进行控制,从而使得柴油机的NOx排放满足国IV排放标准,同时使其具有良好的动力性及燃油经济性能的柴油机优化燃烧控制器。
本发明步骤是:
①激励数据的选取:在柴油机燃烧模型中,将设计的喷油量、喷油时刻、EGR阀的开度、VGT的齿条位置给柴油机,开环运行模型,得到实际的NOx的排放量,发动机曲轴输出的扭矩,燃油消耗率;根据这些数据得到***的输入Hankel矩阵、,输出Hankel矩阵、、:
公式(1)中的表示***在2…N-2M+2时刻***的输出;公式(3)中的表示***在2+M…N-M+2时刻***的输出;其中第i时刻的***输出为,其中、、分别表示第i时刻燃油消耗率、NOx排放量以及柴油机输出的曲轴扭矩 ;
公式(2)中表示***在2……N-M+2时刻***输出的增量,其中第i时刻***输出的增量为,即第i时刻***的输出减去第i时刻的输出,其中、、分别表示第i时刻燃油消耗率、NOx排放量以及柴油机输出的曲轴扭矩;
公式(4)中的表示***在2……N-M+2时刻***输入的增量;(5)中的表示***在2+M…N+2 时刻***输入的增量,其中第i时刻***输入的增量为,即第i时刻***的输出减去第i时刻的输出,其中、、、分别表示第i时刻的喷油量、EGR阀的开度、VGT截面积及喷油时刻;
②通过输入输出Hankel矩阵构造***增量形式的预测输出方程为:
④数据驱动预测控制器设计:对目标函数对应的二次规划问题求解得到优化后的喷油量、喷油时刻、VGT的齿条位置、EGR阀的开度,目标函数的数学表达式如式(9)所示:
其中k1、k2、k3、k4为目标函数的加权系数,k、p分别为当前时刻、预测时域,RT为期望的柴油机输出曲轴扭矩,RNOx由式(10)计算得到
⑤将式(6)进行等量代换得式(11)
⑥通过式(11),得到如下所示的目标函数:
其中
式(13)描述的是对控制量的变化率的约束,其中、分别为控制量变化率的最大值与最小值,式(14)描述的是对控制量的约束,、分别为控制量的最大值与最小值,通过求解带约束(13)、(14)的优化问题(12),就可以得到最优的控制率。
本发明能很好地解决现有技术的三点问题,基于数据驱动预测控制算法能有效表面复杂的***机理建模,
1、使用***的输入输出数据离线辨识得到***面向控制的模型,这大大地减少了传统的机理建模的时间;同时由于辨识得到的面向控制的模型具有线性的结构,这有利于MPC控制器的设计及其优化问题的在线求解。
2、数据驱动MPC控制器能够对***未来的动态进行预测,然后根据其预测的***的动态做出相应的控制,所以使用MPC控制器能够很好克服惯性对柴油机输出的影响,从而对柴油机的燃烧过程实现优化。
3、柴油机的优化燃烧过程是一个多输入、多输出,带耦合的多目标优化问题。数据驱动MPC控制器具有处理带耦合的多目标优化问题的能力。所以使用数据驱动MPC能够同时对柴油机的燃油经济性、动力性及排放性能进行优化。
附图说明
图1 是本发明所述的基于数据驱动MPC的柴油机优化燃烧控制框图;其中控制量是喷油量、喷油时刻、EGR阀的开度、VGT截面积;
图2 是本发明所述的集中式带EGR+VGT的柴油机GT模型的示意图,该模型由环境设置模块、圆管模块、VGT模块、中冷器模块、三通管模块、EGR模块、四缸2升排量的柴油机模块及尾气后处理模块构成;
图3 是本发明在Simulink中搭建的柴油机优化燃烧数据驱动MPC控制器(1),主要是通过MATLAB中的M文件编译程序来实现数据驱动MPC控制算法;
图4是本发明在Simulink中搭建的柴油机优化燃烧数据驱动MPC控制器(2),主要是通过MATLAB中的M文件编译程序来实现数据驱动MPC控制算法;
图5 是本发明对柴油机进行激励时喷油量的输入,单位为mg指的是每循环的喷油量,横坐标为时间,单位为s;
图6 是本发明对柴油机进行激励时EGR阀的开度的输入,横坐标是时间,单位为s;
图7 是本发明对柴油机进行激励时喷油滞后角的输入,单位为度(°),横坐标为时间,单位为s;
图8 是本发明对柴油机进行激励时VGT齿条位置的输入,横坐标为时间,单位为s;
图9 是本发明对柴油机进行激励时燃油消耗率的输出,其单位为g/kwh,横坐标为时间,单位为s;
图10 是本发明对柴油机进行激励时NOx排放的输出,其单位为g/kwh,横坐标为时间,单位为s;
图11 是本发明对柴油机进行激励时曲轴扭矩的输出,其单位为N.m,横坐标为时间,单位为s;
图12 是本发明对辨识模型进行验证时喷油量的输入,单位为mg指的是每循环的喷油量,横坐标为时间,单位为s;
图13 是本发明对辨识模型进行验证时EGR阀的开度的输入,横坐标是时间,单位为s;
图14 是本发明对辨识模型进行验证时喷油滞后角的输入,单位为度(°),横坐标为时间,单位为s;
图15 是本发明对辨识模型进行验证时VGT齿条位置的输入,横坐标为时间,单位为s;
图16 是本发明对辨识模型进行验证时燃油消耗率的输出,其单位为g/kwh,横坐标为时间,单位为s,其中实线代表GT模型的输出,虚线代表辨识模型的输出;
图17 是本发明对辨识模型进行验证时NOx排放的输出,其单位为g/kwh,横坐标为时间,单位为s,其中实线代表GT模型的输出,虚线代表辨识模型的输出;
图18 是本发明对辨识模型进行验证时曲轴扭矩的输出,其单位为N.m,横坐标为时间,单位为s,其中实线代表GT模型的输出,虚线代表辨识模型的输出;
图19 是本发明在driving模式下对柴油机进行控制时喷油量的输入,单位为mg指的是每循环的喷油量,横坐标为时间,单位为s;
图20 是本发明在driving模式下对柴油机进行控制时EGR阀的开度的输入,横坐标是时间,单位为s;
图21 是本发明在driving模式下对柴油机进行控制时喷油滞后角的输入,单位为度(°),横坐标为时间,单位为s;
图22 是本发明在driving模式下对柴油机进行控制时VGT齿条位置的输入,横坐标为时间,单位为s;
图23 是本发明在driving模式下对柴油机进行控制时燃油消耗率的输出,其单位为g/kwh,横坐标为时间,单位为s;
图24 是本发明在driving模式下对柴油机进行控制时NOx排放的输出,其单位为g/kwh,横坐标为时间,单位为s;
图25 是本发明在driving模式下对柴油机进行控制时曲轴扭矩的输出,其单位为N.m,横坐标为时间,单位为s,其中实线代表GT模型的输出,虚线代表辨识模型的输出;
图26 是本发明在load模式下对柴油机进行控制时喷油量的输入,单位为mg指的是每循环的喷油量,横坐标为时间,单位为s;
图27 是本发明在load模式下对柴油机进行控制时EGR阀的开度的输入,横坐标是时间,单位为s;
图28 是本发明在load模式下对柴油机进行控制时喷油滞后角的输入,单位为度(°),横坐标为时间,单位为s;
图29 是本发明在load模式下对柴油机进行控制时VGT齿条位置的输入,横坐标为时间,单位为s;
图30 是本发明在load模式下对柴油机进行控制时燃油消耗率的输出,其单位为g/kwh,横坐标为时间,单位为s;
图31 是本发明在load模式下对柴油机进行控制时NOx排放的输出,其单位为g/kwh,横坐标为时间,单位为s;
图32 是本发明在load模式下对柴油机进行控制时柴油机转速的输出,其单位为rpm,横坐标为时间,单位为s,其中实线代表GT模型的输出,虚线代表辨识模型的输出。
具体实施方式
本发明首先针对柴油机的经济性能、动力性能及燃油经济性,提出了相应的控制指标,并选择了相应的控制输入;其次,根据***的特性设计了适当的激励数据,以保证对***稳态、动态特性的充分激励;接着,根据***的阶跃响应特性,选择合适大小的输入输出矩阵对***进行辨识;然后,考虑执行器的约束,利用模型预测控制算法构造相应的代价函数;最后,通过求解目标函数对应的最优问题,获得控制输入并将其作用于***,从而实现对***的控制。
本发明所述的基于数据驱动MPC控制器的柴油机优化燃烧是通过GT-power与Simulink的联合仿真实现。其中GT-power是一个商用的复杂***建模仿真平台,本发明利用其搭建一个高保真的柴油机燃烧模型,在仿真实验中用以代替一款真实的柴油机;MATLAB/ Simulink 则是用于控制器的仿真模型搭建,即通过 Simulink 编程来完成基于数据驱动MPC的柴油机优化燃烧控制器的搭建。
从功能上说本发明可以包括以下几部分:控制器仿真模块和带EGR-VGT的柴油机燃烧模型。下面详细说明各部分的作用:
带EGR-VGT的柴油机燃烧模型主要用于离线仿真得到能够体现***特性的输入输出激励数据,从而通过子空间辨识的方法得到面向控制的模型。
控制模块的主要作用是对柴油机燃烧模型的输入、输出量进行采集,然后通过求解相应的二次规划问题得到***最优的控制信号(喷油量、喷油时刻、EGR阀的开度、VGT的齿条位置);
本发明中基于数据驱动MPC控制器的柴油机优化燃烧控制的控制框图,如图1所示。图中的控制器是在Simulink中搭建的,控制器的输入是喷油量、喷油时刻、EGR阀的开度及VGT齿条的位置,而实际的燃油消耗率、NOx的排放及发动机曲轴输出的扭矩是实时反馈回控制器的。考虑到实际存在的约束,喷油量的取值范围为0mg~40mg,喷油时刻-10~10,EGR阀开度对应的值范围0~0.3,VGT截面积对应的齿条位置范围0~1。
本发明的控制目标是,柴油机输出的曲轴扭矩跟踪上期望值,同时优化燃油消耗率以及使得柴油机排放的NOx满足国五排放标准。
根据上述步骤可以得到基于PC机的离线柴油机优化燃烧设计试验平台。该平台的搭建及运行过程如下所示:
一、仿真平台搭建
柴油机优化燃烧控制***的被控对象和控制器分别通过MATLAB/Simulink 和GT-power进行搭建,求解器选择分别为ode 3和Explicit-Runge-Kutta。仿真步长为定步长,本发明中选为0.04s。
二、联合仿真设置
MATLAB/simulink与GT-power的联合仿真有两种方式,一种是在GT-power中调用MATLAB/Simulink中的模型;另一种是在MATLAB/Simulink中调用GT-power的模型。本发明中,为了便于控制器的调试使用的是第二联合仿真的方式。要实现MATLAB/Simulink与GT-power的联合仿真,首先要将GT-power的安装路径添到MATLAB中。然后分别在GT-power界面及MATLAB/simulink界面添加相应的通信接口模块,将MATLAB/Simulink和GT-power间需要通信的变量连接到这个模块。最后就是在MATLAB/simulink中设置仿真步长,在进行联合仿真的过程仿真步长一定要设置为定步长。
三、集中式柴油机燃烧过程仿真模型搭建:采用模块化方法,在GT-power中搭建了带EGR-VGT的柴油机模型,其示意图如图2所示。该模型由环境设置模块、圆管模块、VGT模块、中冷器模块、三通管模块、EGR模块、四缸2升排量的柴油机模块及尾气后处理模块构成。首先在GT中新建一个文件,然后将环境设置模块1(EndEnvironment)拉入该文件中,通过该模块对环境温度及压力进行设置,具体的参数如表1 所示,其中Composition属性及Humidity Species属性分别使用GT自带的air对象及h2o-vap对象;然后通过圆管模块1将环境设置模块与环境设置模块与压缩机模块相连接,其中圆管模块1的参数如表1 所示;接着将压缩与中冷器连接在一起,其中中冷器的参数直接使用demo中的参数;接着通过三通管1将中冷器、四缸2L排放的柴油机模块的入口及EGR模块的出口连接在一起,其中四缸2L排放的柴油机模块的参数如表1 所示,其他的参数直接使用demo自带的参数,三通管1的参数如表1 所示;接着通过三通管模块2将四缸2L排放的柴油机模块的出口、涡轮机及EGR模块的入口连接在一起,三通管模块2的参数如表1 所示;然后将涡轮机通过圆管模块2与尾气后处理***相连接,尾气后处理***的参数直接使用demo的参数,圆管模块2的参数如表1所示,最后通过圆管模3块将尾气后处理模块与环境设置模块2连接在一起,其中圆管模块3的参数如表1 所示,环境设施模块2的参数与环境设置模块1的参数是一致的。
表1 带EGR+VGT的柴油机参数列表
四、柴油机优化燃烧控制方案
首先是控制目标的确定:使NOx的排放的均值小于2g/kwh,以满足国五排放标准对NOx排放的要求;同时使得柴油机输出的曲轴扭矩跟踪上期望值,以满足对柴油机动力性的要求;最后,在满足排放性能及动力性能的前提下,尽量降低燃油消耗。
控制量的确定:对于柴油机而言,其输出的扭矩主要由喷油量决定;为了提高柴油机的燃油经济性及动力性,需要优化柴油机的喷油时刻及VGT的截面;另一方面,为了减低柴油机NOx的排放,需要对EGR阀的开度进行控制,所以为了满足上述控制目标,本发明选择喷油量、喷油时刻、EGR率及VGT的截面积作为控制量。
控制器的选择:因为柴油机的燃烧过程复杂且存在较强的非线性,本发明使用增量型子空辨识的方法建模建立面向控制的模型;另一方面,考虑到MPC控制器能够很好的解决带约束、多目标优化问题,所以本发明使用基于增量型数据驱动MPC控制器。
确定控制目标、控制量及控制器之后就能得到本发明的柴油机优化燃烧控制方案,具体的控制框图如图1所示。
五、柴油机优化燃烧数据驱动预测控制器
①激励数据的选取:选取柴油机的喷油量、喷油时刻、EGR阀的开度、VGT的齿条位置(在GT模型中,VGT的截面积是由齿条位置决定的,齿条的位置在0~1的范围内变化)作为激励输入信号,在GT-power的柴油机燃烧模型中,将设计的喷油量、喷油时刻、EGR阀的开度、VGT的齿条位置给柴油机,开环运行模型,得到实际的NOx的排放量,发动机曲轴输出的扭矩,燃油消耗率;根据这些数据得到***的输入Hankel矩阵、,输出Hankel矩阵、、:
公式(1)中的表示***在2…N-2M+2时刻***的输出;公式(3)中的表示***在2+M…N-M+2时刻***的输出;因为公式(3)中的输出是***未来M时刻的值,所以公式(3)中的输出超前公式(1)中对应的输出M时刻。其中第i时刻的***输出为,其中、、分别表示第i时刻燃油消耗率、NOx排放量以及柴油机输出的曲轴扭矩 。
公式(2)中表示***在2……N-M+2时刻***输出的增量,其中第i时刻***输出的增量为,即第i时刻***的输出减去第i时刻的输出,其中、、分别表示第i时刻燃油消耗率、NOx排放量以及柴油机输出的曲轴扭矩;
公式(4)中的表示***在2……N-M+2时刻***输入的增量;(5)中的表示***在2+M…N+2 时刻***输入的增量,其中第i时刻***输入的增量为,即第i时刻***的输出减去第i时刻的输出,其中、、、分别表示第i时刻的喷油量、EGR阀的开度、VGT截面积及喷油时刻;
②通过输入输出Hankel矩阵构造***增量形式的预测输出方程为:
③在得到增量型的预测矩阵和后,利用增量型子空间预测方程(6),即得到***未来输出值;在得到如(6)所示的预测方程之后,本发明使用图5所示的激励输入对辨识模型的精度进行了验证,验证结果如图7所示。由图7可知,预测得到的***的输出与实际值吻合得很好。
④数据驱动预测控制器设计:对目标函数对应的二次规划问题求解得到优化后的喷油量、喷油时刻、VGT的齿条位置、EGR阀的开度,目标函数的数学表达式如式(9)所示:
其中k1、k2、k3、k4为目标函数的加权系数,k、p分别为当前时刻、预测时域,RT为期望的柴油机输出曲轴扭矩,RNOx由式(10)计算得到,由式(9)可得,目标函数中包括了:对柴油机输出曲轴扭矩的跟踪、对燃油消耗的优化以及对排放的跟踪优化。
⑤为了便于控制器的推导,将式(6)进行等量代换得式(11)
⑥通过式(11),得到如下所示的目标函数:
其中
式(13)描述的是对控制量的变化率的约束,其中、分别为控制量变化率的最大值与最小值,的值分别为40、0.3、10、1、-40、-0.3、-10、-1。式(14)描述的是对控制量的约束,、分别为控制量的最大值与最小值,的值分别为40、0.3、10、1、0、0、-10、0。通过求解带约束(13)、(14)的优化问题(12),就可以得到最优的控制率。如图3所示。
实验验证
模型预测控制的控制性能取决于预测模型的精度。为了验证模型的精度,本发明重新设计了一组激励输入,如图11、12、13、14所示,实验结果如图15、16、17所示。由图15、16、17可知,预测模型的NOx、燃油消耗率、曲轴扭矩输出能够很好地跟踪上实际***的输出值。
为了进一步验证柴油机优化燃烧数据驱动MPC控制器的控制性能,本发明使用FTP75工况分别在测功机模式及driving模式下对控制器的控制性能进行了仿真实验。FTP75工况被称为目前最合理的循环工况测试规则,其中美国加州政府起到了至关重要的作用。本着测试最真实数据的原则,美国FTP75工况设计了很多接近现实的试验内容。FTP75共由一个市区循环工况和两个补充循环工况组成。两个补充循环工况分别为SC03高温空调全负荷运转循环和US06高速、高加速度循环。最终试验结果由这三个试验结果通过不同的比例计算而成,因此这样的数据更接近实际使用。
1、测功机模式
GT中的测功机模式是通过柴油机的转速及柴油机的输入计算相应的扭矩输出。本发明使用FTP75工况下的柴油机转速作为给定输入,以验证基于数据驱动MPC控制器的在复杂工况的性能,实验结果如图22、23、24所示,其最优的控制量如图18、19、20、21所示。由图24可知,柴油机实际输出的曲轴扭矩能够很好地跟踪上期望值;同时,通过计算可知,NOx的平均排放量为1.9152g/kwh ,该值能够满足国五排放标准对于NOx排放的要求;燃油消耗率的平均值为262.4231g/kwh。这时本发明认为,在测功机模式下基于数据驱动MPC的柴油机优化燃烧控制策略能够同时提高柴油机的动力性能、经济性能及排放性能。
2、Driving模式
GT中的driving模式是通过柴油机输入的负载扭矩及柴油机的输入计算柴油机的转速。在GT中柴油机的转速可以由下面的式(15)得到。
其中,为柴油机的角加速度,该值可以由期望的柴油机的转速得到;J为柴油机的转动惯量;为柴油机输出的曲轴扭矩;为柴油机输入的负载扭矩。在实验过程中,柴油机的负载扭矩为认为给定的一个常值;期望的柴油机转速为FTP75工况下的柴油机转速;实验结果如图29、30、31所示,其最优的控制输入如图25、26、27、28所示。
由图31可知柴油机在driving模式下分别能够很好地跟踪上期望的柴油机转速;
同时由图30可以得到NOx排放的均值为1.9225g/kwh,该值能够满足国五排放标准对于NOx排放的要求;由图29可以得到燃油消耗率的平均值为262.1950g/kwh。这时本发明认为,在driving模式下基于数据驱动MPC的柴油机优化燃烧控制策略能够同时提高柴油机的动力性能、经济性能及排放性能。
本发明使用基于数据驱动MPC的控制器对四缸2L排放的柴油机的燃烧过程进行控制。首先使用了***的输入输出数据离线辨识得到了面向控制的模型,并通过实验对该模型的精度进行了验证,这一方面减少了传统机理建模的时间,另一方面提高了模型的精度;然后使用上述预测模型设计了MPC控制器,该控制器能够根据***未来***的动态超前产生控制率,从而能够很好地克服惯性对柴油机输出的影响;另一方面,该控制器能够处理多输入、多输出、带耦合的多目标优化问题,所以它能将柴油机的燃油经济性、动力性及排放性能统一到一个目标函数中,从而对柴油机的燃油经济性、动力性及排放性能进行优化。
Claims (1)
1.一种柴油机优化燃烧控制器,其特征在于:
①激励数据的选取:在柴油机燃烧模型中,将设计的喷油量、喷油时刻、EGR阀的开度、VGT的齿条位置给柴油机,开环运行模型,得到实际的NOx的排放量,发动机曲轴输出的扭矩,燃油消耗率;根据这些数据得到***的输入Hankel矩阵ΔIp、ΔIf,输出Hankel矩阵Op、Of、ΔOp:
公式(1)中的O2,…,ON-2M+2表示***在2,…,N-2M+2时刻的***输出;公式(3)中的O2+M,…,ON-M+2表示***在2+M,…,N-M+2时刻的***输出;其中第i时刻的***输出为其中分别表示第i时刻燃油消耗率、NOx排放量以及柴油机输出的曲轴扭矩;
公式(2)中ΔO2.....ΔON-M+2表示***在2……N-M+2时刻***输出的增量,其中第i时刻***输出的增量为即第i-1时刻***的输出减去第i时刻的输出,其中分别表示第i时刻燃油消耗率、NOx排放量以及柴油机输出的曲轴扭矩;
公式(4)中的ΔI2,…,ΔIN-2M+2表示***在2,…,N-2M+2时刻***输入的增量;公式(5)中的ΔI2+M,…,ΔIN-M+2表示***在2+M,…,N-M+2时刻***输入的增量,其中第i时刻***输入的增量为即第i-1时刻***的输入减去第i时刻的输入,其中分别表示第i时刻的喷油量、EGR阀的开度、VGT截面积及喷油时刻;
②通过输入输出Hankel矩阵构造***增量形式的预测输出方程为:
其中ΔOf=Of-Op,两个预测矩阵系数ΔLw和ΔLu的解由下式求得:
③在得到增量型的两个预测系数矩阵ΔLw和ΔLu后,利用增量型子空间预测方程(6),即得到***未来输出值;
④数据驱动预测控制器设计:对目标函数对应的二次规划问题求解得到优化后的喷油量、喷油时刻、VGT的齿条位置、EGR阀的开度,目标函数的数学表达式如式(9)所示:
其中k1、k2、k3、k4为目标函数的加权系数,k、p分别为当前时刻、预测时域,RT为期望的柴油机输出曲轴扭矩,RNOx由式(10)计算得到,J代表目标函数;
其中,SumNOx为前i时刻总的NOx的排放量;
⑤将式(6)进行等量代换得式(11)
其中F项为***自由相应部分,SΔIf项被称为***的控制相应;
⑥通过式(11),得到如下所示的目标函数:
其中
式(13)描述的是对控制量的变化率的约束,其中δImax、δImin分别为控制量变化率的最大值与最小值,式(14)描述的是对控制量的约束,Imax、Imin分别为控制量的最大值与最小值,通过求解带约束(13)、(14)的优化问题(12),就可以得到最优的控制率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710338848.9A CN107269408B (zh) | 2017-05-15 | 2017-05-15 | 柴油机优化燃烧控制器及对仿真模型控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710338848.9A CN107269408B (zh) | 2017-05-15 | 2017-05-15 | 柴油机优化燃烧控制器及对仿真模型控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107269408A CN107269408A (zh) | 2017-10-20 |
CN107269408B true CN107269408B (zh) | 2022-08-05 |
Family
ID=60063977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710338848.9A Active CN107269408B (zh) | 2017-05-15 | 2017-05-15 | 柴油机优化燃烧控制器及对仿真模型控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107269408B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110329258B (zh) * | 2019-07-23 | 2021-04-30 | 吉林大学 | 智能驾驶汽车节能减排协调控制方法 |
CN113448318B (zh) * | 2021-07-07 | 2022-08-16 | 江铃汽车股份有限公司 | 一种车辆下线故障诊断控制方法 |
CN114934848B (zh) * | 2022-06-16 | 2023-04-07 | 哈尔滨工程大学 | 一种面向柴油机燃烧性能优化控制的模糊神经网络建模方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102588129A (zh) * | 2012-02-07 | 2012-07-18 | 上海艾铭思汽车控制***有限公司 | 高压共轨柴油机氮氧化物和颗粒排放的优化协同控制方法 |
CN105629736A (zh) * | 2016-03-22 | 2016-06-01 | 东南大学 | 数据驱动的火电机组scr脱硝扰动抑制预测控制方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3383874B2 (ja) * | 1993-07-19 | 2003-03-10 | 克享 小西 | ディーゼル機関の燃焼シミュレーション方法 |
DE10020448B4 (de) * | 2000-04-26 | 2005-05-04 | Daimlerchrysler Ag | Verfahren und Vorrichtung zur Optimierung des Betriebs eines Verbrennungsmotors |
US7184877B1 (en) * | 2005-09-29 | 2007-02-27 | International Engine Intellectual Property Company, Llc | Model-based controller for auto-ignition optimization in a diesel engine |
DE102006061936A1 (de) * | 2006-12-29 | 2008-07-03 | Robert Bosch Gmbh | Verfahren und System zur Simulation des Betriebs einer Brennkraftmaschine |
CN104401232B (zh) * | 2014-12-21 | 2016-06-22 | 吉林大学 | 基于数据驱动预测控制的电动汽车扭矩优化方法 |
CN105221275B (zh) * | 2015-08-20 | 2017-10-13 | 吉林大学 | 基于模型辨识的涡轮增压汽油机控制*** |
-
2017
- 2017-05-15 CN CN201710338848.9A patent/CN107269408B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102588129A (zh) * | 2012-02-07 | 2012-07-18 | 上海艾铭思汽车控制***有限公司 | 高压共轨柴油机氮氧化物和颗粒排放的优化协同控制方法 |
CN105629736A (zh) * | 2016-03-22 | 2016-06-01 | 东南大学 | 数据驱动的火电机组scr脱硝扰动抑制预测控制方法 |
Non-Patent Citations (1)
Title |
---|
高炉热风炉燃烧输入输出数学模型的应用;丁洪起等;《钢铁》;20160315(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107269408A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jankovic et al. | Robust nonlinear controller for turbocharged diesel engines | |
Xu et al. | Digital twin-driven optimization of gas exchange system of 2-stroke heavy fuel aircraft engine | |
Dong et al. | Experiment and simulation investigation on energy management of a gasoline vehicle and hybrid turbocharger optimization based on equivalent consumption minimization strategy | |
CN107269408B (zh) | 柴油机优化燃烧控制器及对仿真模型控制方法 | |
CN104408271A (zh) | 一种基于模型的汽油机标定方法 | |
Norouzi et al. | Deep learning based model predictive control for compression ignition engines | |
Yin et al. | Model predictive control of an advanced multiple cylinder engine with partially premixed combustion concept | |
CN104514637A (zh) | 动力传动系控制*** | |
Nikzadfar et al. | Investigating a new model-based calibration procedure for optimizing the emissions and performance of a turbocharged diesel engine | |
CN111608790B (zh) | 电辅助涡轮增压柴油机气路优化控制*** | |
Isermann et al. | Mechatronic combustion engines–from modeling to optinlal control | |
Chiara et al. | An exhaust manifold pressure estimator for a two-stage turbocharged diesel engine | |
Sakai et al. | Torque demand control by nonlinear MPC with constraints for vehicles with variable valve lift engine | |
Wang et al. | Adaptive air-fuel ratio control with MLP network | |
Duraiarasan et al. | Real-time embedded models for simulation and control of clean and fuel-efficient heavy-duty diesel engines | |
Keller et al. | Physics-based modeling and mpc for the air path of a two-stage turbocharged si engine with low pressure egr | |
Li et al. | Air path system control of turbocharged gasoline engine based on fuzzy PID | |
Norouzi et al. | Integration of Deep Learning and Nonlinear Model Predictive Control in Emission reduction of Compression Ignition Combustion Engines: A Simulative Study | |
Kuzmych et al. | Nonlinear control for a diesel engine: A CLF-based approach | |
Zhang et al. | Method of Turbocharger Emulation on Engine Test and Application to Turbocompound System Optimisation | |
Sun et al. | Integrated control of torque and emission of a diesel engine based on LPV‐MPC | |
Topaloglou et al. | Transient load share management of a diesel electric hybrid powertrain for ship propulsion | |
Youssef et al. | Model based control of turbochargers: Application to a Diesel HCCI engine | |
Cieslar et al. | Model based control for closed loop testing of 1-D engine simulation models | |
Huan et al. | Airpath prediction model and predictive control of turbocharged gasoline engine. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |