FR3140944A1 - Procédé de détermination de la fiabilité d’un roulement - Google Patents
Procédé de détermination de la fiabilité d’un roulement Download PDFInfo
- Publication number
- FR3140944A1 FR3140944A1 FR2210599A FR2210599A FR3140944A1 FR 3140944 A1 FR3140944 A1 FR 3140944A1 FR 2210599 A FR2210599 A FR 2210599A FR 2210599 A FR2210599 A FR 2210599A FR 3140944 A1 FR3140944 A1 FR 3140944A1
- Authority
- FR
- France
- Prior art keywords
- bearing
- values
- parameters
- lifespan
- rolling body
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000000463 material Substances 0.000 claims abstract description 15
- 238000005096 rolling process Methods 0.000 claims description 46
- 238000004364 calculation method Methods 0.000 claims description 26
- 238000004088 simulation Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 230000007547 defect Effects 0.000 description 11
- 238000004901 spalling Methods 0.000 description 4
- 238000004513 sizing Methods 0.000 description 3
- 208000031968 Cadaver Diseases 0.000 description 2
- 238000005336 cracking Methods 0.000 description 2
- 235000008612 Gnetum gnemon Nutrition 0.000 description 1
- 240000000018 Gnetum gnemon Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000009661 fatigue test Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012066 statistical methodology Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
L’invention concerne un procédé de détermination de la fiabilité d’un roulement consistant à déterminer les paramètres géométriques, de matériaux et de contrainte résiduelle ainsi que les paramètres moteurs pour un roulement donné puis à réaliser une pluralité d’association desdits paramètres précités. Dans un second temps, des modèles d’endommagement du roulement sont utilisés pour déterminer une durée de vie pour chacune desdites association de paramètres. Enfin, à partir desdites durées de vie, on établit un taux de défaillance que l’on compare à un taux de défaillance admissible. Figure de l’abrégé : Figure 2
Description
La présente invention concerne un procédé de détermination de la fiabilité d’un roulement
Dans le domaine aéronautique, on utilise des roulements, tels que des roulements à billes ou à rouleaux, pour permettre la rotation notamment des rotors du compresseur et de la turbine. Classiquement, un roulement est évalué avec selon un critère dit critère L10 développé par Lundberg et Palmberg [Lundberg G, Pamgre, A. Dynamic capacity of rolling bearings. J Appl Mech 1949 ; 16(2) :165-72] qui évalue le temps au bout duquel 10% des roulements testés ont subi une défaillance. Or, un tel taux de défaillance s’avère inacceptable dans le domaine aéronautique où le taux de défaillance doit être de l’ordre 10- 7, ce qui signifie une occurrence de défaillance pour dix millions d’heures de vol. Par ailleurs, le critère L10 impose un recalage par la réalisation d’essai, ce qui conduirait à devoir réaliser cent millions d’heures d’essais pour arriver à évaluer un taux de défaillance de l’ordre de 10-7.
Ainsi, actuellement, on a recours une détermination empirique des paramètres de dimensionnement des roulements et des matériaux de ceux-ci en fonction du moteur sur lequel ils sont destinés à être installés. On a donc recours à des abaques et à la connaissance technique d’opérateurs qui s’avèrent délicats à utiliser lorsque les paramètres dimensionnels et de matériaux des roulements et/ou les paramètres de fonctionnement du moteur sont modifiés en dehors des plages préétablis en base de données.
Il existe donc un besoin pour permettre la validation de l’adéquation des paramètres dimensionnels et de matériaux d’un roulement avec les paramètres de fonctionnement du moteur sur lequel il est destiné à être installé.
Il est ainsi proposé un procédé de test de la fiabilité d’un roulement pour un moteur, tel qu’un moteur de turbomachine, comprenant les étapes suivantes :
E1)
- Déterminer paramètres géométriques, de matériaux et de contrainte résiduelle de roulement nécessaires à la réalisation d’un calcul d’équilibre roulement, étant un indice entier variant entre et et désignant un desdits paramètres, et
- Obtenir valeurs pour chacun des paramètres , étant un indice variant de 1 à et désignant une valeur d’un paramètre ,
E2)
- Déterminer paramètres moteur nécessaires à la réalisation d’un calcul d’équilibre roulement, étant un indice entier variant entre 1 et et désignant un desdits R paramètres, et
- Obtenir valeurs pour chacun des paramètres , étant un indice variant de 1 à et désignant une valeur d’un paramètre ,
E3) Associer aléatoirement une série de valeurs à une série de valeurs de manière à obtenir N couples de séries de valeurs des paramètres / ,
E4) Pour chacun des couples de séries , déterminer une série de valeurs de paramètres de sortie d’un calcul d’équilibre du roulement,
E5) Etablir un modèle d’endommagement d’une sous-couche des composants unitaires du roulement, à savoir la bague externe, la bague interne et le corps roulant puis détermination de N valeurs minimales de durée de vie en sous couche,
E6) Etablir d’un modèle d’endommagement de surface des composants unitaires du roulement, à savoir la bague externe, la bague interne et le corps roulant puis détermination de N valeurs minimales de durée de vie en surface,
E7) Déterminer une valeur de durée de vie globale du minimum parmi et
E8) Comparer chacune des N valeurs de à la durée de vie cible
E9) Déterminer un taux de défaillance.
Dans notre méthode, nous utiliserons une approche de type data-driven (approche dirigée par les données), combinée avec la physique du roulement et une méthodologie Monte-Carlo afin d’évaluer de façon exacte la fiabilité des roulements de ligne arbre. Il s’agit ici de combiner une approche déterministe et une méthodologie statistique de type Monte-Carlo afin d’estimer de façon exacte la fiabilité de roulements sur une ligne d’arbre d’une turbomachine par exemple. On notera que cette approche pourrait être utilisée dans d’autres domaines que celui de l’aéronautique avec les bénéfices exposés ici.
L’étape E1 peut comprendre :
a. Détermination de paramètres de roulement nécessaires à la réalisation d’un calcul d’équilibre roulement, étant un indice entier variant entre et ,
b. Détermination d’une loi de variation de chacun des ,
c. Génération de valeurs pour chacun des paramètres en utilisant la loi de variation obtenue à l’étape précédente.
L’étape E2 peut comprendre :
a. Détermination de paramètres moteur nécessaires à la réalisation d’un calcul d’équilibre roulement, étant un indice entier variant entre 1 et ,
b. Détermination d’une loi de variation de chacun des paramètres moteurs ,
c. Génération de valeurs pour chacun des paramètres en utilisant la loi de variation obtenue à l’étape précédente.
L’étape E4 peut être réalisée par un calcul d’équilibre du roulement permettant d’obtenir :
i. la pression de contact du corps roulant sur les pistes des bagues interne et externe, en MPa
ii. le produit PV qui correspond à la pression multipliée par la vitesse de glissement corps roulant et les pistes des bagues interne et externe : unité en W/mm²,
iii. le tenseur de contrainte en surface,
iv. le tenseur de contrainte en sous-couche,
v. la contrainte de traction liée au frettage et à la force centrifuge sur les bagues interne et externe en MPa, et
vi. l’épaisseur du film d’huile noté h0 entre le corps roulant et les pistes des bagues interne et externe en microns
L’étape E5 peut comprendre :
a) La réalisation d’une estimation de la durée de vie en sous-couche pour chacun desdits composants du roulement à partir du modèle d’endommagement en sous-couche et de paramètres de sortie d’un calcul d’équilibre du roulement pour chaque composant unitaire du roulement à savoir la bague externe, la bague interne et le corps de roulant respectivement noté , et , à savoir :
1. correspondant à la durée de vie en sous-couche pour la bague interne,
2. correspondant à la durée de vie en sous-couche pour la bague externe,
3. correspondant à la durée de vie en sous-couche pour le corps de roulant,
De manière à obtenir N triplés de valeurs [ ],
b) déterminer le minimum de chaque triplé de valeurs de manière à obtenir N valeurs minimales de durée de vie en sous couche.
L’étape E6 peut comprendre :
a) La réalisation d’une estimation de la durée de vie en en surface pour chacun desdits composants du roulement à partir du modèle d’endommagement en sous-couche et de paramètres de sortie d’un calcul d’équilibre du roulement pour chaque composant unitaire du roulement à savoir la bague externe, la bague interne et le corps roulant respectivement noté , et , à savoir :
1. correspondant à la durée de vie en surface pour la bague interne,
2. correspondant à la durée de vie en surface pour la bague externe,
3. correspondant à la durée de vie en surface pour le corps roulant,
De manière à obtenir de N triplés de valeurs [ ],
b) déterminer le minimum de chaque triplé de valeurs de manière à obtenir N valeurs minimales de durée de vie en surface
Le procédé peut encore comprendre :
- Fourniture d’un taux
- Définition de paramètres
-
Les paramètres et peuvent être déterminés de manière à ce que :
Où avec qui représente le taux de défaillance maximum du roulement
La durée de vie cible peut être égale à au moins 20000 heures, préférentiellement égale à 20000 heures.
D’autres caractéristiques, détails et avantages apparaîtront à la lecture de la description détaillée ci-après, et à l’analyse des dessins annexés, sur lesquels :
La représente une turbomachine 10 d’un type connu et comprenant d’amont en aval, c’est à dire selon le sens d’écoulement d’un écoulement d’air au sein de la turbomachine, une soufflante 12 comportant une roue de soufflante réalisant une accélération d’un écoulement d’air entrant à l’intérieur de la turbomachine. Cet écoulement d’air se divise en un écoulement annulaire d’air primaire et un écoulement annulaire d’air secondaire entourant l’écoulement d’air primaire. L’écoulement annulaire d’air primaire traverse un compresseur 14 basse pression puis haute pression, une chambre de combustion 16, une turbine 18 haute pression puis basse pression.
Le compresseur 14 et la turbine 18 comprennent des rotors dont la rotation de chacun est guidée par l’intermédiaire de roulements de différents types, à rouleaux ou à billes par exemple. Classiquement, chaque roulement comprend un bague annulaire interne et une bague annulaire externe, l’une des bagues interne et externe étant solidaire d’un rotor et l’autre des bagues interne et externe étant solidaire d’une pièce statorique ou bien encore d’une autre pièce mobile lorsque le roulement est agencé entre deux pièces mobiles en rotation à des vitesses différentes et/ou des sens différents.
Une turbomachine comporte de nombreux roulements dont on comprend aisément que la maitrise de la durée de vie c’est-à-dire la fiabilité est primordiale au bon fonctionnement de la turbomachine.
En fonctionnement, deux types majeurs de défauts peuvent intervenir sur un roulement sur les composants d’un roulement. Un premier défaut est l’écaillage en profondeur qui consiste en des défauts de sous-couches. En pratique, on observe un écaillage c’est-à-dire une fissuration en profondeur dans la pièce. Un second défaut est l’écaillage de profondeur initié en surface qui consiste en un écaillage ou une fissuration créée au voisinage d’un défaut créés par un indent qui est un défaut consécutif au passage d’une ou plusieurs particules dans le roulement entre le corps roulant et la bague interne et la bague externe.
Principalement, les deux défauts précités sont observés sur la bague interne, la bague externe et le corps roulant c’est-à-dire le corps réalisant roulant entre la bague interne et la bague externe.
Comme évoqué précédemment, la détermination des paramètres de dimensionnement et des matériaux des roulements sont réalisés avec une méthodologie non adaptée aux conditions de durée de vie souhaitées dans le monde aéronautique.
A cette fin, il est proposé une méthode de détermination des paramètres de dimensionnement et de matériaux d’un roulement.
Le procédé est illustré par la séquence d’étapes successives représentées en .
Le procédé comprend neuf étapes E1, E2, E3, E4, E5, E6, E7, E8 et E9 qui ne sont pas nécessairement successives selon la séquence successive de un à neuf. Ainsi, il doit être compris que le terme « étape » en relation avec un numéro « X » donné ne signifie pas ladite étape doit être réalisé nécessairement après l’étape X-1 et avant l’étape X+1.
Plus particulièrement, comme cela apparaitra à la lecture de la description :
- Les étapes E1 et E2 sont des étapes qui peuvent être réalisées simultanément ou bien l’une après l’autre. Les données obtenues à l’issue de ces étapes E1 et E2 sont utilisées pour réaliser l’étape E3. Ainsi, il doit être compris que les étapes les étapes E1 et E2 doivent être réalisées préalablement à l’étape E3.
- L’étape E4 est réalisée après l’étape E3 et avant les étapes E5 et E6.
- Les étapes E5 et E6 peuvent être réalisées simultanément ou bien l’une après l’autre. Les données obtenues à l’issue des étapes E5 et E6 sont utilisées pour réaliser l’étape E7, cette dernière étape E7 étant donc réalisée après les étapes E5 et E6.
- Les étapes E7, E8 et E9 sont réalisées successivement l’une après l’autre et dans cet ordre.
En préalable à la mise en œuvre des étapes E1 à E9, telles qu’illustrées en , le procédé comprend :
- Fourniture d’un taux de défaillance maximum du roulement
- Définition de paramètres , d’une simulation numérique Monte-carlo permettant d’évaluer ledit taux de défaillance maximum, où :
Ces paramètres et sont des paramètres qui peuvent être utilisés dans une simulation de Monte-Carlo.
Où avec qui représente le taux de défaillance maximum du roulement,
A titre d’exemple, lorsque l’on souhaite atteindre un taux de défaillance maximum de 10-7, le nombre d’heures simulées qui est alors égal à 109ne peut être supérieur au produit afin de vérifier l’équation ci-dessus. Ainsi, si on fixe la durée de vie cible à atteindre à 20000 heures, alors il convient d’avoir un paramètre qui doit être d’au moins 50000 couples de séries .
Nous allons maintenant décrire les étapes E1 et E2.
L’étape E1 consiste ainsi à :
- Déterminer tous les
- Obtenir
L’homme du métier connait bien les paramètres géométriques, de matériaux et de contrainte résiduelle nécessaires à la réalisation d’un calcul d’équilibre de roulement permettant d’obtenir en sortie les paramètres de pression.
A titre d’exemple, pour une application de roulement à billes, les paramètres géométriques, matériaux et de contrainte résiduelle peuvent comprendre : le nombre de billes, le diamètre nominal de chaque bille, le diamètre moyen du roulement calculé du centre du roulement au centre d’une bille, le jeu diamétral, la courbure de la face radialement externe de la bague interne, la courbure de la face radialement interne de la bague externe du roulement, le régime de rotation de la bague interne.
A titre d’exemple, pour une application de roulement à rouleaux, les paramètres géométriques, de matériaux et de contrainte résiduelle peuvent comprendre : le nombre de rouleaux, le diamètre externe de la partie cylindrique de chaque rouleau, la longueur totale de chaque rouleau, le ou les rayons de raccordement de la partie cylindrique de chaque rouleau au flancs du rouleau, la longueur de la partie cylindrique du rouleau, le diamètre moyen calculé du centre du roulement au centre d’un rouleau, la courbure de la face radialement externe de la bague interne, la courbure de la face radialement interne de la bague externe du roulement, ou encore les modules d’Young, de poisson, les contraintes résiduelles.
Cette étape E1 comprend les sous-étapes suivantes :
- Détermination de
- Détermination d’une loi de variation de chacun des paramètres
- Génération de
L’étape b) ci-dessus est réalisée en collectant N valeurs des paramètres . Ainsi, en relation avec l’exemple donné précédemment d’un taux de défaillance de 10-7, pour un paramètre donné tel que par exemple le nombre de billes pour un roulement à billes, il conviendra de collecter au moins N=50000 valeurs distinctes de chaque de nombre de billes.
On comprend ainsi qu’il faudra établir une liste de N valeurs de chacun des paramètres . Les données obtenues peuvent être rangées dans un tableau ayant la forme générale suivante, le tableau étant également exemplifié pour un roulement à billes :
(ex : 1 représente le nombre de billes) |
… |
|
… | (ex : S représente le diamètre nominal de chaque bille) |
|
|
|||||
… | |||||
|
|||||
… | |||||
(ex:N=50000) |
L’étape E2 consiste ainsi à :
- Déterminer
- Obtenir
L’homme du métier connait bien les paramètres moteur nécessaires à la réalisation d’un calcul d’équilibre de roulement permettant d’obtenir en sortie les paramètres de pression.
A titre d’exemple, pour une application de roulement à billes ou de roulements à rouleaux, les paramètres moteurs peuvent comprendre par exemple : pour le ou les points de fonctionnement considérés le régime de rotation de la bague interne, le régime de rotation de la bague externe, la charge axiale admissible par le roulement, la charge radiale admissible par le roulement.
Cette étape E2 comprend les sous-étapes suivantes :
- Détermination de
- Détermination d’une loi de variation de chacun des paramètres moteurs
- Génération de
L’étape b) du paragraphe ci-dessus est réalisée en collectant N valeurs des paramètres . Ainsi, en relation avec l’exemple donné précédemment d’un taux de défaillance de 10-7, pour un paramètre donné tel que par exemple le régime de rotation de la bague interne, il conviendra de collecter au moins N=50000 valeurs distinctes de régime de rotation de la bague interne[question pratique : A quoi correspondraient ces N régimes de rotation, ceux sont des valeurs réelles ?].
On comprend ainsi qu’il faudra établir une liste de N valeurs de chacun des paramètres . Les données obtenues peuvent être rangées dans un tableau ayant la forme générale suivante, le tableau étant également exemplifié pour un roulement à billes :
En pratique, on a autant de tableaux que de points de mission. On pourrait utiliser une technique dite de rainflow pour obtenir un seul tableau
(ex : 1 représente la force axiale) | … | | … | (ex : S représente la force radiale) | |
| |||||
… | |||||
| |||||
… | |||||
(ex:N=50000) |
Après l’obtention des données concernant les roulements et relatives au fonctionnement desdits roulements dans le moteur, le procédé consiste à réaliser une étape E3 dans laquelle on associe aléatoirement une série de valeurs à une série de valeurs de manière à obtenir N couples de séries de valeurs des paramètres / .
Dans une étape E4, pour chacun des couples de séries , on réalise un calcul d’équilibre du roulement en résolvant les équations d’équilibre du roulement comme cela est illustré aux figures 2 et 3 ce qui permet d’obtenir les paramètres de sortie dudit calcul d’équilibre du roulement, à savoir
- la pression entre les éléments roulants et les pistes
- le produit PV qui correspond à la pression multipliée par la vitesse entre les éléments roulants et les pistes,
- le tenseur des contraintes en surface,
- le tenseur des contraintes en sous-couche,
- la contrainte de traction liée au frettage et à la force centrifuge, et
- l’épaisseur du film d’huile noté h0 entre éléments roulants et les pistes de roulement
Ainsi, chaque série de paramètres d’entrée est associé à une série de paramètres de sortie (i à v ci-dessus) pour chaque cas d’usage qui varie de 1 à et jusqu’à 50000 dans l’exemple donné précédemment ( ).
On notera que les équations d’équilibre du roulement ainsi que le principe mathématique d’équilibrage d’un roulement sont décrits dans l’ouvrageRolling Bearing Analysis, 5th Edition, de Tedric A. Harris & Michael N. Kotzalas.
Les données de sorties peuvent se présenter ainsi sous la forme du tableau suivant (voir également ) :
Pression | Produit PV | Tenseur de contrainte en surface | Tenseur de contrainte en sous couche | La contrainte de traction liée au frettage et à la force centrifuge | L’épaisseur du film d’huile noté h0 | |
| ||||||
… | ||||||
| ||||||
… | ||||||
(ex:N=50000) |
On notera que le procédé peut faire appel aux connaissances de l’utilisateur de sorte que celui-ci peut aussi utiliser ses connaissances empiriques pour estimer que les paramètres de sortie. Par exemple, il peut être considérer que les paramètres en entête de colonne du tableau immédiatement ci-dessus varie selon une loi normale particulière de sorte que l’estimation des N valeurs de chaque paramètre peut être réalisée plus rapidement. Cette méthode s’avère toutefois moins précise que le calcul des paramètres par l’intermédiaire d’un calcul d’équilibre du roulement mais plus rapide.
Les étapes D1, E2, E3 et E4 précédemment décrites sont réalisées pour chacun des composants unitaires dudit roulement à savoir la bague interne BI, la bague externe BE et le corps roulant CR. On obtient ainsi valeurs de paramètres pour la bague interne noté , valeurs de paramètres pour la bague externe noté et valeurs de paramètres pour le corps roulant noté comme cela est illustré en .
Afin de réaliser un calcul de durée de vie pour chaque composant unitaire de l’organe de roulement, à savoir la bague interne, la bague externe et le corps roulant, on établit un modèle d’endommagement suivant le premier défaut décrit précédemment à savoir un défaut en sous couche. Un modèle d’endommagement en sous couche est établi pour chacun de la bague interne, la bague externe et le corps roulant ou roulant.
Un modèle d’endommagement, en sous-couche ou bien en surface par exemple, est réalisé préalablement. A cette fin, une série d’au moins roulements où est obtenue. Pour chaque roulement, il est réalisé un test de fatigue de contact qui est définie comme une défaillance ou un enlèvement de matière dû à la propagation d'une fissure causée par un champ de contraintes alternatif proche de la surface.
Lors de cette campagne d’essais, pour chaque composant unitaire du roulement, à savoir la bague externe, la bague interne et le corps roulant, on mesure ou calcule les paramètres suivants :
- La contrainte résiduelle dans la direction de roulement et la direction orthogonale,
- La contrainte consécutive au frettage des bagues interne et externe,
- La pression, et
- Le tenseur de contrainte en sous-couche (ou surface)
L’homme du métier pourra utiliser les deux formules suivantes pour calculer dans un premier temps la contrainte totale puis dans un second temps la contrainte équivalente de Tresca.
S’agissant de la contrainte totale pour ce qui concerne l’endommagement en sous-couche, celle-ci est donnée par l’équation suivante :
Où représente le tenseur de contraintes exprimés en MPa.
La première partie de l’équation ci-dessus correspond à la contrainte générée par le problème de contact, la deuxième partie correspond la contrainte consécutive au frettage des bagues interne et externe, la troisième correspond à la contrainte résiduelle.
S’agissant de la contrainte totale pour ce qui concerne l’endommagement en surface, celle-ci est donnée par l’équation suivante :
Où représente le tenseur de contraintes exprimés en MPa,
hb est la hauteur de bourrelet,
h0 l’épaisseur de film d’huile, et
La première partie de l’équation ci-dessus correspond à la contrainte générée par le problème de contact, la deuxième partie correspond à la contribution de l’indent, la troisième partie à la contrainte consécutive au frettage des bagues interne et externe, la quatrième correspond à la contrainte résiduelle.
Par convention tensorielle, l’indice 1 correspond à la direction circonférentielle de roulement, l’indice 2 correspond à une direction perpendiculaire de roulement et l’indice 3 correspond à la direction radiale perpendiculaire aux première et deuxième directions.
L’équation de la contrainte équivalente de Tresca (commune aux modèles d’endommagement en sous-couche et modèles d’endommagement en surface) est obtenue selon :
On remarquera que la contrainte équivalent de Tresca permet d’obtenir une valeur scalaire à partir du tenseur de contraintes .
Pour chaque composant unitaires du roulement parmi les S roulements testés, on détermine le nombre de cycles à défaillance, ce qui permet d’obtenir une loi ou modèle d’endommagement associé à chaque composant unitaire. Par exemple, pour la bague interne, l’obtention de la loi d’endommagement est ainsi réalisée en déterminant pour une valeur de contrainte donnée le nombre de cycles au bout duquel une défaillance intervient. Le nombre de cycles correspondant à une défaillance correspond à une fréquence de sollicitation de composant unitaire à une contrainte donnée. Chaque loi d’endommagement pourra être améliorée au fur et à mesure de compagnes d’essais.
La loi d’endommagement prend la forme d’une courbe ayant la forme suivante représentée en . On comprend qu’au fur et à mesure que la contrainte augmente, le nombre de cycles à défaillance est diminué.
Par suite, ayant établi un modèle d’endommagement en sous-couche pour chacun de la bague interne BI, la bague externe BE et le corps roulant CR, il est possible de déterminer la durée de vie durée de vie en sous-couche pour chacun desdits composants du roulement comme cela est illustré en au travers de l’étape E5a.
On obtient ainsi :
Dans une étape E5b, on détermine le minimum de chaque triplé de valeurs de manière à obtenir N valeurs minimales de durée de vie en sous couche.
De même, ayant établi un modèle d’endommagement en surface pour chacun de la bague interne BI, la bague externe BE et le corps roulant CR, il est possible de déterminer la durée de vie durée de vie en surface pour chacun desdits composants du roulement comme cela est illustré en au travers de l’étape E6a.
On obtient ainsi :
Dans une étape E6b, on détermine le minimum de chaque triplé de valeurs de manière à obtenir N valeurs minimales de durée de vie en surface.
Dans une étape E7 illustré en figure 6, on détermine valeurs de durée de vie global selon : (étape E7, ).
Chacune des N valeurs est comparée à une durée de vie cible enregistrée en base de données (étape E8, figure 6) et on établit une liste de valeurs de défaillance positive (oui) ou négative (non) (étape E9, ).
Ensuite, la fiabilité est obtenue en déterminant le rapport du nombre défaillances sur le nombre d’heures qui correspond à la durée de vie cible. Dans le cas d’un moteur d’avion la durée de vie cible est d’au moins 20000 heures.
En pratique, pour réaliser la simulation, on réalisera :
- Fourniture d’un taux
- Définition de paramètres
-
Où avec qui représente le taux de défaillance maximum du roulement
Si la méthode décrite précédemment peut être mise en œuvre avec les deux seuls modèles d’endommage de surface et modèle d’endommagement de sous-couche, on comprend que l’on pourrait ajouter d’autres modes d’endommagements tels que par exemple un écaillage en surface initiée en surface. Toutefois, on notera que dans l’aéronautique, les deux modes d’endommagement évoqués précédemment s’avèrent les plus importants et doivent donc être pris en compte dans une étude et une définition de la fiabilité d’un nouveau roulement.
Claims (9)
- Procédé de test de la fiabilité d’un roulement pour un moteur, tel qu’un moteur de turbomachine, comprenant les étapes suivantes :
E1)- Déterminer
- Obtenir
- Déterminer
- Obtenir
E4) Pour chacun des
E5) Etablir un – d’endommagement d’une sous-couche des composants unitaires du roulement, à savoir la bague externe, la bague interne et le corps roulant puis détermination de N valeurs minimales
E6) Etablir d’un modèle d’endommagement de surface des composants unitaires du roulement, à savoir la bague externe, la bague interne et le corps roulant puis détermination de N valeurs minimales
E7) Déterminer une valeur de durée de vie globale
E8) Comparer chacune des N valeurs de
E9) Déterminer un taux de défaillance. - Déterminer
- Procédé selon l’une des revendications précédentes, dans lequel l’étape E1 comprend :
a. Détermination de
b. Détermination d’une loi de variation de chacun des
c. Génération de - Procédé selon l’une des revendications précédentes, dans lequel l’étape E2 comprend :
a. Détermination de
b. Détermination d’une loi de variation de chacun des paramètres moteurs
c. Génération de - Procédé selon l’une des revendications précédentes, dans lequel l’étape E4 est réalisée par un calcul d’équilibre du roulement permettant d’obtenir :
i. la pression de contact du corps roulant sur les pistes des bagues interne et externe , en MPa
ii. le produit PV qui correspond à la pression multipliée par la vitesse de glissement corps roulant et les pistes des bagues interne et externe : unité en W/mm²,
iii. le tenseur de contrainte en surface,
iv. le tenseur de contrainte en sous-couche,
v. la contrainte de traction liée au frettage et à la force centrifuge sur les bagues interne et externe en MPa, et
vi. l’épaisseur du film d’huile noté h0 entre le corps roulant et les pistes des bagues interne et externe en microns - Procédé selon l’une des revendications précédentes, dans lequel l’étape E5 comprend :
a) La réalisation d’une estimation de la durée de vie en sous-couche pour chacun desdits composants du roulement à partir du modèle d’endommagement en sous-couche et de paramètres de sortie
b) déterminer le minimum -
- Procédé selon l’une des revendications précédentes, dans lequel l’étape E6 comprend :
a) La réalisation d’une estimation de la durée de vie en en surface pour chacun desdits composants du roulement à partir du modèle d’endommagement en sous-couche et de paramètres de sortie
b) déterminer le minimum -
- Procédé selon l’une des revendications précédentes, dans lequel il comprend :
- Fourniture d’un taux
- Définition de paramètres
-
- Fourniture d’un taux
- Procédé selon la revendication précédente, dans lequel les paramètres
Où - Procédé selon l’une des revendications précédentes, dans lequel la durée de vie cible est égale à au moins 20000 heures, préférentiellement égale à 20000 heures.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR2210599A FR3140944A1 (fr) | 2022-10-14 | 2022-10-14 | Procédé de détermination de la fiabilité d’un roulement |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR2210599 | 2022-10-14 | ||
FR2210599A FR3140944A1 (fr) | 2022-10-14 | 2022-10-14 | Procédé de détermination de la fiabilité d’un roulement |
Publications (1)
Publication Number | Publication Date |
---|---|
FR3140944A1 true FR3140944A1 (fr) | 2024-04-19 |
Family
ID=86468799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR2210599A Pending FR3140944A1 (fr) | 2022-10-14 | 2022-10-14 | Procédé de détermination de la fiabilité d’un roulement |
Country Status (1)
Country | Link |
---|---|
FR (1) | FR3140944A1 (fr) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1780501B1 (fr) * | 2005-10-25 | 2008-12-10 | Snecma | Procédé de contrôle du profil de la zone de raccordement entre la partie cylindrique et la dépouille d'une pièce de turbomachine |
EP3855153A1 (fr) * | 2020-01-23 | 2021-07-28 | Siemens Aktiengesellschaft | Procédé pour déterminer une durée de vie utile restante d'un palier |
-
2022
- 2022-10-14 FR FR2210599A patent/FR3140944A1/fr active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1780501B1 (fr) * | 2005-10-25 | 2008-12-10 | Snecma | Procédé de contrôle du profil de la zone de raccordement entre la partie cylindrique et la dépouille d'une pièce de turbomachine |
EP3855153A1 (fr) * | 2020-01-23 | 2021-07-28 | Siemens Aktiengesellschaft | Procédé pour déterminer une durée de vie utile restante d'un palier |
Non-Patent Citations (3)
Title |
---|
LUNDBERG GPAMGRE, A: "Dynamic capacity of rolling bearings", J APPL MECH, vol. 16, no. 2, 1949, pages 165 - 72 |
LUNDBERG GUSTAF ET AL: "Dynamic Capacity of Rolling Bearings", JOURNAL OF APPLIED MECHANICS., vol. 16, no. 2, 1 June 1949 (1949-06-01), US, pages 165 - 172, XP093082810, ISSN: 0021-8936, Retrieved from the Internet <URL:http://asmedigitalcollection.asme.org/appliedmechanics/article-pdf/16/2/165/6746260/165_1.pdf> DOI: 10.1115/1.4009930 * |
OMAR DJEBILI: "Contribution à la maintenance prédictive par analyse vibratoire des composants mécaniques tournants. Application aux butées à billes soumises à la fatigue de contact de roulement", 26 September 2013 (2013-09-26), XP055254003, Retrieved from the Internet <URL:http://www.theses.fr/2013REIMS030.pdf> [retrieved on 20160229] * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2749214C (fr) | Procede et systeme de surveillance de phenomenes vibratoires survenant dans un moteur a turbine a gaz d'aeronef en fonctionnement | |
EP2496921B1 (fr) | Procédé de détection d'un endommagement d'au moins un roulement de palier d'un moteur | |
Ahmadi et al. | The path of rolling elements in defective bearings: observations, analysis and methods to estimate spall size | |
CA2623617C (fr) | Procede de detection d'un endommagement d'un roulement de palier d'un moteur | |
FR2932266A1 (fr) | Procede de prediction du comportement dynamique d'une structure d'un aeronef | |
FR2882393A1 (fr) | Procede et dispositif pour assembler des machines tournantes | |
FR2584496A1 (fr) | Procede et appareil de detection de fissures axiales dans des rotors de machines tournantes. | |
WO2016075409A1 (fr) | Procédé de surveillance d'un moteur d'aéronef en fonctionnement dans un environnement donné | |
CA2830834C (fr) | Procede d'inspection des impacts observes dans des carters de soufflante | |
CA3116655A1 (fr) | Dispositif et procede de surveillance de duree de vie d'un equipement hydraulique d'un aeronef | |
FR2971870A1 (fr) | Methode de detection de pannes d'une turbomachine au moyen d'un modele theorique du cycle thermodynamique de ladite turbomachine | |
FR3140944A1 (fr) | Procédé de détermination de la fiabilité d’un roulement | |
CA2984522C (fr) | Procede de surveillance d'un moteur d'aeronef en fonctionnement pendant un vol | |
EP2872956A1 (fr) | Méthode de détection d'une dégradation d'une turbomachine par surveillance des performances de ladite turbomachine | |
WO2021044098A1 (fr) | Procédé de relaxation des contraintes par rotation | |
FR3039676A1 (fr) | Procede de determination d'une contrainte alternee maximale admissible pour une piece soumise a des sollicitations cycliques; unite de determination d'une telle contrainte | |
FR3059784A1 (fr) | Installation de diagnostic d'une machine electrique | |
US20200182730A1 (en) | Component inspection using measured and theoretical moment weights | |
FR3125335A1 (fr) | Méthode de réguler les contraintes résiduelles dans les pièces forgées du disque vierge du superalliage par pré-rotation | |
WO2022128717A1 (fr) | Procédé et système de détermination d'un ou plusieurs défauts d'une machine tournante d'aéronef | |
EP3066445A1 (fr) | Procédé d'essai technique | |
Tetter et al. | Researching the Possibility of Determining the Technical Condition of Rolling Bearings by the Kurtosis Factor | |
FR3123718A1 (fr) | Dispositif de determination des modes propres d’un element de turbomachine par emission d’un bruit blanc et dispositif de caracterisation vibratoire de l’element | |
CN114295375A (zh) | 一种滚动轴承极限转速试验方法 | |
FR3138555A1 (fr) | Procédé de détermination de la température des gaz de sortie d’un moteur et procédé de mesure de l’usure d’un moteur |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20240419 |