Brevet FR3096136A1 — publié le 20201120

Description Titre de l'invention : Analyse du vieillissement d'une turbomachine par analyse de modèles prédictifs Domaine technique [0001] La présente invention concerne la mise à jour de modèles prédictifs dans le cadre d'une turbomachine. [0002] En propos liminaire, plusieurs définitions sont données. Comme illustré en figure 1, on se place dans le cadre d'une turbomachine 1 comprenant deux compresseurs 2, 3 successifs (compresseur basse pression 2 et haute pression 3) suivi d'une chambre de combustion 4. Ces définitions sont applicables pour toute la demande. [0003] La Ps3 est la pression statique mesurée ou calculée dans le plan en amont de la chambre de combustion). [0004] Xn12R est le régime du compresseur basse pression 2, réduit sur la température dudit compresseur T12 (pour s'affranchir des variations de température), exprimé en tours par minute. [0005] PCN12R (ou Ni dans le cas d'un entrainement direct ou « direct drive » dans la ter- minologie anglo-saxonne) est le régime du compresseur basse pression 2, réduit sur la T12 (pour s'affranchir des variations de température), exprimé en pourcentage du régime basse pression maximum. 100061 Xn25R est le régime du compresseur haute pression 3, réduit sur la T25 (pour s'affranchir des variations de température), exprimé en tours par minute. 100071 PCN25R (ou N2) est le régime du compresseur haute pression 3, réduit sur la tem- pérature dudit compresseur T25 (pour s'affranchir des variations de température), exprimé en pourcentage du régime haute pression maximum. 100081 PT2 est la pression totale extérieure (fournie par l'avion). 100091 P25 est la pression statique modélisée dans le compresseur haute pression. 100101 Un modèle est une loi mathématique décrivant l'évolution d'une grandeur physique (paramètre) en fonction d'un ou plusieurs variables physiques. Technique antérieure [0011] En exploitation, il arrive que des turbomachines subissent des fausses détections de pompages (décrochage des aubes d'un des cieux compresseurs) en phase de croisière Ces évènements ont un fort impact opérationnel (endoscopie du moteur). [0012] Dans ces deux cas, il a été observé au moment où les évènements ont eu lieu une panne d'écart entre les deux voies Ps3, c'est-à-dire entre les deux voies d'acquisition de la pression statique en amont de la chambre de combustion. [0013] L'impact des fausses détections de pompage a un impact opérationnel important dans 2 le sens où l'avion est immobilisé tant que le moteur n'a pas été endoscope pour vérifier s'il a subi des dommages. [0014] La ligne d'acquisition Ps3 est parfois constituée d'une canalisation qui prélève la pression en amont de la chambre de combustion 4 et de deux capteurs de pression situés directement dans le calculateur de l'avion (FADEC, pour full authority digital engitze control°. [0015] La mesure de Ps3 est effectuée grâce à deux capteurs indépendants. Afin de consolider les informations issues des deux capteurs, une logique de sélection entre les deux capteurs a été mise en place. On suppose ici que les capteurs relèvent des mesures valides (pas de panne électrique et mesure contenue dans une plage de mesures physiquement vraisemblables), et que les deux capteurs relèvent des mesures en écart l'une par rapport à l'autre. Cette configuration provoque une panne d'écart, mais il est en l'état impossible de voter pour l'une ou l'autre des mesures, ne sachant laquelle est la plus proche de la valeur de Ps3 réelle. [0016] Pour pallier ce problème, un modèle de Ps3 basé sur des lois thermodynamiques est calculé. Ce modèle permet théoriquement de lever le doute en fournissant une troisième grandeur (redondance analytique), indépendante des mesures de Ps3, qui permettra de voter pour l'un ou l'autre des relevés via la logique de sélection. La figure 2 illustre ce principe, avec les deux voies d'acquisition V10, V20, le modèle mod Ps3 et la bascule qui se fait lorsque la voie V10 redevient plus proche du modèle mod_Ps3 que la voie V20 qui avait divergé de la voie V10 auparavant. La bascule entraine au niveau du calculateur l'observation d'une variation de pression APs3 importante [0017] Néanmoins, on observe en pratique des valeurs du modèle assez éloignées de la valeur de Ps3 réelle. Cela peut conduire à un arbitrage de voie erroné. Le Déposant s'est aperçu, après des études, que la fausse détection de pompage était due à un brusque changement de sélection de Ps3 : comme les deux mesures de Ps3 étaient en écart, la voie sélectionnée est passée de la mesure la plus forte en Ps3 à la mesure la plus faible en un pas de calcul puisque le modèle était initialement plus proche de la Ps3 la plus importante pour ensuite se rapprocher de la Ps3 la plus faible. C'est ce faux saut APs3 d'au moins 15% de valeur relative qui peut déclencher une fausse détection de pompage alors que la pression n'a pas chuté réellement. [0018] 11 existe donc un besoin de se prémunir contre ce genre d'évènement notamment en améliorant la gestion de l'arbitrage, en particulier au sujet de la pression Ps3, mais pour tout autre paramètre. [0019] Plus généralement, il existe un besoin de mieux traiter les modèles thermody- namiques, afin qu'ils reflètent davantage la réalité, que ce soit pour la Ps3 ou d'autres paramètres. 3 100201 De plus, diverses améliorations ou utilisations du modèle thermodynamique pourraient être apportées pour améliorer la rapidité, l'efficacité et la pertinence des modèles thermodynamiques. Exposé de l'invention [0021] Un but de l'invention est de proposer des solutions aux problèmes mentionnés. [0022] Il est à cet effet. proposé, un procédé d'analyse de vieillissement d'une turbomachine, le procédé consistant à mettre en oeuvres les étapes suivantes : - F1 : Effectuer un recalage d'un modèle d'un paramètre physique de la turbomachine, le modèle étant une loi définissant la valeur du paramètre en fonction d'au moins une variable, - F2 : Sauvegarder le modèle recalé dans une mémoire non-volatile, les étapes FI et F2 étant répétées au moins deux fois, et préférablement plus, - F3 : Comparer les différents modèles recalés pour en déduire unc évolution de l'état de la turbomachine. [0023] Dans un mode de réalisation, le recalage des modèles se fait à l'aide d'une mesure d'une valeur du paramètre par un capteur ou d'une mesure d'un paramètre tiers. [0024] Dans un mode de réalisation, le recalage des modèles se fait à l'aide d'un correcteur PID, préférablement un correcteur uniquement intégral. [0025] Dans un mode de réalisation, les modèles sont définis par segment ou par plan. [0026] Dans un mode de réalisation, au moins deux modèles recalés ont été générés à plus d'un mois d'intervalle. [0027] En complément ou indépendamment, dans un mode de mise en oeuvre possible, le modèle recalé est un modèle de pression statique en amont de chambre de combustion (mod Ps3(PCN25R)), dit « modèle de Ps3», dans une turbomachine comprenant un compresseur, le modèle de Ps3 étant utilisé pour arbitrer entre deux voies d'acquisition de la pression statique en amont de la chambre de combustion (Ps3), dite « pression Ps3», les deux voies d'acquisition faisant intervenir deux capteurs, le procédé utilisant un modèle de Ps3 stocké dans une mémoire, le modèle exprimant la pression Ps3 en fonction au moins du régime du compresseur (3), dit « régime PCN25R » et comprenant les étapes suivantes : El : mesure d'une valeur de la pression Ps3, par un des deux capteurs, E2 : recalage du modèle de Ps3 à l'aide de la mesure de la valeur de Ps3. [0028] En variante, dans un autre mode de mise en oeuvre, le modèle recalé est un modèle de paramètre de fonctionnement (mod_PARAM) de turbomachine ou d'aéronef, le modèle étant défini comme une loi par segment indiquant la valeur dudit paramètre en fonction d'une variable (Mod_PARAM(Var)), ou étant défini comme une loi par plan indiquant la valeur dudit paramètre en fonction de deux variables (mod PARAM(Varl, Var2)' ladite loi étant affine sur chaque segment ou étant affine sur chaque plan, le modèle de paramètre étant stocké clans une mémoire, le procédé de recalage comprenant les étapes suivantes : - obtention d'une valeur du paramètre (Val PARAM), - calcul (l'une erreur (s) par comparaison de ladite valeur du paramètre (Val_PARAM) avec la valeur correspondante du modèle (Val mod PARAM), ladite valeur du modèle (Val mod PARAM) appartenant à un des segments ou des plans du modèle (mod PARAM), - application d'un correcteur en minimisant ladite erreur (e) pour déterminer une correction (corr), - recalage du segment du modèle (mod_PARAM) ou du plan du modèle à l'aide de la correction (corr), pour repositionner ledit segment ou plan et obtenir ainsi un modèle recalé du paramètre physique. Brève description des dessins [0029] D'autres caractéristiques, buts et avantages de l'invention ressortiront de la des- cription qui suit, qui est purement illustrativc et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : [0030] [fig.1] La figure 1 illustre de façon schématique une turbomachine. [0031] [fig.2] La figure 2 illustre un procédé d'arbitrage entre deux voies d'acquisition à l'aide d'un modèle thermodynamique. [0032] [fig.3] La figure 3 illustre graphiquement un procédé dc recalage dc la pression Ps3. [0033] [fig.4] La figure 4 illustre par un schéma-bloc un procédé dc recalage d'un modèle de paramètre, tel que la pression Ps3. [0034] W8.51 La figure 5 illustre un correcteur. [0035] [fig.64 [0036] [fig.611] Les figures 6a et 6h illustrent des procédés de recalage d'un modèle 2D par segment. [0037] [fig.74 La figure 7a illustre, pour un segment, un procédé dc recalage d'un modèle 2D en segment par pondération. [0038] [fig.7b] La figure 7b illustre, pour plusieurs segments, un procédé de recalage d'un modèle 2D en segment par pondération. [0039] [fig.8] La figure 8 illustre un modèle 3D par plan. [0040] [fig.9] La figure 9 illustre par un schéma-bloc un procédé de recalage d'un modèle 3D de paramètre, tel que la pression Ps3, en fonction des pressions PCN12R et PCN25R. [0041] [fig.10a] La figure 10a illustre, pour un plan, un procédé de recalage d'un modèle 3D en seement par pondération. [0042] [fig.1011 La figure 10h illustre le choix d'un triangle parmi le rectangle formant un plan du modèle 3D. [0043] [fig.10c] La figure 10c illustre le choix de la pondération pour un triangle parmi le rectangle formant un plan du modèle 3D. [0044] [fig.1 1] La figure 11 illustre par un schéma-bloc une sélection de modèle en fonction d'une variable, préalablement au recalage du modèle. [0045] [fig.12] La figure 12 illustre un procédé d'analyse de vieillissement de turbomachine. Description des modes de réalisation [0046] Le contexte et les définitions données en introduction sont repris ici. [0047] Dans un premier temps, un procédé de recalage de modèle de pression statique en amont de chambre de combustion va être décrit. Cette pression sera appelée pression Ps3 et ce modèle sera appelé « modèle de Ps3 » et référencé mod_Ps3. Il s'agit d'un modèle thermodynamique. [0048] Lc modèle de Ps3 a notamment pour but final de permettre d'arbitrer entre deux voies d'acquisition V 10, V20, redondantes, dont la fonction est de mesurer la pression Ps3. Chaque voie d'acquisition V 10, V20 comprend un capteur 10, 20. Le capteur 10, 20 est standard et ne sera pas décrit ici. [0049] On décrira par la suite un procédé d'arbitrage entre les deux voies d'acquisition V 10, V20. [0050] Une unité de calcul 100 est prévue, qui comprend un processeur 110 et une mémoire 120. L'unité de calcul 100 peut être un FADEC («full authority digital engine control ») ou bien être un composant distinct, positionné au plus près des voies d'acquisition V 10, V20 pour davantage de réactivité. [0051] La mémoire 120 stocke un modèle mod Ps3, qui permet d'obtenir la valeur de la pression PS3 en fonction d'au moins une variable Var, qui est le régime PCN25R (régime compresseur haute pression) : le modèle mod_ Ps3 est alors écrit sous la forme mod Ps3(PCN25R). En pratique, le modèle mod_Ps3 fait intervenir plusieurs sous-modèles, comme notamment le modèle de Ps3 sur la pression de compresseur haute-pression P25 (cc modèle est appelé mod_Ps3/P25) et le modèle mod_Ps3/P25 est lui-même exprimé en fonction du régime du compresseur haute-pression PCN25R réduit sur sa température T25. Ce modèle est alors écrit sous la forme mod_Ps3/P25(PCN25R/T25). [0052] 11 suffit ensuite de multiplier la valeur de Ps3/P25 par P25 pour obtenir la valeur de Ps3. [0053] Plutôt que de recaler directement le modèle mod_Ps3, il est ainsi préférable de recaler le modèle mod_Ps3/P25. La dénomination de « modèle de Ps3 », sous la forme 6 mod Ps3, inclut les modèles n'exprimant pas directement la pression Ps3 niais permettant de l'obtenir par la suite, tel que le modèle mod_Ps3/P25. [0054] Dans une première étape El, une des deux voies d'acquisition V10, V20, à l'aide de son capteur 10, 20, mesure une valeur Val Ps3 de la pression Ps3 sur la turbomachine (pour une valeur réelle de la grandeur physique qui est utilisée comme variable, c'est-à-dire PCN25R). A ce stade, on suppose que les deux voies d'acquisition V10, V20 sont saines et que les deux capteurs 10, 20 donnent une mesure correcte. En d'autres termes, il n'y a pas de panne de capteurs 10, 20 ou d'écart au-delà d'un seuil prédéterminé entre les deux mesures. [0055] Cette mesure d'une valeur Val Ps3 de la pression Ps3 est ensuite envoyée à l'unité de calcul 100. [0056] Une étape E2 de conversion ou de traitement des données peut être mise en oeuvre : par exemple, Val Ps3 est une valeur de pression statique Ps3, alors que le modèle mod_Ps3/P25 utilise la pression Ps3 réduite sur la P25 : il faut donc diviser la valeur de la pression statique par P25 pour obtenir la valeur Val Ps3/P25. [0057] Ensuite, dans une étape E3, l'unité de calcul 100 recale le modèle de Ps3 stocké dans sa mémoire 120 à l'aide de ladite mesure de la valeur de la pression Ps3. Par recaler, on entend qu'il existe au moins un point du modèle mod_Ps3 (en pratique une pluralité, voire une infinité, si le modèle est continu) dont l'ordonnée a été déplacée (à abscisse constante donc). On note Rmod PS3/P25 le modèle recalé. Par la suite, on simplifiera l'écriture en gardant mod_PS3/F'25 qui désigne un modèle avant et après recalage. [0058] En l'espèce, il existe au moins un point P de la courbe mod_Ps3(Var) dont la valeur Val_mod_Ps3(Var) a changé avant et après le recalage, pour une valeur de la variable donnée. Dans le mode de réalisation préféré, on travaille avec mod_Ps3/P25(PCN25R) et Val_mod_Ps3/P25(PCN25R). [0059] On définit enfin une étape E4 de stockage du modèle de Ps3 recalé dans la mémoire 120. Dans un mode de réalisation, le modèle recalé mod_Ps3 (en l'occurrence mod_Ps3/P25) remplace en supprimant dans la mémoire 120 le modèle précédent. Dans un autre mode de réalisation, il ne le supprime pas. [0060] Préférablement les étapes El, E2 et E3 sont répétées à intervalle régulier, de type à chaque pas de calcul. Le pas de calcul correspond à environ-0,015s. Durant un pas de calcul, les deux étapes El et E3 peuvent être mises en oeuvre ou bien une étape El et en parallèle l'étape E3 à l'aide des données issues del 'étape El du pas précédent sont mises en oeuvre. [0061] Plus le pas de calcul est court et plus le modèle mod_Ps3 pourra être recalé ra- pidement, notamment pour s'adapter à des conditions de vol qui se réalisent rapidement et dont l'influence est critique : on peut citer les prélèvements d'air dans les 7 compresseurs (connus aussi sous l'expression « bleed air » dans la terminologie anglo-saxonne) par exemple. [0062] Comme le modèle mod Ps3 est mis à jour à intervalle régulier, l'arbitrage peut se faire plus rapidement et donc plus correctement, en évitant les sauts APs3 liés au changement de voie V10, V20 intempestif. [0063] Le recalage s'effectue avantageusement à l'aide d'un correcteur 112 qui est intégré dans une boucle de la chainc de contrôle. Cc correcteur sera décrit en détail par la suite. [0064] On définit aussi un procédé d'arbitrage entre deux voies d'acquisitions V10, V20, le procédé d'arbitrage comprenant une étape Al de mise en oeuvre d'un procédé de recalage tel que décrit précédemment et d'une étape A2 de choix de la voie d'acquisition V10, V20, durant laquelle le processeur choisit une voie V10, V20 parmi les deux voies V10, V20. Le choix se fait en fonction de la voie d'acquisition V10, V20 qui est la plus proche du modèle recalé. L'étape A2 est classique et ne sera pas décrite ici. [0065] Dans un deuxième temps, un procédé spécifique de recalage d'un modèle mod PARAM de paramètre de turbomachine ou d'aéronef (par exemple température, pression, en absolu ou en relatif) va être décrit, en référence à la représentation général de la figure 4. On parlera de « paramètre d'intérêt ». Le modèle est à nouveau un modèle thermodynamique. Le modèle décrit l'évolution du paramètre en fonction d'une ou plusieurs variables Var qui sont elles aussi en réalité des paramètres de turbomachine ou d'aéronef (par exemple température, pression, en absolu ou en relatif). 11 est stocké dans la mémoire 120 de l'unité de calcul 100. [0066] Ce procédé est pleinement applicable au procédé de recalage de la pression Ps3 décrit précédemment. On utilisera d'ailleurs la pression Ps3 comme exemple de paramètre PARAM et la pression PCN25R comme variable Var mais le procédé peut être appliqué à tout paramètre physique PARAM d'un aéronef et toute variable Var (par exemple la pression PT2) : par exemple mod_Ps3/P25(PCN25R), mod_Ps3/P25(PCN25R, PCN12R), mod_Ps3/P25(PCN25R, PT2), mod_T25(PCN12R, PT2), mod_Xn25(PCN12R,PT2) où Mach est la vitesse de l'avion, mod_T3(T25), etc. [0067] Un modèle est ici défini comme une loi par segment (dans une configuration dite 2D) ou par plan (dans une configuration dite 3D) indiquant la valeur dudit paramètre d'intérêt en fonction respectivement de une variable Var (2D) ou de deux variables Var l, Var2 (3D). La loi est linéaire respectivement sur chaque segment (ou autrement dit, affine par morceau : c'est-à-dire que son équation est sous la forme générique z=ax+c) ou sur chaque plan (équation sous la forme générique z=ax+by+c). [0068] L'intérêt d'un modèle défini comme une loi par segment (2D) ou par plan (3D) est 8 l'application des principes de l'automatique linéaire. Par exemple, le modèle mod Ps3/P2.5(Xn25r) ou mod Ps3/P25(PCN25R) est non linéaire dans son entièreté. [0069] On se place dans le même cadre que précédemment, avec les deux voies d'acquisition V10, V20. [0070] Dans une étape El, une valeur Val_PARAM du paramètre d'intérêt PARAM est obtenue. Cette obtention peut se faire dans le cadre de l'étape El décrite précédemment, par une mesure d'un capteur 10, 20 d'un ou de plusieurs voies d'acquisitions V10, V20, avec notamment l'acquisition d'un paramètre tiers et on en déduit ledit paramètre intérêt. Alternativement ou complémentairement, le paramètre d'intérêt PARAM peut être obtenu à l'aide d'une simulation. [0071] Les étapes et sous-étapes suivantes sont mises en oeuvre par le processeur 110 et la mémoire 120 de l'unité de calcul 100. [0072] Une étape E2 de conversion des données peut être mise en oeuvre lorsque le paramètre mesuré ne correspond pas au paramètre du modèle : par exemple, comme expliqué précédemment, Val Ps3 est une valeur de pression statique Ps3, alors que le modèle mod Ps3/P25 utilise la pression Ps3 réduite sur la P25. Dans le cas d'un paramètre tiers, ladite unité de calcul 100 calcule une valeur du paramètre d'intérêt Val_PARAM à partir de la valeur du paramètre tiers. [0073] Ensuite, l'étape de recalage E3 est mise en oeuvre. Cette étape de recalage E3 comprend plusieurs sous-étapes. [0074] Dans une sous-étape E31, le processeur 110 récupère la valeur Val_mod_PARAM du modèle mod_PARAM qui correspond à la valeur du paramètre d'intérêt Val_PARAM obtenue à l'étape El. [0075] La valeur du modèle Val_mod_PARAM se trouve donc sur un des segments ou des plans de du modèle mod_PARAM. Cette correspondance peut se faire via la valeur de la variable Var du modèle mod_PARAM : on prend la valeur du modèle Val_mod_PARAM dont l'abscisse correspond à celle de la valeur Val_PARAM du paramètre d'intérêt. Pour cela, il peut être nécessaire d'effectuer en fait deux mesures : une sur le paramètre PARAM et une sur la variable Var, pour avoir un couple de données. [0076] Dans le cas de la pression Ps3, on peut ainsi avoir une mesure de la PCN25R en même temps que la mesure de la Ps3. [0077] Avec les deux valeurs Val_mod_PARAM et Val_PARAM, la sous-étape E31 comprend le calcul d'une erreur E, typiquement par soustraction : e=Val_mod_PARAM-Val_PARAM. Cette erreur E est illustrée en figure 5. [0078] Dans une sous-étape E32, cette erreur E est traitée par un correcteur 122, dont le rôle est de minimiser ladite erreur E. Le correcteur 122 permet de calculer une correction con qui est un écart à appliquer sur les coordonnées des points de la loi corrigée, 9 obtenue via le correcteur PID, à partir de l'erreur (écart entre la mesure et le modèle) et qui doit être apportée au modèle mod PARAM. Du fait de la segmentation (segment ou plan) du modèle m PARÀM, le correcteur n'est implémenté que sur le segment ou le plan considéré lors de la mise en oeuvre de l'étape E3. [0079] Un correcteur particulier sera décrit par la suite. [0080] Enfin, dans une sous-étape E33, la correction corr est utilisée pour recaler le segment ou le plan du modèle mod PARAM. Cette étape consiste à recalculer un segment ou un plan, à partir du modèle mod_PARAM précédent et de la correction corr calculée à la sous-étape E32. En particulier, le recalage consiste à déplacer un nombre minimum de points du modèle mod_PARAM dans une sous-étape E331 et à interpoler le reste du modèle entre ces points dans une sous-étape E332 : deux points pour le modèle par segments et trois points pour le modèle par plan. [0081] Plusieurs modes de réalisation du recalage seront décrits par la suite. [0082] On remarque en outre, par exemple sur la figure 3, que le recalage d'un segment va influer aussi les segments adjacents dans le cas où l'extrémité du segment recalé est déplacée. Une étape d'interpolation des segments adjacents peut en outre être mise en oeuvre. [0083] Le correcteur choisi est un correcteur PID (proportionnel intégral dérivé), illustré en figure 5, où Gp, Gd et Gi sont respectivement le gain du correcteur proportionnel, du correcteur dérivé et du correcteur intégral, s'étant la variable dans le domaine fréquentiel. Le correcteur intégral (le I du PID) permet d'introduire une certaine inertie au système rebouclé, ce qui permet d'éviter une hyper-sensibilité aux perturbations et aux points fous, par rapport à un correcteur tout ou rien. Le correcteur intégral permet également de maîtriser la vitesse de recalage, et éviter une dérive instantanée du modèle m(param) vers la moyenne entre les deux voies VIO, V20 en cas de dérive d'un des capteurs 10, 20. Un correcteur proportionnel (le P du PID) et un correcteur dérivé (le D du PID) sont implémentés pour régler plus finement le correcteur 122 en cas de besoin mais ne sont pas utilisés (l'approche empirique a montré que leur contribution est marginale comparée à celle de l'intégrateur qui transcrit naturellement bien mieux le comportement souhaité pour le recalage). On peut ainsi avoir Gp=Gd=O. Le réglage du correcteur est fait de telle sorte que le modèle mod_PARAM soit recalé suffisamment rapidement pour rendre compte des reconfigurations de la turbomachine (par exemple un changement des niveaux de prélèvements d'air sur le compresseur haute pression). [0084] Modèle par segment (2D) [0085] On se place ici sur le segment du modèle mod_PARAM qui est concerné par la mesure Val_PARAM effectuée en étape E1. Ce segment a deux points extrémaux, à gauche et à droite, notés A et B. [0086] Recalage point par point [0087] La première solution, illustrée en figure 6a et 6b, consiste à rendre compte de la correction en modifiant les coordonnées d'un seul point du segment, par exemple un des points extrémaux A ou B. [0088] Dans ce cas de figure, la sortie du correcteur 122 irnpacte directement le point B (respectivement le point A), et le point A (respectivement le point B) reste figé. Cette solution contraint cependant à figer au moins un des points du modèle mod_PARAM pour servir de référence, à partir de laquelle seront impactés les autres segments du modèle mod PARAM. Ainsi, durant l'étape de recalage E2 et plus précisément durant la sous-étape E231, seul un des deux points extrémaux est déplacé. Ensuite, l'étape d'interpolation E232 est mise en oeuvre. [0089] Cette solution est la plus simple et la plus rapide à calculer. [0090] Recalage pondéré des deux points du segment [0091] La deuxième solution, illustrée en figures 7a et 7b, consiste à répartir la correction de manière pondérée pour permettre de recaler de manière plus représentative et plus efficace le segment sélectionné. Dans un mode de réalisation avantageux, la pondération est effectuée en fonction de la distance entre valeur Val PARAM, ici Val Ps3/P25, et les points A et B du segment. [0092] Les figures 7a et 7b illustrent le recalage sur un intervalle et un pas de calcul : - étape El : la valeur mesurée Val_PARAM est obtenue par une ou deux voies d'acquisition VIO, V20 ; dans l'exemple, il s'agit de Val_Ps3 , - étape E2 (image (a) de la figure 7b): la valeur mesurée Val_PARAM est convertie pour être homogène avec le modèle mod_PARAM ; par simplification, on garde la même référence Val_PARAM, - étape E31 (image (b) de la figure 7b) : on mesure E qui est l'écart entre la valeur mesurée Val_PARAM et valeur du modèle Val_mod_PARAM ; dans l'exemple avec la pression Ps3 : Val_PARAM=Val_PS3/P25, c'est-à-dire la pression Ps3 mesurée divisée par la pression P25 modèle et Val_mod_PARAM = Val_mod_Ps3/P25, la pression Ps3 du modèle recalé (par de précédentes itérations) que l'on divise par P25 modèle, - étape E32 (image (b) de la figure 7b) : l'erreur E est minimisée via le correcteur 122, en l'intégrant, pour calculer une correction coin, - étape E331 (figure 7a): on mesure ensuite (ou avant l'étape E31) la distance du point Val_mod_PARAM, ici Val_mod(Ps3/P25), au point A, qui constitue la borne inférieure de l'intervalle de la variable Var (ici PCN25R) et qui est fonction de la linéarisation du modèle choisie, pal- rapport à la distance entre les points A et B. Enfin, la 11 correction est distribuée sur l'ordonnée des points A (pour donner A') et B (pour donner B'), - étape E332 (image (c) de la figure 7b) : un nouveau segment est interpolé entre les deux points recalés A' et B'. [0093] Le principe de fonctionnement est de distribuer la correction torr du correcteur 122 d'un intervalle sur les ordonnées des points A et B selon le même principe que précédemment : dans un mode de réalisation, on distribue X% de la correction sur l'ordonnée du point B, avec X le rapport entre la distance du point Val mod PARAM au point A sur la distance du point A au point B. On distribue 100-X% de la correction sur l'ordonnée du point A (30% et 70% sur la figure 7a). [0094] Une fois les deux points A' et B' replacés, il suffit dans l'étape E232 d'interpoler le modèle entre ces deux points. La loi étant définie par segment, l'interpolation linéaire (ou affine) est simple. [0095] Alternativement, n'importe quels autres points (distincts) du segment peuvent être déplacés par la correction : il suffit de choisir deux points et l'interpolation linéaire (ou affine) permet de compléter le reste du segment considéré. [0096] Cette méthode permet ainsi un recalage efficace et rapide pour obtenir un modèle recalé [nad PARAM. Néanmoins, ce modèle mod_PARAM ne dépendant que d'une variable Var (PCN25R dans le cas de niod Ps3), il peut être insuffisant pour certaines situations de vol, notamment lorsque le paramètre d'intérêt PARAM dépend de plusieurs variables Varl, Var2. [0097] Modèle par plan (3D) [0098] A cet égard, pour tenir compte de plusieurs variables, le modèle mod_PARAM peut être fonction de deux variables (mod_PARAM(Varl, Var2)) et être exprimé sous la forme d'une loi définie par plans, la loi étant linéaire sur chaque plan comme illustré en figure 8. [0099] La figure 9 illustre l'implémentation d'une méthode de recalage dans le cas d'un modèle par plan. [0100] Par exemple, dans le cas de la pression Ps3, lors de l'activation d'un niveau de pré- lèvement d'air, le modèle mod_Ps3/P25(PCN25R) (c'est-à-dire le modèle Ps3 réduit sur P25 en fonction de PCN25R) est modifié car une partie de l'air compressé par le compresseur haute pression est envoyée dans le système d'air avion). Le correcteur 122 du modèle 2D par segment permet éventuellement de s'adapter à cette reconfiguration si les gains du correcteur 122 sont réglés de sorte que le recalage du modèle soit rapide, mais cela peut poser d'autres difficultés. [0101] Toujours dans l'exemple de la pression Ps3, pour s'affranchir du problème des pré- lèvements d'air, un modèle de Ps3/P25 qui dépend non plus seulement de PCN25R, mais de PCN12R également, est alors implémenté : on définit alors 12 mod Ps3/P25(PCN25R, PCN12R). Lors de l'activation des prélèvements, la loi reliant PCN25R et PCN12R est modifiée, ce qui nous permet de rendre compte de la reconfiguration du système. Le recalage de cette loi nécessite donc un nouveau correcteur « 3D ». [0102] Recalage point par point [0103] La première solution, non illustrée, consiste à rendre compte de la correction en fixant les coordonnées d'un seul point du rectangle, par exemple un des sommets A, B, C ou D du rectangle et en modifiant les coordonnées de deux points du rectangle, par exemple deux des sommets A, B, C ou D. Alternativement, on peut fixer deux points fixant les coordonnées de deux points du rectangle, par exemple deux des sommets A, B, C ou D du rectangle et en modifiant les coordonnées de un point du rectangle, par exemple deux des sommets A, B, C ou D. [0104] Les points concernés sont déplacés pendant la sous-étape E331 puis l'étape d'interpolation E332 sur l'ensemble du rectangle est mise en oeuvre. Comme on travaille sur trois points à chaque fois, on est assuré de l'existence du rectangle interpolé. [0105] Recalage pondéré [0106] Pour permettre un recalage pondéré, pour lequel aucun point n'est fixe, le modèle mod_PARAM est linéarisée en découpant le rectangle ABCD en triangles ABC, ABD, typiquement deux triangles complémentaires (figure 8). En effet, trois points A, B, C sont toujours coplanaires, avant et après recalage, ce qui assure de l'existence de l'interpolation du triangle recalé à la sous-étape d'interpolation E332, une fois la sous-étape E331 de recalage des trois points effectuée. Les trois nouveaux points issus de la correction peuvent ainsi être utilisés pour décrire l'équation cartésienne d'un plan, nous permettant ainsi d'interpoler linéairement le modèle mod_PARAM. [0107] En effet, si une correction pondérée sur trois points de la surface était appliquée sur quatre points, par exemple les quatre sommets ABCD du rectangle, il y aurait une déformation du rectangle si les quatre points du rectangle n'étaient plus coplanaires (impossible d'interpoler les coordonnées du paramètre PARAM grâce à l'équation cartésienne cl' un plan). [0108] Dans la sous-étape E331, il s'agit de d'abord de sélectionner le triangle à recaler en fonction de la valeur de Val_PARAM (appelé point X) obtenue par les étapes El et E2. Pour cela une différence de pente entre le segment AC qui divise le rectangle en deux et le segment AX (figure 101)). On peut utiliser n'importe quel sommet B, C ou D. [0109] En effet, avec les quatre points A, B, C, D formant un rectangle et le point X cor- respondant au point mesuré Val_PARAM, il faut déterminer si X appartient au triangle ABC ou au triangle ACD (on rappelle que ces triangles ont été choisis de manière arbitraire par rapport à ABD et DBC). 13 [0110] Pour cela, une comparaison des valeurs des taux de variation AAC, AAX des droites (AC) et (AX) est effectuée lors de la sous-étape E331. En effet si AAX > AAC alors ACD est sélectionné et si AAX REVENDICATIONS[Revendication 1] Procédé d'analyse de vieillissement d'une turbomachine (1), le procédé consistant à mettre en oeuvre les étapes suivantes : - FI : Effectuer un rccalage d'un modèle (mod_PARAM) d'un paramètre physique de la turbomachine (1), le modèle étant une loi définissant la valeur du paramètre en fonction d'au moins une variable, - F2 : Sauvegarder le modèle recalé (mod_PARAM) dans une mémoire non-volatile (120), les étapes FI et F2 étant répétées au moins deux fois, et préférablement plus, - F3 : Comparer les différents modèles recalés (mod_PARAM) pour en déduire une évolution de l'état de la turbomachine. [Revendication 2] Procédé selon la revendication 1, dans lequel le recalage des modèles se fait à l'aide d'une mesure d'une valeur du paramètre (Val_PARAM) par un capteur (10, 20) ou d'une mesure d'un paramètre tiers. [Revendication 3] Procédé selon la revendication 1 ou 2, dans lequel le recalage des modèles se fait à l'aide d'un correcteur PID, préférablement un correcteur uniquement intégral. [Revendication 4] Procédé selon l'une quelconque des revendications 1 à 3, dans lequel les modèles sont définis par segment ou par plan. [Revendication 5] Procédé selon l'une quelconque des revendications 1 à 4, dans lequel au moins deux modèles recalés (mod_PARAM) ont été générés à plus d'un mois d'intervalle. [Revendication 6] Procédé selon l'une des revendications 1 à 5, dans lequel le modèle recalé dans l'étape F1 est un modèle de pression statique en amont de chambre de combustion (mod_Ps3(PCN25R)), dit « modèle de Ps3». dans une turbomachine comprenant un compresseur (3), le modèle de Ps3 étant utilisé pour arbitrer entre deux voies d'acquisition (V10. V20) de la pression statique en amont de la chambre de combustion (Ps3), dite « pression Ps3», les deux voies d'acquisition (V10, V20) faisant intervenir deux capteurs (10, 20), le procédé utilisant un modèle de Ps3 stocké dans la mémoire (120), le modèle exprimant la pression Ps3 en fonction au moins du régime (PCN25R) du compresseur (3), dit «régime PCN25R » et comprenant les étapes suivantes : El : mesure d'une valeur de la pression Ps3, par un des deux capteurs (10, 20), 17 [Revendication 7] [Revendication 8] [Revendication 9] [Revendication 10] E2 : recalage du modèle de Ps3 à l'aide de la mesure de la valeur de Ps3. Procédé selon la revendication 6, dans lequel le modèle de Ps3 est un modèle de Ps3 sur la pression de compresseur (3), dite « pression P25 », le modèle étant dit « modèle PS3/P25 » (mod PS3/P25). Procédé selon la revendication 7, dans lequel le modèle Ps3/P25 est exprimée en fonction du régime de compresseur, réduit sur sa température, dite «température T25 », dit « régime PCN25R » (mod PS3/P25(PCN25R) ou « régime Xn25R ». Procédé selon l'une des revendications 1 à 5, dans lequel le modèle recalé dans l'étape FI est un modèle de paramètre de fonctionnement (mod_PARAM) de turbomachine (1) ou d'aéronef_ le modèle étant défini comme une loi par segment indiquant la valeur dudit paramètre en fonction d'une variable (mod PARAM(Var)), ou étant défini comme une loi par plan indiquant la valeur dudit paramètre en fonction de deux variables (mod PARAM(Varl , Var2)_ ladite loi étant affine sur chaque segment ou étant affine sur chaque plan, le modèle de paramètre étant stocké dans la mémoire (120), le procédé de recalage comprenant les étapes suivantes : - obtention d'une valeur du paramètre (Val_PARAM) (étape El ), - calcul d'une erreur (E) par comparaison de ladite valeur du paramètre (Val_PARAM) avec la valeur correspondante du modèle (Val_mod_PARAM), ladite valeur du modèle (Val_mod_PARAM) appartenant à un des segments ou des plans du modèle (mod_PARAM) (étape E31), - application d'un correcteur (112) en minimisant ladite erreur (E) pour déterminer une correction (con) (étape E32), - recalage du segment du modèle (mod_PARAM) ou du plan (lu modèle à l'aide de la correction (con-), pour repositionner ledit segment ou plan et obtenir ainsi un modèle recalé du paramètre physique (étape E33). Procédé selon la revendication 9, dans lequel l'étape d'obtention de la valeur du paramètre (Val_PARAM) se fait par : - une mesure directe dudit paramètre (PARAM) à l'aide d'un capteur (10, 20), ou - une mesure d'un paramètre tiers dont dépend ledit paramètre (PARAM), ou - une simulation.