These
These
Etude de la performance
énergétique d’une transmission
de puissance haute vitesse
_____________________________________________________________
1
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Résumé
Etude de la performance énergétique d’une transmission de puissance
haute vitesse
Une des voies d’amélioration des véhicules hybrides et électriques est l’utilisation de moteurs
tournant plus vite, jusqu’à plus de 42.000tr/min. Le but est d’augmenter la densité de puissance et
le rendement des groupes motopropulseurs. Pour utiliser ces moteurs de nouveaux réducteurs
mécaniques doivent être développés. Cela doit se faire sans générer de surcoût important face aux
solutions utilisées à basse vitesse et en assurant un niveau de performance énergétique élevé.
Cette thèse se situe en amont de la phase de conception d’un réducteur haute vitesse lubrifié
par barbotage. Elle a pour but d’identifier les problèmes d’échauffement et de pertes de puissance
ainsi que de proposer des pistes d’amélioration énergétique.
Ce travail propose la modélisation thermomécanique de l’étage grande vitesse (GV) du
réducteur, réalisée à l’aide de la méthode des réseaux thermiques. Ce modèle couple les pertes de
puissance avec les températures. Une attention particulière est portée sur la modélisation des
roulements de l’arbre GV. Un nouveau modèle thermomécanique de roulement est développé.
Les pertes par barbotage deviennent importantes à grande vitesse. Une méthode permettant de
fortement les réduire est caractérisée.
Mots clés : Engrenages, roulement à billes, rendement, modélisation thermomécanique,
barbotage
3
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Abstract
Study of the efficiency of a high-speed mechanical power transmission
A way to improve both electric and hybrid vehicles is to use high speed motors, operating
over than 42.000rpm. The goal is to increase the power density and the efficiency of powertrains.
Using these new motors, new gearboxes should be developed. This must be done without
generating significant additional cost regarding already mastered low speed solutions. High energy
performance level also has to be maintained.
This PhD comes before the design phase of a high-speed oil bath lubricated gearbox. It aims
to identify the warm-up and power loss problems, and propose ways to improve efficiency.
This work proposes a thermomechanical modelling of the gearbox’s first stage, using the
thermal network method. This model links power losses with temperatures. Particular attention is
paid to high speed bearing modelling. A new thermomechanical model of rolling element bearing
is developed. As churning losses being significant at high speeds, a method to greatly reduce this
power loss is characterized.
Key words: Gears, ball bearing, efficiency, thermomechanical modelling, oil bath lubrication
5
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Nomenclature
𝑎 Demi-axe de l’ellipse de contact dans le sens du roulement m
𝑎𝑑 Diffusivité thermique m².s-1
a’ et b’ Facteurs de correction fonction du type de roulement
A Maître couple m²
𝐴𝑣 et Coefficients pour le calcul de la viscosité
𝐵𝑣
𝐴% Aération de l’huile %
𝑏 Demi-axe de l’ellipse de contact dans le sens perpendiculaire du m
roulement
𝐵 Largeur m
[𝐶] Matrice des inerties thermiques J/K
𝐶𝐷 Coefficient de traînée
𝐶𝑏𝑎𝑟𝑏 Couple de barbotage N.m
𝐶𝑓 Couple de frottement d’un palier fluide N.m
𝐶𝑖 Chaleur massique de l’élément i J/kg.K
𝐶𝑚 Couple adimensionnel de trainée
𝐶𝑝 Capacité calorifique massique J/kg.K
𝑑 Diamètre intérieur m
𝑑𝑚 Diamètre moyen m
𝐷 Diamètre extérieur m
𝐷𝑝 Diamètre primitif m
𝑒 Epaisseur m
𝑒𝑐 Epaisseur de la cage m
𝑓 Coefficient de frottement
𝑓0 Facteur fonction du type de roulement et du type de lubrification
𝑓1 Facteur fonction de la charge statique équivalente et statique de base,
et du type de roulement
𝐹12 Facteur de forme de 1 vers 2 pour le rayonnement
𝐹𝑅 Force radiale N
𝐹𝑛 Effort normal N
𝐹𝑛𝑢 Effort normal unitaire N/m
𝑔 Accélération de la pesanteur = 9,81 m/s²
𝑔𝑎 Longueurs de retrait m
𝑔𝑓 Longueurs d’approche m
𝐺∗ Paramètre adimensionné de matériau = 𝛼𝑏 𝐸 ′
𝐺𝑟𝑟 Facteur fonction des charges, des dimensions, et du type de roulement
𝐺𝑠𝑙 Facteur fonction des charges, des dimensions, et du type de roulement
ℎ Hauteur immergée m
ℎ𝑐𝑜𝑛𝑣 Coefficient de convection W/m3
𝐻𝑉 Facteur géométrique de perte denture
𝐻𝑑𝑒𝑛𝑡 Hauteur des dents m
𝐻𝑓 Hauteur de flasque m
𝐾𝑟𝑠 Constante fonction du type de lubrification
𝐾𝑧 Constante géométrique fonction du type de roulement
7
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
𝐽𝑎 Jeu axial m
𝐽𝑟 Jeu radial m
k Conductivité thermique W.K-1.m-1
𝑙ℎ Demi largeur de contact linéaire m
𝑙𝑐𝑚 Longueur moyenne de contact m
L Longueur m
𝑚 Module normal m
𝑚̇ Débit massique kg/s
𝑚𝑖 Masse de l’élément i kg
𝑀0 Couple indépendant de la charge N.m
𝑀1 Couple dépendant de la charge N.m
𝑀𝑑𝑟𝑎𝑔 Couple de trainée dû au mode de lubrification N.m
𝑀𝑟𝑟 Couple causé par le roulement N.m
𝑀𝑠𝑙 Couple causé par le glissement N.m
𝑀𝑡𝑜𝑡𝑎𝑙 Couple de frottement total d’un roulement N.m
𝑛 Nombre d’élément
𝑁 Vitesse de rotation tr/min
𝑝𝑏𝑡 Pas de base apparent m
𝑃 Puissance mécanique W
𝑃0 Pression dans le contact Pa
𝑃1 Combinaison des charges axiale et radiale sur le roulement N
𝑄 Puissance thermique [W]
{𝑄} Vecteur des puissances thermiques W
𝑟𝑐 Rayon de courbure équivalent au niveau du diamètre primitif m
𝑅 Rayon m
𝑅𝑎𝑟 Moyenne arithmétique des rugosités des flancs actifs du pignon et de la μm
roue
𝑅𝑏 Rayon moyen des bulles d’air m
𝑅𝑝 Rayon primitif m
𝑅𝑡ℎ Résistance thermique K/W
𝑅𝑥 Rayon de courbure équivalent dans le sens du roulement m
𝑆 Surface m3
𝑆𝑅𝑅 Taux de glissement
[𝑆] Matrice des conductances thermiques W/K
𝑆𝑚 Surface immergée du pignon m2
𝑆𝑚𝑑 Surface mouillée de la denture
𝑆𝑚𝑓 Surface mouillée des flancs du pignon m²
𝑇 Température K
⃗ 1 (𝑀)
𝑇 Force élémentaire tangentielle de frottement au point de contact M N/m
2 entre les deux flancs
{𝑇} Vecteur des températures K
{𝑇 ′ } Vecteur des températures dérivées par rapport au temps K/s
𝑍
𝑢 Rapport de réduction de la transmission = 𝑍1 𝑎𝑣𝑒𝑐 𝑍1 ≤ 𝑍2
2
𝑈 Vitesses de roulement (𝑈1 + 𝑈2 ) m/s
𝑈∗ Paramètre adimensionné de vitesse = 𝑈. 𝜈/(2. 𝐸 ′ . 𝑅𝑥 )
𝑈1 Vitesse du contact par rapport au solide 1 m/s
𝑈2 Vitesse du contact par rapport au solide 2 m/s
𝑉 Vitesse m/s
𝑉0 Volume d’huile m3
𝑉𝑎𝑖𝑟 Vitesse du flux d’air m/s
𝑉𝑔 Vitesse de glissement (𝑈1 − 𝑈2 )/2 m/s
8
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
𝑉𝑚𝑜𝑦 Vitesse moyenne m/s
𝑉𝑝𝑖𝑔𝑛𝑜𝑛 Volume immergé du pignon m3
𝑊 Charge normale N
𝑊∗ Paramètre adimensionné de charge = 𝑊/(𝐸 ′ . 𝑅𝑥 )
𝑋𝐿 Coefficient dépendant du type d’huile utilisé
⃗⃗⃗⃗
𝑉12 (𝑀) Vitesse de glissement au point M m/s
𝑍 Nombre de dents rad
Symboles grecs
𝛼 Angle de pression normale rad
𝛼𝑏 Coefficient de piézoviscosité 1/Pa
𝛼𝑡 Angle de pression apparent rad
𝛽 Angle d’hélice rad
𝛽𝑏 Angle d’hélice de base
𝛽𝑐 Inclinaison de l’axe de rotation de la bille par rapport à l’axe de rad
rotation du roulement
𝛾 Accélération centrifuge m/s²
𝛤 Critère d’efficacité des flasques
𝛿 Décalage du centre d’application de charge N
𝛥𝑚 Valeur moyenne d’approche normale par rapport à la position de m
corps rigide
∆𝑇 Différence de température K
𝜀 Emissivité du matériau, comprise entre 0 et 1
𝜀𝛼 Rapport de conduite apparent
𝜂 Efficacité d’ailettes
𝜃 Angle d’immersion rad
𝜃𝑝 Angle de projection rad
𝜅0 Position adimensionnel du point de roulement sans glissement le
long de la ligne d’action
𝜇 Viscosité dynamique Pa.s
𝜇0 Viscosité dynamique en entrée de contact Pa.s
𝛬 Facteur de perte
𝜎𝑠 Tension de surface du lubrifiant N/m
𝛶 Profondeur de correction de profil adimensionnée par rapport à Δ𝑚
𝜈 Viscosité cinématique m²/s
𝜌 Masse volumique kg/m3
𝜎 Constante de Stefan-Boltzmann = 5,67. 10−8 W/m².K4
𝜏 Variable de temps normalisée par rapport à la période d’engrènement
𝜒 Effusivité thermique W⋅K-1⋅m-2⋅s1/2
𝜑𝑖𝑠ℎ Facteur de la réduction de chaleur due au cisaillement d'entrée
𝜑𝑟𝑠 Facteur de réduction de l'appoint/insuffisance cinématique
𝜓 Nombre sans dimension dépendant de la géométrie et de la viscosité
de l’huile
𝜔 Vitesse de rotation rad/s
Nombres adimensionnels
𝐹𝑟 Nombre de Froude
𝐺𝑟 Nombre de Grashof
9
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
𝑁𝑢 Nombre de Nusselt
𝑃𝑒 Nombre de Peclet
𝑃𝑟 Nombre de Prandtl
𝑅𝑒 Nombre de Reynolds
10
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Sommaire
11
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
4.4 DEVELOPPEMENT D’UN MODELE ........................................................................112
4.5 APPLICATION........................................................................................................119
4.6 CONCLUSION ....................................................................................................... 122
12
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Introduction générale
14
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Le second chapitre s’attache à analyser les problèmes que pose la conception d’un réducteur
automobile acceptant une vitesse jusqu’à 42 000 tr/min en entrée. Pour répondre au cahier des
charges du réducteur, deux architectures utilisant des engrenages à axes parallèles sont
envisageables ; soit utiliser deux étages en augmentant les rapports de réduction, soit limiter les
rapports de réduction et utiliser 3 étages. Le problème est l’échauffement de l’étage grande
vitesse, il est identifié comme majeur. Ne voyant qu’un faible couple cet étage peut être très
compact, cependant il transmet toute la puissance du moteur. Un modèle thermomécanique de
l’étage grande vitesse est exposé et utilisé pour identifier les sources d’échauffement. Ce modèle
se base sur un projet de banc d’essai qui servira dans le futur à valider les choix technologiques
faits sur l’étage grande vitesse.
L’arbre d’entrée du réducteur pour être guidé en rotation, est porté par deux roulements à
billes à contact oblique. A grande vitesse ces roulements produisent beaucoup de pertes de
puissance. Leur modélisation thermomécanique est donc importante à l’échelle du réducteur
complet.
Le chapitre trois s’attache à la modélisation de ces roulements autant sur le plan de la
prédiction des pertes de puissance que sur le plan de la modélisation des transferts thermiques
aillant cours en leur sein.
Le bain d’huile a pour rôle de lubrifier les différents éléments du réducteur, cependant il sert
aussi à refroidir les mobiles. A grande vitesse les pertes par barbotage deviennent importantes,
pour les limiter le niveau d’huile doit être faible. Il y a donc une compétition entre avoir beaucoup
d’huile pour bien refroidir, et en avoir peu pour limiter les pertes par barbotage. A haute vitesse
et haute température l’huile se charge d’air. Ce phénomène appelé aération augmente les pertes
par barbotage. Ces deux problématiques font que la lubrification par barbotage n’est pas
recommandée pour un réducteur haute vitesse. Le chapitre quatre porte sur l’analyse et
l’optimisation des pertes par barbotage. A l’aide d’un banc d’essai une solution simple et peu
coûteuse est testée. Cette solution permet de diminuer les pertes par barbotage ainsi que de
limiter le phénomène d’aération tout en conservant la fonction thermique du bain d’huile. Elle
permet d’envisager pour un réducteur haute vitesse l’utilisation du barbotage comme mode de
lubrification.
15
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Chapitre 1
Etude bibliographique
1.2 Le lubrifiant
Les lubrifiants pour réducteur mécanique sont composés d’une huile de base à laquelle sont
ajoutés des additifs. Il existe deux grands types d’huile de base [5]. Le premier est constitué par les
huiles issues de la distillation du pétrole brut. Bien que constituées de chaines organiques, elles
sont appelées huiles minérales à cause de leur origine. Ces huiles sont composées d’un mélange
de chaines d’hydrocarbures de structures et de masses molaires différentes.
Le second type d’huile de base est constitué par les huiles synthétiques. Elles sont obtenues
par synthèse chimique. Il est alors possible de créer des huiles par des procédés variés et même de
substituer le carbone couramment utilisé par des silicones ou encore d’utiliser du fluor dans les
formulations. Pour les transmissions mécaniques, les huiles synthétiques sont le plus souvent des
Polyalphaoléfines hydrogénées (dites PAO), des polyglycoles (dits PAG), ou des huiles à base
d’ester. Enfin il existe des huiles de base dites semi-synthétiques qui sont un mélange d’huiles des
deux types précédemment exposés.
Pour améliorer les propriétés de l’huile de base de nombreux additifs sont rajoutés [6], [7].
Leurs actions sont diverses avec par exemple des effets d’amélioration d’index de viscosité, de
protection de surface ou encore d’anti-mousse. Tout ceci fait que les lubrifiants sont très
complexes du point de vue physico-chimique.
L’huile a la double fonction de lubrifier le système et de le refroidir. Même si elle est très
complexe un certain nombre de caractéristiques sont classiquement utilisées pour décrire ces
deux fonctions :
18
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
- La masse volumique 𝜌 qui est pour les fluides utilisés dans les réducteurs mécaniques
comprise entre 800 et 1000 kg/m3
- La viscosité cinématique 𝜈 [m²/s]
- La viscosité dynamique 𝜇 = 𝜈. 𝜌 [Pa.s]
- Le coefficient de piezo-viscosité de Barus 𝛼𝑏 [1/Pa] qui caractérise l’évolution de la
viscosité dynamique en fonction de la pression 𝑃0 [Pa] dans le contact et de la
viscosité 𝜇0 à son entrée par l’équation suivante [8] :
𝜇 = 𝜇0 𝑒 𝛼𝑏 𝑃0 1.1.
- La conductivité thermique k [W/m.K]
- La capacité thermique massique 𝐶𝑝 [J/kg.K]
Parmi ces propriétés certaines évoluent en fonction de la température. La suite du paragraphe
s’attache à exposer les relations couramment utilisées pour caractériser en fonction de la
température les huiles pour la lubrification de réducteurs mécaniques [9], qu’elles soient minérales
ou synthétiques.
La masse volumique est en général mesurée à 15°C. Elle évolue avec la température par
l’intermédiaire du coefficient de dilatation thermique. Celui-ci varie peu d’une huile à l’autre, une
bonne approximation de ce coefficient est 0,00065. La masse volumique est alors obtenue par [9]:
−0,03
𝑘= . (𝑇 − 273) + 0,14 1.4.
400
La capacité thermique 𝐶𝑝 est dans le même cas. La relation suivante permet de l’estimer [9] :
19
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
𝑄=∫ 𝑑𝑇 ⃗⃗⃗⃗2 (𝑀)
⃗ 1 (𝑀). 𝑉 1 1.6.
𝑀∈𝐿(𝜏) 2
Avec
𝑄 = 𝑃. 𝑓. 𝐻𝑣 1.7.
Avec
- 𝑃 la puissance transmise [W]
- 𝑓 le coefficient de frottement moyen le long de la ligne d’action
- 𝐻𝑉 un facteur géométrique
Il ressort de ces études que les pertes par frottement à la denture dépendent uniquement d’un
coefficient de frottement moyen et d’un facteur géométrique 𝐻𝑉 . Chaque étude a fait émerger
une formulation différente pour le facteur 𝐻𝑉 . Ainsi pour Buckingham l’expression de 𝐻𝑉 est :
𝜋 1
𝐻𝑣 𝐵𝑢𝑐𝑘𝑖𝑛𝑔ℎ𝑎𝑚 = (1 + 𝑢). . . 𝜀𝛼 . (2. 𝜅02 − 2𝜅0 + 1) 1.8.
𝑍1 𝑐𝑜𝑠 𝛽𝑏
Pour Niemann l’expression devient :
𝜋 1 1
𝐻𝑣 𝑁𝑖𝑒𝑚𝑎𝑛𝑛 = (1 + 𝑢). . . 𝜀𝛼 . ( − 1 + (2. 𝜅02 − 2𝜅0 + 1)𝜀.𝛼 ) 1.9.
𝑍1 𝑐𝑜𝑠 𝛽𝑏 𝜀𝛼
Avec
𝑍
- 𝑢 le rapport de réduction de la transmission = 𝑍1 𝑎𝑣𝑒𝑐 𝑍1 ≤ 𝑍2
2
- 𝑍1 le nombre de dents du pignon
- 𝑍2 le nombre de dents de la roue
- 𝛽𝑏 l’angle d’hélice de base
- 𝜀𝛼 le rapport de conduite apparent
- 𝜅0 la position adimensionnel du point de roulement sans glissement le long de la ligne
d’action
20
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
1 1 𝑔𝑓2 + 𝑔𝑎2
𝐻𝑣 𝐻𝑒𝑛𝑟𝑖𝑜𝑡 = ( + ). 1.10.
𝑅𝑝1 𝑅𝑝2 2. (𝑔𝑓 + 𝑔𝑎 ). 𝑐𝑜𝑠 𝛼𝑡
Avec
- 𝑅𝑝1 et 𝑅𝑝2 les rayons primitifs du pignon et de la roue [m]
- 𝑔𝑓 et 𝑔𝑎 les longueurs d’approche et de retrait [m]
- 𝛼𝑡 l’angle de pression apparent [rad]
Aucune des trois formules ne prend en compte les corrections de forme des flancs de
denture. Ce type de correction consiste à un faible enlèvement de matière sur la surface des dents,
comme par exemple sur la Figure 1-1. Cet usinage permet d’optimiser de nombreuses
caractéristiques de fonctionnement des dentures [14]. Il est notamment possible de diminuer les
pertes de puissance en enlevant de la matière en tête et en pied de denture [15]. Dans cette zone
la vitesse de glissement relative entre les flancs en contact est maximum. L’enlèvement de matière
permet de diminuer la pression de contact en tête et pied de dent et ainsi de faire diminuer le
produit de la pression par la vitesse.
𝜋 1
𝐻𝑣 𝑉𝑒𝑙𝑒𝑥 = (1 + 𝑢). . .𝜀 .Λ 1.11.
𝑍1 cos 𝛽𝑏 𝛼
Avec :
2. 𝛤𝑝
[1 − 2. 𝜅0 + 2. 𝜅02 − 𝛶. 𝛤𝑝 . [1 − 3 ]]
𝛬=
𝜋 2 2. 𝛤𝑝
𝑡𝑎𝑛𝛼𝑝 . (2. 𝜅0 − 1) − 𝜀𝛼 . (2𝜅 − 2𝜅0 + 1 − 𝛶. 𝛤𝑝 . (1 − )) 1.12.
𝑍1 0 3
1 − 𝛶. 𝛤𝑝 − 𝑓
𝑐𝑜𝑠𝛽𝑏
[ [ ]]
21
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Avec comme représenté sur la Figure 1-2 :
- 𝛶 profondeur de correction de profil adimensionnée par rapport à Δ𝑚
- 𝛤𝑝 longueur de correction de profil adimensionnée par rapport à la longueur du profil
actif
- 𝑝𝑏𝑡 pas de base apparent [m]
- Δ𝑚 valeur moyenne d’approche normale par rapport à la position de corps rigide [m]
- 𝑓 coefficient de frottement moyen
Γ. 𝜀𝛼 . 𝑝𝑏𝑡
Υ. Δ𝑚
𝜀𝛼 . 𝑝𝑏𝑡
1 2 3
Figure 1-3 Courbe de Stribeck [17] évolution du coefficient de frottement en fonction du rapport hauteur du film
d’huile sur la rugosité moyenne quadratique
22
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
3. Le régime de lubrification en film complet. Il n’y a aucun contact entre les surfaces et
le coefficient de frottement est très dépendant de la viscosité de l’huile
1. Le régime de lubrification limite. La charge est portée par les aspérités. Le coefficient
de frottement est très dépendant des additifs et des caractéristiques de surfaces.
2. Le régime de lubrification mixte. La situation est un mélange des deux précédentes.
C’est le régime de lubrification fréquemment rencontré dans les engrenages et les
roulements.
La caractérisation du coefficient de frottement est très complexe. L’ensemble des
phénomènes physico-chimique mis en jeu nécessite de passer par des études expérimentales.
Dans le cas des engrenages plusieurs expressions du coefficient de frottement moyen ont été
proposées par Kuzmin [18], Misharin [19], O’Donoghue et Cameron [20] ou encore Kelley et
Lemanski [21]. Benedict-Kelley [22] propose un formulation qui semble la plus complète par le
nombre de paramètres pris en compte et par le nombre d’essais ayant servi à son élaboration. Il
propose :
291205,8. 10−7
𝑓 = 0,0127. 𝑙𝑜𝑔10 1.13.
𝜌. 𝜈. 𝑉𝑔 . 𝑈 2
[ 𝐹𝑛𝑢 ]
Avec
- 𝜌 la masse volumique du lubrifiant [kg/m3]
- 𝜈 la viscosité cinématique du lubrifiant [m2/s]
- 𝑉𝑔 la vitesse de glissement (différence des vitesses de roulement) [m/s]
- 𝑈 la somme des vitesses de roulement [m/s]
- 𝐹𝑛𝑢 l’effort normal unitaire [N/m]
Cette loi est proche des expérimentations uniquement pour les taux de glissement élevés
(entre 20 et 50%) [15]. Dans un souci de proposer une loi fonctionnant sur une plage de
conditions opératoires plus large, la norme ISO 14179-2 [23] propose une autre formulation :
0,2
𝐹𝑛 ⁄𝐵 0,25
𝑓 = 0,048. [ ] . (𝜇. 103 )−0,05 . 𝑅𝑎𝑟 . 𝑋𝐿 1.14.
𝑟𝑐 . 𝑈. 106
Avec
- 𝑋𝐿 un coefficient dépendant du type d’huile utilisé
- 𝑟𝑐 le rayon de courbure équivalent au niveau du diamètre primitif [m]
- 𝐵 la largeur de denture en contact [m]
- 𝜇 la viscosité dynamique du lubrifiant [Pa.s]
- 𝑅𝑎𝑟 la moyenne arithmétique des rugosités des flancs actifs du pignon et de la roue [μm]
- 𝐹𝑛 l’effort normal [N]
Une autre approche possible est de décrire le coefficient de frottement instantané. C’est
l’approche qu’a utilisé Diab et al. au sein du LaMCoS [24]. En se basant sur les théories de
Greenwood-Tripp [25] et de Mikic [26] une loi de frottement a été construite. Cette loi est établie
pour un couple lubrifiant/matière à partir de nombreux essais expérimentaux réalisés sur une
machine à galets à haute performance [27] du LaMCoS.
23
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Figure 1-4 Comparaison des coefficients de frottement calculés avec les trois lois de frottement en fonction de la
phase de vol considérée [28]
Figure 1-5 Comparaison entre quelques lois de frottement de la littérature en fonction du taux de glissement SRR
[24]
A travers cette étude bibliographique, il ressort qu’un grand nombre de formulations pour
estimer le coefficient de frottement moyen sont disponibles. Par contre aucune ne semble
s’imposer : chacune fonctionnant dans le domaine de validité ayant servi à son élaboration.
La Figure 1-4 montre une comparaison réalisée par Durand [28] pour une transmission de
puissance aéronautique. Il compare deux lois de frottement avec des essais réalisés par Diab pour
différentes phases de vol qui représentent des conditions de couple et de vitesse différentes. La
norme ISO sous-estime le coefficient de frottement. Pour certaines conditions utilisées la
formule de Benedict-Kelley va même jusqu’à prédire des coefficients négatifs (l’expression dans
la fonction 𝑙𝑜𝑔10 devenant inférieure à 1). Les lois utilisées ne sont pas capables de reproduire la
mesure, elles sous estiment le coefficient de frottement. La Figure 1-5 expose une autre
comparaison réalisée par Diab [24], cette fois en faisant varier le taux de glissement SRR défini
par :
24
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
𝑈1 − 𝑈2
𝑆𝑅𝑅 = 1.15.
𝑈1 + 𝑈2
Avec
- 𝑈1 vitesse du contact par rapport au solide 1 [m/s]
- 𝑈2 vitesse du contact par rapport au solide 2 [m/s]
Cette fois les lois qui sont comparées sur estiment le coefficient de frottement. Il ressort que
les lois donnant un coefficient de frottement moyen peuvent être prises en défaut. A l’aide de
nombreux essais il est bien entendu possible d’élaborer une loi de frottement prenant en compte
beaucoup plus de paramètres [24].
Comme ce travail de thèse se situe dans une phase d’étude amont de la conception du
réducteur, il n’est pas envisageable de mener des essais pour caractériser le coefficient de
frottement. Celui-ci sera alors imposé à une valeur moyenne constante sur tout l’engrènement.
Dans ces conditions les pertes aux dentures ne dépendent que de la puissance transmise et de la
géométrie des dentures (voir formule 1.7).
𝑀𝑡𝑜𝑡𝑎𝑙 = 𝑀0 + 𝑀1 1.16.
Où
25
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
𝑀0 = 103 𝑓0 (𝜈𝑁)2/3 𝑑𝑚
3
1.17.
𝑀0 = 16. 𝑓0 . 𝑑𝑚 3 1.18.
Et
𝑀1 = 𝑓1 𝑃1𝑎′ 𝑑𝑚
𝑏′
1.19.
Avec
- 𝑓0 un facteur fonction du type de roulement et du type de lubrification
- 𝜈 viscosité cinématique du lubrifiant [m²/s]
- 𝑓1 un facteur fonction de la charge statique équivalente et statique de base, et du type de
roulement
- 𝑃1 est une combinaison des charges axiale et radiale [N]
- a’ et b’ sont des facteurs de correction fonction du type de roulement
- 𝑑𝑚 diamètre moyen du roulement [m]
- 𝑁 vitesse de rotation du roulement [tr/min]
Ce modèle est celui utilisé par les normes ISO14179 partie 1 et 2 permettant d’estimer les
pertes dans les transmissions mécaniques à engrenage [23], [31]. Le couple 𝑀0 dépend de la
vitesse à la puissance 2/3. La formule montre qu’à grande vitesse les pertes indépendantes de la
charge croissent beaucoup plus vite que les pertes directement liées à la charge. A grande vitesse
le modèle montre qu’il est intéressant de diminuer la viscosité du lubrifiant pour diminuer les
pertes de puissance.
Un autre modèle est couramment utilisé. Il a été proposé par SKF [32] en 2004. A l’inverse
du modèle de Palmgren celui-ci a été développé sur la base de l’utilisation d’un outil numérique.
La méthode se base sur l’addition de trois moments, l’un généré par la résistance au roulement,
l’autre par le frottement de glissement et un dernier représentant les pertes par trainée liées au
mode de lubrification :
26
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
27
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
1
𝜑𝑖𝑠ℎ = 1.27.
1 + 1.84 ∗ 10−9 ((𝑁𝑑𝑚 . 103 )1.28 𝜈 0.64 )
2. La vitesse des billes tend à éjecter le lubrifiant qui se retrouve en moindre quantité sur le
chemin de roulement. SKF prend en compte la possible sous-alimentation des contacts
résultante via le facteur 𝜑𝑟𝑠 .
1
𝜑𝑟𝑠 = 3 1.28.
𝑒 𝐾𝑟𝑠 𝜈𝑁√𝐾𝑧 ⁄(2(𝐷−𝑑).10 )
Avec
- 𝐾𝑟𝑠 constante fonction du type de lubrification
- 𝐾𝑧 constante géométrique fonction du type de roulement
- 𝐷 diamètre extérieur [mm]
- 𝑑 diamètre intérieur [mm]
La vitesse de rotation a une large influence sur les coefficients 𝜑. A froid le produit 𝜈𝑁 est
grand, ces coefficients ont donc un grand impact au démarrage. A chaud la viscosité diminue
fortement et (𝜑𝑖𝑠ℎ , 𝜑𝑟𝑠 ) ont moins d’impact.
Dans ce modèle la vitesse de rotation intervient à de nombreux endroits. La dépendance du
couple à la vitesse de rotation par une fonction de puissance 2/3 se retrouve dans le couple 𝑀𝑟𝑟 ,
la puissance passant à 0,6. De plus pour certains types de roulement 𝐺𝑟𝑟 et 𝐺𝑠𝑙 dépendent eux
aussi de la vitesse au carré, c’est notamment le cas des formulations pour les roulements à
contacts obliques.
Fernandes et al. [33] montrent que le modèle de SKF donne de meilleurs résultats en
comparaison avec le modèle de Harris pour deux butées, l’une à billes et l’autre à rouleaux. Les
vitesses des essais sont inférieures à 3000tr/min. Il est possible d’améliorer encore la prédiction
des pertes par un recalage de 𝜇𝑏𝑙 et de 𝜇𝐸𝐻𝐿 par des essais. Cependant pour une application haute
vitesse il est difficile de sélectionner un modèle en particulier.
28
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
brassage du bain d’huile et le problème que pose ce mode de lubrification à haute vitesse. En
effet à grande vitesse la lubrification est habituellement faite par injection [43] pour des
problématiques de rendement.
𝜌. ω2 . 𝑆𝑚 . 𝑅𝑝 3 . 𝐶𝑚
𝐶𝑏𝑎𝑟𝑏 = 1.29.
2
Avec
- 𝜌 la masse volumique de l’huile [kg/m3]
- ω la vitesse de rotation du pignon [rad/s]
- 𝑆𝑚 la surface immergée du pignon [m²]
- 𝑅𝑝 le rayon primitif du pignon [m]
- 𝐶𝑚 un couple adimensionnel dont la valeur est donnée par des formulations semi-
empiriques
29
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Pignon
Rp Bain
d’huile
θ
h
Avec
2. 𝑍. 𝜃. 𝐻𝑑𝑒𝑛𝑡 . 𝐵
𝑆𝑚𝑑 = 𝐷𝑝 . 𝑏. 𝜃 + 1.31.
𝜋. cos 𝛼 . cos 𝛽
Avec
- 𝐷𝑝 le diamètre primitif du pignon [m]
- 𝐵 la largeur du pignon [m]
- 𝑍 le nombre de dents du pignon
- 𝐻𝑑𝑒𝑛𝑡 la hauteur des dents du pignon [m]
- 𝛼 l’angle de pression normale [rad]
- 𝛽 l’angle d’hélice [rad]
Comme dans l’étude de Daily et Nece le couple adimensionnel de trainée est dépendant du
régime d’écoulement. Changenet a trouvé quatre régimes d’écoulement dépendant de deux
paramètres [54]. Le premier est assimilable à une accélération centrifuge :
1⁄
3
𝛾 = ω2 . (𝑅𝑝 . 𝐵. 𝑚) 1.32.
Avec
- 𝑚 le module normal du pignon [m]
ω. 𝑅𝑝 . 𝐵
𝑅𝑒𝑐 = 1.33.
𝜈
Enfin dans les expressions des couples adimensionnels de trainée, un dernier nombre
adimensionnel apparaît : le nombre de Froude :
30
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
ω2 . 𝑅𝑝
𝐹𝑟 = 1.34.
𝑔
Le couple adimensionnel de trainée de barbotage s’exprime alors :
31
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
vitesses et à haute température. Il relève aussi que l’aération est sensible à plusieurs paramètres,
elle augmente avec l’accroissement de :
La vitesse
L’immersion
L’espace inter-dent
La température
Figure 1-7 Evolution de l'aération en fonction de la température pour différentes vitesses de rotation [58]
Leprince met en évidence le rôle de l’aération sur les pertes de puissance. La Figure 1-8
montre le comportement typique des pertes face à l’aération. A basse température, ou à basse
vitesse, l’aération a un comportement linéaire face à la température. Avec l’échauffement le bain
d’huile se dilate et se charge en air, l’immersion augmente faiblement avec la température. Le
modèle exposé dans le paragraphe précédent (trait noir sur la Figure 1-8) est capable de prédire
les pertes de puissance sur ces plages de dilatation (et plus généralement tant que l’aération est
inférieure à 10%). A partir d’une certaine température dépendante des conditions d’essai et de
l’huile, l’aération s’envole. Les pertes de puissance augmentent sans que le modèle ne puisse
retranscrire cette tendance.
32
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
Figure 1-8 Mise en évidence de la relation entre aération du lubrifiant et perte par barbotage [57]
Peu d’études existent sur la problématique d’aération du lubrifiant dans les transmissions
mécaniques. Des études portant sur l’aération des huiles et son impact dans la lubrification des
paliers ont été conduites [59]–[61]. Dans le cas de la lubrification par barbotage les études se sont
concentrées sur le moussage du lubrifiant. Certaines pointent tout de même que les deux
phénomènes pourraient être liés au travers de l’action d’additifs antimousses [62]–[64]. Seul
Leprince s’est attaché à caractériser les pertes de puissance engendrées par l’aération dans le cas
du barbotage [58].
Leprince fait l’hypothèse que les pertes supplémentaires dues à l’aération sont liées au
fractionnement des bulles d’air du bain par la rotation du mobile. Il élabore alors une modèle dit
heuristique basé sur l’hypothèse que la perte liée à l’aération est due à l’énergie nécessaire pour
fractionner les bulles d’air piégées dans le bain d’huile. Il calcule les pertes supplémentaires dues à
l’aération par :
𝑉𝑝𝑖𝑔𝑛𝑜𝑛 𝐴% − 𝐴%0
𝑃𝑎é𝑟𝑎𝑡𝑖𝑜𝑛 = 𝛼0 𝜎𝑠 𝜔𝐻(𝐴% − 𝐴%0 ) 1.39.
𝑅𝑏 100
Avec
- 𝛼0 coefficient fixé empiriquement
- 𝜎𝑠 tension de surface du lubrifiant [N/m]
- 𝐴%0 l’aération à partir de laquelle le comportement des pertes change
- 𝐻 fonction de Heaviside
- 𝜔 La vitesse de rotation [rad/s]
- 𝑉𝑝𝑖𝑔𝑛𝑜𝑛 volume immergé du pignon [m3]
- 𝑅𝑏 rayon moyen des bulles d’air [m]
Il développe aussi une méthodologie pour calculer 𝐴%0 et 𝑅𝑏 . Son modèle ne permet
cependant pas de prédire 𝐴%, la prédiction des pertes est donc soumise soit à une mesure
d’aération soit à une hypothèse sur celle-ci. Il montre cependant qu’à haute vitesse (vitesse
périphérique > 40m/s) et haute température (>80°C) l’aération dépasse 10% et devient un
phénomène majeur. L’aération du bain d’huile est un problème pour un réducteur mécanique à
haute vitesse et devra par conséquent être prise en compte.
33
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
1.6 Les pertes par frottement aux joints
Pour maintenir le lubrifiant dans le système il est nécessaire d’étanchéifier l’arbre d’entrée et
l’arbre de sortie du réducteur. La solution la plus couramment utilisée en mécanique automobile
est le joint à lèvres.
Ce type de joint est composé d’une lèvre maintenue plaquée sur l’arbre par un ressort de
forme torique. La chaleur produite par le frottement est le produit de la vitesse périphérique avec
la force de frottement tangentielle :
𝑑
𝑄 = 𝐹𝑅 . 𝑓. . 𝜔 1.40.
2
Avec :
- 𝐹𝑅 la force radiale appliquée par le joint sur l’arbre [N]
- 𝑓 le coefficient de frottement de la lèvre sur l’arbre
Cette équation n’est pas utilisable directement. En effet les paramètres 𝐹𝑅 et 𝜇 sont influencés
par de nombreux éléments. Il est possible d’en citer quelques-uns, comme la conception du joint,
les matériaux utilisés, les conditions d’utilisation etc.…
Pour approximer les pertes par frottement plusieurs lois existent dans les normes ISO
traitant des pertes dans les transmissions mécaniques. Les pertes aux joints calculées par la norme
ISO 14179-1 [31] sont représentées sur la Figure 1-9. Cette norme suggère que le couple de
frottement ne dépend que du diamètre de l’arbre, de façon linéaire, et du matériau du joint.
34
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
35
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Figure 1-11 Comparaison entre le rendement d’une boîte de vitesses mesuré expérimentalement, estimé à l’aide d’un
modèle isotherme et calculé à partir d’un modèle basé sur la méthode nodale [42]
Les pertes de puissance, notamment aux roulements et par barbotage, augmentent fortement
avec la montée en vitesse. Dans le cas d’un réducteur mécanique haute vitesse beaucoup de
puissance se retrouve dissipée dans un faible volume. L’échauffement résultant peut être critique
pour le composant. Il est donc important d’estimer les températures internes de fonctionnement
de la transmission.
Figure 1-12 Exemple de modélisation par éléments finis de l'échauffement d'un palier à roulement [69]
36
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
Il est aussi possible d’utiliser l’analyse par éléments finis. Cette méthode permet de décrire
finement le gradient de température dans les pièces d’un réducteur mécanique [69], [70]. Un
exemple d’analyse thermique par la méthode des éléments finis est donné sur la Figure 1-12. Avec
cette méthode il est aussi possible de modéliser l’échange thermique entre les solides et l’huile par
une approche dite « computational fluid dynamics » (CFD). Les calculs sont complexes à mettre
en place et demandent une puissance de calcul importante.
∆𝑇
𝑄= 1.43.
𝑅𝑡ℎ
L’équation ci-dessus permet d’évaluer les échanges de chaleur entre les nœuds. Une fois le
réseau thermique construit et les échanges entre les nœuds retranscrits, le premier principe de la
thermodynamique est appliqué sur chaque élément du système. Il est alors possible d’écrire autant
d’équations de conservation de l’énergie qu’il y a d’éléments dans le réseau thermique :
𝑛
𝑑𝑇𝑖 1 (𝑇𝑖 − 𝑇𝑗 )
= [𝑄 − ∑ ] 1.44.
𝑑𝑡 𝑚𝑖 𝐶𝑖 𝑖 𝑅𝑡ℎ,𝑖𝑗
𝑗=1 𝑒𝑡 𝑗≠𝑖
Avec :
- 𝑚𝑖 représente la masse de l’élément i [kg]
- 𝐶𝑖 chaleur massique de l’élément i [J/kg.K]
𝑑𝑇𝑖
- dérivée de la température par rapport au temps [K/s]
𝑑𝑡
- 𝑄𝑖 valeur de la puissance thermique injectée au nœud i [W]
- 𝑇𝑖 (resp. 𝑇𝑗 ) la température de l’élément i (resp. j) [K]
- 𝑅𝑡ℎ,𝑖𝑗 la résistance thermique entre les éléments i et j [K/W]
- 𝑛 le nombre de nœuds du réseau thermique
37
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Les conditions aux limites des réseaux thermiques sont des nœuds où la température est
imposée. Par exemple l’air ambiant ou le sol ont une température définie, ils sont des puits de
𝑑𝑇𝑖
chaleur. Pour ces nœuds 𝑚𝑖 𝐶𝑖 ∗ 𝑑𝑡
= 0 c’est-à-dire que leur énergie interne n’évolue pas. Qi
représente alors un flux de chaleur absorbé par le nœud limite.
b) Résolution
C. Changenet introduit la notation matricielle suivante [73] :
Il est alors possible d’extraire l’équation différentielle matricielle décrivant l’évolution des
températures du système :
𝐶𝑐 𝑇′𝑖𝑛 = 𝑄𝑐 − 𝑆3 𝑇𝑐 − 𝑆4 𝑇𝑖𝑛 1.50.
38
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
Sachant que pour les nœuds limites 𝐶𝑖𝑛 𝑇′𝑐 = 0 il est possible d’utiliser à tout moment
l’équation 1.47. pour déterminer la chaleur absorbée par les nœuds limites.
𝑑𝑇
𝑄 = −𝑘 1.51.
𝑑𝑥
L’application de cette loi permet de calculer la résistance thermique. Par exemple si l’on
considère un flux de chaleur à travers une paroi plane d’épaisseur e, de surface S, et de
conductivité k soumise à une différence de température il est possible d’écrire :
𝑘𝑆
𝑄= ∆𝑇 1.52.
𝑒
La résistance de conduction dans ce cas est alors par identification :
𝑒
𝑅𝑡ℎ = 1.53.
𝑘𝑆
b) La seconde, la résistance de convection, traduit le transfert de chaleur ayant lieux entre un
solide et un fluide en mouvement. L’évaluation du flux de chaleur peut alors se faire par
l’intermédiaire d’un coefficient de convection ℎ𝑐𝑜𝑛𝑣 :
1
𝑅𝑡ℎ = 1.55.
ℎ𝑐𝑜𝑛𝑣 𝑆
Lors de l’échange de chaleur il y a compétition entre la conduction dans le fluide et la
convection due au mouvement du fluide. Le nombre de Nusselt 𝑁𝑢 quantifie le phénomène :
ℎ𝑐𝑜𝑛𝑣 𝐿
𝑁𝑢 = 1.56.
𝑘
Lorsque qu’il est proche de 1, c’est-à-dire pour un fluide stagnant, les effets de conduction
sont prédominants. S’il y a mouvement du fluide le nombre de Nusselt devient plus important
traduisant un échange convectif élevé. L’application du théorème de Vaschy-Buckingham permet
d’exprimer le nombre de Nusselt comme une fonction de deux nombres adimensionnels [72].
Dans le cas de la convection forcée :
39
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Dans le cas de la convection naturelle :
𝑄 = 𝑆𝜎𝜀𝑇 4 1.59.
Avec
- 𝜀 l’émissivité du matériau, comprise entre 0 et 1
- 𝜎 la constante de Stefan-Boltzmann = 5,67. 10−8 [W/m².K4]
Une partie de la puissance rayonnée est absorbée par les solides environnants. Il y a donc
entre chaque objet un équilibre entre l’émission et la réception. Ainsi le rayonnement thermique
s’échangeant entre deux corps, séparés par un milieu non absorbant, peut s’exprimer de la façon
suivante :
𝜎. (𝑇1 4 − 𝑇2 4 )
𝑄=
1 − 𝜀1 1 1 − 𝜀2 1.60.
𝜀1 𝑆1 + 𝑆1 𝐹12 + 𝜀2 𝑆2
Avec
- 𝐹12 facteur de forme de 1 vers 2 dépendant de la géométrie du système
- 𝑆1 surface rayonnante du solide 1 [m²]
- 𝑆2 surface rayonnante du solide 2 [m²]
1 − 𝜀1 1 1 − 𝜀2
+ + 1
𝜀 𝑆 𝑆1 𝐹12 𝜀2 𝑆2
R th = 1 1 . 1.61.
σ (T1 2 + T2 2 )(T1 + T2 )
1.8 Conclusion
Pour étudier le rendement d’un réducteur mécanique lubrifié par barbotage, quatre sources
de perte de puissance ont été identifiées : (i) frottements à la denture, (ii) pertes dans les
roulements, (iii) barbotage et (iv) frottement des joints.
Dans le cadre de la présente étude, ciblant les réducteurs à haute vitesse pour des applications
automobiles :
Sur l’arbre d’entrée il semble critique d’utiliser un joint à lèvre standard, le
concepteur se dirigera vers d’autres solutions provoquant moins de pertes de puissance. Sur
40
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude bibliographique
l’arbre de sortie l’étanchéité sera assurée par des joints à lèvre. Dans l’application visée cet arbre
ne tourne pas vite et les pertes dues aux joints seront faibles.
Trois autres sources de pertes de puissance sont prépondérantes. Les pertes par
frottement des dentures sont inévitables. L’étude bibliographique relève qu’elles dépendent de la
puissance transmise, de la géométrie des dentures et du coefficient de frottement. Les grandes
vitesses vont agir sur le coefficient de frottement moyen. Plusieurs formules sont disponibles
pour l’estimation de ce coefficient, il est aussi possible par de nombreux essais de le caractériser.
Sans essais il sera fixé dans les modèles à une valeur réaliste en vue des technologies automobiles
utilisées.
Deux phénomènes risquent d’être prépondérants avec l’augmentation de la vitesse
d’entrée. Les pertes indépendantes de la charge dans les roulements sont très sensibles à la vitesse
de rotation. Les pertes par barbotage augmentent elles aussi avec la vitesse. Le phénomène est
aggravé par l’aération du bain d’huile, qui est provoquée par des mobiles tournant à grand vitesse.
Le réducteur mécanique sera conçu pour être le plus compact possible. De ce fait
les pertes de puissance seront dissipées sous forme de chaleur dans un faible volume. Il risque
donc d’y avoir des échauffements problématiques. De plus les pertes de puissance dépendent
fortement de la température de fonctionnement. La méthode des réseaux thermique permet
d’identifier les échauffements et d’affiner la prédiction du rendement en couplant pertes de
puissance et température.
41
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Chapitre 2
Positionnement du problème
44
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
Le Tableau 2-1 expose aussi des contraintes liées à la puissance transmise ainsi qu’au couple
en entrée du réducteur. Il est envisagé d’utiliser un moteur d’une puissance nominale de 50 kW,
celle-ci pouvant monter à 80 kW sur un faible temps (inférieur à 30s). Le dimensionnement en
termes de résistance mécanique est réalisé avec le couple maximum de 55Nm. La Figure 2-1
montre le couple que peut fournir le moteur en fonction de sa vitesse de rotation.
60
50
40
Couple [N-m]
30
20
10
0
0 10000 20000 30000 40000
Vitesse [tr/min]
Nominal Transitoire
45
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Tableau 2-2 Profil de mission du réducteur
Vitesse d’entrée Couple moteur Puissance [kW] Heures
[tr/min] [N.m]
Démarrage 5000 55 28 600
Villes 50 km/h 16000 15 21 1688
Routes 90 km/h 28000 17,1 50 1071
Autoroute 130 km/h 42000 11,4 50 714
Total 4073
46
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
47
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Tableau 2-5 Caractéristiques des dentures de l'architecture 2 trains
Engrène Nombre Module Largeur Diamètre Angle de Angle
ment de dents normal [mm] primitif pression d'hélice
[mm] [mm] [°] [°]
GV Pignon 19 22 1,5 30,2 20 20
Roue 108 20 171,8
PV Pignon 14 3 36 44,7 20 20
Roue 80 31 255,3
Engrènement Rapport de Rapport de Matériau Classe de Rugosité
conduite recouvrement qualité arithmétique
[mm]
GV 1.496 1.452 16MnCr5 7 0.6
cémenté
PV 1.399 1.125 16MnCr5 7 0.6
cémenté
Les réducteurs deux trains et trois trains ont respectivement des rapports de réduction de
32,48 et de 32,37. Pour arriver à ce rapport de réduction les roues du réducteur deux trains sont
plus grandes. Le Tableau 2-7 compare le poids du réducteur et le diamètre de la plus grande pièce
mobile c’est à dire de la roue petite vitesse (PV). La roue de sortie est plus grande sur la solution
2 trains, pour l’accueillir le carter sera plus encombrant. De plus en considérant en première
approximation que les engrenages sont des cylindres pleins, la masse totale hors carter du
réducteur deux trains est plus élevée que celle du réducteur trois trains.
Le carter n’a pas fait l’objet d’un dimensionnement, cependant dans les deux cas son poids
serait comparable. Le carter du réducteur 3 trains comporte deux logements de roulements
supplémentaires mais la partie accueillant la roue petite vitesse sera plus petite.
L’étanchéité du carter sera réalisée avec deux joints à lèvre sur l’arbre petite vitesse, et avec
un joint labyrinthe sur l’arbre d’entrée.
48
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
2 Trains
3 Trains
Figure 2-2 Exemple de configuration à iso entraxe entrée sortie, en bleu l’architecture 3 trains, en gris l’architecture 2
trains
Les caractéristiques du lubrifiant envisagé pour le réducteur sont exposées sur le Tableau 2-8.
Cette huile est équivalente à un lubrifiant de boite de vitesses automobile.
Tableau 2-8 Caractéristiques du lubrifiant
Viscosité à 40°C [cSt] Viscosité à 100°C [cSt] Masse volumique à 15°C [kg/m3]
45,1 7,7 885
49
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
- L’application des normes ISO 14179 [23], [31] donne un coefficient de frottement
denture appartenant à l’intervalle [0,03 :0,04] pour les conditions envisagées. Le
paragraphe 1.3.2 a montré que les formules permettant de calculer ce coefficient peuvent
le sous-estimer. Sans essais de caractérisation du coefficient de frottement denture celui-
ci est pris égal à une valeur de 0,05.
- Les pertes par barbotage sont calculées avec le modèle de Changenet et .al [54]
- Quand il est spécifié une aération du bain d’huile, les pertes par barbotage
supplémentaires dues à l’aération sont calculées avec la méthode de Leprince et .al [74]
voir paragraphe 1.5.2.
- Le niveau d’huile est fixé pour que 25% du rayon de la roue de l’étage GV soit immergée
dans l’huile. Ce niveau est un compromis, à grande vitesse l’immersion doit être faible
pour limiter les pertes par barbotage [43] cependant un niveau d’huile trop petit peut
diminuer fortement les projections d’huile [75] et empêcher la bonne lubrification des
roulements.
- Les arbres sont dans un plan parallèle au plan formé par le bain d’huile à l’arrêt.
- Les pertes des joints sont négligées, elles représentent sur l’arbre de sortie 9W et 6W
pour les conditions autoroute et route.
Le Tableau 2-9 donne les rendements calculés pour les deux architectures. Comme identifié
dans le chapitre 1 l’augmentation de la vitesse à iso puissance entraine la diminution du
rendement. Les deux architectures ont des rendements comparables. L’architecture à trois trains
est désavantagée par deux roulements et un engrènement supplémentaires.
50
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
3 Trains Autoroute
2 Trains Autoroute
3 Trains Route
2 Trains Route
0 250 500 750 1000 1250 1500 1750 2000 2250 2500
Puissance perdue [W]
100
90
Pourcentage des pertes
80
roulements [%]
70
60
50
40
30
20
10
0
3 Trains 2 Trains 3 Trains Route 2 Trains Route
Autoroute Autoroute
Figure 2-4 Part des pertes des roulements de l'arbre GV sur la somme des pertes aux roulements
Enfin pour l’hypothèse de niveau d’huile choisie la solution à trois trains est la plus
intéressante du point de vue du barbotage. Le diamètre extérieur des roues est moins grand,
même s’il y a une roue de plus le barbotage est plus faible. Les pertes par barbotage présentées
sur la Figure 2-3 ont été calculées sans prendre en compte l’aération du bain d’huile. Le chapitre 1
a montré qu’à haute température et à grande vitesse le bain d’huile pouvait s’aérer à des niveaux
de 20 à 30%. Le niveau d’aération ne peut pas être prédit précisément, pour les conditions
autoroute et route une aération de 25% semble être une hypothèse réaliste.
51
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
3 Trains Autoroute
2 Trains Autoroute
3 Trains Route
2 Trains Route
100
Pourcentage des pertes par
80
barbotage [%]
60
40
20
0
3 Trains 2 Trains 3 Trains Route 2 Trains Route
Autoroute Autoroute
Figure 2-6 Part des pertes par barbotage de l’étage GV sur les pertes par barbotage totales
Du point de vue du rendement chaque architecture a un point fort. N’utiliser que deux étages
permet de réduire les pertes aux engrènements. Cependant utiliser un troisième étage de
réduction permet de diminuer le diamètre des roues et par ce biais de limiter les pertes par
barbotage. De plus la solution à trois étages est la plus compacte et la moins lourde, elle semble
être la plus intéressante. Dans les deux cas l’étage GV provoque la majorité des pertes.
52
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
53
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Tableau 2-10 Liste des éléments utilisés dans les réseaux thermiques du réducteur HV
1 Bâti
2 Air
3 Partie inférieure du carter
4 Partie latérale du carter
5 Partie supérieure du carter
6 Bain d’huile
7,8 Roulements de l’arbre du pignon
9,10 Roulements de l’arbre de la roue
11 Arbre du pignon
12 Arbre de la roue
13 Pignon
14 Roue
15 Zone d’engrènement
Les nœuds air et bâti sont utilisés pour fixer les conditions limites. Ces éléments sont des
puits de chaleur, leur température est imposée. D’autres éléments représentent les sources de
chaleur. Ainsi il y a un nœud par roulements, un nœud pour l’engrènement et enfin un nœud
pour l’huile. Ces éléments reçoivent respectivement les pertes aux roulements, les pertes par
frottement à la denture et les pertes par barbotage. Les modèles de pertes de puissance utilisés
sont les mêmes que pour l’approche isotherme.
54
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
Les éléments du réducteur doivent être reliés pas des résistances thermiques afin de traduire
les échanges de chaleur. Cinq types de résistance sont utilisés, ils sont listés dans le Tableau 2-11.
Tableau 2-11 Liste des résistances thermiques utilisées dans le réseau thermique du réducteur HV
R1 Résistance de striction
R2 Résistance de convection avec le bain d’huile
R3 Résistance de transport avec le bain d’huile
R4 Résistance de conduction
R5 Résistance de convection et de rayonnement avec l’air
En reliant les éléments du Tableau 2-10 et les résistances thermiques du Tableau 2-11, on
obtient le réseau thermique du réducteur HV. Ce réseau thermique est présenté sur la Figure 2-8.
55
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette résistance est donnée par l’application de la théorie de la température flash de Block
[76] ramenée à une valeur moyenne (spatialement et temporellement) et s’exprime sous la forme
suivante [73] :
0,767
𝑅𝑠𝑡𝑟𝑖𝑐𝑡𝑖𝑜𝑛 = 2.1.
√2. 𝑙ℎ . 𝑙𝑐𝑚 . 𝜒. √𝑈/2
Avec
- 2. 𝑙ℎ la largeur de la surface d’écrasement dans le contact [m]
- 𝑙𝑐𝑚 longueur moyenne de contact [m]
- 𝜒 l’effusivité thermique du pignon [W⋅K-1⋅m-2⋅s1/2]
- 𝑈 somme des vitesses de roulement [m/s]
1⁄
Si 𝑅𝑒 > 5 × 105 𝑁𝑢 = 𝑃𝑟 3 × (0,037 × 𝑅𝑒 0,8 − 850) 2.3.
Pour calculer le nombre de Reynolds 𝑅𝑒 il est nécessaire de définir une vitesse caractéristique
de l’écoulement d’huile. Il faut alors distinguer deux cas : i) la surface est immergée dans le bain
d’huile ii) la surface est hors du bain mais l’huile ruissèle sur celle-ci.
Dans le premier cas la vitesse caractéristique est imposée par la rotation de la roue. Elle est
prise égale à la vitesse périphérique du mobile c’est-à-dire le produit 𝑅𝑝 . Ω.
Dans le second cas l’huile s’écoule sous l’effet de la gravité. Pour une plaque verticale de
hauteur L, sans prendre en compte les forces visqueuses et en négligeant la trainée
aérodynamique, la vitesse moyenne du fluide peut s’exprimer de la façon suivante [73]:
8. 𝑔. 𝐿 2.4.
𝑉𝑚𝑜𝑦 = √
9
Avec
- 𝑔 l’accélération de pesanteur [m/s²]
56
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
La longueur caractéristique utilisée est le diamètre primitif. La vitesse caractéristique est prise
égale à la vitesse périphérique de la roue [71]. La surface de convection est la surface immergée
du flanc de la roue dentée.
La résistance de projection traduit l’échange de chaleur dû à la projection de l’huile se
trouvant entre les dents. Les forces centrifuges projettent cette huile entre le moment où la dent
sort du bain d’huile et le moment où l’engrènement commence. C’est-à-dire pendant la rotation
sur l’angle θ de la Figure 2-9. Pendant cette période, l’huile refroidit les dents du pignon puis une
partie est centrifugée.
Pignon
Roue
n
Rp
Ω Bain d’huile
𝜃𝑝
Avec
- 𝐵 la largeur de la denture [m]
- 𝑍 le nombre de dents de la roue
- 𝑅𝑝 le rayon primitif de la roue [m]
- 𝑎𝑑 la diffusivité thermique de l’huile [m².s-1]
- 𝜈 la viscosité cinématique de l’huile [m².s-1]
- 𝐻𝑑𝑒𝑛𝑡 la hauteur de dent [m]
57
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
- 𝜃𝑝 l’angle de projection [rad]
- 𝜒 l’effusivité thermique de l’huile [J.K-1.m-2.s-1/2]
- 𝜔 la vitesse de rotation [rad/s]
- 𝜓 est un nombre sans dimension dépendant de la géométrie et de la viscosité de
l’huile. Il est défini par la formule :
1
𝑅𝑝 . 𝑎𝑑 . 𝜃 2 4 2.8.
𝜓=( )
𝜈. 𝐻𝑑𝑒𝑛𝑡
1
𝑅𝑡𝑟𝑎𝑛𝑠𝑝𝑜𝑟𝑡 = 2.9.
𝑚̇𝐶𝑝
Avec
- 𝑚̇ le débit massique d’huile entrant dans le roulement [kg/s]
- 𝐶𝑝 la capacité calorifique massique du lubrifiant [J/kg/K]
rint
rext
L
𝑅
ln ( 𝑅𝑒𝑥𝑡 ) 2.10.
𝑖𝑛𝑡
𝑅𝑡ℎ =
2. 𝜋. 𝑘. 𝐿
Avec comme représenté sur la Figure 2-10
- 𝑅𝑒𝑥𝑡 le rayon extérieur [m]
- 𝑅𝑖𝑛𝑡 le rayon intérieur [m]
58
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
Dans le cas des roues dentées cette résistance sépare la température de masse des dents et
l’élément arbre, le diamètre extérieur est pris égal au diamètre primitif, et le diamètre intérieur à
celui de l’arbre.
Dans le cas des roulements cette résistance correspond à la conduction dans les bagues. Le
diamètre moyen de la piste de roulement est utilisé tantôt pour 𝑟𝑒𝑥𝑡 et tantôt pour 𝑟𝑖𝑛𝑡 selon que
l’on calcule la résistance pour la bague intérieure ou extérieure.
b) Résistances de conduction entre les parties du carter et avec le bâti
Des résistances de conduction relient les différentes parties du carter. Elles représentent la
conduction ayant lieu dans la hauteur du carter latéral. Elles s’expriment de la façon suivante :
l 2.11.
𝑅𝑡ℎ =
𝑆. 𝑘
Avec
- 𝑆 la surface de la section horizontale du carter latéral [m²]
- 𝑙 la moitié de la hauteur du carter latéral [m]
Le carter inférieur est posé sur le bâti. La résistance de conduction prend la même forme. S
devient alors la surface de contact et l une longueur de conduction dans le bâti. En absence
d’information sur la conception du bâti l est prise égal à 100 mm.
1
𝑅𝑡ℎ = 2.12.
𝜀. 𝜎. 𝑆. (𝑇𝑝 + 𝑇𝑎𝑖𝑟 2 ). (𝑇𝑝 + 𝑇𝑎𝑖𝑟 )
2
Avec
- 𝜀 l’émissivité du matériau [-]
- 𝜎 la constante de Stefan-Boltzmann = 5,67. 10−8 [W/m².K4]
- S la surface d’échange de l’élément [m²]
- 𝑇𝑝 la température de l’élément [K]
- 𝑇𝑎𝑖𝑟 la température de l’air [K]
La température du milieu environnant est supposée être la même que celle de l’air. Aucun
nœud ne lui est attribué, la résistance de rayonnement est donc dans le réseau imposée entre le
carter, les bouts d’arbre et l’air.
59
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
b) Résistance de convection entre le carter et l’air
Sur le banc d’essais, le carter sera ventilé par un flux d’air arrivant à la verticale. Ce flux
entraine une convection forcée sur la partie supérieure du carter, ainsi que sur sa partie latérale.
Dans les deux cas l’échange est de type convection forcée sur une plaque.
La partie supérieure du carter est normale au flux d’air. Dans ce cas le coefficient d’échange
convectif est donné par l’expression suivante [72] :
Avec
- 𝑉𝑎𝑖𝑟 la vitesse du flux d’air [m/s]
- 𝐿 la dimension caractéristique égale à la surface frontale de la paroi divisée par son
périmètre [m]
La partie latérale du carter est tangentielle au flux d’air vertical. Dans ce cas le coefficient
d’échange peut être calculé avec :
0,66
ℎ𝑐𝑜𝑛𝑣 = 5,6. 𝐿−0,34 . 𝑉 2.14.
1
𝑅𝑡ℎ = 2.15.
𝜂ℎ𝑐𝑜𝑛𝑣 𝑆
Avec
- 𝜂 l’efficacité d’ailettes
Extrémité de l’arbre
Paroi
𝑟
60
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
Où
𝑉. 2. 𝑅 2.19.
𝑅𝑒 =
𝜈
Avec
- 𝑉 la vitesse périphérique du cylindre [m/s]
- 𝑅 le rayon de l’arbre [m]
- 𝜈 la viscosité cinématique de l’air [m²/s]
Et Où
Avec
- g l’accélération de pesanteur [m/s²]
- 𝛽𝑓 le coefficient de dilatation volumique de l’air [1/K]
𝑇 la température de l’arbre [K]
- 𝑇𝑎𝑖𝑟 la température de l’air [K]
- 𝑅 le rayon de l’arbre [m]
- 𝜈 la viscosité cinématique de l’air [m²/s]
61
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Route Autoroute
9% 9%
15%
25%
66% 76%
62
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Positionnement du problème
10%
9%
81%
Figure 2-13 Répartition des pertes après optimisation de l'engrènement, condition route, vitesse de l’air 20m/s
Cumuler la chute du coefficient de frottement et les corrections de profil permet de diminuer
la part des pertes à l’engrènement, voir Figure 2-13. Pour la condition route le gain représente
0,4% de rendement. Dans les deux cas, ce sont des optimisations qui ont un coût mais le gain
peut être important par rapport aux technologies utilisées dans les boites de vitesse automobiles.
L’optimisation des pertes liées à la denture n’a pas été plus poussée dans cette étude.
b) Echauffement des roulements
Les pertes aux roulements augmentent avec la vitesse, le paragraphe 2.3.2 a montré que la
majorité était imputable aux roulements de l’arbre GV. Le modèle alors utilisé était celui de SKF.
La bibliographie a relevé qu’il existe aussi un second modèle plus ancien, le modèle de Harris.
Dans ce modèle l’augmentation des pertes avec la vitesse est largement conditionnée par le
coefficient de perte hydrodynamique f0. Pour les roulements à contact oblique une fourchette est
donnée : 3,3 < 𝑓0 < 6,6. Le choix de ce coefficient dans l’intervalle est laissé à l’utilisateur, la
valeur moyenne de 𝑓0 = 5 est utilisée. La Figure 2-14 montre la différence de réponse en vitesse
du modèle SKF et Harris. La tendance n’est pas la même suivant le modèle utilisé, le modèle de
SKF prévoit des pertes beaucoup plus importantes à haute vitesse.
1500
Perte de puissance [W]
1000
500 Harris
SKF
0
0 10000 20000 30000 40000
Vitesse [tr/min]
Figure 2-14 Pertes de puissance du roulement 7204c en fonction de la vitesse pour le couple d’entrée de la condition
autoroute pour un roulement à 100°C
Les pertes de puissance et les températures sont couplées. La Figure 2-15 montre la
température du bain d’huile et les pertes de puissance du réducteur HV données par la
modélisation thermomécanique. Le coefficient de frottement denture est fixe et le barbotage est
dans ce cas sur un régime d’écoulement indépendant de la température. La diminution des pertes
de puissance avec la température n’est due qu’aux roulements. Le rendement augmente avec la
température d’huile. Cette tendance est aussi sur la Figure 2-16 qui montre les mêmes données
pour la condition route.
63
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
140
2500
120
140
2500
120
100
80 1500
60
1000
40
500
20
0 0
0 10 20 30 40 50
Vitesse de l'air [m/s]
est plus faible. Cependant la température des roulements GV reste très supérieure à la
température d’huile.
Dans la modélisation thermomécanique du réducteur GV les roulements n’ont qu’une
température représentant la température moyenne des éléments constituant le roulement. Le
carter est plus froid que les arbres et il faut s’attendre à ce que la température du roulement ne
soit pas homogène. La modélisation thermomécanique des roulements doit être améliorée.
Tableau 2-13 Perte de puissance et température des roulements de l'arbre GV, vitesse de l’air 20m/s
Autoroute
7204C 7004C
Perte roulement SKF [W] 934 630
Différence de température avec l'huile [°C] 35,4 26,4
Route
7204C 7004C
Perte roulement SKF [W] 442 281
Différence de température avec l'huile [°C] 17,5 13,6
L’aération augmente largement les pertes par barbotage du système. Plus il y a de puissance
dissipée dans le bain d’huile plus celui-ci voit sa température augmenter. L’aération fait passer de
9 à 12% la participation du barbotage aux pertes du réducteur HV, voir Figure 2-17.
12%
14%
74%
Figure 2-17 Répartition des pertes dans le cas d’une aération de 25%, condition autoroute, vitesse de l’air 20m/s
Dans le cas d’un réducteur à plusieurs étages l’aération provoquée par la roue GV impacte
aussi les autres mobiles augmentant les pertes par barbotage provoquées par ceux-ci. Il serait
65
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
intéressant de limiter l’aération provoquée par la roue de l’étage GV et ainsi diminuer le
barbotage provoqué par l’ensemble des roues.
2.6 Conclusion
La conception d’un réducteur mécanique répond à un cahier des charges. Pour répondre à la
demande deux conceptions reposant sur les mêmes règles de dimensionnement sont proposées.
Pour réaliser le rapport de réduction de 32 en utilisant les technologies classiques deux solutions
s’opposent :
- Utiliser deux trains d’engrenages et avoir des grands rapports unitaires et donc des
grandes roues
- Utiliser trois trains et avoir des roues plus petites au prix d’un engrènement et de deux
roulements supplémentaires
Une première comparaison du rendement des architectures a été réalisée en supposant
isothermes les transmissions mécaniques. L’étude révèle que l’utilisation de deux étages permet de
réduire les pertes aux engrènements et de limiter les pertes aux roulements. A haute vitesse le
bain d’huile peut s’aérer, se charger d’air. Dans ce cas utiliser un troisième étage de réduction
permet de diminuer le diamètre des roues et par ce biais de limiter les pertes par barbotage. Au
final dans le cas le plus réaliste où il y a de l’aération, les deux architectures ont un rendement
similaire. De plus la solution à trois étages est la plus compacte et la moins lourde, elle semble
être la plus intéressante.
La comparaison isotherme montre aussi que l’étage GV pose problème du point de vue du
rendement. Il est pour les deux architectures responsable d’une grande partie des pertes dans les
roulements et par barbotage. Cela a motivé une étude plus poussée de l’étage GV du réducteur
trois trains. Un projet de banc d’essai, limité à l’étage GV, a été modélisé thermo-mécaniquement
à l’aide de la méthode des réseaux thermiques. Le modèle permet de coupler pertes de puissance
et température.
Les pertes à l’engrènement sont très dépendantes du coefficient de frottement denture. Ce
coefficient peut être optimisé en utilisant des technologies qui ne sont pas actuellement utilisées
en automobile comme par exemple la superfinition, ou encore des traitements physico-chimiques
de surface des dentures. De plus les corrections de profil de denture permettent aussi d’optimiser
encore le rendement. Ce sont des solutions coûteuses mais un gain non négligeable est accessible.
L’étude révèle deux problématiques. La première est liée aux roulements de l’arbre GV. A
grande vitesse sur autoroute ils provoquent beaucoup de pertes de puissance quel que soit le
modèle de perte de puissance utilisé. Leur échauffement devient critique. Ils sont dans le modèle
représentés par une seule température. Ce choix de modélisation gomme le gradient de
température dans le roulement. Le chapitre 3 s’attachera à améliorer la prédiction des pertes aux
roulements GV et à investiguer les gradients de température internes à ces paliers très sollicités.
La seconde problématique concerne les pertes par barbotage provoquées par la roue GV.
Elles sont importantes à cause de la grande vitesse de rotation. Cette roue tourne à grande vitesse
dans un bain d’huile à haute température. Elle provoque l’aération de l’huile. Cette aération
dégrade le rendement de la transmission au complet augmentant les pertes par barbotage des
roues moyenne et basse vitesse. Le chapitre 4 s’attachera à proposer et étudier une méthode pour
réduire à la fois les pertes par barbotage mais aussi l’aération créée par la roue GV.
66
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Chapitre 3
Etude des roulements Grande Vitesse
68
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
Figure 3-1 Schéma de l'ellipse de Hertz discrétisée en tranches les vecteurs représentent les forces de frottement
tangentielles en bleu et les vitesses de glissement locales en vert
Enfin grâce à la connaissance des coefficients et des vitesses de frottement ainsi que des
efforts, la puissance dissipée par chaque tranche est calculée. La puissance totale dissipée dans un
contact étant la somme des puissances locales.
69
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
hydrodynamique. Cette source de perte vient impacter l’équilibre mécanique du roulement,
cependant le modèle de NTN-SNR ne prend pas en compte ce phénomène.
La résistance au roulement est caractérisée par une distance ici appelée δ (Figure 3-3). Grâce
à une intégration numérique Tevaarwerk [79] propose pour un contact EHL la formulation
suivante :
Ainsi la perte de puissance par roulement hydrodynamique dans un contact s’exprime sous la
forme :
𝑄𝑟𝑜𝑢𝑙𝑒𝑚𝑒𝑛𝑡ℎ𝑦𝑑𝑟𝑜 = W. (ω1 + ω2 ). δ 3.2.
Avec :
- 𝑄𝑟𝑜𝑢𝑙𝑒𝑚𝑒𝑛𝑡ℎ𝑦𝑑𝑟𝑜 perte de puissance par roulement hydrodynamique [W]
- δ décalage du centre d’application de charge [m]
- 𝑊 charge normale au contact [N]
- ω1 vitesse de rotation du solide 1 par rapport au contact [rad.s-1]
- ω2 vitesse de rotation du solide 2 par rapport au contact [rad.s-1]
- 𝐺 ∗ paramètre adimensionné de matériau [-]
- 𝑈 ∗ paramètre adimensionné de vitesse [-]
- 𝑊 ∗ paramètre adimensionné de charge [-]
- a demi-axe de l’ellipse de contact dans le sens du roulement [m]
- b demi-axe de l’ellipse de contact dans le sens perpendiculaire du roulement [m]
70
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
2𝜇𝐿𝜋𝑅 3 𝛥𝜔 3.3.
𝐶𝑓 =
𝐽𝑟
Où :
- 𝐶𝑓 couple de frottement d’un palier fluide [N.m]
- μ viscosité dynamique du lubrifiant [Pa.s]
- L longueur du palier [m]
- R rayon moyen du palier [m]
- Δω différence de vitesse angulaire au sein du palier [rad.s-1]
- 𝐽𝑟 jeu radial [m]
En prenant l’hypothèse que le jeu entre la cage et les billes est totalement gavé d’huile, cette
approche peut être utilisée pour exprimer une perte de puissance. Ainsi dans le cas du contact
lubrifié cage/bille la perte de puissance associée à la rotation orthogonale de la bille par rapport à
la poche s’exprime sous la forme :
Cette formulation a été appliquée aux roulements de l’arbre GV du réducteur étudié. D’après
la formule la perte de puissance augmente avec la vitesse et la viscosité. Pour les roulements
7204C étudiés à 42000 tr/min, lubrifiés par une huile à 40°C (pour maximiser l’effet de la
viscosité) la perte calculée est de moins de 2W. Il a été décidé de négliger cette source de
dissipation thermique.
𝜌𝑒𝑓𝑓 𝐶𝐷 𝐴 𝜔𝑚 𝑑𝑚 3 3.5.
𝑄𝑡𝑟𝑎î𝑛é𝑒 =𝑍 ( )
2 2
𝜋𝐷 3.6.
𝐴 = 𝐷( − ec )
4
Où :
- 𝑄traînée perte de puissance par traînée aérodynamique [W]
- A maître couple d’une bille auquel est retranché l’épaisseur de la cage [-]
- Z nombre de billes [-]
- D diamètre de bille [m]
- ec épaisseur de la cage [m]
- 𝐶𝐷 coefficient de traînée d’une bille [-]
71
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
- 𝑑𝑚 diamètre primitif du roulement [m]
- 𝜔𝑚 vitesse orbitale de la bille par rapport au carter [rad.s-1]
- 𝜌𝑒𝑓𝑓 densité équivalente du brouillard d’huile [kg.m-3]
Dans [29] Harris utilise le coefficient de traînée d’une sphère isolée dans un milieu infini. Il
ne prend pas en compte la cage, les bagues et l’influence des autres billes. La Figure 3-4 montre
l’évolution du coefficient de traînée en fonction du nombre de Reynolds suivant :
𝑉𝐷 3.7.
𝑅𝑒 =
𝜈
Où :
- 𝑉 vitesse d’écoulement du fluide par rapport au solide [m.s-1]
- 𝜈 viscosité cinématique du fluide [m²/s]
Figure 3-4 Evolution du coefficient de traînée d’une sphère et d’un cylindre isolés en fonction du nombre de
Reynolds
F.Pouly [83] et Y.Marchesse [84] ont étudié l’influence des autres billes sur le coefficient de
trainée. Ils montrent que le rapport de l’espace entre les billes à leur diamètre a un fort impact sur
le 𝐶𝐷 . Pour un rapport proche de un, ils concluent que prendre 20% du coefficient d’une bille
dans un milieu infini est une bonne approximation.
Pour calculer la perte par trainée aérodynamique il faut aussi connaitre les caractéristiques du
fluide se trouvant dans le roulement. Ce fluide est un brouillard d’huile dont les caractéristiques
sont très dépendantes de la fraction d’huile dans l’air. Pouly [83] montre que la relation suivante
donne de bons résultats pour déterminer la fraction d’huile dans le roulement :
Où :
- X fraction d’huile dans le roulement [-]
- 𝐶𝐷 coefficient de traînée d’une bille [-]
- 𝑚̇𝑜𝑖𝑙 débit massique de lubrifiant [kg/s]
72
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
1200
1000
800
Puissance [W]
0
Route Autoroute
Condition opératoire
Figure 3-5 Perte de puissance du roulement 7204C calculées avec la modélisation local, le modèle de SKF et de
Harris pour une huile à 100°C
73
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
catalogue, dans le modèle local il est calculé avec une loi de frottement semi empirique qui prend
en compte les paramètres des contacts.
a)
600
Puissance dissipée [W]
550
b)
250
Puissance dissipée [W]
230
210
Pertes par glissement SKF
190
Pertes par glissement du
170
modèle local
150
20 40 60 80 100 120
Température du lubrifiant [°C]
Figure 3-6 Comparaison des pertes par glissement, roulement 7204 cas a) autoroute b) route
La Figure 3-7 montre la comparaison entre les pertes par roulement calculées avec la formule
de Tevaarwerk et par la méthode SKF. Les courbes montrent la même dépendance à la
température. Cette tendance vient du facteur 𝜈 0.6 présent dans les deux formulations. Les pertes
par roulement hydrodynamique par les deux méthodes sont comparables pour les plages de
température de fonctionnement visées.
1000
Puissance dissipée [W]
800
600
Pertes par roulement SKF
400
Pertes par roulement
200 hydrodynamique
0
60 70 80 90 100 110 120
Température du lubrifiant [°C]
Figure 3-7 Comparaison des pertes par roulement, roulement 7204 cas autoroute
74
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
Enfin, la traînée aérodynamique des billes est calculée par le modèle local. Ces pertes sont
représentées sur la Figure 3-8. Cette source de perte de puissance est relativement faible par
rapport aux pertes aux contacts. Bien que le produit (N*dm) du 7204 de l’étage GV soit de 1,4
million, les éléments roulants ont un faible diamètre (7mm environ). Avec la température, la
densité de l’huile présente dans le brouillard diminue et modifie légèrement les pertes par traînée.
10
Puissance dissipée [W]
8
6
Pertes par trainée
4
2
0
20 70 120
Température du lubrifiant °C
Figure 3-8 Pertes par traînée aérodynamique, roulement 7204 cas autoroute
Les roulements GV sont lubrifiés par un filet d’huile amené par une goulotte, ce mode de
lubrification se rapproche de l’injection. Le modèle de perte de puissance de SKF propose de
calculer la perte engendrée par la lubrification par injection en doublant les pertes par barbotage
d’un roulement faiblement immergé. De plus le débit d’injection n’a pas d’impact dans le calcul.
Les pertes prédites par cette méthode sont très supérieures à celles obtenues avec le modèle local
(cf. Figure 3-9). Aucune borne d’application n’est donnée mais cette méthode ne semble pas
cohérente pour les grandes vitesses.
600
Puissance dissipée [W]
500
400
300 Pertes par injection
200
100
0
20 70 120
Température du lubrifiant °C
Figure 3-9 Pertes par injection SKF, roulement 7204 cas autoroute
Les pertes dans les contacts billes/bagues sont majoritaires. Pour cette source de chaleur le
model local et le modèle SKF ont les mêmes tendances. Chacun calcule des pertes par frottement
et des pertes pas roulement hydrodynamique. Le modèle SKF utilise un coefficient de frottement
qui conditionne les pertes par glissement, celui-ci est donné par le catalogue de manière indicative
selon le type de roulement et le type d’huile (synthétique ou minérale). Ce même coefficient est,
dans le modèle local, calculé avec des essais. En recalant le coefficient de frottement, le modèle
SKF peut donner des totaux de pertes de puissance proches de ceux prédits avec l’approche
locale. Cependant à la différence du modèle SKF cette approche permet de calculer la localisation
des pertes de puissance dans le roulement.
75
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
3.2.4.2 Localisation des sources de chaleur à l’intérieur du roulement
Figure 3-10 Localisation des pertes de puissance dans un roulement à contact oblique
La puissance est dissipée au sein d’un roulement en plusieurs endroits (cf. Figure 3-10). La
perte par trainée aérodynamique est dissipée dans le brouillard interne au roulement de manière
axisymétrique.
Les pertes aux contacts sont par définition localisées à chacun des contacts entre billes et
bagues. Les géométries des pistes de roulement intérieure et extérieure sont différentes. Les
pertes aux contacts ne sont pas équi-réparties entre la bague intérieure et extérieure. Un exemple
est donné Figure 3-11.
1400
1200
Puissance dissipée [W]
1000
800
Bague intérieure
600
Bague extérieure
400
200
0
20 40 60 80 100 120
Température du lubrifiant [°C]
Figure 3-11 Pertes aux contacts billes/bague intérieure et billes/bague extérieure, roulement 7204 condition
autoroute
De plus les roulements GV voient deux types d’effort, l’un radial dû au couple transmis et un
autre axial dépendant de la précharge, de l’angle de contact du roulement, et du couple transmis
par l’intermédiaire de l’angle d’hélice du pignon. Cela est surtout sensible pour les 7004 de l’arbre
GV. Le chargement des contacts dépend alors de la position des éléments roulants, voir Figure
3-12.
76
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
150
50
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Numéro de la bille considérée
Figure 3-12 Chargement des contacts billes/bagues, roulement 7004 condition autoroute
Les chargements étant différents les conditions de contact le sont aussi. Les pertes aux
contacts ne sont donc pas axisymétriquement réparties dans le roulement. Cela a un impact
différent selon la bague considérée.
Dans le cas d’un réducteur, c’est la bague intérieure qui voit tourner l’effort. Cela veut dire
que la piste de roulement intérieure voit passer la zone de charge à chaque révolution. Moyennée
dans le temps, il n’y a donc pas de zone de dissipation préférentielle sur la bague intérieure. A
l’inverse la bague extérieure est fixe par rapport à l’effort. Il y a donc une partie de la piste
extérieure qui est plus chargée. Il y a donc une zone de la bague extérieure qui reçoit plus de
chaleur, les billes de 1 à 4 et 11 à 14 sur la Figure 3-13. Cette zone aura tendance à s’échauffer
plus que la partie moins chargée de la bague extérieure.
35
Puissance dissipée [W]
30
25
20
15
10
5
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Numéro de la bille considérée
Figure 3-13 Pertes aux contacts billes/bague extérieure, roulement 7004 condition autoroute à 90°C
Il y a une forte discontinuité entre les billes 5/6 et 10/11. Pour le cas de charge autoroute, le
7004 voit 4,5 fois plus de charge radiale que de charge axiale. Les billes 6 à 10 décollent de la
bague intérieure et ne transmettent plus de charge. Les contacts entre les billes et la bague
extérieure sont alors légèrement chargés par les efforts centrifuges et non plus par la précharge.
77
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
3.3 Modélisation thermique
0,7
0,6
0,5
0,4
[mm]
0,3 b
a
0,2
0,1
0
1 2 3 4 5 6 7 8 9 10 11
Numéro de la bille considérée
Figure 3-14 Taille des contacts billes/bague extérieure, roulement 7204 condition autoroute
78
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
La liste des éléments considérés dans le réseau thermique d’approche est présentée dans le
Tableau 3-1, les résistances thermiques utilisées dans le Tableau 3-2 et le réseau thermique sur la
Figure 3-15. Les nœuds du roulement sont désignés avec des chiffres. Les nœuds appartenant au
réseau thermique du réducteur le sont avec des lettres, ces nœuds sont communs aux deux
réseaux.
79
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Tableau 3-2 Liste des résistances thermiques utilisées dans le réseau thermique d'approche
R1 Résistance de striction
R2 Résistance de convection avec le brouillard d’huile
R3 Résistance de transport avec le bain d’huile
R4 Résistance de conduction
Les pertes de puissance par traînée aérodynamique sont injectées dans le nœud
correspondant au brouillard d’huile. Les pertes aux contacts sont calculées et imposées aux
nœuds 3, 4 et 6.
𝑉 ∗𝐷
𝑅𝑒 = 3.9.
𝜈
Nombre de Prandtl :
𝜈
𝑃𝑟 = 3.10.
𝑎𝑑
Avec
- 𝑉 Vitesse du fluide [m/s]
- 𝐷 Diamètre de la bille [m]
- ν viscosité cinématique du fluide [m²/s]
- 𝑎𝑑 diffusivité thermique du fluide [m²/s]
0.25
0.54 ) 0.3
𝜇∞
𝑁𝑢 = (1.2 + 0.53 ∗ 𝑅𝑒 ∗ 𝑃𝑟 ∗( ) 3.11.
𝜇𝑝
Avec
- 𝑁𝑢 nombre de Nusselt
- μ∞ viscosité dynamique du fluide hors de la couche limite [Pa.s]
- μp viscosité dynamique du fluide à la paroi de la sphère [Pa.s]
80
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
Nombre de Reynolds :
𝑉 ∗𝐿
𝑅𝑒 = 3.12.
𝜈
Nombre de Prandtl :
𝜈
𝑃𝑟 = 3.13.
𝑎𝑑
Avec
- V∞ vitesse du fluide hors de la couche limite [m/s]
- L longueur de la plaque, ici π fois le diamètre de la surface considérée [m]
- ν viscosité cinématique du fluide [m²/s]
- 𝑎𝑑 diffusivité thermique du fluide [m²/s]
5
Pour
{𝑅𝑒 < 5 ∗ 10 𝑁𝑢 = 0,664 ∗ 𝑅𝑒 1⁄2 ∗ 𝑃𝑟 1⁄3 3.14.
𝑃𝑟 > 0.6
5
Pour
{5 ∗ 10 < 𝑅𝑒 𝑁𝑢 = (0,037 ∗ 𝑅𝑒 0,8 − 871) ∗ 𝑃𝑟 1⁄3 3.15.
𝑃𝑟 > 0.6
Avec
- 𝑁𝑢 nombre de Nusselt
Il est nécessaire de calculer la vitesse hors de la couche limite. Pour ce faire, la vitesse est définie
ainsi :
- Pour la bague extérieure, la vitesse tangentielle des billes au diamètre moyen du roulement par
rapport au carter
-Pour la bague intérieure, la différence de vitesse tangentielle entre les billes et le diamètre moyen
de la piste de roulement intérieure.
𝜋𝑑𝑚 1
R th = ∗ 3.16.
𝑘𝑆 2
Avec
- 𝑑𝑚 Diamètre moyen de la bague [m]
- 𝑘 Conductivité thermique de la bague [m]
- 𝑆 Section de la coupe axiale de la bague [m²]
81
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Pour Pe =
a∗𝑈1
> 10 0.918
𝑎𝑑 R striction = 3.17.
2 ∗ b ∗ χ ∗ √2 ∗ a ∗ √𝑈1
Pour Pe < 1 0,27
R statique = 3.18.
k∗a
Pour 1 < Pe < 10 1 1 1
2 = 2 + 2 3.19.
R th R striction R statique
82
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
Les billes et les contacts billes/bagues sont isothermes : il n’y a plus besoin de
connaitre la répartition des pertes entre les bagues extérieure et intérieure
Le réseau thermique correspondant est représenté sur la Figure 3-17. Il utilise 4 nœuds pour
décrire le roulement (Tableau 3-3). Il n’y a plus qu’un nœud pour la bague extérieure. Les billes et
les contacts étant isothermes un seul nœud les représente. Ce nœud peut recevoir la totalité des
pertes aux contacts.
Tableau 3-3 Liste des noeuds du réseau thermique avec un nœud pour les billes et les contacts
1 Bague Extérieure
2 Billes et contacts avec les bagues
3 Bague Intérieure
4 Brouillard d’huile
A Bain d’huile
B Arbre portant le roulement
C Carter portant le roulement
Figure 3-17 Réseau thermique avec un nœud pour les billes et les contacts
83
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Figure 3-18 Schéma d'un banc FZG
Une autre utilité de la machine FZG est l’étude du rendement des transmissions mécaniques
[67], [71], [91]–[93]. A cette fin le LaMCoS dispose d’un exemplaire qui a été instrumenté. Des
capteurs de température sont installés sur les parties fixes et un couple-mètre permet la mesure
du couple nécessaire à la mise en rotation des deux boites à engrenages. De plus des capteurs de
température ont été installés sur les bagues intérieures des roulements du boitier de test. Ils
permettent de connaitre la différence de température entre les bagues des roulements.
84
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
J.Durand et al. [71] ont modélisé un réducteur de la machine FZG en utilisant la méthode
des réseaux thermiques. Le réseau thermique résultant est présenté sur la Figure 3-19, le label des
éléments se trouve sur le Tableau 3-4. Le réseau est très proche du modèle présenté dans le
chapitre 2. Il y a deux différences : i) il n’y a qu’un nœud limite, l’air ii) il n’y a pas d’échange
direct entre les roulements et le bain d’huile (sur la base des hypothèses retenues par J.Durand et
al. [71]).
Le modèle prend en compte les pertes de puissance présentées dans le chapitre 1 avec les
modèles suivants :
A la denture avec la relation de Velex et Ville [16] et un coefficient de frottement
recalé grâce aux essais.
Dans les roulements avec le modèle SKF avec un coefficient 𝜇𝑠𝑙 recalé.
Le barbotage avec le modèle de C.Changenet [54].
Et enfin des pertes dans les joints avec la formule du catalogue Simrit [65]. Ces pertes
sont imposées directement aux nœuds représentant les arbres.
Dans le cas étudié les deux boites ont des configurations identiques : même couple de
dentures et mêmes roulements. Les roulements utilisés sont des 6406, les dentures sont de type
C, et l’huile est de la Total EP220. Les caractéristiques et l’instrumentation détaillées de la
machine FZG sont présentées à l’annexe 2. Le modèle développé a été validé en utilisant les
températures des parties fixes et les couples mesurés sur le banc FZG du LaMCoS. Les Figure
3-20 et Figure 3-21 montrent que le modèle reproduit bien le comportement thermique global du
banc FZG.
630
620
Perte de puisancce totale [W]
610
Essai
600 Simulation
590
580
570
560
550
540
20 25 30 35 40 45 50
Température d'huile [°C]
Figure 3-20 Comparaison entre la perte de puissance total mesurée et calculée (135 N.m, 2000 tr/min, niveau d’huile
Figure 3-22)
85
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
50
30
Essai
20 Simulation
10
0
0 500 1000 1500 2000
Temps [s]
Figure 3-21 Comparaison entre la température du bain d’huile mesurée et calculée (135 N.m 2000 tr/min, niveau
d’huile Figure 3-22)
Axes
-15mm
Niveau d’huile
86
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
35
Différence de température
30
25
20
[°C]
15
10
5
0
0 10 20 30
Temps [min]
Figure 3-23 Comportement typique de la différence de température entre bague intérieure et extérieure, niveaux
d’huile 15mm sous axe, cas 3000rpm, 135N.m
La Figure 3-23 montre un exemple de différence de température entre la bague intérieure et
la bague extérieure d’un roulement de la FZG. Cette courbe souligne que pendant le démarrage la
bague intérieure est significativement plus chaude. Cela peut être attribué au couple de frottement
des roulements qui est important à froid. Après le pic la différence se stabilise vers les 14°C lors
de cet essai. Cette différence évolue peu avec le temps bien que l’ensemble de la transmission
continue à s’échauffer.
L’investigation ne porte pas sur la caractérisation du démarrage à froid. L’étude de la
modélisation des roulements de la FZG ne se concentre donc pas sur ce point et se concentre sur
des durées où la différence de température entre les bagues est stabilisée.
3.4.3 Résultats
Le modèle de roulement fin à 8 nœuds utilise les pertes de puissance calculées par un modèle
local alors que le modèle original à un nœud et la modélisation simplifiée à 4 nœuds utilisent les
formules de SKF. La Figure 3-24 montre les pertes de puissance des roulements de l’arbre GV,
les puissances calculées sont comparables.
250
Perte de puissance [W]
200
150
100 Modèle SKF
50 Modèle local
0
20 70 120
Température d'huile [°C]
Figure 3-24 Pertes de puissance d’un roulement de l’arbre primaire GV calculées avec différents modèle de pertes
La Figure 3-25 montre la comparaison entre les températures mesurées et calculées de la
FZG après 30 minutes d’essais. En utilisant une approche locale des roulements la température
d’huile est correctement prédite, les changements apportés à la modélisation thermique du banc
FZG impactent peu le comportement thermique prédit. La modélisation fine des roulements est
87
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
capable de reproduire la différence de température de 14°C entre les bagues. Enfin elle reproduit
aussi un gradient dans la bague extérieure, celui-ci est faible.
100
95
Temperature [°C]
90
85
80
75
70
Bain d'huile Bague extérieure, Bague extérieure, Bague intérieure
partie chargée partie non-
chargée
Essai Modélisation fine
100
95
Temperature [°C]
90
85
80
75
70
Bain d'huile Bague Bague Bague Bague
extérieure, extérieure, extérieure intérieure
partie chargée partie non-
chargée
Modélisation fine Modélisation simplifiée
88
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
89
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
90
Température [°C ]
85
80
75
Roulement Bague intérieure Bague extérieure, Bague extérieure,
partie chargée partie moins
chargée
Isotherme 8 nœuds
a)
135
Température [°C ]
130
125
120
Roulement Bague intérieure Bague extérieure, Bague extérieure,
partie chargée partie moins
chargée
Isotherme 8 nœuds
b)
Figure 3-27 Températures du roulement 7204C de l’arbre GV en fonction de la modélisation thermique des
roulements a) route b) autoroute
La température prédite par le modèle isotherme n’est pas la température moyenne des bagues
prédite par l’approche fine. L’approche isotherme sous-estime la température moyenne du
roulement. Cela a un impact sur la prédiction de la température du pignon (voir Figure 3-28).
Discrétiser les roulements GV augmente la température du pignon. Ceci peut s’expliquer ainsi :
A l’engrènement il y a dissipation d’énergie par frottement. Le flux de chaleur est
dissipé en grande partie par la convection avec le bain d’huile.
Cependant une partie de celui-ci s’écoule à travers le pignon, l’arbre et enfin par les
roulements.
Ne plus négliger les gradients de température dans les roulements revient à isoler
thermiquement l’arbre du carter. Les pertes à l’engrènement ayant plus de mal à être
dissipées, la température du pignon augmente.
90 115
Température [°C ]
Température [°C ]
89 114
88 113
87 112
86 111
85 110
Pignon Pignon
90
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
110
Température [°C]
100
90
80
70
60
Bague Bague Bague Bague Huile Pignon
intérieure extérieure extérieure, extérieure,
partie partie moins
chargée chargée
4 nœuds 8 nœuds
a)
130
Température [°C]
120
110
100
90
80
Bague Bague Bague Bague Huile Pignon
intérieure extérieure extérieure, extérieure,
partie partie moins
chargée chargée
4 nœuds 8 nœuds
b)
Figure 3-29 Températures du roulement 7204C de l’arbre GV, de l’huile et du pignon en fonction de la modélisation
thermique des roulements a) route b) autoroute
La température critique du roulement est celle de la bague intérieure. Sur route elle est 22°C
plus élevé que la température d’huile. Sur autoroute par contre il y a potentiellement un problème
avec un deltaT de 45°C. Il y a deux façons de mieux refroidir le roulement, augmenter
l’immersion relative de la roue et de mieux refroidir la ligne d’arbre ou augmenter le débit de
lubrification.
91
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
3.5.3 Influence du niveau d’huile
L’utilisation du modèle de roulement simplifié est pertinente, les études de sensibilité sont
menées avec ce modèle. La Figure 3-30 montre les températures prédites pour plusieurs
immersions relatives de la roue (h/Rp voir Figure 1-6). Pour rappel 25% renvoie au niveau utilisé
jusqu’ici. Le débit arrivant dans les roulements est considéré comme constant, égal à 5ml/s.
L’augmentation de niveau d’huile permet d’augmenter les échanges thermiques entre la roue
et l’huile, elle a donc un impact positif sur la différence de température entre l’huile et le pignon.
Cela dit le gain est faible, en effet une grande partie de l’échange est réalisé par le phénomène de
projection qui est peu sensible à la hauteur d’huile (du moment où celui-ci est plus grand que la
hauteur de dent).
Par contre l’augmentation du niveau d’huile implique plus de perte par brassage. Il y a donc
d’un côté une faible réduction de la différence de température entre l’huile et le pignon, et d’un
autre coté une augmentation de la température d’huile. Cela est d’autant plus sensible à grande
vitesse, donc sur autoroute où les pertes par barbotage passent de 172W à 353W entre les
immersions de 15% et 45%. Il est peu intéressant d’augmenter le niveau d’huile au-delà de 25%.
100
Température [°C]
90
80
70
60
50
Bague extérieure Bague intérieure Pignon Huile
140
Température [°C]
130
120
110
100
90
80
Bague extérieure Bague intérieure Pignon Huile
92
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude des roulements Grande Vitesse
100
Température [°C]
90
80
70
60
50
Bague extérieure Bague intérieure Pignon Huile
140
Température [°C]
120
100
80
60
Bague extérieure Bague intérieure Pignon Huile
93
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
3.6 Conclusion
Les roulements de l’arbre GV produisent beaucoup de pertes de puissance. Le chapitre 2 a
montré qu’ils pouvaient subir un échauffement problématique. Afin de mieux comprendre cet
échauffement ce chapitre expose une modélisation locale des roulements. Pour ce faire un
modèle de perte de puissance local a été développé, il prend en compte :
les pertes par frottement avec des résultats fournis par NTN-SNR
les pertes par roulement hydrodynamique
les pertes par trainée aérodynamique
Cette modélisation des pertes liées au roulement permet de connaitre leur localisation. Elle
vient alimenter un réseau thermique d’approche qui modélise les échanges thermiques dans les
roulements. Ce modèle utilise 8 températures pour décrire le comportement thermique d’un
palier à billes. Appliqué à un réducteur tournant à 3000 tr/min il permet de reproduire les
différences de températures mesurées entre les bagues des roulements.
Le modèle de pertes de puissance local peut s’avérer difficilement ajustable avec les résultats
mesurés sur un banc d’essai roulement. En effet ce type d’équipement ne donne que la perte
globale dans le roulement et non sa répartition. Il est cependant possible de recaler un modèle de
perte global. Le modèle global de SKF a les mêmes tendances que le modèle local. Il peut être
recalé en intervenant sur le coefficient de frottement.
Pour pouvoir utiliser ce type de modèle facilement recalable, le réseau thermique fin est
simplifié. En utilisant quatre nœuds il n’est plus nécessaire de connaître la localisation des pertes
aux contacts. Cette modélisation simplifiée donne des résultats équivalents à l’approche fine. Une
modélisation plus fine des roulements conduit à des températures de fonctionnement des pièces
tournantes plus élevées.
En utilisant cette modélisation simplifiée des roulements une étude de sensibilité sur le
niveau d’huile a été menée. L’augmentation de niveau d’huile augmente les pertes de puissance.
Cependant le pignon est légèrement mieux refroidi. Le gain sur la différence de température entre
l’huile et le pignon est gommé par l’augmentation générale des températures de la transmission. A
grande vitesse il n’y a pas d’intérêt thermique à avoir un niveau d’huile élevé.
Un niveau d’huile haut permet cependant de récupérer plus de projection d’huile pour
lubrifier les roulements. L’étude de sensibilité s’est donc portée le débit de lubrification arrivant
au roulement. Amener plus d’huile au roulement signifie augmenter la densité du brouillard
d’huile et augmenter les pertes par trainée aérodynamique. Cependant ces pertes sont faibles car
les éléments roulants ont un faible diamètre dans l’application visée. Amener plus d’huile au
roulement permet de mieux le refroidir. Les paliers étant plus froids le pignon bénéficie aussi
d’une baisse de température. Pour limiter les échauffements locaux il est intéressant d’amener
beaucoup d’huile au roulement. Ceci ne peut se faire qu’en gardant une immersion significative
de la roue. Cette immersion crée des pertes par barbotage ainsi qu’une aération du bain d’huile
néfaste au rendement de la transmission, elles seront étudiées dans le chapitre suivant.
94
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Chapitre 4
Etude du barbotage
96
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
permettent d’atteindre une vitesse périphérique de 60m/s soit une valeur comparable aux 70m/s
atteinte sur autoroute par le futur réducteur HV. De plus trois huiles de boite de vitesses
automobile ont été utilisées (listées sur le Tableau 4-2).
Le niveau d’huile dans le banc est fixé à l’aide d’une pompe de remplissage et d’une vanne de
vidange. Lors de la mise en rotation de la roue une partie de l’huile est projetée sur les parois du
carter. Il s’établit alors un niveau d’huile dynamique qui est plus faible que le niveau statique dans
le banc. Le niveau d’huile dynamique est quasiment indépendant de la vitesse de rotation à partir
de 1500 tr/min. Dans les essais il est réglé à la vitesse de 2000 tr/min, à cette vitesse le bain
d’huile est calme et la lecture du niveau d’huile aisée. Dans la suite de l’étude, le niveau d’huile est
traduit par l’immersion relative de la roue dentée. De la même manière que dans les autres
chapitres, ce facteur est caractérisé par le rapport entre la hauteur immergée et le rayon primitif
de la roue.
Courroie
Boite à roulements 1
Moteur
Couple-mètre
Boite à roulements 2
97
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Tableau 4-1 Données géométriques des roues dentées
Nom G1,5 G3 G3B G5 P1,5 P3 P5
Module normal [mm] 1,5 3 3 5 1,5 3 5
Nombre de dents 102 53 51 30 64 30 20
Largeur [mm] 14 24 25 25 14 24 25
Rayon primitif [mm] 76,5 79,5 79,2 75 48 45 50
Angle de pression 20 20 20 20 20 20 20
normal [°]
Angle d’hélice [°] 0 0 15 0 0 0 0
La mesure des pertes par barbotage se fait par un couple-mètre FGP CD1140-1a ayant une
précision de ±0,002 N.m sur la pleine échelle de 2 N.m. Le couple généré par la boite à
roulement 2 est dans un premier temps caractérisé sans mettre de roue dans le carter. Cette valeur
est soustraite à la mesure globale du couple pendant les essais pour avoir accès aux seules pertes
par barbotage.
La surface libre est perturbée par la rotation de la roue dentée. Le couple mesuré pendant les
essais n’est pas constant et oscille autour d’une valeur moyenne. Les écarts à la valeur moyenne
sont supérieurs à la précision du capteur, l’incertitude sur la mesure de couple ne sera donc pas
considérée dans la suite du chapitre. Seul le couple de barbotage moyen est considéré, il est
obtenu en moyennant deux secondes d’acquisition.
La Figure 4-3 expose trois essais identiques faits à des dates différentes. Les mesures
montrent les mêmes tendances. La faible différence de valeur moyenne entre les groupes de
points s’explique par une variation de l’immersion. Celle-ci est réglée au mieux grâce à des
graduations, une erreur de plus ou moins 1mm est admise. Les résultats exposés dans la suite du
chapitre sont alors obtenus en moyennant trois essais.
600
500
Perte de puissance [W]
400
300
200
100
0
30 40 50 60 70 80 90 100
Température [C°]
Figure 4-3 Pertes mesurées pour trois essais identiques : Roue G3, 7000 tr/min, Huile 2, h/Rp=0,25
98
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
99
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Fenêtre de Pompe
réception
Pignon
H
Niveau haut
L Niveau bas
100
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
La Figure 4-7 expose le principe de la mesure. Une source radioactive au Cadmium de faible
puissance est placée dans le bain d’huile. Les rayons X vont être alors en partie absorbés par le
bain d’huile puis captés par un détecteur. Plus le bain d’huile est aéré moins il absorbe le
rayonnement.
Huile
Ià vide I
Compteur de rayon X
𝐴% 𝐴%
𝜌𝑒𝑞 = 𝜌ℎ𝑢𝑖𝑙𝑒 (1 − ) + 𝜌𝑎𝑖𝑟 ( ) 4.2.
100 100
L’incertitude de mesure donnée par le constructeur de l’appareil est de ±1 % sur la valeur
d’aération. La Figure 4-8 montre trois essais identiques. Pour les mêmes raisons que pour la
mesure des pertes de puissance, les résultats exposés dans la suite du chapitre sont définis en
moyennant trois essais.
16
14
12
Aération [%]
10
8
6
4
2
0
30 40 50 60 70 80 90 100
Température [C°]
Figure 4-8 Aérations mesurées pour trois essais identiques : Roue G3, Huile 1, 7000 tr/min, h/Rp=0,25
101
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
4.3 Réduction des pertes par barbotage
Plusieurs voies existent pour diminuer les pertes par barbotage. Il est possible de jouer sur la
géométrie des dentures, cependant celle-ci est largement contrainte par le cahier des charges. Il
est aussi possible de sélectionner un lubrifiant ayant une moindre propension à s’aérer, cependant
cela complexifie l’élaboration de celui-ci. Une troisième idée est de modifier les écoulements
d’huile dans le carter en modifiant sa forme ou en insérant des obstacles à l’intérieur de celui-ci.
Déflecteurs
Flasques
Flasques
Ja Jr Ja Roue dentée
Hf Déflecteur
Dp
Jr
102
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
L’étude [95] a été réalisée en utilisant les roues à denture droite de module 1,5, 3, et 5 mm,
pour des vitesses entre 1000 et 6000 tr/min, et pour une immersion de h/Rp=0,5. Les flasques
étaient placées symétriquement entre 1 et 10mm de la roue, les déflecteurs entre 5 et 10mm du
cylindre de tête des dents. Trois huiles de transmission mécanique ont été testées. Les essais se
sont concentrés sur des températures entre 20 et 40°C, les aspects haute température, faible
viscosité, et aération n’ont pas été investigués.
Dans l’article aucune configuration utilisant un seul flasque n’a été analysée, les auteurs ont
cependant testé une configuration avec des jeux différents (un flasque proche à 1mm et un
éloigné à 10mm). Ils mesurent peu de différence par rapport à deux flasques proches. Dans la
suite du chapitre les configurations symétriques feront appel à l’utilisation de deux flasques à des
jeux identiques. Les configurations asymétriques regroupent l’utilisation de jeux différents ou plus
simplement l’utilisation d’un seul flasque, lorsqu’un seul jeu est donné un seul flasque est utilisé.
La Figure 4-11 met en avant un résultat important de l’étude de 2008. Elle montre le rapport
entre Qref la perte de puissance sans flasque et Q la perte mesurée dans le cas d’utilisation
d’obstacles. L’huile est éjectée radialement par rapport à la roue, les déflecteurs se trouvent donc
en aval de l’écoulement. Les déflecteurs ont peu d’influence sur les pertes de puissances. Les
flasques ont par contre un très grand impact sur les pertes. Le gain est d’autant plus grand que le
jeu est faible. Les flasques diminuent fortement l’effet de pompage lié à la denture, pour cela il est
nécessaire que leur hauteur soit supérieure au diamètre de tête de la roue dentée.
Un gain de plus de 50% sur les pertes par barbotage semble atteignable en utilisant des jeux
axiaux de l’ordre de quelques millimètres. Les flasques sont plans et peuvent être réalisés en tôle.
Cette solution est réalisable à bas coût. Cependant des investigations plus poussées sont
nécessaires pour analyser la viabilité de la solution pour un réducteur haute vitesse.
Figure 4-11 Influence du jeu axial et radial sur les pertes de puissance, Roue G3, h/Rp=0,5 [95]
103
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
collés au carter, un sur le côté gauche (2) et un autre sur le fond (1) où une des résistances
chauffantes est collée. Les thermocouples (3) et (4) sont collés, sur les flasques, suffisamment
haut pour ne pas être immergés dans le bain d’huile. Enfin un dernier (5) est suspendu à 10mm
au-dessus du fond.
4
3
1
1
5
1 1
140
120
100
Température [°C]
1 Bas
80 2 Coté
60 3 Flasque devant
4 Flasque derrière
40
5 Huile
20
6 Air ambient
0
0 20 40 60 80 100
Temps [min]
Figure 4-13 Températures mesurées lors d’un essai, Roue G3, Huile 1, 7000 tr/min, h/Rp=0,25, Jeux axiaux
symétriques à 4 mm
Si les roulements ne sont pas partiellement immergés (ce qui est le cas dans le réducteur HV),
ils doivent être lubrifiés par la récupération des projections d’huile générées par les pignons. Afin
d’investiguer l’impact des flasques, des mesures de débit ont été réalisées (voir Tableau 4-3). La
Figure 4-14 montre une partie des résultats. Les flasques ne coupent pas le débit de projection
mesuré. Les débits mesurés avec deux flasques sont du même ordre que ceux mesurés sans
104
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
flasque et ce même en configuration asymétrique. Dans le cas où le jeu est grand et les flasques
placées symétriquement, les projections sont guidées et il est possible de récolter plus d’huile.
Tableau 4-3 Conditions pour les mesures de débit réalisées
Roues G3-G3B-G1,5
h/Rp 0,3-0,5
Jeux axiaux Symétrique à 4, 10mm
Asymétrique à 1, 4, 7mm
Température 25, 40, 60°C
Vitesse 1000 à 7000 tr/min
Huile 3
H 20mm
L 130mm
250
200
Débit en [mL/s]
150
100
50
0
1000 3000 5000 7000
Vitesse [tr/min]
Figure 4-14 Débits mesurés pour plusieurs configurations, Roue G3, h/Rp= 0,3, Huile 3 à 40°C
105
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Les flasques limitent les pertes de puissance en gênant l’effet de pompe centrifuge lié à la
denture [95]. Cependant elles limitent aussi le volume d’huile brassé par la roue, l’huile et l’air se
trouvant dans l’espace inter-dent se mélangent moins. Moins d’air est inséré dans le bain d’huile
et l’aération est plus faible.
20
18
16
Aération [%]
14
12
10
8
6
4
2
0
20 30 40 50 60 70 80 90 100
Température [°C]
Huile 1, Sans flasque Huile 2, Sans flasque
Huile 1, Symétrique 4 mm Huile 2, Symétrique 4 mm
Figure 4-15 Aération pour deux huiles différentes, avec et sans flasques, Roue G3, h/Rp=0,25, 7000tr/min
600
Perte par barbotage [W]
500
400
300
200
100
0
20 30 40 50 60 70 80 90 100
Température [°C]
Huile 1, Sans flasque Huile 2, Sans flasque
Huile 1, Symétrique 4 mm Huile 2, Symétrique 4 mm
Figure 4-16 Pertes de puissance pour deux huiles différentes, avec et sans flasques, Roue G3, h/Rp=0,25,
7000tr/min
Les figures Figure 4-17 et Figure 4-18 présentent le gain apporté par l’utilisation de flasques
pour une roue de module plus faible ainsi qu’à une immersion plus élevée. Elles montrent le
rapport entre les mesures sans flasque prisent comme référence (notées Qref pour la perte par
barbotage et A%ref pour l’aération) et les mesures réalisées avec les flasques (Q et A%). Les
flasques diminuent encore une fois largement l’aération jusqu’à des niveaux où il n’est pas
nécessaire de la prendre en compte (A%<10%).
Les pertes par barbotage sont elles aussi réduites, cependant le gain diminue avec
l’augmentation de la température alors que les niveaux d’aération mesurés sont suffisamment
faibles pour être négligés. Cette diminution du gain est sans doute reliée au remplissage de
l’espace inter flasque. Cet espace est d’autant plus facile à remplir que l’huile est chaude (et donc
fluide) et que l’immersion relative est élevée. L’efficacité des flasques semble découler d’un
équilibre entre remplissage gravitaire de l’espace inter flasque et l’éjection d’huile provoquée par
la denture.
106
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
1
0,9
0,8
A%/A%ref 0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
20 40 60 80 100
Température [°C]
Huile 1 Huile 2
Figure 4-17 Gain d’aération pour deux huiles différentes, Roue G1.5, h/Rp=0,5, 7000tr/min, configuration
symétrique 4mm
1
0,9
0,8
0,7
0,6
Q/Qref
0,5
0,4
0,3
0,2
0,1
0
20 40 60 80 100
Température [°C]
Huile 1 Huile 2
Figure 4-18 Gain de pertes de puissance pour deux huiles différentes, Roue G1.5, h/Rp=0,5, 7000tr/min,
configuration symétrique 4mm
Les huiles 1 et 2 se distinguent par leur propension à s’aérer. En utilisant des flasques
l’aération est diminuée. Le choix du lubrifiant a alors moins d’impact sur les pertes par barbotage
lorsque des flasques sont utilisés.
4.3.3.2 Influence du jeu axial sur l’aération
En éloignant les flasques de la roue le gain sur les pertes diminue. L’espace inter flasque est
plus grand et se remplit plus facilement. Les figures Figure 4-19 et Figure 4-20 comparent les
gains en aération et pertes pour deux jeux différents. A basse température un jeu de 7mm permet
de réaliser un gain de pertes de l’ordre de 50%. La température augmentant, l’espace inter
flasques se remplit et le gain diminue. A haute température les flasques placées à 7mm permettent
un gain faible. Le gain de pertes à haute température est de l’ordre de 15% et correspond au gain
réalisé par la diminution de l’aération.
L’utilisation de flasques placés symétriquement par rapport à la roue permet de réduire
l’aération à des niveaux acceptables (moins de 10%). Les pertes sont doublement impactées.
D’une part l’aération diminue et n’est plus un problème à haute température et haute vitesse.
D’autre part les flasques, en changeant l’écoulement autour de la roue, permettent une large
diminution des pertes de puissance. Cependant ce deuxième effet est dépendant de la
température et de la capacité de lubrifiant à remplir l’espace inter flasque.
107
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
1
0,9
0,8
0,7
A%/A%ref
0,6
0,5
0,4
0,3
0,2
0,1
0
20 30 40 50 60 70 80 90 100
Température [°C]
Symétrique 4 mm Symétrique 7 mm
Figure 4-19 Gain d’aération pour deux jeux différents, Roue G1.5, h/Rp=0,5, 7000tr/min, Huile 1
0,8
0,6
Q/Qref
0,4
0,2
0
20 30 40 50 60 70 80 90 100
Température [°C]
Symétrique 4 mm Symétrique 7 mm
Figure 4-20 Gain de pertes de puissance pour deux jeux différents, Roue G1.5, h/Rp=0,5, 7000tr/min, Huile 1
Les flasques diminuent largement le barbotage cependant elles perdent en efficacité avec la
montée en température. Pour maintenir le gain sur la plus grande plage d’utilisation il est
nécessaire de diminuer le jeu à une valeur inférieure à 4 mm.
4.3.3.1 Utilisation d’un unique flasque
Il n’est pas toujours possible d’insérer un flasque de chaque côté des roues. Les figures
Figure 4-21 et Figure 4-22 montrent la comparaison pour un cas très immergé entre les gains
obtenus avec une configuration symétrique et avec une asymétrique. L’utilisation d’un seul
flasque est beaucoup moins efficace pour réduire l’aération et les pertes. Pour ce cas très
immergé, seule la moitié du gain potentiel est atteinte.
Les flasques ne semblent pas interagir fortement entre elles. Le gain sur les pertes et
l’aération vient de la modification des écoulements sur le côté de la roue et non pas d’un
asséchement de volume piégé entre les deux flasques.
108
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
1
0,9
0,8
A%/A%ref 0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
20 30 40 50 60 70 80 90 100
Température [°C]
Asymétrique 4mm Symétrique 4 mm
Figure 4-21 Gain d’aération pour deux configurations de flasque différentes, Roue G1.5, h/Rp=0,5, 7000tr/min,
Huile 1
1
0,9
0,8
0,7
0,6
Q/Qref
0,5
0,4
0,3
0,2
0,1
0
20 30 40 50 60 70 80 90 100
Température [°C]
Asymétrique 4mm Symétrique 4 mm
Figure 4-22 Gain de pertes de puissance pour deux configurations de flasque différentes, Roue G1.5, h/Rp=0,5,
7000tr/min, Huile 1
Les figures Figure 4-23 et Figure 4-24 montrent le même type de mesures mais réalisées avec
une immersion plus faible ainsi qu’avec une roue de module plus élevé provoquant plus
d’aération. Les mêmes tendances que celles relevées précédemment à l’aide de la roue de module
1,5 sont présentes. Le gain obtenu avec un seul flasque est à mi-chemin entre celui réalisé avec
deux flasques et le cas sans flasque.
L’utilisation d’un seul flasque se révèle donc être moins intéressante pour réduire les pertes
de puissance. Il est donc nécessaire de privilégier l’utilisation de deux flasques pour la roue du
réducteur HV.
109
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
1
0,9
0,8
0,7
A%/A%ref
0,6
0,5
0,4
0,3
0,2
0,1
0
20 30 40 50 60 70 80 90 100
Température [°C]
Symétrique 4 mm Asymétrique 4mm
Figure 4-23 Gain d’aération pour deux configurations de flasque différentes, Roue G3, h/Rp=0,2, 7000tr/min, Huile
1
1
0,9
0,8
0,7
0,6
Q/Qref
0,5
0,4
0,3
0,2
0,1
0
20 30 40 50 60 70 80 90 100
Température [°C]
Symétrique 4 mm Asymétrique 4mm
Figure 4-24 Gain de pertes de puissance pour deux configurations de flasque différentes, Roue G3, h/Rp=0,2,
7000tr/min, Huile 1
110
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
20
18
16
14
Aération [%]
12
10
8
6
4
2
0
20 30 40 50 60 70 80
Température [°C]
Flasque devant Flasque derrière
Figure 4-25 Pertes de puissance selon la position de la flasque, Roue G3b, h/Rp=0,3, 7000tr/min, Huile 1
600
Perte par barbotage [W]
500
400
300
200
100
0
20 30 40 50 60 70 80
Température [°C]
Flasque devant Flasque derrière
Figure 4-26 Aération selon la position de la flasque, Roue G3b, h/Rp=0,3, 7000tr/min, Huile 1
0,8
0,6
Q/Qref
0,4
0,2
0
20 30 40 50 60 70 80 90 100
Température [°C]
Figure 4-27 Gain de pertes de puissance, Roue G3b, h/Rp=0,3, 7000tr/min, Huile 1, configuration symétrique 4
mm
111
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
4.4 Développement d’un modèle
0,8
0,7
0,6
0,5
Q/Qref
0,4
0,3
0,2
0,1
0
0,2 0,4 0,6
h/Rp
Figure 4-28 Gain de perte de puissance apporté par des flasques positionnés à 4mm, en fonction du niveau
d'immersion, Roue G3, huile 1, 30°C, 7000tr/min
Partant du constat que les formulations sur le gain apporté par les flasques proposées en
2007 peuvent être prises en défaut, ce sous chapitre porte sur l’élaboration d’une nouvelle
approche. Cette approche se concentre sur le régime 𝛾 > 1250 𝑚/𝑠² et 𝑅𝑒𝑐 > 4000, régime
qui concerne les hautes vitesses et des températures plus élevées. L’utilisation d’un unique flasque
s’est révélé être moins efficace. La nouvelle formulation se concentrera sur les configurations où
les flasques sont placées symétriquement. L’approche se place dans la continuité du travail exposé
dans le paragraphe 1.5.1.
112
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
Les flasques ont un double effet de réduction des pertes de puissance. D’une part elles
empêchent le bain d’huile de s’aérer. D’autre part en modifiant l’écoulement d’huile elles
diminuent le brassage du lubrifiant. Pour les conditions mêlant grande vitesse et haute
température le rapport Q/Qref comprend donc des facteurs liés à l’aération du bain d’huile.
L’aération sans flasque n’étant pas facilement caractérisable, la définition d’un nouveau facteur
Q/Qref semble complexe.
Une autre façon de faire pour caractériser les pertes de puissance en présence de flasques est
de repartir de l’analyse dimensionnelle du problème de barbotage. Celle-ci a permis d’exprimer le
couple adimensionnel de traînée sous la forme de cinq groupes de paramètres qui utilisent le
diamètre primitif comme paramètre principal associé à la longueur [54] :
ℎ 𝑉0 𝑏
𝐶𝑚 = 𝑓( , 3 , 𝐹𝑟, 𝑅𝑒, ) 4.5.
𝐷𝑝 𝐷𝑝 𝐷𝑝
L’exploitation des essais et l’analyse des phénomènes physiques en jeux permettent ensuite de
définir :
quelles longueurs caractéristiques doivent être prises en compte dans les nombres
adimensionnels
ainsi que la forme de la fonction f :
ℎ 𝜓2 𝑉0 𝜓3 𝑏 𝜓6
𝐶𝑚 = 𝜓1 ∗ ( ) ( 3 ) 𝐹𝑟 𝜓4 𝑅𝑒𝑐 𝜓5 ( ) 4.6.
𝐷𝑝 𝐷𝑝 𝐷𝑝
Les coefficients 𝜓 sont des constantes qui peuvent être fixées grâce aux essais. Par exemple
en réalisant un essai en ne faisant varier que la température, l’évolution de Cm sur un couple
d’essais (a,b) est alors uniquement dépendante de la variation du nombre de Reynolds. Il est alors
possible d’exprimer le coefficient 𝜓5 :
𝐶𝑚
ln(𝐶𝑚𝑎 )
𝑏
𝜓5 = 4.7.
𝑅𝑒
ln(𝑅𝑒𝑎 )
𝑏
Il est ensuite possible de déterminer l’exposant du nombre de Froude en réalisant deux essais
(a,b) où seule la vitesse change. La vitesse intervenant aussi dans l’expression du nombre de
Reynolds il est nécessaire de le prendre son évolution en compte :
𝐶𝑚 𝑅𝑒 𝜓5
ln(𝐶𝑚𝑎 . (𝑅𝑒𝑏 ) )
𝑏 𝑎
𝜓4 = 4.8.
𝐹𝑟
ln(𝐹𝑟𝑎 )
𝑏
113
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
adimensionnel 𝐽𝑎/𝑏 ressort, celui-ci est utilisé pour faire ressortir explicitement que plus la
denture est étroite plus il faut rapprocher les flasques. De plus comme dans l’article de 2008, les
essais révèlent l’influence du rapport 𝐷𝑝 ⁄𝐻𝑓 dans l’exposant du paramètre prenant en compte
l’influence du jeu. Cela conduit à formuler le couple adimensionnel de trainée sous la forme
suivante :
ℎ 𝜓2 𝑉0 𝜓3 𝜓 𝜓5
𝑏 𝜓6 𝐻𝑓 𝜓7 𝐽𝑎 𝜓8 (𝐷𝑝⁄𝐻𝑓 )
𝐶𝑚 = 𝜓1 ∗ ( ) ( 3 ) 𝐹𝑟 𝑅𝑒𝑐 ( ) ∗ ( ) ( )
4
4.9.
𝐷𝑝 𝐷𝑝 𝐷𝑝 𝐷𝑝 𝑏
La Figure 4-29 montre le comportement typique des pertes constaté lors des essais. A partir
d’une certaine température les flasques perdent en efficacité. Jusqu’à cette température les pertes
ne dépendent de la température que par la dilatation du bain d’huile. Au-delà de cette
température, la diminution du gain obtenu sur les pertes est à relier à la viscosité de l’huile et non
à l’aération qui elle est maintenue à des niveaux faibles grâce aux flasques. Les pertes avec les
flasques augmentent et finissent par égaler les pertes par barbotage sans flasque, et sans aération.
La pente de ce régime de transition dépend du jeu axial, de la vitesse, de l’immersion, ainsi que de
l’huile. De manière générale aucun effet des flasques n’a été constaté pour des jeux supérieurs à
10mm. Les jeux qui permettent de maintenir le gain maximum sur une large plage de température
et de vitesse sont inférieurs à 3mm.
1000
Pertes par barbotage [W]
800
600
400
200
0
25 35 45 55 65
Température [°C]
Symétrique à 4 mm Sans Flasques
Figure 4-29 Comportement des pertes de puissance en présence de flasques, Roue G3, 7000tr/min, huile 3,
h/Rp=0,6
3,00E-04
2,50E-04
2,00E-04
Cm
1,50E-04
1,00E-04
5,00E-05
0,00E+00
45 35 25 15 5
Viscosité [cSt]
Figure 4-30 Comportement du Cm en fonction de la viscosité, Roue G3, 7000tr/min, huile 3, h/Rp=0,6
114
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
120
110
100
90
80
5000 7000
Vitesse [tr/min]
25°C 40°C 60°C
Figure 4-31 Débit d'huile en fonction de la température, Roue G3, Huile 3, un flasque à 4mm
115
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Le changement de régime d’écoulement est intimement lié aux phénomènes ayant lieu dans
l‘espace entre la roue dentée et les flasques. Un critère se concentrant sur ces phénomènes a donc
été créé.
Deux nombres adimensionnels sont utilisés. Le premier est le nombre de Froude. Le volume
entre les flasques et la roue dentée est le lieu de phénomènes en concurrence. La rotation de la
roue tend à éjecter l’huile et à vider ce volume. A l’inverse les forces de gravité ramènent l’huile et
tendent à remplir ce volume. Le nombre de Froude est alors formé avec la hauteur immergée et
non avec le rayon primitif comme dans la formule du couple adimensionnel de trainée.
Un second nombre adimensionnel est utilisé pour analyser l’écoulement d’huile qui prend
forme entre les flasques. Il découle d’un nombre de Reynolds formé avec comme longueur
caractéristique le jeu entre le flasque et la roue dentée :
𝜔𝑅𝑝𝐽𝑎
𝑅𝑒 = 4.11.
𝜈
𝐽𝑎
Cependant l’étude des résultats fait ressortir le rôle du rapport 𝑅𝑝 . Celui-ci est alors multiplié
au Reynolds pour faire ressortir un nombre de Taylor. A l’origine ce nombre sert à caractériser la
stabilité de l'écoulement d'un fluide placé entre deux tubes concentriques. Par analogie il est
appliqué au cas de l’espace plan entre le flasque et la roue dentée.
Le rapport des deux nombres adimensionnels constitue le critère créé pour juger de
l’efficacité des flasques :
𝐹𝑟
𝛤= 4.12.
𝑇𝑎
Où :
𝛺2 𝑅𝑝 2
-𝐹𝑟 = = énergie cinétique / énergie potentiel de gravité, nombre de Froude
𝑔ℎ
𝛺𝐽𝑎²
-Ta = = forces d’inertie / forces visqueuses, nombre de Taylor
𝜈
Les essais ont permis de fixer la valeur de ce critère à partir duquel la formule 4.10 est valide.
Ainsi, si 𝛤 est supérieure à 10 les flasques donnent leur plein potentiel. De plus le critère 𝛤 peut
être utilisé comme règle de conception lors du placement des flasques. Il permet de calculer le jeu
maximum à partir duquel les flasques seront moins efficaces.
116
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
117
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Tableau 4-5 Résultats de la nouvelle modélisation en fonction du numéro d’essai
350 350
300 300
250 250
200 200
150 150
100 100
50 50
0 0
1 2 3 4 5 6 7 8 9 10
350 350
300 300
Perte par barbotage [W]
Perte par barbotage [W]
250 250
200 200
150 150
100 100
50 50
0 0
11 12 13 14 15 16 17 18 19 20
350
300
Perte par barbotage [W]
250
200
150
100
50
0
21 22 23 24 25
118
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
4.5 Application
60
50
40
30 Route
𝛤
Autoroute
20
Borne inférieure
10
0
1 2 3 4
Jeu [mm]
Figure 4-32 Critère Γ en fonction du jeu entre flasques et roue, h/Rp=0,5, Huile à 100°C
119
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
90
Température [°C]
80
70
60
50
Bague extérieure Bague intérieure Pignon Huile
130
Température [°C]
120
110
100
90
80
Bague extérieure Bague intérieure Pignon Huile
120
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
90
Température [°C]
80
70
60
50
Bague extérieure Bague intérieure Pignon Huile
130
Température [°C]
120
110
100
90
80
Bague extérieure Bague intérieure Pignon Huile
121
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
4.5.3 Augmentation du rendement attendue sur le réducteur final
Dans le souhait de calculer le gain envisageable sur l’application finale, la démarche exposée
dans le sous-chapitre 2.3 est réutilisée. Le rendement est calculé en supposant une température
d’huile homogène dans le réducteur. Il n’y a pas de couplage thermomécanique.
La Figure 4-35 compare le rendement des deux architectures présentées dans le sous-
chapitre 2.2. En posant l’hypothèse que le bain d’huile est susceptible de s’aérer jusqu’à 25%, les
deux architectures ont des rendements comparables. L’architecture à 3 trains est alors la plus
intéressante car plus compacte.
L’utilisation de flasques permet d’éviter le problème de l’aération ainsi que de diminuer
largement les pertes par barbotage de la roue de l’étage GV. Le rendement des deux architectures
peut être amélioré de plus de 0,4 point.
L’architecture à deux trains est initialement pénalisée par le barbotage à cause de la taille de
ses roues. Provoquant de plus grandes pertes par barbotage, elle bénéficie par contre plus de
l’ajout de flasques. Au final, elle permet ainsi d’obtenir un rendement supérieur à l’architecture à
trois trains. Elle reste cependant la plus encombrante et la plus lourde.
3 Trains Autoroute
2 Trains Autoroute
Aération 25%
3 Trains Route
Flasques 2mm
2 Trains Route
90 91 92 93 94 95 96
Rendement [%]
Figure 4-35 Rendement des architectures 2 et 3 trains avec et sans flasques sur la roue GV
4.6 Conclusion
A grande vitesse, la roue de l’étage GV provoque des pertes par barbotage qui dégradent le
rendement de la transmission. A des températures supérieures à 80°C, elle provoque l’aération du
bain d’huile et encore plus de puissance est perdue. Une solution simple a été identifiée pour
réduire les pertes de puissance. La littérature montre qu’en insérant des flasques parallèles aux
roues et suffisamment proches de ses flancs les pertes par barbotage peuvent être réduites.
Pour investiguer cette solution un banc d’essais dédié a été utilisé. Il permet de mesurer les
pertes par barbotage provoquées par une roue partiellement immergée. Deux systèmes de mesure
supplémentaires viennent équiper ce banc. Le premier est un débitmètre venant réceptionner les
projections d’huile créées par la roue dentée. Le second permet la mesure de l’aération du bain
d’huile in situ.
122
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Etude du barbotage
Pour s’assurer que les flasques ne vont pas gêner le refroidissement du réducteur, des
mesures de températures et de débits ont été réalisées. Les flasques ne gênent pas la circulation de
l’huile dans le bain, de plus, elles n’impactent que faiblement les projections d’huile.
Des mesures de pertes par barbotage et d’aération ont été réalisées avec une vitesse en tête de
denture comparable avec celles du réducteur HV. Grâce à des balayages en température de 30 à
100°C l’utilisation de flasques a été investiguée. Utiliser deux flasques permet de réduire l’aération
à des niveaux inférieurs à 10%. A ce niveau il n’est plus nécessaire de la prendre en compte. De
plus une très grande réduction des pertes par barbotage a été mesurée sur toute la plage de
température. Ces tendances sont aussi valides pour les engrenages hélicoïdaux. Utiliser un seul
flasque a moins d’impact sur les pertes et l’aération. Un gain est réalisable mais celui-ci est
beaucoup plus faible qu’avec deux flasques.
L’étude met en avant deux régimes d’écoulement liés aux flasques. Le premier ne dépend pas
de la température et offre un gain maximal. Ce régime a été caractérisé et une nouvelle
formulation permet de calculer les pertes. Le second régime correspond à une perte d’efficacité
des flasques lorsque la température augmente. L’espace entre la roue dentée et le flasque se
remplit peu à peu et les flasques finissent par ne plus avoir d’impact sur les pertes. En plaçant les
flasques suffisamment proches il est possible d’éviter ce phénomène. Un critère a été créé pour
s’assurer que la nouvelle formulation des pertes est applicable. L’utilisation du critère permet de
dimensionner le jeu entre les flasques et la roue dentée, et de donner un jeu maximal pour des
conditions opératoires.
Le critère appliqué au réducteur HV permet de justifier un jeu de l’ordre de 2mm entre la
roue GV et les flasques. En utilisant le modèle thermomécanique exposé dans le chapitre 3,
l’impact des flasques sur les températures a été analysé. En réduisant les pertes de puissance, les
flasques diminuent globalement les températures de la transmission. Elles permettent aussi
d’augmenter le niveau d’huile dans le réducteur en ne pénalisant que faiblement le rendement. Un
niveau plus grand permet d’augmenter le débit d’huile amené aux roulements et de mieux les
refroidir.
Enfin, les calculs isothermes de rendement réalisés pour comparer deux architectures dans le
chapitre 2 ont été mis à jour. Utiliser des flasques sur la roue de l’étage GV permet de gagner plus
de 0,4 point de rendement. Les architectures deux trains et trois trains ont des rendements
comparables si l’huile s’aère. Avec des flasques l’architecture à deux trains n’est plus pénalisée par
le barbotage et montre un rendement supérieur à l’autre modèle de réducteur. Il reste cependant
plus lourd et plus encombrant.
123
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Conclusion générale
126
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
pertes amplifiées par l’aération du bain d’huile. Cela contribue à l’échauffement général du
réducteur et à la dégradation de son rendement.
Le troisième chapitre propose une amélioration de la modélisation thermomécanique des
roulements de l’arbre GV. En l’absence de mesure de pertes de puissance sur les roulements GV,
un nouveau modèle est développé. Il se base sur une description locale des sources de chaleurs.
Elle vient alimenter un réseau thermique décrivant un roulement par plusieurs températures. Ce
nouveau modèle thermomécanique de roulements est capable de reproduire les différences de
température mesurées entre les bagues des roulements d’un banc d’essai FZG. Appliquée au
modèle de l’étage GV, présenté dans le chapitre deux, cette modélisation des roulements conduit
à des températures de fonctionnement des pièces tournantes plus élevées. Les gradients de
température dans les roulements doivent être pris en compte.
Grâce à ce modèle, deux études de sensibilité sont menées. Il est démontré qu’il n’y a pas
d’intérêt thermique à avoir un niveau d’huile élevé. Même si les pignons sont mieux refroidis le
gain n’est pas suffisant pour compenser la hausse de température liée à la dégradation du
rendement. Il est aussi montré qu’il est intéressant d’amener beaucoup d’huile aux roulements.
Les roulements sont mieux refroidis et la baisse de température se répercute directement sur les
températures dentures. Ceci ne peut se faire qu’en gardant une immersion significative de la roue,
la création de pertes par barbotage importantes semble donc inévitable.
Le quatrième et dernier chapitre porte sur l’étude des pertes par barbotage et sur leur
diminution. Le travail a été réalisé à l’aide d’un banc d’essais dédié à l’étude du barbotage. Ce
moyen d’essais permet de s’approcher des vitesses périphériques visées sur l’étage GV. Il permet
au choix de mesurer le débit d’huile projeté par la denture ou l’aération du bain d’huile. La
littérature montre qu’en insérant des flasques parallèles aux pignons et suffisamment proches de
ses flancs les pertes par barbotage peuvent être réduites. Après s’être assuré que cette solution ne
nuira pas au refroidissement du réducteur, son impact sur les pertes et sur l’aération est investigué
en détail. Utiliser deux flasques permet de réduire les pertes sur de grandes plages de température.
Elles réduisent aussi fortement l’aération au point que la prise en compte de celle-ci ne soit plus
nécessaire. Il n’est alors plus nécessaire d’optimiser le lubrifiant sur des caractéristiques limitant
l’aération. L’utilisation d’un seul flasque est par contre moins intéressante. Un travail
expérimental a permis de caractériser les pertes de puissance à grande vitesse lorsque deux
flasques sont utilisés. Cependant les flasques perdent en efficacité lorsque : la viscosité de l’huile
est trop faible, le jeu est trop grand, la vitesse est trop faible. Pour s’assurer que le gain maximum
est atteint un critère prenant en compte les phénomènes physiques a été créé. En utilisant ce
critère comme règle de conception il est possible d’identifier le jeu optimal.
L’utilisation du modèle thermomécanique exposé dans le chapitre 3 et de la caractérisation
des pertes avec des flasques montrent que : les températures de la transmission sont globalement
diminuées, un niveau d’huile plus haut peut être utilisé sans trop pénaliser le rendement, un
niveau d’huile haut est intéressant pour augmenter le débit d’huile amené aux roulements et ainsi
mieux refroidir paliers et pignons. En rebouclant sur le calcul du rendement des architectures
finales, les flasques permettent d’augmenter le rendement de plus de 0,4% sur route et autoroute.
L’architecture à deux trains se montre plus efficiente, elle reste cependant plus lourde et plus
encombrante que l’alternative à trois trains.
Les modèles analytiques proposées dans ce manuscrit permettent d’estimer le comportement
thermomécanique d’un réducteur haute vitesse. Pour quantifier de manière fiable les pertes à
l’engrènement il sera nécessaire de réaliser des essais de caractérisation du coefficient de
127
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
frottement. Les corrections de profil de denture peuvent largement améliorer le rendement.
Puisqu’elles ont aussi un grand impact sur la réponse vibratoire et sur la résistance mécanique des
engrenages, il parait être intéressant de les optimiser au regard de ces trois critères. Les pertes de
puissance provoquées par l’étanchéité sur l’arbre d’entrée n’ont pas fait l’objet d’une étude
poussée. D’une part les chicanes permettent de répondre au besoin de faible frottement mais
n’apportent pas d’étanchéité statique, d’autre part les joints à lèvre standards ne semblent pas
adaptés à la grande vitesse. Il parait nécessaire de travailler en collaboration avec un fabricant de
joints sur la définition de l’étanchéité de l’arbre GV.
L’étude a pointé le rôle majeur des roulements et de la lubrification sur le comportement
thermomécanique de premier étage d’un réducteur mécanique haute vitesse. Les pistes
d’améliorations identifiées tout au long de ce travail vont être appliquées à un futur banc d’essais
de réducteur haute vitesse (voir Figure C-1). Ainsi il est prévu d’insérer des flasques dans le carter
pour limiter l’aération et les pertes par barbotage. L’intérieur du carter sera conçu pour garder un
volume d’huile suffisamment grand. De plus la récupération des projections d’huile servant à la
lubrification des roulements sera maximisée.
128
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
précision du modèle et le nombre de données nécessaires. Pour cela il apparait indispensable de
réaliser des essais sur un banc d’essai dédié, pour isoler un roulement et l’étudier.
129
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Annexe 1 Caractérisation d’un
brouillard d’huile
𝜇𝑜𝑖𝑙
𝜇𝑒𝑞 = 𝜇
𝑜𝑖𝑙 (1 3.22.
𝜇𝑎𝑖𝑟 ∗ − 𝑋) + 𝑋
Avec
- 𝑋 fraction volumique d’huile dans l’air
- 𝜌 masse volumique [kg/m3]
- 𝑘 conductivité thermique [W/m/K]
- μ viscosité dynamique du fluide [Pa.s]
- 𝐶𝑝 chaleur massique [J/kg/K]
- Indices : - 𝑜𝑖𝑙 caractéristique de l’huile - 𝑎𝑖𝑟 caractéristique de l’air
- 𝑒𝑞 caractéristique du brouillard
132
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Annexe 2 Description du banc FZG
Denture type C
Pignon 16 [-]
Nombre de dents
Roue 24 [-]
Module normal 4.5 [mm]
Angle de pression normal 20 [°]
Largeur de denture 30 [mm]
Coefficient de Pignon 0.1817 [-]
déport Roue 0.1817 [-]
Pignon 72 [mm]
Diamètre primitif
Roue 108 [mm]
Entraxe de fonctionnement 91.6 [mm]
Rapport de conduite 1.449 [-]
La mesure des pertes de puissance provoquées par le banc FZG est réalisée avec un
couplemètre HBM type T22/50Nm, ayant une classe de précision de 0,5. Les mesures de
températures sur les parties fixes et de l’air ambient sont réalisées avec des thermocouples de type
K et des boitiers NI Picolog. Les mesures sur les parties tournantes sont obtenues par l’utilisation
de gauges de contraintes HBM TT-3/100 reliées à des collecteurs tournants HBM SK6. La mise
en œuvre du couplemètre et des collecteurs tournant est représentée sur la figure A2-1
a) b)
Figure A2-1 Accessoires utilisés sur le banc d’essais, a) Couplemètre b) Collecteurs tournants
134
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
Références bibliographiques
136
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
pp. 173–184.
[22] G. H. Benedict and B. W. Kelley, “Instantaneous coefficients of gear tooth friction,”
ASLE Trans., vol. 4, no. 1, pp. 59–70, 1961.
[23] “ISO/TR 14179-2: 2001. Gears – thermal capacity – part 2: thermal load-carrying
capacity.”
[24] Y. Diab, “Analyse des pertes de puissance dans les transmissions par engrenages à grande
vitesse. Applications aux réducteurs industriels et aux machines textiles,” INSA de LYON
- Institut National des Sciences Appliquées de Lyon, 2005.
[25] J. A. Greenwood and J. H. Tripp, “The contact of two nominally flat rough surfaces,”
Proc. Inst. Mech. Eng., vol. 185, no. 1, pp. 625–633, 1970.
[26] B. B. Mikić, “Thermal contact conductance; theoretical considerations,” Int. J. Heat Mass
Transf., vol. 17, no. 2, pp. 205–214, 1974.
[27] F. Ville, D. Nélias, G. Tourlonias, L. Flamand, and P. Sainsot, “On the two-disc machine:
A polyvalent and powerful tool to study fundamental and industrial problems related to
elastohydrodynamic lubrication,” Tribol. Ser., vol. 39, pp. 393–402, 2001.
[28] J. Durand de Gevigney, “Analyse thermomécanique d’un réducteur épicycloïdal:
Application aéronautique.” Lyon, INSA, 2013.
[29] T. A. Harris, Rolling Bearing Analysis, 3° Ed. Wiley Interscience, 1990.
[30] F. Pouly, C. Changenet, F. Ville, P. Velex, and B. Damiens, “Power Loss Predictions in
High-Speed Rolling Element Bearings Using Thermal Networks,” Tribol. Trans., vol. 53,
no. 6, pp. 957–967, Oct. 2010.
[31] “ISO/TR 14179-1:2001. Gears – thermal capacity – part 1: rating gear drives with thermal
equilibrium at 95° sump temperature.,” .
[32] Skf, Rolling bearings. 2013.
[33] C. M. C. G. Fernandes, P. M. T. Marques, R. C. Martins, and J. H. O. Seabra, “Gearbox
power loss. Part I: Losses in rolling bearings,” Tribol. Int., 2014.
[34] Y. Ariura, T. Ueno, T. SunAGA, and S. Sunamoto, “The lubricant churning loss in spur
gear systems,” Bull. JSME, vol. 16, no. 95, pp. 881–892, 1973.
[35] L. S. Akin, J. J. Mross, and D. P. Townsend, “Study of lubricant jet flow phenomena in
spur gears,” J. Tribol., vol. 97, no. 2, pp. 283–288, 1975.
[36] W. Mauz, “Hydraulische Verluste von Stirnradgetriebe bei Umfangsgeschwindigkeit bis 60
m/s,” IMK University of Stuttgart, 1987.
[37] N. E. Anderson and S. H. Loewenthal, “Effect of geometry and operating conditions on
spur gear system power loss,” J. Mech. Des., vol. 103, no. 1, pp. 151–159, 1981.
[38] P. H. Dawson, “Windage loss in larger high-speed gears,” Proc. Inst. Mech. Eng. Part A J.
Power Energy, vol. 198, no. 1, pp. 51–59, 1984.
[39] Y. Diab, F. Ville, P. Velex, and C. Changenet, “Windage Losses in High Speed Gears—
Preliminary Experimental and Theoretical Results,” J. Mech. Des., vol. 126, no. 5, p. 6,
2004.
[40] M. J. Hill, R. F. Kunz, R. W. Noack, L. N. Long, P. J. Morris, and R. F. Handschuh,
“Application and validation of unstructured overset CFD technology for rotorcraft
gearbox windage aerodynamics simulation,” in 64th Annual Forum of the American Helicopter
Society, 2008.
137
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
[41] Y. Marchesse, C. Changenet, F. Ville, and P. Velex, “Investigations on CFD Simulations
for Predicting Windage Power Losses in Spur Gears,” J. Mech. Des., vol. 133, no. 2, p.
024501, 2011.
[42] C. Changenet, X. Oviedo-Marlot, and P. Velex, “Power Loss Predictions in Geared
Transmissions Using Thermal Networks-Applications to a Six-Speed Manual Gearbox,” J.
Mech. Des., vol. 128, no. 3, p. 618, 2006.
[43] C. E. Douglas and A. Thite, “Effect of lubricant temperature and type on spur gear
efficiency in racing engine gearbox across full engine load and speed range,” Proc. Inst.
Mech. Eng. Part J J. Eng. Tribol., vol. 229, no. 9, pp. 1095–1113, Sep. 2015.
[44] J. W. Daily and R. E. Nece, “Chamber dimension effects on induced flow and frictional
resistance of enclosed rotating disks,” J. Fluids Eng., vol. 82, no. 1, pp. 217–230, 1960.
[45] R. W. Mann and C. H. Marston, “Friction Drag on Bladed Disks in Housings as a
Function of Reynolds Number, Axial and Radial Clearance, and Blade Aspect Ratio and
Solidity,” J. Fluids Eng., vol. 83, no. 4, pp. 719–723, 1961.
[46] S. L. Soo and N. J. Princeton, “Laminar flow over an enclosed rotating disk,” Trans.
ASME, vol. 80, pp. 287–296, 1958.
[47] R. J. Boness, “Churning losses of discs and gears running partially submerged in oil,” in
Proc. ASME Int. Power Trans. Gearing Conf, 1989, vol. 1, pp. 255–359.
[48] E. Lauster and M. Boos, “Zum Wärmehaushalt Mechanischer Schaltgetriebe für
Nutzfahrzeuge,” VDI. Ber., no. 488, pp. 45–55, 1983.
[49] A. S. Terekhov, G. Shpitko, and S. Gulezov, “Thermal calculation of transmission
assemblies,” in World congress on gearing and power transmission, 1999, pp. 2247–2257.
[50] A. S. Terekhov, “Hydraulic losses in gearboxes with oil immersion,” Russ. Eng. J., vol. 55,
no. 5, pp. 7–11, 1975.
[51] S. Seetharaman and A. Kahraman, “Load-independent spin power losses of a spur gear
pair: model formulation,” J. Tribol., vol. 131, no. 2, p. 22201, 2009.
[52] S. Seetharaman, A. Kahraman, M. D. Moorhead, and T. T. Petry-Johnson, “Oil churning
power losses of a gear pair: experiments and model validation,” J. Tribol., vol. 131, no. 2, p.
22202, 2009.
[53] C. Changenet and P. Velex, “A Model for the Prediction of Churning Losses in Geared
Transmissions—Preliminary Results,” J. Mech. Des., vol. 129, no. 1, p. 128, 2007.
[54] C. Changenet, G. Leprince, F. Ville, and P. Velex, “A Note on Flow Regimes and
Churning Loss Modeling,” J. Mech. Des., vol. 133, no. 12, p. 121009, 2011.
[55] P. M. T. Marques, C. M. C. G. Fernandes, R. C. Martins, and J. H. O. Seabra, “Power
losses at low speed in a gearbox lubricated with wind turbine gear oils with special focus
on churning losses,” Tribol. Int., vol. 62, pp. 186–197, 2013.
[56] P. Luke and A. V Olver, “A study of churning losses in dip-lubricated spur gears,” Proc.
Inst. Mech. Eng. Part G J. Aerosp. Eng., vol. 213, no. 5, pp. 337–346, 1999.
[57] G. LePrince, “Pertes mécaniques par frottement et lubrification dans une boîte de
vitesses,” INSA de LYON - Institut National des Sciences Appliquées de Lyon, 2011.
[58] G. Leprince, C. Changenet, F. Ville, P. Velex, and F. Jarnias, “Influence of oil aeration on
churning losses,” in JSME 2009 International Motion and Power Transmission Conference, 2009,
pp. 463–468.
[59] T. I. Fowle, “Aeration in lubricating oils,” Tribol. Int., vol. 14, no. 3, pp. 151–157, 1981.
138
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
[60] J. L. Nikolajsen, “The Effect of Aerated Oil on the Load Capacity of a Plain Journal
Bearing©,” Tribol. Trans., vol. 42, no. 1, pp. 58–62, 1999.
[61] S. M. Chun, “Aeration effects on the performance of a turbocharger journal bearing,”
Tribol. Int., vol. 41, no. 4, pp. 296–306, 2008.
[62] P. W. Centers, “Behavior of silicone antifoam additives in synthetic ester lubricants,”
Tribol. Trans., vol. 36, no. 3, pp. 381–386, 1993.
[63] L. T. Dixon and S. Korcek, “Foaming and air entrainment in automatic transmission
fluids,” SAE Technical Paper, 1976.
[64] M. Duncanson, “Controlling Oil Aeration and Foam,” Pract. Oil Anal. Mag., vol. 11, 2001.
[65] Simrit: Radialwellendichtringe, Katalog Nr. 100. .
[66] H. Linke, “Stirnradverzahnung,” in Stirnradverzahnung, Carl Hanser Verlag GmbH & Co.
KG, 2010, pp. 1–14.
[67] R. C. Martins, P. S. Moura, and J. O. Seabra, “MoS2/Ti low-friction coating for gears,”
Tribol. Int., vol. 39, no. 12, pp. 1686–1697, 2006.
[68] R. Martins, J. Seabra, and L. Magalhães, “Austempered ductile iron (ADI) gears: Power
loss, pitting and micropitting,” Wear, vol. 264, no. 9–10, pp. 838–849, 2008.
[69] J. R. Brown and N. H. Forster, “Carbon phenolic cages for high speed bearings. Part III -
Development of numerical models for heat generation and temperature prediction in
lightly lubricated bearings,” Air Force Res. Lab. Rep., pp. 1–23, 2003.
[70] N. Patir and H. S. Cheng, “Prediction of the Bulk Temperature in Spur Gears Based on
Finite Element Temperature Analysis,” A S L E Trans., vol. 22, no. 1, pp. 25–36, Jan.
1979.
[71] J. Durand de Gevigney, C. Changenet, F. Ville, and P. Velex, “Thermal modelling of a
back-to-back gearbox test machine: Application to the FZG test rig,” Proc. Inst. Mech. Eng.
Part J J. Eng. Tribol., vol. 226, no. 6, pp. 501–515, Jan. 2012.
[72] J. P. Holman, Heat transfer, 7th ed. New York: McGraw-Hill Book Company, 1989.
[73] C. Changenet, “Modélisation du comportement thermique des transmissions par
engrenages,” INSA de LYON - Institut National des Sciences Appliquées de Lyon, 2006.
[74] G. LePrince, C. Changenet, F. Ville, P. Velex, C. Dufau, and F. Jarnias, “Influence of
Aerated Lubricants on Gear Churning Losses–An Engineering Model,” Tribol. Trans., vol.
54, no. 6, pp. 929–938, Nov. 2011.
[75] G. Leprince, C. Changenet, F. Ville, and P. Velex, “Investigations on Oil Flow Rates
Projected on the Casing Walls by Splashed Lubricated Gears,” Adv. Tribol., vol. 2012, pp.
1–7, 2012.
[76] H. Blok, “The flash temperature concept.,” J. Wear, vol. 6, pp. 483–493, 1963.
[77] C. Wagner, “Heat Transfer from a Rotating Disk to Ambient Air,” J. Appl. Phys., vol. 19,
no. 9, p. 837, 1948.
[78] H. Winter, K. Michaelis, and G. Funck, “Wärmeabführung bei Getrieben unter
Quasistationären Betriebsbedingungen. Teil II: Untersuchungen zur Wärmeabführung
über Stahlfundamente und Übertragung der Prüfstandsergebnisse auf die Praxis,”
Antriebstechnik, vol. 26, no. 6, pp. 49–55, 1987.
[79] J. L. Tevaarwerk and K. L. Johnson, “The Influence of Fluid Rheology on the
Performance of Traction Drives,” Journal of Lubrication Technology, vol. 101. p. 266, 1979.
139
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
[80] B. J. Hamrock and W. J. Anderson, Rolling-element bearings, vol. 1105. National Aeronautics
and Space Administration, 1983.
[81] P. K. Gupta, Advanced Dynamics of Rolling Elements. New York, NY: Springer New York,
1984.
[82] D. Nelias, P. Sainsot, and L. Flamand, “Power Loss of Gearbox Ball Bearing Under Axial
and Radial Loads©,” Tribol. Trans., vol. 37, no. 1, pp. 83–90, Jan. 1994.
[83] F. Pouly, “Modélisation thermo mécanique d’un roulement à billes grande vitesse,” INSA
de LYON - Institut National des Sciences Appliquées de Lyon, 2010.
[84] Y. Marchesse, C. Changenet, and F. Ville, “Numerical investigations on drag coefficient of
balls in rolling element bearing,” Tribol. Trans., vol. 57, pp. 778–785, 2014.
[85] A. Neurouth, C. Changenet, F. Ville, and A. Arnaudon, “Thermal modeling of a grease
lubricated thrust ball bearing,” Proc. Inst. Mech. Eng. Part J J. Eng. Tribol., Mar. 2014.
[86] J. Padet, “Convection thermique et massique - Nombre de Nusselt: partie 1,” Techniques de
l’ingénieur. 2005.
[87] Y. S. Muzychka and M. M. Yovanovich, “Thermal Resistance Models for Non-Circular
Moving Heat Sources on a Half Space,” J. Heat Transfer, vol. 123, no. 4, p. 624, 2001.
[88] C. M. C. G. Fernandes, P. M. P. Amaro, R. C. Martins, and J. H. O. Seabra, “Torque loss
in cylindrical roller thrust bearings lubricated with wind turbine gear oils at constant
temperature,” Tribol. Int., vol. 67, pp. 72–80, 2013.
[89] A. Neurouth, C. Changenet, F. Ville, and M. Octrue, “Influence de la modélisation des
roulements à billes sur la prédiction du comportement thermique d’un banc FZG,” Congrès
Français de Mécanique, 2015.
[90] K. Michaelis, “Testing procedures for gear lubricants with the fzg test rig,” Industrial
Lubrication and Tribology, vol. 26. pp. 91–94, 1974.
[91] C. M. C. G. Fernandes, R. C. Martins, and J. H. O. Seabra, “Torque loss of type C40 FZG
gears lubricated with wind turbine gear oils,” Tribol. Int., vol. 70, pp. 83–93, 2014.
[92] A. Doleschel, “Method to determine the frictional behavior of gear lubricants using fzg
gear test rig.” FVA Information Sheet n°345, 2002.
[93] R. C. Martins, N. F. R. Cardoso, H. Bock, A. Igartua, and J. H. O. Seabra, “Power loss
performance of high pressure nitrided steel gears,” Tribol. Int., vol. 42, no. 11–12, pp.
1807–1815, 2009.
[94] A. Neurouth, C. Changenet, F. Ville, and M. Octrue, “Influence of ball bearings modeling
on the predicted thermal behavior of a high-speed gear unit,” STLE 70th Annu. Meet.
Exhib., 2015.
[95] C. Changenet and P. Velex, “Housing Influence on Churning Losses in Geared
Transmissions,” J. Mech. Des., vol. 130, no. 6, p. 062603, 2008.
[96] B. Deconninck, T. Delvigne, and G. Videx, “Air-X, an Innovative Device for On-line Oil
Aeration Measurement in Running Engines,” 2003.
[97] A. Neurouth, C. Changenet, F. Ville, and M. Octrue, “Is splash lubrication compatible
with efficient gear units for high-speed applications?,” Int. Gear Conf. Lyon, Conf. Proc., vol.
2, pp. 1060–1068, 2014.
[98] H. S. Isbin, N. C. Sher, and K. C. Eddy, “Void fractions in two‐phase steam‐water flow,”
AIChE J., vol. 3, no. 1, pp. 136–142, 1957.
140
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2016LYSEI017/these.pdf
© [A. Neurouth], [2016], INSA Lyon, tous droits réservés
FOLIO ADMINISTRATIF
Prénoms : Adrien
Ecole doctorale :
MEGA (Mécanique, Energétique, Génie Civile, Acoustique)
Spécialité :
Mécanique
RESUME :
Une des voies d’amélioration des véhicules hybrides et électriques est l’utilisation de moteurs tournant
plus vite, jusqu’à plus de 42.000tr/min. Le but est d’augmenter la densité de puissance et le
rendement des groupes motopropulseurs. Pour utiliser ces moteurs de nouveaux réducteurs
mécaniques doivent être développés. Cela doit se faire sans générer de surcoût important face aux
solutions utilisées à basse vitesse et en assurant un niveau de performance énergétique élevé.
Cette thèse se situe en amont de la phase de conception d’un réducteur haute vitesse lubrifié par
barbotage. Elle a pour but d’identifier les problèmes d’échauffement et de pertes de puissance ainsi
que de proposer des pistes d’amélioration énergétique.
Ce travaille propose la modélisation thermomécanique de l’étage grande vitesse (GV) du réducteur,
réalisée à l’aide de la méthode des réseaux thermiques. Ce modèle couple les pertes de puissance
avec les températures. Une attention particulière est portée sur la modélisation des roulements de
l’arbre GV. Un nouveau modèle thermomécanique de roulement est développé. Les pertes par
barbotage deviennent importantes à grande vitesse. Une méthode permettant de fortement les réduire
est caractérisée.
Directeur de thèse :
CHANGENET Christophe, VILLE Fabrice
Président de jury :
FILLON Michel
Composition du jury :
ROBBE VALLOIRE François, CHANGENET Christophe, FILLON Michel, SEABRA Jorge, VILLE Fabrice,
OCTRUE Michel