Dimensionnement d’un barrage 1

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. Ce problème est largement inspiré de l’examen posé par M. Amestoy en 1991, cf. (Amestoy, 2004) qui
met en valeur la condition de Lévy (1895). On trouvera également une discussion de ce calcul de barrage chez
Bellet (1990).
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

Figure 2: Le barrage en perspective cavalière : repère (X1 , X2 , X3 ), section OAB en X3 = 0. La longueur du

barrage est L.

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.
La résultante des efforts sur le parement amont est le vecteur :
Z Z
R OA = −pn dS = p(ue 1 + e 2 )dX1 dX3

2
avec p = ρeau gX1 . Finalement, R OA = ρeau gL H2 (ue 1 + e 2 ).
Le poids du barrage est
1
P = ρgH 2 (u + v)Le 1
2
La force tangentielle exercée par le barrage sur la fondation est donc donnée par la composante
selon la direction 2 de R OA :
H2
T = ρeau gL
2
La condition de non glissement est
ρeau
T < CP =⇒ <C
ρ(u + v) + ρeau u

où C est le coefficient de frottement de Coulomb.

Si u = 0, l’angle β doit être suffisamment grand :
ρeau
v>
Cρ
L’angle minimal obtenu pour C = 0.8 est β = 27◦ .

2. Etude du barrage dans le cas α = 0

Cette partie du problème est consacrée au cas d’un parement amont vertical : α = u = 0.

3
2.1. Conditions aux limites exhaustives
Donner explicitement l’ensemble des conditions aux limites du problème posé sur la section
plane OAB.
Par ailleurs, on commentera également les conditions sur les faces X3 = ±L/2.
Le problème est posé en conditions mixtes d’efforts et déplacements donnés :
• Sur OA, t = t d = −pn = ρeau gX1 e 2 .
• Sur OB, t = t d = 0.
• Sur AB, u1 = u2 = u3 = 0.
Les conditions en X3 = ±L/2 ne sont pas toutes explicitées. On a u3 = 0 mais rien n’est dit pour
les deux autres directions. Comme la longueur du barrage est prise suffisamment grande par
rapport à sa hauteur, on peut appliquer le principe de Saint-Venant et se contenter d’indiquer
le torseur des efforts selon ces deux directions. On peut en particulier négliger les résultantes
des forces selon 1 et 2 engendrées par l’interaction avec les blocs rocheux latéraux. Ces forces
induiront un effet de bord aux extrémités du barrage.

2.2. Conditions d’équilibre local

Exprimer les conditions d’équilibre en tout point intérieur de la section OAB en fonction des
composantes du tenseur des contraintes.
Il s’agit de la première loi de Cauchy exprimée en coordonnées cartésiennes dans le cas statique
et en présence d’efforts de pesanteur :

σ11,1 + σ12,2 + σ13,3 + ρg = 0, σ21,1 + σ22,2 + σ23,3 = 0, σ31,1 + σ32,2 + σ33,3 = 0 (1)

2.3. Forme proposée pour les contraintes

Dans le contexte des déformations planes, les composantes de contraintes recherchées sont

σ
∼
= σ11 e 1 ⊗ e 1 + σ22 e 2 ⊗ e 2 + σ12 (e 1 ⊗ e 2 + e 2 ⊗ e 1 ) + σ33 e 3 ⊗ e 3 (2)

où l’on a introduit les vecteurs de la base cartésienne orthonormée utilisée.

Ces composantes de contraintes sont recherchées de la forme suivante :

σ11 = aX1 + bX2 (3)

σ12 = cX1 + dX2 (4)
σ22 = eX1 + f X2 (5)
σ33 = h(X1 , X2 ) (6)

où les 6 constantes a, b, c, d, e, f caractérisant les profils linéaires recherchés ont été introduites,
ainsi que la fonction h de deux variables.
Ecrire les conditions d’équilibre local en fonction de ces constantes.
Ecrire les conditions supplémentaires portant sur ces constantes et résultant de l’application
des conditions aux limites sur OA et OB.
Résoudre enfin le système linéaire par rapport à a, b, c, d, e, f .
Constater que les contraintes σ11 , σ12 et σ22 trouvées ne dépendent pas des propriétés
élastiques du matériau constitutif du barrage.
Les deux premières conditions d’équilibre local fournissent :

