FR 2952679 A1 20110520 FR 0905463 A 20091113 10 La présente invention concerne le domaine du contrôle de la phase de combustion d'un moteur à combustion interne. Notamment la présente invention se rapporte à une méthode pour détecter une combustion anormale, du type pré allumage à bas régime et à forte charge, dans une chambre de combustion d'un tel moteur. Elle concerne plus particulièrement mais non exclusivement une telle méthode 15 appliquée à un moteur à allumage commandé « downsizé », fonctionnant sous très fortes charges. Les moteurs à allumage commandé ont l'avantage de limiter les émissions polluantes locales (HC, CO et NON) grâce à l'excellente adéquation entre le mode de fonctionnement (à richesse 1) et leur système de post-traitement simple et à bas coût. Malgré cet avantage 20 fondamental, ces moteurs sont mal positionnés en termes d'émissions de gaz à effet de serre, car les moteurs Diesel qui leur font concurrence peuvent atteindre des émissions de CO2 inférieures de 20 % en moyenne. La combinaison du "downsizing" et de la suralimentation est une des solutions qui se répand de plus en plus pour abaisser la consommation des moteurs à allumage commandé. 25 Malheureusement, le mécanisme classique de combustion dans ces moteurs peut être perturbé par des combustions anormales. 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 30 compression, puis une étape de combustion sous l'effet d'un allumage commandé, par une bougie. Ces étapes sont regroupées sous le vocable de « phase de combustion » dans la suite de la description. 1 Il a pu être constaté que ce mélange carburé peut subir différents types de combustion, et que ces types de combustion sont la source de différents niveaux de pression, ainsi que de contraintes mécaniques et/ou thermiques, dont certaines peuvent endommager gravement le moteur. La première combustion, dite combustion conventionnelle ou combustion normale, est le résultat de la propagation de la combustion d'un mélange carburé comprimé lors d'une étape préalable de compression du moteur. Cette combustion se propage normalement selon un front de flamme à partir de l'étincelle générée à la bougie, et ne risque pas de détériorer le moteur. Un autre type de combustion est une combustion avec cliquetis, qui résulte d'une auto-inflammation indésirable dans la chambre de combustion. Ainsi, après l'étape de compression du mélange carburé, la bougie est actionnée pour permettre l'allumage de ce mélange carburé. Sous l'effet de la pression générée par le piston, et de la chaleur dégagée par le début de la combustion du mélange carburé, il se produit une auto-inflammation brutale et localisée d'une partie du mélange carburé comprimé, avant que n'arrive le front de flamme issu de l'allumage du mélange carburé par la bougie. Ce mécanisme, dénommé cliquetis, conduit à une augmentation locale de la pression et de la température et peut engendrer, en cas de répétitions, des effets destructifs sur le moteur et principalement au niveau du piston. Enfin, un autre type de combustion est une combustion anormale due à un pré-allumage du mélange carburé avant que la bougie n'initie l'allumage du mélange carburé présent dans la chambre de combustion. Cette combustion anormale affecte en particulier les moteurs qui sont le résultat d'une opération de "miniaturisation", plus connu sous le terme anglais de "downsizing". Cette opération vise à diminuer la taille et/ou la cylindrée du moteur tout en conservant la même puissance et/ou le même couple que des moteurs conventionnels. Généralement, ce type de moteurs est principalement de type essence et est fortement suralimenté. Il a pu être constaté que cette combustion anormale se réalise à fortes charges, et généralement lors des bas régimes de fonctionnement du moteur, lorsque le calage de la combustion du mélange carburé ne peut pas être l'optimum à cause du cliquetis. Compte tenu des fortes pressions et des températures élevées atteintes dans la chambre de combustion par la suralimentation, un démarrage de combustion anormale peut se produire, sporadiquement ou de façon continue, bien avant le moment où se réalise l'allumage du mélange carburé par la bougie. Cette combustion se caractérise par une première phase de propagation de flamme qui est calée trop tôt par rapport à celle d'une combustion conventionnelle. Cette phase de propagation peut être interrompue par une auto-inflammation qui va concerner une grande partie du mélange carburé présent dans la chambre de combustion, beaucoup plus grande que dans le cas du cliquetis. Dans le cas où cette combustion anormale se produit de façon répétitive, de cycle- moteur à cycle-moteur, et se réalise à partir d'un point chaud du cylindre par exemple, celle-ci est dénommée "pré allumage". Si cette combustion se produit de manière violente, aléatoire et sporadique, elle est appelée "claquement" ou "rumble" (« pre-ignition »). Cette dernière combustion anormale entraîne des niveaux de pressions très élevés (120 à 250 bars), ainsi qu'une augmentation des transferts thermiques qui peuvent entraîner une destruction partielle ou totale de l'équipage mobile du moteur, comme le piston ou la bielle. Ce type de pré allumage constitue actuellement une véritable limite à la miniaturisation ("downsizing") des moteurs à allumage commandé. Il s'agit d'un phénomène très complexe qui peut avoir de multiples origines. Plusieurs hypothèses ont été évoquées dans la littérature pour expliquer son apparition mais aucune d'entre elles n'a pour l'instant été clairement validée, il semble même plutôt que plusieurs de ces causes potentielles se manifestent simultanément et interagissent entre elles. Cette interaction, la violence du phénomène et son caractère stochastique rendent son analyse extrêmement compliquée. De plus, les différentes études sur le sujet se heurtent toutes au problème de l'identification même de ces combustions anormales. II est en effet difficile de dire si un moteur est plus sensible qu'un autre au pré allumage tant qu'on est incapable de statuer sur la nature de chacune des combustions au sein d'un échantillon donné. Une méthode permettant de détecter et quantifier ces combustions anormales est donc de première nécessité, car elle permet justement d'établir cette hiérarchie et d'identifier les pistes qui permettent d'améliorer la conception et les réglages des moteurs. Cette opération est particulièrement intéressante lors des développements de moteurs au banc moteur. État de la technique La méthodologie générale de traitement de ces combustions anormales est schématisée sur la figure 1, avec dans un premier temps une phase de prévention (PP) pour limiter au maximum les chances d'apparition du phénomène, puis une phase de détection (PD) lorsque la prévention n'a pas suffit à éviter le phénomène, pour déterminer si oui ou non il y a lieu d'intervenir dans le cycle même où le pré allumage a été détecté au moyen d'une phase corrective (PC). La phase de détection comporte une phase d'acquisition de signaux, puis une phase de traitement de signaux permettant de détecter l'apparition du pré allumage à forte charge, de le caractériser et de le quantifier. On connaît par la demande de brevet EP 1.828.737, une méthode pour détecter l'apparition du pré allumage à forte charge, de type rumble. Cette méthode est basée sur la mesure d'un signal relatif au déroulement de la combustion, et une comparaison avec un signal-seuil. La présence d'une combustion anormale, de type "rumble", dans la chambre de combustion, est détectée lorsque l'amplitude du signal dépasse de façon significative celle du signal-seuil. Selon cette méthode, le signal-seuil correspond à l'amplitude du signal produit lors d'une combustion avec cliquetis ou lors d'une combustion normale (conventionnelle) Cependant, selon cette méthode, la détection ainsi réalisée ne permet pas d'agir au cours du cycle même de la détection. Les actions de corrections de ce type de pré allumage ne peuvent donc être réalisée qu'après la réalisation d'un tel phénomène, ce qui peut nuire sérieusement à l'intégrité du moteur. On connaît également la méthode décrite dans le brevet FR 2.897.900. Selon cette méthode, on peut agir plus rapidement après la détection du pré allumage : on est capable d'agir au cours du même cycle que le cycle de détection du phénomène. Pour ce faire, le signal-seuil est préalablement calculé, c'est-à-dire avant le fonctionnement du moteur, puis stocké dans des tables de données du calculateur, appelées cartographies. Cependant, l'utilisation de cartographies, ne permet pas de détecter à tout moment, c'est-à-dire en temps réel, le début d'un tel phénomène. De ce fait, il est toujours possible que la détection se fasse trop tardivement. De plus, aucune quantification de l'évolution du phénomène ne peut être réalisée. Ainsi, la nécessité ou non d'appliquer une phase de correction repose uniquement sur la comparaison de deux amplitudes à un instant donné. Or un tel phénomène peut très bien débuter, puis s'arrêter sans entraîner de dommage pour le moteur, et donc ne pas nécessiter de phase corrective. Ainsi, l'objet de l'invention concerne un procédé alternatif permettant de détecter en temps réel l'apparition d'un phénomène de pré allumage à forte charge (du type rumble), de le caractériser et de le quantifier, avec les dispositifs et systèmes couramment utilisés dans les moteurs, de façon à prendre des mesures permettant de l'éviter dans la suite du fonctionnement du moteur, au cours du même cycle que celui de la détection. La méthode s'appuie sur une modélisation d'indicateurs de l'état de la combustion et l'analyse de l'évolution des coefficients de la loi théorique servant à la modélisation. Le procédé selon l'invention L'objet de l'invention concerne un procédé de contrôle de la combustion d'un moteur à combustion interne à allumage commandé, dans lequel on détecte une combustion anormale dans une chambre de combustion (14) d'au moins un cylindre (12) du moteur, on enregistre au moins un signal représentatif d'un état de la combustion au moyen d'au moins un capteur placé dans le moteur, et on déduit du signal au moins un indicateur de la combustion à chaque cycle moteur. Le procédé comporte les étapes suivantes pour chaque cycle moteur : a- on modélise une distribution de N valeurs d'indicateurs de la combustion acquises sur N cycles précédant le cycle en cours, en déterminant des coefficients d'une loi de distribution théorique; b- on réitère l'étape précédente en retirant des valeurs extrêmes parmi les N valeurs ou, on réitère l'étape précédente pour différents points de fonctionnement du moteur, de façon à obtenir une évolution desdits coefficients; c- on détermine au moins un paramètre caractérisant une évolution d'au moins un des coefficients ; d- on détecte le début d'une combustion anormale en comparant ledit paramètre avec un seuil préalablement défini ; e- on contrôle le déroulement de la combustion anormale détectée dans la chambre de combustion. Selon l'invention, on peut définir un seuil directement à partir d'une tolérance sur une valeur brute de coefficients. On peut également, définir un seuil en calculant les dérivées premières et secondes de chacun des coefficients, et on fixe des bornes de variation maximales et relatives sur chacune des dérivées. On peut également déterminer une fréquence moyenne d'apparition du pré allumage, en analysant séparément chacun des coefficients, et en calculant une moyenne de nombres de pre allumage déterminés séparément grâce à chacun des coefficients. Le signal peut être une mesure continue de pression au sein du cylindre. L'indicateur peut être l'angle vilebrequin où 10 % de l'énergie introduite a été dégagée. Enfin, selon l'invention, on peut transformer les indicateurs en au moins un nouvel indicateur dont une distribution de N valeurs acquises sur N cycles précédant le cycle en 5 cours possède, pour des combustions normales, une dispersion plus faible que celles des indicateurs non transformés. Présentation des figures Les autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description donnée ci-après en se référant aux graphiques annexés où : la figure 1 montre la méthodologie générale de traitement des combustions anormales de type pré allumage. la figure 2 montre un moteur utilisant la méthode de détection selon l'invention. 10 - la figure 3 représente un histogramme de valeurs de CA10 et une loi normale permettant de s'ajuster au mieux à cet histogramme. la figure 4 représente un histogramme de valeurs de CA10 et une loi normale permettant de s'ajuster au mieux à cet histogramme, dans le cas d'apparition de pré allumage (PA). la figure 5 illustre la modélisation des valeurs de CA10 au moyen d'une loi non normale. 15 la figure 6 illustre l'évolution des coefficients d'une loi GEV lorsqu'on retire progressivement tous les échantillons à l'échantillon de données complet de départ. la figure 7 illustre l'évolution des dérivées premières des coefficients de la loi GEV en fonction de la taille de l'échantillon. la figure 8 illustre l'évolution des dérivées secondes des coefficients de la loi GEV en 20 fonction de la taille de l'échantillon. la figure 9 illustre l'évolution de la dérivée seconde du coefficient "Scale" de la loi GEV (moyenne, moyenne +1- l'écart type ajouté) la figure 10 illustre l'évolution des coefficients de la loi GG en fonction de la taille de l'échantillon. 25 - la figure 11 illustre l'évolution des dérivées premières des coefficients de la loi Gamma généralisée en fonction de la taille de l'échantillon. la figure 12 illustre l'évolution des dérivées secondes des coefficients de la loi Gamma généralisée en fonction de la taille de l'échantillon la figure 13 illustre l'évolution des coefficients de la loi GEV en fonction de la charge de 30 fonctionnement. la figure 14 est un exemple d'évolution du délai d'initiation de la combustion sur une acquisition de 1000 cycles la figure 15 illustre le CA10 vs dérivée angulaire de la pression à CA = CA10 la figure 16 illustre l'évolution de la grandeur "réduite" beta sur une acquisition de 1000 5 cycles la figure 17 illustre le CA10 vs dérivée angulaire de la pression à CA = CA10 la figure 18 illustre Y2 en fonction de X2, projection des cycles dans le repère obtenu à partir des directions principales identifiées dans le repère CA10, DP la figure 19 illustre Y2 (en haut) et X2 (en bas) pour chaque échantillon Description détaillée du procédé Sur la figure 2, un moteur à combustion interne 10 suralimenté à allumage commandé, en particulier de type essence, comprend au moins un cylindre 12 avec une chambre de combustion 14 à l'intérieur de laquelle se produit la combustion d'un mélange d'air 15 suralimenté et de carburant. Le cylindre comprend au moins un moyen d'alimentation en carburant sous 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 24 associée à une tubulure d'admission 26 se terminant par un plénum 26b 20 (non représenté sur la figure), au moins un moyen d'échappement des gaz brûlés 28 avec une soupape 30 et une tubulure 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 25 collecteur d'échappement 36 lui-même connecté à une ligne d'échappement 38. Un dispositif de suralimentation 40, par exemple un turbocompresseur ou un compresseur volumétrique. Dans le cas d'un turbocompresseur, celui-ci 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 30 admettre un air d'admission sous pression dans les chambres de combustion 14 par les tubulures d'admission 26. Le moteur comprend des moyens 46a de mesure de la pression cylindre, disposés au sein même du cylindre 12 du moteur. Ces moyens de mesure sont généralement constitués 10 par un capteur de pression qui permet de générer un signal représentatif de l'évolution de la pression dans un cylindre. Le moteur peut également comporter des moyens 46b de mesure de la pression d'admission, 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 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, comme la température de l'eau ou la température de l'huile, pour les traiter par calcul et ensuite commander les organes de ce moteur pour assurer son bon fonctionnement. Ainsi, dans le cas de l'exemple montré à la figure 2, 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 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. Au sein d'un tel moteur, le procédé selon l'invention permet de détecter l'apparition d'un phénomène de pré allumage à forte charge (du type rumble), de le caractériser et de le quantifier. Selon un exemple de réalisation, le procédé comporte les étapes suivantes : 1- on mesure la pression dans le cylindre, et on en déduit au moins un indicateur de la combustion à chaque cycle moteur ; Puis, pour chaque cycle : 2- on modélise des distributions de N valeurs d'indicateurs de la combustion acquises sur N cycles précédant le cycle en cours, en déterminant des coefficients d'une loi de distribution théorique pour différentes configurations : on réitère l'étape de modélisation en retirant des valeurs extrêmes parmi lesdites N valeurs ou, on réitère l'étape de modélisation pour différents points de fonctionnement dudit moteur ; 3- on détermine au moins un paramètre caractérisant une évolution desdits coefficients lorsqu'on s'affranchit progressivement de valeurs extrêmes parmi lesdites N valeurs, et on définit au moins un seuil pour ledit paramètre au moyen desdits coefficients ; 4- on détecte le début d'une combustion anormale en comparant le paramètre avec le seuil 5- on contrôle le déroulement de la combustion anormale détectée dans la chambre de combustion. 1. Détermination d'indicateur de la combustion à chaque cycle 10 On enregistre au moins un signal représentatif de l'état de la combustion, au moyen d'un capteur placé dans le moteur. Selon un mode de réalisation on choisit la pression cylindre. La mesure de la pression cylindre est réalisée à partir des moyens 46a de mesure de la pression cylindre. L'instrumentation des cylindres pour une mesure de pression est de plus en plus courante sur les véhicules. 15 Puis, au cours de chaque cycle, on calcule au moins un indicateur de la combustion, à partir de ce signal. Selon un mode de réalisation on choisit le CA10. Le CA10 correspond à l'angle vilebrequin où 10 % de l'énergie introduite a été dégagée. 2. Modélisation de la distribution de N valeurs d'indicateurs de la combustion. 20 La caractérisation d'une combustion normale, c'est-à-dire sans phénomène de pré allumage, est réalisée en ligne à partir de N valeurs de CA10 calculées sur les N-1 cycles précédent et le cycle en cours. On caractérise une combustion normale, c'est-à-dire sans phénomène de pré allumage, en déterminant un ou plusieurs paramètres caractéristiques de la distribution des 25 N valeurs de CA10 calculées sur les N-1 cycles précédents et le cycle en cours. Selon l'invention, ces paramètres sont déterminés au moyen d'une modélisation de la distribution des N valeurs de CA10. Cette modélisation consiste à déterminer les coefficients d'une loi théorique permettant de reproduire au mieux la distribution observée. Selon un mode de réalisation préféré. on réalise différentes modélisations en 30 s'affranchissant des valeurs dites « extrêmes », par un processus itératif. La dispersion cycle à cycle de la combustion est généralement représentée par une loi normale à laquelle on associe une moyenne, qui sert de référence pour caractériser le comportement du moteur, et un écart type qui sert classiquement à quantifier la stabilité de5 la grandeur observée. Or, il s'avère que l'écart type classique n'est pas toujours représentatif de la stabilité, et que la modélisation de la dispersion cyclique par une loi normale n'est pas optimale dans certaines conditions, comme l'illustre la figure 3. Sur cette figure est représentée l'histogramme de valeurs de CA10 (l'ordonnée représente la densité D) ainsi qu'une loi normale (Lnorm) permettant de s'ajuster au mieux à cet histogramme. Cette instabilité de combustion peut être directement ressentie par l'utilisateur au travers par exemple de ratés d'allumage (problématique de l'agrément de conduite mais aussi des émissions polluantes), et elle peut aussi être liée dans certains cas extrêmes à des combustions anormales pouvant entraîner la destruction du moteur (problématique du pré allumage à forte charge sur les moteurs à allumage commandé). La méthode de modélisation de la distribution, permet de fournir des informations plus pertinentes qu'un couple "moyenne / écart type" calculé au sens de la loi normale. Ceci est important pour détecter une combustion anormale comme du pré allumage. Un exemple de distribution expérimentale de CA10 et de modélisation par une loi normale (Lnorm) dans un cas d'apparition de pré allumage (PA) est donné sur la figure 4. On remarque que les cycles avec pré allumages étirent la distribution vers les faibles valeurs de CA10. Par conséquent, la modélisation effectuée (que ce soit sur la base d'une loi normale ou autre) est inévitablement faussée par cet étirement car celui-ci ne représente en rien une combustion normale. La modélisation des combustions normales seules (c'est à dire des combustions qui ont été initiées par l'étincelle de la bougie) ne peut pas être réalisée à partir de l'échantillon complet de N cycles. La méthode selon l'invention permet ensuite, à partir de ces modélisations de la loi de distribution des indicateurs de combustion, de quantifier la proportion de pré allumage en suivant les évolutions des coefficients des lois de distributions théoriques utilisées. La méthode comporte alors les étapes suivantes : 1. On retire Xvaleurs de CA10, parmi les N valeurs de la distribution à modéliser. De façon préférentielle on retire les valeurs les plus faibles. 2. On modélise la distribution des valeurs de CA10 restantes, en choisissant une loi à priori (loi gamma, « Generalized Extreme Values », Weibull, Nakagami, log-logistic, ...), comme illustré sur la figure 5. 3. on retire encore des valeurs de CA10, parmi les N -X valeurs de CA10. De façon préférentielle on retire les valeurs les plus faibles restantes. Puis on réitère le processus. On obtient ainsi une évolution des coefficients d'une loi théorique modélisant la distribution des valeurs de CA10. Le fait de retirer les X cycles pour lesquels les valeurs de CA10 sont les plus faibles, revient à supprimer de l'échantillon des valeurs qui correspondent potentiellement à des pré allumages. Après quelques itérations, l'algorithme permet d'éliminer tous les pré allumages et donc d'aboutir à un échantillon composé uniquement de combustions normalement initiées à la bougie. Le seuil de détection entre les pré allumages et les combustions normales est donc fixé automatiquement par cet algorithme. Par exemple, si on considère une loi de distribution du type GEV, la fonction de densité a la forme générale suivante : (-'/e)-I Î (x, p, ) = 1 6 1 + I x û exp Les coefficients réels considérés sont donc : ,u E [ 00, mi o.E[O,ao] E [û00,00]Ji û 1+4x6~ -le J appelé coefficient "location" appelé coefficient "scale" appelé coefficient "shape" Les évolutions de ces coefficients au cours du processus itératif montrent des tendances très particulières. Plus précisément, lorsque l'échantillon contient du pré allumage ou lorsqu'il n'en contient pas l'évolution est différente. Il apparaît donc judicieux d'utiliser ces évolutions pour quantifier la proportion de pré allumage. La figure 6 illustre l'évolution des coefficients d'une loi GEV, lorsqu'on retire progressivement tous les échantillons à l'échantillon de données complet de départ. Dans le cas des modélisations itératives par une loi GEV présentées sur la figure 6, on constate que les coefficients 6 (« Scale ») et ,u (« Location ») présentent des évolutions très caractéristiques. II est donc possible de quantifier la fréquence du pré allumage à partir de ces évolutions, en utilisant des valeurs singulières de ces courbes d'évolution, comme des valeurs maximales ou des points d'inflexion. Selon un autre mode de réalisation, les coefficients des lois de distributions théoriques ne sont pas utilisés dans le cadre du processus itératif. II est en effet possible de détecter un 30 début d'apparition de pré allumage en fonction de la charge moteur en appliquant directement les modélisations statistiques à différents points de fonctionnement, et en comparant directement les valeurs des coefficients de la loi utilisée (alors que sur les exemples précédents ils s'agissaient de charge stabilisée). La figure 13 illustre cette comparaison sur une iso régime (i.e. sur différents points enregistrés au même régime de fonctionnement mais à des charges différentes). Elle permet de comparer les valeurs des coefficients de la loi GEV lorsque l'on cherche à modéliser la distribution des CA10 de chacun des points de fonctionnement avec cette loi particulière (les paramètres "k", "sigma" et "mu" qui apparaissent sur cette figure sont appelés "Shape", "Scale" et "Location" sur les figures précédentes). L'apparition du pré allumage se manifeste encore une fois par une évolution particulière de ces coefficients, en particulier sous la forme d'une augmentation importante du coefficient "mu". La méthodologie de détection avec une loi GEV consiste donc, par exemple, à détecter cette franche augmentation en considérant par exemple la valeur moyenne de ce coefficient avant l'apparition du pré allumage à plus faible charge comme étant la référence, et qu'un dépassement de 30 % de cette valeur correspond à l'amorce du pré allumage. 3. Définition de seuils et détection de combustion avec pré allumage Selon un mode de réalisation, on définit un seuil directement à partir d'une tolérance sur la valeur brute des coefficients (p, q. Selon un autre mode de réalisation, on calcule les dérivées premières et secondes de chacun des coefficients (p, a-, 4), et on fixe des bornes de variation maximales et relatives sur chacune de ces dérivées. Les dérivées des trois coefficients présentés sur la figure 6 sont illustrées sur les figures 7 et 8. On peut par exemple calculer sur les coefficients o (« Scale ») et p (« Location ») la valeur moyenne des dérivées premières et secondes sur une plage "bien choisie" et fixer le seuil de détection du pré allumage à une variation correspondant à 2 fois l'écart type autour de la valeur moyenne (illustré sur la dérivée seconde du coefficient 6(« Scale ») sur la figure 9). Bien choisir la zone de calcul consiste à choisir une zone dans laquelle on sait a priori qu'il n'y a pas de pré allumage. Sur la figure 9 par exemple, on a choisi de caractériser la valeur moyenne et l'écart type du coefficient 6 (« Scale ») de la loi GEV sur une plage de 200 cycles uniquement, en évitant ainsi les perturbations : liées aux pré allumages potentiels (on prend une marge confortable en considérant qu'il ne faut pas utiliser les 300 cycles ayant les CA10 les plus faibles étant donné qu'il est extrêmement rare d'avoir 300 cycles en pré allumage sur une acquisition de seulement 1000 cycles) ; - et celles liées aux défauts de convergence des modélisations lors que l'échantillon modélisé ne contient pas suffisamment de données (taille d'échantillon inférieur à 200 cycles par exemple). On peut également déterminer une fréquence moyenne d'apparition du pré allumage en analysant séparément chacun des coefficients, et en faisant ensuite une moyenne des nombres de pre allumage déterminés séparément grâce à chacun desdits coefficients, pour avoir un indicateur moyen de la quantité de pré allumage. La figure 10 illustre les évolutions des coefficients d'une loi Gamma (LG) généralisée sur le même processus itératif. Encore une fois, les coefficients montrent une évolution particulière lorsque les cycles supposés être des pré allumages sont retirés à l'échantillon complet de départ (les coefficients "b" et "c" se stabilisent très clairement lorsqu'une centaine de cycles a été retirée de l'échantillon de départ, le coefficient "a" quant à lui montre lui aussi une évolution particulière puisqu'il atteint un maximum lorsque tous les pré allumages ont été retirés). Les figures 11 et 12 illustrent quant à elles les évolutions des dérivées premières et secondes de ces coefficients pour la loi Gamma généralisée. 4- Contrôle de la combustion anormale Au moyen de cette comparaison, le calculateur moteur peut détecter le début d'une combustion anormale de type "rumble" ou du type "pré allumage" dans la chambre de combustion. En cas de combustion anormale, ce calculateur lance ensuite les actions nécessaires au contrôle de cette combustion afin d'éviter la poursuite d'une telle combustion. Par contrôle de la combustion anormale, il est entendu non seulement la possibilité de maîtriser le déroulement de cette combustion pour éviter les augmentations brutales de pressions destructrices mais aussi d'arrêter complètement une telle combustion, telle que par étouffement. A titre préférentiel, ce contrôle de la combustion est réalisé par une réinjection de carburant à un angle de vilebrequin déterminé par les injecteurs 18. Plus précisément, le calculateur commande les vannes 20 de façon à ce que l'injecteur du cylindre concerné permette d'introduire dans la chambre de combustion une quantité de carburant sous forme liquide. La quantité de carburant réinjectée dépend de la constitution du moteur et peut aller de 10% à 200% de la quantité de carburant initialement introduite dans cette chambre de combustion. De ce fait, le carburant réinjecté sert à contrarier la flamme qui commence à se déployer lors de la combustion anormale. Cette réinjection permet soit de souffler cette flamme, soit d'étouffer cette flamme par augmentation de la richesse du mélange carburé. De plus, le carburant injecté sous forme liquide utilise la chaleur présente autour de cette flamme pour se vaporiser et les conditions de température autour de la flamme vont baisser en retardant la combustion du mélange carburé et surtout son auto-inflammation. Après cette injection de carburant, la pression dans le cylindre augmente mais moins 10 brutalement. Cette pression décroît ensuite pour atteindre un niveau compatible avec le niveau de pression d'une combustion conventionnelle. Par ce mécanisme, tout développement d'une combustion anormale avec une grande vitesse de combustion et des pressions élevées est prohibé. Bien entendu, la mise en oeuvre des moyens pour contrôler la combustion anormale se fait à chaque cycle durant 15 lequel une telle combustion est détectée par le calculateur. Les actions du procédé telles que décrites ci-dessus peuvent être combinées à d'autres actions plus lentes, telles que la fermeture du papillon, pour empêcher que les conditions de pression de la chambre de combustion soient favorables à une combustion anormale dans les cycles qui suivent. 5- Autres modes de réalisation L'invention permet d'utiliser d'autres mesures que la pression cylindre, telle que le couple instantané, le régime instantané, le niveau de vibrations (capteurs accélérométriques), signal d'ionisation,.... 25 L'invention peut également utiliser d'autres indicateurs de combustion : à partir de la pression cylindre : PMI, pression cylindre maximale, angle vilebrequin à la pression maximale, CAxx, maximum de dégagement d'énergie,... à partir du couple instantané : maximum de couple, dérivée maximale de couple, ... à partir du régime instantané : maximum de vitesse, accélération maximale, ... 30 le volume de la chambre de combustion, ou le gradient de volume à certains moments (au CA10 par exemple). 20 L'invention n'est pas limitée à l'utilisation d'un seul indicateur de la combustion. En effet, il peut être avantageux de combiner simultanément plusieurs de ces indicateurs. Dans les procédés actuels, de tels indicateurs sont toujours analysés séparément. Pour analyser simultanément des valeurs de plusieurs indicateurs de l'état d'avancement de la combustion (CA10, PMI, ...), on combine plusieurs de ces indicateurs, pour obtenir d'autres indicateurs permettant de distinguer plus rapidement et plus efficacement des combustions normales et des combustions anormales comme le pré allumage. Indicateurs "directs" : x1, x2, ...xn Indicateurs "réduits" : X1, X2, ...X, En choisissant judicieusement les indicateurs à combiner (CA1o, pression cylindre maximale, dérivée de la pression cylindre maximale,...), il est possible d'obtenir des indicateurs, dits "indicateurs réduits" et notés X,, X2, ..., Xn , qui présentent, dans le cas de combustions normales, des dispersions extrêmement plus faibles que les indicateurs directs x,, x2, ..., xn. Il est ainsi possible de « gommer » la dispersion naturelle associée à la combustion, tout en conservant la même sensibilité au pré allumage. On opère ainsi une réduction de dimension (un indicateur défini une dimension) pour exploiter au mieux les relations entre les différents indicateurs. Pour réaliser cette réduction de dimension, on cherche les directions de l'espace (x1, x2, ..., xn) qui représentent le mieux les corrélations entre les différentes grandeurs. L'hypothèse qui motive cette approche est que les combustions anormales du type pré allumage ont un mode de combustion qui diffère sous certains aspects de celui suivi par une combustion normale. En d'autres termes, il y a une ou plusieurs caractéristiques de ces pré allumages qui présentent une discontinuité, permettant de les distinguer des autres combustions normales. En mettant en évidence les corrélations entre différents indicateurs on met plus facilement en évidence ces changements de modes de combustion. Selon un premier mode de réalisation, on recherche les corrélations entre différents indicateurs au moyen d'une régression entre les indicateurs. La figure 14 illustre l'évolution du CA10 lors d'une acquisition de 1000 cycles moteur en 30 régime stabilisé. La dispersion associée aux combustions normales est de l'ordre de 4 °V. Nous pouvons aussi constater que les pré allumages s'éloignent de la moyenne des CA10 au maximum de l'ordre de 19 degrés vilebrequin (19 °V). La figure 15 présente quant à elle ce même échantillon de points dans un plan CA10 / DP, où DP est la dérivée angulaire de la pression lorsque l'angle vilebrequin est égale à l'angle CA10. Cette représentation permet de constater qu'une relation linéaire lie ces deux variables dans le cas des combustions normales (CN), alors que les combustions anormales (CAN) s'éloignent très clairement de cette tendance. La droite en pointillé sur la figure 15 illustre la droite de régression linéaire entre les deux indicateurs de combustion CA10 et DP pour des combustions normales (CN). L'idée est donc de combiner ces deux indicateurs pour en extraire un troisième potentiellement plus intéressant pour discriminer plus facilement les combustions anormales des normales. On peut par exemple utiliser comme troisième indicateur, indicateur transformé ou combiné, la variable fi définie par l'erreur entre la régression et la CA10 : f=CA10-(a.DP+y) où (a.DP + y) est la droite de régression linéaire entre CA10 et DP pour des combustions normales (CN). La figure 16 illustre l'évolution de la variable f pour les 1000 cycles moteur (NCy). On constate que la dispersion associée aux combustions normales est inférieure à 0.5 ° vilebrequin, alors que les pré allumages restent potentiellement toujours aussi éloignés de la moyenne des valeurs de /3. Ce changement de variable permet donc de resserrer énormément la dispersion normale et d'identifier ainsi plus précisément les combustions anormales. Selon un second mode de réalisation, on recherche les corrélations entre différents indicateurs au moyen d'une analyse en composantes principales, lesdits nouveaux indicateurs correspondants aux composantes principales. On transforme ainsi des indicateurs corrélés en nouveaux indicateurs, indépendants les uns des autres. Ces indicateurs transformés permettent d'expliquer au mieux la variabilité des indicateurs non transformés. L'Analyse en Composantes Principales est une Analyse Factorielle de la famille de l'Analyse des données multivariées, qui consiste à transformer des variables liées entre elles (variables corrélées) en nouvelles variables indépendantes les unes des autres. Cette technique permet de réduire l'information en un nombre de composantes plus limité que le nombre initial de variables. Ces composantes permettant d'expliquer au mieux la variabilité (la variance) des données, c'est-à-dire de "capter" un maximum d'information. Plusieurs degrés de liberté peuvent être utilisés pour optimiser le procédé : le nombre d'indicateurs utilisés ; le choix des indicateurs à utiliser ; la méthode de réduction utilisée pour identifier les corrélations entre les indicateurs retenus (ACP par exemple). Les figures 17 à 19 permettent d'illustrer un exemple simple dans le cas du traitement de données acquises au banc moteur. On effectue un changement de repère, en définissant par ACP deux nouveaux axes (composantes principales) X2 et Y2 (X2 étant "porté" par la zone linéaire et dense des combustions normales, et X2 et Y2 sont orthogonaux), représentés schématiquement sur la figure 17. La représentation du même échantillon de point dans le repère X2, Y2 est fournie sur la figure 18. On analyse ensuite les composantes principales (axes) pour différencier les pré allumages des combustions anormales. Les évolutions de ces composantes sont visibles sur la figure 19 (Y2 en haut, X2 en bas). On peut y constater, comme sur la figure 16, que les pré allumages gardent la même "dynamique" par rapport aux combustions normales qui sont elles beaucoup plus resserrées. Il est en revanche beaucoup plus simple et direct de délimiter des combustions normales et du pré allumage, étant donné que les combustions normales sont toutes très resserrées. Selon ce mode de réalisation, après avoir déterminé ces nouveaux indicateurs, on analyse l'évolution d'au moins un coefficient des lois théoriques modélisant la dispersion de ces nouveaux indicateurs. REVENDICATIONS1. Procédé de contrôle de la combustion d'un moteur à combustion interne à allumage commandé, dans lequel on détecte une combustion anormale dans une chambre de combustion (14) d'au moins un cylindre (12) dudit moteur, on enregistre au moins un signal représentatif d'un état de la combustion au moyen d'au moins un capteur placé dans le moteur, et on déduit dudit signal au moins un indicateur de la combustion à chaque cycle moteur, caractérisé en ce que, pour chaque cycle moteur : a- on modélise une distribution de N valeurs d'indicateurs de la combustion acquises sur N cycles précédant le cycle en cours, en déterminant des coefficients d'une loi de distribution théorique; b- on réitère l'étape précédente en retirant des valeurs extrêmes parmi lesdites N valeurs ou, on réitère l'étape précédente pour différents points de fonctionnement dudit moteur, de façon à obtenir une évolution desdits coefficients; c- on détermine au moins un paramètre caractérisant une évolution d'au moins un desdits coefficients ; d- on détecte le début d'une combustion anormale en comparant ledit paramètre avec un seuil préalablement défini ; e- on contrôle le déroulement de la combustion anormale détectée dans la chambre de combustion. 2. Procédé selon la revendication 1, dans lequel on définit ledit seuil directement à partir d'une tolérance sur une valeur brute de coefficients. 3. Procédé selon la revendication 1, dans lequel on définit ledit seuil en calculant les dérivées premières et secondes de chacun desdits coefficients, et on fixe des bornes de variation maximales et relatives sur chacune desdites dérivées. 4. Procédé selon la revendication 1, dans lequel on détermine une fréquence moyenne d'apparition du pré allumage en analysant séparément chacun desdits coefficients, et en calculant une moyenne de nombres de pre allumage déterminés séparément grâce à chacun desdits coefficients. 5. Procédé selon l'une des revendications précédentes, dans lequel ledit signal est une mesure continue de pression au sein dudit cylindre. 6. Procédé selon l'une des revendications précédentes, dans lequel ledit indicateur est l'angle vilebrequin où 10 % de l'énergie introduite a été dégagée. 7. Procédé selon l'une des revendications précédentes, dans lequel, avant l'étape a, on transforme lesdits indicateurs en au moins un nouvel indicateur dont une distribution de N valeurs acquises sur N cycles précédant le cycle en cours possède, pour des combustions normales, une dispersion plus faible que celles des indicateurs non transformés.