FR 3074542 A1 20190607 FR 1761708 A 20171206 La presente invention concerne le domaine de la commande de la combustion d'un moteur a combustion interne a allumage commande, dans le but de recluire le phenomene de cliquetis. Ce type de moteur comprend au moins un cylindre comportant une chambre de combustion delimitee par la paroi laterale interne du cylindre, par le haut du piston qui coulisse dans ce cylindre et par la culasse. Generalement, un mélange carbure est renferme dans cette chambre de combustion et subit une etape de compression puis une etape de combustion sous l'effet d'un allumage commande, par une Bougie, ces etapes etant regroupees sous le vocable de « phase de combustion » dans la suite de la description. Dans un moteur a explosion, la combustion du mélange air/essence commence normalement apres l'etincelle. 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'elevation de pression et de temperature est parfois suffisamment importante pour que le mélange non brOle coince contre les parois atteigne son point d'auto-allumage et s'auto-enflamme en un ou plusieurs endroits. Ce phenomene est appele « cliquetis ». Ainsi le cliquetis est avant tout un phenomene de combustion anormale dans les moteurs a allumage commande, perceptible exterieurement par un bruit metallique venant du moteur qui resulte de ('apparition d'ondes de pression dans la chambre de combustion. Ces explosions parasites produisent des vibrations dans le domaine acoustique et au- dela (de l'ordre de 5 a 50 KHz). Elles sont fres vives et peuvent rapidement conduire a un echauffement localise. L'accumulation de cliquetis endommage le metal du piston et/ou des parois du cylindre et les segments. Au bout de quelque temps (selon l'intensite) cela conduit a la destruction du piston, des segments ou des parois du cylindre. L'estimation du cliquetis permet un controle de la combustion limitant l'effet du cliquetis et n'endommageant pas le cylindre. Etat de la technique On connait des methodes pour estimer le cliquetis d'un moteur a combustion interne. Elles sont classiquement basees sur l'enregistrement de signaux issus d'un capteur de pression cylindre ou d'un accelerometre. Une premiere methode consiste a appliquer hors ligne, c'est a dire en traitement ulterieur, une technique de transformations de Fourier a ces signaux pour estimer le cliquetis. Une telle methode est decrite dans les documents suivants : 2 - "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) 5 Une seconde methode consiste a appliquer aux signaux une technique de filtrage et de detections de maximum pour estimer le cliquetis hors ligne. Une telle methode est decrite dans le document suivant : - "Modeling of In-cylinder Pressure Oscillations under Knocking Conditions: Introduction to Pressure Envelope Curve", G. Brecq et 0. Le Corre, SAE 2005 10 En outre, on connait une autre methode, egalement basee sur l'enregistrement de signaux issus d'un capteur de pression cylindre. Cette fois, une technique d'ondelettes est appliquee aux signaux pour estimer le cliquetis hors ligne. Une telle methode est decrite dans le document suivant : - "Knocking detection using wavelet instantaneous correlation method", Z. Zhang et , E. 15 Tomita, Journal of SAE Review volume 23 (2002) Cependant, ces methodes ne permettent pas une estimation de l'intensite du cliquetis en temps reel (le calcul doit se faire typiquement en moins de 50 ps pour etre embarque sur un calculateur). 20 On connait par la demande FR 2 949 511 une methode pour estimer en temps reel l'intensite de cliquetis, a partir d'un signal vibratoire. Selon cette methode, on determine en temps reel les coefficients de la decomposition de Fourier de ce signal, on en deduit l'energie contenue dans le signal en sommant les carres de ces coefficients de la decomposition de 25 Fourier. Enfin, on estime l'intensite de cliquetis en determinant un parametre correle a l'intensite de cliquetis egal a la racine carree du maximum de l'energie. II a ete egalement envisage d'utiliser d'autres capteurs comme par exemple une sonde d'ionisation. 30 Concernant la mesure du cliquetis, on obtient de maniere generale un nombre scalaire qui reflete l'intensite du phenomene qui a eu lieu au cours d'un cycle. L'indicateur, quelconque, est appele ici KI (de l'anglais « knock indicator », ce que Ion pourra traduire par indicateur cliquetis). La variabilite des mesures successives de KI, independamment de tout 3 bruit de mesure, est importante au regard des variations de leur moyenne temporelle sur des echelles de temps caracteristiques de revolution d'etat moteur. La distribution des mesures KI, et non des mesures considerees individuellement, contient ('information renseignant sur la propension du moteur a voir apparaitre des cycles cliquetants. II pourrait etre envisage 5 d'utiliser la moyenne des valeurs de KI. Toutefois, la caracterisation de cette distribution par sa seule moyenne est insuffisante pour prevenir ('occurrence de forts KI, susceptibles d'endommager le moteur. En effet, ('experience a montre que la dispersion des mesures n'est pas liee a leur moyenne. Un controleur pilotant correctement le KI moyen sur une cible ne controle pas la probabilite d'occurrence des fortes valeurs de la distribution, ce qui ne 10 permet pas un controle sir du cliquetis. Pour pallier ces inconvenients, la presente invention concerne un procede de commande d'un moteur a combustion interne a allumage commande, dans lequel on commande l'avance a l'allumage au moyen d'une estimation de la distribution des mesures 15 de cliquetis KI. Cette estimation permet de determiner, pour ces mesures, un intervalle de confiance d'un quantile predetermine de la distribution des KI. Ainsi, it est possible de commander, en temps reel, le moteur a combustion interne de maniere a assurer une proportion donnee de cycles consideres a risque, tout en ayant un nombre de mesures limite. De plus, ('estimation, de preference bayesienne, permet la realisation de calculs 20 legers (simples operations scalaires), necessitant peu de temps de calcul, peu de memoire et des ressources materielles moindres. En outre, le procede selon ('invention maximise refficacite energetique moyenne du moteur grace au pilotage fin du niveau de cliquetis. Le procede selon ('invention 25 L'invention concerne un procede de commande d'un moteur a combustion interne a allumage commande, ledit moteur a combustion interne comprenant un capteur de cliquetis agence sur ledit moteur a combustion interne. Ledit procede comporte les etapes suivantes : a) on realise N mesures de cliquetis au moyen dudit capteur de cliquetis ; b) on construit un estimateur d'une distribution desdites N mesures, ledit estimateur 30 determinant un intervalle de confiance d'un quantile predetermine au moyen desdites N mesures ; c) on commande l'allumage au sein dudit moteur a combustion interne en fonction d'une comparaison avec ledit intervalle de confiance : i) si ledit seuil de cliquetis est inferieur a la valeur minimale dudit intervalle de confiance dudit quantile predetermine, alors on diminue l'avance a l'allumage ; 4 ii) si ledit seuil de cliquetis est superieur a la valeur maximale dudit intervalle de confiance dudit quantile predetermine, alors on augmente I'avance a I'allumage ; et iii) si ledit seuil de cliquetis est compris dans ledit intervalle de confiance, alors on ne 5 modifie pas I'avance a I'allumage. Selon un mode de realisation, ledit estimateur est un estimateur bayesien. Conformement a une mise en oeuvre, ladite distribution est une distribution pouvant etre rapportee a une distribution normale par une methode mathematique, de preference 10 ladite distribution est une distribution log-normale. Selon un aspect, ledit estimateur determine en outre l'esperance q dudit quantile predetermine. Avantageusement, on diminue ou on augmente I'avance a I'allumage de maniere proportionnelle a la difference entre ladite esperance dudit quantile predetermine et ledit 15 seuil de cliquetis. Selon une mise en oeuvre, N est compris entre 5 et 20, de preference entre 7 et 15. Conformement a un mode de realisation, ledit intervalle de confiance contient la valeur du quantile reel dans exactement 98% des cas. De maniere avantageuse, ledit capteur de cliquetis est un accelerometre agence sur la 20 culasse dudit moteur a combustion interne. De preference, ladite mesure de cliquetis est le maximum d'amplitude du signal filtre dudit accelerometre. Selon un mode de realisation de l'invention, ledit seuil de cliquetis est obtenu par essais au banc moteur durant la phase de calibration du moteur. 25 Conformement a un aspect, ledit quantile predetermine de ladite distribution est superieur ou egal a 80%. De plus, l'invention concerne un systeme de commande d'un moteur a combustion interne a allumage commande comprenant des moyens pour mettre en oeuvre le procede 30 selon l'une des caracteristiques precedentes. En outre, l'invention concerne un produit programme d'ordinateur telechargeable depuis un reseau de communication et/ou enregistre sur un support lisible par ordinateur et/ou executable par un processeur, comprenant des instructions de code de programme 35 pour la mise en oeuvre du procede selon l'une des caracteristiques precedentes, lorsque ledit programme est execute sur un ordinateur. 5 Presentation succincte des figures D'autres caracteristiques et avantages du procede selon l'invention, apparaitront a la lecture de la description ci-apres d'exemples non limitatifs de realisations, en se referent aux 5 figures annexees et decrites ci-apres. La figure 1 illustre les etapes du procede selon un mode de realisation de l'invention. La figure 2 illustre une distribution des mesures de cliquetis. Les figures 3 a 5 illustrent l'etape de comparaison selon les trois cas de figure : le seuil est compris dans l'intervalle de confiance, le seuil est superieur a l'intervalle de confiance, et 10 le seuil est inferieur a l'intervalle de confiance. Les figures 6a et 6b representent respectivement l'avance a l'allumage et la mesure de cliquetis entouree par ('estimation de son intervalle de confiance obtenu en regime stabilise avec le procede de commande selon un mode de realisation de l'invention. Les figures 6a et 6b representent respectivement l'avance a l'allumage et la mesure de 15 cliquetis obtenu en regime stabilise avec un procede de commande selon l'art anterieur. Description detainee de l'invention Le procede selon l'invention concerne un procede de commande de la combustion d'un moteur a combustion interne prenant en compte les phenomenes de cliquetis. Le 20 procede comporte la comparaison de valeurs obtenues par un estimateur (esperance d'un quantile predetermine et intervalle de confiance) a un seuil de cliquetis, les valeurs stochastiques etant determinees a partir d'un traitement d'un signal issu d'un capteur de cliquetis, ce signal ayant un sens physique pour qualifier directement la combustion et l'eventuelle apparition de cliquetis. Le procede est independent de la nature de la 25 combustion, de la technologie du capteur utilise et de l'indicateur KI. Ce type de moteur comprend au moins un cylindre comportant une chambre de combustion delimitee par la paroi laterale interne du cylindre, par le haut du piston qui coulisse dans ce cylindre et par la culasse. Generalement, un mélange carbure est renferme dans cette chambre de combustion et subit une etape de compression puis une etape de 30 combustion sous l'effet d'un allumage commande, au moins par une Bougie. Le capteur de cliquetis utilise pour le procede de commande selon l'invention, peut-titre de tout type, par exemple un accelerometre ou un capteur de pression. 6 Selon un mode de realisation de ('invention, le capteur de cliquetis est un accelerometre agence sur la culasse du moteur a combustion interne. Le procede selon ('invention comporte les etapes suivantes : 5 1) Mesures de cliquetis 2) Construction d'un estimateur 3) Comparaison avec un seuil 4) Commande de l'allumage au sein du moteur a combustion interne 10 La figure 1 illustre, schernatiquement et de maniere non !imitative, les etapes du procede selon un mode de realisation de l'invention. Le procede comporte une etape de realisation de N mesures (N MES) au moyen du capteur de cliquetis, une etape de construction d'un estimateur (EST) qui determine une esperance q (de maniere optionnelle) et un intervalle de confiance INT d'un quantile predetermine, une etape de comparaison 15 (COMP) avec une valeur seuil 0 et une etape de commande (COM) en temps reel. De maniere optionnelle, et tel qu'illustre sur la figure 1, le procede peut comporter egalement une etape de determination de la difference (DIF) entre ('esperance q du quantile predetermine et la valeur seuil 0, cette difference peut etre prise en compte pour l'etape de commande (COM). 20 Les differentes etapes du procede selon ('invention vont etre decrites ci-dessous en reference a la figure 1. 1) Mesures de cliquetis (N MES) 25 Lors de cette etape, on realise N mesures de cliquetis au moyen du capteur de cliquetis. Chaque mesure est une mesure realisee pour un cycle de combustion. Ces mesures peuvent etre enregistrees dans une memoire interne. Par consequent, pour chaque cycle, on connait les mesures pour les N cycles precedents le cycle en cours. Pour le mode de realisation de l'invention pour lequel le capteur de cliquetis est un 30 accelerometre, la mesure peut consister en le maximum d'amplitude du signal filtre de ('accelerometre. Cette valeur reflete l'intensite du cliquetis (KI). Selon un mode de realisation de l'invention, le nombre N est superieur ou egal a 5, de preference compris entre 5 et 20, et de maniere preferee entre 7 et 15. En effet, au moins 5 35 valeurs permettent de construire un estimateur representatif. De plus, la limite de 20 7 mesures permet de former un estimateur reactif et permet de limiter la memoire vive utilisee pour le procede selon ('invention. Cependant, pour une application ne necessitant pas un estimateur reactif, le nombre N peut etre pris aussi grand que possible. 5 2) Construction d'un estimateur (EST) II s'agit de construire un estimateur qui determine un intervalle de confiance INT autour de l'esperance d'un quantile predetermine au moyen des N mesures realisees a l'etape precedente. L'estimation de ('intervalle de confiance permet d'obtenir une commande robuste. En outre, ('estimateur peut egalement determiner l'esperance du quantile 10 predeterminee. Pour cette etape, on considere que les mesures suivent un loi de distribution stochastique. De preference, la loi stochastique est une distribution pouvant etre rapportee a une distribution normale par une methode mathematique, de preference ladite distribution est une distribution log-normale. 15 Pour cette construction, on considere que les mesures de cliquetis suivent une distribution log-normale. II s'agit d'une hypothese courante dans la litterature. Le document suivant decrit une telle consideration : « Spelina, J. M.; Peyton Jones, J. C. & Frey, J. Characterization of knock intensity distributions: Part 2: parametric models Proceedings of 20 the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, Sage Publications Sage UK: London, England, 2013, 227, 1650-1660 ». En effet, cette distribution est particulierement adaptee, car le cliquetis est un phenomene exponentiel, pour lequel la distribution correspond sensiblement a une distribution normale en echelle logarithmique. Une telle distribution est determinee par les 25 deux parametres suivants : la moyenne et l'ecart type. Avec ces valeurs, it est possible de calculer tout quantile de cette distribution. Selon un aspect de l'invention, le quantile predetermine correspond a au moins 80 %. Ce quantile eleve est un bon indicateur de l'etat cliquetis, les cycles dangereux 30 correspondant aux valeurs elevees de KI. De preference, le quantile predetermine peut valoir 80%, 95%, 99% ou 99,7%. Une telle estimation peut etre realisee rapidement avec un nombre limite de mesures, ce qui permet une commande en temps reel, avec une demande en memoire limitee. De 35 plus, cet estimateur apporte la stabilite au procede selon l'invention, quel que soit le nombre de mesures. 8 Selon un mode de realisation de l'invention, l'estimateur est un estimateur bayesien. Cet estimateur permet de determiner l'intervalle de confiance de ('estimation du quantile de maniere reactive avec un nombre limite de mesures. 5 Conformement a une mise en oeuvre de l'invention, l'intervalle de confiance contient la valeur du quantile reel dans exactement 98% des cas. En d'autres termes, la valeur du quantile reel est inferieure a l'intervalle de confiance dans 1% des cas et superieure a ('intervalle de confiance dans 1% des cas. 10 Un exemple d'estimation bayesienne d'un intervalle de confiance est donne, de maniere non limitative, en annexe de la presente demande (a la fin de la description). 3) Comparaison avec un seuil (COMP) Lors de cette etape, on compare l'esperance et l'intervalle de confiance determines a 15 l'etape precedente avec un seuil de cliquetis. En d'autres termes, it s'agit de verifier si le seuil de cliquetis est dans l'intervalle de confiance. Le seuil de cliquetis peut etre obtenu par essais au banc moteur durant la phase de calibration du moteur. 20 Les figures 2 a 5, illustrent schematiquement et de maniere non limitative, cette etape. La figure 2 illustre une densite de probabilite de l'intensite du signal du capteur de cliquetis KI. II s'agit d'une distribution du type log-normale. Sur cette figure est egalement illustre le seuil de cliquetis 0. La zone hachuree represente la zone pour laquelle le cliquetis est present. 25 La figure 3 illustre la densite de probabilite de la figure 2, sur laquelle on represente l'esperance q, et l'intervalle de confiance, delimite par sa valeur minimale qn-', et sa valeur maximale gmax- L'esperance et l'intervalle de confiance sont obtenus pour le quantile predetermine (par exemple 95%) au moyen des N mesures. Pour cet exemple, le seuil de cliquetis 0 est dans l'intervalle de confiance determine par l'estimateur bayesien. Cette 30 configuration represente 1% des cas de figure. La figure 4 illustre la densite de probabilite de la figure 2, sur laquelle on represente l'esperance q, et l'intervalle de confiance, delimite par sa valeur minimale chi, et sa valeur maximale q'x. L'esperance et l'intervalle de confiance sont obtenus pour le quantile predetermine (par exemple 95%) au moyen des N mesures. Pour cet exemple, le seuil de 35 cliquetis Best superieur a l'intervalle de confiance determine par l'estimateur bayesien. Cette configuration represente 98% des cas de figure. 9 La figure 5 illustre la densite de probabilite de la figure 2, sur laquelle on represente l'esperance q, et l'intervalle de confiance, delimite par sa valeur minimale qm,n et sa valeur maximale q'x. L'esperance et l'intervalle de confiance sont obtenus pour le quantile predetermine (par exemple 95%) au moyen des N mesures. Pour cet exemple, le seuil de 5 cliquetis 0 est inferieur a l'intervalle de confiance determine par l'estimateur bayesien. Cette configuration represente 1% des cas de figure. 3') Difference (DIF) II s'agit d'une etape optionnelle du procede selon ('invention. 10 Lors de cette etape, on determine la difference entre l'esperance q du quantile predetermine determine a l'etape precedente et le seuil de cliquetis 0. La determination de cette difference permet sa prise en compte dans la commande, pour realiser une commande adaptee a l'intensite du cliquetis. 15 4) Commande de l'allumage au sein du moteur (COM) Lors de cette etape, on commande l'avance a l'allumage dans les cylindres du moteur a combustion interne, en fonction du cliquetis, en prenant en compte la comparaison realisee a l'etape 3). Selon l'invention, on commande l'avance a l'allumage du cycle suivant. En d'autres termes, avec les mesures des N cycles precedents le cycle n (c'est-a-dire les cycles 20 n-N, n-N+1, ..., n-2, n-1), on commande l'avance a l'allumage du cycle n. On rappelle que l'avance a l'allumage est un degre de liberte controlable qui determine, en angle vilebrequin, ('angle auquel la combustion est initie. II s'agit d'un actionneur a action instantanee (immediatement appliquee au prochain cycle), qui reduit la propension du moteur a cliqueter. 25 La commande est realisee de la maniere suivante i) si le seuil de cliquetis est inferieur a la valeur minimale de l'intervalle de confiance autour de l'esperance du quantile predetermine, alors on diminue l'avance a l'allumage (cas de la figure 5) ; 30 ii) si le seuil de cliquetis est superieur a la valeur maximale de l'intervalle de confiance autour de l'esperance du quantile predetermine, alors on augmente l'avance a l'allumage (cas de la figure 4) ; et iii) sinon (si ledit seuil de cliquetis est compris dans ledit intervalle de confiance (qmin, q,)) on ne modifie pas l'avance a l'allumage (cas de la figure 3), en d'autres 35 termes on maintient l'avance a l'allumage du cycle precedent. Ainsi, le procede selon !Invention permet de ne prendre que des actions correctives pertinentes : la commande de l'allumage n'est modifiee que si Ion est si- (par exemple a 10 99%, si l'intervalle de confiance represente 98 % des realisations de la distribution) que le quantile est en dehors de l'intervalle de confiance. La commande d'avance a l'allumage est par consequent plus stable, car des changements ne sont imposes qu'un cycle sur cinquante environ en regime stabilise. La reduction de la dispersion de la commande et de la variation 5 de l'etat moteur sur les courtes echelles de temps implique en retour une distribution plus stable des mesures KI, une meilleure estimation de ses parametres, et donc la possibilite de s'approcher de la limite cliquetis, d'obtenir un meilleur rendement tout en limitant le nombre de cycles cliquetants. 10 Selon une mise en oeuvre de ('invention, cette commande peut prendre en compte egalement la difference etablie a l'etape 3'). Pour cela, ('augmentation ou la diminution de l'avance a l'allumage (dans les situations i) et ii) decrites ci-dessus) peut etre proportionnelle a la difference entre I'esperance q du quantile predetermine determine a l'etape precedente et le seuil de cliquetis O. Ainsi, la commande est adaptee a l'intensite du cliquetis. 15 De plus, l'invention concerne un systeme de commande d'un moteur de combustion interne a allumage commande comprenant des moyens pour la mise en oeuvre du procede selon l'une quelconque des combinaisons de variantes du procede decrites precedemment. 20 En particulier, le systeme de commande peut comprendre : - un capteur de cliquetis, - des moyens de traitement du signal du capteur de cliquetis, - une memoire pour enregistrer N mesures (signal de cliquetis des N precedents cycles de combustion), 25 - des moyens de calcul pour construire l'estimateur bayesien et pour realiser la comparaison, et - des moyens de commande de l'avance a l'allumage. Les moyens de traitement du signal, la memoire, les moyens de calcul, et les moyens de commande peuvent etre integres au sein d'un calculateur embarque d'un vehicule. 30 L'invention concerne egalement un moteur de combustion interne equipe d'un tel systeme de commande. En outre, l'invention concerne un produit programme d'ordinateur telechargeable 35 depuis un reseau de communication et/ou enregistre sur un support lisible par ordinateur (calculateur embarque) et/ou executable par un processeur. Ce programme comprend des 11 instructions de code de programme pour la mise en oeuvre du procede tel que decrit cidessus, lorsque le programme est execute sur un ordinateur. Exemple comparatif 5 Les caracteristiques et avantages du procede selon l'invention apparaitront plus clairement a la lecture de l'exemple comparatif ci-apres. L'exemple comparatif est base sur la simulation d'un moteur a combustion interne 1.2 L comprenant trois cylindres a injection directe et a allumage commande. L'exemple comparatif concerne un regime stabilise du moteur a combustion interne a 1750 tr/min. 10 On rappelle que le cliquetis est un phenomene stochastique. Par consequent, merne sous des conditions identiques, chaque test est unique. Dans le but d'obtenir les intensites de cliquetis correspondant a des conditions reelles d'utilisation, les distributions experimentales des intensites de cliquetis sont determinees pour chaque point de fonctionnement, en utilisant tous les cycles precedents. 15 Dans le cadre de l'exemple comparatif, le procede selon l'invention est compare a un procede selon l'art anterieur, pour lequel l'avance a l'allumage augmente a chaque cycle moteur d'une faible valeur, et decroit fortement si la valeur de cliquetis mesuree excede un seuil predetermine. Classiquement, la faible valeur est de 0,01 ° CA (angle vilebrequin), et la forte decroissance correspond a 3° CA (angle vilebrequin). Cette strategie assure un temps 20 de reaction rapide pour tous les cycles cliquetants, tout en conservant 99,7% de cycles avec une intensite de cliquetis inferieure a la valeur seuil. Pour le procede selon l'invention, on utilise N=10 et le quantile est predetermine a 95%. 25 Les figures 6a et 6b illustrent les resultats obtenus par le procede selon l'invention. La figure 6a illustre l'avance a l'allumage A en °CA (angle vilebrequin) en fonction du temps T (s), et Ia figure 6b illustre les valeurs de KI sur Ia merne plage temporelle. Sur Ia figure 6b, le seuil 0 est symbolise par une droite horizontale, chaque cycle KI est represents par un point, et l'intervalle de confiance INT associe a chaque point est represents par un segment de 30 droite vertical. Comme attendu, avec le procede selon l'invention, l'avance a l'allumage n'est modifiee que lorsque le seuil est en dehors de l'intervalle de confiance. On observe sur Ia figure 6a que Ia commande est fres stable (peu de variation de l'avance a l'allumage). Les figures 7a et 7b illustrent les resultats obtenus par le procede selon l'art anterieur. La figure 7a illustre l'avance a l'allumage A en °CA (angle vilebrequin) en fonction du temps 35 T (s), et Ia figure 7b illustre les valeurs de KI sur Ia merne plage temporelle. Sur Ia figure 7b, 12 le seuil S est symbolise par une droite horizontale, et chaque cycle KI est represents par un point. On observe sur la figure 7a que la commande est peu stable avec de nombreuses modifications de l'avance a l'allumage, avec notamment des retraits d'avance brutal et des retours lents vers la zone cliquetante. 5 Ensuite, ces procedes sont compares en termes de consommation de carburant. II ressort de cette comparaison que le procede selon l'invention permet un gain de 1,1 % en consommation par rapport au procede selon l'art anterieur. 10 Des essais identiques ont ete merles pour d'autres regimes stabilises. Pour chaque essai, le procede selon l'invention permet un gain de consommation par rapport au procede selon l'art anterieur, ce gain est compris entre 0,6 et 2,2 %. Le procede selon l'invention permet donc de limiter le cliquetis, tout en recluisant la 15 consommation du moteur a combustion interne. 13 Annexe Estimation bayesienne de I'intervalle de confiance du quantile. Dans un but de simplification, cette annexe concerne ('estimation bayesienne pour une distribution normale. Cette estimation bayesienne peut etre adaptee pour toute distribution 5 pouvant etre rapportee a une distribution normale, en particulier a une distribution log- normale. I. MOTIVATIONS On suppose que Ion dispose de n mesures du phenomene physique etudie, dont les 10 caracteristiques aleatoires ne peuvent etre negligees. On admet par ailleurs que ces mesures sont independantes et normalement distribuees. On dispose alors de n realisations qu'on note X1,... ,Xn. Soit q le vrai quantile 95% de cette distribution. La formule de Bayes permet d'exprimer la distribution de la valeur de q, au regard des mesures observees X1,... ,Xn. Naturellement, q peut etre estime sans incertitude lorsque n 09. Neanmoins, lorsque n 15 REVENDICATIONS1) Procede de commande d'un moteur a combustion interne a allumage commande, ledit moteur a combustion interne comprenant un capteur de cliquetis agence sur ledit moteur a combustion interne, caracterise en ce que ledit procede comporte les etapes suivantes : a) on realise N mesures (MES) de cliquetis au moyen dudit capteur de cliquetis ; b) on construit un estimateur (EST) d'une distribution desdites N mesures, ledit estimateur determinant un intervalle de confiance -inn) qmax) d'un quantile predetermine au moyen desdites N mesures ; c) on commande (COM) I'allumage au sein dudit moteur a combustion interne en fonction d'une comparaison avec ledit intervalle de confiance (611 nn) ,max, n -i i) si ledit seuil de cliquetis 6 est inferieur a Ia valeur minimale qm,n dudit intervalle de confiance dudit quantile predetermine, alors on diminue I'avance a I'allumage ; ii) si ledit seuil de cliquetis 6 est superieur a Ia valeur maximale qmax dudit intervalle de confiance dudit quantile predetermine, alors on augmente I'avance a I'allumage ; et iii) si ledit seuil de cliquetis 6 est compris dans ledit intervalle de confiance (qmin, qmax), alors on ne modifie pas I'avance a I'allumage. 2) Procede de commande selon la revendication 1, dans lequel ledit estimateur (EST) est un estimateur bayesien. 3) Procede de commande selon l'une des revendications precedentes, dans lequel ladite distribution est une distribution pouvant etre rapportee a une distribution normale par une methode mathematique, de preference ladite distribution est une distribution log-normale. 4) Procede de commande selon l'une des revendications precedentes, dans lequel ledit estimateur determine en outre l'esperance q dudit quantile predetermine. 5) Procede de commande selon la revendication 4, dans lequel on diminue ou on augmente I'avance a I'allumage de maniere proportionnelle a la difference (DIF) entre ladite esperance dudit quantile predetermine et ledit seuil de cliquetis. 6) Procede de commande selon l'une des revendications precedentes, dans lequel N est compris entre 5 et 20, de preference entre 7 et 15. 17 7) Procede de commande selon I'une des revendications precedentes, dans lequel ledit intervalle de confiance contient la valeur du quantile reel dans exactement 98% des cas. 8) Procede de commande selon I'une des revendications precedentes, dans lequel ledit 5 capteur de cliquetis est un accelerometre agence sur la culasse dudit moteur a combustion interne. 9) Procede de commande selon la revendication 5, dans lequel ladite mesure de cliquetis est le maximum d'amplitude du signal filtre dudit accelerometre. 10 10) Procede de commande selon I'une des revendications precedentes, dans lequel ledit seuil de cliquetis est obtenu par essais au banc moteur durant la phase de calibration du moteur. 15 11) Procede de commande selon I'une des revendications precedentes, dans lequel ledit quantile predetermine de ladite distribution est superieur ou egal a 80%. 12) Systeme de commande d'un moteur a combustion interne a allumage commande comprenant des moyens pour mettre en oeuvre le procede selon I'une des revendications 20 precedentes. 13) Produit programme d'ordinateur telechargeable depuis un reseau de communication et/ou enregistre sur un support lisible par ordinateur et/ou executable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procede 25 selon I'une des revendications 1 a 11, lorsque ledit programme est execute sur un ordinateur.