Approche lagrangienne Approche eulerienne Conclusion - Questions

Mecanique des milieux continus solides et fluides

Seance 3
Emmanuel Plaut
Jean-Philippe Chateau, Denis Funfschilling, Yann Gunzburger,
Rainier Hreiz, Mathieu Jenny, Marc Poncot & Jean-Baptiste Vaney
Cinematique avancee : deformations dun milieu continu

1 Approche lagrangienne
2 Approche eulerienne
Approche lagrangienne Approche eulerienne Conclusion - Questions

Principe de base : on fait une analyse differentielle locale

Donc en lagrangien, comme le mouvement est decrit par
(.,t) : D0 Dt
domaine initial domaine actuel
X 7 x = (X,t)
position initiale position actuelle
la question est : si X varie de dX, comment varie x ? dx = ?

dx = FdX avec F = X le tenseur gradient de la transformation

Approche lagrangienne Approche eulerienne Conclusion - Questions

Alors : comment mesurer les deformations ?

em ent 1
mouv
00
11
x11
00
00dx
11

X
1
0
0
1
dX

mo
uve
me
nt x00
11
00
2 11
dx
Approche lagrangienne Approche eulerienne Conclusion - Questions

Idee : regarder le produit scalaire entre un d X et un d X

em ent 1
mouv dx
dX 00
11
x11
00
00dx
11

X1
0
0
1
dX

mo
uve dx
me
nt x00
11
00
2 11
dx
Approche lagrangienne Approche eulerienne Conclusion - Questions

Idee : regarder le produit scalaire entre un d X et un d X

Le produit scalaire transporte se calcule a laide du tenseur des
dilatations de Cauchy
C = FT F
suivant la formule
dx dx = C(dX,dX ) = dX C dX .
Le mouvement nest pas deformant autour de X C = 1
e = 0
1 
si e = C 1 tenseur des deformations de Green-Lagrange.
2
Approche lagrangienne Approche eulerienne Conclusion - Questions

Interpretation du
tenseur des deformations de Green-Lagrange
1 
e = C1 C = FT F = 1 + 2e
2
donc ce tenseur permet de calculer un produit scalaire transporte

dx dx = dX C dX = dX dX + 2dX e dX

C et e symetriques sont diagonalisables sur une base orthonormee ;

vecteurs et valeurs propres ont une signification ; les composantes des
matrices representatives ont une signification (exercice 2.1)...
Approche lagrangienne Approche eulerienne Conclusion - Questions

Reformulation en termes du champ de deplacement

u = xX = gradient de la transformation
F = X x = 1 + X u
= tenseur des dilatations de Cauchy
C = FT F = 1 + u + uT + uT u
= tenseur des deformations de Green-Lagrange
1  1 
e = C1 = u + uT + uT u
2 2
Approche lagrangienne Approche eulerienne Conclusion - Questions

Cas de petites transformations

Si    
F 1 1 u 1

alors
1 
e = u + uT + uT u
2
1 
e u + uT =: tenseur des deformations linearise
2
Approche lagrangienne Approche eulerienne Conclusion - Questions

Interpretation du tenseur des deformations linearise

1 
= u + uT
2
tel que
dx dx dX dX + 2dX dX

dx
dX = dX e1 = = 1 + 211 1 + 11
dX
dX = dX e1 et dX = dX e2
\
= = angle de glissement = (dx, dx ) 212
2
Approche lagrangienne Approche eulerienne Conclusion - Questions

Cas dun mouvement plan deformant sans rotation

Pour dX = ex dX et dX = dY = ey dX :

configuration de reference : actuelle :

12 dX
22 dX

dY dY

2 dy

dx
12 dX
11
00
00
11 11
00

2

dX dX 11 dX
X x
Approche lagrangienne Approche eulerienne Conclusion - Questions

Approximation locale dun champ de deplacement

dx = dX + du avec u = u(X+X) u(X) du = u dX

