FR 3064310 A1 20180928 FR 1752457 A 20170324 1 La présente invention concerne les procédés pour détecter et évaluer un niveau de cliquetis dans un cylindre d'un moteur à combustion interne à allumage commandé. Dans les moteurs à allumage commandé, une étincelle de bougie est générée 5 par une impulsion électrique commandée par un calculateur électronique, cette impulsion étant positionnée par rapport au point mort haut avec une certaine avance à l'allumage. Dans certains cas, le mélange compressé peut faire l'objet d'un auto-allumage, phénomène à éviter. Une avance à l'allumage trop prononcée ou un auto-allumage conduisent à un 10 phénomène connu de cliquetis (aussi appelé 'cognement' ou `knock' en anglais). C'est dans ce contexte que les systèmes de contrôle moteur font maintenant appel à une fonction de détection de cliquetis, grâce à un capteur acoustique, notamment piézoélectrique, installé sur le carter cylindre. La présence de signaux significatifs dans une certaine gamme de fréquence, 15 typiquement entre 5 kHz et 25 kHz, est caractéristique de la présence d'un tel phénomène de cliquetis. Si du cliquetis est détecté, alors l'avance à l'allumage doit être réduite. Toutefois, la détection du phénomène de cliquetis peut être polluée par des bruits auxiliaires parasites. Par exemple, il a été découvert que la fermeture d'un injecteur provoque un bruit qui comporte des composantes fréquentielles dans la gamme spectrale 20 d'intérêt pour le cliquetis. Suivant les différentes conditions de fonctionnement du moteur, la position angulaire de la fermeture de chaque injecteur évolue. Pour la détection du cliquetis, pour chaque cylindre, on prévoit généralement une fenêtre temporelle d'observation, pour éliminer les bruits parasites qui ne coïncident 25 pas temporellement avec le moment où le cliquetis peut apparaitre. Mais, pour certaines conditions de fonctionnement du moteur, il s'avère qu'un bruit de fermeture d'injecteur peut 'tomber' dans la fenêtre temporelle d'observation et peut provoquer des non détection ou des détections à tort. Le document US2012192835A1 propose une solution visant à réduire ce 30 problème. Toutefois, dans le cas où la position angulaire de la fermeture de l'injecteur évolue rapidement ou de manière discontinue, il reste des cas de non détection et/ou des cas de détection à tort. Il subsiste donc un besoin d'améliorer les solutions existantes pour proposer un procédé d'évaluation de cliquetis qui ne soit pas perturbé par le fait qu'un bruit 3064310 2 parasite (tel que le bruit de fermeture d'injecteur) puisse 'tomber' dans la fenêtre temporelle d'observation du cliquetis. À cet effet, il est proposé ici un procédé pour gérer (à savoir détecter / évaluer / corriger) du cliquetis dans un cylindre d'un moteur à combustion interne, le procédé étant mis en oeuvre dans un système comprenant au moins un capteur acoustique, un calculateur (unité de calcul), le procédé étant destiné à prendre en compte une pollution acoustique provenant d'un bruit parasite répétable connu et de position temporelle variable, le procédé comprenant : - mettre en forme et numériser les signaux en provenance du capteur acoustique, - appliquer un filtre passe bande pour ne conserver que la gamme de fréquences d'intérêt pour le cliquetis, dans un bruit filtré, - déterminer une fonction de correction à gain ajustable utilisant une courbe de correction de gain, destinée à diminuer l'influence du bruit parasite, - déterminer, en fonction de la position angulaire de fin d'injection, le point de courbe de correction de gain à utiliser pour transformer le bruit filtré en score de cliquetis corrigé, - comparer un score de cliquetis corrigé ainsi obtenu à un seuil de décision de cliquetis, - en déduire une valeur de correction d'avance d'allumage à appliquer au prochain cycle, caractérisé en ce que la courbe de correction de gain est définie par une valeur de calibration et quatre points angulaires, lesdits quatre points angulaires étant obtenus par calcul, à partir des positions de début et fin de la fenêtre d'observation de cliquetis et au moins d'une caractéristique connue du bruit. Grâce aux dispositions ci-dessus, on peut diminuer très fortement l'influence d'un bruit parasite qui tombe dans la fenêtre d'observation, sans modification importante de la structure logicielle préexistante. La correction de gain fonctionne immédiatement sans aucun délai lorsque la fermeture d'injecteur évolue rapidement ou par saut et tombe dans la fenêtre 30 d'observation. Dans divers modes de réalisation de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes. Selon une option, le calcul des quatre points angulaires est effectué à partir d'une caractéristique de durée de bruit et d'une caractéristique de temps de transmission 35 de bruit vers le capteur acoustique ces deux valeurs étant prises comme caractéristiques connues du bruit. 3064310 3 Bénéfice ainsi obtenu : en connaissant seulement deux paramètres relatifs au bruit parasite, on arrive à construire la courbe de correction qui va permettre d'atténuer l'effet de ce bruit parasite. Selon une option, le bruit parasite auquel on s'intéresse est un bruit de 5 fermeture d'injecteur. Il s'avère que ce bruit est parfois prépondérant car il provient d'une source proche, directement fixée sur le carter cylindre En pratique, il s'agit du bruit de fermeture de l'injecteur du cylindre suivant dans le cycle, c'est-à-dire celui qui suit le cylindre en cours d'explosion dans la séquence de fonctionnement du moteur. 10 Selon une option, le temps de transmission du bruit est issu d'un ensemble de paramètres de temps de transmission, un paramètre pour chaque couple (cylindre en fin d'injection, cylindre en cours d'explosion). Avantageusement on prend ainsi en compte la configuration géométrique du moteur avec les temps de transmission qui peuvent être différents pour chaque couple de cylindres successifs dans la séquence. 15 Selon une option, la valeur de calibration représente le plateau bas de la courbe de correction de gain entre les deuxième et troisième points angulaires, et la courbe de correction de gain prend la valeur 1 pour les points antérieurs au premier point angulaire et pour les points postérieurs au quatrième point angulaire. Selon une option, on utilise une interpolation linéaire pour obtenir les valeurs 20 de la courbe de correction entre les premier et deuxième points angulaires, et entre les troisième et quatrième points angulaires. Avantageusement, ce calcul simple est très rapide et mobilise très peu de ressources du microcontrôleur. Selon une option, on utilise, pour calibrer la valeur de calibration, une calibration simple à deux dimensions selon la charge moteur et le régime moteur. 25 Avantageusement, cette calibration simple est très peu gourmande en mémoire. Selon une option, les valeurs de calibration sont comprises entre 0,25 et 1. Selon une option, les valeurs de calibration décroissent au fur et à mesure que la charge moteur augmente. Selon une option, on calcule les abscisses Xa, Xb, Xc, Xd des quatre points 30 angulaires comme suit : Xa = WS1 - (Tbr + Ttr)*6*RPM Xb = WS1 - (Ttr)*6*RPM Xc = WS2 - (Tbr + Ttr)*6*RPM Xd = WS2 - (Ttr)*6*RPM 35 Où : 3064310 4 WS1 et WS2 représentent respectivement les positions de début et de fin de la fenêtre d'observation de cliquetis, exprimées en degrés de rotation du vilebrequin, Tbr représente une caractéristique de durée de bruit exprimée en 5 seconde, Ttr représente une caractéristique de temps de transmission de bruit vers le capteur acoustique exprimée en seconde, RPM représente le régime moteur exprimé en tours par minute. Selon une option, le seuil de décision de cliquetis KSC est déterminé par 10 l'expression KSC = G*KS + 1 - G, où : KS représente le seuil en l'absence de bruit, et G représente la valeur sur la courbe de correction au point de l'angle de fermeture de l'injecteur. Selon une option, le procédé est effectué en temps réel, pour chaque cylindre 15 et à chaque cycle. Avantageusement, il n'y a aucun délai de prise en compte, par exemple si la fermeture injecteur connaît un saut d'angle, la prise en compte de la nouvelle situation est immédiate. Selon une option, le filtre passe bande est de [5 kHz - 25 kHz]. On élimine ainsi avantageusement tous les bruits qui se trouvent en dehors de cette bande. 20 L'invention vise aussi un système pour gérer du cliquetis dans un cylindre d'un moteur à combustion interne, le système comprenant au moins un capteur acoustique et un calculateur caractérisé en ce que ce dernier est configuré pour mettre en oeuvre le procédé tel que décrit ci-dessus. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de 25 la description suivante donnée à titre d'exemple non limitatif, en regard des dessins joints, sur lesquels : les figures 1 et 1A montrent schématiquement le système et les composants dans lesquels le procédé selon l'invention est mis en oeuvre, la figure 2 représente un chronogramme illustratif, la figure 3 illustre la courbe de gain de correction, la figure 4 illustre une cartographie en deux dimensions de correction, la figure 5 montre un bloc-diagramme schématique du processus de calcul de correction, la figure 6 montre un diagramme illustrant les étapes du procédé. Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. Pour des raisons de clarté de l'exposé, certaines éléments temporels ne sont pas représentés à l'échelle. 3064310 5 Sur la figure 1, on a représenté un calculateur 1 pour gérer un moteur à explosion, typiquement un moteur à combustion interne à essence. L'exemple illustré est basé sur un moteur à quatre cylindres, mais le procédé peut tout aussi bien être appliqué à un moteur à trois cylindres, voire un moteur à six 5 cylindres en ligne ou en V. Il faut noter que l'invention pourrait être appliquée à un moteur fonctionnant au GPL ou Gaz de Pétrole Liquéfié, et plus généralement à tout moteur à combustion interne à allumage commandé. Chacun des cylindres est équipé d'un injecteur de carburant avec un nez 10 d'injection donnant directement dans la chambre de combustion : un injecteur 31 dans le premier cylindre CYL C1, un injecteur 32 dans le deuxième cylindre CYL C2, un injecteur 33 dans le troisième cylindre CYL C3, un injecteur 34 dans le quatrième cylindre CYL C4. Sur le carter-cylindre 4 (autrement appelé « bloc-moteur »), est installé un 15 capteur de cliquetis 2, relié par un câble, en général blindé au calculateur moteur. Ce capteur est un capteur acoustique, de préférence de type piézoélectrique ; il est configuré pour être sensible aux vibrations produites dans le carter-cylindre. On note que le bloc moteur 4 peut être en fonte ou en alliage d'aluminium. Sur l'exemple illustré, le carter-cylindre est équipé d'un seul capteur 20 acoustique de cliquetis, mais il pourrait y en avoir plusieurs ; dans le cas des moteurs à deux rangées de cylindres comme par exemple les moteurs à six cylindres en V ou V6 ou à huit cylindres en V ou V8, il y a au moins deux capteurs de cliquetis, un par rangée de cylindres. Comme connu en soi, le fonctionnement du moteur est basé sur une 25 séquence ordonnée de temps (admission, compression, explosion, échappement) pour chacun des cylindres, l'ordre successif des cylindres étant donné par construction. Dans l'exemple illustré du moteur à quatre cylindres, l'ordre habituel d'allumage est CYL Cl -CYL C3-CYL C4-CYL C2. Comme connu en soi, la position angulaire du vilebrequin (repérée 30 génériquement par (p) est connue grâce à une cible codée ou dentée 14, avec un capteur 11 de position du volant moteur. On illustre à la figure 2 un chronogramme qui montre en bas une portion temporelle centrée sur l'explosion du cylindre N (CYL N) et qui montre en haut la phase d'injection d'essence sur le cylindre suivant dans la séquence repéré cylindre N+1. 35 Sur le cylindre N+1, un signal électrique 36 commandé par le calculateur provoque l'ouverture de l'injecteur à la position angulaire (pOp et la fermeture de l'injecteur à la position angulaire (pCI. Cette fermeture provoque des ondes acoustiques, autrement 3064310 6 dit un bruit qui va se propager dans l'ensemble du carter-cylindre. Comme déjà mentionné en introduction, le bruit de cette fermeture peut coïncider temporellement avec la fenêtre d'observation WS du cliquetis sur le cylindre N. Cette fenêtre d'observation commence à la position angulaire WS1 et se 5 termine à la position angulaire WS2. Le début WS1 de la fenêtre d'observation et la fin WS2 de la fenêtre d'observation WS font l'objet d'une calibration préexistante, connue en soi non décrite en détail ici. Selon la commande 36, 37 dictée par le point de fonctionnement moteur (cartographié), le bruit de la fermeture injecteur peut tomber juste avant la fenêtre 10 d'observation (repère 5), il peut tomber à cheval sur le début de la fenêtre d'observation (repère 51), il peut tomber totalement dedans la fenêtre d'observation (repères 52); il peut tomber à cheval sur la fin la fenêtre d'observation (repère 53) ou après la fenêtre d'observation (repère 54). Il faut remarquer que sur la figure 2, les abscisses représentent des angles de 15 position vilebrequin, ce qui correspond aussi au temps qui passe moyennant la prise en compte du facteur régime moteur (noté RPM dans la suite). On comprend que le bruit de fermeture de l'injecteur du cylindre N+1 peut donc représenter un bruit parasite pour l'observation du cliquetis sur le cylindre N. Ce bruit a une durée connue, de l'ordre de la milliseconde, la durée notée Tbr 20 est portée de préférence dans un paramètre de calibration. Il n'est pas exclu d'avoir une durée spécifique par cylindre. Une autre caractéristique importante pour la suite est la connaissance du temps de transmission du bruit à partir de sa source jusqu'au capteur de cliquetis 2. Ce temps de transmission est noté Ttr. On comprend, à l'aune de la figure 1, que les temps 25 de transmission Ttr suivant chaque cylindre sont différents voire très différents, du fait notamment des distances différentes. Par conséquent, il est prévu de stocker un ensemble de paramètres de temps de transmission, préférentiellement un paramètre pour chaque couple (cylindre en fin d'injection, cylindre en cours d'explosion). Dans l'exemple illustré, on peut se limiter à 30 quatre paramètres ; mais on peut avoir un paramétrage plus complexe notamment si on utilise un second capteur de cliquetis. Les inventeurs ont remarqué que le bruit de fermeture injection est connu et répétable ; la durée de ce bruit et le temps de transmission pour qu'il arrive jusqu'au capteur de cliquetis sont des paramètres simples connus à l'avance. 35 En revanche, ce qui n'est pas connu à l'avance mais qui est déterminé en temps réel à chaque cycle d'injection dans un cylindre, c'est la position angulaire (pCI correspondant au moment de fermeture d'injecteur. Cette donnée est calculée en temps 3064310 7 réel à partir d'une calibration complexe préexistante (déjà disponible pour le fonctionnement conventionnel du moteur). Lorsque le bruit injecteur coïncide partiellement ou totalement avec la fenêtre d'observation du cliquetis, cela engendre un accroissement qui n'est pas souhaité du 5 niveau de bruit dans la gamme de fréquences d'intérêt (c'est-à-dire dans la bande [5 kHz-25 kHz]). Les inventeurs ont astucieusement proposé d'utiliser une courbe de correction de gain 7 pour réduire fortement les effets de cet accroissement indésirable. Cette courbe de correction de gain permet de choisir, en fonction de l'angle 10 auquel se produit la fermeture de l'injecteur, un gain G ajusté en temps réel entre la valeur de base 1 est une valeur de gain où la réduction est la plus forte à savoir Gin. Comme illustré à la figure 3, la définition de la courbe de correction repose sur la détermination de quatre points angulaires (A, B, C, D). On voit que cette courbe de correction prend la valeur unité (c'est-à-dire 1) 15 avant un premier point angulaire noté A (segment 71). Ceci correspond au cas où la fermeture d'injecteur provoque un bruit dont les effets se situent intégralement avant la fenêtre d'observation WS. Cette courbe de correction prend également la valeur 1 après un quatrième point angulaire noté D (segment 75). Ceci correspond au cas où la fermeture d'injecteur provoque un bruit dont les effets se situent intégralement après la 20 fenêtre d'observation WS. On voit que cette courbe de correction prend la valeur notée Gin entre le deuxième point angulaire noté B et le troisième point angulaire noté C. Dans ce segment repéré 73, la courbe forme un plateau prenant la valeur Gin. Nous verrons plus loin comment est déterminée la valeur de calibration Gin. 25 Entre le premier point angulaire noté A et le deuxième point angulaire noté B, la courbe 7 est rectiligne 72, autrement dit il s'agit d'une interpolation linéaire entre les premier et deuxième points angulaires. De même entre le troisième point angulaire noté C et le quatrième point angulaire noté D, la courbe 7 est rectiligne 74, autrement dit il s'agit d'une interpolation linéaire entre les troisième et quatrième points angulaires. 30 Bien entendu, au lieu d'une interpolation linéaire, on pourrait avoir une fonction différente qui prenne en compte la dissymétrie du bruit par exemple une différence de puissance entre le début et la fin du bruit. Sur la figure 3, on a représenté une autre instance de courbe de correction de gain en trait mixte avec le repère 7'. Il faut en effet noter qu'à chaque cycle moteur, les 35 bornes WS1, WS2 de la fenêtre d'observation peuvent bouger ainsi que le régime moteur RPM. Autrement dit on fait le calcul en temps réel à chaque cycle d'explosion sur un des cylindre. 3064310 8 La figure 4 illustre la table de calibration qui permet d'obtenir la valeur Gin en fonction du point de fonctionnement moteur ; il s'agit d'une calibration 6 en deux dimensions, à savoir on obtient une valeur scalaire Gin dans un espace à deux dimensions à savoir charge moteur CHGMOT, et régime moteur RPM. Ce type de 5 calibration est très courant et donc non détaillé plus avant ici. On remarque tout d'abord que les valeurs de Gin sont toujours inférieures ou égales à 1 ; par ailleurs, dans l'exemple illustré, les valeurs Gin sont comprises entre 0,25 et 1, ce qui permet d'aller jusqu'à un facteur de réduction de quatre pour atténuer le bruit parasite présent totalement dans la fenêtre observation WS. Des valeurs plus basses 10 pour Gin ne sont pas pour autant exclues. La courbe de correction de gain peut donc être définie très simplement par les quatre points de définition avec leurs abscisses et leurs ordonnées A(Xa, 1), B(Xb, Gin), C(Xc, Gin) et D(Xd, 1). Pour calculer les quatre abscisses d'intérêt, on applique les formules 15 suivantes avec RPM exprimé en tours par minute, Tbr et Ttr exprimés en seconde, et les autres valeurs étant exprimées en degrés de rotation du vilebrequin : Xa = WS1 - (Tbr + Ttr)*6*RPM Xb = WS1 - (Ttr)*6*RPM Xc = WS2 - (Tbr + Ttr)*6*RPM 20 Xd = WS2 - (Ttr)*6*RPM Sur la figure 5, sont illustrés le calcul de Gin qui fait appel aux données charge moteur CHGMOT, et régime moteur RPM ; au-dessus se trouve le bloc de calcul des quatre abscisses Xn comme illustré ci-dessus ; le bloc repéré 16 positionne l'angle de fin d'injection (pCI sur la courbe de correction de gain et identifie le gain G instantané à 25 utiliser pour les corrections du cycle d'observation de cliquetis en cours. Le score de bruit brut constaté noté Bi est multiplié par le gain instantané G, ce qui donne un score de bruit corrigé noté KNK SCORE CORR. Il faut noter que le seuil de référence utilisé en l'absence de bruit, est ici corrigé de manière à conserver une sensibilité suffisante. 30 On fonctionne alors avec un seuil corrigé KNK SEUIL CORR (KSC en bref) qui forme le seuil de décision de cliquetis et qui est déterminé par la formule : KSC = G*KS + 1 - G, dans laquelle KS représente le seuil conventionnel d'absence de bruit. Ce seuil conventionnel KS n'est généralement pas une valeur absolue ; il s'agit d'un multiple de la 35 valeur moyenne glissante du bruit constaté en l'absence de cliquetis, par exemple trois fois la valeur moyenne glissante du bruit constaté en l'absence de cliquetis. 3064310 9 La décision de correction d'avance à l'allumage se prend au final sur la comparaison entre le score de bruit corrigé KNK SCORE CORR et le seuil de décision corrigé KNK SEUIL CORR (KSC). La figure 6 illustre les différentes étapes du procédé, à savoir la mise en forme 5 et numérisation 82 qui transforme les signaux analogiques 20 en données numériques 22 ; le bloc de filtrage passe bande 84 qui fournit des signaux numériques filtrés 24 à partir desquelles est calculé le score de bruit brut Bi mentionné ci-dessus ; le bloc de correction de gain 86 a été détaillé ci-dessus ; il fournit un score de bruit corrigé ; le bloc de comparaison 88 fournit une proposition de correction d'avance d'allumage 90. 10 REVENDICATIONS1. Procédé pour gérer du cliquetis dans un cylindre d'un moteur à combustion interne, le procédé étant mis en oeuvre dans un système comprenant au moins un capteur acoustique (2) et un calculateur (1), le procédé étant destiné à prendre en compte une pollution acoustique provenant d'un bruit parasite (5) répétable connu et de position 5 temporelle variable, le procédé comprenant : - mettre en forme et numériser (82) les signaux (20) en provenance du capteur acoustique, - appliquer un filtre passe bande (84) pour ne conserver que la gamme de fréquences d'intérêt pour le cliquetis, dans un bruit filtré (24), 10 - déterminer une fonction de correction (86) à gain ajustable utilisant une courbe de correction de gain (7), destinée à diminuer l'influence du bruit parasite, - déterminer, en fonction de la position angulaire de fin d'injection, le point de courbe de correction de gain à utiliser pour transformer le bruit filtré en score de cliquetis corrigé, 15 - comparer (88) un score de cliquetis corrigé ainsi obtenu à un seuil de décision de cliquetis, - en déduire (90) une valeur de correction d'avance d'allumage à appliquer au prochain cycle, caractérisé en ce que la courbe de correction de gain est définie par une valeur de 20 calibration (Gin) et quatre points angulaires (A, B, C, D), lesdits quatre points angulaires étant obtenus par calcul, à partir des positions de début et fin (WS1, WS2) de la fenêtre d'observation de cliquetis et au moins d'une caractéristique connue du bruit. 2. Procédé selon la revendication 1, dans lequel le calcul des quatre points angulaires est effectué à partir d'une caractéristique de durée de bruit (Tbr) et d'une 25 caractéristique de temps de transmission (Ttr) de bruit vers le capteur acoustique (2), ces deux valeurs étant prises comme caractéristiques connues du bruit. 3. Procédé selon l'une des revendications 1 à 2, dans lequel le bruit parasite est un bruit de fermeture d'injecteur. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel le 30 temps de transmission (Ttr) du bruit est issu d'un ensemble de paramètres de temps de transmission, un paramètre pour chaque couple (cylindre en fin d'injection, cylindre en cours d'explosion). 3064310 11 5. Procédé selon l'une des revendications 1 à 4, dans lequel la valeur de calibration (Gin) représente le plateau bas de la courbe de correction de gain entre les deuxième et troisième points angulaires (B, C), et la courbe de correction de gain prend la valeur 1 pour les points antérieurs au premier point angulaire (A) et pour les points postérieurs au quatrième point angulaire (D). 6. Procédé selon l'une des revendications 1 à 5, dans lequel on utilise, pour calibrer la valeur de calibration (Gin), une calibration simple (6) à deux dimensions selon la charge moteur et le régime moteur. 7. Procédé selon la revendication 2, dans lequel on calcule les quatre points 10 angulaires (Xa, Xb, Xc, Xd) comme suit : Xa = WS1 - (Tbr + Ttr)*6*RPM Xb = WS1 - (Ttr)*6*RPM Xc = WS2 - (Tbr + Ttr)*6*RPM Xd = WS2 - (Ttr)*6*RPM, 15 20 WS1 et WS2 représentant respectivement les positions de début et de fin de la fenêtre d'observation de cliquetis, exprimées en degrés de rotation du vilebrequin, Tbr représentant une caractéristique de durée de bruit exprimée en seconde, Ttr représentant une caractéristique de temps de transmission de bruit vers le capteur acoustique exprimée en seconde, RPM représentant le régime moteur exprimé en tours par minute. 8. Procédé selon la revendication 1, dans lequel le seuil de décision de cliquetis (KSC) est déterminé par l'expression : KSC = G*KS + 1 - G, KS représentant le seuil en l'absence de bruit, et 25 - G représentant la valeur sur la courbe de correction au point de l'angle de fermeture de l'injecteur. 9. Procédé selon l'une des revendications 1 à 8, caractérisé en ce qu'il est effectué en temps réel, pour chaque cylindre et à chaque cycle. 10. Système pour gérer du cliquetis dans un cylindre d'un moteur à combustion 30 interne, le système comprenant au moins un capteur acoustique, un calculateur caractérisé en ce qu'il est configuré pour mettre en oeuvre le procédé selon l'une des revendications 1 à 9.