a + d = −ρg, c+f =0

4
Les conditions aux limites sur OA, i.e. pour X2 = 0, fournissent, quant à elles :

c = 0, e = −ρeau g

Les conditions sur OB, i.e. pour X2 = vX1 , donnent

− σ11 sin β + σ12 cos β = 0, −σ12 sin β + σ22 cos β = 0 (7)

Elles deviennent :
a + bv − d = 0, −dv 2 + e = 0
Finalement, on obtient :

a = −ρg + ρeau gv −2 , b = ρgv −1 − 2ρeau gv −3 , c=0

d = −ρeau gv −2 , e = −ρeau g, f =0
Ces calculs montrent que l’analyse statique du problème du barrage suffit à déterminer
entièrement les contraintes qui s’y établissent. Il faut toutefois terminer la résolution du
problème avant de statuer sur la pertinence de la forme proposée.

2.4. Calculer finalement σ33

Il faut pour cela utiliser la loi de Hooke et la condition de déformation plane :
1+ν ν
ε33 = σ33 − trace σ
∼
=0
E E
On en déduit
σ33 = ν(σ11 + σ22 ) = ν((a + e)X1 + (b + f )X2 ) (8)

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 ?
Les composantes de déformation sont également linéaires par rapport aux coordonnées X1 et
X2 , en vertu de la loi de Hooke :

ε11 = AX1 + BX2 , ε12 = DX2 , ε22 = F X1 + GX2 (9)

avec
(1 − ν 2 )a − ν(1 + ν)e (1 − ν 2 )b 1+ν
A= , B= , D= d
E E E
(1 − ν 2 )e − ν(1 + ν)a ν(1 + ν)b
F = , G=−
E E
La linéarité du champ assure que les conditions de compatibilité sont toutes satisfaites.

5
2.6. Champ de déplacements
Déterminer les composantes du champ de déplacements.
On ne cherchera pas à ce stade à identifier les constantes d’intégration.
On construit le champ de déplacement par la méthode systématique disponible dans le contexte
cartésien orthonormé. On commence par déterminer les rotations infinitésimales en intégrant
la relation ωij,k = εik,j − εjk,i . On trouve :

ω12 = BX1 + (D + F )X2 , ω23 = ω31 = 0

auxquelles il faut ajouter des rotations infinitésimales de corps rigide.

Le déplacement s’obtient alors en intégrant ui,j = εij + ωij :

A 2 1 B G
u1 = X1 + BX1 X2 + (2D − F )X22 , u2 = − X12 + F X1 X2 + X22 (10)
2 2 2 2
auxquels s’ajoutent les translations et rotations infinitésimales de corps rigide.

2.7. Respect du contexte infinitésimal

En examinant successivement les composantes ε12 et ε11 , établir deux conditions nécessaires
au respect du contexte infinitésimal en fonction des caractéristiques du barrage.
On vérifiera que ces conditions sont remplies avec les valeurs considérées dans le problème.
On doit garantir que
 
1+ν 1 + ν H 
|ε12 | = |DX2 | ≤ |DHv| = |d|Hv =  ρeau g   1
E E v

qui fournit la condition nécessaire :

H 2E

ρeau g
v 3
Cette relation indique que le barrage ne doit pas être trop haut à cause de la pression exercée
par l’eau, et pas trop mince bien entendu. La condition v  3ρeau gH/2E ' 3 × 10−5 est
largement garantie même avec des valeurs très petites de β.
On considère ensuite |ε11 | = |AX1 +BX2 | ≤ (|A|+|Bv|)H. En majorant généreusement chacun
des termes, on obtient la condition :

(ρg − 2ρeau g/v 2 )H

