Mef NL3 2018

Génie Mécanique
Chapitre V
Contact & MEF
M. LAROUSSI HELLARA
Génie Mécanique
Introduction
A m km
Atomes Procédés de fabrication Couches tectoniques

M. LAROUSSI HELLARA Copyright © 2016
Crédit photos : CS – NECS – Altair Engineering – Sirehna - EDF
Contact rails/roues
Mécanique du contact :  Non Linéaire

 Non régulière
Exemple : Loi de coulomb

F
Force de frottement
V
81 Vitesse de glissement
Génie Mécanique
Introduction (suite)
Hypothèses :  Statique
 HPP
 Elasticité linéaire
 div  1  f 1  0
 1
   R
1
:  1
  1
  ( grad u )
1 s

Surfaces potentielles c  CL1 + Contact
de contact
c  div  2  f 2  0
 2
   R
2
:  2
 2

  ( grad u )
2 s
82 
 CL2
Génie Mécanique
Introduction (suite)
But : Résoudre numériquement le problème
 L’écrire sous la forme K .U   F 

Plan :
I- Pb simple : membrane tendue
II- Pb de l’obstacle rigide pour une

membrane tendue
 Modèle de contact
 Algorithme de résolution
83
Génie Mécanique
Partie 1 : Pb de la membrane tendue
e p

e
   R2
e
Modélisation :
Modèle 2D avec : u3  u (u1  u2  0)
1
 13  u3,1   1
2
1
 23  u3, 2   2
2
  2T0  soit  13  T0 u,1   1
84  23  T0 u, 2   2
Génie Mécanique
Partie 1 (suite)
Formulation Locale du pb : F. Forte

Trouver u,  et  /
div   p  0 dans   R 2 (borné et régulier)

  2T0 
1
  grad u
2
u  0 sur 
Formulation Globale du pb :  div   p  w dS  0 w
IPP : On affaiblit la régularité de 
85
 div  w dS    .n w ds    . grad w dS
 
Génie Mécanique
Partie 1 (suite)
Trouver u C.A. /
Pm
  (u ) . grad w dS   p w dS   .n w ds

w CA 0
De la forme : Trouver u C.A. / w CA 0 a(u , w)  l ( w)

Avec ici :
a (u , w)   T0 grad u . grad w dS énergie interne virtuelle


l ( w)   p w dS Travail virtuel des efforts extérieurs


86
Génie Mécanique
Partie 1 (suite)
Equivalence énergétique : Th de stationarité de l’Ep
L’énergie potentielle totale : E p ( w)  W ( w)   ( w)

1
W ( w)   T0 grad w . grad w dS
2
87
 ( w)   p w dS  l ( w)

Génie Mécanique
Partie 1 (suite)
L’énergie potentielle totale : E p ( w)  W ( w)   ( w)

1
W ( w)   T0 grad w . grad w dS
2
 ( w)   p w dS  l ( w)

Trouver u C.A. / w CA 0
Qm
E p (u )  E p ( w)
Et on mq : Pm  Qm
1
Rq : E p ( w)  a( w, w)  l ( w)
2
88
Génie Mécanique
Partie 1 (suite)
Approximation numérique :
Trouver uh CA h /wh CA 0 h
Pmh
a (uh , wh )  l ( wh )
avec CA h  CA / CA h de dimension finie N h
Soit Bh =(N1, N2, …, NNh) une base de CA h

Nh
On écrit alors : uh   N i uhi

i 1
Et on prend : wh  N i : Fonction de forme

Nh
(ou encore wh   N i whi

89 i 1 factorisé puis éliminé)
Génie Mécanique
Partie 1 (suite)
Nh
D’où  a( N
j 1
j , Ni ) u j  l ( Ni )
K .U   F 
Soit
K ij  a ( N i , N j ) Fi  l ( N i )
Ici : K ij   T0 grad N i . grad N j dS


Fi   p N i dS

90
91
Génie Mécanique
Partie 1 (suite)
Traitement des CL cinématiques :
Génie Mécanique
Partie 1 (suite)
 Méthode de relèvement :
u
(valable si contact bilatéral mais pas unilatéral)
u0
u  u0  u *
avec u 0 CA et
*
u CA0 
soit a (u, w)  a(u 0  u * , w)  l ( w)

 a (u * , w)  l ( w)  a(u 0 , w)
Rq :
Cette méthode marche en 1D
en 2D quadrangulaire
mais pas en 2D courbe
en raison de la variation de la normale
92
Génie Mécanique
Partie 1 (suite)
 Méthode de pénalisation (ou de pénalité) :
Trouver uK C.A. réalisant
QK
Inf E p ( w)  EK ( w, u )
w  CA0
1
avec E p ( w)  a( w, w)  l ( w)
2
K
EK ( w, u )    2
et ( w u ) ds
2 
 Méthode approchée dont on contrôle la convergence vers

la solution exacte
93
Génie Mécanique
Partie 1 (suite)
Ou encore :
Trouver uK C.A. / w CA 0
PK
a (u , w)  l ( w)  K  (u  u ) w ds  0

La matrice de rigidité devient alors :
K   CK 
avec CK ij  K  N i N j ds

 Choix de K : suffisamment grand % aux données homogènes

94
Génie Mécanique
Partie 1 (suite)
 Méthode de pénalisation (ou de pénalité) :
Trouver u C.A. réalisant
QK
Inf E p ( w)  EK ( w, u )
w  CA0
1
avec E p ( w)  a( w, w)  l ( w)
2
K
EK ( w, u )    2
et ( w u ) ds
2 
 Méthode approchée dont on contrôle la convergence vers

