Cours Lubrification VES4

Cours de Lubrification industriel
Professeur: R.BENMESSAOUD
Niveau :Deuxième année cycle ingénieur (CI2)
Année scolaire :2020/2021
Powerpoint Templates
Equation de Reynolds Positionnement du problème
Un film de fluide visqueux sépare les deux surfaces 𝑠1 et 𝑠2 d’un système
La charge W que peut supporter le contact,

La force F ou le couple de frottement C ,
Objectifs Calcul
Le débit Q du fluide dans le mécanisme,
La puissance P dissipée dans le contact.
• Déterminer le champ de pression et le champ de vitesse dans le fluide
• La force ou le couple de frottement peuvent être obtenu par l’intégration des contraintes de cisaillement aux
surfaces du contact
• La pression peut être créée par une pompe extérieure au contact
c’est l’hydrostatique :
Certains paliers et des glissières de certaines machines outils.

• La pression peut être créée par le déplacement relatif des
surfaces (hydrodynamique) :
Les butées, les paliers fluides, les garnitures mécaniques….

Pr.BENMESSAOUD,ENSAM-CASA , 2020/2021
Equation de Reynolds Positionnement du problème
• On considère le parallélépipède de volume dxdydz soumis à des contraintes de pression
• Au cas où la distance h entre les 2 plans est suffisamment petite et

l’épaisseur du film ne varie pas significativement suivant (Ox)
La pression peut être supposé constant suivant l’épaisseur
Détermination de la contrainte de cisaillement
Condition d’équilibre
u représente la vitesse du fluide selon 0x

Avec
μ sa viscosité supposée constante
Equation de Reynolds Forme simplifiée
Un développement limité au premier ordre permet d’écrire :
Equation de Navier Stokes simplifiée (écoulement plan)
Ainsi
Conditions au limites sur la vitesse ( pas de glissement des « pillicules » au niveau du paroi)
u=U si y=0
u=0 si y=h
Equation de Reynolds Forme simplifiée
Le débit q(x) par unité de largeur à travers une section droite :
Le débit q (x) est constant si le fluide incompressible et l’ écoulement est unidirectionnel
On obtient une équation très simplifiée de l’équation de Reynolds:
Avec h∗ , l’épaisseur du film pour un gradient de pression nul

Equation de Reynolds Forme générale
Pour un fluide newtonien, les équations de la mécanique des milieux continus conduisent à l’équation de Reynolds :
𝑥𝑖 et t: variables d’espace et de temps
Loi fondamentale de la dynamique 𝑢𝑖 : composantes de la vitesse