1
E
qui exige que le poids du barrage soit suffisant. Ce terme est de l’ordre de 10−5 pour un angle
β de 40◦ .

2.8. Qualité de la solution approchée

Montrer finalement que le champ de déplacement trouvé n’est pas acceptable.
Les contraintes et déformations trouvées représentent toutefois une bonne approximation de la
solution du problème. Pour le voir, des calculs par éléments finis ont été réalisés pour obtenir
une solution numérique du problème aussi précise qu’on le souhaite, pour des caractéristiques
particulières du barrage.
Commenter les cartes de contraintes de la figure 3.

6
 (a) (b)

-1.8 -1.5 -1.2 -0.9 -0.6 -0.3 0 -0.8 -0.6 -0.4 -0.2 0

(c) (d)

-0.4 -0.3 -0.2 -0.1 0 -0.5 -0.42 -0.33 -0.25 -0.17 -0.085 -0.002

(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◦ .

Les contraintes, déformations et déplacements sont tracés le long des lignes des parements
amont et aval sur la figure 4. Discuter la qualité de la solution approchée du problème du
barrage construite précédemment.
Pour identifier complètement le champ des déplacements, il reste à vérifier les conditions
d’encastrement en X1 = H. C’est là que le bât blesse. En effet, s’il est possible de fixer
la translation de corps rigide de telle sorte que le déplacement du point A soit nul, les
rotations infinitésimales de corps rigide ne permettent pas d’annuler les termes quadratiques
des expressions (10) tout le long de la ligne X1 = H. C’est une profonde déception. Toutefois,
il s’avère que la solution trouvée est de bonne qualité. Le principe de Saint-Venant ne peut pas
être invoqué ici pour justifier ce fait car la zone d’encastrement ne possède pas de dimension

7
 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.0005
0

-5e-06
0
-1e-05
déplacement (m)

-1.5e-05 -0.0005
déformation

-2e-05
-0.001
u1 analytique
-2.5e-05 ε11 analytique u1 EF
ε11 EF u2 analytique
ε22 analytique u2 EF
-3e-05 -0.0015
ε22 EF
-3.5e-05
-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) ; (b) 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.

négligeable devant les autres. Seule une analyse numérique permet dès lors de trancher, au
moins pour quelques configurations particulières du barrage (c’est l’inconvénient de l’approche
numérique, même si une étude numérique paramétrique est bien sûr possible).
La déformée de la figure 3(b) illustre bien les effets conjoints du poids du barrage lui-même et
de l’action de l’eau par rapport à la configuration de référence dans un état naturel 3(a). Les
cartes de contraintes de la figure 3 montrent que les composantes σ11 , σ12 et σ22 sont quasi–
linéaires dans un voisinage relativement grand des points O et B. En particulier, au voisinage
de O, la composante σ22 (resp. σ12 ) est proportionnelle à la coordonnée X1 (resp. X2 ), comme
attendu d’après la solution trouvée dans les paragraphes précédents : σ12 = dX2 , σ22 = eX1 .
La différence principale entre les solutions analytique et numérique réside dans la forte
concentration de contrainte située au point A, qui perturbe le champ linéaire dans un large
voisinage du point A. Une singularité des contraintes est présente en A, en conséquence de
la condition d’encastrement de A alors même qu’une pression s’applique sur le parement aval
dans son voisinage. L’existence de cette singularité peut être mise en évidence grâce au calcul
par éléments finis en utilisant des maillages de plus en plus raffinés autour du point A. La

