Chapitre V:

METHODES DE RESOLUTION
DES PROBLEMES DE MMC

I- Introduction :

Soit un corps élastique soumis aux conditions extérieures suivantes :

•Forces de volume fV (N/m3) sur le corps entier.
•Forces de surfaces F (N/m²) sur certaines parties de sa frontière.
•Champs de variation de la température (T-T0) sur tout le corps.
•Déplacements imposés sur certaines parties de sa frontière.

Bilan des équations et des inconnues :

 •15 Equations (respectant les conditions aux limies) :

•(03) Equations d’équilibre (03) : divσ + ρfv = 0

•(06) Relations déformations-déplacements: (
ε = 1 . gradU + gradTU
2
)
•(06) Lois de comportement: 1 +ν
(
1 − 2ν 1 − 2ν
())
σ = E . ε + ν .Tr ε .I − 2E (T − T0 ).I

•15 Inconnues :

•(03) Déplacements Ui .

•(06) Contraintes σij .

•(06) Déformations εij .

On a un système d’équations différentielles de 15 équations à 15

inconnues relativement délicat à résoudre.

Pour la résolution des problèmes d’élasticité, on utilise deux types de formulations :

•Formulation en déplacement : les équations de NAVIER.
•Formulations en contraintes : les équations de BELTRAMI.
 II- Formulation en déplacement : Equation de
NAVIER :

équations d’équilibre nous donnent : () r

div σ + ρ f = 0
Avec la loi de comportement d’élasticité linéaire, nous obtenons :

() ( ) r
2 µ div ε + λ grad trε + ρ f = 0
On obtient alors les formes vectorielles suivantes des équations de
NAVIER :
r
µ ∆ ( u ) + ( µ + λ ) grad ( div( u ) ) + ρ f = 0
r r

En utilisant l’opérateur divergence, on obtient :

r
[ ( )]
(2 µ + λ ) ∆ [div(u )]+ div ρ f = 0
r
Cas particulier : Dans le cas d’un domaine à masse volumique
constante, en équilibre etrplacé dans un champ de force volumique
( r
dérivant d’un potentiel f = − grad p ⇒ div f = 0 )
Les équations de Navier deviennent :

∆(divU ) = 0
On peut écrire les équations de Navier en fonction de E et ν :

E
2(1 + ν )
()
.∆ u + E
2(1 + ν )(1 − 2ν )
[ ( )] + ρ f
.grad div u = 0

En thermo élasticité, les équations de Navier devient :

E
2(1 + ν )
()
.∆ u + E
2(1 + ν )(1 − 2 ν )
.grad div u [ ( )]
+ ρf − E α .grad ( T − T0 ) = 0
1 − 2ν
 Ces équation de Navier doivent vérifier les conditions aux limites en
fonctions des déplacements suivantes :
•Déplacement imposés :
u( N ) = d ∀N ∈ ∂Dd
•Efforts imposés :
( )
T N, n = F ∀N ∈ ∂DF
 Lois de l’élasticité linéaire :
Equations d’équilibre (
ε = 1 . gradU + gradTU
2
)
divσ + ρfv = 0
σ =
1+ν
( 1 − 2ν
())
E . ε + ν .Tr ε .I − 2E (T − T0 ).I
1 − 2ν

Equations de NAVIER :
E
2(1 + ν )
.∆ u + () E
2(1 + ν )(1 − 2ν )
.grad div u [ ( )] + ρ f = 0

Conditions de
symetrie et les Proposi-
conditions aux tion de
limites en U(x).
déplacements.

U(x) est élastostatiquement accéptable

Calcul de :
(
ε = 1 . gradU + gradTU
2
)
σ = E . ε + ν .Tr ε .I  − 2E (T − T0 ).I
1+ν  1 − 2 ν    1 − 2ν

Conditions aux
limites en T = σ .n
contraintes

Toutes les conditions aux limites sont vérifiés :

C’est la solution
 Processus de résolution par
l’approche en déplacement

• Postuler la forme du champs de

déplacement.
• Vérifier les conditions aux limites en
déplacements.
• Vérifier les équations de Lame-Navier.
• Calculer les déformations puis les
contraintes.
• Vérifier les conditions aux limites en
contraintes.
 TD1 :EQUATION DE NAVIER
.
Soit un corps cylindrique repéré dans un trièdre orthonormé direct. Pour
un point M de coordonnées (x1 ; x2 ; x3 ). Les forces de volumes sont
supposées négligeables. On propose le vecteur déplacement suivant :
r r
( 2 2
)
r
U ( M )=(2a x1 x2 )E1 + 7a x1 −3a x2 E2
(a est un nombre infiniment petit)