Ρ: masse volumique du fluide
𝑓𝑖 :forces massiques extérieures
𝜎𝑖 ∶tenseur des contraintes
Loi de comportement rhéologique(fluide newtonien) 𝛿𝑖𝑗 ∶ symbole de Kronecker
(𝛿𝑖𝑗 = 1 si i = j ; 𝛿𝑖𝑗 = 0 si i ≠ j)
λ et μ : coefficients de viscosité de
Navier (fonctions de l’espace et du temps)
θ est le taux de dilatation cubique obtenu par la relation suivante:
𝜀𝑖𝑗 : tenseur des taux de déformations obtenu par la relation suivante:
Remplaçant 𝜎𝑖𝑗 par son expression (comportement rhéologique) dans la loi fondamental de la dynamique (forces
massiques extérieures négligées) ,on obtient :
Equations de Navier
i correspond à la direction considérée et j joue le rôle d'indice muet
En lubrification, la dimension suivant l’épaisseur du film (0, x2 ) est très faible par rapport aux deux autres.
Changement de variable :
x
x1 = L1
u
u1 = V1
x μ
x2 = H2 μത =μ
u2 L 0
x u2 = VH
x3 = L3 ത λ
λ=
u λ 0
u3 = V3
tV
t=ҧ L
Dans ces expressions L et V sont respectivement la dimension et la vitesse caractéristiques selon les axes celles suivant (0,
x2 )
L
, μ0 , λ0 définissent l’ordre de grandeur du temps et des coefficients de viscosité
V
Nous supposerons que ces deux coefficients sont du même ordre de grandeur
H2
On pose donc : pത =p μ0 VL
Ce changement de variable met en évidence les 3 nombres sans dimension suivants :
VHρ0
ℜ= Nombre de Reynolds relatif à l’épaisseur de l’écoulement
μ0
𝐻
ε= Paramètre d’échelle, qui en lubrification est de l’ordre de 10−3
𝐿
𝜆0
Rapport des coefficients de viscosité, généralement de l'ordre de 1.
𝜇0
Compte tenu de ces changements de variable l’équation de Navier s’écrit :
En négligeant les termes multipliés soit par 𝜀 2 , soit par 𝜀 3 𝑅 il vient :
Les termes en ε ℜ seront négligeables devant les termes de viscosité dans le cas ε ℜ << 1 ; ce qui est souvent
vérifié dans les écoulements rencontrés en lubrification:
Dans le cas du contact d’une bille sur un plan (roulement à billes ou engrenage):
L=10−3 𝑚, 𝐻 = 10−6 𝑚, 𝑉 = 10 𝑚𝑠 −1
𝜇0 = 0.05 𝑃𝑎. 𝑠, 𝜌0 = 900 𝑘𝑔/𝑚3 , ε ℜ=1.8 10−4
Equation de Reynolds Ordre de grandeur des paramètres
Dans le cas d’un palier à huile de dimensions moyennes :
L=0.05 m, H = 10−4 m, V = 30 ms −1
μ0 = 0.05 Pa. s, ρ0 = 900 kg/m3 , ε ℜ=1.8 10−2
Quand la viscosité du fluide est faible ou l’épaisseur du film est importante
μ0 = 1 mPa. s, ρ0 = 1000 kg/m3 , ε ℜ=0.6
Dans le cas des paliers de grandes dimensions:
L=0.05 m, H = 10−3 m, V = 30 ms−1
μ0 = 0.05 Pa. s, ρ0 = 900 kg/m3 , ε ℜ=3.6
Equation de Reynolds
Si on néglige les forces d’inertie on aura:
Et, en revenant aux variables dimensionnées en coordonnées cartésiennes (0, x, y, z) :
Equation de Reynolds Résolution
Si μ est indépendant de y, on intègre les équations précédentes
Pour déterminer les constantes A, B, C et D , on utilise les conditions aux limites sur la vitesse avec l’hypothèse
de non glissement du fluide sur les parois
y=0 pour u=𝑈1 ,v=0,w=𝑊1

y= h pour u=𝑈2 ,v= 𝑉2 ,w=𝑊2
Système d'axes en coordonnées cartésiennes
𝑈1 = 𝐵
1 𝜕𝑝 2 𝐴
𝑈2 = ℎ + ℎ+𝐵
2𝜇 𝜕𝑥 𝜇
𝐴 𝑈2 − 𝑈1 1 𝜕𝑝
= − ℎ
𝜇 ℎ 2𝜇 𝜕𝑥
𝑈2 − 𝑈1 1 𝜕𝑝
𝑢= 𝑦+ 𝑦 − ℎ + 𝑈1
ℎ 2𝜇 𝜕𝑥
𝑤2 − 𝑤1 1 𝜕𝑝
𝑤= 𝑦+ 𝑦 𝑦 − ℎ + 𝑤1
ℎ 2𝜇 𝜕𝑧
Ces vitesses sont remplacé par leurs expression dans l’équation de continuité.
En intégrant, à travers l’épaisseur, on obtient :
ℎ
𝜕(𝜌𝑣)
න 𝑑𝑦 = 𝜌𝑉2 𝐶𝑎𝑟 𝑉1 = 0
0 𝜕𝑦
Si ρ est indépendant de y:
Pour intégrer les deux autres termes, il faut intervertir les signes dérivés et sommes en remarquant que :
En regroupant tous les termes on obtient l’équation de Reynolds:
Les équations donnant les vitesses dans le fluide
ℎ−𝑦 1 𝜕𝑝 𝑦
𝑢= 𝑈1 + 𝑦 𝑦 − ℎ + 𝑈2
ℎ 2𝜇 𝜕𝑥 ℎ
ℎ−𝑦 1 𝜕𝑝 𝑧
𝑤= 𝑤1 + 𝑦 𝑦 − ℎ + 𝑤2
ℎ 2𝜇 𝜕𝑧 ℎ
Les contraintes de cisaillement dans le fluide :
𝜕𝑢 1 𝜕𝑝 𝜇
𝜏𝑥𝑦 = 𝜇 𝜕𝑦 = 2 𝜕𝑥 2𝑦 − ℎ + ℎ (𝑈2 - 𝑈1 )
𝜕𝑤 1 𝜕𝑝 𝜇
𝜏𝑦𝑧 = 𝜇 = 2𝑦 − ℎ + (𝑊2 − 𝑊1 )
𝜕𝑦 2 𝜕𝑧 ℎ
Dans le système d'axes (O, r, θ, z) (Cylindrique), les vitesses du fluide sur les surfaces du contact s’écrivent :
z=0 pour u=𝑈1 ,v= 𝑉1 ,w=0

