Brevet FR3087265A1 — publié le 20200417

Arrière-plan de l'invention L'invention concerne le domaine général des machines tournantes. De telles machines sont utilisées dans de nombreux contextes, et notamment dans l'aéronautique. 11 peut s'agir par exemple d'une turbine, d'un compresseur, d'un engrenage ou encore d'un roulement, etc. La surveillance de l'état de telles machines tournantes est cruciale et présente de nombreux intérêts, notamment en termes de prévisibilité des opérations de réparation et/ou de remplacement, de sécurité (pour permettre l'arrêt d'équipements défectueux avant qu'un danger ne se présente), d'optimisation des opérations de maintenance (pour permettre notamment une maintenance prédictive), et de fiabilité. Pour réaliser la surveillance d'une machine tournante, une technique largement appliquée et particulièrement efficace est l'analyse 15 des signaux vibratoires et/ou acoustiques générés par la machine tournante. Ces signaux transportent en effet avantageusement de nombreuses informations sur les défauts susceptibles d'affecter les différentes pièces mécaniques en mouvement se trouvant au sein de la machine tournante, qui peuvent être révélés au moyen d'une analyse 20 spectrale (au premier ordre ou à un ordre supérieur). Plus précisément, le signal vibratoire ou acoustique généré par la machine tournante comprend une multitude de contributions engendrées par différentes sources mécaniques présentes au sein de la machine tournante, celles-ci pouvant résulter notamment d'un 25 fonctionnement régulier des différentes pièces mécaniques en mouvement au sein de la machine tournante (ex. arbres, engrenages, roulements, boîte de vitesse, etc.) ou de défauts affectant ces pièces mécaniques (ex. fissure, usure, etc.). A titre d'exemple, un signal vibratoire ou acoustique mesuré sur une chaîne de transmission mécanique de puissance 30 comprenant un défaut de roulement et des engrenages, comprend notamment une contribution provenant des chocs générés par le défaut de roulement, une contribution provenant du déplacement généré au point du 35 contact des engrenages, et 2 une contribution aléatoire liée à des bruits divers (erreurs de mesures, bruit de fond, etc.). La présence de ces multiples contributions au sein du signal vibratoire ou acoustique réduit le rapport signal-sur-bruit et peut rendre 5 difficile le traitement et l'analyse du signal vibratoire. Pour pallier cet inconvénient, il existe dans l'état de la technique des méthodes de séparation de sources qui permettent de distinguer les contributions vibratoires ou acoustiques dues aux différentes pièces mécaniques (ou « sources ») présentes dans la machine tournante. 10 Une telle méthode est présentée dans l'article de F. Bonnardot et al. intitulé « Extraction of second-order cyclostationnary sources - application to vibration analysis », Mechanical Systems and Signal Processing n°19, 2005, pages 1230-121'1. Il s'agit d'une méthode de séparation de sources aveugle ne requérant aucune connaissance du nombre de sources constituant le signai vibratoire considéré, et exploitant la cyclostationnarité à différents ordres du signal vibratoire. Cette méthode s'appuie notamment sur le filtrage du signal vibratoire par un filtre de Wiener temporel cyclique, dont les coefficients sont calculés à partir de la fonction d'autocorrélation cyclique du signal. Lors de ce calcul, il est nécessaire d'évaluer l'inverse (ou la pseudo-inverse) d'une matrice carrée dont la dimension (c'est-à-dire le nombre de lignes et de colonnes) est égale à la longueur du signal vibratoire (autrement dit, pour un signal vibratoire numérique, au nombre d'échantillons dans le signal vibratoire acquis). Or cette longueur peut être de l'ordre de quelques millions d'échantillons, de sorte que l'inversion requise est extrêmement complexe et peut de surcroît, introduire des erreurs numériques. Pour remédier à ce problème, une réduction du nombre de points du signal vibratoire pris en compte pour calculer le filtre de Wiener temporel cyclique peut être envisagée. Toutefois, une telle solution entraîne des erreurs d'estimation significatives de la fonction d'autocorrélation cyclique du signal et par conséquent, des coefficients du filtre de Wiener temporel cyclique. Objet et résumé de l'invention La présente invention permet notamment de résoudre les inconvénients précités en proposant un procédé de surveillance d'une pièce mécanique en rotation associée à au moins une fréquence cyclique 5 de défaut, ce procédé de surveillance comprenant : une étape d'obtention d'un signal vibratoire ou acoustique temporel x(t) acquis au moyen d'un capteur et incluant une contribution vibratoire ou acoustique de la pièce mécanique en rotation une étape d'estimation de ladite contribution de la pièce mécanique en 10 filtrant le signal x(t) au moyen d'un filtre de Wiener ; le procédé de surveillance étant remarquable en ce qu'il comprend une étape de détermination, dans un domaine fréquentiel, d'un nombre entier K de coefficients du filtre de Wiener, chaque coefficient correspondant à une dite fréquence cyclique de défaut associée à la pièce mécanique, cette 15 étape de détermination comprenant : - une étape d'évaluation, à au moins une fréquence f, d'une matrice 1,(f) dont chaque élément4(f), i et I désignant deux entiers compris entre 1 et K, est égal au spectre de puissance croisé de signaux xi(t) = x(t)exp(-2Pra,t) et x,(t) =x (t) exp(-2jira, t) , a, et al désignant des 20 fréquences cycliques de défaut associées à la pièce mécanique en rotation une étape d'évaluation, à ladite au moins une fréquence f, d'un vecteur b(f) dont chaque élément b',(f), m désignant un entier compris entre 1 et K, est égal au spectre de puissance croisé d'un signal xm(t) = 25 x(t) exp(-2p-cr,,t) et du signal x(t), ex', désignant une fréquence cyclique de défaut associée à la pièce mécanique ; les K coefficients du filtre de Wiener étant déterminés à partir du produit d'une matrice inverse de la matrice L(f) par le vecteur b(f) évalué pour ladite au moins une fréquence f. 30 Corrélativement, l'invention vise aussi un dispositif de surveillance d'une pièce mécanique en rotation, ce dispositif de surveillance comprenant - un module d'obtention, configure pour obtenir un signal vibratoire ou acoustique temporel x(t) acquis au moyen d'un capteur et incluant une contribution vibratoire ou acoustique de la pièce mécanique en rotation 4 un module d'estimation, configuré pour estimer ladite contribution de la pièce mécanique en filtrant le signal x(t) au moyen d'un filtre de Wiener ; le dispositif de surveillance étant remarquable en ce qu'il comprend en 5 outre un module de détermination configuré pour déterminer dans un domaine fréquentiel, un nombre entier K de coefficients du filtre de Wiener, chaque coefficient correspondant à une dite fréquence cyclique de défaut associée à la pièce mécanique, ce module de détermination comprenant 10 un premier module d'évaluation, configuré pour évaluer à au moins une fréquence f, une matrice L(f) dont chaque élément La(f), i et I désignant deux entiers compris entre 1 et K, est égal au spectre de puissance croisé de signaux x, (r) = x(r)exp(-21Tra,t) et x,(t) -- x(t)exp(-2jn-att) , ai et al désignant des fréquences cycliques de défaut 15 associées à la pièce mécanique en rotation - un second module d'évaluation, configuré pour évaluer à ladite au moins une fréquence f, un vecteur h(f) dont chaque élément bin(f), m désignant un entier compris entre 1 et K, est égal au spectre de puissance croisé d'un signal xm(t) = x(t)exp(-2jrca',t) et du signal 20 x(t), a7,1 désignant une fréquence cyclique de défaut associée à la pièce mécanique ; les K coefficients du filtre de Wiener étant déterminés par le module de détermination à partir du produit d'une matrice inverse de la matrice L(f) par le vecteur (f) évalué pour ladite au moins une fréquence f. 25 L'invention propose donc une méthode de séparation de sources aveugle (c'est-à-dire qui ne requiert pas l'utilisation d'un signal de référence pour procéder à la séparation des sources) rapide et efficace, s'appuyant sur un signal vibratoire ou acoustique acquis au moyen d'un unique capteur de signaux vibratoires ou acoustiques (comprenant par 30 exemple un accéléromètre ou un microphone), et permettant d'extraire de façon optimale la contribution d'une pièce mécanique en rotation d'intérêt, embarquée dans une machine tournante pouvant intégrer une pluralité de telles pièces mécaniques. Cette méthode de séparation de sources s'appuie sur le filtrage du signal vibratoire ou acoustique acquis via le 35 capteur au moyen d'un filtre de Wiener, linéaire et périodiquement variable dans le temps. L'invention permet ainsi de débruiter le signal 5 acquis par le capteur, et notamment la contribution liée à la pièce mécanique d'intérêt, en vue de son analyse afin de détecter si celle-ci est affectée ou non d'un défaut. Aucune limitation n'est attachée à la technique utilisée à cet effet pour traiter et analyser la contribution 5 vibratoire ou acoustique extraite grâce à l'invention. Il peut s'agir par exemple d'une technique d'analyse spectrale, d'une technique d'analyse statistique (comprenant le calcul d'indicateurs statistiques tels qu'un kurtosis, un facteur de crête, une énergie, une moyenne quadratique ou RMS (pour Root Mean Square), un coefficient d'aplatissement, etc.), d'une 10 technique d'analyse temps-fréquence, d'une technique d'analyse d'ondelettes, d'une méthode paramétrique, etc. Le procédé objet de l'invention exploite la périodicité des statistiques du signal vibratoire ou acoustique acquis au moyen du capteur, et propose de déterminer les coefficients du filtre de Wiener dans 15 le domaine fréquentiel à partir de spectres de puissance croisés de signaux dérivés du signal vibratoire. De façon connue en soi, le spectre de puissance croisé de deux signaux temporels y(t) et z(t) cyclostationnaires de fréquences cycliques respectives ay et a., (exprimées en hertz) est défini par : 1 yz(f) = v-v- E[Z(f + a7)Yw( f + ay)] 20 où f désigne la fréquence, W la fenêtre temporelle d'observation des signaux y(t) et z(t) et Yw et Zw les transformées de Fourier des signaux y(t) et z(t) sur la fenêtre W. La détermination des coefficients du filtre de Wiener est donc particulièrement simple à mettre en oeuvre et avantageusement moins 25 coûteuse que celle proposée par F. Bonnardot et al. En effet, si par exemple, un estimateur de Welch est utilisé pour déterminer les coefficients du filtre de Wiener à différentes fréquences sur une fenêtre glissante comprenant N points (i.e. N fréquences, N étant par exemple pris égal à 512 ou 1024 pour une fréquence d'échantillonnage du signal de 30 25 ou 50 kHz), la détermination des coefficients du filtre de Wiener conformément à l'invention requiert N inversions de matrices de dimension KxK où K désigne le nombre de coefficients du filtre de Wiener. Au contraire la solution proposée par F. Bonnardot et al. requiert l'inversion d'une matrice de dimension LgxLg où Lg désigne la longueur du 6 signal qui, comme évoqué précédemment, peut atteindre quelques millions d'échantillons. Le choix du nombre K de coefficients du filtre de Wiener résulte d'un compromis entre efficacité du filtrage réalisé pour extraire la source 5 d'intérêt et complexité de calcul. Dans un mode préféré de réalisation, ce nombre de coefficients K du filtre est inférieur ou égal à 4. Les inventeurs ont en effet constaté que cette valeur permet d'extraire de façon aveugle et optimaie ou quasi-optimale un grand nombre de sources quelle que soit la nature de la pièce mécanique considérée (engrenage, roulement, etc.). 10 En d'autres mots, il n'est pas nécessaire de considérer un filtre avec un grand nombre de coefficients pour extraire du signal vibratoire ou acoustique acquis via le capteur la contribution de la pièce mécanique considérée. Il en résulte l'inversion de quelques centaines (ex. 512 ou 1024 dans l'exemple envisagé précédemment) de matrices de faible 15 dimension (4x4 au maximum si K=4) comparativement à la technique proposée par F. Bonnardot et al. qui suppose d'inverser une matrice ayant comme dimension quelques millions d'échantillons. La méthode proposée par l'invention ne nécessite en outre aucune connaissance a priori du nombre de sources présentes dans le 20 signal vibratoire ou acoustique acquis, et ne s'appuie sur aucune hypothèse d'indépendance des sources présentes dans le signal vibratoire ou acoustique (non valables en mécanique où les signaux émis par les différentes sources subissent des interactions convolutives et multiplicatives). Elle requiert uniquement la connaissance de la 25 cinématique de la source dont on cherche à extraire la contribution, autrement dit ici de la pièce mécanique en rotation que l'on cherche à surveiller éventuellement entachée d'un défaut que l'on souhaite détecter. Cette cinématique comprend la ou les fréquences caractéristiques de défaut de la pièce mécanique, aussi appelées 30 fréquences cycliques de défaut en raison de la rotation de la pièce mécanique et de la manifestation cyclique et périodique du défaut affectant cette pièce le cas échéant. Cette fréquence cyclique de défaut est généralement liée de manière linéaire à la rotation de la pièce mécanique que le défaut est susceptible d'affecter. A titre d'exemple, 35 lorsque la pièce mécanique considérée est une roue dentée et que l'on s'intéresse à la détection d'un défaut affectant la denture de cette roue, 7 les fréquences caractéristiques de défaut de la roue comprennent la fréquence de rotation de la roue et ses harmoniques. Ces fréquences sont des fréquences générées par la source que l'on cherche à surveiller et qui sont propres à cette source : elles ne sont pas partagées avec d'autres 5 sources à l'origine de contributions dans le signal vibratoire. Selon un autre exemple, lorsque la pièce mécanique est un roulement affecté d'un défaut, il s'agit de la fréquence d'apparition du défaut lors de la rotation du roulement et de ses harmoniques Les fréquences cycliques de défaut associées à une pièce mécanique en rotation peuvent être, dans la 10 majorité des cas, déterminées de façon théorique, ou de manière alternative, expérimentalement. Dans un mode particulier de réalisation, les éléments de la matrice L(f) vérifient : Lu(f) L((f) Autrement dit, dans ce mode de réalisation, on considère que la 15 matrice L(f) est hermitienne ce qui permet non seulement de réduire la complexité de calcul de cette matrice, mais également d'assurer une meilleure stabilité numérique des calculs (et notamment des calculs requis lors de l'inversion de la matrice L(f)). Dans un mode particulier de réalisation, le filtrage du signal x(t) 20 est réalisé, dans le domaine fréquentiel. Ceci permet d'utiliser directement les coefficients du filtre de Wiener déterminés dans le domaine fréquentiel lors de l'étape de détermination, et de simplifier les opérations mises en oeuvre lors du filtrage. 25 En variante, le filtrage peut être réalisé dans le domaine temporel en dérivant les coefficients de filtre de Wiener à utiliser dans le domaine temporel à partir des coefficients déterminés dans le domaine fréquentiel. Dans un mode particulier de réalisation, le procédé de 30 surveillance comprend en outre une étape de ré-échantillonnage angulaire du signal x(t), la contribution de la pièce mécanique étant estimée lors de l'étape d'estimation en filtrant le signal x(t) ré-échantillonné, et les K coefficients du filtre de Wiener étant déterminés à partir du signal x(t) ré-échantillonné au lieu du signal x(t). 8 Lorsque la pièce mécanique en rotation se trouve intégrée au sein d'une machine tournante, le signal x(t) acquis par le capteur regroupe les contributions des différents éléments mécaniques compris dans la machine tournante et sa périodicité est intrinsèquement liée à 5 l'angle de la rotation de la machine tournante. Lorsque la machine tournante fonctionne en régime stationnaire, cette périodicité se maintient dans le temps. En pratique toutefois, les machines tournantes subissent des fluctuations de leur régime de rotation de sorte que l'analyse temporelle du signal x(t) et des différentes contributions le constituant 10 peut induire des erreurs si ces fluctuations ne sont pas prises en compte. Pour pallier cet inconvénient, un ré-échantillonnage du signal x(t) acquis par le capteur est mis en oeuvre en tenant compte de la position de l'angle de rotation de la machine tournante ou d'un arbre de référence équipant celle-ci et l'entraînant en rotation. Cette position peut être aisément 15 obtenue en utilisant par exemple un tachymètre à proximité de l'arbre de référence. En variante, le capteur, et plus particulièrement sa chaîne d'acquisition, peuvent être configurés pour acquérir des échantillons correspondant directement à un pas angulaire fixe de la machine 20 tournante. Cette solution est toutefois plus complexe à mettre en oeuvre qu'un ré échantillonnage du signal temporel x(t) acquis par le capteur avec une période d'échantillonnage donnée. Dans un autre mode de réalisation, le procédé de surveillance comprend en outre une étape de décomposition du signal x(t) ou du 25 signal x(t) ré échantillonné en une composante déterministe et en une composante aléatoire, la contribution de la pièce mécanique étant estimée lors de l'étape d'estimation en filtrant la composante parmi la composante déterministe et la composante aléatoire qui contient la contribution vibratoire ou acoustique de la pièce mécanique, et les K coefficients du 30 filtre de Wiener étant déterminés à partir de cette composante au lieu du signal x(t). Cette décomposition en composantes déterministe et aléatoire, et la sélection de la composante contenant la contribution vibratoire ou acoustique recherchée, permettent d'améliorer le rapport signal-sur-bruit 35 des traitements effectués. On réduit en effet par ce biais les signaux interférents par rapport à la contribution recherchée. La nature de la 9 composante à préserver dépend de la pièce mécanique à laquelle on s'intéresse. Ainsi par exemple, si on s'intéresse à la contribution liée à un défaut de roulement qui est de nature aléatoire, seule la composante aléatoire sera conservée et utilisée pour déterminer les coefficients du 5 filtre de Wiener et extraire cette contribution. Si en revanche, on vise l'extraction d'une contribution liée à un défaut affectant un engrenage qui est de nature déterministe, c'est la composante déterministe qui sera conservée et utilisée par le procédé de surveillance selon l'invention. Dans un mode particulier de réalisation, les différentes étapes 10 du procédé de surveillance sont déterminées par des instructions de programmes d'ordinateurs. En conséquence, l'invention vise aussi un programme d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en oeuvre dans un dispositif de surveillance ou plus 15 généralement dans un ordinateur, ce programme comportant des instructions adaptées à la mise en oeuvre des étapes d'un procédé de surveillance tel que décrit ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de 20 code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable. L'invention vise aussi un support d'informations lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur 25 tel que mentionné ci-dessus. Le support d'informations peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen 30 d'enregistrement magnétique, par exemple un disque dur. D'autre part, le support d'informations peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier 35 téléchargé sur un réseau de type Internet. 10 Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en question. 5 L'invention vise également un système comprenant : - un capteur, configure pour acquérir un signal vibratoire ou acoustique généré par au moins une pièce mécanique en rotation ; un dispositif de surveillance selon l'invention, configure pour déterminer une contribution vibratoire ou acoustique de la pièce 10 mécanique en rotation à partir du signal vibratoire ou acoustique acquis par le capteur. Le système selon l'invention bénéficie des mêmes avantages cités précédemment que le dispositif et le procédé de surveillance selon l'invention. 15 Le système peut être un équipement portable permettant de surveiller à la volée l'état d'une pièce mécanique en rotation, ou être intégré dans la structure comprenant ladite pièce mécanique ou sur un banc d'essai. Dans un mode particulier de réalisation, le système 20 comprend en outre un module de restitution de la contribution vibratoire ou acoustique déterminée par le dispositif de surveillance. Un tel module de restitution est par exemple un écran ou un haut-parleur. Ce mode de réalisation facilite l'inspection (visuelle ou auditive) de la pièce mécanique et la caractérisation d'éventuels défauts 25 affectant ladite pièce. Par exemple, dans le cas d'un défaut local, la contribution vibratoire ou acoustique de la pièce peut comprendre des signatures transitoires et répétitives générées par des chocs aisément détectables sur un écran. Dans un mode particulier de réalisation, le capteur du système 30 selon l'invention comprend un accéléromètre (délivrant alors un signal vibratoire) ou un microphone (délivrant un signal acoustique). Ce type de capteurs est couramment utilisé dans le contexte des machines tournantes, en particulier lorsque ces machines tournantes sont des moteurs à combustion ou à explosion, de sorte que l'invention est 35 particulièrement simple à mettre en oeuvre. On note que l'invention ne se 11 limite pas à des moteurs d'aéronefs mais peut s'appliquer avantageusement à tout type de machines tournantes. Dans un autre mode de réalisation, le système peut également embarquer un tachymètre permettant de mesurer la vitesse de rotation de 5 la machine tournante ou de la pièce mécanique en rotation considérée, notamment en vue de ré-échantillonner le signal vibratoire ou acoustique acquis par le capteur comme décrit précédemment. On peut également envisager, dans d'autres modes de réalisation, que le procédé de surveillance, le dispositif de surveillance, et 10 le système selon l'invention présentent en combinaison tout ou partie des caractéristiques précitées. Brève description des dessins D'autres caractéristiques et avantages de la présente invention 15 ressortiront de la description faite ci-dessous, en référence aux dessins annexés qui en illustrent un exemple de réalisation dépourvu de tout caractère limitatif. Sur les figures : la figure 1 représente, dans son environnement, un système conforme à l'invention dans un mode particulier de réalisation ; 20 la figure 2 illustre schématiquement l'architecture matérielle d'un dispositif de surveillance conforme à l'invention, inclus dans le système de la figure 1, dans un mode particulier de réalisation ; la figure 3 représente les principales étapes d'un procédé de surveillance selon l'invention tel qu'il est mis en oeuvre par le dispositif 25 de surveillance de la figure 2 ; la figure 4 illustre l'architecture d'un filtre de Wiener et la figure 5 représente les principales étapes mises en oeuvre par le dispositif de surveillance pour déterminer les coefficients du filtre de Wiener dans le domaine fréquentiel conformément à ['invention. 30 Description détaillée de l'invention La figure 1 représente, dans son environnement et dans un mode particulier de réalisation, un système 1 conforme à l'invention permettant la surveillance d'une machine tournante 2 comprenant au 35 moins une pièce mécanique 3 en rotation à laquelle on s'intéresse plus 12 particulièrement pour déterminer si celle-ci est affectée ou non d'un défaut. Dans l'exemple illustré à la figure 1, la machine tournante 2 est une boîte d'engrenages d'un moteur à combustion équipant un aéronef tel 5 qu'un avion, et la pièce mécanique 3 en rotation est une roue dentée d'un engrenage particulier de la boîte d'engrenages. Cette roue dentée lorsqu'elle se trouve en rotation est caractérisée par une cinématique, et plus particulièrement, est associée une ou plusieurs fréquences cycliques de défaut qui caractérisent, lorsque la roue dentée présente un défaut, la 10 répétitivite de la manifestation de ce défaut (qui est périodique dans des conditions stationnaires). Les fréquences cycliques de défaut sont donc propres à chaque pièce mécanique considérée. Dans le cas d'une roue dentée d'un engrenage telle que considérée à titre illustratif ici, ces fréquences cycliques de défaut comprennent typiquement la fréquence de 15 rotation de la roue dentée ainsi que ses harmoniques (i e ses multiples positifs ou négatifs). Préférentiellement on ne considère pour mettre en oeuvre l'invention qu'un nombre faible d'harmoniques, afin d'éviter les harmoniques d'ordre supérieur souvent assez bruitées. De façon connue de l'homme du métier, les fréquences 20 cycliques de défaut associées à une pièce mécanique en rotation peuvent être en général déterminées de façon théorique ou à défaut de façon empirique, et dépendent de la pièce mécanique en rotation à laquelle on s'intéresse et du défaut que l'on souhaite caractériser. Ainsi par exemple, si l'on s'intéresse à une bague externe d'un roulement, les fréquences 25 cycliques de défaut associées à cette bague externe sont la fréquence de passage de la bille sur la bague externe (qui est proportionnelle à la fréquence de rotation du roulement) et ses harmoniques. Dans la suite de la description on désigne par A l'ensemble des fréquences cycliques de défaut considérées associées à la pièce 30 mécanique 3. On suppose que l'ensemble A contient un nombre K de fréquences cycliques de défaut associées â la pièce mécanique 3, K désignant un entier supérieur ou égal â 1; dans la suite de la description, chaque fréquence cyclique de l'ensemble A est indexée par un entier k=1,...K et notée ak. Le choix de K est détaillé plus en détail 35 ultérieurement. 13 Bien entendu, les hypothèses précédentes ne sont en aucun cas limitatives et l'invention s'applique à d'autres pièces mécaniques en rotation (par exemple à un roulement, à une aube, etc.) ainsi qu'à d'autres machines tournantes embarquant de telles pièces mécaniques, et 5 dont on souhaite assurer une surveillance. Elle s'applique par ailleurs à d'autres domaines que l'aéronautique. Dans le mode de réalisation décrit ici, le système 1 est un appareil portable comprenant un capteur 4 de vibrations incluant ici un accéléromètre 4A configure pour mesurer sur une durée déterminée un 10 signal vibratoire généré par la boîte d'engrenages et les pièces mécaniques en rotation embarquées dans celle-ci lors de son fonctionnement (dont la pièce mécanique 3), ainsi qu'un dispositif de surveillance 5 conforme à l'invention, apte à traiter le signal vibratoire fourni par le capteur 4 de vibrations. Le système 1 peut être ainsi 15 aisément utilisé pour surveiller différentes types de machines tournantes, à bord d'un aéronef ou sur un banc d'essai. Il est lors de cette utilisation placé à proximité de la machine tournante que l'on souhaite surveiller (contenant la pièce mécanique en rotation 3), de sorte que le signal vibratoire mesuré par l'accéléromètre 4A comprend une contribution 20 vibratoire due à la pièce mécanique en rotation 3. Il convient de noter que le mode de réalisation décrit ici s'attache à un signal vibratoire acquis par un capteur 4 comprenant un accéléromètre, toutefois l'invention s'applique également à un signal acoustique, acquis par exemple au moyen d'un microphone ou de tout 25 autre capteur adapté à acquérir un tel signal. De façon connue en soi, le capteur 4 de vibrations comprend également une chaîne ou une carte d'acquisition 4B, apte à conditionner le signal acquis par l'accéléromètre 4A, et notamment à numériser (voire à amplifier et à filtrer) le signal (analogique) mesuré par l'accéléromètre 4A. 30 Une telle chaîne d'acquisition est connue en soi, et n'est pas décrite en détail ici. On désigne par x(t) le signai vibratoire numérique disponible à la sortie de la chaîne d'acquisition 48 le signal x(t) est un signal temporel échantillonné avec une période d'échantillonnage notée Te (ou 35 de façon équivalente avec une fréquence d'échantillonnage fe). Par exemple, fe=50kHz. La longueur du signal x(t) dépend de la durée de la 14 mesure effectuée par l'accéléromètre 4A. Préférentiellement, on choisit une durée sur laquelle la pièce mécanique en rotation (et plus généralement la machine tournante dans laquelle elle est intégrée) présente des conditions de fonctionnement stationnaires, c'est-à-dire sans 5 variation significative de couple et/ou de vitesse (par exemple on s'assurera préférentiellement que l'écart type de la vitesse de rotation ne dépasse pas 5% de sa moyenne). Cette durée doit toutefois ne pas être choisie trop courte de sorte à ne pas entacher la précision de la méthode de séparation aveugle proposée par l'invention. Les inventeurs ont 10 constaté expérimentalement qu'une durée de quelques secondes pour une fréquence d'échantillonnage fe=-50 kHz représente un bon compromis. Dans le mode de réalisation décrit ici, le dispositif de surveillance 5 a l'architecture matérielle d'un ordinateur, telle que représentée schématiquement à la figure 2. Il comprend un processeur 15 6, une mémoire vive 7, une mémoire morte 8, une mémoire flash non volatile 9, et des moyens de communication 10 lui permettant de communiquer notamment avec le capteur 4, via par exemple un bus numérique. Dans le mode de réalisation décrit ici, la mémoire morte 8 du 20 dispositif de surveillance 5 constitue un support d'enregistrement conforme à l'invention, lisible par le processeur 6 et sur lequel est enregistré un programme d'ordinateur PROG conforme à l'invention, comportant des instructions pour l'exécution d'une partie des étapes du procédé de surveillance selon l'invention. Ce programme d'ordinateur 25 PROG définit des modules fonctionnels et logiciels du dispositif de surveillance 5 qui utilisent ou s'appuient sur les éléments matériels 6-10 mentionnés précédemment. Ces modules comprennent notamment, comme illustré à la figure 1 -- un module d'obtention 5A, configuré pour obtenir le signal x(t) acquis 30 au moyen du capteur 4 et incluant la contribution vibratoire générée par la pièce mécanique 3 lorsque celle-ci est en rotation - un module d'estimation 5B, configuré pour estimer cette contribution en filtrant le signal x(t) au moyen d'un filtre de Wiener, linéaire et périodiquement variable en temps. Un tel filtre est connu en soi et décrit 35 par exemple dans le document de F. Bonnardot et al. cité précédemment. 15 Il est décrit plus en détail ici ultérieurement en référence à la figure 4 notamment ; et un module de détermination 5C configure pour déterminer dans un domaine fréquentiel, les coefficients du filtre de Wiener. Ces coefficients 5 sont au nombre de K, chaque coefficient correspondant à une fréquence cyclique de défaut 04,, k=1,...,K, de l'ensemble A associée à la pièce mécanique 3. Conformément à l'invention, le module de détermination 5C comprend o un premier module d'évaluation 50., configure pour évaluer à au 10 moins une fréquence f, une matrice carrée 1.(f) de dimension K, dont chaque élément Lii(f), est à la fréquence f un scalaire ; et o un second module d'évaluation 5C2, configure pour évaluer à ladite au moins une fréquence f, un vecteur b(f) dont chaque élément 15 b(f), m=1,...,K est à la fréquence f un scalaire. La matrice L(f) et le vecteur b(f) sont évalués par les modules d'évaluation 50. et 5C2 pour au moins une fréquence f à partir de spectres de puissance croisés de signaux dérivés du signal vibratoire x(t) comme décrit ultérieurement en référence aux différentes étapes du procédé de 20 surveillance selon l'invention. En outre, dans le mode de réalisation décrit ici, le module de détermination 5C comprend o un module d'inversion 5C3, configure pour évaluer la matrice inverse de la matrice carrée L(f) évaluée par le premier module d'évaluation 5C1 ; et 25 o un module de calcul 5C4, configure pour déterminer le produit de la matrice inverse évaluée par le module d'inversion 5C3 et du vecteur b(f) évalué par le second module d'évaluation 5C2, le produit étant évalué à ladite au moins une fréquence f. Dans le mode de réalisation décrit ici, le dispositif de 30 surveillance 5 comprend également un module de ré-échantillonnage angulaire 5D, configure pour réaliser un ré-échantillonnage du signal vibratoire acquis par le capteur 4 permettant de tenir compte de la vitesse de rotation et/ou du couple de la machine tournante 2 et de la pièce mécanique 3 (la vitesse de rotation de la pièce mécanique 3 varie 35 généralement de manière linéaire par rapport à la vitesse de rotation de la machine tournante 2), ainsi qu'un module de décomposition 5E configure 16 pour décomposer le signal vibratoire ré-échantillonné en une composante aléatoire et en une composante déterministe. Il convient de noter que ces deux modules sont optionnels. En variante, aucun ré-échantillonnage n'est appliqué au signal acquis et échantillonné par le capteur 4 (i.e. le dispositif de surveillance ne comprend pas de module 5D). Dans une autre variante encore l'échantillonnage du signal vibratoire mesuré par l'accéléromètre 4A appliqué par la chaîne d'acquisition 4B du capteur 4 tient compte de susceptibles variations des 10 conditions de fonctionnement de la machine tournante 2 et de la pièce mécanique 3. Dans une autre variante encore, aucune décomposition du signal vibratoire échantillonné ou ré-échantillonné n'est réalisée par le dispositif de surveillance (i.e. celui-ci ne comprend pas de module de 15 décomposition 5E). Dans le mode de réalisation décrit ici, le système 1 comprend également un capteur 11 permettant de mesurer l'angle de rotation de la pièce mécanique 3 ou de la machine tournante 2 tel que par exemple un tachymètre ; on note que l'angle de rotation de la pièce mécanique 3 et 20 celui de la machine tournante 2 sont liés par une relation linéaire de sorte qu'il est possible de déduire de la connaissance de l'un, la connaissance de l'autre. Le système comprend aussi ici un module 12 de restitution de la contribution vibratoire de la pièce mécanique 3 estimée par le dispositif de surveillance 5. Ce module 12 de restitution comprend ici un écran 25 permettant de visualiser la contribution vibratoire extraite du signal x(t). En variante, il peut comprendre un haut-parleur, permettant de restituer de façon auditive cette contribution. Nous allons maintenant décrire, en référence aux figures 3 et 4, les principales étapes du procédé de surveillance mis en oeuvre par le 30 système 1, et plus particulièrement par le dispositif de surveillance 5 pour extraire d'un signal vibratoire x(t) acquis par le capteur 4 lors d'un fonctionnement de la machine tournante 2, la contribution vibratoire due à la pièce mécanique 3 afin de permettre son analyse et la détection d'un éventuel défaut affectant la pièce mécanique 3. 35 Comme mentionné précédemment, on suppose que le signal vibratoire x(t) a été acquis par le capteur 4, et plus particulièrement par 17 l'accéléromètre 4A (étape E10), alors que les conditions de fonctionnement de la machine tournante 2 (et corrélativement de la pièce mécanique 3) étaient stationnaires ou quasi-stationnaires, autrement dit à couple ou vitesse de rotation constante ou quasi-constante. De façon 5 connue, et comme décrit précédemment, on suppose ici que le signal vibratoire x(t) désigne le signal mesuré par l'accéléromètre 4A puis échantillonné par la chaîne d'acquisition 413 à la fréquence d'échantillonnage fe. La longueur du signal vibratoire temporel x(t) est de Lg échantillons espacés de la période Te. 10 Le signal vibratoire x(t) ainsi échantillonné est transmis au dispositif de surveillance 5 et plus particulièrement à son module d'obtention 5A (par l'intermédiaire notamment de ses moyens de communication 10) (étape E20). Dans le mode de réalisation décrit ici, le module de ré- 15 échantillonnage angulaire 5D réalise alors un ré-échantillonnage angulaire du signal vibratoire échantillonné x(t), en utilisant une mesure de l'angle de rotation de la machine tournante 2 réalisée par le tachymètre 12 sur la durée du signal vibratoire x(t) (étape E30). Ce ré-échantillonnage est réalisé en fonction de l'angle de rotation de la machine tournante 2 pour 20 tenir compte de variations éventuelles de la vitesse de rotation ou du couple de la machine tournante 2 lors de la durée d'acquisition du signal vibratoire x(t). Un tel ré-échantillonnage est connu en soi et n'est pas décrit plus en détail ici. Il peut être réalisé par exemple de façon identique à ce qui est décrit dans l'article de K.R. Fyle et al. intitulé « Analysis of 25 computed order tracking », Mechanical Systems and Signal Processing, n°11(2), pages 187-205, 1997. Par souci de simplification, on conserve dans la suite de la description la notation x(t) pour désigner le signal vibratoire ré-échantillonné fourni en sortie du module de ré-échantillonnage 5D. 30 Puis, dans le mode de réalisation décrit ici, le signal vibratoire ré-échantillonné x(t) est traité par le module de décomposition 5E du dispositif de surveillance 5. Plus spécifiquement, celui-ci décompose le signal vibratoire ré-échantillonné x(t) en une composante déterministe et en une composante aléatoire (étape E40). Il procède à cet effet comme 35 décrit par exemple dans l'article de J. Antoni et al. intitulé « Unsupervised Noise Cancellation for Vibration Signais : Part II - A Navel Frequency- 18 Domain Algorithm Mechanical Systems and Signal Processing, n°18, pages 103-117, 2004. Le module de décomposition 5E sélectionne ensuite parmi les deux composantes déterministe et aléatoire obtenues, celle qui contient la 5 contribution vibratoire de la pièce mécanique 3 (étape E50). Cette sélection permet de réduire le bruit contenu dans le signal vibratoire qui va être utilisée pour extraire la contribution vibratoire de la pièce mécanique 3. La composante sélectionnée dépend de la nature de la pièce mécanique 3 considérée et du défaut susceptible de l'affecter. Dans 10 l'exemple envisagé ici, la pièce mécanique est une roue dentée d'un engrenage par définition, la contribution vibratoire d'une roue dentée d'engrenage éventuellement entachée d'un défaut est de nature déterministe. Le module de décomposition 5E sélectionne donc la composante déterministe qu'il vient d'extraire du signal ré-échantillonné et 15 conserve uniquement cette composante déterministe. Selon un autre exemple, si la pièce mécanique 3 est un roulement, la contribution vibratoire liée à un défaut de roulement est par nature aléatoire de sorte que le module de décomposition 5E sélectionne la composante aléatoire qu'il vient d'extraire du signal ré-échantillonné et 20 conserve uniquement cette composante aléatoire.. Dans la suite de la description, par souci de simplification, on conserve la notation x(t) pour désigner la composante du signal vibratoire ré-échantillonné sélectionnée par le module de décomposition 5E. Il convient par ailleurs de noter que d'une part, les étapes E30 et E40-E50 25 sont optionnelles, et que d'autre part, elles peuvent être réalisées dans l'ordre inverse de celui qui vient d'être décrit. Conformément à l'invention, le dispositif de surveillance 5 estime la contribution vibratoire de la pièce mécanique 3 contenue dans le signal vibratoire x(t) acquis par le capteur 4 en filtrant la composante 30 aléatoire ou déterministe sélectionnée à l'étape E50 dans le signal vibratoire ré-échantillonné (ou dans la composante sélectionnée dans le signal vibratoire x(t) échantillonné si aucun ré-échantillonnage n'est mis en oeuvre ou encore dans le signal vibratoire x(t) échantillonné ou ré-échantillonné si aucune décomposition du signal n'est mise en oeuvre). Le 35 filtrage est réalisé au moyen d'un filtre de Wiener, connu en soi et tel que décrit par exemple dans l'article de F. Bonnardot et al. cité précédemment. 19 L'invention se distingue toutefois de la technique proposée dans la façon dont les coefficients de ce filtre de Wiener sont déterminés. Afin de mieux comprendre l'invention, l'architecture d'un filtre de Wiener est représentée à la figure 4. Ce filtre est un filtre linéaire périodiquement variable dans le temps qui est optimal pour extraire d'un signal vibratoire des sources présentant des propriétés de cyclostationnarité. Il est défini dans le domaine temporel par la réponse impulsionnelle h(t,u) suivante K h(t,u) = hk(t - u)exp(-2jrraku k-1 où 12=-1, et am, k-=1,...,K désigne les fréquences cycliques associées à la 10 pièce mécanique 3 comprises dans l'ensemble A et indexées par k. Chaque coefficient hk(t) désigne lui-même la réponse impulsionnelle d'un filtre linéaire invariant dans le temps. Il convient de noter que le choix du nombre K de coefficients du filtre de Wiener résulte d'un compromis entre efficacité du filtrage réalisé 15 pour extraire la source d'intérêt et complexité de calcul. Dans un mode préféré de réalisation, ce nombre de coefficients K du filtre est inférieur ou égal à 4. Les inventeurs ont en effet constaté que cette valeur permet d'extraire de façon aveugle et optimale ou quasi-optimale un grand nombre de sources quelle que soit la nature de la pièce mécanique 20 considérée (roue dentée, engrenage, roulement, etc.). Dans l'exemple illustratif envisagé ici on considère que K=4, et que les 4 fréquences cycliques correspondant respectivement à ces 4 coefficients sont la fréquence d'engrènement, la première harmonique de cette fréquence et les deux fréquences opposées correspondant à la fréquence 25 d'engrènement et à la première harmonique. Si in(t) désigne le signal fourni en entrée du filtre de Wiener représenté à la figure 4A, sa sortie out(t) est donnée par K Out (r) = hk(t) Ø ink(t) (1) k =1 où tø désigne l'opérateur de convolution et ink(t) - in(t)exp(-2 jircekt) 20 Conformément à l'invention, pour faciliter l'évaluation des coefficients du filtre de Wiener, on se place dans le domaine fréquentiel, et le module de détermination 5C du dispositif de surveillance 5 évalue les K coefficients notés lik(f), k=1,..,K du filtre de Wiener dans ce domaine 5 fréquentiel, f désignant la fréquence (étape E60). La figure 5 représente plus en détail la façon dont les coefficients Hk(f), k=1,...,K du filtre de Wiener sont déterminés par le module de détermination 5C. Plus spécifiquement, le module de détermination 5C du dispositif 10 de surveillance 5 évalue, via son premier module d'évaluation 5C1,, une matrice carrée L(f) dont chaque élément La(e), i et I désignant deux entiers compris entre 1 et K, est égal au spectre de puissance croisé noté S xi xt(I) de signaux Mt) et xi(t) dérivés du signal vibratoire x(t) obtenu à l'issue de l'étape E50 (étape E61). Les signaux x(t) et xi(t) sont définis 15 de la façon suivante : xi(t) = x(t)exp(-2jirat) aci(t) = x(t)exp(-2frreeit) où a, et al désignent des fréquences cycliques de défaut associées à la pièce mécanique 3, contenues dans l'ensemble A et indexées respectivement par i et L 20 Comme rappelé précédemment, de façon connue en soi, le spectre de puissance croisé de deux signaux temporels y(t) et z(t) cyclostationnaires de fréquences cycliques respectives ay et a, (exprimées en hertz) est défini par : 53,7(f) = IimE[Zcif(f + crz)Vw(f + tiy)1 où f désigne la fréquence, W la fenêtre temporelle d'observation des 25 signaux y(t) et z(t) et Yw et Zw les transformées de Fourier des signaux y(t) et z(t) sur la fenêtre W. En s'appuyant sur une description de Cramer d'un processus aléatoire, on peut montrer que de façon équivalente, le spectre de puissance croisé S1(f) des signaux x(t) et x/(t) (et donc l'élément 30 1,1(f)) est défini par : S x,,x,(f) = E[(1)(,4 (f + a() X i(J + ar)]/df (2) 21 où X,(f) désigne la transformée de Fourier du signal temporel xi(t), Xi'([) désigne le transconjugué (transposé du conjugué) de X,(f), dX,(f) désigne un incrément spectral du signal fréquentiel X,(f) et E[j désigne l'opérateur de moyennage d'ensemble. On note que chaque élément Lit(f) 5 de la matrice L(f), dX,` (net dX i()) sont ici des scalaires pour une fréquence donnée f. Dans le mode de réalisation décrit ici, les signaux considérés étant réels, la matrice carrée L(f) est par définition hermitienne : il en résulte que seuls les éléments de la partie triangulaire supérieure de la 10 matrice LU) (incluant la diagonale) sont calculés à proprement parler par le premier module d'évaluation 5C1, autrement dit, seuls les éléments L,,(f) pour 1,1=1,...,k avec I i sont estimés par le premier module d'évaluation 5C1 à partir de la relation (2) donnée ci-dessus Les éléments de la partie triangulaire inférieure de la matrice L(f) sont ensuite déduits 15 des éléments calculés en utilisant la relation suivante (liée à la propriété hermitienne de la matrice LU') : Lit(f) = (f) pour avec i > I. Pour évaluer les éléments de la partie triangulaire supérieure de la matrice L(f) vérifiant la relation (1), le premier module d'évaluation 5C1 20 utilise ici un estimateur de Welch, connu de l'homme du métier, et permettant d'évaluer facilement numériquement des densités de puissance. Un tel estimateur est connu et exposé par exemple dans le document de P.D. Welch intitulé « The Use of Fast Fourier Transform for the Estimation of Power Spectra: A method based on Time Averaging over 25 Short, Modified Periodograms IEEE Transactions on Audio Electroacoustics, AU-15, 70-73, 1967. Il n'est pas décrit plus en détail ici. Il permet au premier module d'évaluation 5C1 d'évaluer la matrice L(f) pour différentes valeurs de fréquence f, correspondant à N points d'une fenêtre glissante espacés selon une résolution fréquentielle b.f 30 paramétrable. Par exemple, si la fréquence d'échantillonnage est Fe=25 kHz et que l'on désire une résolution fréquentielle de AU = 200Hz, le nombre N doit être de l'ordre de Fe/M= 125, Pour accélérer le temps de calcul requis, on peut choisir un nombre qui est une puissance de 2 proche de 125 soit N=128. 22 En variante, le premier module d'évaluation 5C1 peut utiliser d'autres techniques pour évaluer les coefficients de la partie triangulaire supérieure de la matrice L(f), comme par exemple un estimateur spectral non paramétrique tel qu'un périodogramme, un périodogramme modifié, 5 ou encore un estimateur de Barlett, ou un estimateur spectral paramétrique basé par exemple sur une méthode autorégressive ou sur les techniques MUSIC ou ESPRIT, etc. Dans une autre variante, le premier module d'évaluation 5C1 peut évaluer les coefficients de la partie triangulaire inférieure de la 10 matrice L(f) incluant la diagonale et déduire de ces coefficients ceux de la partie triangulaire supérieure en exploitant la propriété hermitienne de la matrice L(f). Puis le module d'inversion 50 évalue la matrice inverse de la matrice carrée L(f), notée L-1(f), de façon connue en soi (étape E62). Il 15 peut à cet effet, calculer la matrice pseudo-inverse de la matrice L(f) pour parer le cas où L(f) n'est pas inversible. Par ailleurs, le module de détermination 5C, par l'intermédiaire de son second module d'évaluation 5C2, évalue un vecteur b(f) dont chaque élément b,n(f), m désignant un entier compris entre 1 et K, est 20 égal au spectre de puissance croisé noté S ,m,'(f) du signal x(t) = x(t)exp(--2inamt), dérivé du signal vibratoire x(t) obtenu à l'issue de l'étape E50 et du signal x(t), am désignant la fréquence cyclique de défaut associée à la pièce mécanique 3 appartenant à l'ensemble A et indexée par m (étape E63). En d'autres mots 25 brn(f) = sx',,'(1)= Erclx*(f)dx,n(f + (3) Dans le mode de réalisation décrit ici, pour évaluer les éléments bm(f),nn--1,...,K du vecteur b(f) vérifiant la relation (3), le second module d'évaluation 5C2 utilise également un estimateur de Welch sur une fenêtre glissante comprenant N points (i.e. N fréquences espacées de 30 la résolution fréquentielle M précisée a-dessus). On note que chaque élément bm(f), m.=1,...,K du vecteur b(f) est un scalaire. L'étape E63 peut être réalisée indifféremment en parallèle des étapes E61-E62, avant ou après ces étapes. Le module de calcul 5C4 du dispositif de surveillance 5 calcule 35 alors le produit de la matrice inverse 1.,-)-(f) et du vecteur b(f) et obtient 23 les coefficients du filtre de Wiener dans le domaine fréquentiel (étape E64), soit 11(f) = [H1(1) H2(f) = L-1(f) b(f) (4) Cette relation peut être aisément démontrée en considérant l'équivalent de la relation (1) introduite précédemment dans le domaine fréquentiel. Plus précisément, dans le domaine fréquentiel, la relation (1) mène à la relation suivante : dOUT(f) = k(ncliNk(f + a k) k=1 avec ink(t) = in(t) exp(-2 jrra kt) et où OUT(f) et nvk(t) désignent respectivement les transformées de Fourier des signaux temporels out(t) 10 et ink(t). En multipliant la relation ci-dessus par dIN,(f+aï.) et en appliquant l'opérateur de moyennage d'ensemble E[] on obtient la relation suivante K S 11-C(f) S 17t, I it k k-1 où et S,,,,,,,kdésignent respectivement les spectres de puissance croisés entre les signaux in, et in , et entre les signaux tnt et in k. En 15 considérant comme signal d'entrée du filtre, le signal temporel x(t), il en résulte l'équation (4) utilisée par le module de calcul 5C4 pour calculer les coefficients du filtre de Wiener dans le domaine fréquentiel. En référence à la figure 3, une fois les coefficients fik(f), k=1,...,K du filtre de Wiener déterminés par le module de détermination 20 5C sur la fenêtre glissante comprenant N points, le module d'estimation 5B du dispositif de surveillance 5 estime la contribution vibratoire de la pièce mécanique 3 comprise dans le signal x(t) issue de l'étape E50 en filtrant ce signal au moyen du filtre de Wiener ainsi défini (étape E70), Dans le mode de réalisation décrit ici, le module d'estimation 58 réalise ce filtrage 25 dans le domaine fréquentiel en appliquant la relation : C3(f) = [1-ii(f) 112(f).,..HK(f)]lXi(f) X2 (f) --- - xk(f)1T où T désigne ['opérateur de transposition, C3(f) désigne la contribution vibratoire générée par la pièce mécanique 3, et Xk(f) pour k=1,...K désigne la transformée de Fourier du signal xk(t) = x(t)exp(-21ffak t), 24 Dans un autre mode de réalisation, le module d'estimation 5B applique une transformée de Fourier inverse aux coefficients fik(f), k=1,.,,,K évalués dans le domaine fréquentiel et calcule, dans le domaine temporel, la contribution vibratoire de la pièce mécanique 3 en appliquant 5 la relation (1), La contribution vibratoire ainsi extraite du signal x(t) peut ensuite être analysée, de façon connue en soi, pour déterminer si la pièce mécanique 3 est affectée d'un défaut. Dans le mode de réalisation décrit ici, elle est fournie au module 10 de restitution 12 du système 1 pour pouvoir être visualisée par un opérateur, apte à réaliser une inspection visuelle de cette contribution. En variante, on peut envisager une restitution auditive de la contribution vibratoire de la pièce mécanique 3 à l'opérateur, via par exemple un haut-parleur du module de restitution 12. 15 En variante, la contribution vibratoire de la pièce mécanique 3 peut être fournie à un module d'analyse, apte à mettre en oeuvre une détection automatique de défauts à partir du signal vibratoire qui lui est fourni. Un tel module est connu en soi et n'est pas décrit davantage ici. Il met en oeuvre par exemple l'une des techniques de détection décrites 20 dans le document de D-H. Kwak et al. « Fault Detection of Roller-Bearings Using Signal Processing and Optimization Algorithms », Sensors 2014, 14, pages 283-298 ou dans le document de F. Combet et al. « Optimal Filtering of gear signais for early damage detection based on the spectral kurtosis », Mechanical Systems and Signal Processing, vol. 23, issue 3, 25 2009, REVENDICATIONS1. Procédé de surveillance d'une pièce mécanique (3) en rotation associée à au moins une fréquence cyclique de défaut, ledit procédé de surveillance comprenant : - une étape d'obtention (E20) d'un signal vibratoire ou acoustique temporel x(t) acquis au moyen d'un capteur et incluant une contribution vibratoire ou acoustique de la pièce mécanique en rotation - une étape d'estimation (E70) de ladite contribution de la pièce mécanique en filtrant le signal x(t) au moyen d'un filtre de Wiener ; le procédé de surveillance étant caractérisé en ce qu'il comprend une étape de détermination (E60), dans un domaine fréquentiel, d'un nombre entier K de coefficients du filtre de Wiener, chaque coefficient correspondant à une dite fréquence cyclique de défaut de la pièce mécanique, l'étape de détermination comprenant : - une étape d'évaluation (E61), à au moins une fréquence f, d'une matrice LU) dont chaque élément L11([), i et I désignant deux entiers compris entre 1 et K, est égal au spectre de puissance croisé de signaux xi (t) = x(t)exp(-Zjirct,t) et xi (t) = x(t) em.)(-2jrratt), a, et al désignant des fréquences cycliques de défaut associées à la pièce mécanique en rotation ; - une étape d'évaluation (E63), à ladite au moins une fréquence f, d'un vecteur b(f) dont chaque élément bm(f), m désignant un entier compris entre 1 et K, est égal au spectre de puissance croisé d'un signal xm(t) =- x(t) exp(-2jnan,t) et du signal x(t),crm désignant une fréquence cyclique de défaut associée à la pièce mécanique ; les K coefficients du filtre de Wiener étant déterminés (E64) à partir du produit d'une matrice inverse de la matrice LU) par le vecteur b(f) évalué à ladite au moins une fréquence f. 2. Procédé de surveillance selon la revendication 1 dans lequel les éléments de la matrice L(f) vérifient L1(f) ) 3. Procédé de surveillance selon la revendication 1 ou 2 dans lequel le filtrage du signal x(t) est réalisé dans le domaine fréquentiel. 26 4. Procédé de surveillance selon l'une quelconque des revendications 1 à 3 dans lequel le nombre de coefficients K du filtre est inférieur ou égal à 4. 5. Procédé de surveillance selon l'une quelconque des revendications 1 à 4 comprenant en outre une étape (E30) de ré-échantillonnage angulaire du signal x(t), la contribution de la pièce mécanique étant estimée lors de l'étape d'estimation en filtrant le signal 10 x(t) ré-échantillonné, et les K coefficients du filtre de Wiener étant déterminés à partir du signal x(t) ré-échantillonné au lieu du signal x(t). 6. Procédé de surveillance selon l'une quelconque des revendications 1 à 5 comprenant en outre une étape de décomposition 15 (E40) du signal x(t) ou du signal x(t) ré-échantillonné en une composante déterministe et en une composante aléatoire, la contribution de la pièce mécanique étant estimée lors de l'étape d'estimation en filtrant la composante parmi la composante déterministe et la composante aléatoire qui contient la contribution vibratoire ou acoustique de la pièce 20 mécanique, et les K coefficients du filtre de Wiener étant déterminés à partir de cette composante au lieu du signal x(t). 7. Programme d'ordinateur (PROG) comportant des instructions pour l'exécution des étapes du procédé de surveillance selon l'une 25 quelconque des revendications 1 à 6 lorsque ledit programme est exécuté par un ordinateur. 8. Support d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur selon la revendication 7. 9. Dispositif de surveillance (5) d'une pièce mécanique en rotation associée à au moins une fréquence cyclique de défaut, ledit dispositif de surveillance comprenant - un module d'obtention (5A), configure pour obtenir un signal 35 vibratoire ou acoustique temporel x(t) acquis au moyen d'un capteur et 27 incluant une contribution vibratoire ou acoustique de la pièce mécanique en rotation ; - un module d'estimation (5B), configuré pour estimer ladite contribution de la pièce mécanique en filtrant le signal x(t) au moyen d'un 5 filtre de Wiener ; le dispositif de surveillance étant caractérisé en ce qu'il comprend en outre un module de détermination (5C) configuré pour déterminer dans un domaine fréquentiel un nombre entier K de coefficients du filtre de Wiener, chaque coefficient correspondant à une dite fréquence cyclique de 10 défaut associée à la pièce mécanique, ledit module de détermination comprenant - un premier module d'évaluation (bd), configuré pour évaluer à au moins une fréquence f, une matrice L(f) dont chaque élément La(f), i et I désignant deux entiers compris entre 1 et K, est égal au spectre de 15 puissance croisé de signaux x(t) = x(t)exp(-2ftrat) et x(t) = x(t)exp(-2pralt), cri et ai désignant des fréquences cycliques de défaut associées à la pièce mécanique en rotation - un second module d'évaluation (5C2), configuré pour évaluer à ladite au moins une fréquence f, un vecteur b(f) dont chaque élément bm(f), m 20 désignant un entier compris entre 1 et K, est égal au spectre de puissance croisé d'un signal xm(t) x(t)exp(-2jrcumt.) et du signal x(t) ,am désignant une fréquence cyclique de défaut associée à la pièce mécanique ; les K coefficients du filtre de Wiener étant déterminés par le module de 25 détermination (5C,5C4) à partir du produit d'une matrice inverse de la matrice L(f) par le vecteur b(f) évalué à ladite au moins une fréquence f. 10. Système (1) comprenant un capteur (4), configuré pour acquérir un signal vibratoire ou 30 acoustique généré par au moins une pièce mécanique en rotation - un dispositif de surveillance (5) selon la revendication 9, configuré pour déterminer une contribution vibratoire ou acoustique de la pièce mécanique en rotation à partir du signal vibratoire ou acoustique acquis par le capteur. 35 28 11. Système (1) selon la revendication 10 comprenant en outre un module de restitution (12) de la contribution vibratoire ou acoustique déterminée par le dispositif de surveillance. 12. Système (1) selon la revendication 11 dans lequel le module de restitution (12) comprend un écran et/ou un haut-parleur. 13. Système selon l'une quelconque des revendications 10 à 12 dans lequel le capteur (4) comprend un accéléromètre (4A) ou un 10 microphone. 14. Système de surveillance (1) selon l'une quelconque des revendications 10 à 13 dans lequel la pièce mécanique (3) en rotation est une pièce équipant un moteur à combustion ou à explosion.