Travail dirigé de mécanique des milieux continus
UPSSITECH Licence GCGEO - Niveau L3
Université Paul Sabatier - Toulouse 3 - 2014
0.1 méthode de Ritz - barrage
On étudie le barrage poids primatique de section triangulaire représenté sur la gure ci-dessus.
On cherche à établir les champs de contrainte et de déformation par la méthode de Ritz dans
un état plan de déformation. On suppose que les points de la base y = 0 sont bloqués (pas de
déplacement permis).
1. Montrer que la fonction de déplacement → −u (x, y) = Hy
(a→
−x + b→
−y ) est cinématique-
ment admissible.
Le barrage repose sur sa base dont l'équation est y = 0 ; de ce fait la fonction →
−
u est nulle le
long de la base et est donc cinématiquement admissible.
2. Exprimer le potentiel des forces de surfaces le long de [OB].
Le long de OB le potentiel des forces pour une longueur unité de barrage résulte de la
pression de l'eau. Ainsi : ˆ
V =− p p (y) →
−
x ×→
−u (0, y) .dy
OB
soit :
ˆ H ˆ H
y2
y 1
Vp = − $e (H − y) × ady = −$e a y− dy = − $e aH 2
y=0 H 0 H 6
3. Exprimer le potentiel des forces de volume.
Le potentiel gravitaire par unité de longueur de barrage vaut :
ˆ
Vg = − −$b →
−
y ×→ −
u (x, y) .dxdy
barrage
soit :
H (1− B
x
)
ˆB ˆ
b
Vg = $b [Link] dx
H
x=0 y=0
donc :
ˆB
b 1 2 x 2 $b bBH 1 1 1
Vg = $b H 1− dx = 1−2× + = + $b bBH
H 2 B 2 2 3 6
x=0
4. Exprimer le potentiel élastique.
Le potentiel élastique est donné par :
˚
1
Ve = σ : ε.dv
2
barrage
en déformations planes et pour une longueur unitaire de barrage :
¨
1
Ve = σ : ε.dS
2
barrage
il faut préalablement calculer les déformations :
∂u
εxx = =0
∂x
∂v b
εyy = =
∂y H
∂u ∂v a
γxy = + =
∂y ∂x H
notons que :
b
tr ε = Θ =
H
puis les contraintes :
b
σxx = 2µεxx + λΘ = λ
H
b
σyy = 2µεyy + λΘ = (2µ + λ)
H
a
σxy = µγxy = µ
H
d'où le potentiel élastique par unité de volume (elle est uniforme) :
b2 a2
1 1
ve = (σxx × εxx + σyy × εyy + σxy × γxy ) = (2µ + λ) 2 + µ 2
2 2 H H
et le potentiel élastique global par unité de longueur :
¨ ¨
1 1
(2µ + λ) b2 + µa2
Ve = ve .dS = ve dS = BH × 2
2 2H
barrage barrage
soit nalement :
1B
(2µ + λ) b2 + µa2
Ve =
4H
5. Exprimer la fonction potentiel V et calculer les constantes a et b.
Potentiel :
1 1 1B
V = Vp + Vg + Ve = − $e aH 2 + $b bBH + (2µ + λ) b2 + µa2
6 6 4H
Il est stationnaire à l'équilibre de sorte que
∂V = 0
∂a
∂V
=0
∂b
donc :
∂V = − 1 $e H 2 + 1 B µa = 0
∂a 6 2H
∂V 1 1B
= + $b BH + (2µ + λ) b = 0
∂b 6 2H
par voie de conséquence :
1 $e H 3 1 $b
a= et b=− H2
3 µ B 3 2µ + λ
6. Calculer le tenseur des déformations puis des contraintes.
Tenseur des déformations :
1 $e H 2
1 a
ε
xx = 0 ε xy = εxx = 0 ε xy =
ε= 2H 6 µ B
b =
1 $b
εxy εyy = εxy εyy = − H
H 3 2µ + λ
Tenseur des contraintes :
1 H2
1 λ$b
σxx = − 3 2µ + λ H σxy =
3
$e
B
σ=
1
σxy σyy = − $b H
3
7. Le champ de contrainte obtenu est-il statiquement admissible ? Conclusion.
Les conditions de l'équilibre local d'un champ de contrainte uniforme ne sont pas compatibles
avec des forces volumiques gravitaires. Ce tenseur ne convient pas pour représenter l'état de
contrainte réel dans le barrage. A noter que les conditions aux limites ne sont pas davantage
satisfaites.
Université de Toulouse - UPSSITECH GCGEO - [Link] - mmc L3 - travail dirigé.