Détection automatique des défauts de roulements
Détection automatique des défauts de roulements
Université de Toulouse ; LAPLACE ; CNRS, INPT, UPS ; 2 rue C. Camichel, BP 7122, 31071 Toulouse Cedex 7
email : [Link]@[Link]
Résumé—Cet article traite de la détection de défauts de électriques telles que les courants de phase du moteur.
roulements dans les entraînements à base de machine Plusieurs études ont démontré que des signatures relatives
asynchrone par analyse des courants de phase aux défauts de roulements apparaissaient dans les spectres de
statoriques. D’une part, on démontre que les défauts de courant [9] – [10]. Pour cela, deux approches sont
roulements entraînent des oscillations du couple de envisagées : l’une dans laquelle les défauts de roulements
charge à des fréquences caractéristiques du roulement. induisent de l’excentricité [9], l’autre dans laquelle les
D’autre part, un modèle de courant statorique démontre mêmes défauts induisent des oscillations du couple de charge
qu’en présence d’oscillations du couple de charge, il [10]. Les signatures spectrales du courant sont alors
existe des modulations de phase sur ces courants. Ces explorées dans [11]. Cependant, peu d’articles concernent la
modulations de phase se manifestent dans le domaine mise en place d’une détection automatisée à partir de ces
fréquentiel par l’apparition de composantes latérales sur signatures spectrales sur les courants statoriques [12]. Ce
le spectre du courant qui sont utilisées pour la détection travail présente la définition et la validation d’un indicateur
des défauts de roulements. Pour mettre en place un automatique pour la détection de défauts de roulements par
schéma de détection efficace, il est impératif d’étudier les extraction d’énergie dans le spectre des courants statoriques
lois de variations des amplitudes de ces composantes en de la machine asynchrone.
fonction de la fréquence de modulation considérée. La Dans cet article, un récapitulatif des effets mécaniques
fonction de transfert électromécanique de l’entrainement des défauts de roulements est donné dans le paragraphe II. Le
qui relie l’amplitude de ces composantes latérales à la paragraphe III rappelle quant à lui la formulation d’un
fréquence des oscillations de couple, est ainsi étudiée. Elle modèle de courant statorique en cas d’oscillations de couple.
fait apparaître l’existence d’une résonance permettant Par ailleurs, pour plusieurs points de fonctionnement de
d’amplifier les effets des oscillations de couple sur les
l’entraînement, la réponse fréquentielle des harmoniques de
courants statoriques et permet ainsi d’améliorer la
courant relatifs aux oscillations du couple de charge est
qualité de la détection des défauts. Un indicateur de
étudiée en simulation et en expérimentation sur le banc de
détection de défauts de roulements est donc présenté puis
testé sur des roulements artificiellement dégradés. test. La réponse fréquentielle fait apparaître une résonance
Finalement, l’indicateur est validé sur des roulements qui est ensuite utilisée pour favoriser la détection des défauts
usés provenant du service après-vente d’un fabricant de de roulements. Dans le paragraphe IV, un indicateur de
machines. défauts de roulements est présenté, intégrant une extraction
énergétique localisée dans des bandes fréquentielles
I. INTRODUCTION spécifiques aux défauts étudiés. Celui-ci est validé sur des
roulements avec des défauts localisés puis sur un roulement
Les entraînements électriques à base de machine
naturellement dégradé.
asynchrone sont largement utilisés dans les applications
industrielles en raison de leur faible coût, de leurs
II. EFFETS MECANIQUES DES DEFAUTS DE ROULEMENTS
performances et de leur robustesse. Cependant, des modes de
fonctionnements dégradés peuvent apparaître durant la vie de Les défauts de roulements entraînent de manière générale
la machine. Les références [1] – [2] passent en revue un plusieurs effets mécaniques dans les machines tels qu’une
grand nombre de cas de défaillances, ainsi que des techniques augmentation du niveau sonore [5] et l’apparition de
de détection associées. La raison principale des défaillances vibrations. Nous montrerons également que les défauts de
de machine reste les défauts de roulements à billes [3] – [4] roulements induisent des variations dans le couple de charge
qui ont de nombreuses causes telles que la contamination du de la machine asynchrone.
lubrifiant, une charge excessive [5] ou encore des causes
électriques comme la circulation de courants de fuite induits A. Fréquences caractéristiques des roulements à billes
par les onduleurs à MLI [6] – [7]. Des fréquences caractéristiques, provenant de calculs
Dans les cas de défaillances de roulements à billes, une cinématiques, sont associées avec les différents organes
surveillance vibratoire peut être mise en place [8] mais celle- constitutifs d’un roulement à billes : la bague externe, la
ci est souvent onéreuse du fait des capteurs piézoélectriques bague interne, la cage et les billes [3], [13] (1). En cas de
et de leurs conditionneurs. Pour y remédier, la détection et la défaut se présentant sur une partie du roulement, des
surveillance peuvent être basées sur l’analyse de grandeurs vibrations à la fréquence caractéristique associée sont
JCGE'08 LYON, 16 et 17 décembre 2008
engendrées. Des harmoniques supplémentaires, issus de de rotation est nominale, soit proche de 25Hz.
-90
combinaisons entre les fréquences caractéristiques peuvent Defaut bague externe
[Link]-fc [Link]
apparaître dans les spectres vibratoires (cf. Fig. 1). -100 Cas sain
fr D cos(θ )
f be = N B 1 − B
-110
2 D
f D cos(θ ) -130
f bi = r N B 1 + B (1)
2 D -140
p -150
f D cos(θ )
f c = r 1 − B -160
2 Dp
-170
165 170 175 180 185
D cos(θ ) 2
Fréquence (Hz)
N B 1 − B
fr Fig. 2. Spectre du couple de charge – Comparaison cas sain et cas défaillant
f bille =
2 Dp On peut remarquer en Fig. 2, sur les spectres de couple en
présence d’un défaut de roulement, l’apparition de raies
avec : harmoniques à des fréquences caractéristiques des défauts
- fbe est la fréquence de bague externe, bague externe. En revanche, notons que toutes les
- fbi est la fréquence de bague interne, combinaisons possibles des fréquences caractéristiques
- fc est la fréquence de cage, n’existent pas. En conséquence, le contenu spectral du couple
- fbllei est la fréquence de bille, de charge relatif aux défauts de roulements ne peut être prévu
- fr est la fréquence de rotation mécanique,
de manière complète et théorique.
- NB est le nombre de billes du roulement,
Notons aussi que la détection des défaillances de
- DB est le diamètre des billes,
roulements pourrait se faire via l’analyse du contenu spectral
- Dp est le diamètre moyen du roulement,
du couple de charge. Cependant, l’implantation d’un
- θ est l’angle de contact.
-65 couplemètre est bien souvent difficile à réaliser et très
multiples de fbe
onéreuse.
-70
-80
A. Modèle de courant statorique
-85
De nombreux travaux [4], [11], [14] montrent que les
oscillations du couple de charge provoquent
-90
préférentiellement des modulations de phase (PM) sur les
modulations à fc modulations à fr courants statoriques. Comme le montre la Fig. 2, le couple de
-95
10 20 30 40 50 60 70 80
charge peut être représenté par une composante continue
Fréquence (Hz) d’amplitude C0 égale au couple électromagnétique de la
Fig. 1. Spectre vibratoire en cas de défaut bague externe machine Cem et par une série de n harmoniques d’amplitude
Dans le cas des roulements à billes de type 6208 équipant Cn et de pulsation ωn relatifs aux défauts de roulements (2).
notre machine asynchrone, pour une fréquence statorique de Cch arg e (t ) = C0 + ∑ C n cos(ω n t ) (2)
50Hz et donc une fréquence de rotation mécanique proche de n
25Hz, les fréquences de bague externe, interne et de cage Physiquement, en considérant connue la fonction de
valent respectivement fbe=89Hz, fbi=136Hz et fc=9.9Hz. transfert mécanique F(jωn), l’expression de la position
angulaire du rotor θ(t) (3) de la machine asynchrone
B. Oscillations du couple de charge dues aux défauts de comporte également des oscillations aux mêmes fréquences.
roulements
Cn
Lors de la rotation d’un roulement présentant un défaut θ (t ) = θ 0 + Ω 0 t + ∑ cos(ω n t + ϕ n ) (3)
mécanique, il se produit des chocs entraînant des vibrations n ω n F ( jω n )
mais également des à-coups de couple. La répétition La position mécanique permet de calculer la force
périodique de ces à-coups de couple dans le temps induit des magnétomotrice (FMM) rotorique qui se trouve modulée en
oscillations du couple de charge qui sont mesurées sur le phase [11]. Par ailleurs, la FMM statorique est purement
banc de test à l’aide d’un couplemètre. La Fig. 2 montre, sur sinusoïdale. La multiplication de ces FMM par la perméance
une plage fréquentielle réduite du spectre du couple de d’entrefer, considérée comme constante selon la position
charge, une comparaison du contenu harmonique pour une angulaire et le temps, conduit à l’expression du champ
machine saine et pour une machine équipée d’un roulement magnétique. L’intégrale de ce champ sur la surface d’entrefer
présentant un défaut localisé sur la bague externe. La vitesse mène au flux magnétique et ainsi à la force électromotrice
JCGE'08 LYON, 16 et 17 décembre 2008
(FEM) aux bornes d’une phase par dérivation temporelle. En schéma de la Fig. 3 par des inerties J, une raideur
considérant une relation linéaire liant la FEM au courant de d’accouplement kc et des frottements f. La paramétrisation de
phase, on obtient une expression du courant de phase (4) ce modèle est réalisée à partir de mesures pour les inerties et
avec Is l’amplitude du courant statorique, Ir l’amplitude du frottements, la raideur kc est quant à elle ajustée de manière
courant rotorique, ωs la pulsation fondamentale des courants empirique par rapport aux données constructeurs pour être en
statoriques et p le nombre de paires de pôles de la machine. accord avec les mesures. La modélisation du système
i (t ) = I s cos(ωst + φs ) mécanique conduit à une fonction de transfert d’ordre 2 au
numérateur et 3 au dénominateur.
+ I r sin ω s t + ∑ β n cos(ω n t + ϕ n ) (4)
Il est alors possible de comparer les résultats de simulation
n et de mesure en Fig. 4 pour une fréquence d’alimentation de
la MAS de 50Hz.
Cn
avec β n = p .
ω n F ( jω n ) Ωmas, Cem-Ccharge Ωmas, Caxe Ωmcc, Caxe
Jmas Jmcc
Le premier terme de (4) est lié à la contribution de la
kc
FMM statorique et le second terme, modulé en phase, est lié
fmas fmcc
à la contribution de la FMM rotorique. En considérant des
Fig. 3. Modèle du système mécanique
modulations de phase de faible amplitude, on obtient une
expression simplifiée du module de la Transformée de 5
Diagramme mesuré
Fourier (FT) du courant (5) en fonction de la fréquence f avec 0 Diagramme simulé
FT {i (t )} = (I s + I r )δ ( f − f s ) -10
βn
-15
+ Ir ∑ δ ( f − ( f s ± f n ))
Gain (dB)
(5) -20
n 2 -25
de modulation, c'est-à-dire obtenir le ratio βosc/Cosc en Fig. 5. Diagrammes de gain mesurés pour plusieurs fréquences
fonction de fosc. d’alimentation
Dans un second temps, afin de caractériser de manière On constate, en simulation comme par l’expérience, que le
théorique la loi recherchée, une approche de modélisation est comportement fréquentiel du système varie en fonction de la
proposée. La machine asynchrone est modélisée par ses fréquence d’alimentation statorique. La Fig. 5 illustre ce
équations d’état avec comme variables d’état les courants point en montrant, pour différentes fréquences
statoriques et les flux rotoriques dans le repère de Concordia d’alimentation, l’évolution du diagramme de gain.
(α,β) [15]. La charge mécanique est modélisée selon le On peut remarquer que l’amplitude de la résonance décroît
JCGE'08 LYON, 16 et 17 décembre 2008
avec la fréquence statorique mais que sa fréquence est mesures, les transformées de Fourier (FT) de deux courants
constante. La présence de cette résonance permet d’envisager de phase sont multipliées afin de corréler les signatures [16].
son utilisation pour amplifier les effets des oscillations de Cette méthode permet d’augmenter l’efficacité de la
couple sur les signatures fréquentielles des courants lors détection d’harmoniques.
d’une détection.
A. Définition de l’indicateur spectral
C. Exploitation du point de résonance L’étude du couple de charge en cas de défauts de
Il a été démontré que les défauts de roulements à billes roulements montre que les contenus spectraux du couple et
induisaient des oscillations du couple de charge. Ces donc du courant ne peuvent être connus de manière
oscillations sont en général de faible amplitude. Or, nous théorique. Par conséquent, il est nécessaire de rechercher
avons montré l’existence d’une résonance dans les dans les spectres de courant, non pas l’apparition
diagrammes précédents, qu’il pourrait être utile d’utiliser d’harmoniques particuliers, mais plutôt une augmentation
comme un amplificateur naturel des défauts de roulements. d’énergie dans des bandes spectrales où des harmoniques
Sachant que les fréquences caractéristiques des roulements relatifs aux défauts sont susceptibles d’apparaître.
sont linéairement dépendantes de la fréquence de rotation L’indicateur de défaut est donc défini comme une
mécanique et donc de la fréquence d’alimentation de la extraction d’énergie dans les bandes latérales de modulation
MAS, il faut identifier les fréquences statoriques de sorte que à fs±nfdef où fdef peut être la fréquence caractéristique de
les fréquences caractéristiques du roulement soient dans la bague externe ou de bague interne. De plus les bandes de
bande passante d’amplification du diagramme de gain. fréquences sont étendues de façon à tenir compte des
Ainsi, pour détecter via l’analyse des spectres de courant, éventuels harmoniques faisant intervenir la fréquence de
préférentiellement les effets des défauts localisés sur la bague rotation mécanique et la fréquence de cage (8).
externe du roulement (défauts entraînant l’apparition
d’harmoniques de couple de charge relatifs à la fréquence de
[
f s ± nf def − f c ; nf def + f c ]
bague externe) il est nécessaire d’identifier fs telle que fbe≈fres. fs ± [nf def − f r − f c ; nf def − f r + f c ] (8)
± [nf + f ]
Cette condition est réalisée pour une fréquence
d’alimentation de 13.3Hz soit une vitesse de rotation fs def + f r − f c ; nf def + f r c
n ∈ [1;5]
d’environ 400tr/min. Dans ce cas, la fréquence
avec
caractéristique de la bague interne est de fbi≈36Hz. Pour une
telle fréquence statorique, les effets sur les courants des Pour estimer la vitesse de rotation, les harmoniques
oscillations de couple dues à un défaut localisé sur la bague d’encoches rotoriques sont détectés sur les courants de phase
interne sont atténués d’environ 6dB par rapport aux effets [17]. Pour chaque intervalle fréquentiel, l’énergie spectrale
est intégrée et normalisée par la valeur maximum dans la
dus au défaut de la bague externe. Ceci permet bien de
bande. La Fig. 6 décrit le principe de l’extraction
focaliser la détection sur un type de défaut particulier.
énergétique. L’indicateur proposé utilise la différence
Réciproquement, pour détecter préférentiellement les effets
d’énergie relative entre un spectre de courant sain (référence
des défauts localisés sur la bague interne du roulement, il est de mesure) et les spectres en cas de défaut dans chaque
nécessaire d’identifier fs telle que fbi≈fres. Cette condition est intervalle de fréquence indicé par k. Les écarts relatifs pour
réalisée pour une fréquence statorique de 6.7Hz soit une les bandes de fréquences relatives aux fréquences de bague
vitesse de rotation d’environ 200tr/min. Dans ce cas, la externe et interne (∆Efbe et ∆Efbi respectivement) sont sommés
fréquence caractéristique de la bague interne vaut fbe≈12Hz. (9) pour obtenir un écart relatif global ∆Etot.
Les harmoniques sur le courant liés au défaut de la bague
∆Etot (k ) = ∆E fbe (k ) + ∆E fbi (k ) (9)
externe sont donc atténués d’environ 3dB par rapport à ceux
liés au défaut de la bague interne. Par ailleurs, plus la vitesse Enfin, une somme cumulée de ∆Etot est calculée pour
de rotation est basse plus les effets des défauts de roulements construire l’indicateur. Seule la dernière valeur de la somme
sur le couple de charge sont faibles. Ainsi, pour une cumulée définit l’indicateur de défaillance (cf. Fig. 7).
fréquence statorique de 6.7Hz, les défauts de roulements mis
préférentiellement en évidence grâce à la résonance, sont des
défauts localisés sur la bague interne.
Conformément à la définition de l’indicateur, celui-ci doit résonance est essentielle pour la détection des défauts de
être proche de 0, et donc proche de la référence énergétique roulements. Notons que l’utilisation de la résonance permet,
lorsqu’un roulement sain est testé et doit être largement en plus d’assurer la détection des défauts, d’effectuer un
positif lorsqu’un roulement défaillant est testé. diagnostic du défaut en précisant sa localisation.
250
Défaut bague externe
Défaut bague interne
200 Cas sain
150
50
-50
2 4 6 8 10 12 14
Fig. 7. Calcul de l’indicateur Termes des sommes cumulées
(a). fs=13.3Hz
B. Détection de défauts localisés artificiels
250
Tout d’abord, afin d’illustrer la pertinence de l’indicateur Défaut bague externe
pour des défauts de roulements, deux cas de défauts localisés Défaut bague interne
200 Cas sain
sont étudiés. Le premier concerne un défaut constitué d’un
trou de 3mm de diamètre dans la bande de roulement externe Erreur relative (%)
150
et le second d’un trou identique dans la bande de roulement
interne. Enfin un troisième roulement sain est utilisé à titre de 100