Brevet FR2905763A1 — publié le 20080314

Procédé et système de diagnostic d'un aéronef à partir de mesures effectuées sur l'aéronef. La présente invention est relative à un procédé et à un système de diagnostic, i.e. de détection de l'existence d'un défaut au moins, d'une partie au moins d'un aéronef à partir de mesures effectuées sur l'aéronef. Le domaine technique de l'invention est celui de la fabrication d'hélicoptères. La présente invention est notamment relative à un procédé et à un système de détection automatique de l'existence d'un disfonctionnement mécanique à bord d'un giravion, à partir de mesures vibratoires pratiquées sur le giravion pendant un vol. L'invention s'applique notamment au diagnostic d'un mécanisme de transmission inséré entre au moins un moteur et au moins un rotor d'un giravion. Un tel mécanisme comporte habituellement plusieurs arbres équipés d'engrenage (s) et permet de transmettre le couple du(des) moteur(s) au(x) rotor(s) et accessoire(s) du giravion. La surveillance ou détection de l'apparition d'un défaut au sein d'un mécanisme a fait l'objet de nombreuses recherches ; on peut citer le brevet US6301572 relatif à la surveillance de l'état de moteurs ou turbine. La présente invention s'applique notamment au diagnostic d'un mécanisme par analyse des vibrations provoquées, en partie au moins, par le fonctionnement du mécanisme. Pour pratiquer ces mesures, on equipe l'hélicoptère d'accéléromètres qui sont placés (fixés) sur le(s) moteur(s), sur le(s) carter(s) de boîte(s) de transmission, sur des paliers d'arbres, et/ou en d'autres points de la structure de l'hélicoptère. Pendant un vol, les signaux délivrés par ces capteurs peuvent être convertis en données et le cas échéant synchronisés (grâce à des signaux 2905763 2 délivrés par un capteur de rotation) et/ou moyennés , puis enregistrés à bord de l'hélicoptère. De retour au sol, les données enregistrées peuvent être collationnées et analysées. L'interprétation de ces données est complexe : 5 elle nécessite une intervention longue d'un expert. Le document A rapid helicopter drive train fault detection using neuro-fuzzy method (http://erf32.nlr.nl/abstracts/pdf/HUO2.pdf), Bang Tran et al, propose d'appliquer une méthode logique neurofloue à la détection de pannes de systèmes mécaniques à partir de signaux 10 vibratoires. Les outils connus d'analyse automatique de ces données pour diagnostiquer un défaut mécanique dans un mécanisme sont incomplets et imparfaits ; des défauts existants ne sont pas détectés, tandis que des indications injustifiées de défaut sont parfois générées par ces outils 15 d'analyse. Un objectif de l'invention est de proposer un procédé d'analyse de telles données, un programme d'analyse et un dispositif incluant ce programme, qui permettent d'établir rapidement un diagnostic fiable, i.e. maximisant le pourcentage de défauts effectifs détectés et minimisant le 20 pourcentage de défaut non avérés. Un objectif de l'invention est de proposer de tels procédés, programmes, et dispositifs ou systèmes, qui soient améliorés et/ou qui remédient, en partie au moins, aux lacunes ou inconvénients des procédés, programmes, et systèmes connus. 25 Selon un aspect de l'invention, il est proposé un procédé de diagnostic d'un aéronef par analyse d'au moins un spectre obtenu à partir de mesures pratiquées sur l'aéronef, en particulier à partir de mesures vibratoires, dans lequel on applique à l'entrée d'un module de calcul le(s) spectre(s) à analyser ainsi qu'un premier spectre de référence au moins, le 2905763 3 module de calcul déterminant une valeur (unique) de diagnostic en fonction du (des) spectre(s) à analyser, du (des) spectre(s) de référence et de paramètres déterminés à partir d'au moins deux spectres de référence incluant ledit premier spectre de référence. 5 Selon un autre aspect de l'invention, il est proposé un procédé de diagnostic d'un aéronef par analyse d'au moins un spectre obtenu à partir de mesures, en particulier à partir de mesures vibratoires, pratiquées sur l'aéronef, dans lequel : dans une phase préparatoire, on applique au moins deux spectres 10 de référence à l'entrée d'un module de calcul paramétré qui produit en sortie une valeur de diagnostic ; on modifie les paramètres du module de calcul pour obtenir en sortie une valeur de référence ; puis dans une phase de diagnostic, on applique un ou plusieurs 15 spectres à analyser ainsi qu'une partie (certain(s)) des spectres de référence à l'entrée du module de calcul paramétré qui produit en sortie la valeur de diagnostic ; et on diagnostique un défaut lorsque cette dernière valeur dépasse la valeur de référence. En d'autres termes, on détermine d'abord des paramètres communs 20 aux spectres de référence ; on remplace ensuite certain(s) des spectres de référence par un ou plusieurs spectre(s) à analyser - en conservant au moins un spectre de référence - ; puis on détermine une valeur de diagnostic du (des) spectre(s) analysé(s), en fonction des spectres de référence conservés et des paramètres communs. 25 On utilise pour cela un module de calcul comportant une pluralité d'opérateurs structurés en couches et/ou en réseau, incluant des opérateurs paramétriques d'appartenance et des opérateurs paramétriques de pondération des spectres présentés en entrée. 2905763 4 Selon un mode de réalisation, il est proposé un tel procédé dans lequel, dans la phase de diagnostic, on applique au moins un vecteur de données spectrales - ou données fréquentielles - à diagnostiquer et au moins un vecteur de données spectrales de référence à l'entrée du module 5 de calcul qui détermine, pour chaque vecteur en entrée, plusieurs vecteurs de degré d'appartenance respective des composantes du vecteur à plusieurs classes, et qui calcule une valeur (unique) de diagnostic en fonction des vecteurs de degré d'appartenance, et dans lequel on signale une défaillance de l'aéronef lorsque la valeur de diagnostic dépasse un 10 seuil û ou valeur de référence-. Chacun des vecteurs de données spectrales peut résulter d'une transformée de Fourier appliquée à des données de mesures résultant de l'échantillonnage de signaux délivrés par des microphones ou accéléromètres à une fréquence qui peut être au moins égale à 1 kHz 15 environ (notamment une fréquence de l'ordre de 5, 10, 50 ou 250 kHz), chacun des vecteurs de données fréquentielles pouvant comporter plusieurs milliers de coordonnées/composantes/valeurs. Certaines au moins des opérations de la phase préparatoire et/ou de la phase de diagnostic des procédés selon l'invention peuvent être mises en 20 oeuvre par une unité électronique de traitement de données, telle qu'un calculateur, fonctionnant sous la commande d'un programme. Ainsi, selon un autre aspect de l'invention, il est proposé un programme comportant un code fixé sur un support û tel qu'une mémoire û ou matérialisé par un signal, le code étant lisible et/ou exécutable par au 25 moins une unité de traitement de données - telle qu'un processeur embarqué ou embarquable sur un aéronef, pour diagnostiquer le fonctionnement (bon ou mauvais) d'un aéronef, le code comportant des segments de code pour effectuer respectivement les différentes opérations d'un procédé selon l'invention. Selon un autre aspect de l'invention, il est proposé un système de diagnostic d'un aéronef, qui comporte : 2905763 5 - un organe de lecture d'un support de données qui est agencé pour lire des donnés de mesures effectuées sur l'aéronef, - une base de données contenant des données de spectres de référence pour l'aéronef, 5 - un dispositif de transformation des données de mesure du domaine temporel au domaine fréquentiel, qui est connecté à l'organe de lecture pour en recevoir les données de mesures et pour délivrer en sortie des données spectrales à analyser, et - un organe de calcul qui est relié à la base de données et au 10 dispositif de transformation, et qui est programmé pour effectuer les opérations d'un procédé selon l'invention. Quelques heures de vol peuvent suffire à constituer une base de données caractéristiques d'un type déterminé d'aéronef ; cette base comportant une quantité restreinte de données peut facilement être 15 enregistrée dans une mémoire associée à un calculateur pour servir ensuite à estimer, après chaque vol, l'état de bon fonctionnement de tout ou partie des mécanismes d'un aéronef. D'autres aspects, caractéristiques, et avantages de l'invention apparaissent dans la description suivante, qui se réfère aux dessins annexés 20 et qui illustre, sans aucun caractère limitatif, des modes préférés de réalisation de l'invention. La figure 1 illustre schématiquement un système de diagnostic selon un premier mode de réalisation de l'invention et sa mise en oeuvre, dans la phase préparatoire (d'apprentissage), pour la détermination des paramètres 25 du module de calcul paramétré qu'il comporte. La figure 2 illustre schématiquement l'utilisation du système de diagnostic de la figure 1 pour diagnostiquer l'état d'un aéronef à partir d'un nouveau spectre. 2905763 6 La figure 3 illustre schématiquement un système de diagnostic selon un second mode de réalisation de l'invention et sa mise en oeuvre, dans la phase préparatoire (d'apprentissage), pour la détermination des paramètres du module de calcul paramétré qu'il comporte. 5 La figure 4 illustre schématiquement l'utilisation du système de diagnostic de la figure 3 pour diagnostiquer l'état d'un aéronef à partir d'un nouveau spectre. La figure 5 illustre schématiquement la structure d'un module de calcul d'un système selon l'invention. 10 La figure 6 illustre schématiquement cinq courbes d'appartenance qui sont respectivement caractéristiques de cinq modules/opérateurs d'appartenance faisant partie d'un module de calcul tel que celui illustré figures 7 et 8. La figure 7 illustre un autre exemple de structure d'un module de 15 calcul paramétré d'un système selon l'invention et son utilisation dans la phase préparatoire. La figure 8 illustre l'utilisation du module de calcul de la figure 7 pour diagnostiquer le bon/mauvais fonctionnement d'un aéronef à partir d'un spectre vibratoire. 20 La figure 9 est un graphe illustrant schématiquement un exemple de spectre issu de mesures vibratoires et la figure 10 illustre, avec la figure 9, les coordonnées de ce vecteur. La figure 11 est un graphe illustrant schématiquement les variations au cours du temps de la valeur de diagnostic calculée pour un hélicoptère à 25 partir de spectres issu de mesures vibratoires au cours de vols successifs. Selon un aspect de l'invention, on détermine si (au moins) un spectre issu de mesures sur un aéronef est - ou non -significatif du bon fonctionnement global de l'aéronef, ou du bon fonctionnement d'un 2905763 7 mécanisme particulier de l'aéronef, en comparant l' allure ou la forme de ce spectre avec l'allure ou forme commune à (au moins) deux autres spectres - dits spectres de référence. Plus particulièrement, on compare l'allure du spectre à analyser à 5 certain(s) seulement des spectres de référence, à l'aide de paramètres préalablement déterminés à partir de tous les spectres de référence. Dans une certaine mesure, on peut considérer que ces paramètres sont caractéristiques de l'allure commune à tous les spectres de référence. Comme il ressort de la description suivante, le nombre de ces paramètres 10 est réduit, le nombre d'opérations nécessaires à la comparaison des spectres est également réduit, et ces opérations sont simples. Les procédés, programmes et systèmes selon l'invention permettent ainsi un diagnostic rapide ; il est en outre apparu que les résultats obtenus sont étonnamment fiables. 15 En pratique, un aéronef déterminé est généralement équipé de capteurs à demeure, notamment d'accéléromètres. Pendant chaque vol, des signaux délivrés par chaque accéléromètre sont numérisés et les données de mesure correspondantes sont enregistrées ; le procédé et système selon l'invention permettent de 20 diagnostiquer, après un vol, à partir de ces données enregistrées, si l'état de fonctionnement d'une partie au moins de l'aéronef est satisfaisant, par rapport à des mesures précédemment enregistrées et représentatives d'un état de fonctionnement satisfaisant du même aéronef. Pour que les résultats de cette comparaison puissent être 25 valablement utilisés aux fins d'un diagnostic, il importe non seulement que l'aéronef considéré soit le même, mais également que les capteurs à l'origine des mesures dont les spectres sont analysés, soient ceux ayant servi à l'élaboration des spectres de référence. 2905763 8 En pratique, il importe donc de constituer, au fur et à mesure de l'utilisation de l'aéronef, une base de données de spectres représentatifs du bon fonctionnement des différentes parties de l'aéronef et d'utiliser ensuite ces spectres comme spectres de référence. 5 Sauf indication explicite ou implicite contraire, dans la présente demande, les expressions donnée(s) spectrale(s) et donnée(s) fréquentielle(s) sont employées indifféremment pour désigner les valeurs des composantes d'un spectre obtenu à partir de mesures, par transformation de Fourier par exemple. 10 Sauf indication explicite ou implicite contraire, dans la présente demande, les expressions vecteur(s) de données spectrales et spectre(s) sont utilisées indifféremment pour désigner une suite finie de données spectrales qui est appliquée en entrée des modules de calcul des systèmes selon l'invention. 15 A titre d'exemple, on a illustré figure 9 un courbe 20 représentative d'un spectre obtenu à partir de mesures ; sur ce graphe, l'axe des abscisses représente la fréquence (f), et l'axe des ordonnées représente un niveau d'accélération qui peut être exprimé en décibel (dB), en mètre par seconde par seconde (m/s'), ou en g (gz9,81. m/s'). 20 Cette courbe est assimilée à une suite de points d'abscisses respectives f1 à f8 et d'ordonnées respectives v1 à v8. La figure 10 montre un vecteur Y de données spectrales [v1, v2, v3, v4, v5, v6, v7, v8] correspondant aux huit points en question de la courbe 20 figure 9. 25 En pratique, on utilise des vecteurs comportant un nombre plus élevé de valeurs spectrales ; le nombre de ces valeurs, qui dépend généralement de la valeur de la fréquence d'échantillonnage, peut par exemple être de l'ordre de 1000 valeurs - ou composantes du vecteur - ou de l'ordre de 10 000 valeurs. 2905763 9 Par référence aux figures 1 et 2, on utilise un module 23 paramétré de calcul qui délivre en sortie une valeur DIAG de diagnostic ; par référence à la figure 1, cette valeur résulte de l'application aux deux entrées du module 23, de deux spectres VR1 et VR2 de référence issus de 5 deux mémoires 21 et 22 (i.e. d'une base de données). Dans une phase préparatoire de détermination des paramètres de calcul du module 23, les deux spectres VR1 et VR2 sont présentés une seule fois en entrée du module, et le module fournit en sortie la valeur DIAG calculée, en fonction des valeurs initiales des paramètres du module 23, à 10 partir des deux vecteurs VR1 et VR2. Un additionneur 25 compare la valeur DIAG obtenue à une valeur Valref de référence enregistrée dans une mémoire 24, et délivre en sortie une valeur CorParam de correction des paramètres du module 23 ; cette valeur de correction est utilisée pour modifier les valeurs des paramètres 15 du module ; un nouveau calcul de la valeur DIAG est ensuite effectué à partir des mêmes spectres VR1 et VR2. Par itérations et modifications successives, on détermine ainsi les valeurs des paramètres du module 23 qui sont telles que la valeur DIAG obtenue en sortie est égale à la valeur Valref de référence. 20 Les valeurs ainsi obtenues pour les paramètres du module 23 sont alors figées (enregistrées) et le module 23 peut ensuite être utilisé pour établir un diagnostic à partir d'un spectre Vtest à analyser. A cet effet, comme illustré figure 2, le spectre à analyser est enregistré dans une mémoire 26 alimentée par des données issues du 25 capteur ayant auparavant délivré les données correspondant aux spectres VR1 et VR2 . 2905763 10 Le spectre de référence VR1 seulement et le spectre Vtest à analyser sont appliqués aux entrées du module 23 qui fournit en sortie une valeur DIAG de diagnostic ; cette valeur DIAG est comparée à la même valeur Valref de référence que précédemment, par l'additionneur 25 qui délivre en 5 sortie un indicateur Idiag. En pratique, lorsque la valeur DIAG de diagnostic dépasse la valeur Valref de référence de plus d'une valeur Vseuil, l'indicateur Idiag peut prendre la valeur 1 pour indiquer que le spectre Vtest est représentatif d'un dis fonctionnement d'une partie de l'aéronef ; au contraire, tant que la 10 différence (DIAG-Valref) reste inférieure à cette valeur Vseuil, l'indicateur Idiag prend la valeur nulle pour indiquer que le spectre Vtest est représentatif d'un bon fonctionnement de la partie considérée de l'aéronef. Dans le mode de réalisation correspondant aux figures 3 et 4, quatre vecteurs VR1, VR2, VR3, VR4 de référence sont utilisés (figure 3) pour 15 déterminer les paramètres du module 23. De la même façon que décrit ci avant, ces quatre vecteurs sont lus dans la mémoire 21, 22, 210, 220 et sont appliqués aux quatre entrées du module 23 qui fournit en sortie une valeur DIAG. La différence CorParam entre la valeur DIAG et la valeur Valref lue 20 dans la mémoire 24, est utilisée pour modifier, de façon itérative, les paramètres du module 23 jusqu'à obtenir en sortie une valeur DIAG qui soit égale à la valeur Valref - qui est par exemple égale à 1 - . Les valeurs des paramètres ainsi obtenues sont enregistrées et le module 23 peut alors être utilisé pour diagnostiquer globalement l'état 25 de composants de l'aéronef représentés par deux vecteurs Vtestl et Vtest2 à analyser. 2905763 I1 Ces deux vecteurs sont lus dans la mémoire 26, 260 et sont appliqués (figure 4), ainsi que les deux vecteurs VR1 et VR2 de référence aux entrées du module 23. Le module 23 fournit en sortie une valeur DIAG qui est comparée, 5 par le comparateur/additionneur 25, à la valeur Valref ; lorsque la valeur Idiag fournie en sortie est élevée, on peut générer un signal/message d'alarme, puis analyser séparément les vecteurs Vtestl et Vtest2, par comparaison avec un seul vecteur de référence (cf. figures 1 et 2) ou par comparaison avec trois des quatre vecteurs VR1 à VR4 (cf. figures 3 et 4). 10 Lorsqu'au contraire la valeur Idiag est inférieure à la valeur seuil, on peut générer un signal/message indiquant que les deux spectres analysés reflètent un comportement satisfaisant de l'aéronef considéré. Par référence aux figures 5, 7 et 8 notamment, on observe que le module 23 comporte des opérateurs regroupés en cinq couches CO1 à COS 15 et interconnectés pour former un réseau. Les modules représentés sur ces figures comportent deux entrées 27, 28 et une sortie 31. La première couche CO1 est constituée d'opérateurs A, à A5 et B, à B5 d'appartenance qui sont chacun caractérisés par une fonction 20 d'appartenance. Chaque opérateur d'appartenance - tel que A, - est relié à une entrée - telle que 27 du module 23 par laquelle un vecteur - tel que VR1 - est appliqué en entrée ; l'opérateur fournit en sortie un vecteur - dit vecteur de degré d'appartenance - de même dimension que le vecteur appliqué en 25 entrée, dont chaque composante résulte de la transformation de la composante de même indice du vecteur en entrée, par la fonction d'appartenance de l'opérateur considéré. 2905763 12 Les opérateurs A, et A2 du module 23 fournissent ainsi, lorsque le vecteur VR1 est appliqué à l'entrée 27, deux vecteurs de degré d'appartenance (cf. figure 5), et les opérateurs A1-A5 fournissent cinq vecteurs de degré d'appartenance (cf. figures 7 et 8). 5 Le fonctionnement des opérateurs Bj (B1 - B5) reliés à l'entrée 28 est identique à celui des opérateurs Ai (A, - A5). La seconde couche CO2 est constituée d'opérateurs Ili,j de multiplication qui reçoivent chacun deux vecteurs de degré d'appartenance délivrés par les opérateurs de la couche CO1, et fournissent en sortie une 10 valeur Wi,j de produit de ces couples de vecteurs ; plus précisément, cette valeur est égale au produit d'un des deux vecteurs par la transposée du second vecteur. La troisième couche CO3 est constitué d'opérateurs N de normalisation qui reçoivent chacun un produit d'un couple de vecteurs 15 d'appartenance et fournissent en sortie une valeur normalisée de produit. Le nombre de ces opérateurs est identique au nombre des opérateurs de multiplication ; ce dernier nombre est égal au produit des nombres d'opérateur d'appartenance respectivement associés aux différentes entrées, soit 25 opérateurs FI et 25 opérateurs N dans le cas illustré figures 7 et 8. 20 La quatrième couche CO, est constituée de 25 opérateurs Pi,j (P1 , à P51,... ,P, 5 à PS 5) de pondération respective des vecteurs d'entrée (tels que VR1 et VR2) par les produits normalisés déterminés par les opérateurs de la couche CO3, et la cinquième couche CO5 est constituée d'un seul opérateur de sommation des résultats délivrés par les opérateurs Pi,j de pondération. 25 Les sorties des couches CO, et CO5 étant des vecteurs de même dimension que celle des vecteurs d'entrée VR1 et VR2, la valeur Diag peut être définie comme la norme du vecteur délivré par le sommateur de la couche CO5. 2905763 13 La structure et le fonctionnement du module 23 sont similaires à ceux des réseaux adaptatifs multicouches - ou réseaux neuroflous - décrits dans le document ANFIS : Adaptative - Network - Based Fuzzy Inference System , Jyh-Shing Roger Jang, IEEE Transactions on Systems, Man and 5 Cybernetics, Vol. 23, 1993, pp 665-684. En particulier, le module 23 correspondant aux figures 7 et 8 s'apparente à un système d'inférence floue à deux entrées et une sortie, et dont la base de règle comporte cinq règles de type Takagi-Sugéno de forme : _-) -h - 10 Si VR1 est Ai et VR2 est Bj alors DIAGi, j = ai, j VR1 + /Ji,j VR2 + Si ,j. Ces règles sont combinées (pondérées et sommées) pour obtenir le vecteur DIAG dont la norme (moyenne quadratique des composantes par exemple) est la valeur Diag de diagnostic qui est comparée à la valeur Valref de référence. 15 Plus précisément, une fonction paramétrique d'appartenance ,um (x) est respectivement associée à chaque noeud (ou opérateur) de la première couche CO,, la fonction d'appartenance pouvant être définie par l'équation : 1+(x-cm) 26m am 20 dont la représentation graphique CCI-CC5 présente, comme illustré figure 6, une forme de cloche. Comme mentionné ci avant, cette fonction est appliquée à chaque composante d'un vecteur présenté en entrée, et l'application de cette fonction à une composante produit la composante de même indice du 25 vecteur correspondant en sortie, dit vecteur de degré d'appartenance. 2905763 14 Chacune de ces fonctions d'appartenance se caractérise par la valeur particulière des trois paramètres am, bm et c,,,, ces trois valeurs étant déterminées par itérations dans la phase préparatoire. En pratique, on peut par exemple déterminer (par itérations) une 5 valeur aä commune aux cinq paramètres a,-a5, et déterminer (par itérations) une valeur bä commune aux cinq paramètres b,-b5, pour ce qui concerne les opérateurs d'appartenance Al-A5 reliés à l'entrée 27 ; on peut faire de même pour les opérateurs B1-B5. Comme illustré figure 6, les valeurs des paramètres et-c5 10 correspondent à l'abscisse, sur ce graphe, des sommets respectifs des courbes en cloches CC1-CC5 caractéristiques des fonctions ,11(x) - p5(x) d'appartenance respective des opérateurs Al-A5. Dans la phase préparatoire, la détermination itérative des valeurs des paramètres am, bm et cm des opérateurs d'appartenance, peut se faire par 15 une méthode du gradient, par rétropropagation d'erreurs ; la détermination itérative des valeurs des trois coefficients ai,j ,Bi,j et U1 ,j de pondération caractérisant respectivement chacun (les opérateurs de pondération Pi,j peut se faire par une méthode directe, en particulier par la méthode des moindres carrés. 20 Contrairement aux opérateurs Ai, Bj, et Pi,j, les opérateurs II (couche CO2) et N (couche CO3) réseau d'opérateurs du module 23, ne sont pas paramétrés. Chaque opérateur FI délivre en sortie le produit des deux vecteurs qui lui sont présentés en entrée, i.e. le produit d'un vecteur de degré 25 d'appartenance issu d'un premier vecteur d'entrée, par la transposée d'un vecteur de degré d'appartenance issu d'un second vecteur d'entrée, soit : Wi, j = VR1 i . VR2 j '' 2905763 15 Ce simple opérateur multiplicateur peut être remplacé par un opérateur de type t-norme ou t-conorme, i.e. un opérateur borné monotone, commutatif, et associatif. Chaque opérateur N de normalisation de la couche CO, délivre en 5 sortie une valeur normalisée Wi, j définie par Wi, j Wi,j - Wk, t k,l Ces valeurs normalisées sont utilisées par les opérateurs Pi,j de pondération des entrées qui délivrent une sortie Pi, j définie par la relation : 10 Pi, i = Wi, j (ai,j vl + f3i,j v2 + bi, j ) vl et v2 étant les deux vecteurs (tels que VR1, Vtestl, VR2, Vtest2) présentés en entrée du module 23. Les vecteurs Pi, j délivrés en sortie de la couche CO, sont additionnés par le sommateur de la couche CO5 qui délivre en sortie un 15 vecteur DIAG, ces vecteurs présentant la même dimension que les vecteurs d'entrée. La norme DIAG du vecteur DIAG constitue alors la valeur DIAG de diagnostic qui est comparée à la valeur de référence. Le procédé selon l'invention a été appliqué à des données de 20 mesures enregistrées après des vols successifs d'un hélicoptère et a permis d'obtenir des valeurs DIAG respectives de diagnostic à l'issue de ces vols ; l'évolution de l'amplitude de cette valeur de diagnostic, qui est illustrée à titre d'exemple figure 11, a montré que la valeur de référence Valref, 2905763 16 permet de détecter le premier vol \Tm après lequel est apparu un disfonctionnement d'un mécanisme de la boîte de transmission principale de l'hélicoptère. REVENDICATIONS 1. Procédé de diagnostic d'un aéronef par analyse d'au moins un spectre (Vtest, Vtestl, Vtest2) obtenu à partir de mesures pratiquées sur l'aéronef, en particulier à partir de mesures vibratoires, dans lequel on applique à l'entrée d'un module (23) de calcul le(s) spectre(s) à analyser ainsi qu'un premier spectre (VR1, VR2, VR3, VR4) de référence au moins, le module de calcul déterminant une valeur (Diag) de diagnostic en fonction du (des) spectre(s) à analyser., du (des) spectre(s) de référence et de paramètres (am, bRä cm, ai,j, (3i,j, Si,j) déterminés à partir d'au moins deux spectres de référence incluant ledit premier spectre de référence. 2. Procédé selon la revendication 1 dans lequel : - dans une phase préparatoire, on applique au moins deux spectres de référence à l'entrée du module (23) de calcul paramétré qui produit en sortie une valeur de diagnostic, puis on modifie par itérations les paramètres du module de calcul pour obtenir en sortie une valeur (Valref) de référence ; puis -dans une phase de diagnostic, on applique un ou plusieurs spectres à analyser ainsi qu'une partie (certain(s)) des spectres de référence à l'entrée du module de calcul paramétré qui produit en sortie la valeur de diagnostic ; et on diagnostique un défaut lorsque cette dernière valeur dépasse la valeur de référence. 3. Procédé selon la revendication 1 ou 2, dans lequel on utilise un module de calcul comportant une pluralité d'opérateurs structurés en couches et en réseau, incluant des opérateurs (Ai, Bj) paramétriques d'appartenance et des opérateurs (Pi,j) paramétriques de pondération des spectres présentés en entrée. 4. Procédé selon l'une quelconque des revendications 1 à 3 dans lequel, dans la phase de diagnostic, on applique au moins un vecteur de 2905763 18 données spectrales - ou données fréquentielles - à diagnostiquer et au moins un vecteur de données spectrales de référence à l'entrée du module de calcul qui détermine, pour chaque vecteur en entrée, plusieurs vecteurs de degré d'appartenance respective des composantes du vecteur à plusieurs 5 classes, et qui calcule une valeur (unique) de diagnostic en fonction des vecteurs de degré d'appartenance, et dans lequel on signale une défaillance de l'aéronef lorsque la valeur de diagnostic dépasse un seuil û ou valeur de référence-. 5. Procédé selon l'une quelconque des revendications 1 à 4 dans 10 lequel chacun des vecteurs de données spectrales résulte d'une transformée de Fourier appliquée à des données de mesures résultant de l'échantillonnage de signaux délivrés par des microphones ou accéléromètres à une fréquence au moins égale à 1 kHz environ (notamment une fréquence de l'ordre de 5, 10, 50 ou 250 kHz), chacun des 15 vecteurs de données comportant de l'ordre de un ou plusieurs milliers de composantes. 6. Procédé selon l'une quelconque des revendications 1 à 5 dans lequel on transforme des séries de données temporelles de mesure en une série/vecteur de données fréquentielles de mesure, la transformation 20 comportant les opérations successives suivantes : - synchronisation de plusieurs séries de données temporelles de mesure, - calcul d'une série moyenne de données temporelles de mesure à partir de plusieurs séries de données temporelles de mesure 25 synchronisées, - transformation du domaine temporel au domaine fréquentiel pour obtenir le vecteur de données fréquentielles de mesure à partir de la série moyenne de données temporelles. 2905763 19 7. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel le module (23) de calcul comporte une structure en couches incluant une couche (COI) d'opérateurs d'appartenance, et dans lequel on applique un vecteur en entrée des opérateurs d'appartenance du module (23) pour 5 obtenir en sortie (le chaque opérateur d'appartenance un vecteur de degrés d'appartenance qui présente une dimension identique à celle du vecteur en entrée. 8. Procédé selon la revendications 7, dans lequel le module (23) de calcul inclue une couche (CO2) d'opérateurs (17) de multiplication, et dans 10 lequel on applique chaque couple de vecteurs de degrés d'appartenance en entrée des opérateurs de multiplication du module (23) pour obtenir en sortie de chaque opérateur de multiplication une valeur de produit du couple de vecteurs correspondant, ce produit représentant le degré d'activation d'une règle associée au couple de classes d'appartenance ??? 15 correspondant au couple de vecteurs. 9. Procédé selon la revendication 8, dans lequel le module (23) de calcul inclue une couche (CO3) d'opérateurs (N) de normalisation, et dans lequel on applique chaque valeur de produit de vecteurs en entrée des opérateurs de normalisation du module (23) pour obtenir en sortie de 20 chaque opérateur de normalisation une valeur normalisée de produit du couple de vecteurs correspondant. 10. Procédé selon la revendication 9, dans lequel le module (23) de calcul inclue une couche (CO4) d'opérateurs (Pi,j) de pondération des vecteurs d'entrée, et dans lequel on applique les valeurs normalisées de 25 produit de vecteurs et les vecteur d'entrée en entrée des opérateurs de pondération pour obtenir en sortie de chaque opérateur de pondération un vecteur pondéré. 11. Procédé selon la revendication 10, dans lequel le module (23) de calcul inclue une couche (CO5) comportant un opérateur (~) de sommation des vecteurs pondérés, et dans lequel on applique chaque vecteur pondéré en entrée de l'opérateur de sommation pour obtenir en sortie de 2905763 20 l'opérateur de sommation un vecteur global, et on calcule la valeur de diagnostic en fonction de (égale à) une norme du vecteur global. 12. Procédé selon l'une quelconque des revendications 7 à 11, dans lequel les opérateurs d'appartenance de la lère couche (CO1) présentent 5 une courbe (CC1-CC5) d'appartenance continue et dérivable en forme générale de cloche ou de gaussienne. 13. Programme comportant un code fixé sur un support ù tel qu'une mémoire ù ou matérialisé par un signal, le code étant lisible et/ou exécutable par au moins une unité de traitement de données - telle qu'un 10 processeur, pour diagnostiquer le fonctionnement (bon ou mauvais) d'un aéronef, le code comportant des segments de code pour effectuer les opérations d'un procédé selon l'une quelconque des revendications 1 à 12. 14. Système de diagnostic d'un aéronef, qui comporte : - un organe de lecture d'un support de données qui est agencé pour 15 lire des donnés de mesures effectuées sur l'aéronef, - une base de données contenant des données de spectres de référence pour l'aéronef, - un dispositif de transformation des données de mesure du domaine temporel au domaine fréquentiel, qui est connecté à l'organe de 20 lecture pour en recevoir les données de mesures et pour délivrer en sortie des données spectrales à analyser, et - un organe de calcul qui est relié à la base de données et au dispositif de transformation, et qui est programmé pour effectuer les opérations d'un procédé selon l'une quelconque des revendications 1 25 à 12.