Brevet FR3005180A1 — publié le 20141031

B12549 - B232 - DI 6180-01 1 SYSTEME DE DETECTION D'ANOMALIE Domaine de la présente description La présente description concerne un procédé et un système pour détecter une anomalie dans un système mécanique ou électromécanique, et en particulier un procédé et un système pour identifier des composantes spectrales harmoniques et/ou de modulation dans un signal provenant d'au moins un capteur détectant des vibrations et/ou des fluctuations électriques dans le système. Arrière-plan Il est connu d'utiliser une analyse de signal sur un signal provenant d'un capteur de vibrations pour surveiller et diagnostiquer des anomalies dans des machines. Par exemple, un procédé pour la surveillance de l'état d'hélicoptères est décrit dans la publication intitulée "Condition Monitoring Diagnosis Methods of Helicopter Units", L.M. Gelman et al., Mechanical Systems and Signal Processing 14.4 (July 2000). En outre, la surveillance de l'état d'éoliennes est décrite dans la publication intitulée "Condition monitoring of wind turbines : Technique and methods", Pedro Fausto Garcia Marquez et al' Renewable Energy 46 (2012), pages 169-175. Dans de telles techniques, on utilise un ou plusieurs capteurs pour mesurer des vibrations, et le signal provenant de B12549 - B232 - DI 6180-01 2 chaque capteur est converti en un spectre de fréquence afin d'analyser la distribution fréquentielle. En particulier, .1a présence de fréquences fondamentales et harmoniques, et/ou de bandes latérales de modulation à certaines fréquences, peut indiquer une anomalie dans le système, comme l'usure ou le dysfonctionnement de certaines parties. Toutefois, l'efficacité des techniques courantes pour réaliser une telle détection d'anomalies est limitée en raison d'un problème technique pour identifier avec précision les 10 fréquences fondamentales, harmoniques et de bande de modulation dans le signal provenant des capteurs. Résumé Un objet de modes de réalisation de la présente description est de résoudre au moins partiellement un ou 15 plusieurs problèmes de l'art antérieur. Selon un aspect, on prévoit un procédé pour détecter une anomalie dans un système mécanique ou électromécanique, comprenant les étapes suivantes : identifier par un dispositif de traitement, à partir d'un ensemble de composantes spectrales 20 généré sur la base d'un signal provenant d'au moins un capteur détectant des vibrations et/ou des fluctuations électriques dans le système, au moins des première et deuxième composantes spectrales, chaque composante représentant un harmonique ou une modulation associé à une composante candidate : en déterminant, 25 sur la base d'une fréquence centrale et d'une marge d'erreur de la composante candidate, des limites de fréquence inférieure et supérieure d'une première fenêtre de recherche pour identifier la première composante ; en identifiant la première composante sur la base de la première fenêtre de recherche ; en modifiant 30 la fréquence centrale et la marge d'erreur de la composante candidate sur la base d'une fréquence centrale et d'une marge d'erreur de la première composante ; et en déterminant, sur la base de la fréquence centrale modifiée et de la marge d'erreur modifiée de la composante candidate, des limites de fréquence B12549 - B232 - DI 6180-01 3 inférieure et supérieure d'une deuxième fenêtre de recherche pour identifier la deuxième composante. Selon un mode de réalisation, le procédé comprend en outre l'identification de la deuxième composante dans la deuxième fenêtre de recherche. Selon un mode de réalisation, le procédé comprend en outre la définition d'une série de composantes spectrales comprenant au moins la composante candidate et les première et deuxième composantes, et la détection de l'anomalie au moins 10 partiellement sur la base de ladite série. Selon un mode de réalisation, le procédé comprend en outre la génération d'un signal de sortie pour contrôler une partie du système mécanique ou électromécanique sur la base de l'anomalie détectée. 15 Selon un mode de réalisation, les première et deuxième composantes sont des harmoniques de la composante candidate, et l'anomalie est détectée sur la base de paramètres d'une ou plusieurs composantes de modulation associées à la composante candidate ou à la première ou la deuxième composante. 20 Selon un mode de réalisation, chacune des limites inférieure et supérieure des première et deuxième fenêtres de recherche est déterminée sur la base des équations rk(Vi-Ei) et rk(Vi+Ei), où rk est l'ordre harmonique des première et deuxième composantes respectivement, Vi est la fréquence centrale 25 associée à la composante candidate, et Ei est la marge d'erreur associée à la composante candidate. Selon un mode de réalisation, la première composante est considérée comme tombant dans la première fenêtre de recherche si la condition suivante est vérifiée : 30 [a; = - Ei;Vi + EAn h(vi - El); k(Vi + = 0 où Vj est la fréquence centrale de la première composante, Ej est la marge d'erreur de la première composante, Vi est la fréquence centrale de la composante candidate, Ei est la marge d'erreur de la composante candidate et rk est l'ordre de 35 l'harmonique. 2549 - 5232 - DI 6180-01 Selon un mode de réalisation, à la suite de l'identification de la première composante, la fréquende centrale Vi et la marge d'erreur Ei sont mises à jour de la façon suivante : = (a + b)Ei (b -a) 2 x r ; 2 x r Selon un mode de réalisation, le procédé comprend en outre la détermination d'au moins un élément parmi une valeur de densité Di de la série, une valeur de richesse Ri de la série et un taux de distorsion THDi de la série sur la base des équations 10 suivantes : Di = card(Hi) NrraX ill = rua x max r. 2 with MaX Ni AI + A2 + ...+ AN THDi = -2 Ao où Di est la densité de la série, rimax est le plus grand 15 rapport rk dans la série Hi, card(Hi) est la cardinalité (taille) de la série Hi, Ri est la richesse de la série, Nimax est la taille maximum de la série sur la base de la fréquence fondamentale Vi de la série et de la fréquence maximum Vmax dans le spectre, Lxi est l'entier le plus grand qui n'est pas 20 supérieur à x, THDi est le taux de distorsion dans la série Hi, sur la base des amplitudes Al à AN des composantes 1 à N dans la série. Selon un autre aspect, on prévoit un support lisible par un ordinateur mémorisant un programme d'ordinateur qui, 25 lorsqu'il est exécuté par un dispositif de traitement, amène la mise en oeuvre du procédé susmentionné. Selon un autre aspect, on prévoit un dispositif de calcul pour détecter une anomalie dans un système mécanique ou électromécanique comprenant un dispositif de traitement agencé 2549 - 3232 - DI 6180-01 pour : identifier, à partir d'un ensemble de composantes spectrales généré sur la base d'un signal provenant d'au moins un capteur détectant des vibrations et/ou des fluctuations électriques dans le système, au moins des première et deuxième 5 composantes spectrales, chaque composante représentant un harmonique ou une modulation associé à une composante candidate : en déterminant, sui la base d'une fréquence centrale et d'une marge d'erreur de la composante candidate, des limites de fréquence inférieure et supérieure d'une première fenêtre de recherche pour identifier la première composante ; en identifiant la première composante sur la base de la première fenêtre de recherche ; en modifiant la fréquence centrale et la marge d'erreur de la composante candidate sur la base d'une fréquence centrale et d'une marge d'erreur de la première composante ; et en déterminant, sur la base de la fréquence centrale modifiée et de la marge d'erreur modifiée de la composante candidate, des limites de fréquence inférieure et supérieure d'une deuxième fenêtre de recherche pour identifier la deuxième composante. Selon encore un autre aspect, on prévoit un système comprenant : un système mécanique ou électromécanique ; au moins un capteur agencé pour détecter lesdites vibrations et/ou fluctuations électriques dans le système ; et le dispositif de signal provenant calcul susmentionné agencé pour traiter un 25 dudit au moins un capteur. Brève description des dessins Les fonctionnalités et avantages susmentionnés, et d'autres, apparaîtront clairement avec la description détaillée suivante de modes de réalisation, faite à titre d'illustration 30 et non de limitation, en référence aux dessins joints dans lesquels : la figure 1 illustre un système mécanique ou électromécanique muni d'une détection d'anomalies selon un exemple de réalisation ; B12549 - B232 - DI 6180-01 6 de spectre de de la figure 1 la figure 2 illustre un exemple fréquence d'un signal capturé par un capteur selon un exemple de réalisation ; représentant des la figure 3 est un -diagramme spectral paramètres dans un procédé de détection d'anomalie selon un exemple de réalisation ; la figure 4 est un organigramme illustrant des étapes dans un procédé de détection d'une anomalie dans un système mécanique ou électromécanique selon un exemple de réalisation de 10 la présente description ; et la figure 5 illustre un dispositif pour mettre en oeuvre un procédé de détection d'anomalie selon un exemple de réalisation de la présente description. Description détaillée 15 La figure 1 illustre un système 100 à surveiller. Ce système est par exemple un système mécanique ou électromécanique, comme une éolienne ou une hydrolienne, un moteur à combustion ou à réaction, un moteur électrique, une machinerie industrielle, ou d'autres types de systèmes comprenant des 20 parties mobiles dans lesquelles on veut surveiller des anomalies. L'anomalie à détecter est par exemple l'usure graduelle ou la défaillance d'une ou plusieurs parties du système. L'anomalie pourrait également être une opération non optimale de 25 la totalité ou d'une partie du système. On utilise un capteur 102 pour détecter des vibrations et/ou des fluctuations électriques en un point du système 100. Par exemple, le capteur 102 est un capteur de vibrations, comme un capteur MENS (capteur micro-électromécanique), un accélé- 30 romètre ou similaire. Dans un tel cas, le capteur 102 est par exemple en contact physique avec la partie du système à surveiller, ou en contact physique avec un élément du système qui est mécaniquement relié à une ou plusieurs parties du système à surveiller. A titre de variante, le capteur pourrait 35 être adapté à détecter des fluctuations électriques dans le B12549 - B232 - DI 6180-01 7 câblage électrique du système, comme sur une borne d'entrée ou de sortie d'une partie à surveiller, par exemple un Moteur électrique ou similaire. Un tel capteur est par exemple un détecteur de courant ou de tension, soit câblé électriquement soit couplé de façon électromécanique à des chemins conducteurs dans le système. Ici encore, le capteur peut être positionné de façon à pouvoir détecter des fluctuations provenant de plusieurs parties du système à surveiller. Bien qu'un seul capteur 102 soit illustré en figure 1, 10 dans des variantes de réalisation, il peut y avoir plusieurs capteurs pour surveiller divers points dans le système 100. La sortie de chaque capteur est-par exemple traitée séparément de la même manière que celle qui va être décrite plus en détail ci-après. 15 Le capteur 102, et tous les autres capteurs présents dans le système, fournissent un signal de sortie sur une ligne 104 à un dispositif de calcul 106. Le dispositif de calcul 106 réalise une analyse de signal sur le signal de sortie. En particulier, le dispositif de calcul 106 convertit le signal de 20 sortie variable dans le temps en un spectre de fréquence, par exemple en appliquant une transformation de Fourier discrète (DFT). Le dispositif de calcul 106 détecte ensuite des pics dans le spectre de fréquence- pour fournir une liste de composantes spectrales. Le dispositif de calcul 106 analyse ensuite ces 25 composantes spectrales pour identifier des séries de compo- santes, chaque série comprenant une composante de fréquence fondamentale, un ou plusieurs harmoniques de la fréquence fondamentale, et/ou une ou plusieurs composantes formant des bandes latérales de modulation des composantes fondamentales ou 30 harmoniques. Le dispositif de calcul 106 détecte par exemple la présence d'une anomalie sur la base de l'une quelconque de telles composantes spectrales identifiées. Dans certains modes de réalisation, le dispositif de calcul 106 fournit une indication de série identifiée de compo35 santes spectrales, et/ou des détails d'une anomalie détectée, B12549 - B232 - DI 6180-01 8 un module de sortie 108, qui est par exemple un afficheur, une interface avec un site web, ou un autre type de moyen de sortie informant ou alertant un expert en maintenance du système des résultats de l'analyse. Dans certains modes de réalisation, le dispositif de calcul 106 fournit en plus ou à la place un signal de sortie 110 renvoyé vers le système 100, qui par exemple altère le fonctionnement du système sur la base de l'anomalie détectée. Par exemple, le signal de sortie 110 peut être agencé pour immobiliser la totalité ou une partie du système dans le cas où une anomalie sérieuse est détectée, comme une pièce qui est sur le point de présenter une défaillance. En variante, le signal 110 peut être utilisé pour altérer le fonctionnement d'une ou plusieurs parties du système, comme réduire la vitesse de fonc- tionnement, jusqu'à ce qu'une réparation ou un remplacement de la pièce défectueuse puisse être réalisé. La figure 2 est un graphique illustrant un exemple de composantes spectrales générées par le dispositif de calcul 106 de la figure 1 sur la base du signal de sortie du capteur 102. Chacun des pics représentés en figure 2, qui dans cet exemple sont dans une plage allant d'une fréquence d'environ 35 Hz à une fréquence d'environ 105 Hz, correspond à une composantes spectrale. Bien sûr, les composantes de fréquence peuvent aussi être présentes à des fréquences plus basses ou plus élevées non représentées en figure 2. Chaque composante spectrale correspond à un pic dans le spectre de fréquence du signal de sortie du capteur 102, et est par exemple associée à une fréquence centrale correspondant à la fréquence au centre du pic, et à une marge d'erreur indiquant la précision de la valeur de la fréquence centrale. Par exemple, les composantes spectrales sont calculées sur la base du procédé décrit dans la publication intitulée "Vers une carte d'identité spectrale", Nadine Martin et al., GRETSI 2005, et/ou dans la demande PCT W02006/051242. En effet, cette demande de brevet, qui a des inventeurs en commun avec la B12549 - B232 - DI 6180-01 9 présente demande, concerne un procédé pour l'identification automatique de composantes spectrales d'un signal, et fournit, pour chaque pic d'un spectre de fréquence, une forme de carte d'identité indiquant des attributs de la composante spectrale correspondante, comme la fréquence centrale et l'erreur. Des techniques permettant de convertir le signal de sortie d'un capteur en un spectre de fréquence sur la base d'une transformation de Fourier discrète, et de détecter dans le spectre de fréquence 10 fréquence centrale capacités de l'homme cations mentionnées les composantes spectrales associées à une et une marge d'erreur, sont dans les de l'art, au moins sur la base des publici-avant, et ne seront pas décrites en détail ici. Par exemple, il sera clair pour l'homme de l'art que la marge d'erreur pour une composante de fréquence donnée peut 15 être déterminée en fonction de la fréquence de cette composante, et en fonction du nombre de coefficients utilisés lors de l'application de la transformation de Fourier discrète pour générer le spectre de fréquence. Est considéré comme une composante spectrale un pic 20 correspondant dans le spectre de fréquence qui dépasse par exemple un niveau de seuil minimum. La génération de nouveaux ensembles de composantes de fréquences est. par exemple réalisée périodiquement, par exemple avec une périodicité de quelques minutes ou de quelques heures. De cette façon, dès qu'un pic 25 dans le spectre de fréquence dépasse le niveau de seuil minimum, une composante spectrale correspondante va être générée, celle-ci pouvant indiquer la présence d'une nouvelle anomalie dans le système. 30 2, Parmi les composantes spectrales illustrées en figure certaines de ces composantes peuvent correspondre à des fréquences fondamentales présentes dans le système, alors que d'autres peuvent correspondre à des premiers, deuxièmes, troisièmes harmoniques ou des harmoniques supérieurs de ces B12549 - B232 - DI 6180-01 10 fréquences fondamentales. En outre, autour de toute fréquence fondamentale ou harmonique, des bandes latérales de modulation peuvent être présentes, aussi représentées par des composantes spectrales correspondantes. Par exemple, la composante située 5 autour de 50 Hz en figure 2 peut correspondre à une fréquence fondamentale, et une composante située juste au-dessus de 100 Hz peut correspondre à une première harmonique de la composante à 50 Hz. Une ou plusieurs des composantes de chaque côté des composantes à 50 et 100 Hz peuvent correspondre à des bandes 10 latérales de modulation. Une série de composantes spectrales, et en particulier la fréquence et l'amplitude d'une composante fondamentale et la présence et l'amplitude d'harmoniques et/ou de bandes de modulation associés, sont des indicateurs clés de la présence de 15 certaines anomalies au niveau de certaines parties du système mécanique ou électromécanique, et peuvent être analysés pour déterminer non seulement le type et l'étendue d'une anomalie, mais aussi la partie ou la pièce du système au niveau de laquelle l'anomalie prend sa source. Il y a une difficulté 20 technique pour identifier des séries de composantes spectrales. En effet, dans la sortie du capteur d'un système électromécanique typique, il peut y avoir des centaines ou des milliers de composantes spectrales identifiées, dont au moins certaines peuvent être relativement proches les unes des autres. Une 25 erreur dans la sélection d'une telle composante spectrale peut conduire au fait qu'une anomalie soit manquée ou mal interprétée. La figure 3 est un schéma représentant une série de composantes spectrales identifiées par un procédé selon un 30 exemple de réalisation. Comme cela est illustré, une composante spectrale Ci est sélectionnée parmi l'ensemble de composantes spectrales en tant que fréquence fondamentale candidate. La composante Ci est associée à une fréquence centrale Vi et à une marge d'erreur E. La marge d'erreur Ei correspond à la distance B12549 - B232 - DI 6180-01 11 positive ou négative par rapport à la valeur Vià laquelle la fréquence centrale Vi peut en réalité être située. L'existence d'une ou plusieurs composantes harmoniques associées à la composante Ci, et/ou d'une ou plusieurs compo- santes de bande de modulation associées à la composante Ci selon un rapport r est déterminée de la façon suivante. Le rapport r va dépendre de la composante harmonique ou de modulation à identifier. Par exemple, en supposant que le premier harmonique est à identifier, r est égal à 2. Une fenêtre de recherche est définie pour servir à la localisation d'une composante harmonique ou de modulation potentielle, ayant une limite inférieure égale à r(Vi-E) et une limite de fréquence supérieure égale à r(Vi+Ei). Ensuite, toute composante spectrale qui a une 15 fréquence centrale tombant dans la fenêtre de recherche après avoir pris en compte sa marge d'erreur est considérée comme étant une composante potentiellement associée. Par exemple, en figure 3, une composante de fréquence Cj a une fréquence centrale Vj tombant dans la plage de fréquence [r(Vi-E), rVi] 20 et une composante de fréquence Cj+i a une fréquence centrale V1+1, qui, lorsqu'elle est diminuée de la marge d'erreur Ej+i de la composante Cj+l, tombe dans la plage de fréquence [rVi, r(Vi+Ei)]. La condition peut être formalisée de la façon suivante pour la composante Cj, qui a une marge d'erreur Ej, où 25 [a, b] est la plage de fréquence de chevauchement : [a; b] = [v; - Ei; .17 Eii n [r(Vi - Ei); r(Vi + Es)] # 0 -1 Ensuite, l'une des composantes potentiellement associée est sélectionnée comme étant la composante la plus vraisemblablement associée à la composante candidate C. Par exemple, 30 la composante la plus proche du centre de la fenêtre de recherche, en d'autres termes la plus proche de rVi, est sélectionnée. Dans l'exemple de la figure 3, cela correspond à la composante C. Avant de répéter l'opération décrite précédemment pour 35 une nouvelle valeur de r afin d'identifier un autre élément de B12549 - B232 - DI 6180-01 12 la série comprenant la composante candidate Ci, la fréquence centrale Vi et la marge d'erreur Ei de la composante Ci sont misesn à jour sur la base des paramètres de la dernière composante identifiée, en d'autres termes la composante Cj dans l'exemple de la figure 3. En particulier, la composante Cj est supposée tomber dans la plage de fréquences de chevauchement [a, b] de la fenêtre de recherche avec la fréquence centrale Vj de la composante Cj plus et moins la marge d'erreur E. Ainsi, la limite inférieure "a" de cette plage est la fréquence r(Vi-E), et la limite supérieure "b" de cette plage en figure 3 est la fréquence Vj+Ej. La nouvelle fréquence centrale et la nouvelle marge d'erreur de la composante Ci deviennent par conséquent : = (a + b) (5, - a) Và = -2 2 x r ; 2 x r Ainsi la nouvelle fréquence centrale Vi est basée sur 15 la composante Cj se trouvant à mi-chemin entre a et b, et la nouvelle marge d'erreur est basée sur la moitié de l'intervalle de fréquence entre a et b. Le rapport r d'une composante harmonique ou de modulation suivante à identifier est par exemple une valeur 20 supérieure, par exemple égale à 3 dans le cas où la deuxième harmonique doit être identifiée. Par conséquent, plutôt qu'une nouvelle fenêtre de recherche soit 2r fois plus grande que la marge d'erreur initiale Ei de la composante Ci, la nouvelle fenêtre va avoir de 2r fois la taille de la valeur modifiée de 25 la marge d'erreur, ce qui est notablement plus petit. Une deuxième composante dans la nouvelle fenêtre de recherche peut être identifiée de la même façon que cela a été décrit précédemment, avant que le processus soit répété pour l'identification d'autant de composantes associées addition- 30 nelles qu'on le souhaite. La figure 4 est un organigramme illustrant un exemple de procédé pour détecter une anomalie dans un système mécanique ou électromécanique selon un exemple de réalisation. B12549 - B232 - DI 6180-01 13 On suppose en figure 4 qu'un ensemble de I composantes spectrales basé sur un signal détecté par un capteur a déjà été généré, et que chacune des composantes spectrales Ci est associée à une fréquence centrale Vi et à une marge d'erreur Ei, pour i compris entre 1 et I. Dans une opération 401, des paramètres i et k sont initialisés, par exemple à la valeur 1. Comme cela a été indiqué précédemment, le paramètre i désigne l'une des composantes Ci qui est appelée ici composante candidate. Le paramètre k désigne un rapport rk correspondant à une composante harmonique ou de bande de modulation courante à identifier, pour k égal à 1 à K. Le nombre total K de rapports va dépendre du nombre de composantes harmoniques/de modulation à identifier, et pourrait par exemple être un nombre quelconque entre 2 et 50. Il sera clair pour l'homme de l'art que ce nombre va dépendre de la taille du spectre, et de la fréquence fondamentale. Par valeurs r1 à rk sont respectivement égales à 2, Dans une opération suivante 402, candidate Ci est sélectionnée, 20 composante Cl. Dans une opération suivante 403, une fenêtre de correspondant initialement à 3, la exemple, les 4, ... K+1. fondamentale la recherche est déterminée pour la composante associée au rapport rk. Initialement, le rapport est rl, qui correspond par exemple à un rapport égal à- 2. La fenêtre de recherche est égale à 25 [rk(Vi-E), rk(Vi+Ei)]. Dans le cas d'une recherche de composantes spectrales de bande de modulation, on peut utiliser les mêmes valeurs de.r que celles utilisées pour les composantes harmoniques. Toutefois, en supposant que des bandes de modulation situées à 30 droite, c'est-à-dire ayant une fréquence supérieure, d'une composante Cm ayant une fréquence centrale Vm doivent être identifiées, la valeur Vm est par exemple soustraite de la fréquence centrale Vi de la composante candidate Ci avant de déterminer la fenêtre de recherche, et les limites inférieure et 35 supérieure de la fenêtre de recherche vont être égales à 2549 - 3232 - DI 6180-01 14 [Vm+rk(Vi-Ei), Vm+rk(Vi+Ei)]. Si les bandes de modulation situées à gauche, c'est-à-dire de fréquence inférieure, de la composante Pm doivent être identifiées, la fréquence centrale Vi de la composante candidate est par exemple soustraite de la valeur Vm avant de déterminer la fenêtre de recherche, et les limites inférieure et supérieure de la fenêtre de recherche vont être égales à [Vm-rk(Vi+Ei), Vm-rk(Vi-Ei)]. Dans une opération suivante 404, on détermine si une ou plusieurs composantes ont été identifiées dans la fenêtre de recherche, en prenant en compte la marge d'erreur Ej de chaque composante C. En particulier, la condition de la formule 1 précédente est appliquée pour les composantes C. Si au moins une composante est trouvée dans la fenêtre de recherche, l'opération suivante est l'opération 405, dans laquelle une composante est sélectionnée parmi les composantes trouvées. Comme cela a été mentionné précédemment, la composante la plus proche du niveau de fréquence rkVi est par exemple sélectionnée. L'opération après l'opération 405, ou après l'opération 404 si aucune composante n'est trouvée dans la fenêtre de recherche, est l'opération 406, dans laquelle on vérifie si k est égal à K. Dans la négative, l'opération suivante est l'opération 407, dans laquelle k est incrémenté, puis dans une opération suivante 408, la valeur Vi de la fréquence centrale et l'erreur Ei de la composante Ci sont redéfinies sur la base de la composante identifiée dans l'opération 404. En particulier, les formules 2 décrites précédemment sont par exemple appliquées. Le procédé revient ensuite à l'opération 403 pour déterminer une nouvelle fenêtre de recherche sur la base de la nouvelle valeur de rk et des valeurs modifiées de Vi et E. Si dans l'opération 407, on a trouvé que k est égal à K, l'opération suivante est l'opération 409, dans laquelle la série Hi des composantes harmoniques et/ou de modulation associées à la composante Ci est définie. Une telle famille Hi est par exemple définie seulement si au moins une composante B12549 - 3232 - DI 6180-01 15 harmonique ou de modulation associée à la composante Ci a été identifiée. Après l'opération 409, l'opération suivante est l'opération 410, dans laquelle on vérifie si i est égal à I. 5 Dans la négative, i est incrémenté dans une opération 411, et le procédé revient à l'opération 402, où une nouvelle composante candidate est sélectionnée sur la base de la nouvelle valeur de i. Dans l'autre cas, si dans l'opération 410 on a trouvé que i est égal à I, le processus par exemple se termine au niveau 10 d'une opération 412, dans laquelle une anomalie est identifiée sur la base de la série Hi de composantes qui ont été identifiées. Ainsi le procédé continue jusqu'à ce que toutes les composantes Ci à CI aient été évaluées, pour tous les rapports 15 r1 à rK ou jusqu'à ce qu'aucune composante ne puisse être identifiée dans la fenêtre de recherche. On notera que certains éléments de la série Hi peuvent comprendre la même composante, si par exemple une telle composante résulte de la somme d'harmoniques de fréquences fondamentales différentes. 20 L'identification d'une anomalie réalisée dans l'opé- ration 412 est par exemple obtenue en déterminant au moins l'une des caractérisations suivantes d'une série identifiée. = card(Hi) Di -3 max max r. I vmax 1 -4 R» = Nnwith Nrx [v J / riax 25 THDi = -5 -2 A0 où Di est la densité de la série, rimax est le rapport rk le plus grand dans la série Hi, card(Hi) est la cardinalité (taille) de la série Hi, Ri est la richesse de la série, Nimax est la taille maximum de la série sur la base de la fréquence 30 fondamentale Vi de la série et de la fréquence maximum Vmax dans le spectre, Lxi est l'entier le plus grand qui n'est pas B12549 - B232 - DI 6180-01 16 supérieur à x, THDi est le taux de'distorsion dans la série Hi, sur la base des amplitudes Al à AN des composantes 1 à N dans la série. La figure 5 illustre un exemple d'un dispositif de 5 traitement 500 du dispositif de calcul 106 de la figure 1. Le dispositif de traitement 500 est agencé pour mettre en oeuvre le procédé de détection d'anomalie décrit ici. Le dispositif 500 comprend par exemple un ou plusieurs processeurs 502 (P) sous le contrôle d'instructions mémorisées dans une mémoire d'ins- 10 tructions 504. En outre, le dispositif 500 comprend par exemple une mémoire 506 par exemple configurée pour mémoriser la liste de composantes spectrales Ci à CI, accompagnée d'au moins les fréquences centrales Vi correspondantes et de la marge d'erreur Ei pour chaque composante, comme cela a été décrit précédemment, 15 et de toutes les séries Hi identifiées. En outre, la mémoire 506, peut mémoriser le signal variable dans le temps provenant du capteur 102, qui est par exemple reçu par le dispositif de traitement 502 par l'intermédiaife d'une interface d'entrée/sortie 508 couplée à la ligne d'entrée 104. L'interface d'entrée/sortie 20 508 peut comprendre un convertisseur analogique-numérique pour échantillonner le signal d'entrée et le numériser. L'interface d'entrée/sortie 508 peut aussi être couplée à la ligne de sortie 110 pour contrôler le fonctionnement du système mécanique OU électromécanique dans le cas où une anomalie est détectée. 25 Le dispositif 500 comprend par exemple en plus un afficheur 510, qui par exemple fournit une interface d'utilisateur, ainsi que des moyens pour alerter un utilisateur en indiquant si et quand une anomalie a été détectée. Un avantage des modes de réalisation décrits ici est 30 qu'une anomalie dans un système peut être détectée avec précision même dans des systèmes relativement complexes dans lesquels le spectre de fréquence calculé par un capteur comprend de nombreuses composantes spectrales. B12549 - B232 - DI 6180-01 17 Avec la description ainsi faite d'au moins un mode de réalisation illustratif, diverses altérations, modifications et améliorations apparaîtront facilement à l'homme de l'art. Par exemple, il sera clair pour l'homme de l'art que, bien qu'on ait décrit un procédé Utilisant des formules spécifiques, il s'agit simplement d'exemples d'une façon de mettre en oeuvre le procédé, et que diverses modifications pourraient être appliquées à ces formules. REVENDICATIONS1. Procédé pour détecter une anomalie dans un système mécanique ou électromécanique, comprenant les étapes suivantes : identifier par un dispositif de traitement, à partir d'un ensemble de composantes spectrales généré sur la base d'un signal provenant d'au moins un capteur détectant des vibrations et/ou des fluctuations électriques dans le système, au moins des première et deuxième composantes spectrales, chaque composante représentant un harmonique ou une modulation associé à une composante candidate (Ci) : en déterminant, sur la base d'une fréquence cen- trale (Vi) et d'une marge d'erreur (Et) de la composante can- didate, des limites de fréquence inférieure et supérieure (rl(Vi±Ei)) d'une première fenêtre de recherche pour identifier la première composante ; en identifiant la première composante sur la base de la première fenêtre de recherche ; en modifiant la fréquence centrale et la marge d'erreur de la composante candidate sur la base d'une fréquence centrale et d'une marge d'erreur de la première composante ; et en déterminant, sur la base de la fréquence centrale modifiée (Vi) et de la marge d'erreur modifiée (Et) de la composante candidate, des limites de fréquence inférieure et supérieure (r2(V1±Ei)) d'une deuxième fenêtre de recherche pour identifier la deuxième composante. 2. Procédé selon la revendication 1, comprenant en outre l'identification de la deuxième composante dans la deuxième fenêtre de recherche. 3. Procédé selon la revendication 1 ou 2, comprenant en outre la définition d'une série de composantes spectrales (Hi) comprenant au moins la composante candidate et les première et deuxième composantes, et la détection de l'anomalie au moins partiellement sur la base de ladite série. 4. Procédé selon l'une quelconque des revendications 1 à 3, comprenant en outre la génération d'un signal de sortieB12549 - B232 - DI 6180-01 19 pour contrôler une partie du système mécanique ou électromécanique sur la base de l'anomalie détectée. 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel les première et deuxième composantes sont des harmoniques de la composante candidate; et dans lequel l'anomalie est détectée sur la base de paramètres d'une ou plusieurs composantes de modulation associées à la composante candidate ou à la première ou la deuxième composante. 6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel chacune des limites inférieure et supérieure des première et deuxième fenêtres de recherche est déterminée sur la base des équations rk(VI-E) et rk(Vi+Ei), où rk est l'ordre harmonique des première et deuxième composantes respectivement, Vi est la fréquence centrale associée à la composante candidate, et Ei est la marge d'erreur associée à la composante candidate. 7. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel la première composante est considérée comme tombant dans la première fenêtre de recherche si la condition 20 suivante est vérifiée : [a; b] = lv - j; V + Eii n [ric(vi Ei); rk(vi + El)] o où [a ; b] est une plage de fréquence de chevauchement, Vj est la fréquence centrale de la première composante, Ej est la marge d'erreur de la première composante, Vi est la fréquence centrale 25 de la composante candidate, Ei est la marge d'erreur de la composante candidate et rk est l'ordre de l'harmonique. 8. Procédé selon la revendication 7, dans lequel à la suite de l'identification de la première composante, la fréquence centrale Vi et la marge d'erreur Ei sont mises à jour 30 de la façon suivante : (a+= t) (be - a) v ; 2 x r 2 x rB12549 - B232 - DI 6180-01 20 où a est la limite inférieure de la plage de fréquence de chevauchement et b est la limite supérieure de la plage de fréquence de chevauchement. 9. Procédé selon l'une quelconque des revendications 1 à 8, comprenant en outre la détermination d'au moins un élément parmi une valeur de densité Di de la série, une richesse Ri de la série et un taux de distorsion THDi de la série sur la base des équations suivantes : carel(111) D. = max r. 1 max r. with eax - max Ni THDi = Â02 où Di est la densité de la série, rimax est le plus grand rapport rk dans la série Hi, card(Hi) est la cardinalité (taille) de la série Hi, Ri est la richesse de la série, Nimax 15 est la taille maximum de la série sur la base de la fréquence fondamentale Vi de la série et de la fréquence maximum Vmax dans le spectre, [x] est l'entier le plus grand qui n'est pas supérieur à x, THDi est le taux de distorsion dans la série Hi, sur la base des amplitudes Ai à AN des composantes 1 à N dans la 20 série. 10. 10. Support lisible par un ordinateur mémorisant un programme d'ordinateur qui, lorsqu'il est exécuté par un dispositif de traitement, amène la mise en oeuvre du procédé de l'une quelconque des revendications 1 à 9. 25 11. Dispositif de calcul pour détecter une anomalie dans un système mécanique ou électromécanique comprenant un dispositif de traitement agencé pour : identifier, à partir d'un ensemble de composantes spectrales généré sur la base d'un signal provenant d'au moins 30 un capteur détectant des vibrations et/ou des fluctuations électriques dans le système, au moins des première et deuxième10B12549 - 5232 - DI 6180-01 21 composantes spectrales, chaque composante représentant un harmonique ou une modulation associé à une composante candidate (Ei) en déterminant, sur la base d'une fréquence 5 centrale (Vi) et d'une marge d'erreur (Ei) de la composante candidate, des limites de fréquence inférieure et supérieure (ri(Vi±Ei)) d'une première fenêtre de recherche pour identifier la première composante ; en identifiant la première composante sur la base 10 de la première fenêtre de recherche ; en modifiant la fréquence centrale et la marge d'erreur de la composante candidate sur la base d'une fréquence centrale et d'une marge d'erreur de la première composante ; et en déterminant, sur la base de la fréquence 15 centrale modifiée (Vi) et de la marge d'erreur modifiée (Et) de la composante candidate, des limites de fréquence inférieure et supérieure (r2(V1±Ei)) d'une deuxième fenêtre de recherche pour identifier la deuxième composante. 12. Système comprenant : 20 un système mécanique ou électromécanique ; au moins un capteur agencé pour détecter lesdites vibrations et/ou fluctuations électriques dans le système ; et le dispositif de calcul de la revendication 11 agencé pour traiter un signal provenant dudit au moins un capteur.