Brevet FR3000206A1 — publié le 20140627

Domaine de l'invention La présente invention concerne les essais sur les moteurs à combustion interne visant à caractériser leur fonctionnement, par exemple à des fins de calibration, lesdits essais étant par exemple réalisés sur un banc moteur de tout type - stationnaire, dynamique ou haute- dynamique - ou réalisés sur le véhicule en fonctionnement. Contexte général Les essais moteurs consistent à mesurer et enregistrer les réponses d'un moteur placé dans des conditions particulières de fonctionnement, ces conditions pouvant être simulées ou réelles, et par exemple correspondre à des conditions de conduite particulières d'un véhicule en terme environnemental (température d'air extérieur, altitude, pente de la route) ou en terme d'utilisation (roulage urbain ou autoroutier, comportement sportif ou économique, masse remorquée par le véhicule). Les informations obtenues permettent alors de caractériser le fonctionnement du moteur en relation avec les conditions particulières appliquées, et peuvent par exemple être utilisées pour calibrer les variables des stratégies de contrôle du moteur programmées dans le calculateur électronique afin notamment de répondre à des normes et critères divers, tel que des normes concernant les émissions de polluants, le bruit, la souplesse de la conduite, etc. Ces essais sont traditionnellement réalisés sur des bancs de tests moteur, aussi appelés bancs d'essais moteurs ou bancs moteurs, par exemple lors des phases de conception des moteur, de réglage avant mise dans le commerce, ou tout au long de la production industrielle du moteur. Chaque essai se compose généralement de nombreuses mesures de différents paramètres de fonctionnement du moteur ou de ses organes, ces paramètres reflétant leur état de fonctionnement, en particulier l'état de la combustion du moteur et les conséquences thermiques, mécaniques, cinétiques et chimiques de cette réaction. Les paramètres de fonctionnement sont par exemple des concentrations de composés chimiques dans les gaz d'échappement, par exemple des substances polluantes, la pression et la température de ces gaz, le bruit généré par la combustion dans le moteur, la vitesse de rotation de l'arbre du turbo-compresseur, etc. Lors d'une séquence d'essais moteur, on fait classiquement varier à plusieurs reprises les réglages du moteur (pression d'injection, phasage d'injection, pression d'air, etc), et/ou le point de fonctionnement du moteur (régime, charge). Le point de fonctionnement du moteur correspond à un couple (régime, charge), et sera également dénommé point de puissance dans la présente description, en relation avec la puissance du moteur qui s'exprime commé étant le produit de la vitesse de rotation du moteur (régime) par le couple (charge). Parmi les mesures effectuées lors de l'essai moteur, on peut distinguer les mesures continues des mesures ponctuelles, les mesures continues représentant la très grande majorité des mesures généralement effectuées. Les mesures continues consistent à mesurer la réponse du moteur suite à chaque variation des réglages et/ou du point de puissance du moteur. Les mesures continues sont des mesures réalisées par des dispositifs dits synchrones, réalisant des mesures à intervalles réguliers et présentant un taux de disponibilité élevé. Ce type de mesures est réalisé avec un fréquence d'échantillonnage relativement rapide (typiquement de 1 à 10Hz): les observations sont rafraîchies dans le temps toutes les secondes ou même plus souvent. Les mesures continues sont par exemple des mesures permanentes, qui sont soit moyennées sur une fenêtre temporelle glissante, soit brutes, telles que les mesures de pression, de température, de vitesse, de tension, d'opacité des gaz de combustion, de pression cylindre permettant d'établir des grandeurs décrivant la combustion. D'autres mesures continues, telles que la mesure de concentration en espèces chimiques des gaz de combustion, la mesure instantanée de la consommation de carburant, peuvent parfois être interrompues pendant des opérations de purge, permettant d'éliminer certains composés chimiques imprégnant les lignes de prélèvement de gaz afin de restaurer la qualité de mesure et d'éviter les dérives, ou pendant des opérations de remplissage (certains appareils sont constitués d'un seau posé sur une balance pour peser le carburant qui débite dans le moteur). Les mesures ponctuelles sont quant à elles réalisées par des dispositifs dits asynchrones qui réalisent une observation à un moment donné et sur une plage de temps limitée. Après chaque observation, la valeur n'évolue pas dans le temps. Il peut s'agir de la mesure d'indice de fumée par exposition d'un filtre à un volume de gaz de combustion prélevé à l'échappement, de la mesure intégrée de la consommation de carburant, effectuée en masse sur une durée écoulée ou en temps pour une masse donnée de carburant consommée. L'ensemble des mesures lors d'une séquence d'essais moteurs permet alors de caractériser le fonctionnement du moteur dans les conditions particulières des essais. Traditionnellement, une mesure continue lors d'un essai moteur, c'est-à-dire la mesure continue d'une grandeur pour des paramètres de réglages et/ou un point de puissance donné du moteur, se décompose en trois étapes, tel qu'illustré à la figure I. La mesure débute avec une première phase 11 de modification rapide (en quelques secondes au plus, représentée par l'intervalle de temps t1-t2 dans la figure 1) des réglages du moteur (PR), par exemple le débit d'air frais à l'admission (Qaid ou la quantité de carburant injectée pour l'injection pilote (Qin; Pi!), ou de son point de puissance, par exemple le régime et/ou la charge. En réaction, les signaux acquis des paramètres de fonctionnement PF du moteur varient en général fortement. La phase suivante est une phase d'attente de la stabilisation du signal 12 (intervalle de temps t243), qui doit être la même pour toutes les mesures effectuées car généralement l'objectif est de mesurer des phénomènes qui se déroulent au même moment. Elle s'achève lorsque tous les paramètres de fonctionnement du moteur ont une variation plus faible qu'un seuil propre à chaque paramètre (quelques pourcents) pendant une durée donnée (environ 10 secondes). Dans le cadre d'un essai sur banc moteur, cette deuxième phase dure typiquement 45 secondes pour obtenir une bonne qualité de mesure. Une troisième phase correspond à une phase de moyennage du signal 13. Cette troisième phase de mesure 13 permet de calculer une valeur représentative de la mesure continue complètement stabilisée, en effectuant la moyenne de tous les échantillons enregistrés dans l'intervalle de temps correspondant t3-t4. Il est courant, dans le cadre d'un essai au banc moteur, que cette phase dure 30 secondes afin d'obtenir un résultat fiable. La mesure continue effectuée selon les différentes phases décrites lors d'un essai moteur est connue dans le domaine des essais moteurs sous le nom de prise de point. Une séquence d'essais moteur permet ainsi l'enregistrement d'une information synthétique, par exemple dans un fichier informatique, où à chaque mesure continue réalisée lors de l'essai moteur correspond une valeur unique représentant la moyenne de tous les échantillons de cette mesure continue sur l'intervalle de temps de la phase de moyennage 13. Les valeurs sont indicées suivant un numéro de point. Cette information synthétique peut alors contribuer à établir une base de données relative à la séquence d'essais moteur réalisée, cette base de donnée pouvant être traitée avec des outils numériques de manière à fournir des informations supplémentaires concernant l'état de fonctionnement du moteur, tels que des modèles comportementaux statistiques (régis par exemple par des fonctions simples telles que des polynômes) ou physiques (régis par des systèmes d'équations complexes à résoudre). Un problème majeur lié à cette prise de point traditionnelle réside dans le fait que les phases de stabilisation et de moyennage sont très consommatrices de temps lors de l'essai moteur. Comme cela est représenté à la figure 1, ces deux phases comptent pour plus de la moitié de la durée d'une séquence de mesure continue lors de l'essai moteur, et souvent plus de 80 % de la durée totale de l'essai moteur. Il est très difficile de réduire la durée de la phase de stabilisation sans prendre le risque d'exposer la phase suivante de moyennage à des phénomènes thermo-mécaniques ayant des inerties et/ou des dynamiques différentes pouvant conduire à invalider cette phase de moyennage, et au final aboutir à devoir répéter la mesure et, par conséquent, allonger la durée de l'essai. De même, réduire la durée de la phase de moyennage peut introduire un biais dans le résultat. Un autre problème rencontré lors de cette prise de point est que la phase de moyennage, déterminée de manière unique pour toutes les mesures continues, ne fournit pas toujours une valeur représentative du signal stabilisé. En effet, il peut y avoir certaines mesures qui varient encore, là où d'autres se sont stabilisées. C'est par exemple le cas de la mesure de la concentration en hydrocarbures imbrulés (HC) dans la figure 1, qui continue d'évoluer lors de la phase de moyennage alors que la mesure de la concentration de monoxyde de carbone (CO) est stabilisée depuis beaucoup plus longtemps. Ce genre de phénomène est relativement fréquent au cours des essais moteurs. La présente invention vise de manière générale à réduire la durée des essais moteurs, et plus particulièrement le délai d'attente lié aux phases de stabilisation et de moyennage du signal, traditionnellement requis pour calculer une valeur représentative d'un signal une fois stabilisé. L'invention a également pour objectif de fournir un procédé fiable pour caractériser le fonctionnement d'un moteur à combustion interne lors d'une séquence d'essais moteur, en ce qu'on fournit pour chaque mesure continue une valeur représentative du signal stabilisé. Pour cela, l'invention propose un procédé de caractérisation du fonctionnement d'un moteur à combustion interne, à partir d'une séquence d'essais moteur, chaque essai moteur comprenant la modification d'au moins un paramètre de réglage du moteur et/ou de son point de puissance, comprenant les étapes suivantes: a) pour chaque essai moteur: i. on acquiert au moins un signal issu de la mesure continue d'un paramètre de fonctionnement PR dudit moteur sur une durée d'acquisition d, la durée d'acquisition d étant choisie de manière à pouvoir déterminer au moins une fonction sigmoïde pour modéliser le dit au moins un signal ; ii. on détermine ladite au moins une fonction sigmoïde modélisant ledit au moins un signal ; iii. on déduit de ladite au moins une fonction sigmoïde une valeur représentative dudit au moins un signal stabilisé ; b) on analyse l'ensemble desdites valeurs représentatives des signaux stabilisés pour caractériser le fonctionnement dudit moteur à combustion interne. Avantageusement, on fait varier rapidement, de préférence entre 0 et 5 secondes, au moins un paramètre de réglage PR et/ou le point de puissance du moteur après le déclenchement de l'acquisition du signal. De préférence, la mesure continue est effectuée selon une fréquence d'échantillonnage comprise entre 1 Hz et 10 Hz. De manière préférée, ladite au moins une fonction sigmoïde est une fonction de Gompertz. La représentation graphique de ladite au moins une fonction sigmoïde comprend de préférence une première asymptote horizontale dont la valeur est déterminée à un instant initial te, du signal et une deuxième asymptote horizontale dont la valeur est déterminée à un instant infini t.. Selon un mode de réalisation, on modélise le signal par une fonction sigmoïde s'écrivant selon l'équation (I) suivante: y(t) = yo + ysiep x exp(-b x exp(-c x (t -0)) (I) comprenant les cinq paramètres à déterminer yo, Ystep, b, c et tb YO étant la valeur du paramètre de fonctionnement PF à un instant initial t=to; ystep étant l'écart entre la valeur de y(t) à un temps infini t=t et to; b et c étant des constantes positives, avec b caractérisant le décalage de la courbe dans le temps et c caractérisant la pente de la courbe ; ti étant étant une constante de temps positive. Selon un mode de réalisation, on modélise le signal par deux fonctions sigmoïdes. De manière préférée, les deux fonctions sigmoïdes sont combinées sous l'équation (Il) suivante: y(t)= yo + y 'epi x exp(-b1x exp(-c1x (t - t11)))± Ystep2x exP(- b2 x exp(-c2x (t - ti2))) (11) comprenant les neuf paramètres à déterminer Yo, v stepl, Ystep2, b1, c1, ti1, b2, C2e ti2 yo étant la valeur du paramètre de fonctionnement PF à un instant initial t=t0; Ystepl étant l'écart entre la valeur de y(t) à un temps infini t=t- et celle à to si un seul premier phénomène thermo-mécanique et/ou thermo-chimique existait, et représentant la variation du signal liée à ce un premier phénomène ; Ystep2 étant l'écart entre la valeur de y(t) à un temps infini t=t si seul le premier phénomène existait et la valeur de y(t) à un temps infini t=t- si un seul deuxième phénomène thermo- mécanique et/ou thermo-chimique antagoniste au premier phénomène existait, Ystep2 représentant la variation du signal liée à ce deuxième phénomène ; Ystepl et Ystep2 étant de signes opposés; !y, et cl étant des constantes positives, avec bl caractérisant le décalage de la première courbe en S dans le temps, et cl caractérisant la pente de la première courbe en S; étant une constante de temps positive relative à l'inertie du premier phénomène ; b2 et c2 étant des constantes positives, avec b2 caractérisant le décalage de la deuxième courbe en S dans le temps, et c2 caractérisant la pente de la deuxième courbe en S; ti2 étant une constante de temps positive relative à l'inertie du deuxième phénomène. L'étape ii) de détermination de ladite au moins une fonction sigmoïde modélisant le signal peut comprendre les étapes suivantes: - on sélectionne le signal ayant un critère de dynamique ou un critère d'évolution temporelle supérieur ou égal à une valeur seuil de dynamique prédéterminée ou une valeur seuil d'évolution prédéterminée (QUAL-S); - on détermine le paramètre yo correspondant à la valeur de la première asymptote horizontale (Val-AS1); - on calcule la dérivée temporelle du signal (CH-FS) ; - on définit la fonction sigmoïde selon l'équation (I) ou les deux fonctions sigmoïdes selon l'équation (II) selon que le signe de la dérivée temporelle du signal reste le même ou change (CH-FS) ; - on calcule les autres paramètres de ladite/desdites fonctions sigmoïdes de manière à minimiser au sens des moindres carrés l'écart entre le signal acquis et le signal estimé selon ladite/desdites fonctions sigmoïdes (CAL-PFs) ; et l'étape iii) de détermination de la valeur représentative du signal stabilisé peut comprendre les étapes suivantes: - on détermine la valeur de la deuxième asymptote y. (Val-AS2); - on calcule un coefficient de corrélation entre le signal acquis et la/les fonctions sigmoïdes et on sélectionne y, en tant que valeur représentative du signal stabilisé si le coefficient de corrélation est supérieur à une valeur de corrélation prédéterminée (QUAL-corr, CAL-PF). Avantageusement, on augmente la durée d'acquisition d à l'issue de l'étape iii) de détermination de la valeur représentative du signal stabilisé de manière à ce que le coefficient de corrélation soit supérieur à une valeur de corrélation prédéterminée. Le procédé peut en outre comprendre un étape de filtrage du signal (FIL) lors de l'étape ii) de détermination de la fonction sigmoïde modélisant le signal, ledit filtrage comprenant un filtre de préférence choisi dans la liste des filtres suivants: moyenne mobile, filtre de Butterworth, filtre passe-bas de premier ordre. Avantageusement, l'essai moteur est réalisé sur un banc de test moteur ou sur un véhicule en fonctionnement. Avantageusement, on calibre le moteur à combustion interne à partir de l'analyse de 35 l'ensemble des valeurs représentatives des signaux stabilisés caractérisant le fonctionnement dudit moteur. Le procédé est de préférence utilisé pour caractériser le fonctionnement d'un moteur à allumage par compression. L'invention propose également un produit programme informatique permettant de conserver dans une mémoire d'une unité d'un processeur ou sur un support de mémoire amovible approprié pour coopérer avec ladite unité du processeur, le produit-programme comprenant des instructions pour mettre en oeuvre le procédé selon l'invention. L'invention propose encore un dispositif de test d'un moteur à combustion interne permettant de réaliser une séquence d'essais moteur, comprenant: - une instrumentation de mesure pour la mesure continue d'au moins un paramètre de 10 fonctionnement PF dudit moteur générant un signal temporel brut lors d'un essai moteur ; - un système d'automatisation pour superviser la réalisation de l'essai moteur et pour afficher et/ou enregistrer les signaux acquis avec l'instrumentation de mesure, - un programme informatique selon l'invention pour déterminer une valeur représentative dudit signal stabilisé ; 15 ledit système d'automatisation pouvant afficher et/ou enregistrer l'ensemble des valeurs représentatives déterminées par le programme informatique, en vue de leur analyse pour caractériser le fonctionnement dudit moteur à combustion; et- ledit programme informatique pouvant être intégré audit système d'automatisation ou à ladite instrumentation de mesure, ou encore à un ordinateur distinct dudit système d'automatisation 20 et de ladite instrumentation de mesure, ledit ordinateur étant en communication avec ledit système d'automatisation. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures 25 annexées suivantes. La figure 1 est un diagramme illustrant schématiquement une séquence de mesure continue traditionnelle (prise de point) de paramètres de fonctionnement PF d'un moteur à combustion interne. 30 La figure 2 est un diagramme représentant l'évolution temporelle du CO2 émis dans les gaz d'échappements d'un moteur à combustion interne lors d'une séquence d'essais moteur. Chaque point du diagramme correspond à une valeur représentative de la concentration en CO2 lors d'un essai moteur (Pts) calculée selon l'art antérieur l'invention (points symbolisés par des cercles gris) et selon l'invention (points symbolisés par des carrés noirs). 35 La figure 3 correspond au signal temporel brut (quantité de CO2 mesurée pendant 220 secondes) permettant d'estimer la valeur représentative (a) de la figure 2, et fait apparaître le signal modélisé selon l'invention prenant en compte les 60 premières secondes du signal brut. Les figures 4A, 4B et 4C représentent un signal temporel brut correspondant à la mesure continue de la concentration en hydrocarbures imbrulés (HC) lors de la même séquence d'essais moteur que celle des figures 2 et 3, et représentent la modélisation de ce signal par une double courbe sigmoïde selon l'invention à partir des 60 premières secondes de signal. La figure 5 est un organigramme décrivant le calcul d'une valeur représentative d'un signal stabilisé selon l'invention. Les figures 6A et 6B sont des exemples de signaux temporels non qualifiés pour la modélisation par une fonction sigmoïde selon l'invention. Les figures 7A et 7B illustrent deux modes de réalisation de l'invention: la figure 7A illustre une modélisation du signal selon une fonction sigmoïde et la figure 7B illustre la modélisation du signal par deux fonctions sigmoïdes. Les figures 8A, 8B et 8C sont des diagrammes comparant l'évolution de la quantité de CO2 en fonction du temps calculée selon l'art antérieur et selon l'invention (Figure 8A) et représentant la performance du procédé selon l'invention (Figures 8B et 8C). Les figures 9A, 9B et 9C sont des diagrammes comparant l'évolution de la concentration en hydrocarbures imbrulés (HC) en fonction du temps, calculée selon l'art antérieur et selon l'invention (Figure 9A), et représentant la performance du procédé selon l'invention (Figures 9B et 9C). La figure 10 est un schéma illustrant la mise en oeuvre du procédé selon l'invention dans un banc moteur. On entend par paramètre(s) de fonctionnement du moteur, notés PF dans la présente description, des paramètres pouvant correspondre à une grandeur physique ou chimique permettant d'obtenir des informations sur le fonctionnement du moteur et reflétant ainsi un état de fonctionnement du moteur ou de ses organes. Ce peut être une concentration d'un composé dans les gaz d'échappement, comme le monoxyde de carbone (CO), le dioxyde de carbone (CO2), les hydrocarbures imbrûlés (HC), les différents oxydes d'azote (N0x), la température ou la pression des gaz à l'admission ou à l'échappement, la pression maximale des gaz dans la chambre de combution, un angle du vilebrequin auquel se situe une fraction de masse de carburant brûlée, un bruit de combustion calculé d'après un signal de pression du cylindre, un débit de carburant, une vitesse de rotation d'un arbre de turbo-compresseur, la richesse des gaz à l'admission ou à l'échappement, etc. On entend par paramètre(s) de réglage du moteur, noté(s) PR dans la présente description, des paramètres tels qu'un débit d'air, une pression des gaz à l'admission, une quantité de carburant injectée pour chaque injection, une pression d'injection de carburant, un temps de séparation entre différentes injections, un niveau de swirl dans la chambre de combustion, ou tout autre paramètre permettant de modifier le déroulement de la combustion et par conséquence le comportement du moteur. Selon l'invention, il est possible de caractériser le fonctionnement d'un moteur à combustion interne, en particulier un moteur à allumage par compression, à partir d'une séquence d'essais moteur, chaque essai comprenant la modifications d'au moins un paramètre de réglage du moteur PR et/ou de son point de puissance. Pour cela, le demandeur propose un procédé comprenant les étapes suivantes: a) pour chaque essai moteur: I. on acquiert au moins un signal issu de la mesure continue d'un paramètre de fonctionnement PF du moteur sur une durée d'acquisition d, la durée d'acquisition d étant choisie de manière à pouvoir déterminer au moins une fonction sigmoïde pour modéliser ledit au moins un signal ; ii. on détermine ladite au moins une fonction sigmoïde modélisant ledit au moins un signal ; iii. on déduit de ladite au moins une fonction sigmoïde une valeur représentative dudit au moins un signal stabilisé ; b) on analyse l'ensemble des valeurs représentatives des signaux stabilisés pour caractériser le fonctionnement dudit moteur à combustion interne. L'ensemble des valeurs représentatives des signaux stabilisés issus d'une séquence d'essais moteur comprenant des mesures continues d'un ou plusieurs paramètres de fonctionnement PF du moteur permet d'établir l'évolution temporelle du ou des paramètres lors de l'essai. C'est l'analyse de cette évolution du ou des paramètres durant toute la durée de la séquence des essais moteur (somme des durées des essais moteur) qui permet de caractériser le fonctionnement dudit moteur à combustion. La figure 2 illustre une telle évolution temporelle pour la concentration en CO2 dans les gaz d'échappement lors d'une séquence d'essais moteur. L'évolution de la concentration en en CO2 est représentée par la valeur de la concentration en CO2 (en pourcentage) en fonction des essais moteurs indicés en abscisse, aussi appelés "points" (points 9330 à 9520). Chaque point de mesure, symbolisé par un carré, correspond à la valeur représentative d'un signal stabilisé calculée selon l'invention, c'est-à-dire à la valeur estimée stable du signal modélisé par au moins une fonction sigmoïde selon l'invention. Les cercles gris symbolisent les valeurs calculées selon une prise de point traditionnelle, où la valeur représentative correspond a une moyenne réalisée sur une phase de moyennage du signal. Ainsi, pour un nombre donné n de signaux acquis relativement à un paramètre de fonctionnement du moteur PR lors d'une séquence de n essais moteur, n étant suffisant à caractériser le fonctionnement du moteur, on calcule n valeurs représentatives (ou points de mesure) permettant alors de déterminer l'évolution du paramètre de fonctionnement PR en fonction du temps, chaque signal acquis, aussi communément appelés "point", ayant une durée d'acquisition d. Sur la figure 2, l'évolution du CO2 dans le temps est réalisée grâce à 191 signaux acquis lors de 191 essais moteurs (indice 9520 - indice 9330), représentés en 10 abscisses sous la nomination "points", donnant 191 points de mesure. La séquence d'essais moteur illustrée à la figure 2 comprend les caractéristiques suivantes: la séquence d'essais enchaîne 191 variations rapides et simultanées du point de puissance et des réglages en plan d'expériences. La durée totale de la séquence d'essais est de 9h45, soit environ 184 s par essai moteur unitaire. Le moteur testé est un moteur 15 diesel (PSA DV6 Turbo-compresseur Diesel Injection directe Common-Rail). Les essais sont réalisés sur un banc moteur dynamique. L'acquisition de tous les signaux est réalisée à une fréquence de 1Hz. Les paramètres de fonctionnement suivants sont mesurés: - polluants : CO2, CO, HC, NOx, 02, CO2 recirculé à l'admission; - tension délivrée par une sonde de richesse; 20 - calculs instantanés : richesse "5gaz" (richesse calculée à partir des concentrations de cinq gaz), bruit de combustion (système Noisemeter), angle de la fraction de gaz brûlés à 50%, maximum de pression cylindre; cinétique : vitesse de rotation du turbo-compresseur; températures : gaz à l'admission, gaz à l'échappement; 25 pressions : gaz à l'admission, gaz à l'échappement; débits : air frais à l'admission, carburant (via un système KMA4000). La représentation graphique de la fonction ou des fonctions sigmoïdes modélisant le signal comprend de préférence une première asymptote horizontale dont la valeur est 30 déterminée à un instant initial 1.0 du signal et une deuxième asymptote horizontale dont la valeur est représentative d'un instant qui serait situé à l'infini to 0. L'instant initial to correspond aux premiers échantillons du signal temporel, lorsque les paramètres de réglages PR et/ou le point de puissance du moteur n'ont pas encore été modifiés. L'instant infini L correspond aux échantillons du signal stabilisé. Ainsi, la deuxième asymptote correspond à la partie 35 stabilisée du signal, c'est-à-dire la partie du signal où les variations ne sont plus significatives, c'est-à-dire qu'elles sont plus faibles qu'un seuil propre à chaque paramètre, par exemple quelques pourcents, pendant une durée donnée, par exemple pendant environ 10 secondes. Selon l'invention, le signal temporel est modélisé par au moins une fonction sigmoïde. De manière préférée, le signal est modélisé par au moins une fonction sigmoïde qui est une fonction de Gompertz. Selon un mode de réalisation, le signal est modélisé par une fonction sigmoïde s'écrivant selon l'équation (I) suivante, qui est une fonction de Gompertz : y(t) = yc, + Ystep x exp(-b x exp(-c x (t - t1))) (I) avec: yo étant la valeur du paramètre de fonctionnement PF à un instant initial t=to; Ystep étant l'écart entre la valeur de y(t) à un temps infini t=t«, et celle à un temps initial to; b et c étant des constantes positives, avec b caractérisant le décalage de la courbe dans le temps, et c caractérisant la pente de la courbe, et; t étant une constante de temps positive, et traduisant que l'inertie du phénomène physique (impact de la modification d'un paramètre de réglage du moteur et/ou de son point de puissance sur un paramètre de fonctionnement du moteur) n'est pas nécessairement liée à to. Dans ce cas, la fonction comprend cinq paramètres à déterminer yo, Ystep, b, c et ti; La figure 3 illustre la représentation graphique d'une fonction sigmoïde selon l'équation (I), modélisant le signal issu de la mesure de la concentration en CO2 dans les gaz d'échappement au cours d'un essai moteur (indicé Pt 9382) de la séquence d'essais moteur de la figure 2. La figure 3 fait apparaître le signal temporel brut acquis sur 220 secondes, et la courbe en S modélisant un tel signal déterminée selon l'invention, permettant de fournir la valeur représentative (a) symbolisé par un carré dans la figure 2. Une durée de 60 secondes suffit à établir la fonction sigmoïde modélisant le signal acquis. Selon un autre mode de réalisation, le signal est modélisé par deux fonctions sigmoïdes, dont la représentation graphique peut comprendre deux courbes sigmoïdes de sens de variation opposés. Une telle double courbe en S avec des sens de variation opposés est illustrée par la figure 3 et les figures 4A, 4B et 4C. La courbe FS de la figure 4A se superpose à un signal représentant la mesure continue dans le temps de la concentration en HC lors d'un essai moteur (essai moteur indicé par le point 9501) de la séquence d'essais décrite ci-dessus en relation avec les figures 2 et 3. De préférence, les deux fonctions sigmoïdes sont des fonctions de Gompertz et s'expriment sous l'équation (Il) suivante: y(t) =- yo + y stepi x exp(-b1x exp(-c1x (t - ti1)))+ Ystep2 x exp(-b2x exp(-c2x (t - t12))) (Il) avec: yo étant la valeur du paramètre de fonctionnement PF à un instant initial t=to ; ystepi étant l'écart entre la valeur de y(t) à un temps infini t=t., et celle à un temps initial to si un seul premier phénomène thermo-mécanique et/ou thermo-chimique existait et représentant la variation du signal liée à ce premier phénomène ; Ystep2 étant l'écart entre la valeur de y(t) à un temps infini t=t., si seul le premier phénomène existait et la valeur de y(t) à un temps infini t=tu si un seul deuxième phénomène thermo- mécanique et/ou thermo-chimique antagoniste au premier phénomène existait, Ystep2 représentant la variation du signal liée à ce deuxième phénomène ; Ystepl et Ystep2 étant de signes opposés ; 1)1 et c, étant des constantes positives, avec bl caractérisant le décalage de la première courbe en S dans le temps, et c1 caractérisant la pente de la première courbe en S; t11 étant une constante de temps positive relative à l'inertie du premier phénomène ; b2 et c2 étant des constantes positives, avec b2 caractérisant le décalage de la deuxième courbe en S dans le temps, et c2 caractérisant la pente de la deuxième courbe en S; t12 étant une constante de temps positive relative à l'inertie du deuxième phénomène ; Cette deuxième fonction comprend neuf paramètres à déterminer v , , v stepl, Ystep2, 131, Cl, b2, c, ti2 ; La modélisation du signal par une telle double fonction sigmoïde avec une double courbe en S, les deux courbes ayant des sens de variation opposés, est particulièrement bien adaptée au cas où il existe une compétition entre plusieurs phénomènes thermomécaniques et/ou thermo-chimiques avec une inertie et/ou une dynamique différentes et antagonistes. Cela peut par exemple être le cas des mesures continues de la concentrations en HC, telle que représenté dans les diagrammes des figures 4A, 4B et 4C. Les paramètres à déterminer de ces fonctions permettent d'obtenir une courbe ou plusieurs courbes en S qui s'adaptent bien à la forme des signaux temporels enregistrés lors d'essais sur un moteur à combustion interne. D'autres fonctions sigmoïdes aptes à modéliser le signal peuvent être utilisées dans le procédé selon l'invention, telle que la fonction sigmoïde s'exprimant de la manière suivante: y(t)= yo + ystep x exp(-b xi) , où yo est la valeur du paramètre de fonctionnement PF à un instant initial t=to, ystep est l'écart entre la valeur de y(t) à un temps infini t=t., et celle à un temps initial to, et b est une constante positive caractérisant le décalage de la courbe dans le temps. La figure 5 présente schématiquement ce qui est réalisé pour un signal donné relatif à un paramètre de fonctionnement PF à chaque essai moteur, afin de déterminer une valeur représentative du signal stabilisé relatif à un paramètre de fonctionnement donné, acquis sur une durée d'acquisition d. Comme l'illustre la figure 5, huit étapes principales sont réalisées: - Étape 1: Acquisition d'un signal temporel brut à partir d'une mesure continue (ACQS). - Étape 2: Qualification du signal temporel brut (QUAL-S). - Étape 3: Calcul de la valeur yo (valeur de la première asymptote) (Val-AS1). - Étape 4: Étape optionnelle de filtrage du signal brut temporel (FIL). - Étape 5: Choix de la(des) fonction(s) sigmoïde(s) et du nombre de paramètres associés (CH-FS). - Étape 6: Détermination des paramètres de la(des) fonction(s) sigmoïde(s) (CAL- PFs). - Étape 7: Calcul de la valeur y. (valeur de la deuxième asymptote) (Val-AS2). - Étape 8: Qualification du résultat avec le calcul d'un coefficient de corrélation entre le signal et la(les) fonction(s) sigmoïde(s) (QUAL-corr), permettant de déterminer (CAL- PF) la valeur représentative du signal stabilisé (point de mesure). Les étapes 2 à 6 permettent de déterminer la (les) fonction(s) sigmoïde(s) modélisant le signal (brut ou filtré), et les étapes 7 et 8 permettent la détermination d'une valeur représentative dudit signal stabilisé à partir de ladite modélisation par la(les) fonction(s) sigmoïde(s) effectuée dans les étapes 2 à 6.30 Étape 1 : Acquisition d'un signal temporel brut à partir d'une mesure continue Lors de l'essai moteur, on fait de préférence varier rapidement, par exemple entre 0 et 5 secondes, au moins un paramètre de réglage PR du moteur et/ou son point de puissance. Le point de puissance du moteur peut notamment être caractérisé par le régime du moteur ou la charge du moteur. La charge du moteur peut être caractérisée par différentes grandeurs, comme par exemple les quantités d'air introduit et de carburant injecté, le couple moyen effectif (ou indiqué si on prend en compte les frottements) ou la pression moyenne effective (ou indiquée). Cette dernière grandeur correspond à l'effort constant qui s'exercerait sur chaque piston au cours d'un cycle complet et qui peut être reliée au couple par la 10 cylindrée du moteur. Ces variations sont assimilables à un échelon quand elles sont instantanées et à une rampe quand elles durent quelques secondes. Cela a pour effet de changer les conditions initiales du mélange introduit dans la chambre de combustion ainsi que le déroulement de cette réaction et en conséquence tout ou partie des paramètres de fonctionnement du 15 moteur PF. Lorsqu'un paramètre de fonctionnement PF évolue de manière significative, on observe que l'enregistrement temporel du signal issu de sa mesure présente le comportement caractéristique suivant: un palier précède la variation qui commence lentement puis s'accélère puis ralentit avant de s'achever par un autre palier ; entre l'accélération et le 20 ralentissement, on distingue un point d'inflexion. Ce comportement est le plus souvent à l'ceuvre même quand le bruit de mesure masque sa présence. La variation brutale du paramètre de réglage PR du moteur et/ou son point de puissance est telle que le signal acquis est identifiable par au moins une fonction sigmoïde. Une étape d'acquisition dans le temps de signaux bruts est réalisée préalablement à 25- toute étape utilisant de tels signaux tel que décrit ci-après. Cette étape permet d'enregistrer les effets de la variation d'au moins un paramètre de réglage PR du moteur et/ou de son point de puissance sur les paramètres de fonctionnement PF du moteur. Cette étape d'acquisition consiste à récupérer au moins un signal issu de la mesure continue d'un paramètre de fonctionnement PF du moteur à combustion interne, cette 30 mesure étant effectuée à l'aide de moyens de mesure adaptés, et à enregistrer ce signal dans une unité d'acquisition. On entend, par mesure continue lors d'un essai moteur, toute mesure effectuée à intervalles de temps réguliers et courts, de préférence effectuée à une fréquence d'échantillonnage f comprise entre 1 Hertz et 10 Hertz. Cette gamme de fréquence est 35 adaptée à la durée de variation des paramètres de fonctionnement PF qui s'établit de préférence sur quelques dizaines de secondes. Plus la fréquence est élevée et plus le signal offre une réserve de filtrage pour en éliminer le bruit. Ces mesures peuvent être des mesures permanentes de pression et température des gaz, de vitesse de rotation d'un arbre de turbo-compresseur, de tension de sonde de richesse. Certaines mesures permanentes, telles que les mesures rapides de pression cylindre, peuvent également permettre de calculer des grandeurs décrivant la combustion, par exemple un angle du vilebrequin pour lequel 50% de la fraction de carburant a brûlé pour chaque cylindre, une pression maximale instantanée dans chaque cylindre, un bruit de combustion dans chaque cylindre. Certaines mesures peuvent être interrompues par des opérations de purge ou de remplissage, comme la concentration en espèces chimiques des gaz de combustion, l'opacité des gaz d'échappement ou le nombre de particules qu'ils contiennent, la mesure instantanée de la consommation de carburant. Ces mesures sont effectuées par des dispositifs synchrones présentant un taux de disponibilité élevé. Chaque mesure continue est effectuée sur une durée d'acquisition donnée d, cette durée étant choisie de manière à pouvoir déterminer au moins une fonction sigmoïde pour modéliser le signal. Cette durée peut être déterminée empiriquement, sur la base d'essais moteurs antérieurs, et/ou être calculée lors du procédé, de manière à ce que le signal puisse être modélisé par une fonction sigmoïde. Un mode de réévaluation de la durée d lors du procédé est détaillé plus loin dans la description, en relation avec l'étape 8. Dans la figure 8A, chaque point de mesure calculé selon la présente invention, et symbolisé par un carré noir, est par exemple issu de la modélisation du signal sur les 60 premières secondes d'acquisition avec une unique fonction sigmoïde selon l'équation (I). Une durée d'acquisition d de 60 secondes suffit donc dans cet exemple à permettre la modélisation du signal par une fonction sigmoïde selon l'invention. La durée d'acquisition d peut être fixe ou variable d'une acquisition de signal à un autre, lors d'un essai moteur ou d'un essai moteur à un autre. Au cours d'un même essai moteur, lors duquel l'acquisition de plusieurs signaux peut être réalisée, une durée d'acquisition d variable permet de s'adapter aux différentes grandeurs mesurées, certaines pouvant être plus réactives que d'autres. C'est le cas par exemple pour la richesse des gaz obtenue par l'intermédiaire d'une sonde à oxygène située à l'échappement, ou le bruit de combustion qui est un calculé à partir de la mesure réalisée par un capteur de pression cylindre placé directement dans la chambre de combustion. Ces paramètres sont issus de mesures réalisées grâce à des capteurs situés très près de la source qu'ils mesurent, et leur temps de réponse est élevé, de l'ordre de la seconde à quelques secondes. Dans le cas où la durée d est fixe pour un même essai moteur, on détermine de préférence cette durée d sur la base du paramètre mesuré qui nécessite la plus longue durée d'acquisition. Par exemple, ce paramètre de fonctionnement PF peut être un polluant tel que les hydrocarbures imbrûlés HC qui présentent généralement l'évolution du signal la plus lente. Étape 2: Modélisation par au moins une fonction sigmoïde: qualification du sianal temporel brut Cette étape a pour but de déterminer si le signal acquis peut être traité au moyen d'une modélisation par au moins une fonction sigmoïde afin de calculer une valeur représentative du signal stabilisé. Lors de cette étape, on peut calculer deux critères permettant de qualifier ou non le signal : - un premier critère evol d'évolution du signal sur la durée d, exprimant l'évolution du signal dans le temps. Le premier critère évol s'écrit selon la formule suivante: evol = 1 moyenne(signal brut sur les premiers échantillons) moyenne(signal brut sur la durée d) - un deuxième critère dyn sur la dynamique du signal brut par rapport à sa variabilité intrinsèque, ce qui revient à évaluer un rapport amplitude du signal sur bruit de mesure. Ce deuxième critère s'écrit de la manière suivante: dyn = max(signal filtré) - min(signal filtré) écart - type(signal brut - signal filtré) Le filtrage du signal effectué pour le calcul de ce critère peut être différent de celui éventuellement réalisé à l'étape 4 décrite plus bas. Avantageusement, cette formule permet d'éviter de définir une valeur de bruit pour chaque mesure. Si le premier critère evol est inférieur à une valeur seuil AA prédéterminée ou si le deuxième critère de dynamique dyn est inférieur à une valeur seuil prédéterminée de dynamique A, traduisant le cas d'un signal qui varie peu ou qui est très bruité, alors on estime que le signal n'est pas assimilable à une courbe sigmoïde, et on ne poursuit pas avec les étapes suivantes 3 à 8. Dans un tel cas, on détermine de préférence une valeur représentative du signal qui n'évolue pas, et considéré par conséquent comme stabilisé, comme étant une valeur alternative ALT-PF représentant le signal à un temps infini, par exemple la valeur médiane de tous les échantillons de l'acquisition. Les figures 6A et 6B sont des exemples de signaux non qualifiés sur la base respectivement du critère d'évolution du signal evol (signal évoluant trop peu) et du critère de dynamique du signal dyn (signal trop bruité). Si les premier et deuxième critères sont satisfaits (supérieurs respectivement aux valeurs seuils AA et A), alors le signal est qualifié, i.e. il contient de l'information pouvant être traitée, et les étapes suivantes du procédé peuvent être réalisées. Étaie 3: Modélisation ar au moins une fonction si moïde: calcul de valeur valeur de la première asymptote) 10 Lors de cette troisième étape, on utilise les premiers échantillons temporels du signal pour déterminer définitivement la valeur yo de la première asymptote. La valeur yo correspond à la valeur stable enregistrée avant la modification des paramètres de réglages PR et/ou le point de puissance du moteur. Ce paramètre peut être calculé quelque soit la(les) fonction(s) sigmoïde(s) utilisée(s) pour modéliser le signal car il fait nécessairement partie de 15 l'équation de la ou des fonctions sigmoïdes modélisant le signal. Il est par exemple commun aux équations (I) et (II). Étape 4 (optionnelle): Modélisation par au moins une fonction sigmoïde: filtrage du signal 20 Une étape optionnelle de filtrage du signal brut acquis, représentée en pointillés dans la figure 5, peut être réalisée. Certains signaux peuvent en effet présenter un bruit de mesure important, qui peut être atténué par lissage du signal en appliquant un filtre, lequel est choisi en fonction de sa fréquence d'échantillonnage. 25 De préférence, le filtre a un effet de déphasage minimum, c'est-à-dire qu'il produit un décalage temporel minimum entre le signal filtré et le signal brut. Parmi les filtres pouvant être appliqués, on peut citer de manière non limitative les filtres suivants bien connus de l'homme du métier, qui peuvent être choisis en prenant en compte leur degré de simplicité d'implémentation dans un programme d'ordinateur et le 30 temps de calcul requis: - une moyenne mobile, qui a l'avantage d'être simple d'implémentation et peu coûteuse en temps de calcul. La moyenne mobile est de préférence centrée afin de ne pas produire de déphasage du signal ; 3 0002 06 18 - un filtrage passe-bas de I er ordre, qui a l'avantage d'être peu coûteux en temps de calcul. L'utilisation d'un tel filtre requiert de préférence la compensation du déphasage créé ; - un filtre de Butterworth. L'utilisation d'un tel filtre requiert de préférence la compensation du déphasage créé ; Étape 5: Modélisation par une au moins fonction sigmoïde: choix de la(des) fonction(s) siomoïde(s) Cette cinquième étape a pour objectif de déterminer la(les) fonction(s) sigmoïde(s) utilisée(s) pour modéliser le signal. Pour cela, le choix entre la modélisation par une fonction sigmoïde selon l'équation (I) ou par deux fonctions sigmoïdes selon l'équation (II) est basé sur le calcul de la dérivée temporelle du signal, brut ou filtré. Selon si la dérivée temporelle conserve ou change de signe, on choisit respectivement la modélisation selon l'équation (I) ou (II). Si la modélisation par une fonction sigmoïde selon l'équation (I) est choisie, les inconnues à déterminer sont constituées par les quatre paramètres ystep, b, c et t. Si la modélisation par deux fonctions sigmoïdes selon l'équation (II) est choisie, ce sont les huit paramètres v stepl, Ystep2, b1, C1, 1.11, b2, C2, ti2 qui doivent être déterminés. Étaie 6: Modélisation ar une fonction si moïde: calcul des aramètres P s de la des fonction(s) sipmoïde(s) choisie(s) Cette étape vise à déterminer tous les paramètres non encore déterminés de la(des) fonction(s) sigmoïde(s) choisie(s) à l'étape précédente. Pour cela, on réalise une régression statistique mettant en oeuvre une méthode des moindres carrés pour déterminer les valeurs des quatre paramètres ou des huit paramètres de la(des) fonction(s) sigmoïde(s), le paramètre yo ayant déjà été déterminé à l'étape 3, à partir du signal brut ou filtré. La minimisation au sens des moindres carrés effectuée vise à minimiser un critère qui est égal à la somme quadratique des écarts entre les valeurs du signal acquis (brut ou filtré) et les valeurs du signal modélisé par la (les) fonction(s) sigmoïde(s): E (échantillon du signal brut ou filtré - échantillon du signal identifie échantillons Lorsque les valeurs des quatre (ou huit) paramètres PFs de la (des) fonction(s) sigmoïde(s) évoluent en faisant converger le critère vers un minimum, on obtient alors des valeurs optimales pour ces quatre (ou huit) paramètres PFs qui permettent une estimation du signal acquis. Dans le cas où le processus de minimisation ne converge pas vers une solution (N conv dans la figure 5), on peut augmenter la durée d'acquisition d et reprendre à l'étape 1 et/ou changer de filtre et reprendre à l'étape 4, jusqu'à obtenir une convergence, ou effectuer une prise de point traditionnelle. Étape 7: Détermination d'une valeur représentative du signal stabilisé: calcul de la valet valeur de la deuxième asymptote) Cette étape permet de déterminer la valeur représentative du signal stabilisé pour le paramètre de fonctionnement PF . Pour cela, on calcule la valeur de la deuxième asymptote de la représentation graphique de la(des) fonction(s) sigmoïde(s), y. Si le signal est modélisé par au moins une fonction sigmoïde selon les équations (I) ou (Il), la valeur de la deuxième asymptote est calculée respectivement selon les équations (III) et (IV) ci-dessous: Y = Y 0 ± Y step (III) Yoe = YO + Y stepl Y step2 (IV) Les figures 7A et 7B illustrent le calcul de la valeur de la deuxième asymptote selon les équations (III) et (IV). Connaissant les quatre ou huit paramètres de la fonction sigmoïde, on connait alors la valeur de y. qui représente la valeur du signal temporel stabilisé. Étape 8: Détermination d'une valeur représentative du signal stabilisé: qualification du résultat A partir des valeurs obtenues pour les paramètres de la(des) fonction(s) sigmoïde(s) PFS, on construit la représentation graphique de la(des) fonction(s) sigmoïde(s) modélisant le signal, et on calcule un coefficient de corrélation C', entre ce signal estimé, et le signal acquis, brut ou filtré. Le coefficient de corrélation peut s'écrire selon la formule suivante: c = cor \IEni.1(x 57)2 liEni.1(Yi _ 17)2 Le coefficient de corrélation varie entre -1 (signaux anti-corrélés) et 1 (signaux parfaitement corrélés). Lorsqu'il prend la valeur nulle, les signaux ne sont pas corrélés. De préférence, le coefficient de corrélation est calculé sur toute la durée d'acquisition du signal (n représentant alors tous les échantillons du signal acquis sur la durée d). Alternativement, le coefficient de corrélation est calculée sur une fenêtre temporelle située vers la fin du signal (n représentant alors le nombre d'échantillons du signal sur cette fenêtre temporelle). Si ce coefficient de corrélation est élevé, par exemple supérieur ou égal à une valeur de corrélation prédéterminée B (par exemple 8 = 0,98), on considère que le résultat est fiable, et on sélectionne la valeur de y. (deuxième asymptote) en tant que valeur représentative du signal stabilisé CAL-PF. On peut passer à l'étape 1 si un autre signal est acquis lors de l'essai moteur. Dans le cas contraire (Ccor Avantageusement, on augmente la durée d'acquisition d du signal à l'issue de cette étape de qualification du résultat, dans le cas où le résultat n'est pas qualifié en raison d'un coefficient de corrélation inférieur à la valeur de corrélation prédéterminée 8, jusqu'à ce que le coefficient de corrélation soit supérieur à la valeur B. Pour cela, l'ordre de poursuivre la mesure est par exemple envoyé à l'instrumentation de mesure qui fournit la mesure continue du paramètre de fonctionnement, et on répète au moins les étapes 1, et 4 à 8, un nombre de fois suffisant pour atteindre un coefficient de corrélation supérieur à la valeur B. Il peut être fixé un nombre de fois maximum n max pour lequel on augmente la durée d à l'issue de cette huitième étape, par exemple n max égal à 2 fois. L'augmentation de durée d peut être fixée par défaut à une valeur donnée d_plus, par exemple d_plus égal 30 secondes. On s'assurera que la durée maximum d + (n_max * d_plus) est inférieure à la durée moyenne d'une prise de point classique afin de ne pas à dégrader la durée globale de la séquence d'essais par rapport à une approche traditionnelle. Le temps de calcul des étapes 2 à 8 assurant le traitement de l'ensemble des paramètres de fonctionnement PF est suffisamment court pour permettre de réaliser cette augmentation optionnelle de la durée d'acquisition d du signal lors du procédé. Ce temps de calcul est de préférence inférieur à la valeur de d_plus afin de permettre la réalisation de l'augmentation optionnelle de d, sous réserve que l'unité d'acquisition du signal est apte à poursuivre son traitement de manière indépendante de l'unité de calcul qui exécute les étapes 2 à 8. Cela est le cas d'un système d'automatisation qui réunit au moins ces deux unités dans un même ensemble, et tel que décrit à la figure 10. Une fois déterminées toutes les valeurs représentatives des signaux stabilisés acquis lors de la séquence d'essais moteur, concernant au moins un paramètre de fonctionnement du moteur, on analyse l'ensemble de ces valeurs pour caractériser le fonctionnement du moteur. L'évolution dans le temps de chaque paramètre de fonctionnement, donnée par l'ensemble des valeurs représentatives lors de la séquence d'essais moteur, est enregistrée et/ou affichée dans un système comprenant une unité d'affichage, une mémoire, et éventuellement une unité de calcul. L'analyse peut être réalisée automatiquement par l'unité de calcul, qui compare par exemple l'évolution temporelle des paramètres de fonctionnement du moteur testé avec des valeurs types, par exemple dans le but d'assurer la sécurité de fonctionnement du moteur ou de ses organes en respectant son cahier des charges d'utilisation, ou par un utilisateur du procédé, qui peut formuler des conclusions sur le fonctionnement du moteur, par exemple dans le but d'indiquer le niveau de polluants émis, à partir des valeurs représentatives déterminées et de leur évolution dans le temps. Certaines séquences d'essais moteur visent par exemple à régler le fonctionnement du moteur en pleine charge. Au cours d'une séquence de ce type, on place le moteur à un régime donné et à sa puissance maximale telle que renseignée au cahier des charges, et on fait varier les réglages du moteur par saut à partir d'un réglage de départ. Ces sauts produisent une évolution des paramètres de fonctionnement PF du moteur, dont certains sont liés à la sécurité d'exploitation et sont de préférence analysés. Ces paramètres sont par exemple la pression des gaz à l'intérieur du cylindre, qui ne doit pas dépasser une certaine valeur au risque de provoquer un flambement d'une bièle, ou encore la température d'échappement, qui ne doit pas dépasser une certaine valeur au risque de provoquer la fusion des aubes de la turbine pour un moteur suralimenté. D'autres séquences d'essai visent par exemple à régler le fonctionnement du moteur en zone dépollution, c'est-à-dire une zone correspondant à une plage en régime /couple du moteur dans laquelle on maîtrise les émissions de polluants pour répondre aux contraintes réglementaires. Au cours d'une séquence de ce type, on place le moteur à un régime et une puissance fixés représentatifs des conditions d'utilisation que le véhicule est susceptible de rencontrer en suivant le profil de conduite d'un cycle d'homologation des émissions polluantes. A partir d'un réglage de départ, on fait varier les réglages par saut ; ces sauts vont produire une évolution des paramètres de fonctionnement PF du moteur et ceux qu'il convient d'analyser sont liés aux émissions polluantes : par exemple la concentration de CO2 dans les gaz d'échappement, qui constitue un polluant global ayant un impact sur le climat que l'on cherche à limiter pour répondre aux normes en vigueur, ou bien la concentration des hydrocarbures imbrûlés (famille des HC) dans ces mêmes gaz d'échappement, qui constituent un polluant local ayant un impact sur la qualité de l'air. Les hydrocarbures imbrûlés doivent être également limités afin de respecter les normes en vigueur (ex: normes Euro5 sur les émissions polluantes s'appliquant sur les véhicules neufs et jusqu'à ce qu'ils atteignent 100000km en conditions d'utilisation normales). La figure 8A montre l'évolution d'un paramètre de fonctionnement d'un moteur à combustion interne, la concentration en CO2 dans les gaz d'échappement, au cours d'une séquence d'essais moteur. Lors de cette séquence d'essais, 189 variations rapides et simultanées des réglages en plan d'expériences sur un point de puissance fixe à (1750tr/min, 2bar) sont enchaînées. La durée totale de la séquence d'essais est de 7h32, soit environ 143 s par essai moteur unitaire. Cette séquence d'essais est réalisée sur un banc moteur à frein. Le moteur testé est un moteur diesel (PSA DV6 Turbo-compresseur Diesel Injection directe Common-Rail), alimenté avec un carburant standard. La fréquence d'acquisition de tous les signaux réalisés est de 1Hz. Les paramètres de fonctionnement mesurés sont les suivants: - polluants : CO2, CO, HC, NOx, 02, CO2 recirculé à l'admission; - tension délivrée par une sonde de richesse; - calculs instantanés : richesse "5gaz"; - températures : gaz à l'admission, gaz à l'échappement; - pressions : gaz à l'admission, gaz à l'échappement; - débits : air frais à l'admission. Dans la figure 8A, les points symbolisés par des carrés noirs sont les valeurs représentatives de 189 signaux (189 points) acquis lors de la séquence d'essais moteur, calculées selon le procédé de l'invention, avec une modélisation de chaque signal par une seule fonction sigmoïde selon l'équation (I) et une durée d'acquisition d de 60s. A titre comparatif, le diagramme de la figure 8A comprend des points symbolisés par des cercles gris représentant les valeurs calculées selon une prise de point traditionnelle, comprenant une phase de moyennage de 30 secondes servant à déterminer les points de mesure. On observe clairement une superposition des deux types de valeurs, ce qui traduit la bonne performance du procédé selon l'invention permettant d'obtenir en un temps réduit des valeurs représentatives des signaux stabilisés équivalentes à celles calculées de manière traditionnelles en un temps réduit. La bonne cohérence entre les valeurs calculées selon l'invention et selon l'art antérieur est également nettement visible sur le diagramme de la figure 8B, qui représente le résidu, c'est-à-dire la différence entre la valeur estimée selon l'invention et la valeur moyenne calculée selon l'art antérieur, entre les deux types de valeurs en fonction de la quantité de CO2 mesurée. L'ensemble des résidus a une forme de tuyère divergente. La figure 8C est un diagramme de parité, où les valeurs de CO2 estimées selon l'invention (CO2cAL en ordonnée) sont représentées en fonction des valeurs calculées selon l'art antérieur (CO2TRAD en abscisse), est une autre manière de représenter les résidus: on voit que le résidu est faible, les points étant proches sur la droite oblique y=x. Les figures 9A, 9B et 9C représentent, de la même manière que pour les figures 8A à 8C, l'évolution de la quantité de HC lors de la même séquence d'essais moteur, déterminée selon l'invention et selon l'art antérieur (figure 9A), ainsi que les performances du procédé selon l'invention (figure 9B et 9C). On observe, à l'instar de la quantité de 002, une très bonne correspondance entre les valeurs calculées selon une prise de point traditionnelle (crecles gris) et selon l'invention (carrés noirs). Les résidus représentés à la figure 9B sont faibles, bien qu'un peu plus dispersés que pour le 002, et reflètent une très bonne performance du procédé pour une mesure de HC qui est généralement plus difficile à réaliser que celle du 002. La figure 9C montrent une très bonne correspondance, et donc un faible résidu, entre les valeurs de HC estimées selon l'invention et les valeurs de HC calculées selon une prise de point traditionnelle. Le procédé selon l'invention permet un gain de temps significatif lors de l'essai moteur comparativement aux prises de points effectuées de manière classique. Ce gain de temps est estimé à environ 50 % du temps traditionnellement requis pour effectuer l'essai moteur. Dans les exemples présentés aux figures 2, 3 et 4A, 4B, 4C, la durée de la séquence d'essais moteur est réduite de 5 heures et 15 minutes, soit près de 54 %, en considérant une durée d'acquisition fixe de 60 secondes pour chaque point suivie de 25 secondes en moyenne permettant d'obtenir les valeurs des systèmes de mesure asynchrones, par rapport à une séquence d'essais traditionnelle qui aurait requis 9h45 avec une durée moyenne d'acquisition de 184 secondes pour chaque point. Dans les exemples présentés aux figures 8A, 8B, 8C, 9A, 9B et 9C, la durée de la séquence d'essais moteur est réduite de 3 heures, soit près de 41 %, en considérant une durée d'acquisition fixe de 60 secondes pour chaque point suivie de 25 secondes en moyenne permettant d'obtenir les valeurs des systèmes de mesure asynchrones, par rapport à une séquence d'essai traditionnelle qui aurait requis 7h32 avec une durée moyenne d'acquisition de 143 secondes pour chaque point. Le procédé selon l'invention permet de prédire la valeur représentative du signal stabilisé, ce qui autorise à anticiper le passage à la mesure continue suivante lors de l'essai sans avoir à attendre la stabilisation complète du signal. Il est alors fait l'économie du délai d'attente lié au moins à la phase de moyennage traditionnelle de la prise de point. Une telle économie de temps permet non seulement une réduction significative de la durée totale de l'essai, mais permet également d'éviter au moteur d'atteindre de manière prolongée une limite thermomécanique, et en cela offre une meilleur protection du moteur et de ses organes. L'invention permet donc de réduire la durée de l'essai moteur tout en améliorant la sécurité d'exploitation du moteur. Le gain de temps réalisé peut contribuer à accélérer les cadences d'essais sur bancs moteurs ou permettre d'obtenir plus rapidement des informations sur le fonctionnement du moteur dans un véhicule en fonctionnement, mais aussi contribuer à limiter l'usure des moteurs et/ou réduire ou alléger les opérations de maintenance sur les moteurs (vidange/renouvellement de l'huile moteur toutes les 100 heures) et les équipements de mesure (intervention bisannuelle pour la baie d'analyse des gaz d'échappement). Cela a également un autre effet vertueux qui est l'amélioration de la qualité de mesure dans la mesure où les niveaux de répétabilité et de reproductibilité des mesures s'en trouvent réduits puisque les séquences d'essais s'enchaînent plus rapidement et que les conditions de réalisation sont plus homogènes. Le procédé selon l'invention est robuste et fiable en ce que, pour chaque mesure continue, on fournit une valeur représentative du signal stabilisé, contrairement à une prise de point traditionnelle où l'on observe une plus grande variabilité du résultat lorsque le signal présente une instabilité. Le procédé de l'invention permet en outre de prédire quand une valeur donnée du signal sera atteinte dans le futur. Dans le cas où le signal représentatif d'un paramètre de fonctionnement PF présente une dynamique assez lente, c'est-à-dire une durée de transition importante entre les deux asymptotes, et est de même ordre que la durée d'acqusition d, on sera en mesure d'estimer une valeur stable qui est assez différente des derniers échantillons acquis après cette durée d. Après avoir identifié les paramètres de l'équation (I) par exemple, on peut prédire le moment dans le futur qui sépare les derniers échantillons acquis d'une ( ( 1 -1 Ytgt YO valeur cible ytgt donnée sur le chemin vers yinfini trgt = --x ln - x ln Ystep Ainsi, il est possible soit de corriger les réglages appliqués soit d'appliquer d'autres réglages afin, dans un objectif de sécurité de fonctionnement, d'éviter de rencontrer une limite ou bien, dans un objectif de contrôle avancé, de piloter d'autres actionneurs pendant l'intervalle de temps disponible. Les étapes principales du procédé selon l'invention sont de préférence réalisées de manière automatique, avantageusement à l'aide d'un programme informatique. Ainsi, l'invention porte également sur un produit programme informatique permettant de conserver dans une mémoire d'une unité d'un processeur ou sur un support de mémoire amovible approprié pour coopérer avec ladite unité du processeur, le produit-programme comprenant des instructions pour mettre en oeuvre le procédé selon l'invention. L'invention a également pour objet un dispositif de test d'un moteur à combustion interne (banc de test moteur) permettant de réaliser des essais moteur dans le cadre d'une séquence d'essais. La figure 10 illustre un tel dispositif, permettant la mise en oeuvre du procédé selon l'invention. Le dispositif 100 selon l'invention comprend: - une instrumentation de mesure (102, 104, 105) pour la mesure continue d'au moins un paramètre de fonctionnement dudit moteur à combustion 101 générant un signal temporel brut (Si., Sib, S2a, S2b) lors d'un essai moteur. Cette instrumentation peut comprendre une ou plusieurs sondes de mesure (102, 104), par exemple fixées sur le moteur, sur ses organes (par exemple une vanne EGR ou un turbo-compresseur), ou sur des lignes de prélèvement. Les lignes de prélèvement sont des tuyaux souples chauffés, afin d'éviter toute condensation, qui sont chargés de convoyer les gaz d'échappement dont une partie du débit est aspiré, grâce à une pompe, depuis la ligne d'échappement vers les différents analyseurs, souvent plus d'une dizaine, de la baie d'analyse. Cette instrumentation comprend également le ou les conditionneurs 105 (ou convertisseurs) chargés de convertir les signaux acquis par les capteurs (qui peuvent délivrer par exemple une charge électrique ou une très faible tension) en un signal adapté au traitement par l'unité d'acquisition (par exemple une tension dans une plage allant de 0 à 10V) ; - un système d'automatisation 106 pour superviser la réalisation de l'essai moteur et pour afficher et/ou enregistrer les signaux acquis avec l'instrumentation de mesure (102, 104, 105). Ce superviseur 106 comprend donc un dispositif d'affichage et/ou une mémoire (unité d'enregistrement) pour l'affichage et/ou l'enregistrement des signaux; - un programme informatique selon l'invention pour déterminer une valeur représentative dudit signal stabilisé ; ledit système d'automatisation 106 pouvant afficher et/ou enregistrer l'ensemble des 10 valeurs représentatives déterminées par le programme informatique, en vue de leur analyse pour caractériser le fonctionnement dudit moteur à combustion 101; et ledit programme informatique pouvant être intégré audit système d'automatisation 106, ou à ladite instrumentation de mesure (102, 104, 105), ou encore à un ordinateur 107 distinct dudit système d'automatisation 106 et de ladite instrumentation de mesure (102, 104, 105), 15 ledit ordinateur 107 étant alors en communication avec ledit système d'automatisation 106. Le procédé selon l'invention peut être mis en oeuvre dans le cadre d'essais moteur réalisés sur bancs moteurs, mais également dans le cadre d'essais réalisés sur un véhicule en fonctionnement, c'est-à-dire circulant, comprenant ledit moteur. Dans ce dernier cas, le 20 véhicule comprend notamment une instrumentation embarquée permettant de réaliser une mesure continue d'un paramètre de fonctionnement du moteur. Le procédé selon l'invention donne alors accès à des informations sur le fonctionnement du moteur en temps réel ou en différé, pouvant potentiellement être utilisées pour calibrer automatiquement le moteur et/ou être utilisées par le conducteur afin de modifier sa conduite et atteindre un fonctionnement 25 optimal du véhicule (limitation des émissions de polluants, etc.). De manière avantageuse, la présente invention est utilisée pour calibrer le moteur à combustion interne. On calibre le moteur à combustion interne à partir de l'analyse de l'ensemble des valeurs représentatives des signaux stabilisés caractérisant le 30 fonctionnement dudit moteur. Il est ainsi possible de modifier des paramètres de réglage du moteur de manière à optimiser le fonctionnement du moteur, selon les spécifications prédéterminées. Ainsi, il est par exemple possible de modifier des paramètres de réglage du moteur liés à l'injection, à la masse d'air dans le cylindre, à la recirculation des gaz d'échappement (EGR), à la suralimentation, etc., en fonction de paramètres de 35 fonctionnement déterminés selon l'invention, de manière à limiter la consommation du moteur et/ou les émissions polluantes. Ces modifications des paramètres de réglage du moteur en fonction des informations issues du procédé selon l'invention peuvent également être effectuées en dehors du cadre des essais sur bancs moteurs, par exemple lors du fonctionnement du véhicule, comme mentionné plus bas. Le procédé selon l'invention peut par exemple être utilisé pour déterminer les dizaines de paramètres d'un modèle statistique de comportement pour chaque paramètre de fonctionnement PF du moteur, qui peuvent ensuite permettre de calibrer le moteur. Dans ce contexte, des séquences de plusieurs centaines d'essais organisés en plan d'expériences sont réalisées. Des sauts de point de puissance et de réglages du moteur d'un point à l'autre sont obtenus par la recherche d'un critère optimal (maximum du déterminant de la matrice des régresseurs pour un modèle polynômial ou maximum des distances minimales inter- point pour un modèle de type interpolatoire). Ce type de séquence, effectué de manière classique avec des prises de point traditionnelles, est très coûteux en temps (plusieurs dizaines d'heures, parfois 2 jours en continu). L'utilisation du procédé selon l'invention, pouvant permettre de diviser par deux la durée de la séquence, est donc particulièrement bien adaptée pour cette application. Les étapes suivantes peuvent être réalisées pour déterminer les modèles de comportement des paramètres de fonctionnement et caractériser le moteur, en vue par exemple de sa calibration: 1.1 lister les actionneurs du moteur à piloter; 1.2 définir les limites minimum et maximum de variation du réglage de chaque actionneur. Le procédé selon l'invention peut être mis en oeuvre à cette étape; 1.3 identifier les interdépendances entre les réglages qui ne peuvent pas varier indépendamment les uns des autres (par exemple le débit d'air frais aspiré et la pression de suralimentation). Le procédé selon l'invention peut être mis en oeuvre à cette étape; 1.4 calculer le plan d'expériences; 1.5 réaliser la séquence d'essais de manière plus rapide grâce au procédé selon l'invention en faisant se succéder les combinaisons de réglages inscrites dans le plan; 1.6 caler les paramètres de chaque modèle statistique de paramètres de fonctionnement du moteur en utilisant les valeurs finales estimées par le procédé selon l'invention; 1.7 optimiser les modèles statistiques et déterminer les réglages des actionneurs qui permettent de satisfaire les contraintes de fonctionnement du moteur; 1.8 réunir l'ensemble de ces réglages optima dans une table (une table par réglage) qui sera inscrite dans le calculateur électronique du moteur qui assure le contrôle du moteur équipant chaque véhicule produit en série ; l'ensemble de ces tables peut etre utilisée pour la calibration du moteur. L'invention peut également être utilisée pour caractériser de manière exhaustive l'évolution de chaque paramètre de fonctionnement PF du moteur dans sa plage d'utilisation en régime et en puissance. Pour cela, on effectue par exemple des variations croissantes par saut de couple (ou de la pression) moyen effectif (ou indiqué) sur une série croissante de niveaux de régime moteur. Ce type de séquence d'essais dure une dizaine d'heures mais se répète plusieurs fois au cours d'une étude ; c'est pourquoi diviser par deux la durée d'exécution est intéressant. Les étapes pour ce type d'essai peuvent être les suivantes : 2.1 charger dans le calculateur un ensemble de tables optimales (voir l'étape 1.8 décrite précédemment); 2.2 réaliser la séquence d'essai de manière plus rapide grâce au procédé selon l'invention en commençant à faible régime et faible couple puis en réalisant des sauts de couple à iso-régime jusqu'au couple maximum puis en augmentant le régime après avoir réduit le couple au minimum; 2.3 inscrire dans un tableau indexé en niveaux de régime et de couple les valeurs des paramètres de fonctionnement du moteur en utilisant les valeurs finales estimées par le procédé selon l'invention ; l'ensemble de ces tableaux forme une carte d'identité du moteur lorsqu'il est réglé avec ces tables optimales des réglages. Le procédé selon l'invention peut également être utilisé pour tester le moteur une fois celui-ci mis dans le commerce, afin d'évaluer son bon fonctionnement. En effet certains programmes qui s'exécutent dans le calculateur électronique ont pour but de vérifier que certains organes remplissent toujours leur rôle dans la dépollution des gaz d'échappement (ex. la vanne EGR ou le filtre à particules). Ces programmes détectent lorsque le véhicule en cours d'utilisation s'inscrit dans une suite logique d'états afin d'effectuer un relevé d'un paramètre de fonctionnement du moteur (ex. une recopie de position ou une différence de pression). Lorsque ce relevé n'est pas conforme à une valeur pré-enregistrée, suivant la gravité du dysfonctionnement du système de dépollution, soit un défaut est enregistré dans le calculateur ou bien un voyant peut être allumé au tableau de bord. Dans ce dernier cas, l'utilisateur est alors alerté d'un dysfonctionnement qui nécessite un retour en atelier. L'invention peut faciliter la réalisation du relevé et réduire les risques de détecter de manière erronée un dysfonctionnement. REVENDICATIONS1. Procédé de caractérisation du fonctionnement d'un moteur à combustion interne à partir d'une séquence d'essais moteur, chaque essai moteur comprenant la modification d'au moins un paramètre de réglage du moteur PR et/ou du point de puissance du moteur, comprenant les étapes suivantes: a) pour chaque essai moteur : I. on acquiert au moins un signal issu de la mesure continue d'un paramètre de fonctionnement PR dudit moteur sur une durée d'acquisition d, la durée d'acquisition d étant choisie de manière à pouvoir déterminer au moins une fonction sigmoïde pour modéliser ledit au moins un signal ; ii. on détermine ladite au moins une fonction sigmoïde modélisant ledit au moins un signal ; iii. on déduit de ladite au moins une fonction sigmoïde une valeur représentative dudit au moins un signal stabilisé ; b) on analyse l'ensemble desdites valeurs représentatives des signaux stabilisés pour caractériser le fonctionnement dudit moteur à combustion interne. 2. Procédé selon la revendication 1, dans lequel on fait varier rapidement, de préférence entre 0 et 5 secondes, au moins un paramètre de réglage PR et/ou le point de puissance du moteur après le déclenchement de l'acquisition du signal. 3. Procédé selon l'une des revendications précédentes, dans lequel la mesure continue est effectuée selon une fréquence d'échantillonnage comprise entre 1 Hz et 10 Hz. 4. Procédé selon l'une des revendications précédentes, dans lequel ladite au moins une fonction sigmoïde est une fonction de Gompertz. 5. Procédé selon l'une des revendications précédentes, dans lequel la représentation graphique de ladite au moins une fonction sigmoïde comprend une première asymptote horizontale dont la valeur est déterminée à un instant initial 1'0 du signal et une deuxième asymptote horizontale dont la valeur est déterminée à un instant infini t.. 6. Procédé selon la revendication 5, dans lequel on modélise le signal par une fonction sigmoïde s'écrivant selon l'équation (I) suivante:y(t) = yo + y step x exp(-b x exp(-c x (t - ))) (I) comprenant les cinq paramètres à déterminer yo, Ystep, b, c et yo étant la valeur du paramètre de fonctionnement PF à un instant initial t=t0; ystep étant l'écart entre la valeur de y(t) à un temps infini t=t., et 1.0; b et c étant des constantes positives, avec b caractérisant le décalage de la courbe dans le temps et c caractérisant la pente de la courbe ; ti étant étant une constante de temps positive 7. Procédé selon la revendication 5, dans lequel on modélise le signal par deux fonctions sigmoïdes. 8. Procédé selon la revendication 7, dans lequel les deux fonctions sigmoïdes sont combinées sous l'équation (II) suivante: y(t) = yo + y stepi x exp(-b1 x exp(-c1 x (t - t11))) + Ystep2 x exp(-b2 x exp(-c2 x (t - ti2))) (Il) comprenant les neuf paramètres à déterminer y0, 0, , V step1, Ystep2, b1, C1, ti1, b2, C2, t12; YO étant la valeur du paramètre de fonctionnement PF à un instant initial t=t0; Ystepl étant l'écart entre la valeur de y(t) à un temps infini t=t..> et celle à t, si un seul premier phénomène thermo-mécanique et/ou thermo-chimique existait, et représentant la variation du signal liée à ce premier phénomène ; Ystep2 étant l'écart entre la valeur de y(t) à un temps infini t=t.., si seul le premier phénomène existait et la valeur de y(t) à un temps infini t=t., si un seul deuxième phénomène thermomécanique et/ou thermo-chimique antagoniste au premier phénomène existait, Ystep2 représentant la variation du signal liée à ce deuxième phénomène ; Ystep1 et Ystep2 étant de signes opposés; bl et c1 étant des constantes positives, avec b1 caractérisant le décalage de la première courbe en S dans le temps, et cl caractérisant la pente de la première courbe en S; étant une constante de temps positive relative à l'inertie du premier phénomène !; b2 et c2 étant des constantes positives, avec b2 caractérisant le décalage de la deuxième courbe en S dans le temps, et c2 caractérisant la pente de la deuxième courbe en S; ti2 étant une constante de temps positive relative à l'inertie du deuxième phénomène. 9. Procédé selon la revendication 6 ou 8, dans lequel l'étape ii) de détermination de ladite au moins une fonction sigmoïde modélisant le signal comprend les étapes suivantes:- on sélectionne le signal ayant un critère de dynamique ou un critère d'évolution temporelle supérieur ou égal à une valeur seuil de dynamique prédéterminée ou une valeur seuil d'évolution prédéterminée (QUAL-S); - on détermine le paramètre yo correspondant à la valeur de la première asymptote horizontale (Val-AS1); - on calcule la dérivée temporelle du signal (CH-FS) ; - on définit la fonction sigmoïde selon l'équation (I) ou les deux fonctions sigmoïdes selon l'équation (II) selon que le signe de la dérivée temporelle reste le même ou change (CH-FS) ; - on calcule les autres paramètres de ladite/desdites fonctions sigmoïdes de manière à minimiser au sens des moindres carrés l'écart entre le signal acquis et le signal estimé selon ladite/desdites fonctions sigmoïdes (CAL-PFs) ; et dans lequel l'étape iii) de détermination de la valeur représentative du signal stabilisé comprend les étapes suivantes: - on détermine la valeur de la deuxième asymptote y.. (Val-AS2); - on calcule un coefficient de corrélation entre le signal acquis et la/les fonctions sigmoïdes et on sélectionne y.. en tant que valeur représentative du signal stabilisé si le coefficient de corrélation est supérieur à une valeur de corrélation prédéterminée (QUAL-corr, CAL-PF). 10. Procédé selon la revendication 9, dans lequel on augmente la durée d'acquisition d à l'issue de l'étape iii) de détermination de la valeur représentative du signal stabilisé de manière à ce que le coefficient de corrélation soit supérieur à une valeur de corrélation prédéterminée. 11. Procédé selon l'une des revendications précédentes, comprenant en outre une étape de filtrage du signal (FIL) lors de l'étape ii) de détermination de la fonction sigmoïde modélisant le signal, ledit filtrage comprenant un filtre de préférence choisi dans la liste des filtres suivants: moyenne mobile, filtre de Butterworth, filtre passe-bas de premier ordre. 12. Procédé selon l'une des revendications précédentes, dans lequel l'essai moteur est réalisé sur un banc de test moteur ou sur un véhicule en fonctionnement. 13. Procédé selon l'une des revendications précédentes, dans lequel on calibre le moteur à combustion interne à partir de l'analyse de l'ensemble des valeurs représentatives des signaux stabilisés caractérisant le fonctionnement dudit moteur. 14. Procédé selon l'une des revendications précédentes, utilisé pour caractériser le fonctionnement d'un moteur à allumage par compression. 15. Produit programme informatique permettant de conserver dans une mémoire d'une unité d'un processeur ou sur un support de mémoire amovible approprié pour coopérer avec ladite unité du processeur, le produit-programme comprenant des instructions pour mettre en oeuvre le procédé selon l'une quelconque de revendications 1 à 14. 16. Dispositif de test (100) d'un moteur à combustion interne (101) permettant de réaliser une séquence d'essais moteur, comprenant: - une instrumentation de mesure (102, 104, 105) pour la mesure continue d'au moins un paramètre de fonctionnement PF dudit moteur (101) générant un signal temporel brut (S1a, S1b, S2a, S2b) lors d'un essai moteur ; - un système d'automatisation (106) pour superviser la réalisation de l'essai moteur et pour afficher et/ou enregistrer les signaux acquis avec l'instrumentation de mesure (102, 104, 105) - un programme informatique selon la revendication 15 pour déterminer une valeur représentative dudit signal stabilisé ; ledit système d'automatisation (106) pouvant afficher et/ou enregistrer l'ensemble des valeurs représentatives déterminées par le programme informatique, en vue de leur analyse pour caractériser le fonctionnement dudit moteur à combustion (101); et ledit programme informatique pouvant être intégré audit système d'automatisation (106) ou à ladite instrumentation de mesure (102, 104, 105), ou encore à un ordinateur (107) distinct dudit système d'automatisation (106) et de ladite instrumentation de mesure (102, 104, 105), ledit ordinateur (107) étant en communication avec ledit système d'automatisation (106).