La présente invention concerne le domaine de la régulation d'un moteur à combustion interne, en particulier du type HCCI (« Homogeneous Charge Compression Ignition » en anglais, pour « allumage par compression d'une charge homogène ») utilisant un capteur de pression cylindre.

La spécificité de certains types de moteurs comme les moteurs HCCI par exemple où la combustion est spontanée et se fait par auto-allumage lors de la compression homogène du mélange carburant/comburant, ou les moteurs à très forte stratification rendent nécessaire un contrôle accru de la combustion.

La détermination de l'instant de combustion n'étant pas faisable en temps réel, le contrôle de la combustion d'un moteur se fait par exemple par la détermination de l'angle du vilebrequin (ou l'instant) CA_X correspondant à un pourcentage MFB_X donné (Masse Fraction Brûlée égale à X% de la masse totale brûlée dans le cylindre lors d'un cycle complet). Selon l'invention CA_X correspond donc indifféremment à l'angle vilebrequin ou à l'instant pour lequel X% du mélange est brûlé.

Le paramètre CA_X permet par exemple de piloter ensuite la combustion en jouant sur des paramètres tels que l'instant et la durée d'injection, le réglage des soupapes et autres éléments de contrôle, etc.

La fraction de la masse des gaz brûlés MFB_X n'étant pas mesurable, elle est déterminée à l'aide de modèles thermodynamiques basés sur l'information de la pression cylindre, par exemple celui de Rassweiler-Withrow.

Classiquement, la pression cylindre est acquise par des capteurs de pressions actifs, c'est à dire à électronique intégrée, du type piézo-électrique, montés sur chaque bougie de chaque cylindre ou dans la culasse.

Du fait d'une électronique intégrée, par exemple de type ASIC (« Application Spécifie Integrated Circuit » en anglais, pour « circuit intégré pour application spécifique »), les inconvénients de ces capteurs sont notamment un coût élevé, et un procédé de fabrication compliqué, car l'électronique intégrée dans chaque capteur doit compenser le gain et l'offset du capteur mesuré lors des tests de production.

La présente invention a pour but de remédier à ces inconvénients en proposant une solution moins onéreuse, mettant en œuvre des capteurs de pression passifs, c'est-à- dire ne comportant pas d'électronique intégrée.

Avec cet objectif en vue, l'invention concerne plus précisément, selon un premier de ses objets, un procédé de détermination de la fraction de la masse des gaz brûlés (MFB) lors d'un cycle de moteur à combustion interne de véhicule automobile. Selon l'invention, le procédé est essentiellement caractérisé en ce qu'il comprend des étapes consistant à :

- effectuer une première (Y(k)) et une deuxième (Y(k+1)) mesure de la pression dans un cylindre, respectivement en un premier instant (k) et un deuxième instant (k+1), par un capteur de pression passif possédant un offset (Y _Off ),

- déterminer la valeur V(k) et V(k+1) du volume dudit cylindre auxdits premier (k) et deuxième (k+1) instants respectifs, et

- calculer la fraction de la masse des gaz brûlés (MFB) à un instant (j) selon l'équation :

Où Y est le coefficient polytropique d'expansion des gaz, et m l'instant de fin du cycle moteur complet.

Grâce à l'invention, on peut déterminer le paramètre CA_X sans avoir à déterminer la valeur du gain K du capteur mais seulement son offset, celui-ci étant avantageusement déterminé entre la phase de fermeture des soupapes et la phase d'allumage.

Dans un mode de réalisation, le procédé comprend en outre une étape consistant à calculer l'offset du capteur de pression passif.

De préférence, pour le calcul de l'offset Y _Off du capteur de pression passif, la première Y(k) et la deuxième Y(k+1) mesure de la pression dans le cylindre sont effectuées au cours d'une période de temps correspondant à une compression hors combustion.

Dans un mode de réalisation, y, le coefficient polytropique d'expansion des gaz, est assimilé à une constante dont la valeur est comprise dans l'ensemble [1 ,3 ; 1 ,4].

Dans un autre mode de réalisation, le procédé comprend en outre une étape consistant à calculer y, le coefficient polytropique d'expansion des gaz, selon l'équation :

γ = (1 - MFB).γ _Q + MFB.γ _x

Où MFB est la fraction de la masse des gaz brûlés lors du ou d'un cycle précédent, et χ ₀ et γ _χ des valeurs prédéterminées.

Avantageusement, le procédé comprend en outre une étape de filtrage de la mesure Y(k) de la pression dans le cylindre.

De préférence, l'étape de filtrage est mise en œuvre par un filtre de Kalman. Dans un mode de réalisation, le procédé comprend en outre une étape consistant à calculer la valeur de l'instant) pour lequel la fraction de la masse des gaz brûlés MFB est égale à une valeur prédéterminée, en particulier 50%.

Dans un mode de réalisation, le procédé comprend en outre une étape consistant à asservir, pour le ou un cycle moteur donné, tout ou partie des actuateurs du moteur selon une loi dont une variable d'entrée est la valeur de l'instant y pour lequel la fraction de la masse des gaz brûlés MFB est égale à une valeur prédéterminée, calculée lors du ou d'un cycle précédent.