1-Ce champ de déplacement est-il élastostatiquement acceptable.

2-Déterminer le tenseur symétrique ε et le tenseur anti-symétrique ω
de la transformation.
3-Ecrire les équations d’équilibre. Quelle est la valeur du coefficient de
Poisson qui permet de les satisfaire en l’absence de forces de volume?
4-Avec les conditions précédentes, déterminer le tenseur des contraintes en
M dans la base (E1 , E2 , E3 )
r r r
 III- Formulation en contraintes : Equation de
BELTRAMI :
Les équations de BELTRAMI proviennent des équations de
compatibilité :

() () [ ( )]
T
   
grad div ε +  grad div ε  − grad  grad tr ε  − ∆ ε = 0
   
La loi de comportement nous permet d’écrire :

ε=
1 +ν
E
ν
σ −  tr σ
E
( ) I ⇒ ∆ε =
1 +ν
E
ν
( ( ))

∆σ −  ∆ tr σ  I
E 
D’autre part nous avons :

tr σ  = I I = E tr ε  = E θ
  1 − 2ν   1 − 2ν

∆ II = E ∆θ
1 − 2ν

Avec la forme scalaire des équations de NAVIER :

1+ν
( 2 µ + λ ) ∆θ = div[ ρ ( γ ) ] div[ ρ ( γ ) ]
r r
⇒ ∆ II =
1 −ν
De plus nous pouvons écrire :
1 +ν ν
div ε = div σ − grad (I1 )
E E
  1 +ν
grad  div ε  =
  E
  ν
[
grad  div σ  − grad grad (I1 )
  E
]
  1 +ν  ν
[ ]
T T
   
 grad  div ε   =  grad  div σ   − grad grad (I1 )
   E    E

Enfin nous pouvons utiliser les équations d’équilibre : () r

div σ + ρ f = 0

L’équation de compatibilité donne les équation de BELTRAMI :

( ) [ ( ) ]I = 0
T
       r
grad div σ  + grad div σ   − 1 grad grad ( I1 ) − ∆σ +  ν  div ρ
 f
  1 +ν  1 −ν 
 
   

Sous forme indicielle on obtient :

∂ [ ρ ( fi ) ] ∂ [ ρ ( f j ) ]
∂x j
+
∂xi
− 1 ∂ 2 I1
1 + ν ∂xi ∂x j
− ∆ σ ij + ν
1+ν
r
[ ( )]
div ρ f δij = 0
Cas particulier : Dans le cas d’un domaine à masse volumique
constante, en équilibre etr placé dans un champ de force volumique
( r
dérivant d’un potentiel f = − grad p ⇒ div f = 0 )
, on obtient :

1
1+ν
grad (grad (I 1 ))+ ∆σ = 0

1 ∂ 2 I1
+ ∆σ =0
1 + ν ∂xi ∂x j
ij

Les conditions aux limites :

•Les efforts imposés : ( )
T N, n = F ( N ) ∀N ∈ ∂DF

•Les déplacements imposés : u( N ) = d ∀N ∈ ∂Dd

En résolvant les équations de BELTRAMI, on détermine les

contraintes et on utilise les conditions aux limites en contraintes pour
lever les indéterminations restantes.
 Equations de compatibilité :
Equations d’équilibre T
   
divσ + ρfv = 0 grad div ε  +  grad div ε   − grad  grad  tr ε    − ∆ ε = 0
           


Equations de BELTRAMI :

( )
T
    
grad divσ  +  grad divσ   − 1 grad grad( I1 ) − ∆σ +
 ν 
 [ ( ) ]I =0
r
 divρ f
  1 +ν ν
   1 − 
   

Conditions de
symetrie et les Proposi-
conditions aux tion de σ
limites en (x).
contraintes.

σ(x) est élastostatiquement acceptable

Calcul de :
(
ε = 1 . gradU + grad TU
2
)
1 +ν
ε= σ −  ν tr σ   I
E E   

Conditions aux
u( N ) = d
limites en
déplacements

Toutes les conditions aux limites aux limites sont vérifiées

C’est la solution
 Processus de résolution par
l’approche en contraintes

• Postuler la forme du champs de contrainte.

