FR 3070761 A1 20190308 FR 1758238 A 20170907 L'invention concerne un procédé et un dispositif de mesure d'une marge de poussée d'une turbomachine. Un domaine d'application de l'invention est les turbomachines d'aéronefs, comme par exemple turboréacteurs d'avion ou turbopropulseurs d'avion. La principale caractéristique d'une turbomachine est sa poussée. Ainsi, il est fondamental de mesurer sa poussée avec précision. Dans l'état de la technique, la calibration des accessoires bancs est réalisée 10 sur des périodes très éloignées. Elle est habituellement déclarée faite avec une précision de +/-0,25% en termes de marge de poussée. Le document FR-A-1 416 112 décrit un dispositif d'étalonnage d'un banc de mesure de poussée d'un moteur de fusée. Le document FR-A-2 535 847 décrit une cellule d'essais de moteurs d'avion, 15 prévoyant une correction aérodynamique de la poussée mesurée. Le document US-A-2016/0 202 131 décrit la structure mécanique d'un capteur de mesure de poussée. Les inventeurs se sont aperçus de ce que les tests de réception pour la vérification des fonctions de la turbomachine sont effectués sur des bancs de test 20 différents, ce qui impacte les mesures de poussée. Un processus de certification impose le calcul d'une marge de poussée à partir de la poussée mesurée. On a ainsi constaté qu'il y avait une dispersion de la marge de poussée en raison des conditions de test. La dispersion de la marge de poussée résulte de la complexité des conditions de l'essai tels que : les ATC (Adaptations techniques 25 complémentaires), les bancs et les sites. Ainsi, il a été constaté que la mesure de la marge de poussée y est particulièrement influencée par certains composants importants tels que la nacelle secondaire et le banc. En outre, la mesure de la marge de poussée peut être influencée par le fait que les différentes parties de la turbomachine viennent de fournisseurs différents. 30 L'invention vise à obtenir un procédé et un dispositif de mesure, permettant d'effectuer des mesures de poussée depuis différents bancs d'essai avec différents équipements et de faire en sorte que ces mesures soient comparables. 2 A cet effet, un premier objet de l'invention est un procédé de mesure d'une marge de poussée d'une turbomachine, dans lequel on acquiert des données comportant la marge de poussée laquelle est déterminée en fonction d'une poussée spécifiée et d'une poussée mesurée, la poussée 5 mesurée étant déterminée sur un banc de mesure, qui comporte au moins un équipement de banc et sur lequel se trouve la turbomachine, caractérisé en ce que l'on modélise une évolution temporelle de la marge de poussée par au moins une fonction linéaire ou affine, qui est calculée pour au moins un intervalle de temps 10 déterminé de la marge de poussée, on calcule au moins un biais de l'équipement de banc par rapport à la au moins une fonction linéaire ou affine ayant été calculée, on corrige la marge de poussée en lui soustrayant le au moins un biais de l'équipement de banc dans le au moins un intervalle de temps déterminé. 15 Grâce à l'invention, la marge de poussée est rendue indépendante des conditions d'essai et sa dispersion est réduite. L'invention permet un meilleur contrôle de la production en réduisant l'incertitude sur la mesure de la poussée. En outre, l'invention permet de repérer les tendances d'efficacité de production de chaque fournisseur, ce qui permet d'alerter le fournisseur quand la poussée de la 20 turbomachine produite diminue et se rapproche de la marge limite. Suivant un mode de réalisation de l'invention, la fonction linéaire ou affine est calculée comme approchant un point prescrit ou passant par un point prescrit, le point prescrit ayant été calculé en fonction de la marge de poussée sur l'intervalle de temps déterminé. 25 Suivant un mode de réalisation de l'invention, le point prescrit est calculé comme étant un point moyen de la marge de poussée et au milieu de l'intervalle de temps déterminé. Suivant un mode de réalisation de l'invention, le point moyen de la marge de poussée est calculé comme étant au milieu de la marge de poussée et au milieu de 30 l'intervalle de temps déterminé. 3 Suivant un mode de réalisation de l'invention, le point prescrit a une valeur temporelle déterminée dans le au moins un intervalle de temps et une valeur de marge de poussée, la fonction linéaire ou affine est calculée comme minimisant, sur le au moins 5 un intervalle de temps, l'écart entre la marge de poussée et la fonction linéaire ou affine, rapporté à l'écart entre la valeur de marge de poussée du point prescrit et la valeur de la fonction linéaire ou affine à la valeur temporelle déterminée du point prescrit. Suivant un mode de réalisation de l'invention, la fonction linéaire ou affine 10 q(t) est calculée comme minimisant, sur le au moins un intervalle de temps, l'expression suivante : E Ptmoy)-(9,(t,)-9,(tptmoy))12 T 1(0 1,/ OÙ y, est la marge de poussée, Ti(i) est une fonction binaire indicatrice de chaque intervalle de temps (INT1), 15 tptmoy est la valeur temporelle déterminée du point prescrit dans le au moins un intervalle de temps, ptmoy est la valeur de marge de poussée du point prescrit dans le au moins un intervalle de temps. Suivant un mode de réalisation de l'invention, le au moins un biais est calculé 20 comme minimisant, sur le au moins un intervalle de temps, la marge de poussée à laquelle sont soustraits le biais et la fonction linéaire ou affine qui est calculée à l'instant de mesure ou jour de mesure correspondant à la marge de poussée. Suivant un mode de réalisation de l'invention, le au moins un biais bk est calculé comme minimisant, sur le au moins un intervalle de temps, l'expression 25 suivante : 2 y. - 9 (t k(i) OÙ y, est la marge de poussée, 4 8k(1) est une fonction binaire indicatrice du au moins un équipement k du banc de mesure, q(t) est la fonction linéaire ou affine. Suivant un mode de réalisation de l'invention, on tient compte de la fonction 5 yi(ti) indicatrice de l'origine des données depuis la source j de données ou du fournisseur j du au moins un équipement k de banc, et le ou les biais bk est calculé comme minimisant, sur le au moins un intervalle de temps, l'expression suivante : 2 y, - 9 ( )-Ibk kW OÙ y, est la marge de poussée, 1 0 15k(i) est une fonction binaire indicatrice du au moins un équipement k de banc, q(t) est la fonction linéaire ou affine, y(t) est la fonction indicatrice de l'origine des données depuis la source j de données ou du fournisseur j du au moins un équipement k de banc. 15 Suivant un mode de réalisation de l'invention, le au moins un intervalle de temps correspond à au moins un intervalle de temps de production et/ou à au moins un intervalle de temps, pendant lequel le au moins un équipement de banc est présent sur le banc de mesure, et/ou à l'intersection entre le au moins un intervalle de temps de production et le au moins un intervalle de temps, pendant lequel le au moins un 2 0 équipement de banc est présent sur le banc de mesure. Un deuxième objet de l'invention est un dispositif de mesure d'une marge de poussée d'une turbomachine, le dispositif comportant un banc de mesure, qui comporte au moins un équipement de banc et sur lequel est destiné à être montée la turbomachine, le banc de mesure comportant des moyens d'acquisition de données 25 comportant la marge de poussée de la turbomachine, caractérisé en ce que le dispositif comporte en outre une unité de calcul automatique configurée pour 5 modéliser une évolution temporelle de la marge de poussée par au moins une fonction linéaire ou affine, qui est calculée pour au moins un intervalle de temps déterminé de la marge de poussée, calculer au moins un biais de l'équipement de banc par rapport à la au moins 5 une fonction linéaire ou affine ayant été calculée, corriger la marge de poussée en lui soustrayant le au moins un biais de l'équipement de banc dans le au moins un intervalle de temps déterminé. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif en référence aux dessins annexés, 1 0 sur lesquels : - les figures 1, 8 et 11 représentent des mesures de marges de poussée en fonction du temps, - la figure 2 représente un changement de pente d'une fonction linéaire modélisant la marge de poussée en fonction du temps, suivant un mode de réalisation 15 du procédé et du dispositif de mesure suivant l'invention, - la figure 3 représente des intervalles de temps pouvant être utilisés pour définir la fonction linéaire de modélisation de la marge de poussée en fonction du temps suivant un mode de réalisation du procédé et du dispositif de mesure suivant l'invention, 20 - la figure 4 représente des points moyens pouvant être utilisés pour définir la fonction linéaire de modélisation de la marge de poussée en fonction du temps suivant un mode de réalisation du procédé et du dispositif de mesure suivant l'invention, - la figure 5 est un synoptique modulaire d'un mode de réalisation du procédé 25 et du dispositif de mesure suivant l'invention, - les figures 6, 9 et 12 représentent la marge de poussée débiaisées en fonction de l'instant de mesure dans le cas respectivement des mesures des figures 1, 8 et 11, - les figures 7, 10 et 13 représentent un résidu de la marge de poussée 30 débiaisée par rapport à son modèle de fonction linéaire en fonction de l'instant de mesure dans le cas respectivement des mesures des figures 1, 8 et 11. 6 A la figure 5, la turbomachine TUR est montée sur le banc de mesure BM pour mesurer la poussée exercée par la turbomachine TUR (pouvant être un turboréacteur d'avion ou un turbopropulseur d'avion ou autre) dans les conditions d'essai. Le banc BM de mesure comporte un ou plusieurs équipements k de banc et 5 des moyens d'acquisition de données de mesure (pouvant comprendre un ou plusieurs instruments de mesure, un ou plusieurs capteurs de mesure ou autres), permettant de mesurer un ou plusieurs paramètres de la turbomachine TUR, dont la poussée PST procurée par la turbomachine TUR lors de son fonctionnement sur le banc BM de mesure (par exemple par un ou plusieurs capteurs de poussée exercée 1 0 par la turbomachine TUR lors de son fonctionnement). Ci-dessous, les composants et équipements du banc BM de mesure sont désignés par la lettre k. Le banc BM de mesure peut comporter par exemple comme équipement k une nacelle secondaire et/ou autres. La turbomachine TUR comporte des aubes rotatives montées sur un rotor de celle-ci. Ces aubes rotatives peuvent provenir de différents fournisseurs. 15 Le dispositif de mesure met en oeuvre les étapes décrites ci-dessous du procédé de mesure suivant l'invention. Au cours d'une étape de mesure, on acquiert la marge de poussée de la turbomachine TUR sur le banc BM de mesure. Le banc BM de mesure comporte des moyens d'acquisition de données comportant la marge de poussée y, y, de la 2 0 turbomachine TUR. Par exemple, pour ce faire, on acquiert des données de mesure, comportant une poussée mesurée PST de la turbomachine TUR, par les moyens d'acquisition du banc BM de mesure. Pour chaque essai, la poussée PST est mesurée dans les conditions du jour t, de l'essai. Suivant un mode de réalisation, la poussée PST de la turbomachine TUR, mesurée par le banc BM de mesure, peut avoir été en 25 plus corrigée pour être ramenée dans des conditions standard, par exemple selon une procédure CFM, et étant alors appelée poussée standard PST. Puis, on calcule la marge de poussée y en fonction d'une poussée spécifiée PSP et de la poussée mesurée ou standard PST. Ce calcul est effectué par des moyens de calcul, comme par exemple l'unité de calcul automatique UC décrite ci- 30 dessous. La poussée spécifiée PSP est par exemple imposée selon un cahier des charges. 7 La marge de poussée y (ou marge de poussée y en réception) est formulée pour mesurer l'écart de la poussée mesurée ou standard PST rapporté à la poussée spécifiée PSP. La marge de poussée y (ou y, dans ce qui suit) peut être par exemple 5 calculée de la manière suivante PST - PSP Y - L'invention utilise un modèle MOD moyen d'évolution temporelle de la marge de poussée yi, ce modèle MOD étant linéaire ou affine par intervalles INT, de temps. Chaque intervalle de temps INT, peut couvrir une durée de plusieurs jours. 1 0 La marge de poussée y, est donnée en plusieurs instants de mesure ou jours de mesure, désignés par t,. A titre d'exemple, les figures 1 et 4 représentent la marge de poussée y, en ordonnées en fonction de l'instant t, de mesure ou jour t, de mesure, pour plusieurs équipements k du banc BM de mesure, comme par exemple la nacelle secondaire (en 15 anglais : cowl) N4, N7, N2, N5, N6, N9 du banc BM pour un premier fournisseur d'aubes rotatives de la turbomachine TUR. A titre d'exemple, la figure 8 représente la marge de poussée y, en ordonnées en fonction de l'instant t, de mesure ou jour t, de mesure, pour plusieurs équipements k de banc, comme par exemple la nacelle secondaire (en anglais : cowl) N4, N7, N2, 20 N5, N6, N9 du banc BM pour un deuxième fournisseur d'aubes rotatives de la turbomachine TUR. A titre d'exemple, la figure 12 représente la marge de poussée y, en ordonnées en fonction de l'instant t, de mesure ou jour t, de mesure, pour plusieurs équipements k de banc, comme par exemple la nacelle secondaire (en anglais : cowl) 25 N4, N7, N2, N5, N6, N9 du banc BM pour un troisième fournisseur d'aubes rotatives de la turbomachine TUR. On modélise une évolution temporelle de la marge de poussée y ou y, par au moins une fonction linéaire ou affine (I) ou cp,(t,), qui est calculée pour au moins un intervalle INT, de temps déterminé de la marge de poussée y ou y,. PSP 8 On calcule au moins un biais bk de l'équipement k de banc par rapport à la au moins une fonction linéaire ou affine (I) ou q(t) ayant été calculée. On corrige la marge de poussée y ou yi en soustrayant à cette marge de poussée y ou yi le au moins un biais bk de l'équipement k de banc dans le au moins 5 un intervalle INTI de temps déterminé. Pour une production donnée, on applique un modèle MOD selon lequel entre les instants correspondant à un intervalle de temps INTI identifié d'un composant k du banc BM on a une fonction linéaire ou affine (1), q(t) de la marge de poussée yi dépendant du temps t1. Il peut y avoir plusieurs composants k du banc BM, qui peuvent être ou comprendre par exemple le banc BM et la nacelle. Chaque entrée ou mesure est définie par un triplet (y, t, k), où y est la marge de poussée, ayant été mesurée de la manière décrite ci-dessus, t est la date d'essai, et k le composant du banc BM utilisé entre deux maintenances. Suivant un mode de réalisation de l'invention, on construit une fonction 8k(i) binaire pour chaque composant k du banc BM, cette fonction 8k(i) étant égale à 1 si l'observation est faite sur le composant k du banc BM et étant égale à 0 sinon. Cette fonction 8k(i) binaire est donc une fonction binaire indicatrice du composant k du banc BM. Suivant un mode de réalisation de l'invention, on résout le modèle MOD 20 suivant y, = b .6 k(i) + (I)(t ,) + £ où yi est la marge de mesure, ayant été calculée à partir du banc BM de mesure, l'indice i est l'indice de chaque donnée observée de cette marge de mesure yi 25 à un instant ou jour t1, le nombre de i étant fini et désigne chaque banc BM de mesure, q(t) est une fonction linéaire (ou affine) par intervalles INTI de temps sur le temps t,, qui représente l'évolution temporelle de la marge de poussée yi, bk est le biais introduit par le composant k du banc BM, 30 c est l'erreur de la mesure. 9 La fonction (1), q(t) est donc constituée de segments de droite en fonction du temps t,. Suivant un mode de réalisation, la fonction (1), q(t) linéaire ou affine par intervalles INT, de temps sur t, est continue, en pouvant changer de pente, au passage 5 d'un intervalle INT, de temps à l'intervalle INT, de temps suivant. Il peut y avoir une indétermination initiale des seuils de la marge de poussée, qui fait que le modèle MOD est positionné à une constante additive près. Suivant, un mode de réalisation, on prend une condition initiale sur le biais moyen bk, pouvant être par exemple : 10 uk -NL,/vk-uk OÙ Nk est le nombre de mesures faites avec le composant k du banc BM. Par exemple, Nk = k (i) Cela permet de résoudre cette indétermination initiale des seuils de la marge de poussée. 15 Suivant un mode de réalisation, on peut avoir des jeux différents de fonctions linéaires ou affines (1), cp,(ti) par intervalles INT, de temps sur t1, appelés ci-dessous (I),, pour lesquels les tendances sont spécifiques mais pour lesquels les biais bk doivent rester identiques, comme par exemple dans le cas de fournisseurs différents des composant k pour des tests similaires. On peut définir les fonctions (Di pour chaque tendance en gardant les biais bk et en tenant compte des fonctions indicatrices y,(ti) de sources j de données, par exemple indiquant le fournisseur j des composants k. La fonction y,(ti) indicatrice est égale à un pour la source j de données ou le fournisseur j du composant k et est égale à zéro pour les autres sources j de données ou autres fournisseurs j du composant k. La fonction (I) linéaire ou affine par intervalles INT, de temps sur t, peut être spécifique à une source de données j ou fournisseur j et tenir compte de la fonction y(t) indicatrice de l'origine des données depuis la source j de données ou du fournisseur j du composant k, par exemple selon l'équation suivante : 10 (Kt, ) = cl) (t )y (t ) La marge de poussée y, est alors égale à : y = bk .6k (i)+ cl) (t, )y (t, ) + £ Cela rajoute une nouvelle indétermination par les valeurs initiales des 5 fonctions T. Suivant un mode de réalisation, on calcule la ou les fonction (1), q(t) linéaire ou affine par intervalles INTI de temps sur t1. Suivant un mode de réalisation, une technique de résolution comprend une première étape de résolution, au cours de laquelle on décrit l'évolution (fonction (1), 1 0 q(t) linéaire ou affine par intervalles INTI de temps sur ti) de la marge de poussée y, indépendamment du niveau absolu résultant d'un ou plusieurs composants k du banc BM. Suivant un mode de réalisation, par exemple au cours de cette première étape de résolution, on associe chaque mesure y, de la marge de poussée à un ou plusieurs 15 premiers intervalles de temps INTI où les conditions d'essai sont constantes, appelé intervalles de temps INTI constant. On a par exemple L premiers intervalles de temps INTI ou L intervalles de temps INTI constant, pour 1 étant un entier variant de 1 à L, où L est un entier naturel non nul. Suivant un mode de réalisation de l'invention, le au moins un premier 2 0 intervalle de temps INTI correspond à au moins un deuxième intervalle de temps IPR de production et/ou à au moins un intervalle de temps ITC d'utilisation du composant k sur le banc BM de mesure. Suivant un mode de réalisation, on définit un ou des deuxièmes intervalles IPR de temps de production. Sur chacun de ces deuxièmes intervalles IPR de temps 25 de production, l'évolution (fonction (1), q(t) linéaire ou affine par intervalle de temps sur ti) de la marge de poussée est considérée différente. Cela est illustré à titre d'exemple à la figure 2, montrant en abscisse un intervalle de temps de production IPR1, où la fonction (I) est égale à une première fonction linéaire ou affine (Di, suivi d'un intervalle de temps de production IPR2, où la fonction (1) est égale à une 11 deuxième fonction linéaire ou affine e1e2 et un changement C de pente de la fonction (I) linéaire ou affine entre la première fonction linéaire ou affine e1e2 et la deuxième fonction linéaire ou affine e1e2 au passage P12 entre l'intervalle IPR1 et l'intervalle IPR2 (par exemple pour N4 précité). 5 Suivant un mode de réalisation, on définit un ou des troisièmes intervalles ITC de temps du composant k du banc BM, définis comme le ou les intervalles de temps où le composant k est présent sur le banc BM et a été utilisé pour effectuer des mesures y,. L'exemple de la figure 3 illustre l'exemple de définition de troisièmes intervalles de temps ITC1, ITC2, ITC3 d'un composant k étant l'une des nacelles, à 1 0 savoir la nacelle N9 précitée. Suivant un mode de réalisation, le premier intervalle de temps INT, ou intervalle de temps INT, constant est défini comme étant l'intersection entre un deuxième intervalle IPR de temps de production et un troisième intervalle ITC de temps du composant k. 15 Suivant un mode de réalisation, la fonction linéaire ou affine (1), q(t) est calculée comme approchant un point prescrit pt ou passant par un point prescrit pt. On calcule le point prescrit pt sur le premier intervalle de temps INT, ou intervalle de temps INT, constant déterminé, et ce en fonction de la marge de poussée y ou y, présente sur celui-ci. Le point prescrit pt a une valeur temporelle déterminée tptmoy 2 0 (selon l'axe des temps t, en abscisses aux figures 1 à 4) dans le premier intervalle de temps INT, ou intervalle de temps INT, constant et une valeur de marge de poussée ptmoy (selon l'axe des ordonnées y, aux figures 1 à 4). Un exemple de plusieurs points moyens pt est illustré à la figure 4 dans le cas des différentes nacelles N4, N7, N2, N5, N6 et N9 de la figure 1. 25 Suivant un mode de réalisation, on calcule le point prescrit pt comme étant un point moyen de la marge de poussée y ou yi, dont la valeur temporelle déterminée tp,moy est située sensiblement au milieu du premier intervalle de temps INT, ou intervalle de temps INT, constant déterminé. Suivant un mode de réalisation, on calcule le point moyen de la marge de 30 poussée y ou y, comme étant au milieu de la marge de poussée y ou y, et au milieu du premier intervalle de temps INT, ou intervalle de temps INT, constant déterminé. 12 On décrit ci-dessous une première étape de calcul de la fonction linéaire ou affine (I), (1),(0. Suivant un mode de réalisation, on calcule la valeur (pt(tptmoy) de la fonction linéaire ou affine (I), q(t) à la valeur temporelle déterminée tptmoy du point prescrit pt. 5 Suivant un mode de réalisation, la fonction linéaire ou affine (1), q(t) est calculée comme minimisant, sur le premier intervalle de temps INT1 ou intervalle de temps INT1 constant déterminé, l'écart entre la marge de poussée y ou yi et cette fonction linéaire ou affine (1), ((t), rapporté à l'écart entre la valeur de marge de poussée ptmoy du point prescrit pt et la valeur (pt(tptmoy) de la fonction linéaire ou 1 0 affine (1), q(t) à la valeur temporelle déterminée tptmoy du point prescrit pt, c'est-à- dire minimisant l'écart (y, - ptmoy) - (t,) -p1 (t pimoy)) sur le premier intervalle de temps INT1 ou intervalle de temps INT1 constant déterminé. Pour ce faire, suivant un mode de réalisation, la fonction linéaire ou affine (1), q(t) est calculée comme minimisant, sur le premier intervalle de temps INT1 ou 15 intervalle de temps INT1 constant déterminé, l'expression suivante : E=IkY,-Ptmoy)-(9 (t)-( (t ptmoy))12 T 1(1) 1,/ où Ti(i) est une fonction binaire pour chaque premier intervalle de temps INT1 ou intervalle de temps INT1 constant déterminé, qui est égale à un pour les mesures de marges de poussée yi faites pendant cet intervalle INT1 et 0 pour les mesures de 2 0 marges de poussée yi faites en dehors de cet intervalle INT1. La fonction binaire Ti(i) est donc une fonction binaire indicatrice du premier intervalle de temps INT1 ou intervalle de temps INT1 constant déterminé. On caractérise ainsi l'évolution temporelle (marge de poussée mesurée yi - point moyen pt) qui est indépendante des conditions de mesure. 25 Suivant un mode de réalisation, la fonction linéaire (1), q(t) est calculée comme minimisant, sur le premier intervalle de temps INT1 ou intervalle de temps INT1 constant déterminé, l'expression suivante : E =r; T 1(1) *[(); - Pt moy)- (9 j (t 1)- (t ptMoy)12 13 Suivant un mode de réalisation, la fonction linéaire ou affine d), q(t) est calculée comme étant égale à Cp, ayant l'expression suivante : 1Y (i)* T (1) ky, - (9, 9,( ptmoy)12} 01= Arg min ''j où Argmin désigne l'argument minimum sur les fonctions (pj linéaires ou 5 affines. Il peut y avoir une indétermination initiale des seuils, qui fait que le modèle MOD est positionné à une constante additive près. Lorsque la fonction linéaire ou affine (1), q(t) a été calculée au cours de la première étape, on calcule le ou les biais bk au cours d'une deuxième étape, décrite 1 0 ci-dessous. Suivant un mode de réalisation, le ou les biais bk est calculé comme minimisant, sur le premier intervalle de temps INT, ou intervalle de temps INT, constant déterminé, la marge de poussée y, à laquelle sont soustraits le biais bk et la fonction linéaire ou affine (1), q(t) qui est calculée à l'instant t, de mesure ou jour de 15 mesure t, correspondant à la marge de poussée yi, c'est-à-dire minimisant l'écart Yi-(p,(ti)-Ibk.6k(i) sur le premier intervalle de temps INT, ou intervalle de temps INT, constant déterminé. Pour ce faire, suivant un mode de réalisation, le ou les biais bk est calculé comme minimisant, sur le premier intervalle de temps INT, ou intervalle de temps 2 0 INT, constant déterminé, l'expression suivante : - 2 y1-9 (t)-Ibk.6k(i) - Suivant un mode de réalisation, on tient compte de la fonction y,(ti) indicatrice de l'origine des données depuis la source j de données ou du fournisseur j du composant k, et le ou les biais bk est calculé comme minimisant, sur le premier 25 intervalle de temps INT, ou intervalle de temps INT, constant déterminé, l'expression suivante : 2 14 y. - ço. (t 1) -lb k k (i) Suivant un mode de réalisation, le ou les biais bk est calculé comme étant égal à b ayant l'expression suivante : 2} bk= Arg min y , - (t ,) k .6 (i) bk 5 où Argmin désigne l'argument minimum sur le ou les biais bk. Après la deuxième étape, on corrige la marge de poussée yi pour calculer la marge de poussée débiaisée yidébiaisée, qui est égale à la marge de poussée yi à laquelle on a soustrait le ou les biais bk valable dans le ou les intervalles de temps INTI ou une moyenne des biais bk dans le ou les intervalles de temps INTI. 10 Suivant un mode de réalisation, on calcule la marge de poussée débiaisée yidébiaisée comme étant égale à yidébiaisée = où Fk est le biais moyen du composant k du banc BM. Suivant un mode de réalisation, on identifie les biais moyens Fk introduits par les composants k du banc 15 BM et éventuellement, dans le cas de présence de plusieurs fournisseurs j, les écarts entre les fournisseurs j (ou usineurs j). L'identification du biais moyen Fk est relative à une référence pour laquelle son biais associé est supposé nul. Suivant un mode de réalisation, on calcule la marge de poussée corrigée, débiaisée c'est-à-dire débiaisée, yi e comme étant égale à 20 débiaisée = Yi - 6k (/) A la figure 5, le dispositif 1 de mesure suivant l'invention prend en entrée un jeu d'essais et permet de débiaiser la marge yi de poussée et la rendre indépendante des composants k du banc BM. On décrit ci-dessous en référence à la figure 5 une utilisation d'un dispositif 1 25 de mesure suivant l'invention, au cours de laquelle la marge de poussée yi est 15 mesurée et corrigée, suivant un mode de réalisation de l'invention. Ce dispositif 1 prend en entrée un jeu d'essais et permet de débiaiser la marge y, de poussée et la rendre indépendante des composants k du banc BM. Suivant un mode de réalisation de l'invention, l'opérateur OP récupère les 5 données de mesure correspondant aux dernières périodes de production. Ces données ont été automatiquement archivées dans une base de données BDD. Suivant un mode de réalisation de l'invention, une interface homme-machine IHM (écran de visualisation ou autre) propose une interface de visualisation des mesures de poussée permettant d'extraire les moteurs dont les aubes (en anglais : 10 fan) proviennent d'un fournisseur donné et qui sont mesurés par un banc BM spécifique. Cette interface IHM permet par exemple de superposer de tels affichages en utilisant des couleurs différentes par fournisseur et par banc BM/nacelle, par exemple selon les figures 6 à 13. Suivant un mode de réalisation de l'invention, l'opérateur utilise l'interface 15 IHM pour identifier les intervalles de temps INTI durant lesquels les équipements de banc sont constants et les enregistrer dans la base de données BDD. Pour cela il utilise par exemple des journaux de bord (en anglais : logbooks) de maintenance et la visualisation des courbes d'évolution temporelle des mesures de poussée. Les instants délimitant les intervalles de temps INTI sont par exemple choisis par clics 20 avec la souris et sont automatiquement enregistrés. Une fois les instants de maintenance de bancs enregistrés, l'interface IHM permet de les visualiser sur les courbes. Suivant un mode de réalisation de l'invention, l'opérateur OP définit pour chaque fournisseur k des instants de changement de tendance initiaux, c'est-à-dire 25 les passages P12 entre intervalles de temps de production IPR1 et IPR2, par exemple par une initialisation manuelle. Pour cela il utilise par exemple des journaux de bord (en anglais : logbooks) de maintenance et la visualisation des courbes d'évolution temporelle des mesures de poussée. L'interface IHM permet par exemple de visualiser ces instants de changement de tendance initiaux P12, par exemple par une visualisation légèrement différente (par exemple des traits verticaux continus au lieu de traits pointillés pour les changements de configuration). 16 Suivant un mode de réalisation de l'invention, une unité de calcul automatique UC, telle que par exemple un ou plusieurs calculateur(s) et/ou un ou plusieurs ordinateur(s), et/ou un ou plusieurs processeur(s) et/ou un ou plusieurs serveur(s) et/ou une ou plusieurs machine(s), pouvant être programmées à l'avance 5 par un programme informatique préenregistré, exécute la modélisation automatique par le modèle MOD linéaire ou affine par intervalles de temps de la marge de poussée Yi pour un échantillon ECH d'un certain nombre n d'essais la correction automatique de la marge de poussée Yi et les étapes du procédé de mesure décrit ci-dessus. Cette unité de calcul automatique UC fournit la marge de poussée corrigée 10 Yi débiaisée. Suivant un mode de réalisation de l'invention, cette marge de poussée corrigée yidébiaisée est enregistrée dans la base de données BDD ou dans une autre mémoire, par exemple en association avec les instants ou jours ti, auxquels la mesure de poussée a été faite. Par exemple, l'unité de calcul automatique UC calcule les tendances de production des fournisseurs k et identifie les biais bk des bancs/nacelles 15 ainsi que décrit ci-dessus. La marge de poussée corrigée yidébiaisée et/ou eoles biais bk calculés sont par exemple enregistrés dans la base de données BDD ou dans une autre mémoire, par exemple en association avec les instants ou jours ti, auxquels la mesure de poussée a été faite, et/ou sont par exemple affichés et/ou communiqués sur l'interface IHM et/ou transmis vers l'extérieur sur au moins un port de sortie. Un 20 dispositif de transmission peut être prévu pour transmettre ces biais bk aux opérateurs OP de bancs qui peuvent corriger leurs mesures. Des itérations des étapes décrites ci-dessus peuvent être effectuées en optimisant la position et le nombre des instants de changement de tendance. Le calcul des biais correspond à une optimisation de la vraisemblance du 25 modèle MOD linéaire ou affine par morceaux. Il est donc possible de choisir automatiquement par des méthodes d'optimisation les instants P12 de changement de tendance de production. Suivant un mode de réalisation de l'invention, selon une première approche possible de choix des instants P12 de changement de tendance de production, une 30 approche locale consiste, à partir d'une définition pour chaque fournisseur k des instants de changement de tendance initiaux P12 telle que décrite ci-dessus, à 17 modifier localement ces instants de changement de tendance initiaux P12 par une descente de gradient. Cela ne change pas le nombre de points P12 de changement de tendance, mais optimise le modèle MOD de pentes. Suivant un mode de réalisation de l'invention, selon une deuxième approche 5 possible de choix des instants P12 de changement de tendance de production, on utilise une approche globale par une optimisation stochastique. Cette méthode consiste à tirer aléatoirement une population de sélection d'instants P12 (une liste de listes d'instants P12) et de calculer la qualité des modèles produits. Les meilleurs modèles des instants P12 sont conservés et une nouvelle population d'instants P12 à 1 0 tester est générée automatiquement à partir de ceux-ci, par mutation locale (quelques pas de descente de gradient ou des variations locales aléatoires des positions P12 avec une variance de plus en plus petite) et par combinaison de listes (croisement). Au bout de quelques itérations, les listes survivantes correspondent aux meilleurs modèles. On a ainsi un mécanisme de type « génétique » permettant d'identifier les 15 meilleurs instants P12 de changement de tendance ainsi que leur nombre. La figure 6 représente la marge de poussée débiaisée yidébiaisée en fonction de l'instant ti de mesure ou jour ti de mesure, pour plusieurs équipements k de banc, comme par exemple la nacelle secondaire (en anglais : cowl) N4, N7, N2, N5, N6, N9 pour le premier fournisseur d'aubes rotatives de la turbomachine TUR, à partir 2 0 de la marge de poussée yi des figures 1 et 4. La figure 7 représente le résidu R de la marge de poussée débiaisée yidébiaisée de la figure 6 par rapport à son modèle MOD de fonction linéaire ou affine (1) par intervalles de temps, qui possède un écart-type de 0.29, alors que l'écart-type des marges de poussée yi des figures 1 et 4 est de 0.49. La figure 9 représente la marge de poussée débiaisée yidébiaisée en fonction de 25 l'instant ti de mesure ou jour ti de mesure, pour plusieurs équipements k de banc, comme par exemple la nacelle secondaire (en anglais : cowl) N4, N7, N2, N5, N6, N9 pour le deuxième fournisseur, à partir de la marge de poussée yi de la figure 8. La figure 10 représente le résidu R de chaque marge de poussée débiaisée yidébiaisée de la figure 8 par rapport à son modèle MOD de fonction linéaire ou affine (1) par 30 intervalles de temps, qui possède un écart-type de 0.29, alors que l'écart-type des marges de poussée yi de la figure 8 est de 0.54. 18 La figure 12 représente la marge de poussée débiaisée yidébiaisée en fonction de l'instant ti de mesure ou jour ti de mesure, pour plusieurs équipements k de banc, comme par exemple la nacelle secondaire (en anglais : cowl) N4, N7, N2, N5, N6, N9 pour le troisième fournisseur, à partir de la marge de poussée yi de la figure 11. 5 La figure 13 représente le résidu R de chaque marge de poussée débiaisée yidébiaisée de la figure 12 par rapport à son modèle MOD de fonction linéaire ou affine (1) par intervalles de temps, qui possède un écart-type de 0.28, alors que l'écart-type des marges de poussée yi de la figure 11 est de 0.47. Ainsi que le montrent ces figures 6, 7, 9, 10, 12 et 13 par comparaison aux 1 0 figures 1, 8 et 11, le procédé et le dispositif suivant l'invention ont permis de réduire de plus de 40% environ l'écart-type et donc la dispersion de la marge de poussée yi. L'invention permet d'assurer l'indépendance du traitement de l'évolution de la marge de poussée et la quantification des biais des composants du banc. En outre, l'invention permet de déterminer l'évolution de la marge de poussée associée à des 15 aubes fabriquées par des fournisseurs différents et de quantifier les biais introduits par les composants de bancs sans interaction entre ces 2 traitements. L'invention permet d'attribuer de façon quasi-certaine un résultat de poussée s'écartant par rapport à la spécification du moteur soit à une dérive des aubes fabriquées soit à une particularité du système d'essai. L'invention permet de contrôler la calibration des 2 0 bancs BM et de leurs accessoires k d'une façon précise et continue en se basant sur les mesures opérationnelles. L'invention permet également de caractériser les tendances qualité des aubes fournies par les différents fabricants. L'invention permet de rendre les mesures de poussée plus robuste vis-à-vis des évolutions des conditions d'essai. L'invention permet un meilleur contrôle de la production en réduisant 25 l'incertitude sur la mesure de la poussée. Bien entendu, les modes de réalisation, caractéristiques et exemples ci-dessus peuvent être combinés l'un avec l'autre ou être sélectionnés indépendamment l'un de l'autre. REVENDICATIONS1. Procédé de mesure d'une marge de poussée d'une turbomachine (TUR), dans lequel on acquiert des données comportant la marge de poussée (y, yi) laquelle est déterminée en fonction d'une poussée spécifiée (PSP) et d'une poussée mesurée (PST), la poussée mesurée (PST) étant déterminée sur un banc (BM) de mesure, qui comporte au moins un équipement (k) de banc et sur lequel se trouve la turbomachine (TUR), 1 0 caractérisé en ce que l'on modélise une évolution temporelle de la marge de poussée (y, yi) par au moins une fonction linéaire ou affine (4), ((t)), qui est calculée pour au moins un intervalle (INT1) de temps déterminé de la marge de poussée (y, yi), on calcule au moins un biais (bk) de l'équipement (k) de banc par rapport à la 15 au moins une fonction linéaire ou affine (4), ((t)) ayant été calculée, on corrige la marge de poussée (y, yi) en lui soustrayant le au moins un biais (bk) de l'équipement (k) de banc dans le au moins un intervalle (INT1) de temps déterminé. 2. Procédé de mesure suivant la revendication 1, caractérisé en ce que la 2 0 fonction linéaire ou affine (4), ((t)) est calculée comme approchant un point prescrit (pt) ou passant par un point prescrit (pt), le point prescrit (pt) ayant été calculé en fonction de la marge (y, yi) de poussée sur l'intervalle (INT1) de temps déterminé. 3. Procédé de mesure suivant la revendication 2, caractérisé en ce que le point prescrit (pt) est calculé comme étant un point moyen de la marge de poussée 25 (y, yi) et au milieu de l'intervalle (INT1) de temps déterminé. 4. Procédé de mesure suivant la revendication 3, caractérisé en ce que le point moyen de la marge de poussée (y, yi) est calculé comme étant au milieu de la marge de poussée (y) et au milieu de l'intervalle (INT1) de temps déterminé. 5. Procédé de mesure suivant l'une quelconque des revendications 2 à 4, 30 caractérisé en ce que le point prescrit (pt) a une valeur temporelle déterminée (tptmoy) 20 dans le au moins un intervalle de temps (INT1) et une valeur de marge de poussée (Ptmoy), la fonction linéaire ou affine (4), ((t)) est calculée comme minimisant, sur le au moins un intervalle de temps (INT1), l'écart entre la marge de poussée (y, yi) et la 5 fonction linéaire ou affine (4), q(0), rapporté à l'écart entre la valeur de marge de poussée (ptmoy) du point prescrit (pt) et la valeur (q),(tptmoy)) de la fonction linéaire ou affine (4), ((t)) à la valeur temporelle déterminée (tptmoy) du point prescrit (pt). 6. Procédé de mesure suivant l'une quelconque des revendications 2 à 5, caractérisé en ce que la fonction linéaire ou affine q(t) est calculée comme 10 minimisant, sur le au moins un intervalle de temps (INT1), l'expression suivante : E=IkY,-Ptmoy)-(9,(t)-( (tptmoy))12 - T 1(1) 1,/ où yi est la marge de poussée, Ti(i) est une fonction binaire indicatrice de chaque intervalle de temps (INT1), tptmoy est la valeur temporelle déterminée du point prescrit (pt) dans le au 15 moins un intervalle de temps (INT1), ptmoy est la valeur de marge de poussée du point prescrit (pt) dans le au moins un intervalle de temps (INT1). 7. Procédé de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que le au moins un biais (bk) est calculé comme 20 minimisant, sur le au moins un intervalle de temps (INT1), la marge de poussée (yi) à laquelle sont soustraits le biais (bk) et la fonction linéaire ou affine (4), ((t)) qui est calculée à l'instant (t) de mesure ou jour de mesure (t) correspondant à la marge de poussée (yi). 8. Procédé de mesure suivant l'une quelconque des revendications 25 précédentes, caractérisé en ce que le au moins un biais bk est calculé comme minimisant, sur le au moins un intervalle de temps (INT1), l'expression suivante : - 2 y1 - (t ) - b k k (i) où yi est la marge de poussée, 21 8k(1) est une fonction binaire indicatrice du au moins un équipement k du banc (BM) de mesure, q(t) est la fonction linéaire ou affine. 9. Procédé de mesure suivant l'une quelconque des revendications 5 précédentes, caractérisé en ce que on tient compte de la fonction y(t) indicatrice de l'origine des données depuis la source j de données ou du fournisseur j du au moins un équipement k de banc, et le ou les biais bk est calculé comme minimisant, sur le au moins un intervalle de temps (INT1), l'expression suivante : 2 y, -9,(t,)-Ibk.6k(i) OÙ yi est la marge de poussée, 8k(i) est une fonction binaire indicatrice du au moins un équipement k de banc, q(t) est la fonction linéaire ou affine, y(t) est la fonction indicatrice de l'origine des données depuis la source j de données ou du fournisseur j du au moins un équipement k de banc. 10. Procédé de mesure suivant l'une quelconque des revendications précédentes, caractérisé en ce que le au moins un intervalle de temps (INTI) correspond à au moins un intervalle de temps (IPR) de production et/ou à au moins un intervalle de temps (ITC), pendant lequel le au moins un équipement (k) de banc est présent sur le banc (BM) de mesure, et/ou à l'intersection entre le au moins un intervalle de temps (IPR) de production et le au moins un intervalle de temps (ITC), pendant lequel le au moins un équipement (k) de banc est présent sur le banc (BM) de mesure. 11. Dispositif de mesure d'une marge de poussée d'une turbomachine (TUR), le dispositif comportant un banc (BM) de mesure, qui comporte au moins un équipement (k) de banc et sur lequel est destiné à être montée la turbomachine (TUR), le banc (BM) de mesure comportant des moyens d'acquisition de données comportant la marge de poussée (y, yi) de la turbomachine (TUR), 22 caractérisé en ce que le dispositif comporte en outre une unité (UC) de calcul automatique configurée pour modéliser une évolution temporelle de la marge de poussée (y, yi) par au moins une fonction linéaire ou affine (4), q(0), qui est calculée pour au moins un 5 intervalle (INTI) de temps déterminé de la marge de poussée (y, yi), calculer au moins un biais (bk) de l'équipement (k) de banc par rapport à la au moins une fonction linéaire ou affine (4), cp,(0) ayant été calculée, corriger la marge de poussée (y, yi) en lui soustrayant le au moins un biais (bk) de l'équipement (k) de banc dans le au moins un intervalle (INTI) de temps 1 0 déterminé.