8
figure 5 donne le profil de σ11 le long du parement amont pour deux discrétisations différentes.
On voit que la contrainte σ22 augmente et devient positive lorsqu’on s’approche suffisamment
près de A. La contrainte de traction n’a pas de limite en A. Cela signifie d’une fissure apparaı̂t
nécessairement à cet endroit. La singularité mérite une étude théorique spécifique. Elle peut
être relâchée dans la pratique en évitant la présence d’un angle droit en A. Au contraire, on
n’observe pas de singularité en B.
Les courbes de la figure 4 confirment la bonne performance de la solution analytique proposée.
Celle–ci constitue une bonne approximation de la solution réelle excepté au voisinage du point
A et du point B. Les contraintes le long du parement aval sont surestimées, en valeur absolue,
dès qu’on s’éloigne du point O.
La figure 4(d) montre que le déplacement u1 coı̈ncide avec le champ trouvé pour lequel le
mouvement de corps rigide a été fixé à zéro. Par contre la composante u2 trouvée numériquement
n’est pas une parabole centrée en O. Elle présente un terme linéaire en X1 supplémentaire qui
correspond à la rotation infinitésimale rX1 qui n’a pas pu être fixée par les conditions aux
limites dans la solution analytique.
5
50 éléments
70 éléments
4

2
σ11 (MPa)

-1

-2
0 20 40 60 80 100
H − X1 (m)

Figure 5: Profil de la contrainte σ11 le long du parement amont pour deux discrétisations par éléments finis
différentes. Le maillage raffiné à 70 éléments possède plus d’éléments proches de A et des valeurs plus grandes
en traction sont atteintes.

2.9. Critère de Mohr–Coulomb

Le béton est un matériau fragile dont le domaine d’élasticité peut être décrit par le critère de
Mohr–Coulomb :
f (σn , τ ) = |τ | + σn tan φ − C < 0 (11)
où σn et τ sont respectivement la contrainte normale et la contrainte tangentielle sur une
facette de normale donnée au point considéré. Ce critère d’endommagement ou de rupture fait
intervenir deux paramètres : la cohésion C et l’angle de frottement interne φ. Ces paramètres
sont déterminés expérimentalement. Tant que f (σn , τ ) < 0 pour toutes les facettes possibles,
la réponse du matériau reste élastique linéaire. La limite d’élasticité est atteinte lorsque, dans
le plan (σn , τ ), dit plan de Mohr, le grand cercle de Mohr vient tangenter la droite d’équation
τ + σn tan φ − C = 0 (avec τ > 0).

9
 Comparer ce critère aux critères de plasticité et/ou de rupture que vous connaissez.
Calculer la limite d’élasticité en tration selon ce critère en fonction des paramètres du modèle.
Interpréter ce critère à la lumière de ce qu’évoque pour vous le nom de Coulomb.
Exprimer le critère de Mohr–Coulomb uniquement en fonction des contraintes principales
σ1 ≥ σ2 ≥ σ3 .
Tracer les cercles de Mohr associés au tenseur des contraintes au point X1 = H, X2 = 0 à
la base du parement amont en utilisant la solution analytique trouvée précédemment. Pour
cela, on calculera les valeurs numériques des composantes de contraintes à l’aide des formules
établies au paragraphe 2.3, et des données suivantes : u = 0, v = 2, H = 100 m.
Dans ce même diagramme, tracer la frontière du domaine d’élasticité décrit par la fonction
(11). 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 vis–à–vis du critère de Mohr–
Coulomb.
Le critère proposé fait intervenir l’effet conjoint de la contrainte tangentielle et de la contrainte
normale pour une facette donnée. Lorsque tan φ = 0, il s’identifie avec le critère de plasticité de
Tresca. A l’inverse, si tan φ devient très grand, le critère s’approche d’un critère de contrainte
normale critique, typique de la rupture des matériaux fragiles. Il est représenté par une droite
dans le plan de Mohr, dont l’inclinaison est donnée par φ. L’intersection avec l’axe des abscisses
fournit la limite de rupture en traction : σn = C/ tan φ. La combinaison des contraintes normale
et tangentielle dans le critère (11) fait irrésistiblement penser au critère de frottement de
Coulomb entre corps matériels. Du point de vue physique, il s’agit du frottement intérieur
entre les grains d’un matériau granulaire tels qu’on les rencontre dans le domaine du génie
civil. La résistance au cisaillement est d’autant plus grande que la compression est forte. Le
béton est un matériau granulaire dont la limite d’élasticité est atteinte par des effets combinés
de contraintes normale et de cisaillement.
Le grand cercle de Mohr a pour centre (σ1 + σ3 )/2 et pour rayon (σ1 − σ3 )/2. Il est tangent à
la droite de Mohr–Coulomb lorsque
σ1 − σ3
tan φ =
2C/ tan φ − (σ1 + σ3 )
ce qui fournit l’expression alternative du critère de Mohr–Coulomb :

