FR 2947864 A1 20110114 FR 0903369 A 20090708 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 20 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 25 Specific 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 30 une solution moins onéreuse, mettant en oeuvre 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 (Yoff), - 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 (V(kû1)y x(Yn ûY(k-1)) Y(k) û Yof. + V(k) k=0 MFB= 10 Y(k) û Yoff + (V(kû1) V(k) x(Yoff -Y(kû1)) k=0 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 15 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 Yoff du capteur de pression passif, la 20 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]. 25 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 : y = (1û MFB)•yo + MFB.y Où MFB est la fraction de la masse des gaz brûlés lors du ou d'un cycle précédent, et ro et y des valeurs prédéterminées. 30 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 oeuvre par un filtre de Kalman. 3 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] 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 10 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 15 - 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 brulée MFB en fonction de la position angulaire du vilebrequin CA. Quelques événements notables sont repris sur cette figure en abscisses. Le 20 point I, tout d'abord, correspond à la fermeture des soupapes et au début de la phase de compression, repérée Il, du mélange carburant/ comburant. Le point III correspond pour sa part à l'allumage dudit mélangea Enfin, le point CA50 correspond au point de l'abscisse auquel 50% de la masse totale à brûler a bien été consumée. 25 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 brulée MFB à un instant j du cycle est définie par 30 l'équation MFB = mb(j) / mb(total) Où mb(j) représente la fraction de la masse des gaz brûlés à l'instant j 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 brulée MFB à un instant j du cycle est approximée par l'équation l 1AI'c(k) MFB = k-° (1) 1 AP (k) k=0 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 Al" 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(j+1)=Pv(j+1)+APc(j+1) (1') 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 : Pv(j+1)=P(j)x V( j) ~~ (1") V(j+1)) Dans laquelle • P est la pression réelle, • Pv est la pression de compression, • Va) et V(j+1) sont les volumes du cylindre aux instants j et j+1, calculés selon la position du piston à ces deux instants respectifs, et • y est le coefficient polytropique d'expansion des gaz. On a donc, d'après les équations (1') et (1") : APc(j+1)=P(j+1)ùP(j)x~ V(j) ~~ (2) y(�fi)) 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 Yoff est ajouté, soit : Y(j+1)=KxP(j+l)+Yoe. (3) Les capteurs actifs permettent de s'affranchir des contraintes de gain et d'offset. Toutefois, l'invention met en oeuvre 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 Yoff, 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 j respectivement les équations (~ YU + 1) û Ya~ (4) P '+1) = K YU) û Yo P(j) _ K Des équations (2) et (4), on obtient : 10 APc(j+1) Y(j+1)ùYotf _P(j)X V(j) K y(j+1)) Et des équations (5) et (4'), on obtient : Y(j+1)ûYoff Y(j)ûYoff r V(j) m' APc(J+1)= K - K x Soit i V(j+1)l (4') (5) (5') APc(j+1)=K Y(j+1)ûYaf +(Yo ûY(j))x V(j) V(j+1)) i (6) 15 Des équations (1) et (6), on obtient donc l'équation : (V(kû1) x(YoJf ûY(kû1)) V(k) x(Yçff-Y(kû1)) MFB = Y(k) û Yo + YV(k 1)m' V(k) (7) Y(k) - Y) + 1 k=0 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. 20 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 Yoff est connu et stocké dans une mémoire. 25 Il peut par exemple être assimilé à une constante. Dans un autre mode de réalisation, l'offset Yoff du capteur de pression passif est calculé, et éventuellement stocké dans une mémoire. 6 Le calcul (l'estimation) de l'offset Yoff est avantageusement mis en oeuvre 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 effeten posant C11+1 le ratio entre deux valeurs (de préférence successives) de la pression, on a (Y(j)ù1'of) 10 Soit Y' (Y(j+ 1) ù Ch;+IY(j)) (8) (1ù C1,1+1 ) 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 15 d'admission et l'allumage), on a sur cette période particulière APc = O. Grâce à cette caractéristique, on déduit de l'équation (2) que V(j) " (8') Comme vu précédemment, pour chaque instant j les valeurs de VO) sont connues. On peut donc calculer le ratio C11+1 et la valeur de l'offset, grâce aux 20 équations (8) et (8'). De préférence, la valeur calculée de l'offset Yoff 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. 25 Dans ce cas, on peut réécrire l'équation (8) sous la forme : Y(j + 1) = C l,l+1 •Y(j) + (1ù C;,1+1)Yoff (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 30 taille importante, il est possible d'utiliser un filtre prédicteur-correcteur pour une estimation en continue basée sur les dernières mesures. P(j+1) Ainsi, d'après les équations (4) et (4'), on a : (Y(j+1)-Ynrr) C1,1+1 = Duns ce but, en utilisant un filtre de Kalman, on résout l'équation suivante j + 1) .Y( j) +�l ù C~,~+~ )Yo Où @ et @, 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 e observée au temps (j+1) est égale à e(j+1) = Y(j+1)ù@(j+1) Et l'offset utilisé dans r étape (cycle) suivante est corrigé selon l'équation : Y0.f(j+1)=Y0ff(j)+Gxe(j+1) 10 Où G est le gain de la résolution de la méthode numérique, tel que G = (1ùCi,i+1)/g avec g le gain du filtre de Kalman compris entre 0 et 1. La mise en oeuvre 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 15 à 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 y le coefficient polytropique d'expansion des gaz 20 Dans un mode de réalisation, la valeur de y 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 25 masse brûlée supérieures à 50%, on calcule la valeur de y à partir d'un (du) cycle précédent, selon l'équation y = (1- MFB).yo + MFB.y1 Où ro et yl sont des données connues/ estimées par le contrôle moteur. Grâce aux valeurs de l'offset Yoff et de y, on peut alors calculer la fraction de 30 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 oeuvre par exemple par interpolation. L'invention est avantageusement mise en ceuvre 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] 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. REVENDICATIONS1. 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, 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 (Yoff), - déterminer la valeur V(k) et V(k+1) du volume dudit cylindre auxdits premier (k) et deuxième (k+1) instants respectifs, - calculer la fraction de la masse des gaz brûlés (MFB) à un instant j selon l'équation Y(k) k ûY +~V(kû1)~~ x -Y(kû1 or~V(k)i )~ ~ Y(k) - off f + V(k û (k) 1) x (Yo - Y(k -1) V Où y est le coefficient polytropique d'expansion des gaz, et m l'instant de fin du cycle moteur complet. 2. Procédé selon la revendication 1, comprenant en outre une étape consistant à 15 calculer l'offset (Y0ff) du capteur de pression passif. 3. Procédé selon la revendication 2, dans lequel pour le calcul de l'offset (Yoff) 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. 20 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel y, le coefficient polytropique d'expansion des gaz, est assimilé à une constante dont la valeur est comprise dans l'ensemble [1,3 ; 1,4]. 5. Procédé selon l'une quelconque des revendications 1 à 3, comprenant en outre une étape consistant à calculer y, le coefficient polytropique d'expansion des gaz, 25 selon l'équation y= (1ûMFB).yo +MFB. y, Où MFB est la fraction de la masse des gaz brûlés lors du ou d'un cycle précédent, et yo et y, des valeurs prédéterminées. 6. Procédé selon l'une quelconque des revendications précédentes, comprenant 30 en outre une étape de filtrage de la mesure (Y(k)) de la pression dans le cylindre. 7. Procédé selon la revendication 6, dans lequel l'étape de filtrage est mise en oeuvre par un filtre de Kalman. MFB k=0 ) 8. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre une étape consistant à calculer la valeur de l'instant (j) pour lequel la fraction de la masse des gaz brûlés (MFB) est égale à une valeur prédéterminée, en particulier 50%. 9. Procédé selon la revendication 8, comprenant en outre une étape consistant à 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 (j) 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. 10. Programme d'ordinateur, comprenant des instructions de code de 10 programme pour l'exécution des étapes du procédé selon l'une quelconque des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur.