EP1729001A1 - Estimation method using a non-linear adaptive filter of air-fuel mixture richness in a cylinder of a combustion engine - Google Patents
Estimation method using a non-linear adaptive filter of air-fuel mixture richness in a cylinder of a combustion engine Download PDFInfo
- Publication number
- EP1729001A1 EP1729001A1 EP06290558A EP06290558A EP1729001A1 EP 1729001 A1 EP1729001 A1 EP 1729001A1 EP 06290558 A EP06290558 A EP 06290558A EP 06290558 A EP06290558 A EP 06290558A EP 1729001 A1 EP1729001 A1 EP 1729001A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- air
- richness
- cylinders
- exhaust
- exhaust manifold
- 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
Images
Classifications
-
- 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/008—Controlling each cylinder individually
-
- 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/1454—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 oxygen content or concentration or the air-fuel ratio
- F02D41/1458—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 oxygen content or concentration or the air-fuel ratio with determination means using an estimation
-
- 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/1413—Controller structures or design
- F02D2041/1415—Controller structures or design using a state feedback or a state space representation
- F02D2041/1416—Observer
-
- 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/1413—Controller structures or design
- F02D2041/143—Controller structures or design the control loop including a non-linear model or compensator
-
- 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/1413—Controller structures or design
- F02D2041/1431—Controller structures or design the system including an input-output delay
-
- 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
- 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/1402—Adaptive control
-
- 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/1454—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 oxygen content or concentration or the air-fuel ratio
Definitions
- the present invention relates to a method for estimating the fuel richness of each cylinder of an internal combustion engine injection, from a measurement of the wealth downstream of the collector and an adaptive nonlinear filter.
- the knowledge of wealth characterized by the ratio of the mass of fuel on the air mass, is important for all vehicles, whether they are petrol engines because it conditions a good combustion of the mixture when it is close of 1, or for vehicles with diesel engines for which the interest of the knowledge of the wealth is different since they work with poor mixture (wealth lower than 1).
- catalysts using a NOx trap lose their effectiveness over time. In order to return to optimum efficiency, the richness must be kept close to 1 for a few seconds, then return to normal operation at a lean mixture. Depollution by DeNOx catalysis therefore requires precise control of the cylinder-by-cylinder richness.
- a probe placed at the outlet of the turbine (turbocharged engine) and upstream of the NOx trap, gives a measure of the average richness by the exhaust process. This measurement, being very filtered and noisy, is used for the control of the masses injected into the cylinders during the phases of richness equal to 1, each cylinder then receiving the same mass of fuel.
- An engine control can thus, from the reconstructed wealth, adapt the fuel masses injected into each of the cylinders so that the wealth is balanced in all the cylinders.
- the object of the present invention is to model the exhaust process more finely so as, on the one hand, to dispense with the identification step, and on the other hand to bring more robustness to the wealth estimation model. , and this for all operating points of the engine.
- the invention also makes it possible to perform a measurement every 6 ° of rotation of the crankshaft and thus to have a high frequency information of the measurement of richness, without falling into the measurement noise.
- the physical model may comprise at least the following three types of variables: the total mass of gas in the exhaust manifold ( M T ), the fresh air mass in the exhaust manifold ( M air ) and the riches in each cylinder ( ⁇ i ).
- This mode can then comprise at least the two following types of output data: the total mass of gas in the exhaust manifold ( M T ) and the mass flow rates leaving said cylinders ( d i ).
- the measured richness ( ⁇ ) can be estimated as a function of the total mass of gas in the exhaust manifold ( M T ) and the fresh air mass in the exhaust manifold ( M air ).
- the estimation of the value of the richness in each of the cylinders may then comprise a real-time correction of the estimate of the total mass of gas in the exhaust manifold ( M T ), the estimation of the mass of the fresh air in the exhaust manifold ( M air ) and estimating the value of the richness in each of the cylinders ( ⁇ i ).
- the method can be applied to an engine control to adapt the fuel masses injected into each of the cylinders to adjust the richness in all the cylinders.
- the composition of the exhaust gas depends on the amount of fuel and air introduced into the combustion chamber, the fuel composition and the development of the combustion.
- the richness probe measures the concentration of O 2 inside a diffusion chamber, connected to the exhaust pipe by a diffusion barrier made of porous materials. This configuration may induce differences depending on the location of the chosen probe, in particular because of temperature variations and / or pressures in the vicinity of the richness probe.
- the measured wealth ( ⁇ ) is connected to the mass of air (or to the air flow) around the probe and to the total mass ( or at the total rate).
- the model is based on a three-gas approach: air, fuel and flue gas.
- the lean mixture richness formula is used in the estimator, at the level of integration of richness in equation (7), neglecting a very small portion of the air ( ⁇ 3%).
- the invention is not limited to this mode, in fact, the formula is continuous in the vicinity of a richness equal to 1, and its inversion does not pose a problem for rich mixtures.
- AMESim is a 0D modeling software, particularly well suited to thermal and hydraulic phenomena. It allows to model volumes, behaviors or restrictions.
- the basic tubing, restriction and volume modeling blocks are described in the AMESim "Thermal Pneumatic Library" user manual. Standard equations are used to calculate a flow through a restriction and energy and mass conservations. In addition, the model takes into account gas inertia, which is important for studying the dynamics of gas composition.
- a unique real-time physical model is defined for modeling the overall system, that is to say the entire path of the exhaust gases, from the cylinders to the downstream exhaust from the turbine, through the collector.
- the exhaust manifold is modeled according to a volume in which there is conservation of the mass. It is assumed that the temperature is substantially constant, and determined from an abacus function of the load and the engine speed.
- Model to determine the flow passing through the turbine model of the turbine
- the turbine is modeled according to a flow passing through a flow restriction.
- the flow rate in the turbine is generally given by mapping (abacus) as a function of the turbine speed and the upstream / downstream pressure ratio of the turbine.
- the parameters of the function f are optimized by correlation with the mapping of the turbine.
- the first equation contains one unknown: M T.
- the second contains two: M air and ⁇ i . This leads to the additional assumptions described below.
- the unknowns of the physical model are ultimately M T , M air and ⁇ i .
- the output data of the physical model is M T and d i .
- the physical model (5) is non-linear, and it is impossible to solve such a system in real time. It is therefore necessary to use an estimator, rather than seeking to directly calculate the unknowns of the system.
- the choice of the estimator according to the invention is based on the fact that the structure of the system is linear as a function of the wealth in the cylinders ⁇ i (the air mass variation is linear as a function of ⁇ i ).
- a particularly suitable technique is to use an adaptive filter.
- the method according to the invention proposes to construct an estimator based on an adaptive filter. This estimator ultimately allows an estimation of the cylinder to cylinder richness from the measurement of wealth by the sensor located behind the turbine.
- the principle of the estimator is to converge the physical model (5), and consequently the riches ⁇ i towards reality.
- the model (5) outputs M T and M Air , and we also have input parameters Y.
- the estimator therefore compares the output values of the RTM model with the input values, then make the appropriate corrections.
- L 1 , L 2 , L ⁇ adjustment parameters, making it possible to control the speed of convergence of the solution to the three unknowns. These are strictly positive real parameters. These parameters are set manually to obtain a good compromise between the speed of convergence and the low sensitivity to measurement noise.
- the estimator thus constructed makes it possible to correct in real time M T , M air and ⁇ , from a first value of M T provided by the RTM model and from the measurement of richness made by the probe.
- the system (8) is numerically solved in real time, the calculator using an explicit Euler discretization, well known to those skilled in the art.
- FIGS. 2A and 2B show below the wealth of references ⁇ i ref AMESIM data by a function of time (T) and above the results of the estimator ( ⁇ i) in function of time (T).
- the four curves correspond to each of the four cylinders.
- the performance of the estimator based on the adaptive filter is very good. However, there is a slight difference in phase, due to the inertia of the gas which is not taken into account in the present model. It is therefore proposed to complete the model and the estimator by an estimator of the exhaust delay time.
- the estimator implemented as described above does not allow the estimation method to take into account the delay time between the cylinder exhaust and the signal acquired by the probe.
- the delay time comes from several sources: transport time in the pipes and through the volumes, dead time of the measuring probe.
- the penalty is given by ⁇ . If there is a positive variation in the estimated wealth value for cylinder 2, then the delay time between the estimator and the measurements is positive. If there is a variation on cylinder 3, the delay is negative and the penalty is negative. A variation of the cylinder 4 can be considered as a consequence of a positive or negative delay.
- the delay D applied to the output variables of the RTM model is an additive delay, it is computed by least squares by minimizing J k .
- the criterion J k is controlled to zero by a PI (Proportional Integral) controller on the delay of the estimator.
- PI Proportional Integral
- FIGS. 4A and 4B illustrate the estimation of the cylinder to cylinder richness by the estimator previously described at 1500tr / min average load. These figures show up the wealth of references ⁇ i ref versus time (T) and below the results of the estimator ( ⁇ i) in function of time (T). The four curves correspond to each of the four cylinders.
- the present invention relates to an estimation method comprising the construction of an estimator, making it possible, from the measurement of the richness of the probe ( ⁇ ) and the total mass of gas information inside the collector ( M T ), to estimate the wealth at the output of the four cylinders ( ⁇ i ).
- the estimator thus produced is efficient, and above all does not require any additional adjustment in the case of change of the operating point. No identification phase is necessary, only a measurement and model noise adjustment must be made once and only once.
- a delay time controller is put in parallel with the estimator, making it possible to reset the delay time following a step of injection time on a cylinder. This allows optimal calibration of the estimator, for example before a rich phase equal to 1.
- the invention also makes it possible to perform a measurement every 6 ° of rotation of the crankshaft and thus to have a high frequency information of the measurement of richness, without falling into the measurement noise.
- the high frequency representation makes it possible to take into account the pulsating effect of the system.
- the modeled system is periodic and makes it possible to obtain an estimator with a better dynamics: one anticipates the pulsation of the escapement.
- the invention makes it possible to reduce the calculation time by a factor of about 80 compared to the previous methods.
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
- Exhaust Silencers (AREA)
- Testing Of Engines (AREA)
Abstract
Description
La présente invention concerne une méthode d'estimation de la richesse en carburant de chaque cylindre d'un moteur à combustion interne à injection, à partir d'une mesure de la richesse en aval du collecteur et d'un filtre non linéaire adaptatif.The present invention relates to a method for estimating the fuel richness of each cylinder of an internal combustion engine injection, from a measurement of the wealth downstream of the collector and an adaptive nonlinear filter.
La connaissance de la richesse, caractérisée par le rapport de la masse de carburant sur la masse d'air, est importante pour tous les véhicules, qu'ils soient à motorisation essence puisqu'elle conditionne une bonne combustion du mélange lorsqu'elle est proche de 1, ou pour les véhicules à motorisation diesel pour lesquels l'intérêt de la connaissance de la richesse est différent puisqu'ils fonctionnent à mélange pauvre (richesse inférieure à 1). En particulier, les catalyseurs utilisant un piège à NOx perdent de leur efficacité avec le temps. Afin de revenir à une efficacité optimale, la richesse doit être maintenue proche de 1 pendant quelques secondes, pour revenir ensuite en fonctionnement normal à un mélange pauvre. La dépollution par catalyse DeNOx nécessite donc un pilotage précis de la richesse cylindre par cylindre.The knowledge of wealth, characterized by the ratio of the mass of fuel on the air mass, is important for all vehicles, whether they are petrol engines because it conditions a good combustion of the mixture when it is close of 1, or for vehicles with diesel engines for which the interest of the knowledge of the wealth is different since they work with poor mixture (wealth lower than 1). In particular, catalysts using a NOx trap lose their effectiveness over time. In order to return to optimum efficiency, the richness must be kept close to 1 for a few seconds, then return to normal operation at a lean mixture. Depollution by DeNOx catalysis therefore requires precise control of the cylinder-by-cylinder richness.
Pour ce faire, une sonde, placée en sortie de la turbine (moteur turbocompressé) et en amont du piège à NOx, donne une mesure de la richesse moyennée par le processus d'échappement. Cette mesure, étant très filtrée et bruitée, est utilisée pour le contrôle des masses injectées dans les cylindres lors des phases de richesse égale à 1, chaque cylindre recevant alors la même masse de carburant.To do this, a probe, placed at the outlet of the turbine (turbocharged engine) and upstream of the NOx trap, gives a measure of the average richness by the exhaust process. This measurement, being very filtered and noisy, is used for the control of the masses injected into the cylinders during the phases of richness equal to 1, each cylinder then receiving the same mass of fuel.
Afin de contrôler d'une manière plus précise, et surtout individuelle, l'injection des masses de carburant dans les cylindres, une reconstruction de la richesse dans chaque cylindre est indispensable. L'implantation de sondes de richesse en sortie de chaque cylindre n'étant pas envisageable sur véhicule étant donné leur prix de revient, la mise en place d'un estimateur fonctionnant à partir des mesures d'une seule sonde, permet avantageusement de connaître séparément les richesses de chaque cylindre.In order to control in a more precise way, and especially individual, the injection of the masses of fuel in the cylinders, a reconstruction of the wealth in each cylinder is essential. Since the implantation of wealth probes at the outlet of each cylinder is not possible on a vehicle because of their cost price, the implementation of an estimator operating from the measurements of a single probe advantageously makes it possible to know separately the riches of each cylinder.
Un contrôle moteur pourra ainsi, à partir des richesses reconstruites, adapter les masses de carburant injectées dans chacun des cylindres afin que les richesses soient équilibrées dans tous les cylindres.An engine control can thus, from the reconstructed wealth, adapt the fuel masses injected into each of the cylinders so that the wealth is balanced in all the cylinders.
On connaît le document
L'estimation de la richesse dans les cylindres est uniquement conditionnée par les coefficients d'une matrice, coefficients qui doivent être identifiés hors ligne grâce à un algorithme d'optimisation. De plus, à chaque point de fonctionnement (régime/charge) correspond un réglage différent de la matrice, donc une identification de ses paramètres. Cet estimateur nécessite donc de mettre en place de lourds moyens d'essais (avec 5 sondes de richesse) d'acquisition, et s'avère peu robuste dans le cas d'un changement de moteur.The estimation of the richness in the cylinders is only conditioned by the coefficients of a matrix, coefficients which must be identified offline thanks to an optimization algorithm. In addition, at each operating point (speed / load) corresponds a different setting of the matrix, so an identification of its parameters. This estimator therefore requires the establishment of heavy means of testing (with 5 wealth probes) acquisition, and is not robust in the case of a change of engine.
La présente invention a pour objet de modéliser plus finement le processus d'échappement afin, d'une part de se passer de l'étape d'identification, et d'autre part d'apporter plus de robustesse au modèle d'estimation de richesse, et ceci pour tous les points de fonctionnement du moteur. L'invention permet en outre d'effectuer une mesure tous les 6°de rotation du vilebrequin et donc d'avoir une information haute fréquence de la mesure de richesse, sans pour autant tomber dans le bruit de mesure.The object of the present invention is to model the exhaust process more finely so as, on the one hand, to dispense with the identification step, and on the other hand to bring more robustness to the wealth estimation model. , and this for all operating points of the engine. The invention also makes it possible to perform a measurement every 6 ° of rotation of the crankshaft and thus to have a high frequency information of the measurement of richness, without falling into the measurement noise.
Ainsi, la présente invention concerne une méthode pour estimer la richesse en carburant dans chacun des cylindres d'un moteur à combustion interne comprenant un circuit d'échappement de gaz comprenant au moins des cylindres reliés à un collecteur et un capteur de mesure de la richesse (λ) en aval dudit collecteur. La méthode est caractérisée en ce qu'elle comporte les étapes suivantes :
- on établit un modèle physique représentant en temps réel l'éjection des gaz dans chacun desdits cylindres et leur parcours dans ledit circuit d'échappement jusqu'audit capteur;
- on définit une estimation de ladite richesse (λ) mesurée par ledit capteur à partir d'au moins une variable dudit modèle ;
- on couple ledit modèle avec un estimateur non linéaire de type adaptatif dans lequel on prend en compte ladite estimation de la mesure de richesse mesurée ;
- on réalise une estimation en temps réel de la valeur de la richesse dans chacun des cylindres, à partir dudit estimateur non linéaire de type adaptatif.
- establishing a physical model representing in real time the gas ejection in each of said cylinders and their course in said exhaust circuit to said sensor;
- an estimate of said richness (λ) measured by said sensor is defined from at least one variable of said model;
- said model is coupled with an adaptive-type non-linear estimator in which said estimate of the measure of measured wealth is taken into account;
- a real-time estimate of the value of the richness in each of the cylinders is made from said adaptive-type nonlinear estimator.
On peut également évaluer un temps de retard dû au temps de transit des gaz et au temps de réponse du capteur, en effectuant une perturbation test dans un cylindre déterminé et en mesurant son effet au capteur.It is also possible to evaluate a delay time due to the transit time of the gases and to the response time of the sensor, by performing a test perturbation in a given cylinder and measuring its effect on the sensor.
Selon un mode de réalisation, le modèle physique peut comporter au moins les trois types de variables suivants : la masse totale de gaz dans le collecteur d'échappement (M T ), la masse d'air frais dans le collecteur d'échappement (M air ) et les richesses dans chacun des cylindres (λ i ). Ce mode peut alors comporter au moins les deux types de données de sortie suivants : la masse totale de gaz dans le collecteur d'échappement (M T ) et des débits massiques sortant desdits cylindres (d i ).According to one embodiment, the physical model may comprise at least the following three types of variables: the total mass of gas in the exhaust manifold ( M T ), the fresh air mass in the exhaust manifold ( M air ) and the riches in each cylinder (λ i ). This mode can then comprise at least the two following types of output data: the total mass of gas in the exhaust manifold ( M T ) and the mass flow rates leaving said cylinders ( d i ).
La richesse (λ) mesurée peut être estimée en fonction de la masse totale de gaz dans le collecteur d'échappement (M T ) et la masse d'air frais dans le collecteur d'échappement (M air ).The measured richness (λ) can be estimated as a function of the total mass of gas in the exhaust manifold ( M T ) and the fresh air mass in the exhaust manifold ( M air ).
L'estimation de la valeur de la richesse dans chacun des cylindres peut alors comporter une correction en temps réel de l'estimation de la masse totale de gaz dans le collecteur d'échappement (M T ), de l'estimation de la masse d'air frais dans le collecteur d'échappement (M air ) et de l'estimation de la valeur de la richesse dans chacun des cylindres (λ i ).The estimation of the value of the richness in each of the cylinders may then comprise a real-time correction of the estimate of the total mass of gas in the exhaust manifold ( M T ), the estimation of the mass of the fresh air in the exhaust manifold ( M air ) and estimating the value of the richness in each of the cylinders (λ i ).
Enfin, la méthode peut être appliquée à un contrôle moteur pour adapter les masses de carburant injectées dans chacun des cylindres afin de régler la richesse dans tous les cylindres.Finally, the method can be applied to an engine control to adapt the fuel masses injected into each of the cylinders to adjust the richness in all the cylinders.
La présente invention sera mieux comprise et ses avantages apparaîtront plus clairement à la lecture de la description suivante d'un mode de réalisation, nullement limitatif, illustré par les figures ci-après annexées, parmi lesquelles :
- la figure 1 présente schématiquement les éléments descriptifs du processus d'échappement ;
- les figures 2A et 2B illustrent les richesses de références
- la figure 3 illustre la structure de l'estimateur ;
- les figures 4A et 4B illustrent les richesses de références
- Figure 1 schematically shows the descriptive elements of the exhaust process;
- FIGS. 2A and 2B illustrate the wealth of references
- Figure 3 illustrates the structure of the estimator;
- FIGS. 4A and 4B illustrate the wealth of references
Les intérêts d'une estimation de la richesse dans chacun des cylindres individuellement sont nombreux par rapport à une estimation de la richesse moyenne de l'ensemble des cylindres :
- gain sur le prix de revient si l'estimation est effectuée à partir d'une seule sonde de richesse en sortie turbine ;
- réduction des émissions polluantes ;
- amélioration de l'agrément de conduite (régularisation du couple délivré) ;
- réduction de la consommation de carburant
- diagnostic du système d'injection (détection de la dérive d'un injecteur ou de la défaillance du système d'injection).
- gain on the cost price if the estimate is made from a single turbine output wealth probe;
- reduction of polluting emissions;
- improvement of the driving pleasure (regularization of the delivered torque);
- reduced fuel consumption
- diagnosis of the injection system (detection of the drift of an injector or the failure of the injection system).
Le processus d'échappement comprend le parcours des gaz entre la soupape d'échappement jusqu'à l'air libre, en sortie du pot d'échappement. Le moteur de la présente mise en oeuvre en exemple est un 4 cylindres de 2200 cm3. Il est équipé d'un turbocompresseur à géométrie variable. Le schéma de la figure 1 présente les éléments descriptifs du processus d'échappement, dans la quelle :
- λ 1 à λ 4 représentent les richesses dans chacun des quatre cylindres ;
- SR représente la sonde de richesse ;
- CE correspond au collecteur d'échappement
- T correspond à la turbine du turbocompresseur ;
- DS1 à DS4 représentent les débits en sorties des cylindres.
- λ 1 to λ 4 represent the wealth in each of the four cylinders;
- SR represents the wealth probe;
- CE is the exhaust manifold
- T corresponds to the turbine of the turbocharger;
- DS1 to DS4 represent the output flows of the cylinders.
La sonde de richesse (SR) est située juste après la turbine (T). Les gaz, après la combustion dans le cylindre, subissent les actions suivantes :
- passage à travers la soupape d'échappement. Cette dernière étant commandée par un arbre à came, la loi de levée est en forme de cloche. Les débits passeront d'une valeur élevée, lors de l'ouverture de la soupape, à une valeur plus faible lorsque les pressions cylindre et collecteur s'égaliseront, pour enfin ré-augmenter lorsque le piston commencera à remonter pour éjecter les gaz d'échappement.
- passage dans une courte tubulure reliant le collecteur à la sortie de la culasse.
- phase de mélange dans le collecteur d'échappement (CE) où les débits (DS1 à DS4) des quatre cylindres se rejoignent. C'est ici que se produit le mélange des bouffées, fonction du type de collecteur (symétrique ou asymétrique), de l'AOE (Avance Ouverture Echappement) et du RFE (Retard Fermeture Echappement), qui détermineront la proportion de recouvrements des débits.
- passage à travers la turbine qui fournit le couple nécessaire au compresseur, situé en amont de l'admission. Bien que son action sur les débits soit peu connue, on peut penser qu'elle va mélanger un peu plus encore les bouffées provenant des différents cylindres.
- mesure par la sonde de type UEGO.
- passage through the exhaust valve. The latter being controlled by a camshaft, the law of emergence is bell-shaped. The flow rates will increase from a high value, when opening the valve, to a lower value when the cylinder and manifold pressures equalize, and finally increase when the piston starts to rise to eject the gas exhaust.
- passage in a short pipe connecting the manifold to the output of the cylinder head.
- mixing phase in the exhaust manifold ( CE ) where the flows ( DS1 to DS4) of the four cylinders meet. This is where the mixing of puffs occurs, depending on the type of collector (symmetrical or asymmetric), the AOE (Avance Opening Exhaust) and the RFE (Exhaust Closure Delay), which will determine the proportion of flow recoveries.
- passage through the turbine which provides the torque necessary for the compressor, located upstream of the intake. Although his action on flows is little known, we can think that it will mix a little more puffs from different cylinders.
- measurement by the UEGO type probe.
La composition des gaz d'échappement dépend des quantités de carburant et d'air introduites dans la chambre de combustion, de la composition du carburant et du développement de la combustion.The composition of the exhaust gas depends on the amount of fuel and air introduced into the combustion chamber, the fuel composition and the development of the combustion.
En pratique, la sonde de richesse mesure la concentration en O2 à l'intérieur d'une chambre de diffusion, mise en relation avec la conduite d'échappement par une barrière de diffusion fabriquée en matériaux poreux. Cette configuration pourra induire des différences suivant l'emplacement de la sonde choisit, notamment à cause des variations de températures et/ou de pressions à proximité de la sonde de richesse.In practice, the richness probe measures the concentration of O 2 inside a diffusion chamber, connected to the exhaust pipe by a diffusion barrier made of porous materials. This configuration may induce differences depending on the location of the chosen probe, in particular because of temperature variations and / or pressures in the vicinity of the richness probe.
Ce phénomène de variation de richesse dépendant de la pression ou de la température a cependant été négligé, puisque l'on s'intéresse à détecter des disparités de richesse entre les cylindres, la valeur moyenne étant normalement conservée par l'estimateur.This phenomenon of variation of richness depending on the pressure or the temperature has however been neglected, since one is interested in detecting disparities of richness between the cylinders, the average value being normally preserved by the estimator.
Dans le modèle physique temps réel utilisé par l'estimateur selon l'invention, on relie la richesse mesurée (λ) à la masse d'air (ou au débit d'air) se trouvant autour de la sonde et à la masse totale (ou au débit total). Le modèle est basé sur une approche trois gaz : air, carburant et gaz brûlés. Ainsi, on considère, qu'à mélange pauvre, la totalité du gaz restant après combustion est un mélange d'air et de gaz brûlés. Pour un mélange riche, le carburant étant en excès, on retrouve donc du carburant imbrûlé et des gaz brûlés après combustion, tandis que tout l'air a disparu. En réalité, la combustion n'est jamais à 100% complète, mais pour notre estimateur, la combustion sera considérée comme complète.In the real-time physical model used by the estimator according to the invention, the measured wealth (λ) is connected to the mass of air (or to the air flow) around the probe and to the total mass ( or at the total rate). The model is based on a three-gas approach: air, fuel and flue gas. Thus, it is considered that at lean mixture, all of the gas remaining after combustion is a mixture of air and flue gas. For a rich mixture, the fuel is in excess, so there is unburned fuel and burnt gases after combustion, while all the air is gone. In reality, the combustion is never 100% complete, but for our estimator, the combustion will be considered complete.
On définit une formulation reliant la richesse aux masses des trois espèces citées. Dans le cas d'un mélange pauvre : l'air est en excès, et il ne reste pas de carburant après la combustion. Avant combustion, on suppose que sont présentes les masses suivantes dans le cylindre, avec M air , la masse d'air, M carb , la masse de carburant et M gazB , la masse de gaz brûlés :
Sachant qu'il faut 14,7 fois plus d'air que de carburant pour être à la stoechiométrie, on peut construire ce tableau indiquant les masses de chaque espèce avant et après combustion :
La richesse λ représentant le rapport
Pour un mélange riche, la formule est la suivante :
Cependant, ces formules sont valables dans le cas où le mélange ne contient pas d'EGR, puisque la présence de gaz brûlés à l'admission modifie les concentrations des trois gaz à l'échappement.However, these formulas are valid in the case where the mixture does not contain EGR, since the presence of flue gases on admission modifies the concentrations of the three exhaust gases.
Dans le présent mode de réalisation, seule la formule de richesse pour mélange pauvre est utilisée dans l'estimateur, au niveau de l'intégration de la richesse dans l'équation (7), en négligeant une très faible partie de l'air (<3%). Cependant, l'invention ne se limite pas à ce mode, en effet, la formule est continue au voisinage d'une richesse égale à 1, et son inversion ne pose pas de problème pour les mélanges riches.In the present embodiment, only the lean mixture richness formula is used in the estimator, at the level of integration of richness in equation (7), neglecting a very small portion of the air ( <3%). However, the invention is not limited to this mode, in fact, the formula is continuous in the vicinity of a richness equal to 1, and its inversion does not pose a problem for rich mixtures.
Afin de mieux appréhender la manière dont se mélangent les gaz dans les conduites d'échappement, un modèle de moteur diesel a été utilisé sous le logiciel AMESim de la Société IMAGINE (France). Ce modèle, qui ne peut pas être inversé servira de référence pour valider le modèle selon l'invention.In order to gain a better understanding of the way gases mix in the exhaust pipes, a diesel engine model was used under the AMESim software of IMAGINE (France). This model, which can not be reversed, will serve as a reference for validating the model according to the invention.
AMESim est un logiciel de modélisation 0D, particulièrement bien adapté aux phénomènes thermiques et hydrauliques. Il permet notamment de modéliser des volumes, conduites ou restriction.AMESim is a 0D modeling software, particularly well suited to thermal and hydraulic phenomena. It allows to model volumes, behaviors or restrictions.
Le modèle d'échappement comprend :
- les tubulures d'échappement représentées par un volume et un tube ;
- le collecteur d'échappement avec échanges thermiques ;
- la turbine et la vanne de by-pass ;
- un volume à la confluence des débits turbine et vanne ;
- un tube entre la turbine et la sonde de mesure ;
- un volume et un tube pour la ligne d'échappement.
- the exhaust pipes represented by a volume and a tube;
- the exhaust manifold with heat exchange;
- the turbine and the bypass valve;
- a volume at the confluence of the turbine and valve flows;
- a tube between the turbine and the measuring probe;
- a volume and a tube for the exhaust line.
Les blocs élémentaires de modélisation des tubulures, restrictions et volumes sont décrits dans le manuel d'utilisation AMESim "Thermal Pneumatic Library". On utilise les équations standards pour calculer un débit à travers une restriction et les conservations d'énergie et de la masse. De plus, le modèle prend en compte les inerties des gaz, ce qui est important pour étudier la dynamique de la composition des gaz.The basic tubing, restriction and volume modeling blocks are described in the AMESim "Thermal Pneumatic Library" user manual. Standard equations are used to calculate a flow through a restriction and energy and mass conservations. In addition, the model takes into account gas inertia, which is important for studying the dynamics of gas composition.
Comme ce modèle est 0D, la dimension temps n'est pas prise en compte, et il n'est pas possible de modéliser un temps de retard avec une approche physique. Si une variable d'entrée est modifiée, la sortie est immédiatement changée. Le temps de transport est ainsi négligé. Cette limitation est importante lorsque l'on essaye de travailler sur des acquisitions en temps réel.Since this model is 0D, the time dimension is not taken into account, and it is not possible to model a delay time with a physical approach. If an input variable is changed, the output is immediately changed. The transport time is neglected. This limitation is important when trying to work on acquisitions in real time.
Selon l'invention, l'on définit un modèle physique temps réel unique pour modéliser le système global, c'est-à-dire l'ensemble du parcours des gaz d'échappement, depuis les cylindres jusqu'à l'échappement en aval de la turbine, en passant par le collecteur.According to the invention, a unique real-time physical model is defined for modeling the overall system, that is to say the entire path of the exhaust gases, from the cylinders to the downstream exhaust from the turbine, through the collector.
Dans le présent mode de réalisation, on considère que la variation de la température est faible sur un cycle moteur, et que son action est limitée sur les variations de débits. Ce sont en effet les variations de pressions qui sont primordiales dans le processus, puisque directement reliées aux débits. On s'impose donc une température fixe pour chaque élément : cylindres, collecteur et turbine. Les échanges de chaleur ne sont donc pas non plus modélisés. Cette hypothèse de simplification se révèle sans grande incidence.In the present embodiment, it is considered that the variation of the temperature is low on a motor cycle, and that its action is limited on the variations of flow rates. It is indeed the pressure variations that are essential in the process, since they are directly related to flow rates. We therefore impose a fixed temperature for each element: cylinders, collector and turbine. Heat exchanges are therefore not modeled either. This simplification hypothesis is not very important.
Dans une première approche, deux gaz sont considérés: de l'air frais et des gaz brûlés. Les équations classiques décrivent l'évolution de la masse totale des gaz dans les volumes, et de la masse de l'air frais. Les gaz brûlés peuvent ensuite en être déduit. Cette démarche est valable dans le cas d'un fonctionnement en mélange pauvre, mais des équations similaires peuvent être écrites pour le carburant et les gaz brûlés, dans le cas de mélange riche.In a first approach, two gases are considered: fresh air and flue gases. The classical equations describe the evolution of the total mass of gases in volumes, and the mass of fresh air. The flue gas can then be deduced. This approach is valid in the case of a lean mixture operation, but similar equations can be written for the fuel and the flue gases, in the case of rich mixture.
Le collecteur d'échappement est modélisé selon un volume dans lequel il y a conservation de la masse. On suppose que la température est sensiblement constante, et déterminée à partir d'abaque fonction de la charge et du régime moteur.The exhaust manifold is modeled according to a volume in which there is conservation of the mass. It is assumed that the temperature is substantially constant, and determined from an abacus function of the load and the engine speed.
Selon l'invention, on a choisi de relier la richesse mesurée à la masse d'air se trouvant autour de la sonde et à la masse totale. Ainsi, la conservation de la masse totale dans le collecteur exprime le fait que la masse de gaz d'échappement dans le collecteur est égal à la masse de gaz d'échappement entrante dans le collecteur (débit de sortie des cylindres) diminuée de la masse sortante du collecteur. On suppose que la composition du débit dans la turbine est la même qu'à la sortie du collecteur. Ainsi la masse sortante du collecteur est égale au débit passant par la turbine. On a ainsi la formule suivante pour la masse totale :
avec :
- Ne
- : Régime moteur
- α :
- Angle du vilebrequin
- M T
- : Masse totale dans le collecteur d'échappement
- d i :
- Débit massique sortant du cylindre i
- d T :
- Débit total passant par la turbine
with:
- Born
- : Engine speed
- α :
- Angle of the crankshaft
- M T
- : Total mass in the exhaust manifold
- d i :
- Mass flow coming out of cylinder i
- d T :
- Total flow through the turbine
De même pour la conservation de la masse d'air, nous avons :
avec
- Ne :
- Régime moteur
- α :
- Angle du vilebrequin
- M air :
- Masse d'air frais dans le collecteur d'échappement
- λ i :
- Richesse dans chacun des cylindres
- d i :
- Débit massique sortant du cylindre i
- d air :
- Débit d'air passant par la turbine
with
- Do not :
- Engine speed
- α :
- Angle of the crankshaft
- M air :
- Fresh air mass in the exhaust manifold
- λ i :
- Wealth in each of the cylinders
- d i :
- Mass flow coming out of cylinder i
- of air:
- Air flow through the turbine
On décrit maintenant les modèles physiques permettant de déterminer le débit de sortie des cylindres et les débits passant par la turbine.The physical models for determining cylinder output and flow through the turbine are now described.
Le débit de sortie des gaz à la sortie des cylindres peut être modéliser grâce un modèle physique décrivant le débit en sortie des soupapes d'échappement. Pour ce modèle, d'éjection des gaz par les soupapes, on utilise trois variables :
- l'angle du vilebrequin (α);
- le débit aspiré par le cylindre d asp (variable estimée par le contrôle moteur en amont) ;
- la valeur moyenne de richesse mesurée par la sonde λ̅ sur un cycle.
- the angle of the crankshaft ( α );
- the flow rate sucked by the cylinder asp (variable estimated by the engine control upstream);
- the mean value of richness measured by the probe λ̅ on a cycle.
Le débit sortant moyen est connu à partir du débit aspiré et du débit d'essence injectée. La valeur instantanée du débit sortant, quant à elle, est basée sur un gabarit dépendant du débit aspiré. Ce gabarit est un modèle physique (une courbe) reposant sur une loi empirique permettant d'estimer un débit moyen pour un cylindre en fonction de l'angle de vilebrequin à partir du régime moteur, de l'angle du vilebrequin, du débit aspiré par le cylindre et de la valeur moyenne de richesse mesurée par la sonde sur un cycle. Cette loi physique a pour seule contrainte de respecter le débit sortant moyen (aire de la courbe) et de fournir une courbe rendant compte des deux phénomènes suivants :
- l'équilibre de pression cylindre/échappement qui se traduit par un pic de débit en fonction de l'angle du vilebrequin ;
- un débit qui se calque sur la section efficace de la soupape d'échappement qui se traduit par un second pic de débit d'amplitude plus faible.
- the cylinder / exhaust pressure equilibrium which results in a flow peak as a function of crankshaft angle;
- a flow which is traced on the effective section of the exhaust valve which results in a second peak of lower amplitude flow rate.
Ce gabarit (d̅) fournit en sortie une allure (courbe) du débit massique en sortie des soupapes d'échappement d i , c'est-à-dire une estimation commune (d̅) du débit pour tous les cylindres. Il est obtenu en corrélation avec les mesures de banc moteur. En fonction du débit aspiré et de la richesse moyenne mesurée par la sonde λ̅, on effectue ensuite une homothétie du gabarit et l'on effectue un déphasage pour chaque cylindre en fonction de l'angle de vilebrequin, pour déduire le débit de sortie des gaz à la sortie de chaque cylindre :
avec :
- d i (α) : le débit de sortie des gaz à la sortie du cylindre i ;
- d̅(α) : le gabarit, c'est-à-dire une estimation du débit en sortie des cylindres ;
- α i : l'angle de déphasage pour le cylindre i.
with:
- d i ( α ): the output flow rate of the gases at the outlet of the cylinder i;
- d̅ ( α ): the template, that is to say an estimate of the output flow of the cylinders;
- α i : the phase shift angle for the cylinder i.
On peut observer schématiquement le déphasage de la courbe du gabarit sur la figure 1 (DS1 à DS4). It is possible to observe schematically the phase shift of the curve of the template in FIG. 1 ( DS1 to DS4 ) .
On décrit maintenant les modèles physiques permettant de déterminer le débit passant par la turbine.The physical models for determining the flow through the turbine are now described.
La turbine est modélisée selon un débit passant à travers une restriction du débit. Le débit dans la turbine est généralement donné par cartographie (abaque) en fonction du régime de la turbine et du rapport de pression amont/aval de la turbine.The turbine is modeled according to a flow passing through a flow restriction. The flow rate in the turbine is generally given by mapping (abacus) as a function of the turbine speed and the upstream / downstream pressure ratio of the turbine.
Le débit passant par la turbine d T est une fonction de la masse totale (M T ) dans le collecteur d'échappement, de la température dans le collecteur d'échappement, du régime du turbocompresseur et de la géométrie du turbocompresseur. Les données d'entrée de ce modèle sont donc :
- La masse totale de gaz d'échappement (M T ) ;
- La masse d'air (M air ) ;
- Le régime moteur (N e ) ;
- Le régime de la turbine (du turbocompresseur).
- The total mass of exhaust gas ( M T );
- The air mass ( M air );
- The engine speed ( N e );
- The speed of the turbine (turbocharger).
Ce débit peut être estimé à partir d'une fonction concave de la masse totale M T . Cette fonction est notée p. Le débit dans la turbine s'écrit alors : d T (M T )=M T .p(M T ).This flow rate can be estimated from a concave function of the total mass MT. This function is noted p. The flow rate in the turbine is then written: d T ( M T ) = M T .p ( M T ).
La fonction p est une fonction de type racine qui s'exprime en fonction d'une part du régime de la turbine, et d'autre part en fonction du rapport entre la masse totale dans le collecteur d'échappement (M T ) et la masse dans le collecteur dans les conditions atmosphériques (M 0). Ainsi la cartographie donne p(M T ) en fonction du rapport
où :
- f est une fonction polynomiale
- g est une constante
or :
- f is a polynomial function
- g is a constant
Les paramètres de la fonction f sont optimisés par corrélation avec la cartographie de la turbine.The parameters of the function f are optimized by correlation with the mapping of the turbine.
De plus la composition en air est supposée la même que dans le collecteur d'échappement. Le débit d'air passant à travers la turbine est donc :
Ainsi, à l'aide des modèles physiques d'éjection des gaz et de la turbine, les équations (1) et (2) s'écrivent :
Ce système d'équation (3) constitue le modèle physique du collecteur d'échappement.
Les données d'entrées de ce modèle sont :
- Ne :
- Régime moteur
- α :
- Angle du vilebrequin
- d i :
- Débit massique sortant du cylindre i
- d T :
- Débit total passant par la turbine
- d air :
- Débit d'air passant par la turbine
- λ i :
- Richesse dans chacun des cylindres
- M T :
- Masse totale dans le collecteur d'échappement
- M air :
- Masse d'air frais dans le collecteur d'échappement
The input data of this model are:
- Do not :
- Engine speed
- α :
- Angle of the crankshaft
- d i :
- Mass flow coming out of cylinder i
- d T :
- Total flow through the turbine
- of air:
- Air flow through the turbine
- λ i :
- Wealth in each of the cylinders
- M T :
- Total mass in the exhaust manifold
- M air :
- Fresh air mass in the exhaust manifold
La première équation contient une inconnue : M T . La seconde en contient deux : M air et λ i . Ceci conduit aux hypothèses supplémentaires décrites ci-après.The first equation contains one unknown: M T. The second contains two: M air and λ i . This leads to the additional assumptions described below.
Pour compléter le modèle physique temps réel (RTM) du collecteur d'échappement, on suppose que les richesses sorties cylindres sont constantes sur un point de fonctionnement, on a donc :
En fait, le calcul s'effectuant en temps réel, l'on estime des constantes λ i . To complete the real-time physical model (RTM) of the exhaust manifold, it is assumed that the cylindrical outputs are constant on an operating point, so we have:
In fact, the computation taking place in real time, one estimates constant λ i .
Finalement, le modèle physique temps réel RTM peut se mettre sous la forme matricielle suivante :
Les inconnues du modèle physique sont au final M T , M air et les λ i .
Les données de sorties du modèle physique sont M T et d i . Finally, the RTM real-time physical model can be in the following matrix form:
The unknowns of the physical model are ultimately M T , M air and λ i .
The output data of the physical model is M T and d i .
Les temps de retard dus au transport du gaz dans les tubulures et les différents volumes, ainsi que du temps « mort » de la sonde de mesure, n'est pas pris en compte dans le modèle physique décrit ci-dessus (système d'équations 5). Cependant le modèle a été construit de manière linéaire vis-à-vis de ce ces retards car on néglige le transport dans les tubulures. Aussi, les retards peuvent être compilés en un seul délai de retard pour l'ensemble du processus d'échappement, et le modèle physique peut être inversé tel quel, l'influence du délai de retard pouvant être considéré ultérieurement, comme explicité ci-après.The delay times due to the gas transport in the pipes and the different volumes, as well as the "dead" time of the measuring probe, are not taken into account in the physical model described above (system of equations 5). However, model was built linearly with respect to these delays because we neglect the transport in the pipes. Also, the delays can be compiled in a single delay time for the entire escape process, and the physical model can be reversed as is, the influence of the delay delay can be considered later, as explained below. .
Le modèle physique (5) décrit ci-dessus que la richesse en aval de la turbine (que l'on considère identique à la richesse dans le collecteur) est exprimée comme une fonction de la composition du débit de gaz à l'entrée du collecteur d'échappement.The physical model (5) described above that the wealth downstream of the turbine (which is considered identical to the richness in the collector) is expressed as a function of the composition of the gas flow at the inlet of the collector exhaust.
Les données mesurées sont :
- Richesse mesurée par la sonde: λ
- Wealth measured by the probe: λ
Les autres données connues du système sont :
- Régime moteur : Ne
- Angle du vilebrequin : α
- Régime de la turbine (du turbocompresseur)
- Débit aspiré par le cylindre
- Engine Speed: Ne
- Angle of the crankshaft: α
- Turbocharger (Turbocharger) Regime
- Flow sucked by the cylinder
Les données modélisées du système sont :
- Débit massique sortant du cylindre i : d i
- Débit total passant par la turbine : d T
- Débit d'air passant par la turbine : d air
- Masse totale dans le collecteur d'échappement : M T
- Mass flow rate out of cylinder i : d i
- Total flow through the turbine: d T
- Air passing through the turbine flow rate: of air
- Total mass in the exhaust manifold: M T
Les inconnues sont donc :
- Richesses dans chacun des quatre cylindres : λ i
- Masse d'air frais dans le collecteur d'échappement : M air
- Riches in each of the four cylinders: λ i
- Fresh air mass in the exhaust manifold: M air
Le modèle physique (5) est non linéaire, et il est impossible de résoudre un tel système en temps réel. Il est donc nécessaire d'utiliser un estimateur, plutôt que de chercher à calculer directement les inconnues du système. Le choix de l'estimateur selon l'invention repose sur le fait que la structure du système est linéaire en fonction des richesses dans les cylindres λ i (la variation de masse d'air est linéaire en fonction des λ i ). Dans ce cadre, une technique particulièrement bien adaptée consiste à utiliser un filtre adaptatif. Ainsi, pour estimer les inconnues à partir du modèle physique RTM, la méthode selon l'invention propose de construire un estimateur basé sur un filtre adaptatif. Cet estimateur permet, au final, une estimation de la richesse cylindre à cylindre à partir de la mesure de richesse par le capteur situé derrière la turbine.The physical model (5) is non-linear, and it is impossible to solve such a system in real time. It is therefore necessary to use an estimator, rather than seeking to directly calculate the unknowns of the system. The choice of the estimator according to the invention is based on the fact that the structure of the system is linear as a function of the wealth in the cylinders λ i (the air mass variation is linear as a function of λ i ). In this context, a particularly suitable technique is to use an adaptive filter. Thus, to estimate the unknowns from the physical model RTM, the method according to the invention proposes to construct an estimator based on an adaptive filter. This estimator ultimately allows an estimation of the cylinder to cylinder richness from the measurement of wealth by the sensor located behind the turbine.
De façon générale les filtres adaptatifs sont des systèmes appliqués sur des données bruitées pour obtenir une information utile à un certain instant t, ces systèmes étant mis en oeuvre dans trois configurations :
- le filtrage, c'est-à-dire l'extraction de l'information utile au temps t à partir des données bruitées mesurées jusqu'au temps t inclus ;
- le lissage, qui utilisera aussi les données postérieures au temps t ;
- la prédiction, qui ne se sert que des données jusqu'au temps t-τ pour déduire l'information qui nous intéresse au temps t.
- the filtering, that is to say the extraction of the useful information at time t from the noisy data measured up to the time t included;
- smoothing, which will also use the data after time t;
- prediction, which only uses data up to time t - τ to derive the information of interest at time t.
Dans notre cas, le but est de reconstruire, à partir des deux éléments considérés comme des mesures y 1 = M T et y 2 = λ, les données suivantes :
- M T :
- Masse totale dans le collecteur d'échappement
- M air :
- Masse d'air frais dans le collecteur d'échappement
- λ i :
- Richesse dans chacun des cylindres
λ est mesuré par la sonde, et M T est estimé à partir du modèle RTM (5).In our case, the goal is to reconstruct, from the two elements considered as measures y 1 = M T and y 2 = λ , the following data:
- M T :
- Total mass in the exhaust manifold
- M air :
- Fresh air mass in the exhaust manifold
- λ i :
- Wealth in each of the cylinders
λ is measured by the probe, and M T is estimated from the RTM model (5).
En posant :
- M̂ T :
- Estimateur de la masse totale dans le collecteur d'échappement
- M̂ air :
- Estimateur de la masse d'air frais dans le collecteur d'échappement
- λ̂ i :
- Estimateur de la richesse dans chacun des cylindres
En d'autres termes, il s'agît du modèle RTM (5) dans lequel on applique un terme correctif pour chaque paramètre à estimer : C MT , C Mair , C λi . By asking :
- M T :
- Estimator of the total mass in the exhaust manifold
- M air :
- Estimator of the fresh air mass in the exhaust manifold
- λ i :
- Estimator of wealth in each of the cylinders
In other words, it is the RTM model (5) in which a correction term is applied for each parameter to be estimated: C MT , C Mair , C λi .
Le principe de l'estimateur est de faire converger le modèle physique (5), et par conséquent les richesses λ i vers la réalité. En effet, le modèle (5) fournit en sortie M T et M Air , et l'on a par ailleurs les paramètres d'entrée Y. L'estimateur compare donc les valeurs de sorties du modèle RTM avec les valeurs d'entrée, puis effectue les corrections adéquates. Par exemple, les richesses λ i doivent s'adapter en fonction de l'erreur sur M T et λ : si l'erreur entre les valeurs d'entrée, M T et λ , et les valeurs estimées correspondantes,M̂ T etλ̂ , est négative, alors il faut augmenter les valeurs estimées et réciproquement. C'est pourquoi on écrit :
En appliquant le même principe pour les corrections de M T et M air , on obtient :
où L 1, L 2, L λ sont des paramètres de réglages, permettant de contrôler la vitesse de convergence de la solution aux trois inconnues. Ce sont des paramètres réels strictement positifs. Ces paramètres sont réglés manuellement afin d'obtenir un bon compromis entre la vitesse de convergence et la faible sensibilité au bruit de mesure.
Enfin, en utilisant l'équation de la richesse précédemment décrite, il est possible d'estimer la richesse mesurée à partir de M T et M air :
On en déduit : M air = (1-λ)M T
Et en utilisant les notations y 1 =M T et y 2 = λ, l'estimateur du modèle physique RTM, basé sur le filtre adaptatif et sur l'estimation de la richesse à partir de M T et M air , s'écrit alors :
avec δ T = 1The principle of the estimator is to converge the physical model (5), and consequently the riches λ i towards reality. Indeed, the model (5) outputs M T and M Air , and we also have input parameters Y. The estimator therefore compares the output values of the RTM model with the input values, then make the appropriate corrections. For example, the wealth λ i must adapt according to the error on M T and λ : if the error between the input values, M T and λ , and the corresponding estimated values, M T and λ , is negative, then the estimated values must be increased and vice versa. This is why we write:
Applying the same principle for the corrections of M T and M air , we obtain:
where L 1 , L 2 , L λ are adjustment parameters, making it possible to control the speed of convergence of the solution to the three unknowns. These are strictly positive real parameters. These parameters are set manually to obtain a good compromise between the speed of convergence and the low sensitivity to measurement noise.
Finally, using the previously described wealth equation, it is possible to estimate the wealth measured from M T and M air :
We deduce: M air = (1 -λ ) M T
And by using the notations y 1 = M T and y 2 = λ , the physical model RTM estimator, based on the adaptive filter and the estimation of the richness from M T and M air , is written :
with δ T = 1
L'estimateur ainsi construit, permet de corriger en temps réel M T , M air et λ, à partir d'une première valeur de M T fournie par le modèle RTM et à partir de la mesure de richesse effectuée par la sonde.The estimator thus constructed makes it possible to correct in real time M T , M air and λ , from a first value of M T provided by the RTM model and from the measurement of richness made by the probe.
On résout numériquement en temps réel le système (8), le calculateur faisant appel à une discrétisation d'Euler explicite, bien connu de l'homme du métier.The system (8) is numerically solved in real time, the calculator using an explicit Euler discretization, well known to those skilled in the art.
A partir des richesses individuelles connues, on estime par la modélisation de référence AMESim une richesse au niveau de la sonde (λ). Cette valeur de richesse (λ) est utilisée en entrée de l'estimateur. La dynamique de la sonde n'a pas été prise en compte. On applique des déséquilibres d'injection et on observe l'estimation cylindre à cylindre de la richesse (λ i ) à partir de la mesure de richesse derrière la turbine (λ).From the known individual wealth, one estimates by reference modelization AMESim a richness on the level of the probe ( λ ). This wealth value (λ) is used at the input of the estimator. The dynamics of the probe has not been taken into account. Injection imbalances are applied and the cylinder-to-cylinder estimation of the richness ( λ i ) is observed from the wealth measurement behind the turbine ( λ ).
Pour la simulation, on déséquilibre successivement les 4 cylindres en introduisant 80µs d'injection en plus sur le cylindre, puis on déséquilibre le cylindre 1 et 4 de la même façon. Les figures 2A et 2B montrent en bas les richesses de références
L'estimateur implémenté comme décrit plus haut, ne permet pas à la méthode d'estimation de prendre en compte le temps de retard entre l'échappement cylindre et le signal acquis par la sonde. Dans la réalité, le temps de retard provient de plusieurs sources : temps de transport dans les tubulures et à travers les volumes, temps mort de la sonde de mesure.The estimator implemented as described above, does not allow the estimation method to take into account the delay time between the cylinder exhaust and the signal acquired by the probe. In reality, the delay time comes from several sources: transport time in the pipes and through the volumes, dead time of the measuring probe.
En appliquant un temps de retard D à l'entrée de l'estimateur sur les variables provenant du modèle, on peut synchroniser l'estimateur avec les mesures de richesse. La structure de l'estimateur avec retard est illustrée sur la figure 3, dans laquelle :
- N e et α sont les données d'entrées du modèle RTM décrit par les équations (5) ;
- MMBO est le Modèle de Masse Boucle Ouverte (modèle RTM) ;
- D est le retard appliqué aux variables de sorties du modèle RTM (MMBO) ; ce retard est issu de l'équation (9) ;
- SR est la sonde de mesure de la richesse en aval de la turbine utilisée dans l'estimateur via l'équation (7) ;
- ERFA est l'Estimateur de Richesse basé sur un Filtre Adaptatif et décrit par l'équation (8) ;
- λ i est la richesse du cylindre i estimée par l'ERFA.
- N e and α are the input data of the RTM model described by equations (5);
- MMBO is the Open Loop Mass Model (RTM model);
- D is the delay applied to the output variables of the RTM model ( MMBO ); this delay comes from equation (9);
- SR is the measurement of the wealth downstream of the turbine used in the estimator via equation (7);
- ERFA is the Adaptive Filter-based Wealth Estimator and described by Equation (8);
- λ i is the richness of the cylinder i estimated by the ERFA.
Le temps de retard dépend des conditions de fonctionnement : vitesse du moteur, charge, pression au collecteur d'échappement, etc. Comme le retard est difficile à modéliser, on a développé une méthode d'identification, de manière à calculer en temps réel le retard entre l'estimateur et les mesures sans appel à une instrumentation supplémentaire. Le principe consiste à appliquer un petit échelon dans le voisinage du point d'injection du cylindre 1, et à calculer les variations estimées de richesse pour chacun des cylindres. Puis, un critère d'identification J k est construit de manière à pénaliser les variations des cylindres 2, 3, et 4.
avec :
with:
La pénalisation est donnée par β. S'il y a une variation positive de la valeur de richesse estimée pour le cylindre 2, alors le temps de retard entre l'estimateur et les mesures est positif. S'il y a une variation sur le cylindre 3, le retard est négatif et la pénalisation est négative. Une variation du cylindre 4 peut être considérée comme une conséquence d'un retard positif ou négatif. Le retard D appliqué aux variables de sorties du modèle RTM, est un délai additif, il est calculé par moindre carré en minimisant J k . The penalty is given by β . If there is a positive variation in the estimated wealth value for
Le critère J k est contrôlé à zéro par un contrôleur PI (Proportionnel Intégral) sur le retard de l'estimateur. Quand le contrôleur est stabilisé, la variation de la richesse estimée est maximale sur le cylindre 1, et minimale sur le cylindre 4. L'estimateur est alors en phase avec les mesures.The criterion J k is controlled to zero by a PI (Proportional Integral) controller on the delay of the estimator. When the controller is stabilized, the variation of the estimated richness is maximum on the
Les figures 4A et 4B illustrent l'estimation de la richesse cylindre à cylindre par l'estimateur décrit précédemment à 1500tr/min moyenne charge. Ces figures montrent en haut les richesses de références
La présente invention concerne une méthode d'estimation comportant la construction d'un estimateur, permettant, à partir de la mesure de richesse de la sonde (λ) et de l'information de masse totale de gaz à l'intérieur du collecteur (M T ), d'estimer les richesses en sortie des quatre cylindres (λ i ). L'estimateur ainsi réalisé est performant, et surtout ne nécessite aucun réglage supplémentaire dans le cas de changement du point de fonctionnement. Aucune phase d'identification n'est nécessaire, seul un réglage des bruits de mesure et de modèle doit être effectué, une seule et unique fois.The present invention relates to an estimation method comprising the construction of an estimator, making it possible, from the measurement of the richness of the probe (λ) and the total mass of gas information inside the collector ( M T ), to estimate the wealth at the output of the four cylinders (λ i ). The estimator thus produced is efficient, and above all does not require any additional adjustment in the case of change of the operating point. No identification phase is necessary, only a measurement and model noise adjustment must be made once and only once.
Pour rendre plus robuste l'estimation selon l'invention, quelles que soient les conditions de fonctionnement, un contrôleur de temps de retard est mis en parallèle de l'estimateur, permettant de recaler le temps de retard à la suite d'un échelon de temps d'injection sur un cylindre. Ceci permet un calage optimal de l'estimateur, par exemple avant une phase de richesse égale à 1.To make the estimation according to the invention more robust, whatever the operating conditions, a delay time controller is put in parallel with the estimator, making it possible to reset the delay time following a step of injection time on a cylinder. This allows optimal calibration of the estimator, for example before a rich phase equal to 1.
L'invention permet en outre d'effectuer une mesure tous les 6° de rotation du vilebrequin et donc d'avoir une information haute fréquence de la mesure de richesse, sans pour autant tomber dans le bruit de mesure. De plus, la représentation haute fréquence permet de tenir compte de l'effet pulsatoire du système. Le système modélisé est périodique et permet d'obtenir un estimateur avec une meilleure dynamique : on anticipe la pulsation de l'échappement.The invention also makes it possible to perform a measurement every 6 ° of rotation of the crankshaft and thus to have a high frequency information of the measurement of richness, without falling into the measurement noise. In addition, the high frequency representation makes it possible to take into account the pulsating effect of the system. The modeled system is periodic and makes it possible to obtain an estimator with a better dynamics: one anticipates the pulsation of the escapement.
Par ailleurs l'invention permet de réduire le temps de calcul d'un facteur 80 environ par rapport aux méthodes antérieures.Moreover, the invention makes it possible to reduce the calculation time by a factor of about 80 compared to the previous methods.
Claims (7)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0505442A FR2886345B1 (en) | 2005-05-30 | 2005-05-30 | METHOD OF ESTIMATING AN ADAPTIVE NON-LINEAR FILTER OF WEALTH IN A CYLINDER OF A COMBUSTION ENGINE |
Publications (2)
Publication Number | Publication Date |
---|---|
EP1729001A1 true EP1729001A1 (en) | 2006-12-06 |
EP1729001B1 EP1729001B1 (en) | 2008-03-26 |
Family
ID=35058371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP06290558A Expired - Fee Related EP1729001B1 (en) | 2005-05-30 | 2006-04-03 | Estimation method using a non-linear adaptive filter of air-fuel mixture richness in a cylinder of a combustion engine |
Country Status (5)
Country | Link |
---|---|
US (1) | US7483782B2 (en) |
EP (1) | EP1729001B1 (en) |
JP (1) | JP4964503B2 (en) |
DE (1) | DE602006000790T2 (en) |
FR (1) | FR2886345B1 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008135312A1 (en) * | 2007-05-07 | 2008-11-13 | Continental Automotive Gmbh | Method and device for determining the combustion lambda value of an internal combustion engine |
WO2009011191A1 (en) * | 2007-07-19 | 2009-01-22 | Toyota Jidosha Kabushiki Kaisha | Abnormality detection device for internal combustion engine and air/fuel ratio control apparatus for internal combustion engine |
FR2929650A1 (en) * | 2008-04-04 | 2009-10-09 | Bosch Gmbh Robert | METHOD AND DEVICE FOR ADAPTING A DYNAMIC MODEL OF AN EXHAUST GAS PROBE. |
EP2687709A4 (en) * | 2011-03-16 | 2017-01-11 | Toyota Jidosha Kabushiki Kaisha | Internal combustion engine control apparatus |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2886346B1 (en) * | 2005-05-30 | 2010-08-27 | Inst Francais Du Petrole | METHOD OF ESTIMATING EXTENDED KALMAN FILTER OF WEALTH IN A CYLINDER OF A COMBUSTION ENGINE |
DE102005057975A1 (en) * | 2005-12-05 | 2007-06-06 | Robert Bosch Gmbh | Method for controlling fuel or air flow to individual cylinder of internal-combustion engine, involves evaluation of signal which is influenced by combustion or affects value which has influence on combustion |
ES2758794T3 (en) * | 2009-12-23 | 2020-05-06 | Fpt Motorenforschung Ag | Method and apparatus for measuring and controlling the EGR rate in a combustion engine |
US7987840B2 (en) * | 2010-04-14 | 2011-08-02 | Ford Global Technologies, Llc | Delay compensated air/fuel control of an internal combustion engine of a vehicle |
JP2013253593A (en) * | 2012-05-11 | 2013-12-19 | Denso Corp | Cylinder-by-cylinder air fuel ratio control device for internal combustion engine |
US9291526B2 (en) | 2012-11-12 | 2016-03-22 | Kerdea Technologies, Inc. | Oxygen sensing method and system |
US10030593B2 (en) * | 2014-05-29 | 2018-07-24 | Cummins Inc. | System and method for detecting air fuel ratio imbalance |
JP6800799B2 (en) * | 2017-04-05 | 2020-12-16 | オムロン株式会社 | Control devices, control programs, control systems, and control methods |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0688945A2 (en) * | 1994-06-20 | 1995-12-27 | Honda Giken Kogyo Kabushiki Kaisha | Air/fuel ratio detection system for multicylinder internal combustion engine |
EP0724073A2 (en) * | 1995-01-27 | 1996-07-31 | Matsushita Electric Industrial Co., Ltd. | Air-fuel ratio control system |
US5839415A (en) * | 1995-02-24 | 1998-11-24 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio control system having function of after-start lean-burn control for internal combustion engines |
US5911682A (en) * | 1996-08-29 | 1999-06-15 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio control system for internal combustion engines |
FR2834314A1 (en) | 2001-12-31 | 2003-07-04 | Peugeot Citroen Automobiles Sa | Procedure for estimating fuel richness of mixture burned in i.c. engine with fuel injection uses Kalman estimator and set equations |
US20050022797A1 (en) * | 2003-07-30 | 2005-02-03 | Denso Corporation | Cylinder-by-cylinder air-fuel ratio calculation apparatus for multi-cylinder internal combustion engine |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4903668A (en) * | 1987-07-29 | 1990-02-27 | Toyota Jidosha Kabushiki Kaisha | Fuel injection system of an internal combustion engine |
JP2705165B2 (en) * | 1988-12-09 | 1998-01-26 | トヨタ自動車株式会社 | Fuel injection amount control device for multi-cylinder internal combustion engine |
US5535135A (en) * | 1993-08-24 | 1996-07-09 | Motorola, Inc. | State estimator based exhaust gas chemistry measurement system and method |
JP3149714B2 (en) * | 1994-12-28 | 2001-03-26 | 日産自動車株式会社 | Catalyst deterioration diagnosis device for internal combustion engine |
JP3749971B2 (en) * | 1995-02-27 | 2006-03-01 | 本田技研工業株式会社 | Fuel injection control device for internal combustion engine |
JP3340058B2 (en) * | 1997-08-29 | 2002-10-28 | 本田技研工業株式会社 | Air-fuel ratio control system for multi-cylinder engine |
JP3304844B2 (en) * | 1997-08-29 | 2002-07-22 | 本田技研工業株式会社 | Plant control equipment |
JP3304845B2 (en) * | 1997-08-29 | 2002-07-22 | 本田技研工業株式会社 | Plant control equipment |
JPH1185719A (en) * | 1997-09-03 | 1999-03-30 | Matsushita Electric Ind Co Ltd | Parameter estimation device |
FR2773847B1 (en) * | 1998-01-19 | 2000-03-24 | Sagem | INJECTION SYSTEM RICHNESS ESTIMATING DEVICE FOR INTERNAL COMBUSTION ENGINE |
JP3842528B2 (en) * | 2000-07-07 | 2006-11-08 | 株式会社日立製作所 | Air-fuel ratio feedback control device for internal combustion engine |
JP3960106B2 (en) * | 2001-07-12 | 2007-08-15 | 日産自動車株式会社 | Engine air-fuel ratio control device |
US6882929B2 (en) * | 2002-05-15 | 2005-04-19 | Caterpillar Inc | NOx emission-control system using a virtual sensor |
JP4144405B2 (en) * | 2003-04-02 | 2008-09-03 | 日産自動車株式会社 | Deterioration judgment device for exhaust aftertreatment device |
JP4357863B2 (en) * | 2003-04-14 | 2009-11-04 | 株式会社デンソー | Multi-cylinder internal combustion engine cylinder-by-cylinder air-fuel ratio calculation device |
FR2867232B1 (en) | 2004-03-05 | 2006-05-05 | Inst Francais Du Petrole | METHOD OF ESTIMATING FUEL WEALTH IN A CYLINDER OF A COMBUSTION ENGINE |
US7015725B1 (en) * | 2004-12-06 | 2006-03-21 | Faraday Technology Corp. | Delay-locked loop device capable of anti-false-locking |
JP4420288B2 (en) * | 2005-04-25 | 2010-02-24 | 株式会社デンソー | Cylinder-by-cylinder air-fuel ratio control apparatus for internal combustion engine |
FR2886346B1 (en) | 2005-05-30 | 2010-08-27 | Inst Francais Du Petrole | METHOD OF ESTIMATING EXTENDED KALMAN FILTER OF WEALTH IN A CYLINDER OF A COMBUSTION ENGINE |
-
2005
- 2005-05-30 FR FR0505442A patent/FR2886345B1/en not_active Expired - Fee Related
-
2006
- 2006-04-03 DE DE602006000790T patent/DE602006000790T2/en active Active
- 2006-04-03 EP EP06290558A patent/EP1729001B1/en not_active Expired - Fee Related
- 2006-05-22 US US11/437,702 patent/US7483782B2/en not_active Expired - Fee Related
- 2006-05-30 JP JP2006149259A patent/JP4964503B2/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0688945A2 (en) * | 1994-06-20 | 1995-12-27 | Honda Giken Kogyo Kabushiki Kaisha | Air/fuel ratio detection system for multicylinder internal combustion engine |
EP0724073A2 (en) * | 1995-01-27 | 1996-07-31 | Matsushita Electric Industrial Co., Ltd. | Air-fuel ratio control system |
US5839415A (en) * | 1995-02-24 | 1998-11-24 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio control system having function of after-start lean-burn control for internal combustion engines |
US5911682A (en) * | 1996-08-29 | 1999-06-15 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio control system for internal combustion engines |
FR2834314A1 (en) | 2001-12-31 | 2003-07-04 | Peugeot Citroen Automobiles Sa | Procedure for estimating fuel richness of mixture burned in i.c. engine with fuel injection uses Kalman estimator and set equations |
US20050022797A1 (en) * | 2003-07-30 | 2005-02-03 | Denso Corporation | Cylinder-by-cylinder air-fuel ratio calculation apparatus for multi-cylinder internal combustion engine |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008135312A1 (en) * | 2007-05-07 | 2008-11-13 | Continental Automotive Gmbh | Method and device for determining the combustion lambda value of an internal combustion engine |
US8364377B2 (en) | 2007-05-07 | 2013-01-29 | Continental Automotive Gmbh | Method and device for determining the combustion lambda value of an internal combustion engine |
WO2009011191A1 (en) * | 2007-07-19 | 2009-01-22 | Toyota Jidosha Kabushiki Kaisha | Abnormality detection device for internal combustion engine and air/fuel ratio control apparatus for internal combustion engine |
US8050852B2 (en) | 2007-07-19 | 2011-11-01 | Toyota Jidosha Kabushiki Kaisha | Abnormality detection device for internal combustion engine and air/fuel ratio control apparatus for internal combustion engine |
FR2929650A1 (en) * | 2008-04-04 | 2009-10-09 | Bosch Gmbh Robert | METHOD AND DEVICE FOR ADAPTING A DYNAMIC MODEL OF AN EXHAUST GAS PROBE. |
EP2687709A4 (en) * | 2011-03-16 | 2017-01-11 | Toyota Jidosha Kabushiki Kaisha | Internal combustion engine control apparatus |
Also Published As
Publication number | Publication date |
---|---|
DE602006000790D1 (en) | 2008-05-08 |
JP2006336644A (en) | 2006-12-14 |
FR2886345A1 (en) | 2006-12-01 |
EP1729001B1 (en) | 2008-03-26 |
US20060271271A1 (en) | 2006-11-30 |
JP4964503B2 (en) | 2012-07-04 |
FR2886345B1 (en) | 2010-08-27 |
US7483782B2 (en) | 2009-01-27 |
DE602006000790T2 (en) | 2008-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1729001B1 (en) | Estimation method using a non-linear adaptive filter of air-fuel mixture richness in a cylinder of a combustion engine | |
EP1571318B1 (en) | Method for estimating richness of air-fuel mixture in a cylinder of a combustion engine | |
EP1729000B1 (en) | Method for estimating the air/fuel ratio in a cylinder of an internal combustion engine using an extended Kalman filter | |
FR3065990A1 (en) | METHOD FOR REALIZING A DYNAMIC OF ADAPTING A WEALTH VALUE TO A SET IN A MOTOR | |
EP2148979B1 (en) | Method of controlling the combustion of a diesel engine | |
FR3074524A1 (en) | SYSTEM AND METHOD FOR CONTROLLING AN INTERNAL COMBUSTION ENGINE EQUIPPED WITH A SELECTIVE CATALYSIS-TYPE EXHAUST GAS POST-PROCESSING SYSTEM | |
EP1931868B1 (en) | Method for estimating instantaneous speed produced by each of the cylinders of an internal combustion engine | |
WO2010112719A1 (en) | Control system and method for estimating the flow rate of recycled exhaust gases in an internal combustion engine | |
WO2015059034A1 (en) | System and method for estimating the flow of nitrogen oxides in the exhaust gases of an internal combustion engine for a motor vehicle | |
EP2195519B1 (en) | Engine state parameter estimation comprising the measurement of the internal pressure of a cylinder | |
EP2430298B1 (en) | Estimation of the nitrogen oxide (nox) concentration in an internal combustion engine | |
EP1890024A1 (en) | Determination of combustion start in an internal combustion engine | |
EP2829711A1 (en) | Method and device for determining the concentration of exhaust gas recirculating at the entrance to the intake distributor of a heat engine | |
FR2835281A1 (en) | Method for estimating mass of air admitted into engine combustion chamber consists of modeling air mass as function of pressure variation in combustion chamber from polytropic gas compression law | |
EP2507491B1 (en) | System and method for estimating the mass of particles stored in a particle filter of a motor vehicle | |
EP1662121A1 (en) | Method of controlling an intake system of an internal combustion engine and motor vehicle for carrying out said method | |
WO2014095052A1 (en) | Method for determining the recycled air flow rate and the quantity of oxygen available at the inlet of an internal combustion engine cylinder | |
WO2008043952A2 (en) | System for determining the mass flow rate of nitrogen oxides emitted in the exhaust gases of an internal combustion engine | |
EP1647692A1 (en) | Air inlet control method for an internal combustion engine and automotive vehicle for applying this method | |
FR2833649A1 (en) | Procedure for estimating and predicting mass of air entering i.c. engine combustion chamber includes measurement of inlet manifold pressure and calculation of Kalman gain | |
EP2811140A1 (en) | System and method for determining the mass fraction of fresh gas in the intake manifold of an internal combustion engine of a motor vehicle | |
FR2936015A1 (en) | State variables estimation system for diesel engine of vehicle, has multilayer type neuron network for estimating variables of state of engine from tables of intermediate temporal variables | |
FR3053117A1 (en) | METHOD FOR DETECTING AN AIR LEAK DOWNWARD OF A COMPRESSOR OF A TURBOCHARGER | |
FR2999648A1 (en) | Method for determining e.g. nitrogen dioxide concentration on outlet side of diesel engine, involves determining nitrogen oxide concentration from maximum temperature depending on engine angular position when maximum temperature is reached | |
FR2909762A1 (en) | Cylinder combustion noise determining system for internal combustion engine of motor vehicle, has determining module determining combustion noise in cylinder from selected pressure values in chamber of cylinder in real time |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LI LT LU LV MC NL PL PT RO SE SI SK TR |
|
AX | Request for extension of the european patent |
Extension state: AL BA HR MK YU |
|
17P | Request for examination filed |
Effective date: 20070606 |
|
AKX | Designation fees paid |
Designated state(s): DE GB IT |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): DE GB IT |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D Free format text: NOT ENGLISH |
|
REF | Corresponds to: |
Ref document number: 602006000790 Country of ref document: DE Date of ref document: 20080508 Kind code of ref document: P |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed |
Effective date: 20081230 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R081 Ref document number: 602006000790 Country of ref document: DE Owner name: IFP ENERGIES NOUVELLES, FR Free format text: FORMER OWNER: INSTITUT FRANCAIS DU PETROLE, RUEIL-MALMAISON, FR Effective date: 20110331 Ref country code: DE Ref legal event code: R081 Ref document number: 602006000790 Country of ref document: DE Owner name: IFP ENERGIES NOUVELLES, FR Free format text: FORMER OWNER: INSTITUT FRANCAIS DU PETROLE, RUEIL-MALMAISON, HAUTS-DE-SEINE, FR Effective date: 20110331 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 20110421 Year of fee payment: 6 |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20120403 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20120403 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 20160426 Year of fee payment: 11 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: IT Payment date: 20160429 Year of fee payment: 11 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R119 Ref document number: 602006000790 Country of ref document: DE |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20171103 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IT Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20170403 |