g(σ1 , σ3 ) = σ1 − σ3 + (σ1 + σ3 ) tan φ − 2C (12)

Les contraintes le long du parement amont dans le cas u = 0, v = 2 valent :

σ11 = aX1 , σ12 = 0, σ22 = eX1 , σ33 = ν(a + e)X1

avec a = (ρeau /4 − ρ)g et e = −ρeau g. Le tenseur des contraintes y est donc diagonal et les
contraintes principales valent, dans l’ordre décroissant, en X1 = H = 100 m :
σ1 = σ33 = −0.64 MPa, σ2 = σ22 = −0.98 MPa, σ3 = σ11 = −2.21 MPa. On constate qu’elles
sont toujours négatives. C’est d’ailleurs un résultat qui s’étend aux contraintes principales dans
tout le barrage si l’on en croit les résultats du calcul par éléments finis illustrés par la figure 6, à
l’exception toutefois du voisinage du point singulier A où une étude spécifique est indispensable.
Etant donné la faible cohésion du béton (limite d’élasticité en traction C/ tan φ de l’ordre de
0.1 MPa), les contraintes de traction doivent être évitées absolument, c’est une consigne de
dimensionnement en génie civil.
Les cercles de Mohr correspondant aux valeurs précédentes des contraintes principales sont

10
représentées sur la figure 7. La frontière du domaine d’élasticité est la droite d’équation
f (σn , τ ) = 0. Elle coupe l’axe des abscisses en σn = C/ tan φ, i.e. la limite d’élasticité en
traction simple. On observe que les cercles de Mohr sont situés sous cette droite, dans le
domaine f (σn , τ ) < 0 assurant un fonctionnement élastique du point matériel considéré.
Plus généralement, Albert Caquot (1881–1976), professeur de Matériaux à l’Ecole des Mines
de Paris à partir de 1922, a proposé de considérer une courbe intrinsèque plus générale que la
droite de Mohr–Coulomb, issue des observations expérimentales, pour caractériser le domaine
d’élasticité du béton. La frontière du domaine d’élasticité est atteinte lorsque le grand cercle
de Mohr vient tangenter la courbe intrinsèque de Caquot. C’est le critère de Mohr–Caquot.
Pour finir, il faut indiquer qu’un critère d’endommagement supplémentaire pour les fortes
valeurs négatives de σn est nécessaire. La rupture alors a lieu sous forme de bandes de
localisation de l’endommagement conduisant à de multiples fissures.

-0.5 -0.42 -0.34 -0.26 -0.18 -0.1 0 -0.8 -0.67 -0.54 -0.41 -0.28 -0.15 -0.02

Figure 6: Champs de deux contraintes principales : σ1 à gauche et σ2 à droite. Ils montrent que les contraintes
principales sont partout négatives, sauf pour σ1 au voisinage du point A, d’après la figure 5. Les valeurs sont
données en MPa. Les caractéristiques du barrage sont H = 100 m et u = 0, v = 2, correspondant à un angle
β = 63◦ .

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.

3.1. Détermination des contraintes

Rappeler les conditions aux limites sur les parements amont et aval en tenant compte du fait
que α 6= 0.
On recherche à nouveau un champ de contraintes de la forme donnée par les équations (3) à
(6).
Etablir le système d’équations permettant de déterminer de manière unique les constantes
a, b, c, d, e, f . Pour gagner du temps, on ne demande pas de résoudre ce système. Les solutions

11
 τ (MPa)
2

1.5

0.5

σ3 0
σn (MPa) σ2 σ1
-2 -1.5 -1 -0.5 0

Figure 7: Cercles de Mohr des contraintes en X1 = H, X2 = 0 ; comparaison avec le critère de Mohr–Coulomb.

