Examen 2017
Examen 2017
Examen 2017
On se propose d’étudier, dans l’hypothèse des petites perturbations, l’équilibre d’un barrage
constitué d’un matériau élastique, homogène et isotrope. Plusieurs critères de dimensionnement
seront envisagés pour aboutir finalement à une optimisation des caractéristiques du barrage.
La géométrie de cette structure est décrite sur les figures 1 et 2. Il s’agit d’un barrage poids 2 de
forme prismatique de génératrice parallèle à la direction 3. Le repérage est cartésien orthonormé.
La direction X1 est la verticale orientée vers le bas dans le sens de la gravité. Les sections,
perpendiculaires à l’axe X3 , sont toutes équivalentes et la figure 1 montre la section centrale à
la cote X3 = 0, avec le repère de la figure 2. Le barrage est de hauteur H et de longueur L. La
forme triangulaire de la section est caractérisée par les angles α et β, choisis positifs. A toutes
fins utiles, on notera :
u = tan α, v = tan β
Le chargement est associé aux deux facteurs suivants : d’une part, la pesanteur propre
s’exerçant sur le matériau constitutif de masse volumique ρ ; d’autre part, l’effet de pression
dû à la mise en eau (masse volumique ρeau ) jusqu’à affleurement. La pression atmosphérique
agissant sur le parement aval est, quant à elle, négligée. L’état initial (hors pesanteur !) est
naturel et l’évolution isotherme.
On ne prendra pas en compte les effets des écoulements internes de l’eau par infiltration dans
les interstices du barrage, sauf mention du contraire 3 .
On considère un encastrement à la base du barrage (fondation), correspondant à des
déplacements donnés nuls en X1 = H. On suppose également qu’en raison de la raideur
importante des blocs rocheux latéraux sur lesquels s’appuie le barrage, et de la longueur L
suffisamment grande par rapport à H, les faces X3 = ±L/2 ne subissent aucun déplacement
latéral u3 = 0. Cette condition nous place dans le contexte des déformations planes, c’est–à–dire
que la condition u3 = 0 s’étend à tous les points du barrage et que les champs de contraintes
et de déformations ne dépendent pas de X3 .
Pour les applications numériques demandées dans le problème, on considérera que le
barrage est en béton renforcé de comportement élastique linéarisé, avec un module de Young,
E = 50 GPa, un coefficient de Poisson, ν = 0.2 et une masse volumique ρ = 2500 kg.m−3 . De
plus ρeau = 1000 kg.m−3 . Enfin, l’accélération de la pesanteur est supposée constante égale à
g = 9.8 m.s−2 .
On prendra la hauteur H = 100 m.
1. Quelques considérations sur la construction des grands barrages, par M. Maurice Lévy, Comptes Rendus
des séances de l’Académie des Sciences, vol. 121, pp. 288–300, 1895. Cet article est paru peu après la rupture
du barrage de Bouzey (avril 1895, département des Vosges).
2. par opposition aux barrages voûtes dont la forme n’est pas prismatique.
3. Leur modélisation représente une difficulté majeure dans le calcul des barrages poids. Cette mécanique
des milieux poreux n’est pas abordée dans ce problème.
1
O X2
β
H
α
A B
X1
Figure 1: Coupe bidimensionnelle du barrage poids étudié. Le parement amont (resp. aval) désigne la partie
OA (resp. OB) du barrage. Les lignes pointillées ondulées représentent l’eau du barrage. La hauteur du barrage
est H.
X2
O
X3
B
A
L
X1
2
1. Critère de non glissement du barrage
Le barrage poids transmet la poussée de l’eau à la fondation par cisaillement de celle-ci,
principalement par frottement. Le poids de l’ouvrage doit être suffisant pour assurer le non
glissement du barrage sur sa fondation. Cette condition représente une vérification fondamentale
à effectuer : le glissement sur la fondation est la première cause de rupture des barrages poids.
Calculer la résultante des forces agissant sur le parement amont.
Calculer le poids du barrage.
Calculer la force tangentielle résultante s’exerçant sur la fondation.
Exprimer la condition de non glissement en faisant intervenir le coefficient de frottement de
Coulomb. Il est question ici de la condition globale portant sur les résultantes des efforts. Une
étude du glissement local n’est pas envisagée dans ce problème.
Les résultats seront donnés en fonction des paramètres géométriques du barrage et de
l’accélération de la pesanteur.
Lorsque α = 0, en déduire une condition sur v. Evaluer cette condition pour un coefficient de
frottement égal à 0.8.
3
2.4. Calculer finalement σ33
2.5. Calcul des déformations
Calculer les composantes des déformations au sein du barrage engendrées par le champ de
contraintes précédemment déterminées.
S’agit–il d’un champ de déformations compatible ?
4
Dans ce même diagramme, tracer la frontière du domaine d’élasticité décrit par la fonction
(6). On prendra les caractéristiques : C = 0.1 MPa et φ =46◦ , typiques du béton.
Conclure sur l’admissibilité des contraintes locales en ce point A vis–à–vis du critère de
Mohr–Coulomb.
(a) (b)
(c) (d)
(e) (f)
Figure 3: Résultats de calculs par éléments finis du barrage : (a) maillage par éléments finis (état naturel
sans pesanteur), (b) déformée du barrage (déplacements × 10000), (c) champ de contrainte σ11 , (d) champ de
contrainte σ22 , (e) champ de contrainte σ12 , (f) champ de contrainte σ33 . Les valeurs sont données en MPa. Les
caractéristiques du barrage sont ici H = 100 m et u = 0, v = 2, correspondant à un angle β = 63◦ .
5
0
0
-0.2
-0.5
contrainte (MPa)
contrainte (MPa)
-1 -0.4
σ11 analytique
-1.5 σ11 EF -0.6 σ11 analytique
σ22 analytique σ11 EF
σ22 EF σ22 analytique
σ33 analytique σ22 EF
-2 σ33 EF -0.8 σ12 analytique
σ12 EF
-2.5 -1
0 20 40 60 80 100 0 50 100 150 200
X1 (m) X2 (m)
(a) (b)
0 0
-0.0002
-5e-06
-0.0004
-1e-05
-0.0006
déplacement (m)
-1.5e-05
déformation
-0.0008
-0.001 u1 analytique
-2e-05
u1 EF
-0.0012
-2.5e-05 ε11 analytique -0.0014
ε11 EF
-3e-05 ε22 analytique -0.0016
ε22 EF
-3.5e-05 -0.0018
-0.002
0 20 40 60 80 100 0 20 40 60 80 100
X1 (m) X1 (m)
(c) (d)
Figure 4: Résultats de calculs par éléments finis du barrage : (a) tracé des contraintes le long du parement amont
(droite d’équation X2 = 0) ; (c) tracé des contraintes le long du parement aval (droite d’équation X2 = 2X1 ) ; (c)
tracé des déformations le long du parement amont ; (d) profils des déplacements le long du parement amont (la
translation de corps rigide a été fixée de telle sorte que le déplacement du point O soit nul). Les caractéristiques
du barrage sont les mêmes que pour la figure 3.
6
3. Etude du barrage dans le cas général
On revient au cas général de la figure 1 où, a priori, u 6= 0. L’objectif de cette partie est
d’optimiser les valeurs de u et v au regard du critère de Maurice Lévy associé à la présence
d’une fissure dans le parement amont.
7
O
X2
X1
Figure 5: Zoom sur une fissure perpendiculaire au parement amont OA à une certaine cote X1 = X0 . La fissure
est représentée légèrement ouverte pour l’illustration.
* * *