z= h pour u=𝑈2 ,v= 𝑉2 ,w=𝑊2
L'équation de Reynolds en coordonnées cylindriques est donnée par :
Les composantes de la vitesse dans le film s'écrivent :
ℎ−𝑧 1 𝜕𝑝 𝑧
𝑢= 𝑈1 + 𝑧 𝑧 − ℎ + 𝑈2
ℎ 2𝜇 𝜕𝑟 ℎ
ℎ−𝑧 1 𝜕𝑝 𝑧
v= ℎ
𝑉1 + 2𝜇𝑟 𝜕𝜃 𝑧 𝑧 − ℎ + ℎ 𝑉2
Les contraintes de cisaillement dans le film sont :
𝜕𝑢 1 𝜕𝑝 𝜇
𝜏𝑟𝑧 = 𝜇 𝜕𝑧 = 2 𝜕𝑟 2𝑧 − ℎ + ℎ (𝑈2 - 𝑈1 )
𝜕𝑣 1 𝜕𝑝 𝜇
𝜏𝜃𝑧 =𝜇 = 2𝑧 − ℎ + (𝑉2 − 𝑉1 )
𝜕𝑧 2𝑟 𝜕𝜃 ℎ
Hypotheses (Equation de reynolds)
• L’épaisseur du film est très faible devant les autres dimensions du contact
• Le milieu est continu (Pas de rupture du film)
• Le fluide est newtonien (proportionnalité entre le taux de cisaillement et les contraintes de cisaillement)
• Le régime est laminaire
• Une des surfaces du contact est parfaitement plane et ne présente pas de rugosité
• Il n’y a pas de glissement entre le fluide et les parois du contact (Vitesse du fluide au paroi est égale à celle de la
paroi)
• Les forces d’inertie sont négligeables.
• Les forces massiques extérieures sont négligées
• La viscosité et la masse volumique du fluide ne varient pas à travers l'épaisseur du film
• La température ne varie pas à travers l’épaisseur du film :
T=T(𝑥1 , 𝑥3 , 𝑡) Pr.BENMESSAOUD,ENSAM-CASA , 2020/2021
μ (T,ρ) et ρ (T, p) ne sont fonction que de la température et de la pression
Equation de Reynolds Contacts hydrostatiques
Il existe des formes plus générales de l'équation de Reynolds qui s'affranchissent de plusieurs hypothèses
Elle s ’appliquent à différents lubrifiants et à différents types de contact.
Pour les contacts hydrostatiques les seuls contacts pour lesquels , l’équation de Reynolds devient:
Coordonnée cartésiennes
Coordonnées cylindriques
Equation de Reynolds Exemple de problèmes
L’exemple du contact formé d’un cylindre tournant autour de son axe , à la vitesse
angulaire ω, en regard d’un plan se déplaçant parallèlement à lui même, à la vitesse
U, illustre ce problème (figure ci-contre).
Dans le cas d’un écoulement incompressible, bidimensionnel, l’équation de Reynolds se
réduit à :
𝜕 ℎ3 𝜕𝑝 𝜕ℎ
= 6(𝑈1 −𝑈2 ) 𝜕𝑥 + 12𝑉2
𝜕𝑥 𝜇 𝜕𝑥
Dans le système d’axes choisi, la vitesse des surfaces s’écrit :
𝑈1 = 𝑈
𝑈2 = 𝜔𝑅 cos 𝛼
𝑉2 = 𝜔𝑅 sin 𝛼
Equation de Reynolds Exemple de problèmes
et en regroupant tous les termes il vient :
𝜕ℎ
tan 𝛼 ≈ sin 𝛼 ≈ 𝛼 ≈
𝜕𝑥
cos 𝛼 ≈1
𝜕ℎ
𝑈2 = 𝜔𝑅 𝑒𝑡 𝑉2 = 𝜔𝑅
𝜕𝑥
Ainsi l’équation de Reynolds est donnée par :
𝜕 ℎ3 𝜕𝑝 𝜕ℎ
= 6(𝑈 +𝜔𝑅)
𝜕𝑥 𝜇 𝜕𝑥 𝜕𝑥
Le second membre de l’équation de Reynolds est proportionnel à U + ωR (vitesse de roulement) et c’est la somme des
vitesses tangentielles du contact qui génère la portance dans le film.
Equation de Reynolds Problème 1 (Patin de Rayleigh )
Soit le patin de Rayleigh (Figure ci-dessous), de largeur infinie selon Oz.
Le patin est constitué de deux régions : une région d’entrée où l’épaisseur du film est constante et égale à ℎ1 , une
région de sortie où l’épaisseur est également constante et égale à ℎ2 (ℎ1 < ℎ2 ).
Les régions sont séparées par une discontinuité dans l’épaisseur du film, située au point x = 𝐵0
Pour étudier le problème, on isole la discontinuité et on résout l’équation de Reynolds séparément dans les régions
d’entrée et de sortie
Dans la région d’entrée définie par O < x < B0 ,’équation de Reynolds s’écrit :
𝑑2𝑝
=0
𝑑𝑥 2
D’où
𝑝 = 𝐶1 𝑥 + 𝐶2
𝐶1 et 𝐶2 sont calculées à partir des conditions aux limites qui s’écrivent :
P=0 pour x=0