• Verifier les conditions d’equilibre en volume.
• Vérifier les conditions aux limites en
contraintes.
• Vérifier les équations de Beltrami.
• Calculer les déformations puis les
déplacements.
• Vérifier les conditions aux limites en
déplacements.
 TD2 : EQUATION DE BELTRAMI
Un barrage est constitué par un mur vertical à loi de comportement
élastique linéaire, de largeur 2h et de hauteur L. Il est défini dans le repère
(O; X , X , X ) selon la figure ci-dessous. L’axe (O; Xr 1 ) est vertical descendant.
r r r
1 2 3

La face supérieure est notée Ss, la face inférieure en appui sur le sol Si et les
deux parois verticales Sh et S-h. Sur toute la hauteur de la paroi Sh s’exerce la
poussée d’un fluide de poids volumiqueϖ = ρ g selon l’axe (O; X .) Les forces de
r
2

volume dans le corps seront négligées. On considérera un état plan de

déformation selon l’axe (O; X ) c’est à dire un état déterminé par ε 13 = ε 23 = ε 33 = 0.
r
3

On propose un champ du tenseur des contraintes dans le barrage de

la forme:
r
X2
 ϖ x13 x 2 ϖ  6 2 3
σ
 11 = 3
+ 3 
h x x
1 2 − 2 x x
1 2
h h
 4 h 4 h  5 
L  ϖ x1  x 23 3 x 2 
σ 22 = − + ϖ x1  3 − 
 2  4 h 4 h 
3ϖ x1 2 ϖ 4 4 3ϖ 2

[ ] [ ] [ ]
2
σ
 12 = h − x 2
2 − h − x 2 + h − x 2
2
 8 h3 8 h3 20 h
r
X1
r
X2 1-Ce tenseur des contraintes est-il élastostatiquement
acceptable ?
h h
2-Vérifier que les conditions aux limites sur les deux
L parois verticales Sh et S-h sont bien respectées.
3-Si la face supérieure Ss est libre, montrer que ce
tenseur ne rend pas exactement compte des conditions
aux limites sur la face mais que l’effort tranchant
r résultant des contraintes de cisaillement est nul.
X1
Le liquide est de l’eau de poids volumique ϖ = 9810 N / m .
3

Les dimensions du barrage sont h = 0,5 m et L = 5 m. Le

coefficient de Poisson du matériau est υ. = 0,25 .
4- Tracer le tricercle de Mohr de l’état de contrainte pour les points M1(L,0) et
M2(L,h).
5- En quel point est situé la contrainte de cisaillement maximale et quelle est sa
valeur ?
 x1
TD3 :POUTRE EN FLEXION

O’
x3
2a O
x2
2a
L

Section S0 S’0 S1 S’1 S2 S’2

position X3 =0 X3 =L X1=a X1=-a X2 =a X2 =-a

Soit une poutre longue de longueur ‘L’ et de section carrée de coté

‘2a’. Le matériaux de la poutre est de module d’Young E et de
coefficient de Poisson ‘v’. La poutre est encastrée en S0
On suppose qu’il existe une fonction F(x1 , x2 , x3 ), telle que les
composantes du tenseur des contraintes sont:
σ 11 = σ 22 = σ 12 = 0 ; σ 33 = m.x1.x3 − M .x2 .x3

σ = ∂ F m . x 2
∂F M .x22
− 1
; σ 23 =− +
 13
∂x2 2 ∂x1 2
Les forces de volume sont néligeables.
1- Montrer que pour que la forme proposées de σ vérifie les équations
d’équilibre , il faut que la fonction F(x1 , x2 , x3 ) soit tel que:
F(x1 , x2 , x3 ) = G(x1 , x2 ) + H(x3 )
2- On donne maintenant à F la forme:
F = − (x12 + x22 ) + (α .x12 + β .x1.x22 )
C D
2 d
dans laquelle d est une longueur donnée et C, α et β sont trois
constantes.
Donner les expressions auxquelles doivent satisfaire α et β en fonction
des autres données ou constantes du problème pour que le champs de
contraintes soit élastostatiquement acceptable.
3- Si la surface (S1 ) est libre , calculer m et M. En déduire alors l’état de
chargement sur les faces (S’1 ) , (S2 ) et (S’2 ) en fonction du paramètre
de charge C. Expliciter les résultats à l’aide de schémas.
4-Quelles sont en O la résultante R et le moment M0 des efforts dus à
l’encastrement ?
5- En supposant que les conditions d’encastrement satisfassent à
()
rot u O = 0 et ()u O = 0 , déterminer le champ du vecteur déplacement
u(x).