Brevet FR3096031A1 — publié le 20201120

Description Titre de l'invention : Recalage de modèle par segment ou plan dans une turbomachine 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 direct drive) 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. [0006] 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 deux 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 le sens où l'avion est immobilisé tant que le moteur n'a pas été endoscopé pour vérifier 2 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 awhority digital engine 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. [0020] De plus, diverses améliorations ou utilisations du modèle thermodynamique 3 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é de recalage d'un modèle de paramètre de fonc- tionnement 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, ou étant défini comme une loi par plan indiquant la valeur dudit paramètre en fonction de deux variables de fonctionnement, ladite loi étant affine sur chaque segment ou étant affine sur chaque plan, le modèle de paramètre étant stocké dans une mémoire. [0023] Les paramètres de fonctionnement et les variables sont par exemple relatif à une tem- pérature ou une pression, ou encore à un régime de compresseur (typiquement les régimes Xn12 et Xn25 du corps basse pression et du corps haute pression. Plus généralement, ils peuvent être tout paramètre de fonctionnement pour lequel on dispose d'une mesure et d'un modèle permettant une redondance analytique [0024] Lc procédé de recalage comprend les étapes suivantes : - obtention d'une valeur du paramètre, - calcul d'une erreur par comparaison de ladite valeur du paramètre avec la valeur correspondante du modèle, ladite valeur du modèle appartenant à un des segments ou des plans du modèle, - application d'un correcteur en minimisant ladite erreur pour déterminer une correction, - recalage du segment du modèle ou du plan du modèle à l'aide de la correction, pour repositionner ledit segment ou plan et obtenir ainsi un modèle recalé du paramètre physique. [0025] Dans un mode de réalisation, l'étape d'obtention de la valeur du paramètre se fait par : - une mesure directe dudit paramètre à l'aide d'un capteur, ou - une mesure d'un paramètre tiers dont dépend ledit paramètre, ou - une simulation. [0026] Dans un mode de réalisation, le correcteur est un correcteur PID ou un correcteur intégral. [0027] Dans un mode de réalisation, lorsque le modèle est une loi par segment, le recalage se fait en figeant un point du segment et en déplaçant un autre point du segment à l'aide de la correction, les deux points étant préférablement des extrémités du segment. [0028] Dans un mode de réalisation, lorsque le modèle est une loi par segment, le recalage se fait en ne gardant aucun point du segment fixe, par exemple en déplaçant les deux extrémités du segment à l'aide de la correction. [0029] Dans un mode de réalisation, le déplacement des extrémités du segment se fait en fonction de leur distance respective avec ladite valeur correspondante du modèle. [0030] Dans un mode de réalisation, la répartition de la correction à appliquer une extrémité du segment est égal au ratio de la distance de la valeur correspondante du modèle à l'autre extrémité du segment, sur la longueur du segment. [0031] Dans un mode de réalisation, l'étape de recalage du segment du modèle comprend une interpolation linéaire entre deux points recalés. [0032] Dans un mode de réalisation, lorsque le modèle est une loi par plan, le plan a la forme d'un rectangle qui est découpé en triangles, et le recalage se fait en figeant un ou deux sommets du triangle et en déplaçant les deux derniers ou le dernier sommet du triangle à l'aide de la correction. [0033] Dans un mode de réalisation, lorsque le modèle est une loi par plan, le plan est découpé en triangles, et le recalage se fait en déplaçant les trois sommets du triangle. [0034] Dans un mode de réalisation, le déplacement de chaque sommet du triangle se fait en fonction de l'aire du sous-triangle défini par les deux autres sommets et ladite valeur correspondant du modèle. [0035] Dans un mode de réalisation, la répartition de la correction à appliquer à un sommet du triangle est égal au ratio de l'aire du sous-triangle défini par les autres sommets et ladite valeur correspondante du modèle, sur l'aire du triangle. [0036] Dans un mode de réalisation, l'étape de recalage du triangle comprend une inter- polation linéaire à partir des points recalés. [0037] Dans un mode de réalisation, le paramètre est la pression Ps3 ou la pression Ps3 divisée par la pression P25 et dans lequel : - la variable est, lorsque le modèle est une loi par segment, le régime PCN25R et - les variables sont, lorsque le modèle est une loi par plan, la PCN25R et la PCN12R, ou la PCN25R et la PT2. [0038] Dans un mode de réalisation, le modèle à recaler est choisi en fonction d'une variable, la mémoire stocke une pluralité de modèles exprimés en fonction du prélèvement d'air avion, la variable pouvant être le niveau de prélèvement d'air avion dans les compresseurs. [0039] Dans un mode de réalisation, les gains du correcteur sont différents pour différents segments ou plans du modèle. [0040] 11 est aussi proposé un procédé d'analyse de vieillissement d'une turbomachine, le procédé consistant à mettre en oeuvre les étapes suivantes : - Fl : Mise en oeuvre d'un procédé de recalage tel que décrit précédemment, - F2 : Sauvegarder le modèle recalé dans une mémoire non-volatile, les étapes Fl 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 une évolution de l'état de la turbomachine. Brève description des dessins [0041] D'autres caractéristiques, buts et avantages de l'invention ressortiront de la des- cription qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : [0042] [fie.1] La figure 1 illustre de façon schématique une turbomachine. [0043] [fie.2] La figure 2 illustre un procédé (l'arbitrage entre deux voies d'acquisition à l'aide d'un modèle thermodynamique. [0044] [fig.3] La figure 3 illustre graphiquement un procédé dc recalage dc la pression Ps3. [0045] [fig.4] La figure 4 illustre par un schéma-Hoc un procédé dc recalage d'un modèle de paramètre, tel que la pression Ps3. [0046] W8.51 La figure 5 illustre un correcteur. [0047] [fig.6a] [0048] [fig.611] Les figures 6a et 6h illustrent un procédé de recalage d'un modèle 2D par segment. [0049] [fig.74 La figure 7a illustre, pour un segment, un procédé dc rccalage d'un modèle 2D en segment par pondération. [0050] [fig.711] La figure 7h illustre, pour plusieurs segments, un procédé de recalage d'un modèle 2D en segment par pondération. [0051] [fig.8] La figure 8 illustre un modèle 3D par plan. [0052] [fig.9] La figure 9 illustre par un schéma-Hoc un procédé dc recalage d'un modèle 3D de paramètre, tel que la pression Ps3, en fonction des pressions PCN12R et PCN25R. [0053] [fig.104 La figure 10a illustre, pour un plan, un procédé dc recalage d'un modèle 3D en segment par pondération. [0054] [fig.1011] La figure 10h illustre le choix d'un triangle parmi le rectangle formant un plan du modèle 3D. [0055] [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. [0056] [fig.11] 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. [0057] [fig.12] La figure 12 illustre un procédé d'analyse de vieillissement de turbomachine. Description des modes de réalisation 100581 Le contexte et les définitions données en introduction sont repris ici. 100591 Dans un premier temps, un procédé de recalage de modèle de pression statique en 6 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. [0060] Le modèle de Ps3 a notamment pour but final de permettre d'arbitrer entre deux voies d'acquisition V10, V20, redondantes, dont la fonction est de mesurer la pression Ps3. Chaque voie d'acquisition V10, V20 comprend un capteur 10, 20. Le capteur 10, 20 est standard et ne sera pas décrit ici. [0061] On décrira par la suite un procédé d'arbitrage entre les deux voies d'acquisition V10, V20. [0062] 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 englue control ») ou bien être un composant distinct, positionné au plus près des voies d'acquisition V10, V20 pour davantage de réactivité. [0063] 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 (ce 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). [0064] Il suffit ensuite de multiplier la valeur de Ps3/P25 par P25 pour obtenir la valeur de Ps3. [0065] 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 mod_Ps3, inclut les modèles n'exprimant pas directement la pression Ps3 mais permettant de l'obtenir par la suite, tel que le modèle mod_Ps3/P25. [0066] Dans une première étape El, une des deux voies d'acquisition VIO, 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 VIO, 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. [0067] Cette mesure d'une valeur Val_Ps3 de la pression Ps3 est ensuite envoyée à l'unité de calcul 100. [0068] Une étape E2 de conversion ou de traitement des données peut être mise en oeuvre : 7 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. [0069] 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 unc pluralité, voire une infinité, si le modèle est continu) dont l'ordonnée a été déplacée (à abscisse constante donc). On note Rtnod PS3/P25 le modèle recalé. Par la suite, on simplifiera l'écriture en gardant mod PS3/P25 qui désigne un modèle avant et après recalage. [0070] 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 unc 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). [0071] 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. [0072] 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 de 1 'étape El du pas précédent sont mises en oeuvre. [0073] 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 VIO, V20 intempestif. [0074] Le recalage s'effectue avantageusement à l'aide d'un correcteur 112 qui est intégré dans une boucle de la chaine de contrôle. Ce correcteur sera décrit en détail par la suite. [0075] On définit aussi un procédé d'arbitrage entre deux voies d'acquisitions VIO, 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 VIO, V20 qui est la plus proche du modèle recalé. L'étape A2 est classique et ne sera pas décrite ici. [0076] 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, 8 pression, en ahsolu 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 ahsolu ou en relatif). Il est stocké dans la mémoire 120 de l'unité de calcul 100. [0077] 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 tnod Ps3/P25(PCN25R), mod Ps3/P25(PCN25R, PCN12R), tnod Ps3/P25(PCN25R, PT2), mod T25(PCN12R, PT2), mod Xn25(PCN12R, PT2) où Mach est la vitesse de mod T3(n5), etc. [0078] 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 Varl, 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). [0079] L'intérêt d'un modèle défini comme une loi par segment (2D) ou par plan (3D) est l'application des principes de l'automatique linéaire. Par exemple, le modèle mod_Ps3/P25(Xn25r) ou mod_Ps3/P25(PCN25R) est non linéaire dans son entièreté. [0080] On se place dans le même cadre que précédemment, avec les deux voies d'acquisition V10, V20. [0081] 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. [0082] 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. [0083] 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 9 Val_PARAM à partir de la valeur du paramètre tiers. [0084] Ensuite, l'étape de recalage E3 est mise en oeuvre. Cette étape de recalage E3 comprend plusieurs sous-étapes. [0085] Dans une sous-étape E31, le processeur 110 récupère la valeur Val [nad PARAM du modèle mod_PARAM qui correspond à la valeur du paramètre d'intérêt Val_PARAM obtenue à l'étape E1. [0086] 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. [0087] 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. [0088] Avec les deux valeurs Val mod PARAM et Val PARAM, la sous-étape E31 comprend le calcul d'une erreur E, typiquement par soustraction : mod PARAM-Val PARAM. Cette erreur F. est illustrée en figure 5. [0089] Dans une sous-étape E32, cette erreur F. 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 con-igée, obtenue via le correcteur MD, à 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_PARAM, 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. [0090] Un correcteur particulier sera décrit par la suite. [0091] Enfin, dans une sous-étape E33, la correction con 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 con 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. [0092] Plusieurs modes de réalisation du recalage seront décrits par la suite. [0093] 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. [0094] 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 (variable de Laplace). Le correcteur intégral (le I du PM) permet d'introduire une certaine inertie au système rehouclé, 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 V10, 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 P1D) 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=0. Le réglage du correcteur est fait de telle sorte que le modèle mod_PARAM soit recalé suffisamment rapidement pour rendre compte des rcconfigurations de la turbomachine (par exemple un changement des niveaux de prélèvements d'air sur le compresseur haute pression). [0095] Modèle par segment (2D) [0096] On se place ici sur le segment du modèle mod_PARAM qui est concerné par la mesure Val_PARAM effectuée en étape El. Ce segment a deux points extrémaux, à gauche et à droite, notés A et B. [0097] Recalage point par point [0098] La première solution, illustrée en figures 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. tandis que l'autre est figé. [0099] Dans ce cas de figure, la sortie du correcteur 122 impacte 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. [0100] Cette solution est la plus simple et la plus rapide à calculer. [0101] Recalage pondéré des deux points du segment [0102] La deuxième solution, illustrée en figures 7a et 7b, consiste à répartir la con-ection de manière pondérée pour permettre de recaler de manière plus représentative et plus 11 efficace le segment sélectionné. Dans un mode de réalisation avantageux, la pon- dé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. [0103] 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 V10, 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 (h) 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 (h) de la figure 7b) : l'erreur E est minimisée via le correcteur 122, en I' intégrant, pour calculer une correction torr, - é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 home 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 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'. [0104] Le principe de fonctionnement est de distribuer la correction con- 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). [0105] 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. [0106] 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é. [0107] Cette méthode permet ainsi un recalage efficace et rapide pour obtenir un modèle 12 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. [0108] Modèle par plan (3D) [0109] 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. [0110] La figure 9 illustre l'implémentation d'une méthode de recalage dans le cas d'un modèle par plan. [0111] 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. [0112] 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 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 ». [0113] Recalage point par point [0114] 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. [0115] 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é. [0116] Recalage pondéré 13 [0117] 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. [0118] 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 d'un plan). [0119] 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 10b). On peut utiliser n'importe quel sommet B, Cou D. [0120] 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). [0121] 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 I] Procédé de recalage d'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é dans une 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 du 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). [Revendication 2] Procédé selon la revendication 1, 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. [Revendication 3] Procédé selon la revendication 1 ou 2, dans lequel le correcteur (112) est un correcteur PID ou un correcteur intégral. [Revendication 4] Procédé selon l'une quelconque des revendications 1 à 3, dans lequel, lorsque le modèle (mod_PARAM(Var)) est une loi par segment, le recalage se fait en figeant un point du segment (A, B) et en déplaçant un autre point du segment (A, B) à l'aide de la correction (corr), les deux points (A, B) étant préférablement des extrémités du segment. [Revendication 5] Procédé selon l'une quelconque des revendications 1 à 4, dans lequel, lorsque le modèle est une loi par segment (mod PARAM(Var)), le 18 recalage se fait en ne gardant aucun point du segment fixe, par exemple en déplaçant les deux extrémités (A, B) du segment à l'aide de la correction (corr). [Revendication 6] Procédé selon la revendication 5, dans lequel le déplacement des ex- trémités du segment (A, B) se fait en fonction de leur distance respective avec ladite valeur correspondante du modèle (Val mod PARAM). [Revendication 7] Procédé selon la revendication 5 ou 6, dans lequel la répartition de la correction (corr) à appliquer une extrémité du segment (A, B) est égal au ratio de la distance de la valeur correspondante du modèle (Val mod PARAM) à l'autre extrémité du segment (B, A), sur la longueur du segment (AB). [Revendication 8] Procédé selon l'une quelconque des revendications 4 à 7 dans lequel l'étape de recalage du segment du modèle comprend une interpolation linéaire entre deux points recalés. [Revendication 9] Procédé selon l'une quelconque des revendications 1 à 4, dans lequel, lorsque le modèle (mod PARAM(Varl, Var2)) est une loi par plan, le plan a la forme d'un rectangle (ABCD) qui est découpé en triangles (ABC, ABD), et le recalage se fait en figeant un ou deux sommets du triangle et en déplaçant les deux derniers ou le dernier sommet du triangle à l'aide de la correction (corr). [Revendication 10] Procédé selon l'une quelconque des revendications 1 à 4, dans lequel, lorsque le modèle (mod_PARAM(Varl, Var2)) est une loi par plan, le plan est découpé en triangles, et le recalage se fait en déplaçant les trois sommets du triangle (A, B, C). [Revendication 11] Procédé selon la revendication 10, dans lequel le déplacement de chaque sommet (A, B, C) du triangle se fait en fonction de l'aire du sous-triangle (XBC, XAC, XAB) défini par les deux autres sommets (BC, AC, AB) et ladite valeur correspondant du modèle (Val_mod_PARAM). [Revendication 12] Procédé selon la revendication 11, dans lequel la répartition de la correction (con) à appliquer à un sommet du triangle est égal au ratio de l'aire du sous-triangle défini par les autres sommets et ladite valeur correspondante du modèle, sur l'aire du triangle (ABC). [Revendication 13] Procédé selon l'une quelconque des revendications 9 à 12, dans lequel l'étape de recalage du triangle comprend une interpolation linéaire à partir des points recalés (ABC). [Revendication 14] Procédé selon l'une quelconque des revendications 1 à 13, dans lequel le paramètre est la pression Ps3 ou la pression Ps3 divisée par la pression P25 (Ps3/P25) et clans lequel 19 [Revendication 15] [Revendication 16] [Revendication 17] - la variable (Var) est, lorsque le modèle est une loi par segment (mod Ps3/P25(Var), le régime PCN25R et - les variables (Varl, Var2) sont, lorsque le modèle est une loi par plan (mod Ps3 (Varl, Var2)), la PCN25R et la PCN12R, ou la PCN25R et la PT2. Procédé selon l'une quelconque des revendications 1 à 14, dans lequel le modèle (mod_PARAM) à recaler est choisi en fonction d'une variable, la mémoire (120) stocke une pluralité de modèles (mod_PARAM) exprimés en fonction du prélèvement d'air avion, la variable pouvant être le niveau de prélèvement d'air avion dans les compresseurs. Procédé selon l'une quelconque des revendications 1 à 15, dans lequel les gains du correcteur sont différents pour différents segments ou plans du modèle (mod PARAM). Procédé d'analyse de vieillissement d'une turbomachine (1), le procédé consistant à mettre en oeuvres les étapes suivantes : - Fl : Mise en oeuvre d'un procédé de recalagc selon l'une quelconque des revendications 1 à 16, - F2 : Sauvegarder le modèle recalé (mod_PARAM) dans une mémoire non-volatile (120), les étapes Fl 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.