P= 𝑝𝑚 pour x=B
où 𝑝𝑚 est la pression inconnue qui existe au niveau de la discontinuité.
Et en regroupant tous les termes il vient :
𝑝𝑚
Donc : 𝑝= 𝑥
𝐵0
De même dans la zone de sortie définie par B0 < x < B , l’équation de Reynolds s’écrit :
𝑑2𝑝
=0
𝑑𝑥 2
𝑝 = 𝐶3 𝑥 + 𝐶4
P=0 pour x=B

Les conditions aux limites
P= 𝑝𝑚 pour x=B
𝑝𝑚
𝑝= 𝑥
𝐵0
𝑝𝑚
𝑝= (B-x)
𝐵−𝐵0
Ainsi, la variation de pression est linéaire dans le contact (figure). La pression inconnue Pm est calculée en écrivant
l’égalité des débits à l’entrée et à la sortie.
Le débit pour une largeur L s’écrit à l’entrée ou à la sortie :
ℎ
𝑄𝑘 = 𝐿 න 𝑢𝑘 𝑑𝑦
0
1 𝑑𝑝𝑘 ℎ𝑘 −𝑦
𝑢𝑘 = 2𝜇 y(y-ℎ𝑘 )+ 𝑢
𝑑𝑥 ℎ𝑘
Avec
𝐿ℎ𝑘3 𝑑𝑝𝑘 𝐿𝑈ℎ𝑘
𝑄𝑘 = − 12𝜇 𝑑𝑥 + 2
𝑑𝑝1 𝑃𝑚
=
𝑑𝑥 𝐵0
Et
𝑑𝑝2 −𝑃𝑚
=
𝑑𝑥 𝐵−𝐵0
L’égalité des débits Q1 et Q2 permet d’écrire

ℎ −ℎ2
𝑃𝑚 =6μUℎ3 1 3
1 + ℎ2
𝐵0 𝐵−𝐵0
Pour une largeur L, la charge est:
𝐵𝑃𝑚
𝑊=𝐿
2
Avec
𝐵0 ℎ1
𝑠= 𝑎=
𝐵 ℎ2
6μU 𝑠 1−𝑠 𝑎−1

𝑃𝑚 = ℎ2 𝑎3 1−𝑠 +𝑠
2
Les valeurs de a et s qui rendent la charge

maximale sont :
a = 1,866 ; s = 0,718
On considère un patin à échelon, ou patin de Rayleigh, semblable à celui schématisé sur la figure ci-dessous. La forme
du contact et les conditions de fonctionnement sont telles que :
Longueur totale : B = 0,5 m

Longueur de la zone d’entrée : 𝐵0 = 0,36 m
Largeur: L = 1 m
Epaisseurs du film : ℎ1 =0,2 mm ; ℎ2 =0,1 mm
Vitesse de translation : U = 6 m/ s
Viscosité dynamique du lubrifiant : μ = 0,01Pa.s
Si on néglige l’écoulement axial, les relations présentées précédemment permettent de calculer les caractéristiques de
fonctionnement :
La pression maximale : P= 1,226 106 Pa