en sont données ci-dessous :

u2 + v 2 u + 2uv 2 − v u−v u2 + 3uv − 2
a = −ρg − ρeau g , b = −ρg + ρeau g (13)
(u + v)2 (u + v)3 (u + v)2 (u + v)3
uv(u − v) uv(uv − v 2 − 2) 2uv u + 2uv 2 − v
c = −ρg + ρeau g , d = −ρg + ρ eau g (14)
(u + v)2 (u + v)3 (u + v)2 (u + v)3
2u2 v 2 v 2 (2u2 v − 3u − v)
e = −ρg + ρeau g , f = −c (15)
(u + v)2 (u + v)3

Les conditions sur OB (surface libre où X2 = vX1 ) sont toujours données par les équations (7).
Sur la surface OA d’équation X2 = −uX1 , le vecteur contrainte est imposé :

−σ11 sin α − σ12 cos α = p sin α, −σ12 sin α − σ22 cos α = p cos α

où p = ρeau gX1 . Ces conditions fournissent les 4 équations suivantes :

(a + bv)v − (c + dv) = 0, (c + dv)v − (e + f v) = 0 (16)

(a − bu)u + (c − du) + ρeau gu = 0, (c − du)u + (e − f u) + ρeau g = 0 (17)

Les deux dernières équations résultent des conditions d’équilibre local (1) :

a + d + ρg = 0, c+f =0 (18)

On peut vérifier que la résolution du système précédent conduit aux valeurs (13) à (15).

12
3.2. Contraintes et directions principales sur le parement amont
Déterminer les directions principales (vecteurs propres du tenseur des contraintes) puis les
contraintes principales (valeurs propres) des points matériels du parement amont OA. Cette
détermination ne nécessite pas de calcul fastidieux.
La forme (2) montre que le vecteur e 3 est une direction principale des contraintes. La contrainte
principale correspondante est
σIII = σ33
Sur le parement amont, les contraintes sont de la forme :

σ11 = (a − bu)X1 , σ12 = (c − du)X1 , σ22 = (e − f u)X1

La condition aux limites sur le parement amont indique que le vecteur normal
n = − sin α e 1 − cos α e 2 est un vecteur propre de σ
∼
pour la valeur

σI = −p = −ρeau gX1

La troisième direction principale est donc m = n ∧ e 3 = − cos α e 1 + sin α e 2 . La valeur

principale associée s’obtient à partir des deux autres et de la trace du tenseur des contraintes :

σII = σ11 + σ22 + σ33 − σI − σ33 = σ11 + σ22 + ρeau gX1 (19)

Par suite,
σII = (a − bu + e − f u + ρeau g)X1 (20)

O
X2

X1

Figure 8: 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.

13
3.3. Présence d’une fissure perpendiculaire au parement amont
On considère que, en raison d’un défaut ou de l’endommagement du béton, une petite fissure
plane apparaı̂t perpendiculairement au parement amont, à une certaine cote X1 = X0 donnée,
comme sur la figure 8. Il s’agit d’une amorce de fissure dont la profondeur est négligeable devant
les caractéristiques du barrage. La fissure est perpendiculaire au parement aval et parallèle à
la direction 3. On se doit d’étudier la nocivité d’un tel défaut.
La sécurité du barrage sera assurée si le chargement mécanique a pour effet de refermer
cette amorce de fissure. Formuler cette condition en utilisant les contraintes trouvées le long
du parement amont. On supposera d’abord que l’eau ne s’infiltre pas dans la la fissure.
Formuler ensuite cette condition en considérant que l’eau a pénétré par la fissure et que le
fluide interstitiel exerce une pression sur les lèvres de la fissure.
Les critères seront exprimés en fonction des constantes a à f déterminées précédemment et
des données du problème.
Montrer que les critères obtenus ne font pas intervenir la position X0 de l’amorce de fissure.
Le champ de contrainte le long du parement n’est pas modifié tant que la fissure reste fermée
puisque les lèvres en contact transmettent l’intégralité des charges. La contrainte agissant sur
le plan de la fissure est la contrainte principale σII déterminée au paragraphe 3.2. La fissure
reste fermée tant que cette contrainte principale est négative :