Approche lagrangienne Approche eulerienne Conclusion - Questions

Decomposition locale dun champ de deplacement

dx = dX + du avec du = u dX = dX}
| {z + {zdX}
|
deformation rotation

= +

  1 
= vd avec = u uT
2
1
= rot(u) signification de rotationnel de u
2
Approche lagrangienne Approche eulerienne Conclusion - Questions

Decomposition locale dun champ de deplacement

dx = dX + du avec du = u dX = dX}
| {z + {zdX}
|
deformation rotation

= +

En general la dilatation volumique

d 3x  
J = = det F = det 1 + u
d 3X
ce qui donne en petite transformation (cf. lexercice de CT 1.10)
J = 1 + divu signification de divu = tru = tr
Approche lagrangienne Approche eulerienne Conclusion - Questions

Approche lagrangienne des deformations : bilan (1)

Tenseur gradient de la transformation :

F = X : dX 7 F dX = dx .

Tenseur des dilatations de Cauchy :

C = FT F : (dX,dX ) 7 C(dX,dX ) = dx dx

Tenseur des deformations de Green-Lagrange :

1  1
e = C 1 : (dX,dX ) 7 e(dX,dX ) = (dx dx dX dX )
2 2
Petite transformation :

F 1 = u 1 avec u = x X le champ de deplacement

Tenseur des deformations linearise :

  1 
= e a lordre u = u + uT
2
Approche lagrangienne Approche eulerienne Conclusion - Questions

Approche lagrangienne des deformations : bilan (2)

Decomposition locale du champ de deplacement :

du = u dX =
| {zdX} + dX}
| {z
rotation deformation

Rotationnel de u :
1
= rot(u)
2
Divergence de u :
divu = tru = tr
Dilatation volumique :

d 3x
J = = det F
d 3X
ce qui donne en petite transformation

J = 1 + divu
Approche lagrangienne Approche eulerienne Conclusion - Questions

Approche eulerienne : bilan (cf. poly)

Tenseur gradient de vitesse :

d(dx)
K = x v : dX 7 K dx = .
dt
Divergence de la vitesse :

1 d(d 3 x)
divv = = taux de dilatation volumique
d 3 x dt
Tenseur des taux de deformation :
1  1 d(dx dx )
D = K + KT : (dx,dx ) 7 D(dx,dx ) =
2 2 dt
Decomposition locale du champ de vitesse :
1
dv = v dx = |
{zdx} + D dx} avec = rot(v)
| {z 2
rotation deformation
Approche lagrangienne Approche eulerienne Conclusion - Questions

En conclusion...
Il importe detre capable de faire le lien entre les deux approches
voir le document de cours
voir le probleme 2.1 que lon traitera en TD,
et qui consiste a etudier une transformation homogene
!
F homogene i.e. ne depend pas de X
x = x0 + F X avec F indep. de X

de cisaillement pur :

dX dx
00
11 0
1
1
0
X11
00
00
11 x 0
1
dX dx
Approche lagrangienne Approche eulerienne Conclusion - Questions

Attention :
toutes les transformations ne sont pas homogenes !

Cest le cas par exemple de lextrusion :

X x

[ Zienkiewicz, O. C. & Taylor, R. L. 2000 The finite elements method ]

Approche lagrangienne Approche eulerienne Conclusion - Questions

Rappel : equipe pedagogique pour le TD :

Charge de TD Laboratoire Specialite Groupe(s) Salle

Y. Gunzburger Laego Geomecanique X1 337
M. Jenny Lemta Meca. des fluides X2 338
R. Hreiz LRGP Meca. des fluides X3 & Y3 339
J.-B. Vaney IJL Meca. des materiaux X4 & Y4 341
J.-P. Chateau IJL Meca. des materiaux X5 342
D. Funfschilling LRGP Meca. des fluides Y1 337
M. Poncot IJL Meca. des materiaux Y2 338

Avez-vous des questions ?