La portance : W = 306 103 N
Le débit : Q = 3,7 10−4 𝑚3 / s .
La force de frottement est obtenue par intégration du taux de cisaillement :
𝑑𝑢 1 𝑑𝑝 𝑈
𝜏𝑥𝑦 = μ = (2y-h) −μ
𝑑𝑦 2 𝑑𝑥 ℎ
sur la surface inférieure : y = 0. On obtient :
𝐿 𝐵 𝑃𝑚 𝐵0 𝐵−𝐵0
F=-‫׬‬0 ‫׬‬0 (𝜏𝑥𝑦 )𝑦=0 𝑑𝑥𝑑𝑧 =𝐿 ℎ1 − ℎ2 + μ𝑈 ( + )
2 ℎ1 ℎ2
La puissance dissipée par cisaillement est : P = F.U = 1520 W
et le coefficient de frottement : 𝐹
f= 𝑊 ≈ 8.3 10−4

Cours Lubrification VES4

Transféré par

Droits d'auteur :

Formats disponibles

Cours Lubrification VES4

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours Lubrification VES4

Transféré par

Droits d'auteur :

Formats disponibles

Cours de Lubrification industriel

Un film de fluide visqueux sépare les deux surfaces 𝑠1 et 𝑠2 d’un système

La charge W que peut supporter le contact,

Certains paliers et des glissières de certaines machines outils.

Les butées, les paliers fluides, les garnitures mécaniques….

• On considère le parallélépipède de volume dxdydz soumis à des contraintes de pression

• Au cas où la distance h entre les 2 plans est suffisamment petite et

La pression peut être supposé constant suivant l’épaisseur

Détermination de la contrainte de cisaillement

u représente la vitesse du fluide selon 0x

Un développement limité au premier ordre permet d’écrire :

Equation de Navier Stokes simplifiée (écoulement plan)

Le débit q(x) par unité de largeur à travers une section droite :

Le débit q (x) est constant si le fluide incompressible et l’ écoulement est unidirectionnel

On obtient une équation très simplifiée de l’équation de Reynolds:

Avec h∗ , l’épaisseur du film pour un gradient de pression nul

𝑥𝑖 et t: variables d’espace et de temps

Loi fondamentale de la dynamique 𝑢𝑖 : composantes de la vitesse

θ est le taux de dilatation cubique obtenu par la relation suivante:

𝜀𝑖𝑗 : tenseur des taux de déformations obtenu par la relation suivante:

i correspond à la direction considérée et j joue le rôle d'indice muet

Ce changement de variable met en évidence les 3 nombres sans dimension suivants :

𝜇0 = 0.05 𝑃𝑎. 𝑠, 𝜌0 = 900 𝑘𝑔/𝑚3 , ε ℜ=1.8 10−4

μ0 = 0.05 Pa. s, ρ0 = 900 kg/m3 , ε ℜ=1.8 10−2

Quand la viscosité du fluide est faible ou l’épaisseur du film est importante

μ0 = 1 mPa. s, ρ0 = 1000 kg/m3 , ε ℜ=0.6

Dans le cas des paliers de grandes dimensions:

L=0.05 m, H = 10−3 m, V = 30 ms−1

μ0 = 0.05 Pa. s, ρ0 = 900 kg/m3 , ε ℜ=3.6

Et, en revenant aux variables dimensionnées en coordonnées cartésiennes (0, x, y, z) :

y=0 pour u=𝑈1 ,v=0,w=𝑊1

Système d'axes en coordonnées cartésiennes

En intégrant, à travers l’épaisseur, on obtient :

Les équations donnant les vitesses dans le fluide

z=0 pour u=𝑈1 ,v= 𝑉1 ,w=0

Les composantes de la vitesse dans le film s'écrivent :

Les contraintes de cisaillement dans le film sont :

Dans le système d’axes choisi, la vitesse des surfaces s’écrit :

𝐶1 et 𝐶2 sont calculées à partir des conditions aux limites qui s’écrivent :

P=0 pour x=0

où 𝑝𝑚 est la pression inconnue qui existe au niveau de la discontinuité.

P=0 pour x=B

L’égalité des débits Q1 et Q2 permet d’écrire

Pour une largeur L, la charge est:

6μU 𝑠 1−𝑠 𝑎−1

Les valeurs de a et s qui rendent la charge

Longueur totale : B = 0,5 m

La pression maximale : P= 1,226 106 Pa

La force de frottement est obtenue par intégration du taux de cisaillement :

sur la surface inférieure : y = 0. On obtient :

La puissance dissipée par cisaillement est : P = F.U = 1520 W

Vous aimerez peut-être aussi