σII ≤ 0

Si l’eau s’infiltre à travers la fissure et exerce une pression p = ρeau gX0 tendant à ouvrir la
fissure, il faut que la contrainte principale σII puisse contrecarrer cet effet d’ouverture :

σII + p ≤ 0 (21)

C’est une condition plus sévère que la précédente. En utilisant le résultat (20), le critère de non
traction sur la fissure s’écrit
a − bu + e − f u + 2ρeau g ≤ 0 (22)
On remarque que la condition est indépendante de la cote X0 . Alors que les contraintes
augmentent, en valeur absolue, avec X1 , la condition précédente reste la même.

3.4. Condition de Maurice Lévy

Expliciter alors la dernière condition précédente à l’aide des quantités u, v, ρ et ρeau . Avec
l’hypothèse (bien sûr réalisée en pratique), ρ > ρeau , établir finalement son écriture sous la
forme : p
−uρ + u2 ρ2 + 4ρeau (ρ − ρeau )
v≥ (23)
2(ρ − ρeau )
dite condition de Maurice Lévy.
En tenant compte des équations (16) et (17), on note que
1
a − bu + ρeau = − (c − du), e − f u + ρeau = u(c − du)
u
En substituant ces relations dans la condition (22), on obtient
1
−(c − du)(u + ) ≤ 0
u

14
qui devient c − du ≥ 0 puisque α ∈ [0, π2 ]. On utilise alors les expressions de c et d données par
(14) pour obtenir la condition

ρ(u + v)v − ρeau (1 + v 2 ) ≥ 0

qui se met aussi sous la forme

v 2 (ρ − ρeau ) + ρuv − ρeau ≥ 0

Le discriminant ∆ = ρ2 u2 + 4(ρ − ρeau )ρeau est toujours positif dès que ρ ≥ ρeau . Le binôme
précédent est alors positif seulement si la condition cherchée (23) est également remplie (avec
v ≥ 0). Cette condition fournit un angle minimum β pour u fixé.

3.5. Optimisation du barrage par rapport à la condition de Lévy

On cherche alors à dimensionner au mieux le barrage, c’est–à–dire à respecter la condition de
sécurité précédente à l’aide d’une masse minimale, la hauteur H et la longueur L étant fixées.
En utilisant la condition, ρ > 2ρeau , remplie en pratique, trouver la valeur optimale du
paramètre u, puis celle du paramètre v.
Faire l’application numérique pour un barrage en béton renforcé, avec les propriétés
correspondantes indiquées au début du problème.
La masse du barrage vaut :
LH 2
m=ρ (u + v)
2
La condition de Lévy permet de trouver la masse minimale :
p p
LH 2 −uρ + u2 ρ2 + 4ρeau (ρ − ρeau ) LH 2 u(ρ − 2ρeau ) + u2 ρ2 + 4ρeau (ρ − ρeau )
m≥ρ (u+ )=ρ
2 2(ρ − ρeau ) 2 2(ρ − ρeau )

Lorsque ρ ≥ 2ρeau , le dernier terme est minimal pour

u=0

On obtient alors la masse minimale et l’angle minimal :

LH 2
r r
ρeau ρeau
mmin = ρ , vmin = (24)
2 ρ − ρeau ρ − ρeau

qui fournissent les caractéristiques optimales du barrage.

Dans le cas du béton renforcé, on trouve v = 0.82 et β ' 40◦ .

Références
Amestoy, M., 2004. Mécanique des milieux continus déformables : recueil de problèmes. Cours
de l’Ecole des Mines de Paris.
Bellet, D., 1990. Problèmes d’élasticité. Cepadues Editions, Toulouse, France.
Lévy, M., 1895. Quelques considérations sur la construction des grands barrages. Comptes
Rendus des séances de l’Académie des Sciences 121, 288–300.

15