FR 3086391 A1 20200327 FR 1858512 A 20180920 La présente invention concerne le domaine de la détermination d'un indicateur du cliquetis pour un moteur à combustion interne, en particulier un moteur à allumage commandé. Ce type de moteur comprend au moins un cylindre comportant une chambre de 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 la description. Dans un moteur à explosion, la phase de 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 est parfois suffisamment importante pour que le mélange non brûlé coincé contre les parois atteigne son point d'auto-allumage et s'auto- enflamme en un ou 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 qui résulte de l'apparition d'ondes de pression dans la chambre de combustion. Ces explosions parasites produisent des vibrations dans le domaine audible et au-delà (de l'ordre de 5 à 50 KHz). Elles sont très vives et peuvent rapidement conduire à un échauffement des parois. L'accumulation de cliquetis endommage le métal du piston et/ou des parois du cylindre et les segments. Au bout de quelque temps (selon l'intensité) cela conduit à la destruction du piston, des segments ou des parois du cylindre. L'estimation du cliquetis permet un contrôle de la combustion limitant l'effet du cliquetis. État de la technique Généralement, la phase de détection du cliquetis comporte une phase d'acquisition de signaux, puis une phase de traitement de ces signaux permettant de détecter l'apparition du cliquetis, de le caractériser et/ou de le quantifier. 2 Une méthode répandue consiste en une mesure de la pression dans le cylindre, suivie par une détermination d'un indicateur nommé MAPO (de l'anglais « Maximum Amplitude of Pressure Oscillations » pouvant être traduit par amplitude maximum des oscillations de pression). 5 D'autres indicateurs ont été développés, ces indicateurs sont généralement basés sur une mesure de la pression locale dans le cylindre ou sur la mesure au moyen d'un accéléromètre. Selon un premier exemple, la demande de brevet JP2017218936 décrit une détermination du cliquetis sur la base du taux de production de chaleur dans le cylindre. 10 Conformément à un deuxième exemple, les demandes de brevet JP2017207015, JP2017214917, WO17126304 et US2017226956 décrivent la détection du cliquetis par une mesure de la pression dans le cylindre. En particulier, la demande de brevet JP2017214917 décrit la détection du cliquetis par une mesure de la pression dans le cylindre, et par la considération de différents modes de 15 vibrations, ainsi que par le calcul d'un ratio des amplitudes entre les modes de vibrations. De plus, la demande de brevet W017126304 divulgue la décomposition du signal par fenêtrage temporel, et détection du cliquetis via une comparaison des pics dans chaque fenêtre. Les documents listés ci-dessous répertorient et comparent différents indicateurs de 20 cliquetis, différents types de capteurs, et leurs performances: Dues, S.; Adams, J.; Shinkle, G.; Div, D. Combustion Knock Sensing : Sensor Selection and Application Issues SAE Trans. J. Engines, 1990. Gao, X.; Stone, R. ; Hudson, C. ; Bradbury, I. The Detection and Quantification of Knock in Spark Ignition Engines SAE Technical Paper, 1993. 25 Naber, J.; Blough, J.; Frankowski, D.; et al. Analysis of Combustion Knock Metrics in Spark-Ignition Engines SAE Trans. J. Engines, 2006. Shahlari, A.; Ghandhi, J. A Comparison of Engine Knock Metrics SAE Technical Paper, 2012. 30 Toutefois, les limites suivantes sont rencontrées par ce type d'indicateurs de cliquetis : ces indicateurs de cliquetis sont dépendants de la position du capteur de pression dans la chambre de combustion, car la pression n'est pas homogène dans le cylindre. En effet, la mesure d'un indicateur de cliquetis à des positions différentes 3 génèrent différentes évaluations de l'intensité du cliquetis, ceci est particulièrement vrai pour l'indicateur MAPO, et pour les indicateurs de cliquetis qui sont fondés sur l'amplitude des oscillations de la pression locale mesurée, le choix de la technique de traitement du signal influe sur la détection du cliquetis. 5 En effet, deux filtres ayant des gammes de fréquences différentes génèrent différentes évaluations du cliquetis, ces indicateurs ne représentent pas le risque d'endommagement du cylindre, en particulier, l'indicateur MAPO est purement qualitatif et non représentatif d'un endommagement du cylindre. 10 Afin de pallier ces inconvénients, la présente invention concerne un procédé de détermination d'un indicateur de cliquetis de la combustion dans un cylindre d'un moteur à combustion interne. Le procédé est mis en oeuvre au moyen d'une mesure de pression locale, dont la décomposition modale permet d'avoir des caractéristiques de la pression 15 globale (représentation de la pression dans toute la chambre de combustion) dans le cylindre. Cette pression permet de déterminer, directement ou indirectement, un indicateur de cliquetis représentatif de l'endommagement du cylindre. L'invention concerne en outre un procédé de contrôle d'un moteur à combustion interne mettant en oeuvre un tel procédé de détermination d'un indicateur de cliquetis. 20 Le procédé selon l'invention La présente invention concerne un procédé de détermination d'un indicateur de cliquetis de la combustion dans au moins un cylindre d'un moteur à combustion interne, notamment à allumage commandé, ledit cylindre étant pourvu d'un capteur de pression. Pour ce procédé, 25 on met en oeuvre les étapes suivantes : a) on mesure une pression locale dans ledit cylindre en un point de mesure dudit cylindre au moyen dudit capteur de pression ; b) on réalise une décomposition modale dudit signal de mesure de ladite pression locale pour déterminer des modes d'oscillation de ladite pression locale ; 30 c) on détermine une pression globale dans ledit cylindre au moyen de ladite décomposition modale dudit signal de mesure de ladite pression locale ; et d) on détermine un indicateur de cliquetis, ledit indicateur de cliquetis étant directement ou indirectement fonction de ladite pression globale dans le cylindre. 4 Selon un mode de réalisation, on détermine ladite pression globale, en mettant en oeuvre les étapes suivantes : i) On définit un modèle de ladite pression globale fonction de l'amplitude des modes de ladite pression globale ; 5 ii) on détermine ladite amplitude desdits modes de ladite pression globale, par identification des oscillations de chaque mode de ladite décomposition modale de ladite pression locale au point de mesure de ladite pression locale. Conformément à une mise en oeuvre, procédé comporte une étape supplémentaire de 10 détermination de la vitesse des gaz dans ledit cylindre, au moins à proximité des parois dudit cylindre, en appliquant un modèle acoustique de la conservation de mouvement de gaz dans ledit cylindre à ladite pression globale dans ledit cylindre, et on détermine ledit indicateur de cliquetis en fonction de ladite vitesse des gaz dans ledit cylindre. 15 Selon une caractéristique, le procédé comporte en outre les étapes suivantes : i) on détermine une vitesse des gaz dans ledit cylindre, au moins à proximité des parois dudit cylindre, en appliquant un modèle acoustique de la conservation de mouvement de gaz dans ledit cylindre à ladite pression globale dans ledit cylindre 20 ii) on détermine un coefficient de transfert thermique auxdites parois dudit cylindre en appliquant un modèle de coefficient de transfert thermique à l'aide de ladite pression locale mesurée et de ladite vitesse des gaz dans ledit cylindre ; et iii) on détermine un indicateur de cliquetis au moyen dudit coefficient de transfert thermique convectif aux parois dudit cylindre. 25 Avantageusement, ledit modèle de coefficient de transfert thermique s'écrit : h= C x Pa x uP x Dy x T5 avec h le coefficient de transfert thermique, C une constante, D une longueur caractéristique de l'écoulement dudit cylindre, P la pression dans ledit cylindre, u ladite vitesse des gaz à proximité des parois du cylindre, T la température des gaz dans le 30 cylindre, a ,)3,y, et S des coefficients constants prédéterminés. De préférence, ladite température T est déterminée par une application de la loi des gaz parfaits dans ledit cylindre (12) Selon un aspect, on détermine ledit indicateur de cliquetis au moyen d'une étape de détermination du transfert thermique Q, de telle sorte que Q = h x (T1 - T2) avec h le 5 coefficient de transfert thermique, et T1 - T2 le gradient de température entre les gaz dans le cylindre et la paroi du cylindre. De manière avantageuse, ledit modèle acoustique de la conservation du mouvement 5 de gaz dans ledit cylindre s'écrit : su po = st +VP1 = 0 avec po la densité des gaz, u1 la vitesse des gaz associée à l'onde acoustique et P1 la pression dans ledit cylindre associée à l'onde acoustique. Selon une mise en oeuvre, on détermine ladite vitesse des gaz dans ledit cylindre en 10 prenant en compte la vitesse nominale des gaz dans ledit cylindre. Conformément à un mode de réalisation, ladite vitesse nominale des gaz dans ledit cylindre est estimée au moyen de caractéristiques de l'admission des gaz au sein dudit cylindre, notamment les caractéristiques de l'admission générant des tourbillons dudit gaz 15 dans ledit cylindre. Selon un aspect, ledit procédé comporte une étape supplémentaire de détermination d'une combustion anormale par analyse dudit indicateur de cliquetis. 20 En outre, l'invention concerne un procédé de contrôle de la combustion d'un moteur à combustion interne, notamment à allumage commandé, dans lequel on détecte un cliquetis de la combustion dans au moins un cylindre dudit moteur à combustion interne. Pour ce procédé, on met en oeuvre les étapes suivantes : a) on détermine un indicateur de cliquetis au moyen du procédé selon l'une des 25 caractéristiques précédentes ; et b) on contrôle la combustion dudit moteur à combustion interne en fonction dudit indicateur de cliquetis. Présentation succincte des figures 30 D'autres caractéristiques et avantages du procédé 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. La figure 1 illustre un moteur à combustion interne adapté à la mise en oeuvre du procédé selon l'invention. 6 La figure 2 illustre les étapes du procédé selon un premier mode de réalisation de l'invention. La figure 3 illustre les étapes du procédé selon un deuxième mode de réalisation de l'invention. 5 La figure 4 illustre les étapes du procédé selon un troisième mode de réalisation de l'invention. Description détaillée de l'invention La présente invention concerne un procédé de détermination d'un indicateur de 10 cliquetis de la combustion dans un cylindre d'un moteur à combustion interne. La présente invention est particulièrement adaptée pour un moteur à allumage commandé (moteur essence). Le moteur à combustion interne peut être de tout type : il peut notamment être suralimenté ou non, à injection de carburant directe ou indirecte, et équipé d'un système de 15 recirculation des gaz d'échappement ou non, etc. La figure 1 illustre, schématiquement et de manière non limitative, un moteur à combustion interne 10 adapté à la mise en oeuvre du procédé selon l'invention. Le moteur à combustion interne 10 illustré est de type suralimenté à allumage commandé, en particulier 20 de type essence, comprend au moins un cylindre 12 (quatre cylindres 12 pour le mode de réalisation de la figure 1) avec une chambre de combustion 14 à l'intérieur de laquelle se produit la combustion d'un mélange d'air suralimenté et de carburant (caractéristiques optionnelles pour la mise en oeuvre du procédé selon l'invention). Le cylindre 12 comprend au moins un moyen d'alimentation en carburant sous 25 pression 16, par exemple sous la forme d'un injecteur de carburant 18 contrôlé par une vanne 20, qui débouche dans la chambre de combustion, au moins un moyen d'admission d'air 22 avec une soupape d'admission 24 associée à une tubulure d'admission 26 se terminant par un plénum 26b (non représenté sur la figure), au moins un moyen d'échappement des gaz brûlés 28 avec une soupape d'échappement 30 et une tubulure 30 d'échappement 32 et au moins un moyen d'allumage 34, comme une bougie, qui permet de générer une ou plusieurs étincelles permettant d'enflammer le mélange carburé présent dans la chambre de combustion. Les tubulures 32 des moyens d'échappement 28 de ce moteur sont raccordées à un collecteur d'échappement 36 lui-même connecté à une ligne d'échappement 38. Un dispositif 35 de suralimentation 40, par exemple un turbocompresseur ou un compresseur volumétrique, 7 est placé sur cette ligne d'échappement et comprend un étage d'entraînement 42 avec une turbine balayée par les gaz d'échappement circulant dans la ligne d'échappement et un étage de compression 44 qui permet de faire admettre un air d'admission sous pression dans les chambres de combustion 14 par les tubulures d'admission 26. 5 Le moteur comprend des moyens 46a de mesure de la pression locale dans le cylindre, disposés au sein même du cylindre 12 du moteur. Ces moyens de mesure comprennent généralement un capteur de pression qui permet de générer un signal représentatif de l'évolution de la pression locale dans un cylindre ou un accéléromètre permettant d'obtenir une représentation indirecte des oscillations de pression dans la 10 chambre de combustion. Le moteur peut également comporter des moyens 46b de mesure de la pression d'admission (caractéristique optionnelle pour le procédé selon l'invention), disposés dans le plénum 26b. Ces moyens de mesure sont généralement constitués par un capteur de pression absolue, de type piézoélectrique, qui permet de générer un signal représentatif de 15 l'évolution de la pression d'admission dans le plénum d'admission. Le moteur comprend également une unité de calcul et de commande 48, dénommée calculateur moteur, qui est reliée par des conducteurs (pour certains bidirectionnels) aux différents organes et capteurs du moteur de façon à pouvoir recevoir les différents signaux émis par ces capteurs, et ensuite commander les organes de ce moteur pour assurer son 20 bon fonctionnement. Ainsi, dans le cas de l'exemple montré à la figure 1, les bougies 34 sont reliées par des conducteurs 50 au calculateur moteur 48 de façon à commander le moment de l'allumage du mélange carburé, le capteur de pression cylindre 46a est connecté par une ligne 52 à ce même calculateur moteur pour lui envoyer les signaux représentatifs de l'évolution de la 25 pression dans le cylindre, et les vannes 20 de commande des injecteurs 18, sont raccordées par des conducteurs 54 au calculateur 48 pour commander l'injection de carburant dans les chambres de combustion. Les moyens 46b sont également connectés par une ligne 53 au calculateur moteur 48. 30 Pour la mise en oeuvre du procédé selon l'invention, au moins un cylindre du moteur à combustion interne comporte un capteur de pression, mesurant une pression locale dans le cylindre. Il s'agit de la pression des gaz dans le cylindre. Les gaz forment un comburant ou un mélange carburé (cas de l'injection indirecte), et peuvent comprendre notamment de l'air à pression ambiante, de l'air suralimenté, un mélange d'air (suralimenté ou non) et des gaz 35 brûlés. 8 Selon un premier mode de réalisation l'invention et tel qu'illustré en figure 2, le procédé de détermination d'un indicateur de cliquetis peut comporter les étapes suivantes : 1) mesure de la pression locale dans le cylindre (MES) 5 2) décomposition modale du signal de mesure de la pression locale (DEC) 3) détermination de la pression globale dans le cylindre (PGL) 6) détermination de l'indicateur de cliquetis (IND) 7) contrôle de la combustion (CON) 10 Selon un deuxième mode de réalisation de l'invention et tel qu'illustré en figure 3, le procédé de détermination d'un indicateur de cliquetis peut comporter les étapes suivantes : 1) mesure de la pression locale dans le cylindre (MES) 2) décomposition modale du signal de mesure de la pression locale (DEC) 3) détermination de la pression globale dans le cylindre (PGL) 15 4) détermination de la vitesse des gaz dans le cylindre (Vg) 6) détermination de l'indicateur de cliquetis (IND) 7) contrôle de la combustion (CON) Selon un troisième mode de réalisation de l'invention et tel qu'illustré en figure 4, le 20 procédé de détermination d'un indicateur de cliquetis peut comporter les étapes suivantes : 1) mesure de la pression locale dans le cylindre (MES) 2) décomposition modale du signal de mesure de la pression locale (DEC) 3) détermination de la pression globale dans le cylindre (PGL) 4) détermination de la vitesse des gaz dans le cylindre (Vg) 25 5) détermination du transfert thermique (TTH) 6) détermination de l'indicateur de cliquetis (IND) 7) contrôle de la combustion (CON) Les étapes 1) à 3) sont communes pour les trois modes de réalisation de l'invention. 30 Les étapes 2) à 6) peuvent être mises en oeuvre par le calculateur du moteur ; ainsi la détermination de l'indicateur du cliquetis est réalisée en ligne. L'étape facultative 7) représente une variante de réalisation de l'invention. Toutes les étapes du procédé sont réalisées, de préférence, en temps réel (c'est-à-dire 35 de cycle à cycle), de manière à déterminer l'indicateur de cliquetis en temps réel, afin de pouvoir éventuellement agir sur le cliquetis avant le prochain allumage. 9 Dans la suite de la demande, on appelle « pression locale » la pression au point de mesure du capteur de pression, et on appelle « pression globale » une représentation de la pression dans toute la chambre de combustion du cylindre. 5 Les sept étapes du procédé sont détaillées dans la suite de la description. 1) Mesure de la pression locale dans le cylindre Il s'agit d'une étape de mesure de la pression locale dans la chambre de combustion du cylindre au moyen du capteur de pression équipant le cylindre. Le capteur de pression 10 génère un signal de capteur représentatif de la pression locale en un point dans le cylindre. 2) Décomposition modale du signal de mesure de la pression locale Cette étape concerne la décomposition du signal de mesure de la pression locale, afin de déterminer les modes d'oscillation de pression dans la chambre de combustion. 15 Selon un mode de réalisation de l'invention, cette décomposition modale du signal de mesure peut être mise en oeuvre par un filtrage du signal de mesure de la pression locale autour des fréquences d'intérêt correspondant aux modes d'oscillation. Les fréquences d'intérêt peuvent être obtenues par la résolution de l'équation d'onde acoustique pour la géométrie spécifique du cylindre. Selon un exemple de réalisation, 1 8t22 20 l'équation d'onde acoustique peut s'écrire : OP = X avec P la pression, c la vitesse c2 8P du son. Cette équation aux dérivées partielles peut être résolue par des solutions de Bessel, qui contiennent les fréquences recherchées, ainsi que la forme d'onde dans la chambre. Selon un mode de réalisation de l'invention, la vitesse du son peut être calculée par une formule du type : c = yrT avec y le coefficient adiabatique, r la constante 25 thermodynamique spécifique des gaz et T la température. 3) Détermination de la pression globale dans le cylindre Lors de cette étape, on détermine des caractéristiques de la pression globale dans le cylindre au moyen de la décomposition modale de la pression locale dans le cylindre. 30 Une fois que les oscillations (les modes) de la pression locale sont décomposées mode par mode, la mesure de pression locale est associée avec l'oscillation de pression globale (mode par mode) de toute la chambre de combustion du cylindre. Pour cela, la somme des oscillations (modes) de chaque mode de la pression locale permet de reconstruire les oscillations (modes) de la pression globale, par exemple au moyen des solutions de Bessel. 35 Ainsi, la mesure de pression locale est exploitée pour retrouver la description complète de la pression à l'intérieur de la chambre de combustion. 10 Cette détermination de la pression globale est permise car les oscillations de pression dans la chambre de combustion peuvent être complètement décrites au moyen d'une décomposition modale. Ainsi, grâce à seulement un capteur de pression locale, il est possible de remonter à la 5 pression globale dans le cylindre. Selon un mode de réalisation de l'invention, on peut mettre en oeuvre les étapes suivantes pour réaliser la détermination de la pression globale ; a) au moyen de la résolution de l'équation d'onde à l'étape 1), on construit un modèle de la pression globale P par une équation du type : P(x, y, z, t) = Xi(At(t) x 10 h(x, y, z) x g i(0) avec A l'amplitude qui dépend du cycle analysé, les fonctions f et g définies pour chaque mode i (ce sont les solutions de l'équation d'onde de l'étape 2), t le temps et x, y, z les coordonnées d'un point dans la chambre de combustion, b) on calcule le modèle de pression globale aux coordonnées du capteur (x0,y0,z0), pour identifier les fonctions Ai(t) par corrélation avec la mesure. Cette identification 15 est possible grâce à la décomposition modale permettant de corréler modèle et mesure mode par mode. Cette identification permet de complètement définir le modèle de pression globale. Pour le premier mode de réalisation (figure 2), le procédé continue directement avec 20 l'étape 6). Pour les deuxième et troisième modes de réalisations (respectivement figures 3 et 4)), le procédé continue avec l'étape 4. 4) Détermination de la vitesse des gaz dans le cylindre Lors de cette étape, on détermine la vitesse des gaz dans au moins une portion du 25 cylindre, notamment au moins à proximité des parois latérales du cylindre et/ou du piston et/ou de la culasse. Alternativement, on peut déterminer la vitesse des gaz dans la totalité du cylindre. Cette étape est mise en oeuvre au moyen d'un modèle acoustique de la conservation de quantité de mouvement des gaz. Le modèle acoustique est appliqué à la pression globale dans le cylindre déterminée à l'étape précédente. Ainsi, il est possible de 30 déterminer de manière représentative le comportement physique des gaz dans le cylindre. Selon un mode de réalisation de l'invention, le modèle acoustique de la conservation de mouvement de gaz dans le cylindre peut s'écrire par une équation de la forme po su, + VPl = 0 avec po la densité des gaz, u1 la vitesse des gaz associée à l'onde acoustique et P1 35 la pression dans le cylindre associée à l'onde acoustique. Cette équation acoustique de la conservation de mouvement de gaz permet de calculer la vitesse des gaz d'après le gradient 11 de pression dans la chambre de combustion. Une décomposition de la pression globale en termes spatial et temporel permet d'exploiter l'équation acoustique. Grâce à la décomposition modale de la pression globale (calculée dans les étapes 2 et 3), la contribution des ondes acoustiques à la vitesse des gaz dans la chambre de combustion du 5 cylindre peut être directement calculée. Conformément à une mise en oeuvre de l'invention, la combinaison de la vitesse nominale avec les oscillations de vitesse calculées par le modèle acoustique peut permettre de reconstruire la vitesse totale des gaz dans la chambre. La vitesse nominale peut résulter 10 des caractéristiques de l'admission des gaz dans le cylindre. Ces caractéristiques peuvent être des caractéristiques de l'admission permettant le swirl, le tumble ou le swumble (combinaison de swirl et de tumble). Le swirl, qui est un mouvement macroscopique de rotation du mélange carburé autour d'un axe colinéaire à l'axe du cylindre, se caractérise par une bonne conservation au cours 15 du processus d'admission et plus spécifiquement au cours de la remontée du piston. C'est un mouvement macroscopique aérodynamique qui est généralement utilisé pour les moteurs à combustion interne à allumage par compression où il est un bon moyen d'homogénéiser le mélange carburé. Le tumble est lui aussi un mouvement macroscopique de rotation du mélange carburé 20 mais autour d'un axe globalement perpendiculaire à l'axe du cylindre. Il a la particularité de se transformer en mouvements aérodynamiques microscopiques qui créent de la turbulence lors de la remontée du piston. C'est un mouvement macroscopique aérodynamique qui est généralement utilisé pour les moteurs à combustion interne à allumage commandé où il est un bon moyen d'obtenir une vitesse de combustion appropriée. Par ailleurs, ce mouvement 25 est assez sensible à la géométrie de la chambre de combustion ainsi qu'à la loi de levée, aussi bien en termes d'étalement que de hauteur de levée maximale. L'utilisation du swumble permet de bénéficier des avantages des deux structures aérodynamiques détaillées ci-dessus et donc de bénéficier d'une excellente homogénéisation et d'une meilleure vitesse de combustion grâce à un niveau de turbulence 30 plus élevée lors de la phase de compression que ce qu'on observe avec les meilleurs moteurs à allumage commandé actuels. Pour le deuxième mode de réalisation (figure 3), le procédé continue avec l'étape 6). Pour le troisième mode de réalisation (figure 4), le procédé continue avec l'étape 5). 35 5) Détermination du transfert thermique aux parois du cylindre 12 Lors de cette étape, on détermine le transfert thermique convectif aux parois latérales du cylindre et/ou du piston et/ou de la culasse, en appliquant un modèle de transfert thermique à la pression locale et à la vitesse des gaz dans la chambre de combustion déterminée lors de l'étape 4). Ainsi, en sortie de cette étape, on obtient un coefficient de 5 transfert thermique aux parois du cylindre, représentatif des phénomènes physiques mis en oeuvre dans le cylindre, qui peuvent générer un endommagement du cylindre. Selon un aspect de l'invention, le modèle de transfert thermique peut dépendre de la pression locale, de la température des gaz dans le cylindre, et de la vitesse des gaz à 10 proximité des parois du cylindre. Selon une mise en oeuvre de réalisation de l'invention, le modèle de transfert thermique peut s'écrire : h = C x Pa x ufl x DY x T8 avec h le coefficient de transfert thermique, C une constante, D une longueur caractéristique de l'écoulement dudit cylindre, P la pression dans le cylindre déterminée à l'étape 3), u la vitesse des gaz à proximité des 15 parois du cylindre déterminée à l'étape 4), T la température des gaz dans le cylindre, a y, et â des coefficients constants prédéterminés. Par exemple, les coefficients utilisés peuvent être a = 0.8, fi = 0.8, y = -0.2 et 8 = -0.53 avec D l'alésage cylindre (tels qu'utilisés dans le modèle de Woschni). Conformément à une mise en oeuvre de ce mode de réalisation de l'invention, la 20 température n'est pas mesurée, mais peut être déterminée par une application de la loi des gaz parfaits dans la chambre de combustion du cylindre. La loi des gaz parfaits peut s'écrire : PV = mrT avec P la pression locale mesurée, V le volume du cylindre (connu), m la masse des gaz dans la chambre de combustion, r la constante thermodynamique spécifique des gaz et T la température. Cette mise en oeuvre permet de se passer d'un capteur de 25 température dans le cylindre. Pour cette mise en oeuvre, la masse de gaz peut être estimée par d'autres capteurs externes au cylindre. Alternativement, la température peut être mesurée au moyen d'un capteur de température. 30 6) Détermination d'un indicateur de cliquetis Lors de cette étape, on déduit un indicateur de cliquetis de la grandeur mesurée à l'étape précédente (pression globale pour le premier mode de réalisation, vitesse des gaz pour le deuxième mode de réalisation, et coefficient de transfert thermique pour le troisième 35 mode de réalisation). De cette manière, l'indicateur de cliquetis dépend directement (cas du premier mode de réalisation) ou indirectement (cas des deuxième et troisième modes de 13 réalisation) de la pression globale du cylindre. Ainsi, l'indicateur de cliquetis est indépendant de la position du capteur de pression, et est indépendant d'un quelconque filtrage du signal. Pour le premier mode de réalisation, l'indicateur de cliquetis dépend directement de la 5 pression globale du cylindre. Dans ce cas, l'indicateur de cliquetis peut notamment être une valeur maximale de l'amplitude d'oscillation de la pression globale, la valeur maximale (ou moyenne) du gradient de la pression, etc. Ces méthodes permettent de mieux qualifier la sévérité du cliquetis en utilisant la globalité des conditions dans toute la chambre. 10 Pour le deuxième mode de réalisation, l'indicateur de cliquetis dépend de la vitesse des gaz qui est fonction de la pression globale du cylindre. Dans ce cas, l'indicateur de cliquetis peut notamment être la valeur maximale de la vitesse des gaz dans le cylindre, etc. La vitesse des gaz informe de l'interaction des gaz brulés avec les parois et, permet ainsi de qualifier la susceptibilité d'endommagement du moteur. 15 Pour le troisième mode de réalisation, l'indicateur de cliquetis est dépendant du coefficient du transfert thermique qui est fonction de la vitesse des gaz, elle-même fonction de la pression globale du cylindre. Dans ce cas, l'indicateur de cliquetis peut être la valeur maximale au cours du temps de la moyenne spatiale (dans le volume du cylindre) du 20 coefficient de transfert thermique h sur les parois. Toutefois, l'indicateur de cliquetis peut prendre différentes formes dépendantes de la valeur du coefficient de transfert thermique h sur les parois, par exemple sa moyenne, sa médiane, etc. On peut aussi réaliser un calcul du transfert thermique aux parois. Dans ce cas, on peut déterminer ledit indicateur de cliquetis au moyen d'une étape de détermination 25 du transfert thermique Q, de telle sorte que Q = h x (T1 - T2) avec h le coefficient de transfert thermique, et Ti - T2 le gradient de température. Ainsi, l'indicateur de cliquetis est corrélé aux risques d'endommagement du cylindre. En effet, ce transfert thermique traduit la sévérité de la combustion au niveau des parois du cylindre, qui sont les éléments qui s'endommagent lors du cliquetis. 30 7) Contrôle de la combustion On rappelle que cette étape est facultative. Il s'agit de contrôler la phase de combustion au sein du cylindre du moteur à combustion interne, en fonction de l'indicateur de cliquetis déterminé à l'étape précédente, afin de réduire l'apparition du cliquetis et donc 35 les risques d'endommagement du cylindre. Le contrôle de la combustion peut être mis en oeuvre par le contrôleur du moteur. 14 En cas de cliquetis détecté au moyen du procédé selon l'invention, le calculateur peut lancer les actions nécessaires au contrôle des combustions suivantes afin d'éviter la poursuite de l'occurrence du cliquetis lors des combustions suivantes. Selon un mode de réalisation, le contrôle peut être fondé sur une comparaison de 5 l'indicateur de cliquetis avec un seuil de cliquetis. Avantageusement, on peut contrôler la phase de combustion par un contrôle de la variation de l'avance à l'allumage en fonction de l'indicateur de cliquetis. La présente invention n'est REVENDICATIONS1) Procédé de détermination d'un indicateur de cliquetis de la combustion dans au moins un cylindre (12) d'un moteur à combustion interne (10), notamment à allumage commandé, ledit cylindre (12) étant pourvu d'un capteur de pression (46a), caractérisé en ce qu'on met en oeuvre les étapes suivantes : a) on mesure une pression locale dans ledit cylindre (12) en un point de mesure dudit cylindre au moyen dudit capteur de pression (46a) ; b) on réalise une décomposition modale dudit signal de mesure de ladite pression locale pour déterminer des modes d'oscillation de ladite pression locale ; c) on détermine une pression globale dans ledit cylindre (12) au moyen de ladite décomposition modale dudit signal de mesure de ladite pression locale ; et d) on détermine un indicateur de cliquetis, ledit indicateur de cliquetis étant directement ou indirectement fonction de ladite pression globale dans le cylindre (12). 2) Procédé selon la revendication 1, dans lequel on détermine ladite pression globale, en mettant en oeuvre les étapes suivantes : i) On définit un modèle de ladite pression globale fonction de l'amplitude des modes de ladite pression globale ; ii) on détermine ladite amplitude desdits modes de ladite pression globale, par identification des oscillations de chaque mode de ladite décomposition modale de ladite pression locale au point de mesure de ladite pression locale. 3) Procédé selon l'une des revendications précédentes, dans lequel le procédé comporte une étape supplémentaire de détermination de la vitesse des gaz dans ledit cylindre (12), au moins à proximité des parois dudit cylindre (12), en appliquant un modèle acoustique de la conservation de mouvement de gaz dans ledit cylindre (12) à ladite pression globale dans ledit cylindre (12), et dans lequel on détermine ledit indicateur de cliquetis en fonction de ladite vitesse des gaz dans ledit cylindre (12). 4) Procédé selon l'une des revendications 1 ou 2, dans lequel le procédé comporte en outre les étapes suivantes : i) on détermine une vitesse des gaz dans ledit cylindre (12), au moins à proximité des parois dudit cylindre (12), en appliquant un modèle acoustique de la conservation de mouvement de gaz dans ledit cylindre (12) à ladite pression globale dans ledit cylindre (12) ; ii) on détermine un coefficient de transfert thermique auxdites parois dudit cylindre (12) en appliquant un modèle de coefficient de transfert thermique à l'aide de 16 ladite pression locale mesurée et de ladite vitesse des gaz dans ledit cylindre (12) ; et iii) on détermine un indicateur de cliquetis au moyen dudit coefficient de transfert thermique convectif aux parois dudit cylindre (12). 5 5) Procédé selon la revendication 4, dans lequel ledit modèle de coefficient de transfert thermique s'écrit h = C x Pa x uP x DY x T8 avec h le coefficient de transfert thermique, C une constante, D une longueur caractéristique de l'écoulement dudit cylindre, P la pression dans ledit cylindre, u ladite vitesse des gaz à proximité des parois du cylindre, T 10 la température des gaz dans le cylindre, a, fl, y, et S des coefficients constants prédéterminés. 6) Procédé selon la revendication 5, dans lequel ladite température T est déterminée par une application de la loi des gaz parfaits dans ledit cylindre (12) 15 7) Procédé selon l'une des revendications 4 à 6, dans lequel on détermine ledit indicateur de cliquetis au moyen d'une étape de détermination du transfert thermique Q, de telle sorte que Q = h x (T1 - T2) avec h le coefficient de transfert thermique, et T1 - T2 le gradient de température entre les gaz dans le cylindre (12) et la paroi du cylindre (12). 20 8) Procédé selon l'une des revendications 3 à 7, dans lequel ledit modèle acoustique de la conservation du mouvement de gaz dans ledit cylindre s'écrit : 8u po= +VP1 = 0 avec po la densité des gaz, u1 la vitesse des gaz associée à l'onde acoustique et P1 la pression dans ledit cylindre associée à l'onde acoustique. 25 9) Procédé selon l'une des revendications 3 à 8, dans lequel on détermine ladite vitesse des gaz dans ledit cylindre (12) en prenant en compte la vitesse nominale des gaz dans ledit cylindre (12). 30 10) Procédé selon la revendication 9, dans lequel ladite vitesse nominale des gaz dans ledit cylindre (12) est estimée au moyen de caractéristiques de l'admission des gaz au sein dudit cylindre (12), notamment les caractéristiques de l'admission générant des tourbillons dudit gaz dans ledit cylindre (12). 35 11) Procédé selon l'une des revendications précédentes, dans lequel ledit procédé comporte une étape supplémentaire de détermination d'une combustion anormale par analyse dudit indicateur de cliquetis. 17 12) Procédé de contrôle de la combustion d'un moteur à combustion interne (10), notamment à allumage commandé, dans lequel on détecte un cliquetis de la combustion dans au moins un cylindre (12) dudit moteur à combustion interne (10), caractérisé en ce qu'on 5 met en oeuvre les étapes suivantes : a) on détermine un indicateur de cliquetis au moyen du procédé selon l'une des revendications précédentes ; et b) on contrôle la combustion dudit moteur à combustion interne (10) en fonction dudit indicateur de cliquetis. 10