FR 2949511 A1 20110304 FR 0904084 A 20090826 10 La présente invention concerne le domaine du contrôle de la combustion d'un moteur à combustion interne. L'invention concerne notamment une méthode pour estimer l'intensité de cliquetis d'un tel moteur au moyen d'un capteur placé dans le moteur. Ce type de moteur comprend au moins un cylindre comportant une chambre de 15 combustion délimitée par la paroi latérale interne du cylindre, par le haut du piston qui coulisse dans ce cylindre et par la culasse. Généralement, un mélange carburé est renfermé dans cette chambre de combustion et subit une étape de compression puis une étape de combustion sous l'effet d'un allumage commandé, par une bougie, ces étapes étant regroupées sous le vocable de � phase de combustion � dans la suite de 20 la description. Dans un moteur à explosion, la combustion du mélange air/essence commence normalement après l'étincelle. Le front de flamme se propage et son souffle repousse une partie du mélange contre les parois du cylindre et le sommet du piston. L'élévation de pression et de température devient tellement importante que le combustible coincé 25 contre les parois atteint son point d'auto-allumage et s'auto enflamme à plusieurs endroits. Ce phénomène est appelé � cliquetis �. Ainsi le cliquetis est avant tout un phénomène de combustion anormale dans les moteurs à allumage commandé, perceptible extérieurement par un bruit métallique venant du moteur. Il résulte de l'apparition d'ondes de choc dans la chambre de combustion. 30 Les micro explosions qui en résultent produisent des vibrations dans le domaine acoustique (de l'ordre de 5 à 20 KHz). Elles sont très vives et peuvent rapidement créer des points chauds qui accentueront encore plus le problème. L'accumulation de micro explosions va arracher ou faire fondre une petite quantité de métal sur le 1 sommet du piston et/ou sur les parois du cylindre et des segments. Au bout de quelques temps (selon l'intensité) cela conduira à la destruction du piston, des segments ou des parois du cylindre. L'estimation de l'intensité de cliquetis permet un contrôle de la combustion limitant l'effet du cliquetis et n'endommageant pas le cylindre. État de la technique On connaît des méthodes pour estimer le cliquetis d'un moteur à combustion interne. Elles sont basées sur l'enregistrement de signaux issus d'un capteur de pression cylindre. Une première méthode consiste à appliquer hors ligne, c'est à dire en traitement ultérieur, une technique de transformations de Fourier à ces signaux pour estimer le cliquetis. Une telle méthode est décrite dans les documents suivants : "Application of the Reallocated Smoothed Pseudo Wigner-Ville Distribution to Knock Detection'; Olivier BOUBAL et Jacques OKSMAN, Traitement du Signal volume 15, 1998 "Knock measurement for fuel evaluation in spark ignition engines", C. Hudson, X Gao, R. Stone, Fuel volume 80 (2001) Une seconde méthode consiste appliquer aux signaux une technique de filtrage et de détections de maximum pour estimer le cliquetis hors ligne. Une telle méthode est décrite dans le document suivant : "Modeling of ln-cylinder Pressure Oscillations under Knocking Conditions: Introduction to Pressure Envelope Curve", G. Brecq et O. Le Corre, SAE 2005 Enfin, on connaît une autre méthode, également basée sur l'enregistrement de signaux issus d'un capteur de pression cylindre. Cette fois, une technique d'ondelettes est appliquée aux signaux pour estimer le cliquetis hors ligne. Une telle méthode est décrite dans le document suivant : "Knocking detection using wavelet instantaneous correlation method", Z. Zhang et , E. Tomita, Journal of SAE Review volume 23 (2002) Cependant, c'est méthodes ne permettent pas une estimation de l'intensité du cliquetis en temps réel. L'objet de l'invention concerne une méthode pour déterminer l'intensité du cliquetis dans un moteur à combustion interne. Cette estimation est réalisée en temps réel de façon à contrôler au mieux la combustion et préserver le moteur. La méthode est basée sur un traitement d'un signal vibratoire représentatif des vibrations du moteur, d'un filtrage en temps réel de ce signal vibratoire à l'aide d'un filtre, et dune détermination en temps réel des coefficients de la décomposition de Fourier du signal. La méthode selon l'invention L'invention concerne une méthode pour estimer l'intensité de cliquetis d'un moteur à combustion interne comprenant au moins un capteur vibratoire, à partir duquel on acquiert en continu et en fonction de l'angle vilebrequin un signal vibratoire représentatif des vibrations du moteur, caractérisée en ce que l'on réalise les étapes suivantes : on réalise en temps réel un filtrage dudit signal vibratoire à l'aide d'un filtre ; - on construit un modèle de la dynamique dudit filtre que l'on couple avec un estimateur ; on détermine en temps réel des coefficients de la décomposition de Fourier dudit signal vibratoire acquis, en inversant la dynamique dudit filtre à l'aide dudit estimateur ; - on déduit en temps réel un paramètre corrélé à l'intensité de cliquetis à partir desdits coefficients. Selon l'invention, le signal vibratoire peut être un signal issu d'un capteur de pression cylindre ou un signal issu d'un capteur accéléromètrique. Le paramètre peut être fonction de l'énergie contenue dans le signal. Cette énergie est calculée en sommant les carrés des coefficients de la décomposition de Fourier. Ce paramètre est par exemple égal à la racine carrée du maximum de cette énergie. Selon un mode de réalisation, le signal vibratoire comporte plusieurs harmoniques et l'on estime lesdits coefficients de la décomposition de Fourier uniquement pour l'harmonique ou les harmoniques corrélées avec la combustion. On peut sélectionner cette harmonique ou ces harmoniques corrélées avec la combustion en réalisant une analyse spectrale ou une analyse temps/fréquence, ou en choisissant la bande de fréquence 5000Hz - 20000Hz. Selon un autre mode de réalisation, on détermine une plage angulaire incluant l'intégralité de la réponse vibratoire liée à la combustion, et l'on restreint le signal 5 vibratoire à cette plage angulaire avant de réaliser ledit filtrage. Selon l'invention, le filtre peut être un filtre passe-bande, et l'estimateur peut être un estimateur non linéaire de type adaptatif. Enfin, selon l'invention, on peut modifier des paramètres moteurs en fonction de du paramètre d'intensité de cliquetis, pour contrôler la combustion du moteur. D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. 15 Présentation succincte des figures la figure 1 montre le schéma du principe du contrôle de la combustion selon l'invention. la figure 2 montre un schéma fonctionnel décomposant l'estimation de l'intensité du cliquetis via le paramètre Ic,. Description détaillée de la méthode La méthode selon l'invention concerne une méthode de contrôle de la combustion d'un moteur à combustion interne prenant en compte les phénomènes de cliquetis. La méthode comporte l'estimation de l'intensité de cliquetis du moteur, à 25 partir d'un traitement d'un signal issu d'un capteur, ce signal ayant un sens physique pour qualifier directement la combustion et l'éventuelle apparition de cliquetis. La méthode est indépendante de la nature de la combustion et de la technologie du capteur de pression cylindre utilisé. L'invention est décrite selon un mode de réalisation particulier, dans lequel on 30 décrit la méthode appliquée à un cylindre unique, sans intégrer de particularismes associés au mode de combustion du moteur. Selon cet exemple de réalisation, on utilise un capteur de pression cylindre. La nature du moteur à combustion interne ainsi 10 20 que la technologie du capteur de pression cylindre utilisé n'influence pas le principe de la méthode. La figure 1 illustre le principe général de la méthode de contrôle selon l'invention. Dans un premier temps on acquiert un signal de pression cylindre (p(t)) issu d'un capteur de pression cylindre (A). Ce signal de pression correspond à la combustion dans la chambre de combustion (CC). On sait, que seuls certains harmoniques de ce signal contiennent de l'information qualifiant le cliquetis. L'idée est donc de calculer les différentes harmoniques de ce signal vibratoire en temps réel, puis de sélectionner celles qui correspondent à la combustion. Enfin, à partir de ces harmoniques portant une information sur la combustion, on calcule (CAL) l'intensité de cliquetis qualifiant directement la combustion, toujours en temps réel, et permettant l'ajustement (CTRL) des paramètres du moteur pour réguler la combustion. Comme l'illustre la figure 2, la méthode de contrôle moteur de la combustion peut se décomposer en cinq étapes principales : Étape préliminaire : Étape 1- Sélection des harmoniques qualifiant la combustion (SLCT) Étapes réalisées en temps réel : Étape 2- Acquisition d'un signal issu d'un capteur de pression cylindre (ACQ) Étape 3- Détermination en temps réel des harmoniques du signal (HAR) Étape 4- Détermination en temps réel de l'intensité du cliquetis pour qualifier directement la combustion (CAL) Étape 5- Contrôle moteur : Ajustement des paramètres du moteur pour réguler la combustion et éviter le cliquetis (CTRL) Étape 1 : Sélection des harmoniques qualifiant la combustion (SLCT) Toute grandeur oscillante périodique (ou grandeur alternativement croissante et décroissante), comme le signal de pression cylindre (ou le signal issu d'un capteur accéléromètrique), peut être considérée comme la somme de composantes sinusoïdales dont les fréquences sont multiples entiers de la fréquence fondamentale. Cette dernière, appelée aussi fréquence de base, est égale à l'inverse de la période de la grandeur considérée. Les fréquences multiples sont des composantes sinusoïdales dont la fréquence est un multiple entier de la fréquence fondamentale. Ces fréquences multiples sont appelées les harmoniques. Chacune des harmoniques portent une ou plusieurs informations : bruit, combustion, claquement des injecteurs, retournement du piston au point mort haut, vibrations induites par le système de distribution, etc. L'objectif est, dans un premier temps, de sélectionner, parmi ces composantes, celles qui contiennent de l'information qualifiant la combustion. Cette opération revient à extraire la partie utile du signal issu du capteur de pression cylindre tout en s'affranchissant des phénomènes parasites. Pour déterminer les harmoniques permettant de qualifier la combustion, deux méthodes peuvent être envisagées : l'analyse temps/fréquence ou l'analyse spectrale. Ces méthodes sont appliquées sur des signaux issus de capteurs de pressions cylindre dans le cadre d'expériences préalables à la mise en place des étapes permettant un contrôle de la combustion du moteur en temps réel. L'analyse temps/fréquence du signal peut être interprétée comme une fenêtre glissante dans laquelle on effectue une transformée de Fourier rapide (FFT), bien connue des spécialistes, à chaque pas de temps. Cette forme de transformée de Fourier est aussi connue comme la transformée de Fourier à temps court, et sa représentation dans le plan temps fréquence est donnée par un spectrogramme. Un spectrogramme est un diagramme associant à chaque instant t d'un signal, son spectre de fréquence. Dans son format le plus courant, l'axe horizontal représente le temps et l'axe vertical la fréquence. Chaque point à l'intérieur du graphique est doté d'une certaine intensité qui indique l'amplitude (souvent en décibels) d'une fréquence particulière à un temps donné. Dans le domaine de l'industrie automobile, il est courant de remplacer l'axe temporel par un axe angulaire, correspond à la valeur de l'angle vilebrequin. Sur un tel spectrogramme, la périodicité des combustions apparaît clairement sur l'axe angulaire. Ces événements correspondent à des phénomènes périodiques et synchrones avec les événements intervenant lors d'un cycle moteur dont une partie est liée à la combustion. Classiquement, dans le cadre d'analyse de signaux vibratoires issus du capteur de pression cylindre placé dans un moteur à combustion interne, on distingue clairement des zones d'étalement angulaire dans certaines bandes de fréquence entre 5kHz et 20kHz. En effet, entre un signal avec cliquetis et un signal sans cliquetis, on observe une différence très importante sur l'énergie dans les zones de fréquences de quelques kHz (le plus souvent entre 7 et 9, c'est pour cela que l'on étend le filtre entre 5 et 20). Cette fréquence dépend principalement de la géométrie du modèle. On obtient expérimentalement cette fréquence en réalisant quelques mesures sur le système de combustion. Ainsi, selon la méthode, on peut soit fixer la bande de fréquence d'intérêt pour qualifier la combustion entre 5kHz et 20kHz, soit déterminer, en fonction du type de moteur par exemple, la bande de fréquence adéquate à l'aide d'une analyse spectrale. Dans la suite de la description, on travaille dans la bande de fréquence 5kHz û 20kHz pour estimer l'intensité de cliquetis. Les harmoniques (composantes fréquentielles) du signal sélectionnées sont notées : wk avec k e [1, nh ] . L'entier nh est choisi de façon à ce que ces nh harmoniques appartiennent à la bande de fréquence 5kHz û 20kHz. Les étapes 2 à 5 sont réalisées en temps réel. Étape 2 : Acquisition d'un signal issu d'un capteur de vibration (ACQ) La méthode de contrôle de la combustion repose sur le traitement d'un signal vibratoire, ayant un sens physique pour qualifier directement la combustion et l'éventuelle apparition de cliquetis. Ce signal est délivré en continu par un capteur : un capteur de pression cylindre selon cet exemple de réalisation. Ce type d'acquisition est bien connu des spécialistes. La technologie du capteur de pression cylindre utilisé n'influence pas non plus le principe de la méthode. Le signal acquis par le capteur de pression cylindre (A) est noté p(t) . La variable t représente le temps. Le signal délivré par le capteur de pression cylindre, p(t) , constitue une des deux entrées de la chaîne de traitement. Le second signal, a(t) , est un signal fournissant la valeur de l'angle vilebrequin en fonction du temps. L'angle vilebrequin est noté a. Ce signal est délivré par un codeur angulaire (ENC). Étape 3 : Détermination en temps réel des harmoniques du siqnal L'objectif est, dans un premier temps, d'estimer en temps réel toutes les composantes harmoniques du signal vibratoire issu du capteur de pression cylindre, ou tout au moins, les composantes harmoniques utiles sélectionnées à l'étape 1. Le fait de reconstruire ces composantes en temps réel signifie que l'on n'enregistre pas un signal en continue pour le traiter ultérieurement pour obtenir les composantes harmoniques, par une décomposition de Fourier par exemple. Au contraire, ces composantes sont estimées au fur et à mesure que le signal est mesuré. Pour ce faire on utilise un observateur de Fourier. La méthode peut comporter les étapes suivantes : - Extraction de la partie utile du signal : prétraitement du signal vibratoire (PRE) - Estimation des harmoniques par un observateur de Fourier (HAR) Extraction de la partie utile du siqnal : prétraitement du siqnal vibratoire (PRE) Le traitement du signal issu du capteur de pression cylindre constitue un élément primordial dans la chaîne de contrôle moteur puisqu'il délivre des variables exploitables pour le contrôle bouclé de la combustion. La première étape consiste en un traitement du signal vibratoire pour tirer un maximum d'information utile. Ce prétraitement (avant la reconstruction des harmoniques et des paramètres de combustion) pour extraire la partie utile du signal vibratoire peut comporter les étapes suivantes : échantillonnage de p(t) dans le domaine angulaire, détermination d'une fenêtre angulaire incluant la combustion. Échantillonnage de p(t) dans le domaine angulaire Le signal p(t) est échantillonné dans le domaine angulaire. Ce signal est alors noté p(a) . Le domaine angulaire est, par opposition au domaine temporel, le domaine dans lequel toutes les données sont exprimées en fonction de l'angle vilebrequin. Cette étape permet de travailler dans un repère invariant comme pour toutes les machines tournantes. Le second intérêt à travailler dans le domaine angulaire, vient du fait que l'avance à l'allumage est pilotée de manière angulaire, par rapport au point mort haut (PMH), ce qui facilite le contrôle de l'avance pour la combustion. Mais il faut évidemment s'assurer du bon calage de l'encodeur angulaire qui définit l'origine du repère dans le domaine angulaire comme étant le point mort haut. Détermination d'une fenêtre angulaire [al ; ad incluant la combustion Le signal p(t) est mesuré en continu par une carte d'acquisition, puis transformé dans le domaine angulaire (p(a) ). Cependant, la suite du traitement est avantageusement effectuée uniquement sur une plage angulaire bornée par des angles al et a2. Ces deux paramètres font partie des éléments de calibration de l'observateur de cliquetis (observateur de Fourier selon l'invention). Les limites al et a2 dépendent de la plage angulaire d'occurrence de la combustion. Celle-ci peut être précoce (plusieurs degrés avant le PMH) ou, au contraire, extrêmement tardive (20 degrés après le PMH) selon le mode de combustion. Pour un moteur donné, les limites a, et a2 doivent être choisies afin d'inclure l'intégralité de la réponse vibratoire de la combustion, et ceci sur la totalité de la plage de fonctionnement du moteur. Ces angles sont définis par rapport au point mort haut qui représente l'origine du repère. Donc, l'angle a, est de signe négatif et a2 de signe positif. En somme, la méthode propose de traiter le signal vibratoire sur une fenêtre de temps déterminée du signal issu du capteur, dont la position correspond à une position angulaire déterminée du vilebrequin du moteur. Le signal du capteur contenant, à l'intérieur de cette fenêtre de temps ou de la fenêtre angulaire [a, ; a2] correspondante, des informations concernant le comportement vibratoire du moteur, permettant de déterminer des paramètres de combustion pour apprécier la qualité de la combustion. Estimation des harmoniques par un observateur de Fourier (HAR) L'idée est de filtrer le signal issu du capteur de pression cylindre p(a), puis d'inverser le filtrage par reconstruction en temps réel des composantes harmoniques utiles, en utilisant un observateur de Fourier. On peut reconstruire les composantes d'un seul ou d'une combinaison de plusieurs signaux de pression cylindre issus d'un ou plusieurs capteurs montés sur un moteur à combustion interne (bloc moteur et/ou culasse). Filtrage passe bande du signal p(a) sur la fenêtre angulaire fia, ; a2] (FILT) Après avoir été fenêtré, le signal p(a) est filtré à l'aide d'un filtre passe bande qui permet d'éliminer des composantes spectrales inutiles déterminées par analyse spectrale. Ce signal filtré sur la fenêtre [al ; a2] est noté p f (a). Les fréquences de coupure (F, et F2) du filtre constituent, elles aussi, des paramètres de réglage de l'observateur de combustion. Celles-ci sont notamment liées à la structure du moteur qui conditionne le contenu fréquentiel de la réponse vibratoire ainsi qu'au capteur qui est caractérisé par une bande passante dépendante de sa technologie. Une analyse spectrale du signal p(a) , telle que celle décrite à l'étape 2, permet de choisir la bande de fréquence utile (qualifiant la combustion) et de régler les fréquences de coupure du filtre passe bande (F, et F2). Reconstruction des harmoniques utiles par inversion du filtre passe bande (REC) On recherche à reconstruire (estimer), à partir de la mesure p f (a) du signal filtré issu du capteur de pression cylindre, les composantes fréquentielles (harmoniques) de ce signal. La méthode comporte principalement les trois étapes suivantes : définition d'un modèle de la dynamique du filtre passe-bande (ou du capteur vibratoire), caractérisation du signal vibratoire par des paramètres quasi invariants au cours du temps, couplage avec un estimateur non linéaire de type adaptatif, et éventuellement, vérification de la qualité de l'estimation. Modèle de la dynamique du filtre passe-bande Les vibrations engendrées par le moteur constituent une excitation, noté a), 5 enregistrée par un capteur vibratoire (le capteur de pression cylindre) sous la forme d'un signal, noté p(a). Cette mesure est ensuite filtrée p f (a) à l'aide d'un filtre passe-bande. Ce filtre peut se mettre, dans le domaine angulaire, sous la forme du système d'état suivant : dx = A.x+Ao.w da y = C.x (1) 10 avec : dpf(a) dnpf(a) da dan, appelé vecteur d'état a : l'angle vilebrequin w l'excitation du moteur y : la sortie du modèle, égale à C.x, c'est-à-dire le signal pf(a) mesuré 15 (A,A0,C) : les paramètres définissant le filtre passe-bande utilisé, avec : ùal 1 0 . 0- - a2 0 1. b, 0 ùan_, 0 . 0 1 ù an 0 0 A= , Ao = bn_1 bn , C = [1 0 . 0] de dimension lxn Les coefficients ao...an et bo...bn sont définis en fonction des fréquences FI et F2. Ils correspondent aux coefficients numériques du filtre passe-bande utilisé. 20 Cette équation (1) constitue le modèle représentant en temps réel la dynamique du filtre passe-bande. Caractérisation du siqnal w par des paramètres quasi invariants au cours du temps On recherche à estimer, à partir de ce modèle de la dynamique du filtre et de la mesure y (égale à p f (a) ), le signal w, c'est-à-dire la reconstruction du signal vibratoire engendré par le moteur p(a) avec les fréquences utiles (contenant l'information relative à la combustion). Au cours de la méthode selon l'invention, on caractérise ce signal w par des paramètres quasi invariants au cours du temps, ce qui permet de réaliser cette estimation en temps réel. En d'autres termes, on définit le signal w à l'aide de paramètres qui, à un instant donné, sont des constantes. Pour ce faire, on exploite le fait que le signal co est mécaniquement périodique. Ainsi, au lieu de réaliser une estimation du signal fortement variable w, on peut estimer les coefficients de la décomposition de Fourier de ce signal. On peut également utiliser tous paramètres permettant de décrire le signal w en relation avec son caractère périodique. La décomposition en coefficients de Fourier du signal co, développée en complexe pour la clarté de l'exposé, s'écrit ainsi : nh w(a) = Ecke(ia'ka) (2) k=1 Les Ck représentent les coefficients de Fourier de la décomposition du signal co. On définit ainsi un signal traduisant le signal vibratoire, en fonction des paramètres ck invariants au cours du temps. En n'utilisant qu'un nombre fini d'harmoniques, le modèle représentant en temps réel la dynamique du capteur vibratoire s'écrit alors : dx da dck da Y nh A.x+Ao.(Ecke(i'ka)) k=1 0 , Vk E [1, nh ] C.x (3) Selon l'invention, le but est d'estimer les coefficients (ck) de Fourier de la décomposition du signal w. Couplage avec un estimateur non linéaire de type adaptatif 12 On souhaite définir un estimateur, que l'on appelle � observateur de Fourier �, pour reconstruire l'excitation périodique du signal accéléromètrique à partir de sa mesure filtrée p f (a) . Pour cela, la dynamique du filtre est inversée et les composantes fréquentielles du signal accéléromètrique sont reconstruites. A partir du modèle décrit par le système (3), on définit un estimateur non linéaire de type adaptatif comportant d'une part, un terme lié à la dynamique et d'autre part, un terme de correction : dx nh ù=Az+A0~cke('wk") ûL(Cxûsf(a)) J da k=1 dck _ ûe( ,wka) Lk (Cz û s f (a)), Vk E [1, nh (4) da =Cz avec : estimée de x ck estimée de Ck L une matrice à calibrer Lk des matrices à calibrer pour assurer la convergence Un choix de matrices L et Lk assurant la convergence de l'estimateur (observateur) est : L= 1 Lk k2 +1 VkE[1,nn1 1 Le système d'équations (4) représente un estimateur non linéaire de type adaptatif permettant d'estimer les coefficients Ck de la décomposition en coefficients de Fourier du signal co. Cette méthode assure la reconstruction de l'excitation co à travers les coefficients Ck. Le signal accéléromètrique reconstruit w est donné par la relation suivante : �a) = > Ck (a)erwka Etl'ona: =Cz=pf(a) Choix des harmoniques On a considéré le cas où l'on utilisait les harmoniques choisies à l'étape 1 (wk ) pour reconstruire le signal : k varie entre 1 et nh, on utilise donc nh harmoniques, et ces harmoniques portent de l'information qualifiant la combustion. 10 Cependant, l'estimation w(a) de l'excitation peut être réalisée à partir d'un nombre d'harmoniques plus important. En effet, plus l'on utilise de coefficients de Fourier, plus la reconstruction de l'excitation w est précise. On ne sortirait donc pas du cadre de l'invention en choisissant un nombre d'harmoniques (nh) différent du nombre d'harmoniques sélectionnées à l'étape 1. En revanche, il est évident que plus le 15 nombre d'harmoniques choisi pour reconstruire l'excitation w est grand, plus le temps de calcul est important. II faut donc faire un compromis entre la rapidité de calcul et la précision. L'objectif est donc de déterminer à la fois les wk pertinentes mais également le paramètre nh adapté, c'est-à-dire, le choix du nombre d'harmoniques nécessaires à la 20 reconstruction du signal. Ce choix dépend des énergies respectives des différentes harmoniques. En général, les trois premières harmoniques sont les plus significatives pour qualifier la combustion et elles suffisent pour reconstruire l'excitation. On prendra dans la majeure partie des cas nh =1. Le signal ainsi reconstruit est noté 13(a) . 25 Exemple de calibration de l'observateur. Le tableau 1, où No, représente le nombre de cylindres du moteur, présente des exemples de valeurs des paramètres nécessaires à la calibration de l'observateur de Fourier.5 Paramètre Unité Description Valeur al degré borne inférieure de la fenêtre angulaire 0 a2 degré borne supérieure de la fenêtre angulaire 80 Fl hertz fréquence de coupure basse du filtre passe bande 5000 F2 hertz fréquence de coupure haute du filtre passe bande 20000 n" _ nombre de coefficient de Fourier nécessaires à la 1 reconstruction du signal L, - matrice de l'observateur périodique 1 ù 2 w, hertz pulsation du premier harmonique 7500 Tableau 1 : Récapitulatif des paramètres de réglages adaptés au moteur d'essai. Vérification de la qualité de l'estimation (VAL) Afin d'éviter l'envoi de paramètres erronés au contrôle moteur, une étape optionnelle permet de valider la pertinence de la reconstruction du signal vibratoire. Un seuil sur le carré du signal d'erreur (différence entre la mesure et l'estimée) assure cette fonction. Si ce seuil e est trop important, supérieur à un seuil S fixé, l'envoi de paramètres 10 au contrôle moteur n'est pas réalisé (NS). Sinon, on réalise l'étape 4 (CAL). Étape 4 : Détermination de paramètres pour qualifier directement l'intensité de cliquetis (CAL) Dès lors qu'est validée l'estimation des composantes fréquentielles du signal 15 vibratoire, l'opération de contrôle moteur peut être menée. Le contrôle bouclé en régime transitoire de la combustion requière l'estimation en temps réel de paramètre corrélé à l'intensité de cliquetis, afin d'atteindre un contrôle cycle à cycle. Le principe de la méthode repose sur l'estimation en temps réel d'un paramètre, noté Ic,, permettant de caractériser l'intensité du cliquetis dans le cylindre d'un moteur. 20 De façon préférentielle, le calcul de ce paramètre repose sur l'analyse de l'énergie contenue dans les harmoniques (composantes spectrales) sélectionnées. L'énergie E(a) est donnée par la somme des nh premiers coefficients de la décomposition de Fourier délivrés par l'observateur périodique et mis au carré : E(a) = E (ck (a))2 En se basant sur cette énergie, on donne un sens physique au paramètre Ic,, puisque l'énergie E(a) augmente avec l'énergie libérée pendant la combustion. On peut par exemple définir le paramètre suivant, dont la valeur est représentative de l'intensité du cliquetis : IC~ = mai~[CE[[Cl,a2i(E(a))) '/z 10 Étape 5 : Contrôle moteur de la combustion (CTRL) La mise en oeuvre d'une stratégie de commande en boucle fermée est possible en se basant sur ce paramètre Ic,. Ce paramètre d'intensité est utilisé pour compenser des variations des caractéristiques du carburant ou des déséquilibres cylindre à cylindre. A titre d'exemple, voici une stratégie de contrôle à partir du paramètre Ic,. 15 Réglage du phasage de l'allumage. Le phasage de l'allumage est une composante primordiale car elle entre en ligne de compte dans l'optimisation du rendement de la combustion et des émissions polluantes. Dans ce cas précis, le Ic, est un candidat idéal pour envisager un bouclage de l'avance à l'allumage afin de ne pas dépasser une intensité impliquant la destruction ou la détérioration moteur. 20 En effet, plus on retarde l'allumage, plus la probabilité de cliqueter devient importante. De ce fait, on contrôle l'avance à l'allumage avec un contrôleur proportionnel intégral, afin de contrôler le cliquetis vers un niveau de référence qui permet une bonne production de couple et de ne pas endommager le moteur. Cette valeur de référence peut être déterminée expérimentalement car elle dépend du 25 moteur et de la géométrie du cylindre. Variantes On a décrit l'invention en utilisant un capteur de pression cylindre. On ne sortirait pas du cadre de l'invention en utilisant tout autre capteur capable de fournir un signal 30 corrélé avec les vibrations du cylindre. On peut par exemple utiliser le signal issu d'un capteur accéléromètrique. En effet, un tel signal correspond à une grandeur oscillante périodique et peut donc être considéré comme la somme de composantes sinusoïdales dont les fréquences sont multiples entiers de la fréquence fondamentale. Ce signal d'accélération a(t) correspond à la vibration du bloc moteur qui est fortement corrélé à la combustion. En appliquant le même traitement qu'au signal de pression cylindre, on peut estimer l'intensité du cliquetis au moyen de ICI. Avantage de l'invention. La méthode proposée présente l'avantage, contrairement aux méthodes antérieures, d'accéder en temps réel à une estimation de l'intensité du cliquetis. Ainsi, à l'issue d'un cycle moteur, le paramètre de cliquetis (Ic,) est disponible, et permet de réaliser un diagnostic du cycle écoulé. En fonction de ce diagnostic, si l'intensité est jugée trop forte, les paramètres d'allumage sont modifiés afin de garantir la stabilité de la combustion suivante et sa conformité par rapport au point de référence. On rappelle que l'analyse fréquentielle d'un signal peut se révéler coûteuse en terme de temps de calcul et inadéquate pour le traitement en temps réel. Ici, la méthode proposée est à la fois précise et rapide. De plus, la méthode selon l'invention ne dépend pas de la nature de la combustion ou de la technologie du capteur utilisé. Les paramètres de réglages assurent une bonne adaptabilité de la méthode à diverses configurations moteur/capteur. REVENDICATIONS1. Méthode pour estimer l'intensité de cliquetis d'un moteur à combustion interne comprenant au moins un capteur vibratoire, à partir duquel on acquiert en continu et en fonction de l'angle vilebrequin un signal vibratoire représentatif des vibrations dudit moteur, caractérisée en ce que l'on réalise les étapes suivantes : on réalise en temps réel un filtrage dudit signal vibratoire à l'aide d'un filtre ; on construit un modèle de la dynamique dudit filtre que l'on couple avec un estimateur ; - on détermine en temps réel des coefficients de la décomposition de Fourier dudit signal vibratoire acquis, en inversant la dynamique dudit filtre à l'aide dudit estimateur ; on déduit en temps réel un paramètre corrélé à l'intensité de cliquetis à partir desdits coefficients. 2. Méthode selon la revendication 1, dans laquelle ledit signal vibratoire est un signal issu d'un capteur de pression cylindre ou un signal issu d'un capteur accéléromètrique. 3. Méthode selon l'une des revendications précédentes, dans laquelle ledit paramètre est fonction de l'énergie contenue dans ledit signal, et ladite énergie est calculée en sommant les carrés desdits coefficients de la décomposition de Fourier. 4. Méthode selon la revendication 3, dans laquelle ledit paramètre est égale à la racine carrée du maximum de ladite énergie. 5. Méthode selon l'une des revendications précédentes, dans laquelle ledit signal vibratoire comporte plusieurs harmoniques et l'on estime lesdits coefficients de la décomposition de Fourier uniquement pour l'harmonique ou les harmoniques corrélées avec la combustion. 6. Méthode selon la revendication 5, dans laquelle on sélectionne lesdites harmoniques corrélées avec la combustion en réalisant une analyse spectrale ou une analyse temps/fréquence, ou en choisissant la bande de fréquence 5000Hz - 20000Hz. 7. Méthode selon l'une des revendications précédentes, dans laquelle on détermine une plage angulaire incluant l'intégralité de la réponse vibratoire liée à lacombustion, et l'on restreint ledit signal vibratoire à cette plage angulaire avant de réaliser ledit filtrage. 8. Méthode selon l'une des revendications précédentes, dans laquelle ledit filtre est un filtre passe-bande. 9. Méthode selon l'une des revendications précédentes, dans laquelle ledit estimateur est un estimateur non linéaire de type adaptatif. 10. Méthode selon l'une des revendications précédentes, dans laquelle on modifie des paramètres moteurs en fonction dudit paramètre corrélé à l'intensité de cliquetis, pour contrôler la combustion dudit moteur.