FR 3016191 A1 20150710 FR 1450114 A 20140108 PROCÉDÉ ET PROGRAMME D'ORDINATEUR POUR LA SURVEILLANCE D'UN INVERSEUR DE POUSSÉE À ACTIONNEURS HYDRAULIQUES DESCRIPTION DOMAINE TECHNIQUE L'invention se rapporte au domaine de la surveillance d'un inverseur de poussée de moteur d'aéronef. En particulier, l'invention concerne un procédé et un programme d'ordinateur de surveillance d'un inverseur de poussée pour en détecter et en pronostiquer des anomalies. ÉTAT DE LA TECHNIQUE ANTÉRIEURE De façon connue, une opération de maintenance préventive consiste à effectuer une inspection d'une machine, par exemple un inverseur de poussée pour moteur d'aéronef, avant qu'une panne empêchant son fonctionnement ne se déclare, cette inspection pouvant conduire éventuellement au remplacement d'une ou de plusieurs pièces de la machine. Bien entendu, pour qu'elle soit efficace d'un point de vue préventif et ne pas engendrer de coûts inutiles, une telle opération ne doit pas être réalisée de façon intempestive mais à un moment opportun de la vie de la machine. On cherche donc d'une manière générale à réaliser une surveillance efficace d'un inverseur de poussée à même de délivrer un avis de maintenance fiable permettant ainsi d'optimiser les opérations de maintenance sur cet inverseur. On connaît ainsi de la demande de brevet FR 2 943 732 Al un procédé de surveillance d'un inverseur de poussée à actionneur commandé par un moteur électrique, procédé qui utilise différents types de paramètres acquis pendant une phase de fonctionnement de l'inverseur de poussée : - des paramètres représentatifs de la durée de cette phase, - des paramètres représentatifs de l'énergie vue par le moteur électrique, - des paramètres représentatifs du couple vu par le moteur électrique, - des paramètres représentatifs de la durée pendant laquelle la vitesse de rotation de consigne du moteur électrique diffère de la vitesse de rotation mesurée. Ce procédé s'applique à un inverseur de poussée électrique, n'incluant pas de commandes hydrauliques. Ainsi, pour un inverseur de poussée à actionneurs hydrauliques, seuls les paramètres représentatifs de la durée de cette phase de fonctionnement pourraient être utilisés. Toutefois, de tels indicateurs de durée peuvent être fortement influencés par le contexte extérieur du moteur (température extérieure en particulier). Or les interactions entre données endogènes et données exogènes ne sont pas bien connues. Il en découle que ce procédé connu ne permet pas de réaliser de manière satisfaisante la surveillance d'un inverseur de poussée à actionneurs hydrauliques. EXPOSÉ DE L'INVENTION L'invention vise à proposer une surveillance efficace d'un inverseur de poussée à actionneurs hydrauliques et propose pour ce faire, un procédé de surveillance d'un inverseur de poussée à portes escamotables de moteur d'aéronef, le moteur comprenant un calculateur configuré pour réaliser des mesures d'un paramètre représentatif de la position de contacteurs agencés pour renvoyer une information sur la position des portes, caractérisé en ce qu'il comprend un calcul d'un ou plusieurs indicateurs statistiques du paramètre mesuré et une analyse de l'évolution temporelle du ou des indicateurs statistiques calculés. Certains aspects préférés mais non limitatifs de procédé sont les suivants : - le calcul d'un ou plusieurs indicateurs statistiques du paramètre mesuré comprend la détermination d'une moyenne des mesures du paramètre et/ou d'une médiane des mesures du paramètre et/ou d'un écart-type des mesures du paramètre et/ou d'un type de distribution des mesures du paramètre ; - l'analyse de l'évolution temporelle d'un indicateur statistique calculé comprend le calcul d'une pente de l'indicateur statistique sur une fenêtre temporelle, ou encore la comparaison de l'indicateur statistique calculé avec un indicateur de référence préalablement calculé à partir de mesures du paramètre réalisées lors d'un fonctionnement considéré comme sain de l'inverseur de poussée ; - le résultat de l'analyse de l'évolution temporelle d'un indicateur statistique calculé est comparé à un seuil, et une alerte de détection d'anomalie est émise en cas de dépassement du seuil ; - il comprend, en cas de dépassement du seuil, le calcul d'une signature d'anomalie détectée à partir du ou des indicateurs statistiques calculés et la comparaison de la signature d'anomalie détectée à une ou plusieurs signatures de référence chacune caractéristique d'un type de défaillance ; - il comprend le calcul, à partir du résultat de l'analyse de l'évolution temporelle du ou des indicateurs statistiques calculés, d'une probabilité d'occurrence d'une anomalie ; - le paramètre mesuré par le calculateur est un ratio de tension entre une tension divisée par un circuit diviseur de tension et une tension d'alimentation du circuit diviseur de tension. L'invention porte également sur un produit programme d'ordinateur comprenant des instructions de code pour l'exécution des étapes du procédé selon l'invention, lorsque ledit programme est exécuté sur un ordinateur BRÈVE DESCRIPTION DES DESSINS D'autres aspects, buts, avantages et caractéristiques de l'invention apparaîtront mieux à la lecture de la description détaillée suivante de formes de réalisation préférées de celle-ci, donnée à titre d'exemple non limitatif, et faite en référence aux dessins annexés sur lesquels : - la figure 1 illustre un circuit diviseur de tension permettant de réaliser la mesure d'un paramètre représentatif de la position de contacteurs agencés chacun pour renvoyer une information sur la position d'une porte de l'inverseur de poussée ; - la figure 2 est un schéma illustrant un exemple de liaison entre les différents contacteurs d'un inverseur de poussée et le calculateur moteur au moyen de circuits diviseur de tension du type de celui de la figure 1 ; - la figure 3 est un organigramme illustrant les principales étapes du procédé de surveillance selon l'invention. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Un moteur d'aéronef comprend un calculateur de régulation connu sous le nom de FADEC (« Full Authority Digital Engine Control ») qui agit en temps réel sur différents actionneurs par l'analyse et l'exploitation continue de données collectées par de multiples capteurs tels que des contacteurs. Un inverseur de poussée pour un tel moteur comprend de manière connue une pluralité de portes escamotables qui sont déployées lors du freinage et fermées lors des autres phases de vol. Chaque porte peut être équipée de trois contacteurs à deux états, ouvert ou fermé, deux contacteurs étant agencées de manière à renvoyer une information sur le niveau de fermeture de la porte et un contacteur étant agencé de manière à renvoyer une information sur le niveau d'ouverture de la porte. Afin d'éviter toute ouverture intempestive des portes, plusieurs lignes de sécurité sont prévues. Prenant l'exemple d'un inverseur à deux portes et à actionneurs hydrauliques, on peut alors retrouver : Une première ligne de sécurité constituée de crochets en forme de S (un sur chaque porte) dont la position est détectée par des contacteurs d'escamotage (« stow switch ») au nombre de deux par crochets. Les crochets sont actionnés par un verrou hydraulique primaire (HPL pour « Hydraulic Primary Lock ») et un verrou électro-hydraulique primaire (EHPL pour « Electro-Hydraulic Primary Lock ») ; Une deuxième ligne de sécurité constituée d'actionneurs hydrauliques PDA (pour « Pivot Door Actuator ») et où l'on retrouve des contacteurs de déploiement (« deploy switch ») pour détecter la position des portes (un contacteur par porte) ; Une troisième ligne de sécurité comprenant un verrou tertiaire à pilotage électrique indépendant des autres contrôles électriques. Les informations des quatre contacteurs d'escamotage et des deux contacteurs de déploiement sont relayées vers le calculateur FADEC afin que celui-ci puisse déterminer la position exacte des portes de l'inverseur de poussée. Afin de réduire la quantité de câbles d'interconnexion entre le calculateur et les contacteurs positionnés au niveau des portes de l'inverseur de poussée, et ainsi de réduire le poids du système global, on a recours à un ou plusieurs circuits de détection des positions individuelles de contacteurs. Ces circuits agissent comme des multiplexeurs pour chacun fournir au calculateur sur une même voie de transmission une information sur l'état de l'ensemble des contacteurs dont ils surveillent la position. Un tel circuit, dont on pourra trouver une description détaillée dans la demande de brevet WO 2011/050289 Al, prend la forme d'un circuit de résistances formant un circuit diviseur de tension. Comme représenté sur la figure 1, un tel circuit diviseur de tension 1 peut comprendre une résistance de sortie Rlect en série avec des modules Ml, M2 reliés en parallèle entre eux, chaque module comprenant un contacteur Sa, Sb relié en série avec une première résistance de module Rai, Rbl ou une seconde résistance de module Rat, Rb2 en fonction de sa position, ouverte ou fermée. En référence à la figure 2, on notera que l'on peut notamment retrouver trois circuits la, lb, lc du type de celui illustré sur la figure 1 pour relayer les informations des différents contacteurs vers le calculateur FADEC (dont une seule voie est ici représentée): un premier circuit la associé aux contacteurs d'escamotage 3a, 4a de la première porte, un deuxième circuit lb associé aux contacteurs de déploiement 5a, 5b de la première et de la deuxième porte, et un troisième circuit associé aux contacteurs d'escamotage 3b, 4b de la deuxième porte. La valeur du ratio de tension entre la tension Ulect aux bornes de la résistance de sortie Rlect et la tension d'alimentation Uexc du circuit diviseur de tension 1 de la figure 1 varie en fonction des positions des contacteurs Sa et Sb. Ces ratios sont connus avec une précision estimée de 1%, et ont des valeurs caractéristiques lorsque l'inverseur est ouvert ou fermé (en phase stabilisée et non transitoire). Ils permettent en outre de donner des statuts de validité aux contacteurs Sa et Sb. Ainsi en mesurant un paramètre représentatif de la position de contacteurs agencés chacun pour renvoyer une information sur la position d'une porte, comme par exemple ce ratio de tension Ulect/Uexc du circuit diviseur de tension de la figure 1, le calculateur FADEC peut déterminer la position exacte des portes de l'inverseur de poussée. Dans le cadre de l'invention, et en référence à la figure 3, il est proposé de réaliser la surveillance d'un inverseur de poussée à portes escamotables de moteur d'aéronef sur la base de telles mesures, réalisées au cours d'une étape El, d'un paramètre représentatif de la position de contacteurs agencés chacun pour renvoyer une information sur la position d'une porte de l'inverseur de poussée, en particulier sur la base des mesures du ratio de tension présenté ci-dessus sans que cela ne soit pour autant limitatif. L'invention porte ainsi sur un procédé de surveillance d'un inverseur de poussée, en particulier un inverseur de poussée à actionneurs hydrauliques, par exemple un inverseur composé de deux portes escamotables. Ce procédé comporte une étape E2 selon laquelle on procède au calcul d'un ou plusieurs indicateurs statistiques du paramètre représentatif mesuré au cours de l'étape El, et une étape E3 selon laquelle on procède à une analyse de l'évolution temporelle du ou des indicateurs statistiques calculés. Les mesures du paramètre représentatif peuvent être acquises au cours d'une ou plusieurs phases de fonctionnement de l'inverseur parmi une phase transitoire d'ouverture, une phase d'ouverture stabilisée, une phase transitoire de fermeture et une phase de fermeture stabilisée. D'une manière générale, la phase transitoire d'ouverture débute à la commande d'ouverture et termine à l'expiration d'une première temporisation. S'en suit la phase d'ouverture stabilisée, puis la phase transitoire de fermeture qui débute à la commande de fermeture et termine à l'expiration d'une seconde temporisation, puis la phase de fermeture stabilisée. Le calcul d'indicateurs statistiques et l'analyse de leur évolution temporelle est spécifique à une phase de fonctionnement. Ainsi, si les quatre phases mentionnées ci-dessus sont surveillées, on procède au calcul d'indicateurs statistiques du paramètre mesuré représentatifs de chacune des phases, et à l'analyse de l'évolution temporelle de ces indicateurs statistiques représentatifs pendant chacune des phases. Selon le procédé de surveillance de l'invention, le calcul d'un ou plusieurs indicateurs statistiques du paramètre mesuré représentatif d'une phase de fonctionnement surveillée de l'inverseur de poussée peut comprendre la détermination d'une moyenne des mesures du paramètre et/ou d'une médiane des mesures du paramètre et/ou un écart-type des mesures du paramètre et/ou d'un type de distribution des mesures du paramètre (on peut ainsi vérifier que les mesures présentent une distribution gaussienne, par exemple par mise en oeuvre du test du chi2), ces mesures étant réalisées au cours de la phase de fonctionnement surveillée. L'analyse de l'évolution temporelle d'un indicateur statistique calculé peut quant à elle comprendre le calcul, typiquement par régression linéaire ou polynomiale, d'une pente de l'indicateur statistique sur une fenêtre temporelle, par exemple une fenêtre de durée paramétrable. Cette analyse peut également comprendre la comparaison de l'indicateur statistique calculé avec un indicateur de référence préalablement calculé à partir de mesures du paramètre réalisées lors d'un fonctionnement considéré comme sain de l'inverseur de poussée. Cette comparaison peut prendre la forme d'un calcul de score standardisé tel qu'un Z-score. Cette comparaison peut également être réalisée pour plusieurs indicateurs statistiques calculés de manière à fournir un score global. L'intérêt du calcul de score est de caractériser l'état de l'inverseur de poussée surveillé, sous forme de distance par rapport à un fonctionnement sain. Cette analyse peut également comprendre le calcul de probabilités conditionnelles d'occurrence d'anomalies, par exemple en parcourant un réseau bayésien à partir des indicateurs calculés, un tel réseau présentant l'avantage de permettre l'introduction de connaissances a priori d'experts dans le graphe en compléments de l'expérience acquise grâce à l'observation des indicateurs calculés. Le résultat de l'analyse de l'évolution temporelle d'un indicateur statistique calculé peut être comparé à un seuil, et une alerte de détection d'anomalie peut être émise en cas de dépassement du seuil comme représenté par l'étape E4 de la figure 3. A titre d'exemple illustratif, la pente de l'indicateur statistique sur une fenêtre temporelle peut être comparée à une pente seuil, une alerte étant émise en cas de dépassement de seuil. Egalement à titre d'exemple illustratif, l'écart entre un indicateur statistique mesuré et son indicateur de référence peut être comparé à un seuil, une alerte étant émise lorsque cet écart dépasse le seuil. Encore à titre d'exemple illustratif, une alerte de détection d'anomalie peut également être émise lorsqu'une probabilité conditionnelle d'occurrence d'une anomalie dépasse un certain seuil. Dans un mode de réalisation possible, le procédé selon l'invention comprend une étape de diagnostic E5 visant à classifier l'anomalie détectée en cas de dépassement de seuil et à localiser la panne correspondante. Cette étape E5 peut comprendre le calcul d'une signature d'anomalie détectée à partir du ou des indicateurs statistiques calculés et la comparaison de la signature d'anomalie détectée ainsi calculée à une ou plusieurs signatures de référence chacune caractéristique d'un type de défaillance. Les signatures de référence sont typiquement enregistrées dans une base de données. Les signatures de référence peuvent comprendre des signatures construites a priori par expertise, conformément à des comportements caricaturaux des indicateurs en cas d'anomalie. Les signatures de référence peuvent également être construites a posteriori une fois une anomalie détectée et le ou les composants fautifs identifiés. Dans un mode de réalisation possible, le procédé selon l'invention comprend une étape de prédiction E6 visant à pronostiquer l'occurrence future d'une anomalie. Cette étape E6 peut notamment comprendre le calcul d'une probabilité d'occurrence a posteriori d'une anomalie en fonction d'une probabilité d'occurrence a priori pour chaque anomalie définie par expertise et les écarts entre une signature de l'évolution courante déterminée à partir du ou des indicateurs statistiques calculés et des signatures de référence mentionnées ci-dessus chacune caractéristique d'une défaillance. Cette probabilité d'occurrence a posteriori permet ainsi de calculer une échéance avant de rencontrer une anomalie donnée (l'échéance se traduisant en une durée de fonctionnement restante ou en un nombre restant de cycles d'actionnement de l'inverseur avant de devoir réaliser une opération de maintenance), ou encore de calculer la probabilité de rencontrer une anomalie donnée sur un horizon temporel donné. On relèvera que d'une manière générale, les étapes E3-E6 du procédé selon l'invention peuvent être mises en oeuvre selon les techniques d'analyse détaillées dans le brevet FR 2939924 B1. Et selon une implémentation préférée, les différentes étapes du procédé selon l'invention sont exécutées au moyen d'instructions de code de programme. En conséquence, l'invention vise aussi un produit programme d'ordinateur comprenant des instructions de code pour l'exécution des étapes du procédé selon l'invention tel que décrit ci-dessus, lorsque ledit programme est exécuté sur un ordinateur. En particulier, le procédé peut être mis en oeuvre aussi bien avec un programme d'ordinateur embarqué, par exemple au sein du FADEC, qu'avec un programme d'ordinateur au sol exploitant les mesures du paramètre représentatif, typiquement les mesures du ratio de tension d'un circuit diviseur de tension, acquises et enregistrées en embarqué, par exemple par le FADEC. REVENDICATIONS1. Procédé de surveillance d'un inverseur de poussée à portes escamotables de moteur d'aéronef, le moteur comprenant un calculateur (3) configuré pour réaliser des mesures (El) d'un paramètre représentatif de la position de contacteurs (3a, 4a, 5a, 3B, 4b, 5b, Sa, Sb) agencés pour renvoyer une information sur la position des portes, caractérisé en ce qu'il comprend un calcul (E2) d'un ou plusieurs indicateurs statistiques du paramètre mesuré et une analyse (E3) de l'évolution temporelle du ou des indicateurs statistiques calculés. 2. Procédé selon la revendication 1, dans lequel le calcul (E2) d'un ou plusieurs indicateurs statistiques du paramètre mesuré comprend la détermination d'une moyenne des mesures du paramètre et/ou d'une médiane des mesures du paramètre et/ou d'un écart-type des mesures du paramètre et/ou d'un type de distribution des mesures du paramètre. 3. Procédé selon l'une des revendications 1 et 2, dans lequel l'analyse (E3) de l'évolution temporelle d'un indicateur statistique calculé comprend le calcul d'une pente de l'indicateur statistique sur une fenêtre temporelle. 4. Procédé selon l'une des revendications 1 et 2, dans lequel l'analyse (E3) de l'évolution temporelle d'un indicateur statistique calculé comprend la comparaison de l'indicateur statistique calculé avec un indicateur de référence préalablement calculé à partir de mesures du paramètre réalisées lors d'un fonctionnement considéré comme sain de l'inverseur de poussée. 5. Procédé selon l'une des revendications 1 à 4, dans lequel le résultat de l'analyse (E3) de l'évolution temporelle d'un indicateur statistique calculé est comparé à un seuil,et dans lequel une alerte de détection d'anomalie (E4) est émise en cas de dépassement du seuil. 6. Procédé selon la revendication 5, comprenant, en cas de dépassement du seuil, le calcul (E5) d'une signature d'anomalie détectée à partir du ou des indicateurs statistiques calculés et la comparaison de la signature d'anomalie détectée à une ou plusieurs signatures de référence chacune caractéristique d'un type de défaillance. 7. Procédé selon l'une des revendications 1 à 4, comprenant le calcul (E6), à partir du résultat de l'analyse de l'évolution temporelle du ou des indicateurs statistiques calculés, d'une probabilité d'occurrence d'une anomalie. 8. Procédé selon l'une des revendications 1 à 7, dans lequel le paramètre mesuré par le calculateur est un ratio de tension entre une tension divisée par un circuit diviseur de tension (1) et une tension d'alimentation du circuit diviseur de tension, le circuit diviseur de tension comprenant une résistance de sortie (Rlect) en série avec des modules (M1, M2) reliés en parallèle entre eux, chaque module comprenant un desdits contacteurs (Sa, Sb) relié en série avec une première (Ra1, Rb1) ou une seconde (Ra2, Rb2) résistance de module en fonction de la position dudit contacteur (Sa, Sb). 9. Procédé selon l'une des revendications 1 à 8, dans lequel l'inverseur de poussée est un inverseur à actionneurs hydrauliques 10. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution des étapes du procédé selon l'une des revendications 1 à 9, lorsque ledit programme est exécuté sur un ordinateur.30