Brevet FR3153425A1 — publié le 20250328

[Procédé de filtrage d’un signal] Ce procédé de filtrage d’au moins un signal à filtrer, le filtrage étant déterminé à partir de mesures représentatives du mouvement d’un arbre en rotation, comprend les étapes suivantes : - définition d’au moins un indicateur de caractérisation de l’arbre en rotation, l’indicateur de caractérisation étant obtenu à partir d’une analyse fréquentielle du comportement de l’arbre en rotation (étape 120) ; - détermination d’un gabarit de filtrage en fonction de l’indicateur de caractérisation (étape 130) ; et - application du gabarit de filtrage audit signal à filtrer (étape 140). Figure pour l’abrégé : Fig. 3 Procédé de filtrage d’un signal Domaine technique de l’invention La présente invention concerne la caractérisation d’un arbre tournant autour d’un axe de rotation. La présente invention concerne plus particulièrement le filtrage d’un indicateur de caractérisation des effets gyroscopiques et des déformations d’un arbre tournant afin d’en faciliter le contrôle, la modélisation, le monitoring, le diagnostic et la maintenance. Etat de la technique antérieure Un système tournant, tel qu’un arbre en rotation, parfaitement équilibré tourne autour d’un axe de rotation de sorte que ledit axe de rotation est immobile par rapport à son positionnement théorique, et donc par rapport à l’environnement proche de l’arbre en rotation. Cependant, des défauts, des déformations, des effets gyroscopiques, ou encore des effets dus à la force de gravité peuvent induire un décalage de l’axe de rotation par rapport à sa position théorique lors de la rotation dudit arbre en rotation. C’est le cas par exemple de déformations de l’arbre en rotation en raison de modes propres flexibles, ces modes déformant l’arbre en fonction de la vitesse de rotation dudit arbre. Lors de la présence de tels défauts, l’axe de rotation n’est pas immobile et vibre. Il s’éloigne de sa position théorique selon une trajectoire particulière. Plus précisément, chaque point de l’axe de rotation peut s’éloigner de sa position théorique, deux points de l’axe de rotation à deux extrémités différentes de l’arbre pouvant avoir une trajectoire différente. Les propriétés d’un arbre en rotation sont classiquement appréhendées par une analyse latérale de paramètres. Les paramètres sont des valeurs de grandeurs radiales telle que la position ou la vitesse radiale de l’arbre, ayant une composante non nulle dans un plan orthogonal à l’axe de rotation. Il est souvent nécessaire de caractériser les évolutions de ces grandeurs selon que l’axe soit en rotation ou non. Cependant, de multiples causes peuvent influencer le comportement en rotation d’un arbre en rotation et il est complexe de caractériser précisément un tel comportement à l’aide d’indicateurs fiables. La présente invention a donc pour but de pallier les inconvénients précités et de proposer un filtrage d’un signal à filtrer caractérisant un arbre en rotation. La présente invention a donc pour objet un procédé de filtrage d’au moins un signal à filtrer, le filtrage étant déterminé à partir de mesures représentatives du mouvement d’un arbre en rotation, le procédé comprenant les étapes suivantes : - définition d’au moins un indicateur de caractérisation de l’arbre en rotation, l’indicateur de caractérisation étant obtenu à partir d’une analyse fréquentielle du comportement de l’arbre en rotation ; - détermination d’un gabarit de filtrage en fonction de l’indicateur de caractérisation ; et - application du gabarit de filtrage audit signal à filtrer déterminé à partir de mesures représentatives du mouvement de l’arbre en rotation. Ainsi, ce procédé permet de rendre l’indicateur de caractérisation plus fiable et d’isoler certaines informations facilitant la compréhension et la rectification de défauts sur un arbre en rotation. Avantageusement, l’indicateur de caractérisation est obtenu à partir des étapes suivantes : - a) acquisition sur une période prédéterminée d’au moins deux signaux d’entrée, les signaux d’entrées comprenant des mesures représentatives du mouvement de l’arbre en rotation, les mesures, liées à une grandeur vectorielle, étant non colinéaires entre elles et optionnellement coplanaires avec un plan sécant à l’axe de rotation de l’arbre ; - b) détermination de deux composantes des signaux d’entrée, les composantes étant orthogonales entre elles et optionnellement à l’axe de rotation ; - c) fusion des deux composantes en un signal bivarié ; - d) décomposition temps-fréquence du signal bivarié en éléments spectraux ; - e) détermination d’une ellipse équivalente pour chaque élément spectral ; et - f) extraction de l’indicateur de caractérisation de l’arbre en rotation à partir de paramètres de l’ellipse équivalente. Les effets gyroscopiques et les déformations du rotor, notamment les modes flexibles, sont plus simples à considérer sous une décomposition fréquentielle des signaux d’entrée. En effet, les effets tournants de ces déformations sont asynchrones, autrement dit leurs fréquences ne sont pas directement liées à la fréquence de rotation de l’arbre, et sont donc plus facilement caractérisés à l’aide d’une décomposition fréquentielle qui peut mettre en évidence une orbite tournante, également appelée ellipse équivalente, pour chaque fréquence, chaque orbite tournant à une vitesse propre correspondant à la fréquence considérée. Cette caractérisation d’un arbre en rotation permet des applications de type contrôle, identification, modélisation, monitoring ou diagnostic. Dans un mode de mise en œuvre, le gabarit de filtrage est appliqué à la décomposition temps-fréquence du signal bivarié de manière à obtenir un signal bivarié filtré. Dans divers modes de mise en œuvre, la décomposition temps-fréquence du signal bivarié sur laquelle est appliqué le gabarit de filtrage est une décomposition temps-fréquence obtenue par transformée de Fourier, ou par ondelette, ou par filtrage, ou par transformée de Fourier quaternionique. Avantageusement, le procédé comprend une étape d’obtention de signaux d’entrée filtrés par projection du signal bivarié filtré sur les axes colinéaires aux grandeurs mesurées. Dans un mode particulier de mise en œuvre, le procédé comprend une étape de calcul d’un signal dual du signal bivarié filtré. Dans un mode de mise en œuvre, l’au moins un indicateur comprend un taux de polarisation, et/ou un indicateur de circularité et/ou d’amplitude et/ou de direction et/ou d’orientation et/ou de synchronisation d’une ellipse équivalente. Avantageusement, le gabarit de filtrage est défini par rapport à un seuil prédéfini d’un taux de polarisation, et/ou d’un indicateur de circularité d’une ellipse équivalente. Avantageusement, le gabarit de filtrage est défini à partir d’au moins un indicateur normalisé compris entre -1 et 1 ou entre 0 et 1. Dans un mode particulier de mise en œuvre, le gabarit de filtrage est un filtre bivarié linéaire et/ou est obtenu par multiplication d’un masque déterminé en fonction de l’indicateur de caractérisation et d’une pondération fréquentielle. Dans un mode de mise en œuvre, les étapes a) à f) sont effectuées au moins une deuxième fois pour un plan additionnel sécant à l’axe de rotation de l’arbre et distinct dudit plan, le procédé comprenant en outre une étape de comparaison des indicateurs de caractérisation des plans. La présente invention a également pour objet une méthode mise en œuvre par ordinateur comprenant les étapes du procédé telle que définies précédemment. L’implémentation peut se faire sur une cible temps réel, notamment dans une armoire de contrôle de palier magnétique. Brève description des figures FIG. 1 est une représentation schématique des différentes étapes d’un procédé d’obtention d’un indicateur pour le procédé selon l’invention ; FIG. 2 est une représentation d’une ellipse et de ses paramètres d’Euler ; et FIG. 3 est une représentation schématique de différentes étapes de filtrage applicables à un signal à filtrer suite aux étapes de la FIG. 1 . Description détaillée de l’invention On a représenté schématiquement sur la FIG. 1 différentes étapes d’obtention d’un indicateur pouvant, dans un mode de mise en œuvre particulier, être incorporées dans le procédé selon l’invention. Ces étapes permettent de caractériser le comportement d’un arbre en rotation. On entend par arbre en rotation un arbre mobile en rotation autour d’un axe de rotation. L’arbre en rotation est par exemple un arbre d’une machine tournante. Pour la mise en œuvre du procédé, on effectue en premier lieu une étape 10 d’acquisition de signaux d’entrée. Au moins deux signaux d’entrées sont acquis sur une période prédéterminée. La période prédéterminée est par exemple de plusieurs secondes. Les signaux d’entrée comprennent des mesures représentatives du mouvement de l’arbre en rotation. Par exemple, les signaux d’entrée comprennent des mesures de position, ou de vitesse, ou d’accélération, ou de force, ou de champ magnétique, ou encore d’induction magnétique de l’arbre en rotation. L’acquisition des signaux d’entrée s’effectue par exemple à l’aide de capteurs configurés pour mesurer ces grandeurs. D’une manière générale, une grandeur vectorielle est mesurée selon plusieurs composantes réparties dans l’espace sensibles à la rotation de l’arbre. La répartition dans l’espace de ces composantes permet d’évaluer le mouvement du rotor. Un cas particulier est celui de mesures non colinéaires entre elles et coplanaires avec un plan sécant à l’axe de rotation de l’arbre. A partir de deux mesures non colinéaires, une caractérisation du comportement de l’arbre en rotation est possible dans un plan sécant de l’axe de rotation de l’arbre. Des mesures supplémentaires permettent une plus grande précision de ladite caractérisation. Des mesures non coplanaires avec un plan sécant à l’axe de rotation de l’arbre peuvent néanmoins être ramenées par projection dans un plan de référence sécant à l’axe de rotation de l’arbre. Pour M mesures avec M≥2, on obtient M signaux d’entrées à , t désignant le temps durant l’acquisition des mesures, dont l’ensemble des valeurs est désigné par T, autrement dit dates, que peut prendre t. On entend par plan sécant un plan dans lequel sont faites les mesures qui coupe l’axe de rotation de l’arbre en rotation. Le plan sécant est par exemple orthogonal à l’axe de rotation ou non orthogonal. Dans le cas où le plan sécant n’est pas orthogonal à l’axe de rotation, on effectue une étape 20 de détermination de deux composantes des signaux d’entrée, les composantes étant orthogonales entre elles et optionnellement à l’axe de rotation. Autrement dit, cette étape 20 est un changement de base facilitant la suite des calculs. Dans le cas où le plan sécant est orthogonal à l’axe de rotation, l’étape précédente est intrinsèquement effectuée. Dans le cas où l’axe de rotation n’est pas rigoureusement et précisément défini, l’arbre n’étant pas toujours parfaitement cylindrique et/ou équilibré, et/ou indéformable et/ou connu, ledit plan sécant orthogonal est assimilé à un plan d’étude supposé orthogonal à l’axe de rotation, ou au minimum sensible au mouvement de rotation de l’arbre. En pratique, on effectue l’étape 20 de détermination des deux composantes orthogonales et en moyennant des projections des signaux d’entrées à sur deux axes orthogonaux x et y, ces deux axes étant idéalement orthogonaux à l’axe de rotation. Ainsi, Les composantes et sont donc désormais dans le plan d’étude, idéalement orthogonal à l’axe de rotation. Puis, on effectue une étape 30 de fusion des deux composantes orthogonales et en un signal bivarié afin de réunir ces deux composantes de manière équivalente sous forme d’un signal complexe, telle que : On peut également l’exprimer de la manière suivante : Ou encore sous forme vectorielle : On effectue ensuite une étape 40 de décomposition fréquentielle du signal bivarié en éléments spectraux caractérisés par leur couple temps-fréquence , d’ailleurs également appelés cellules temps-fréquence, et ce périodiquement afin d’obtenir une décomposition à la fois temporelle et fréquentielle, dite décomposition temps-fréquence telle que : Il existe différentes manières d’effectuer cette étape 40 de décomposition temps-fréquence. Ces différents modes de mise en œuvre sont présentés ci-après, les décompositions temps-fréquence , ou , ou , ou étant données sous forme continue, autrement dit analogique et non discrètes. Une implémentation numérique mise en œuvre par ordinateur comprend alors une étape d’échantillonnage en temps et en fréquence. Dans un premier mode de mise en œuvre, noté 41 sur la FIG. 1 , l’étape 40 de décomposition temps-fréquence est réalisée avec un banc de filtre autour de chaque fréquence f. Ainsi, la décomposition temps-fréquence s’exprime : avec la réponse impulsionnelle d’un filtre passe-bande centré autour de la fréquence f, étant l’ensemble des fréquences supérieures ou égales à 0. Dans un deuxième mode de mise en œuvre 42, l’étape 40 de décomposition temps-fréquence est réalisée par transformée de Fourier, en particulier transformée de Fourier simple face à court terme. De base, la décomposition temps-fréquence s’exprime : Dans la pratique, on utilise de préférence des estimateurs plus performants reposant sur un découpage en intervalles temporels, avec pondérations temporelles, moyennage sur plusieurs intervalles temporels, chevauchements entre intervalles successifs. La formulation est alors : Avec et pour définissant un découpage en p intervalles, avec chevauchement éventuel, autour d’une date donnée t . Classiquement, ce découpage ne dépend pas de t , identique pour chaque cellule. L’étendue temporelle de l’intervalle k est . L’étendue temporelle globale de la cellule ( t,f ) correspondante est : pour pour définissant une fenêtre de pondération temporelle pour chacun de ces intervalles, par exemple de type Hanning. Classiquement, ne dépend pas de k , identique pour chaque intervalle. Selon la grandeur spectrale estimée, il est nécessaire d’introduire les compensations suivantes : Dans un cas particulier, l’absence de pondération correspond à =1, auquel cas . De plus, on introduit un sous-ensemble de dates : Dans un troisième mode de mise en œuvre 43, l’étape 40 de décomposition temps-fréquence est réalisée par transformée de Fourier, en particulier transformée de Fourier double face à court terme, également appelée « Fullspectrum » en termes anglo-saxons, autrement dit pour un ensemble F de fréquences prenant des valeurs dans l’ensemble des nombres réels. Ainsi, de base la décomposition temps-fréquence s’exprime : Ou avec un découpage en blocs temporels : Dans un quatrième mode de mise en œuvre 44, l’étape 40 de décomposition temps-fréquence est réalisée par transformée de Fourier, en particulier transformée de Fourier quaternionique à court terme, faisant intervenir les nombres imaginaires unitaires vérifiant 1i 2 = 1j 2 = 1k 2 = 1i.1j.1k = -1. Ainsi, de base la décomposition temps-fréquence s’exprime : Ou avec un découpage en blocs temporels : Les expressions des transformées de Fourier données dans les modes de mise en œuvre 41 à 44 ci-dessus constituent des formules applicables, réalisant des estimateurs performants notamment dans leur formulation avec pondérations temporelles, moyennage sur plusieurs intervalles temporels, chevauchements entre intervalles successifs, dont on peut donner une expression plus générale : étant la pondération temporelle, et étant défini tel que . Finalement, L’indice k appliqué à w k (ci-dessus) ou p k (ci-dessous) signifie que la pondération peut différer d’un intervalle à l’autre. De plus, d’autres transformées fréquentielles peuvent être utilisées, telles que les décompositions en ondelettes. Une expression plus générale est : Plus généralement encore, il est possible de faire varier le découpage temporel et la pondération en fonction de la cellule temps-fréquence : La trajectoire (orbite) s’obtient par transformation inverse, de type : L’étendue fréquentielle de la cellule étant alors , l’étendue temporelle de la cellule (t, f) étant . De chaque élément spectral , à une date t et une fréquence f fixés, on peut déduire de l’orbite instantanée une ellipse équivalente. On peut la qualifier de locale au sens temporel et fréquentiel, autrement dit à la date t et à la fréquence f. On effectue donc une étape 50 de détermination de l’ellipse équivalente pour chaque élément spectral. Quelle que soit la méthode de décomposition temps-fréquence, l’orbite instantanée peut s’exprimer sous la forme d’une équation paramétrique telle que, pour chaque composante x et y, et avec , et les amplitudes et phases instantanées : Pour une fréquence f fixée, parcourir le temps t permet d’obtenir la trajectoire à cette fréquence. Cette orbite instantanée prend la forme d’une ellipse équivalente, pour peu que , ), , varient lentement en fonction du temps, par exemple durant une révolution complète à la fréquence f , et dans tous les cas pendant la durée dt de l’intervalle temporel définissant l’élément spectral. Dans ce cas, on peut distinguer la date t de l’élément spectral considéré du temps qui évolue dans cet élément spectral et. De plus, en introduisant la fréquence f dans l’expression des phases supposées quasi-linéaires, on obtient l’ellipse équivalente sous forme d’équation paramétrique temporelle pour chaque composante x et y : Une autre écriture de ces équations paramétriques est par exemple : Pour obtenir ces équations paramétriques, notamment dans le cas d’une décomposition fréquentielle par banc de filtres, les paramètres , , , de l’ellipse équivalente peuvent être obtenus par estimations, tels que moyennage direct ou quadratique ou par valeur efficace, à partir des grandeurs instantanées . Dans le cas d’une décomposition fréquentielle par transformée de Fourier, on obtient l’équation paramétrique de l’ellipse équivalente par transformée de Fourier inverse de l’élément fréquentiel : En formant la quantité complexe associée , autrement dit la représentation bivariée, on met en évidence la décomposition d’une équation paramétrique de l’ellipse équivalente en deux cercles de directions opposées : Avec les paramètres suivants obtenus suite à une étape 59 de calcul : La caractérisation de l’ellipse équivalente est naturelle avec la transformée de Fourier double face (« fullspectrum ») qui peut directement s’écrire sous forme polaire : De manière générale, une ellipse équivalente peut être déterminée et tracée par l’obtention de ses paramètres d’Euler. On a représenté sur la figure 2 un exemple d’une ellipse ainsi que ses paramètres d’Euler , , , qui définissent complètement une ellipse équivalente, avec : Par ailleurs, b est le demi-grand axe, autrement dit le grand rayon, et c est le demi-petit axe, autrement dit le petit rayon, dont le signe donne le sens de rotation de l’ellipse équivalente. L’équation paramétrique de la trajectoire de l’ellipse équivalente est plus commode sous forme complexe, car : - on exprime une équation de cercle, notamment en cos et sin de la même phase, - on fait subir une dilatation des axes pour obtenir une ellipse, - on lui applique une rotation pour obtenir la direction voulue, ainsi : On entend par direction l’angle θ de l’ellipse équivalente entre le grand axe et l’axe x. On entend par orientation le sens de rotation décrivant la trajectoire de ladite ellipse. Cette caractérisation de l’ellipse équivalente est naturelle avec la décomposition temps-fréquence par transformée de Fourier quaternionique qui peut directement se mettre sous forme d’Euler : Cette équation paramétrique a l’intérêt particulier de rejeter la paramétrisation temporelle vers un seul terme exponentiel lorsque l’on applique la transformée de Fourier inverse à l’élément fréquentiel afin d’obtenir l’équation paramétrique temporelle de la trajectoire de l’ellipse équivalente : avec une projection dans l’ensemble des nombres complexes . De manière plus précise, l’étape 50 de détermination de l’ellipse équivalente comprend une étape de calcul de paramètres d’Euler de ladite ellipse équivalente. En partant de la décomposition temps-fréquence par banc de filtre, on obtient une étape de calcul 51 de et telle que : étant un sous-ensemble de T. En partant de la décomposition temps-fréquence par transformée de Fourier simple face à court terme, on obtient une étape 52 de calcul de , , et telle que : En partant de la décomposition temps-fréquence par transformée de Fourier double face à court terme, on obtient une étape 53 de calcul de , , et telle que : Par une étape 54 ou 55 supplémentaire de calcul, on obtient les paramètres d’Euler suivants : L’étape 56 de calcul inverse est également possible, tel que : Le procédé comprend optionnellement une étape 60 de calcul de la densité spectrale de la décomposition temps-fréquence. Les densités spectrales de chaque composante et croisées entre composantes permettent d’estimer la répartition énergétique entre les composantes x et y . En partant de la décomposition temps-fréquence par banc de filtre, on obtient une étape 61 de calcul de la densité spectrale telle que : En partant de la décomposition temps-fréquence par transformée de Fourier simple face à court terme, on obtient une étape 62 de calcul de la densité spectrale telle que : En partant de la décomposition temps-fréquence par transformée de Fourier double face à court terme, on obtient une étape 63 de calcul de la densité spectrale telle que : Nous utiliserons aussi : En partant de la décomposition temps-fréquence par transformée de Fourier quaternionique à court terme, on obtient une étape 64 de calcul de la densité spectrale telle que : A partir de ces densité spectrales , , et , on effectue une étape 70 de calcul des paramètres de Stokes , , , et d’au moins une ellipse équivalente. Les paramètres de Stokes sont quatre nombres réels et sont classiquement utilisés pour décrire l’état de polarisation d’une onde électro-magnétique. S 0 est communément l’intensité totale d’une onde électro-magnétique et représente ici par extension l’intensité totale du signal bivarié, S 0 >0 . S 1 et S 2 donnent l’intensité de polarisation rectiligne ainsi que la direction de la polarisation de l’ellipse . S 3 donne l’intensité de polarisation circulaire , son signe donnant le sens de rotation. On définit également les paramètres de Stokes réduits : En partant de la densité spectrale obtenue pour la décomposition temps-fréquence par banc de filtre, on obtient une étape 71 de calcul des paramètres de Stokes d’au moins une ellipse équivalente telle que : En partant de la densité spectrale obtenue pour la décomposition temps-fréquence par transformée de Fourier simple face à court terme, on obtient une étape 72 de calcul des paramètres de Stokes d’au moins une ellipse équivalente telle que : En partant de la densité spectrale obtenue pour la décomposition temps-fréquence par transformée de Fourier double face à court terme, on obtient une étape 73 de calcul des paramètres de Stokes d’au moins une ellipse équivalente telle que : En partant de la densité spectrale obtenue pour la décomposition temps-fréquence par transformée de Fourier quaternionique à court terme, on obtient une étape 74 de calcul des paramètres de Stokes d’au moins une ellipse équivalente telle que : En variante, les paramètres de Stokes , , et sont calculés directement à l’aide de la décomposition temps-fréquence par transformée de Fourier quaternionique : A partir de ces paramètres de Stokes, on effectue une étape 75 de calcul du taux de polarisation ainsi que des paramètres de Stokes réduits précédemment évoqués, ainsi que la partie polarisée de l’intensité . La polarisation exprimée par le taux de polarisation permet de mettre en évidence la cyclo-stationnarité d’une ellipse équivalente, à savoir la régularité de la trajectoire de ladite ellipse par rapport à une période de révolution. La polarisation reflète la structuration de ladite trajectoire relativement à des fréquences définies. Le taux de polarisation peut ainsi être calculé pour chaque cellule temps-fréquence issue de la décomposition temps-fréquence afin de mettre en évidence la cyclo-stationnarité de l’arbre en rotation pour une fréquence donnée. Par ailleurs, une étape 57 de calcul peut être mise en œuvre pour déterminer des paramètres d’Euler , , et à partir de la densité spectrale de polarisation et des paramètres de Stokes réduits : Un quatrième paramètre d’Euler donnant la phase à l’origine de la trajectoire peut également être déduit par l’étape 58 de calcul suivante : Suite à l’étape 50 de détermination d’une ellipse équivalente, notamment par le calcul de paramètres d’Euler de ladite ellipse équivalente, on effectue une étape 80 d’extraction d’au moins un indicateur de caractérisation de l’arbre en rotation à partir de paramètres de l’ellipse équivalente, par exemple des paramètres d’Euler, de Stokes ou des paramètres dérivés de ces derniers. En outre, on effectue optionnellement les étapes précédentes au moins une deuxième fois pour un plan additionnel sécant à l’axe de rotation de l’arbre et distinct du plan pour lequel ces étapes ont été réalisées en premier lieu. Dans ce mode de mise en œuvre, on peut ainsi effectuer une étape de comparaison des indicateurs de caractérisation des plans. Par exemple, l’au moins un indicateur de caractérisation comprend le taux de polarisation précédemment décrit, et/ou un indicateur de circularité (lié à ), et/ou un indicateur d’intensité et/ou un indicateur d’amplitude, et/ou un indicateur de sens de rotation (lié au signe de ), et/ou un indicateur de direction (lié à ), et/ou de synchronisation (lié à ) de l’ellipse équivalente. Chaque indicateur élémentaire contribue à une caractéristique concrète et utile du système en rotation, par exemple un fort taux de polarisation et une forte circularité sont probablement du fait de l’arbre en rotation. Distinguer cette partie du signal, le complémentaire étant assimilé à des perturbations, peut être utile pour par exemple identifier ledit arbre (à des fins de modélisation) ou mieux le contrôler dans le cadre d’un arbre équipé de paliers magnétiques. Optionnellement, l’au moins un indicateur de caractérisation est normalisé, autrement dit, sa valeur est comprise entre -1 et 1 ou entre 0 et 1. L’indicateur de caractérisation est donné pour un plan, il est alors noté A, ou par comparaison entre deux plans, il est alors noté I. L’indicateur de caractérisation est qualifié de critère, noté C, lorsqu’il est construit à partir d’une grandeur caractéristique quelconque réévaluée à chaque instant t. L’indicateur de caractérisation est qualifié de stationnarité, noté S, lorsqu’il est basé sur la stationnarité dans le temps d’une grandeur (maximum lorsque la variation dans le temps de ladite grandeur est nulle). Ainsi, on note l’indicateur du taux de polarisation : Pour un plan, l’indicateur d’intensité est : Avec la fonction de normalisation : TT et FF désignent des ensembles de dates et de fréquences génériques, prenant leur valeur dans l’ensemble des réels R. Pour un plan, l’indicateur d’amplitude est : Pour un plan, l’indicateur de signe est-1 ou 1 selon le signe de : Pour un plan, l’indicateur de circularité est : Pour un plan, l’indicateur de stationnarité de signe, autrement dit de sens de rotation, autrement dit d’orientation, est : Avec la fonction de stationnarité : Pour un plan, l’indicateur de stationnarité de circularité est : Pour un plan, l’indicateur de stationnarité de direction, autrement dit d’inclinaison géométrique de l’ellipse, est : Pour un plan, l’indicateur de stationnarité de synchronisation est : Pour deux plans, on redéfinit chaque paramètre ou indicateur pour un plan additionnel, en notant ces nouveaux paramètres à l’aide d’une apostrophe, par exemple . L’indicateur de caractérisation comparative de l’amplitude est : L’indicateur de caractérisation comparative du signe est un booléen entre 0 et 1 selon le résultat du test suivant : L’indicateur de caractérisation comparative de circularité est : Avec la fonction d’erreur relative telle que : Et avec la fonction division avec exception telle que : L’indicateur de caractérisation comparative de direction, autrement dit d’inclinaison géométrique de l’ellipse, est : L’indicateur de caractérisation comparative de synchronisation est : L’indicateur de caractérisation de stationnarité et comparative de circularité est : L’indicateur de caractérisation de stationnarité et comparative de direction est : L’indicateur de caractérisation de stationnarité et comparative de synchronisation est : D’autres indicateurs peuvent être construits à partir des grandeurs construites dans les étapes 70 et 50. Optionnellement, on effectue une étape 90 d’insertion de l’au moins un indicateur de caractérisation dans un réseau de neurones afin de classifier ledit indicateur en fonction de la rotation de l’arbre et le relier à une anomalie préalablement répertoriée. Les anomalies similaires à la base d’apprentissage sont ensuite détectées. Optionnellement, on effectue une étape 100 d’affichage de l’au moins un indicateur de caractérisation extrait pour être lisible par un opérateur et permettre sa compréhension de potentielles anomalies de l’arbre en rotation. Une étape 110 d’affichage d’une image de la trajectoire d’au moins une ellipse équivalente peut également être effectuée. On a représenté schématiquement sur la FIG. 3 les étapes du procédé de filtrage selon un exemple de l’invention. Dans un mode particulier de mise en œuvre, les étapes suivantes sont réalisées à la suite de l’étape 80 d’extraction d’un indicateur de caractérisation. Le procédé de filtrage est destiné à filtrer au moins un signal à filtrer, par exemple un signal déterminé à partir de mesures représentatives du mouvement d’un arbre en rotation. Le procédé est réalisé en temps réel ou non. Ledit signal peut être l’un des signaux d’entrée déterminés à l’étape 10, et/ou l’une des composantes des signaux d’entrée déterminée à l’étape 20, et/ou encore le signal bivarié déterminé à l’étape 30. Afin d’effectuer ce procédé de filtrage, on effectue en premier lieu une étape 120 de définition d’au moins un indicateur de caractérisation de l’arbre en rotation, l’indicateur de caractérisation étant obtenu à partir d’une analyse fréquentielle du comportement de l’arbre en rotation. Dans un mode particulier de mise en œuvre, l’indicateur de caractérisation est un indicateur extrait durant l’étape 80 précédemment décrite. Ainsi, l’au moins un indicateur de caractérisation comprend le taux de polarisation précédemment décrit, et/ou un indicateur de circularité, et/ou un indicateur d’intensité et/ou un indicateur d’amplitude, et/ou un indicateur de sens de rotation, et/ou un indicateur de direction, et/ou de synchronisation de l’ellipse équivalente. On effectue ensuite une étape 130 de détermination d’un gabarit de filtrage en fonction de l’au moins un indicateur de caractérisation. L’étape 130 de détermination d’un gabarit de filtrage comprend par exemple une étape de calcul 131 d’un masque M, tel que : Où combine est une fonction qui prend en entrée un ou plusieurs indicateurs de caractérisation à valeurs dans [0,1] ou [-1,1] et restitue en sortie une valeur réelle dans [0,1]. Suite à l’étape 131 de calcul du masque M, on effectue une étape de calcul 132 du gabarit de filtrage P et son dual P* tels que : Avec N une pondération fréquentielle à valeurs dans [0,1]. N est une pondération fréquentielle arbitraire, indépendante de M et s’apparente à un filtrage fréquentiel classique. N permet par exemple de se focaliser sur une bande de fréquence autour de la vitesse de l’arbre. Et le masque dual (complémentaire) : L’intérêt de considérer M* et par suite P* est que selon l’application, on peut s’intéresser à une certaine partie du signal, obtenue selon les indicateurs de caractérisation utilisés, donnée par M et/ou au reliquat assimilé aux perturbations donné par M*. Par exemple, on peut imaginer un contrôleur de paliers magnétiques agissant différemment sur la partie tournante et sur les perturbations extérieures. Dans un mode de mise en œuvre, le masque M, et par extension le gabarit de filtrage P, est défini par rapport à un seuil prédéfini d’un indicateur de caractérisation, par exemple un seuil prédéfini d’un taux de polarisation, et/ou d’un indicateur de circularité d’une ellipse équivalente telle que décrite précédemment. La constance d’un indicateur de caractérisation permet également de définir un gabarit de filtrage par rapport audit indicateur de caractérisation. On effectue finalement une étape 140 d’application du gabarit de filtrage P au signal à filtrer, par exemple un signal déterminé à partir des mesures représentatives du mouvement de l’arbre en rotation. Dans un mode de mise en œuvre, le gabarit de filtrage est appliqué à la décomposition temps-fréquence du signal bivarié de manière à obtenir un signal bivarié filtré R et son dual R* calculés conjointement. L’étape 140 est alors effectuée selon différents modes de mise en œuvre. En partant de la décomposition temps-fréquence par transformée de Fourier simple face à court terme, on obtient une étape 141 de calcul d’un signal bivarié filtré intermédiaire , de son dual , et des composantes des signaux d’entrée déterminées à l’étape 20 filtrées , , , et . et En partant de la décomposition temps-fréquence par transformée de Fourier double face à court terme, on obtient une étape 142 de calcul du signal bivarié filtré intermédiaire , de son dual , tels que : En partant de la décomposition temps-fréquence par transformée de Fourier quaternionique à court terme, on obtient une étape 143 de calcul du signal bivarié filtré intermédiaire , de son dual , tels que : Enfin, on l’égalité 144 suivante permet de déterminer le signal bivarié filtré R et son dual R* en reconstituant l’axe temporel complet par morceaux : Optionnellement, on effectue une étape 150 d’obtention de signaux d’entrée filtrés à , par projection du signal bivarié filtré R sur les axes colinéaires aux grandeurs mesurées. Le gabarit de filtrage permet ainsi de sélectionner des plages de fréquences spécifiques sur les signaux d’entrée. Dans un mode particulier de mise en œuvre, on applique plusieurs gabarits de filtrage successivement et/ou en combinaison. De plus, le procédé de filtrage peut être appliqué à des mesures effectuées dans des plans différents ou pour des grandeurs différentes. Dans un mode particulier de mise en œuvre, toutes les étapes du procédé sont mises en œuvre par ordinateur, avec possibilité d’exécution en temps réel, notamment dans une armoire de contrôle de palier magnétique. Certains termes et fonctions utilisés dans les calculs précédemment détaillés sont redéfinis ci-après : L’ensemble des fréquences strictement positives est : L’ensemble des fréquences strictement négatives est : On définit également : On définit également : L’ensemble complet des fréquences est : La fonction générique de stationnarité est : La fonction de maximum avec exception est : La fonction de normalisation est : La fonction de division avec exception est : La fonction d’erreur relative est : Rappelons également la définition de la fonction arctan2, formulation équivalente à l’argument du nombre complexe : Procédé de filtrage d’au moins un signal à filtrer, le filtrage étant déterminé à partir de mesures représentatives du mouvement d’un arbre en rotation, caractérisé en ce qu’il comprend les étapes suivantes : définition d’au moins un indicateur de caractérisation de l’arbre en rotation, l’indicateur de caractérisation étant obtenu à partir d’une analyse fréquentielle du comportement de l’arbre en rotation (étape 120) ; détermination d’un gabarit de filtrage en fonction de l’indicateur de caractérisation (étape 130) ; et application du gabarit de filtrage audit signal à filtrer (étape 140). Procédé selon la revendication 1, dans lequel l’indicateur de caractérisation est obtenu à partir des étapes suivantes : a) acquisition sur une période prédéterminée d’au moins deux signaux d’entrée, les signaux d’entrées comprenant des mesures représentatives du mouvement de l’arbre en rotation, les mesures étant non colinéaires entre elles (étape 10) ; b) détermination de deux composantes des signaux d’entrée, les composantes étant orthogonales entre elles (étape 20) ; c) fusion des deux composantes en un signal bivarié (étape 30) ; d) décomposition temps-fréquence du signal bivarié en éléments spectraux (étape 40) ; e) détermination d’une ellipse équivalente pour chaque élément spectral (étape 50) ; et f) extraction de l’indicateur de caractérisation de l’arbre en rotation à partir de paramètres de l’ellipse équivalente (étape 80). Procédé selon la revendication 2, dans lequel le gabarit de filtrage est appliqué à la décomposition temps-fréquence du signal bivarié de manière à obtenir un signal bivarié filtré. Procédé selon la revendication 3, dans lequel la décomposition temps-fréquence du signal bivarié sur laquelle est appliqué le gabarit de filtrage est une décomposition temps-fréquence obtenue par transformée de Fourier, ou par ondelette, ou par filtrage, ou par transformée de Fourier quaternionique. Procédé selon l’une des revendications 3 et 4, comprenant en outre une étape (150) d’obtention de signaux d’entrée filtrés par projection du signal bivarié filtré sur les axes colinéaires aux grandeurs mesurées. Procédé selon l’une quelconque des revendications 3 à 5, comprenant en outre une étape de calcul d’un signal dual du signal bivarié filtré. Procédé selon l’une quelconque des revendications 2 à 6, dans lequel l’au moins un indicateur comprend un taux de polarisation, et/ou un indicateur de circularité et/ou d’amplitude et/ou de direction et/ou d’orientation et/ou de synchronisation d’une ellipse équivalente. Procédé selon l’une quelconque des revendications 2 à 7, dans lequel le gabarit de filtrage est défini par rapport à un seuil prédéfini d’un taux de polarisation, et/ou d’un indicateur de circularité d’une ellipse équivalente. Procédé selon l’une quelconque des revendications 1 à 8, dans lequel le gabarit de filtrage est défini à partir d’au moins un indicateur normalisé compris entre -1 et 1 ou entre 0 et 1. Procédé selon l’une quelconque des revendications 1 à 9, dans lequel le gabarit de filtrage est un filtre bivarié linéaire et/ou est obtenu par multiplication (étape 131) d’un masque déterminé en fonction de l’indicateur de caractérisation, et d’une pondération fréquentielle.