la solution exacte
95
Génie Mécanique
Partie 1 (suite)
 Méthode de Lagrange :
 Méthode alternative exacte MAIS introduction des réactions
96
Génie Mécanique
Partie 2 : Pb. de l’obstacle rigide
e p

e
e
ddl : flêche u ( x1 , x2 )
Hypothèse : HPP f c  f c e 3 action extérieure obstacle/membrane
Inconnues : - efforts de contact

- étendue de la zone de contact effective
 CL non standards!
97
Génie Mécanique
Partie 2 (suite)
Formulation Locale du pb :
Trouver u et  /
div   p  f c  0
  T0 grad u
+ Loi de contact : Modèle de Signorini (1933)
u0 sur 
Formulation Globale :
FV+IPP :  T
0 grad u . grad w dS   p w dS   f c w dS
 
98
Génie Mécanique
Partie 2 (suite)
Equation de surface de l’obstacle :
O  ( x1 , x2 , x3 )  R 3 ; x3   ( x1 , x2 ) ( x1 , x2 )  c 
Où c est la zone potentielle de contact (majorant la zone de contact)
Rq :  ( )  0 Condition de cohérence x3  
Modèlede Signorini : Loi de contact unilatéral sans frottement

- Condition cinématique de non inter-pénétration
(principe de Fermy sur les atomes) portant sur u % 
d  u   0 sur 
- Condition statique portant sur f c : compressive et non attractive
f c  0 sur 
+ Condition d’exclusion : Condition alternative d f c  0 sur 
99 si pas de contact f c  0 si contact f c  f c (u )
Génie Mécanique
Partie 2 (suite)
Contacts admissibles :
C   w; d ( w)  w   0 ( x1 , x2 )  
w  w'u / w' CA  C u  CA  C
Soit
Alors  T
0 grad u . grad ( w'u ) dS   p ( w'u ) dS   f c ( w'u ) dS

0 
( w' )  (  u )
0
0 si f c  0
Soit  T
0 grad u . grad ( w'u ) dS   p ( w'u ) dS

100
Génie Mécanique
Partie 2 (suite)
Le pb devient :
Trouver u  CA  C / w  CA  C
PIV
 T
0 grad u . grad ( w  u ) dS   p ( w  u ) dS
Equivalence énergétique :
Trouver u  CA  C / w  CA  C
PIV QIV
E p (u )  E p ( w)
1
avec E p ( w)  a( w, w)  l ( w)
2
Rq : C n’est pas un espace vectoriel mais c’est un ensemble convexe
101
Génie Mécanique
Partie 2 (suite)
 
Soit :
K

Inf  E p ( w)   ( w  ) dS 
w CAC
 2

 2 
D’où :
Trouver u  CA/ w  CA
QIVK
Inf  E p ( w)  EK ( w)
w  CA
Et on mq :
Trouver u  CA/ w  CA NL
QIV
K
PIVK
a (u , w)  l ( w)  K  (u  )  w dS  0

102
Génie Mécanique
Partie 2 (suite)
Algorithme de linéarisation numérique :
Newton-Raphson? Ici, la partie + est une fonction non dérivable!
  0 : Je colle la S à  i
On a : K   CK  avec ui  i 
 0 : J' accepte la S
 Méthode des statuts : Pour chaque point, on définit un statut :

1 si contactant 0 si pas de contact
Initialisation : Non contactant

puis on applique la pénalisation et on étudie la convergence
 0  0
  1   
    2 
( St )  K   K  1  .U   F  K  3 
    
  1    4
103   0   0 


Mef NL3 2018

Transféré par

Droits d'auteur :

Formats disponibles

Mef NL3 2018

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

Mef NL3 2018

Transféré par

Droits d'auteur :

Formats disponibles

Génie Mécanique

Atomes Procédés de fabrication Couches tectoniques

Mécanique du contact :  Non Linéaire

Exemple : Loi de coulomb

 L’écrire sous la forme K .U   F 

II- Pb de l’obstacle rigide pour une

Formulation Locale du pb : F. Forte

div   p  0 dans   R 2 (borné et régulier)

Formulation Globale du pb :  div   p  w dS  0 w

IPP : On affaiblit la régularité de 

De la forme : Trouver u C.A. / w CA 0 a(u , w)  l ( w)

a (u , w)   T0 grad u . grad w dS énergie interne virtuelle

l ( w)   p w dS Travail virtuel des efforts extérieurs

L’énergie potentielle totale : E p ( w)  W ( w)   ( w)

L’énergie potentielle totale : E p ( w)  W ( w)   ( w)

avec CA h  CA / CA h de dimension finie N h

Soit Bh =(N1, N2, …, NNh) une base de CA h

On écrit alors : uh   N i uhi

Et on prend : wh  N i : Fonction de forme

(ou encore wh   N i whi

Ici : K ij   T0 grad N i . grad N j dS

soit a (u, w)  a(u 0  u * , w)  l ( w)

 Méthode approchée dont on contrôle la convergence vers

La matrice de rigidité devient alors :

 Choix de K : suffisamment grand % aux données homogènes

 Méthode approchée dont on contrôle la convergence vers

Hypothèse : HPP f c  f c e 3 action extérieure obstacle/membrane

Inconnues : - efforts de contact

Rq :  ( )  0 Condition de cohérence x3  

Modèlede Signorini : Loi de contact unilatéral sans frottement

 Méthode des statuts : Pour chaque point, on définit un statut :

Initialisation : Non contactant

Vous aimerez peut-être aussi