L'invention concerne également un programme d'ordinateur, comprenant des instructions de code de programme pour l'exécution des étapes d'un procédé tel que défini ci-avant lorsque ledit programme est exécuté sur un ordinateur.

D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description suivante donnée à titre d'exemple illustratif et non limitatif et faite en référence aux figures annexées dans lesquelles :

- la figure 1 illustre une vision synoptique d'un mode de réalisation du procédé selon l'invention, et

la figure 2 illustre l'évolution de la pression cylindre P et de la fraction de masse brûlée MFB en fonction de la position angulaire du vilebrequin CA. Quelques événements notables sont repris sur cette figure en abscisses. Le point I, tout d'abord, correspond à la fermeture des soupapes et au début de la phase de compression, repérée II, du mélange carburant / comburant. Le point IiI correspond pour sa part à l'allumage dudit mélange. Enfin, le point CA50 correspond au point de l'abscisse auquel 50% de la masse totale à brûler a bien été consumée.

L'un des buts de la présente invention est de permettre le contrôle d'un moteur de véhicule automobile à combustion interne par la connaissance dans un cycle donné de l'instant ou de l'angle vilebrequin pour lequel la fraction de la masse des gaz brûlés MFB est égale à une valeur prédéterminée.

La fraction de masse brûlée MFB à un instant y du cycle est définie par l'équation :

MFB = mbG) / mb(total)

Où mb(j) représente la fraction de la masse des gaz brûlés à l'instant y du cycle, et mb(total) la masse totale des gaz brûlés durant le cycle moteur complet. Selon l'invention, la fraction de masse brûlée MFB à un instant ; du cycle est approximée par l'équation

MFB = ^ (1 )

∑ΔP _e (*)

Jt=O

Dans l'équation (1) :

- k représente un instant ou une unité d'angle vilebrequin,

m correspond à l'instant ou à l'angle vilebrequin à la fin d'un cycle moteur complet, typiquement 720° pour un moteur à combustion interne de type « 4 temps », et

ΔP _C représente la variation de la pression dans le cylindre dans l'intervalle de temps considéré, due uniquement à la combustion.

Or, la pression P dans le cylindre à l'instant j+ 1 est la résultante de la pression de combustion Pc à l'instant j+1 et de la variation de pression de compression Pv, due au mouvement du piston sur l'intervalle entre l'instant; et l'instant j+1, soit :

p(y + i) = p _κ (y + i) + ΔP _c (y + i) d ¹ )

La composante de la pression Pv venant de la compression à l'instant (j+1) est donnée par l'équation polytropique d'expansion des gaz, soit :

Dans laquelle :

• P est la pression réelle,

• Pv est la pression de compression,

• V(j) et V(j+1) sont les volumes du cylindre aux instants; et j+1, calculés selon la position du piston à ces deux instants respectifs, et

• γ est le coefficient polytropique d'expansion des gaz.

On a donc, d'après les équations (1 ') et (1") :

AP _Γ C UKJ + i) J = P ( KjJ + \) J - P( KjJ) J x (2)

Or la pression mesurée Y par un capteur est différente de la pression réelle P. En effet, pour tout capteur, la pression mesurée Y à un instant j+1 est égale à la pression réelle P à l'instant j+1 multipliée par un gain K et à laquelle un offset Y _off est ajouté, soit :

Y(j + \) = K ^χ PU + \) + Y _off (3)

Les capteurs actifs permettent de s'affranchir des contraintes de gain et d'offset. Toutefois, l'invention met en œuvre de préférence des capteurs passifs. Dans ce contexte, l'invention permet avantageusement de s'affranchir de la connaissance du gain K et d'estimer la valeur de l'offset Y _Off , comme décrit ci-après.

A cet effet, à partir de l'équation (3), on obtient pour l'instant j+1 et pour l'instant y respectivement les équations :

Y(j + \) - Y _off

^/1 O ^' + !) = - (4)

K

et

PU) = (4 ¹ )

K

Des équations (2) et (4), on obtient :

Et des équations (5) et (4'), on obtient :

Soit

(6)

Des équations (1) et (6), on obtient donc l'équation :

qui permet de calculer, au moyen d'un calculateur, la fraction de la masse des gaz brûlés. La fraction exprimée dans l'équation (7) permet avantageusement de s'affranchir du gain K du capteur.

Dans cette équation, les grandeurs Y(k) sont les mesures effectuées par le capteur de préférence dans un même cycle moteur. Les grandeurs V(k) sont connues (données géométriques du moteur calculables ou stockées dans une mémoire), elles sont par exemple déterminées par un calculateur.

Dans un mode de réalisation, l'offset Y _Off est connu et stocké dans une mémoire.

Il peut par exemple être assimilé à une constante.

Dans un autre mode de réalisation, l'offset Y _off du capteur de pression passif est calculé, et éventuellement stocké dans une mémoire. Le calcul (l'estimation) de l'offset Y _0K est avantageusement mis en œuvre lorsque la première Y(k) et la deuxième Y(k+1) mesure de la pression dans le cylindre sont effectuées au cours d'une période de temps correspondant à une compression pure, c'est-à-dire hors combustion (voir figure 2).

En effet, en posant C _{} J+i} le ratio entre deux valeurs (de préférence successives) de la pression, on a :

_ P(j + \)

J ^{J+ι ~} PU)

Ainsi, d'après les équations (4) et (4'), on a :

Soit

Or ₁ lorsque le calcul, c'est-à-dire la première Y(k) et la deuxième Y(k+1 ) mesure de la pression dans le cylindre, sont effectuées au cours d'une période de temps correspondant à une compression hors combustion (entre la fermeture des soupapes d'admission et l'allumage), on a sur cette période particulière ΔP _C = 0. Grâce à cette caractéristique, on déduit de l'équation (2) que

c - (JX-LT m

Comme vu précédemment, pour chaque instant j les valeurs de V(j) sont connues. On peut donc calculer le ratio C _{J J+ι} , et la valeur de l'offset, grâce aux équations (8) et (8').

De préférence, la valeur calculée de l'offset Y _Off pour un cycle donné est stockée dans une mémoire et utilisée pour le pilotage / contrôle du moteur pour le cycle suivant.

On peut prévoir également en outre une étape de filtrage (du signal) de la mesure Y(k) de la pression dans le cylindre, par exemple par un filtre de Kalman.

Dans ce cas, on peut réécrire l'équation (8) sous la forme :

Y(J + 1) = C _{7 7+1} .Y(J) + (1 - C _wtl ) Y _off (9)

Pour l'ensemble des intervalles de la période de temps considérée, l'équation (9) peut s'exprimer sous forme matricielle, qui peut être résolue par la méthode des moindres carrés. Mais pour ne pas avoir à stocker n mesures et inverser a posteriori une matrice de taille importante, il est possible d'utiliser un filtre prédicteur-correcteur pour une estimation en continue basée sur les dernières mesures. Dans ce but, en utilisant un filtre de Kalman, on résout l'équation suivante :

YU + l) = C _JJ+ι .Y(j) + (l - C _JJ+ι )Y _off

Où Y et Y _off désignent respectivement la quantité prédite et la valeur de l'offset utilisée dans le modèle (valeur déterminée à l'instant j, ou à défaut égale au cycle précédent par exemple pour une initiation).

L'erreur ε observée au temps (j+1) est égale à :

ε ( j + l) = Y (j + l ) - Ϋ ( j + l)

Et l'offset utilisé dans l'étape (cycle) suivante est corrigé selon l'équation :

Ϋ _off U + l) = Y _off (j) + G ^χ εU + l)

Où G est le gain de la résolution de la méthode numérique, tel que

G = (l - C, _J+l )/g ,

avec g le gain du filtre de Kalman compris entre 0 et 1.

La mise en œuvre d'une telle étape de filtrage est avantageusement optimisée par le choix des plages de temps traités. Cette estimation, réalisable de manière récursive à l'aide d'un filtre de Kalman, permet de s'affranchir d'un besoin de grande capacité de stockage des mesures et d'inversion de matrice volumineuse.

Le calcul de la fraction de la masse des gaz brûlés (équation (7)) implique également la connaissance de la valeur de γ le coefficient polytropique d'expansion des gaz.

Dans un mode de réalisation, la valeur de γ est constante, en particulier au cours d'une période de temps correspondant à une compression hors combustion, et comprise dans l'ensemble [1 ,3 ; 1 ,4], en particulier pour des valeurs de fraction de masse brûlée inférieures ou égales à 50%.

Dans un autre mode de réalisation, en particulier pour des valeurs de fraction de masse brûlée supérieures à 50%, on calcule la valeur de γ à partir d'un (du) cycle précédent, selon l'équation :

γ = (1 - MFB).γ _Q + MFB.γ _λ

Où γ ₀ et /, sont des données connues / estimées par le contrôle moteur.

Grâce aux valeurs de l'offset Y _Off et de γ , on peut alors calculer la fraction de masse brûlée MFB selon l'équation (7).

La détermination de l'angle du vilebrequin (ou l'instant) CA_X correspondant à un pourcentage MFB_X donné peut alors être mis en œuvre par exemple par interpolation. L'invention est avantageusement mise en œuvre pour la détermination de l'angle du vilebrequin (ou l'instant) CA_50 correspondant à un pourcentage MFB_50 de 50% de la masse brûlée.

Cette valeur pour un cycle donné permet ensuite de piloter divers éléments du moteur au cycle suivant, par exemple pour asservir, pour le cycle moteur, tout ou partie des actuateurs du moteur selon une loi dont une variable d'entrée est la valeur de l'instant y pour lequel la fraction de la masse des gaz brûlés (MFB) est égale à la valeur MFB_X prédéterminée, en l'espèce MFB_50, calculée lors du ou d'un cycle précédent.