17

dilatance, dissipation d’énergie

et critère de rupture tridimensionnel
sous grandes déformations dans les matériaux granulaires
dilatancy, energy-dissipation
and tridimensional sailure criterion
under large strains in granular materials
E. FROSSARD
Ingénieur, Coyne & Bellier*

Résumé
Cet article propose un schéma physique simple décrivant les aspects essentiels du
comportement irréversible. A partir d'une analogie formelle entre le critère de
Coulomb et les relations contraintes-dilatance de Rowe est mis à jour un principe
énergétique élémentaire, sousjacent à ces équations, établies en symétrie axiale ou
en déformation plane.
La traduction objective de ce principe mène alors à une équation de dissipation,
tridimensionnelle. On détaille des résultats expérimentaux variés qui portent la trace
de cette équation de dissipation.
On examine les conséquences théoriques de cette équation de dissipation sur le
critère de rupture tridimensionnel sous grandes déformations, à l’état critique. La
pyramide de Coulomb apparaît alors comme minimum de résistance et comme lieu
des chemins de moindre énergie.

Abstract
A simple physical scheme is proposed in the paper, which describes essential
features of irreversible behaviour. Proceeding from a formal analogy between
Coulomb's criterion and Rowe's stress-dilatancy relations, an elemental energetic
principle appears, underlying these equations, which are established for axisymetric
conditions or place strain.
The objective expression of this principle leads to an energy-dissipation equation
valid for tridimensional conditions. Various experimental results following the line of
this energy-dissipation equation are detailed.
Theorical consequences of this energy-dissipation equation on critical state failure
criterion are examined.
Coulomb's pyramid then appears both as a strengh minimum and as a locus of
minimum energy paths.

* 5, rue d’Héliopolis. 75017 Paris.

18 N° 34 REVUE FRANÇAISE DE GÉOTECHNIQUE

INTRODUCTION Coulomb est largement reconnue, bien qu’il soit néces­

saire de faire apparaître des angles de frottement
Une certaine approche énergétique du comportement internes distincts pour ces trois situations, dès que le
mécanique des sols a été développée, pour les maté­ matériau exhibe une certaine dilatance ; cependant,
riaux fins, par les auteurs du modèle CAM-CLAY et on considère ici les conditions à volume constant de
de ses dérivés. l'état critique.
Pour les matériaux granulaires, et dans un autre ordre Les relations de Rowe ont été établies à l’origine dans
d’idées, ROWE (1962) et HORNE (1965) ont obtenu le cas de la symétrie de révolution, et plus tard
un résultat très convaincant à partir d’une modélisation étendues au cas de la déformation plane sur une
microstructurale et d’un principe de minimisation base expérimentale (Rowe (1969)). Ces relations sont
d’énergie, en symétrie de révolution. L’usage d’un en principe valides pour toute déformation à volume
principe de minimisation d’énergie, parfois discuté, a variable.
conduit de JOSSELIN DE JONG (1976) à proposer
une autre origine aux intrigantes équations de ROWE, 1.1. Les cas de symétrie axiale
fondée cette fois sur les seules lois de la friction,
toujours dans le cas de la symétrie de révolution en et de déformation plane
contraintes. - Considérons d’abord le cas du « triaxial » compres­
L’objet de cet article est d’aller plus loin dans cette sion. Les contraintes de compression et déformations
approche physique énergétique, sur les chemins tridi­ en contraction sont comptées positivement, la coaxialité
mensionnels, et en la reliant à la notion de critère des tenseurs o et ε est supposée réalisée, on se place
de rupture. dans le repère formé par les directions principales,
numérotées dans l’ordre décroissant des contraintes
Le schéma original qui est présenté ici permet de principales σ1 > σ2 > σ3. Les vitesses de déformation
décrire clairement nombre de caractères essentiels du principale, ainsi numérotées, n’ont donc pas d’ordre
comportement de ces matériaux. Comme il ne défini.
comporte qu’un seul paramètre intrinsèque, il est clair —Les matériaux granulaires ne résistant guère à la
qu’il ne constitue pas un modèle complet. traction, on suppose que les contraintes principales
sont toutes des compressions.
Le développement des méthodes de calcul numérique • Le critère de Coulomb
en mécanique des sols a suscité la floraison de modèles
de comportement, parfois sophistiqués, qui tentent (σ1 - σ3) = sin ϕ (1)
d’en saisir l’intégralité.
s'écrit aussi :
Il semble à l’usage que certains calculs d’ouvrages σ1/σ3, = tg2 (ji/4 + ϕ/2)
réels pourraient souvent être aussi bien conduits avec
des modèles simplifiés, plus proches de la caricature
efficace que de la description minutieuse, dont la
précision est souvent rendue superflue par les incer­
titudes et inconnues qui peuvent subsister jusqu’à un
stade avancé de la construction.
C’est dans cet esprit de recherche de schémas simples
que s’inscrit cet article.

1. UNE ANALOGIE FORMELLE

ET SA TRADUCTION PHYSIQUE
— Dans ce chapitre est détaillé le parallélisme qui relie
le critère de rupture de Coulomb aux relations con­
traintes dilatance de Rowe, ainsi que la signification
physique précise que l’on peut en tirer.
Le critère de Coulomb est une relation portant sur
les seules contraintes, tandis que les relations de Rowe
font intervenir à la fois les contraintes et vitesses de
déformation ; aussi, afin d’établir ce parallélisme, consi­
dérera-t-on le critère de Coulomb lorsque le matériau
atteint un régime stationnaire de déformation à volume
constant comme cela se produit après de larges
déformations, à l’état critique.
—Dans les conditions du « triaxial » axisymétrique en
compression, de la déformation plane, ou du « triaxial »
axisymétrique en extension, la validité du critère de Essai triaxial en compression.
DILATANCE. DISSIPATION D'ÉNERGIE ET CRITÈRE DE RUPTURE TRIDIMENSIONNEL 19

Dans les conditions de volume constant, cette dernière Le parallélisme est frappant, non seulement entre les
expression peut être transformée, en tenant compte colonnes du tableau I, mais aussi entre les lignes. Phy­
de la condition ε v = ε1+ 2ε3 = 0, qui s’écrit aussi siquement, ces expressions indiquent que le rapport de
ε1./2ε3 = - 1 certaines puissances développées par les contraintes est
constant.
On obtient ainsi par multiplication : —Examinons d’abord les constantes :
σ1, ε1/2σ3 ε3 = - tg2 (π/4 + ϕJ2) où est l’angle
de frottement interne à volume constant de l’état • d’après les résultats de CORNFORTH (1964), repris
critique. par LAMBE et WITHMAN (1969). les frottements internes
à volume constant seraient très voisins pour les trois
• La relation de Rowe s’écrit ici : chemins, à l’état critique.
σ1, ε3/2σ3 ε3 = tg2 (π/4 + ϕ/2) ϕ cvl= ϕ cv2 = ϕCV3
• d’après les résultats de HORNE (1969) et ceux de
où ϕ1, est une constante expérimentale, intermédiaire l’auteur (1983), l’angle de frottement interne à volume
entre l’angle de friction physique grain sur grain ϕµ constant ϕcv au triaxial compression excède l’angle
et l’angle ϕcV. Les valeurs mesurées pour ϕµ, sont de friction physique ϕµ de quelques degrés. Or ϕ1, est
voisines de ϕ dans le cas du triaxial compression sur compris entre ϕcv et ϕµ. La distinction entre tous ces
matériau dense. angles suppose une précision qui peut paraître un
peu superflue dans l’optique de certains calculs d’ou­
— Les situations de la déformation plane et du « triaxial vrages. L’auteur propose donc de remplacer tous ces
extension » s’étudient sans peine de façon similaire. angles par une seule constante, notée i)> et appelée
En tenant compte des particularités de ces diverses « angle d’écoulement ».
situations, on obtient le tableau suivant : 0 < ψ < π/2

Tableau I. - L'analogie formelle entre le critère de Coulomb à l'état critique et les relations contraintes-
dilatance de Rowe

Chemin de Critère de Coulomb Relations de Rowe

sollicitation rupture à volume constant déformation à volume variable
1. Axisymétrique
en compression σ1,ε1,/2σ3ê3 = - tg2(π/4 + ϕCV,/2) σ1ε1,/2σ3ε3 = - tg2(π/4 + ϕn/2)
o, axial
σ2 = σ3 radial
2. Déformation plane σ1ε1/σ3ε3, = - tg2(π/4 + ϕ 3/2 ) σ1,ε1,/σ3,É3 = - tg2( π/4 ϕf2/2)
ε2 = 0
3. Axisymétrique
en extension 2 σ1,ε1,/σ3É3 = - tg2(π/4 + 2 σ1é1/σ3É3 = - tg2(π/4 + ϕ13/2)
σ1 = σ2 radial
σ3 axial

• Le rapport de puissances qui apparaît au premier —A partir de ces définitions, l’unité formelle du
membre des expressions du tableau I exhibe : tableau 1 s’exprime simplement : les puissances reçues
P+ et et retransmises P_ sont en rapport constant,
—au numérateur, la puissance développée dans les qui s’écrit :
directions principales pour lesquelles la vitesse de
déformation est une contraction, c’est-à-dire la puis­ P+/P. = - tg2 (n/4 + ψψ/2)
sance développée dans les directions principales par
lesquelles le matériau reçoit de l’énergie du système
de sollicitation, que l’on propose de nommer « puis­ On sait relier ce qui précède aux caractères physiques
sance reçue », notée P+ ; élémentaires de la microstructure du matériau et des
- au dénominateur, la puissance développée dans les actions mutuelles entre ses éléments, c’est-à-dire des
contacts à frottement sec, dans le cas d’une sollicitation
directions principales par lesquelles le matériau rend axisymétrique (HORNE, 1965; DE JOSSELIN DE
de l’énergie au milieu extérieur (vitesses de déformation JONG, 1976) ou dans le cas d’un matériau bidimen­
en extension), que l’on propose de nommer « puis­ sionnel (CAMBOU. 1985). Cependant, dans le cas
sance retransmise » notée P_. plus général d’une sollicitation tridimensionnelle, on ne
On note de plus que la somme algébrique des sait guère pour l’instant relier ce qui précède à des
puissances P+ et P_, dans les expressions du tableau I, causes physiques élémentaires.
est égale à la puissance des efforts intérieurs :
Aussi l’auteur propose de poser comme principe la
constance du rapport entre P+ et P_.
20 N° 34 REVUE FRANÇAISE DE GÉOTECHNIQUE

1.2. La transposition tridimensionnelle 2. LA TRACE EXPÉRIMENTALE

Tout ce qui précède a trait à des chemins de
sollicitation bien particuliers, bidimensionnels en fait.
Pour étendre le principe précédent aux chemins tri­
dimensionnels, il convient d’abord de trouver une Les résultats qui fondent l’analyse théorique du cha­
définition objective des puissances P+ et P_, qui est pitre précédent ont été acquis par de nombreuses
donnée par les formulations suivantes : expériences ; on en propose ici une illustration, à
partir d’éléments en majorité publiés dans la littérature,
et reproduits sous leur forme originale. Ces résultats
ont tous trait aux conditions drainées.

2.1. En symétrie axiale, sous sollicitation

monotone
—La relation de Rowe pour l’essai triaxial en com­
Ces définitions vérifient bien sûr la condition : pression s’écrit aussi :
σ1/σ3 = (1 - εv/ε1 ). tg2 (π/4 + ψ/2) (4)

Cette expression prévoit que, dans le diagramme

La formulation tridimensionnelle du principe physique contraintes-dilatance, c’est-à-dire le plan de coordon­
proposé est donc : nées :
(x = 1 — εv/ε1 ; y = σ1/σ3), la trace de l'essai est une
droite (y = kx), indépendante de l’état initial. La
figure 1**, qui décrit des essais sur un sable calcaire
concassé pour trois densités initiales différentes, sous
un même confinement σ3 = 0.05 MPa, traduit cette
propriété : les traces des essais dans le diagramme
Ceci s’écrit encore, de façon plus simple : contraintes-dilatance (fig. 1b) sont effectivement ras­
(3) semblées sur une même droite.
- Si la densité initiale est maintenue constante, mais
que la contrainte de confinement varie d'un essai à
• Cette équation, déjà établie par l’auteur par une l’autre, on sait que l’effet de pic s’affaiblit lorsque le
autre voie (1983), écrit l'égalité de la puissance des confinement augmente, à cause de l’affaiblissement de
efforts intérieurs, au premier membre, avec un second la dilatance.
membre toujours positif. Elle prévoit donc que le
matériau absorbe toujours de l'énergie dès qu’il y a Notons ϕpic. l’angle de frottement interne au pic, défini
déformation (e ≠ 0) sous contraintes (σ ≠ 0 ), ce qui est par sinϕpic = [(σ1, - σ3)/(σ1, + σ1,)]pic. La relation (4)
concevable si cette énergie est dissipée au fur et à relie ϕ au taux de dilatance observé au pic, par :
mesure. Bien que le terme « flow rule » ou « loi
d'écoulement », défini par certains mécaniciens des tg2 (π/4 + ϕpic/2) = (1 - ε v /ε1 .tg (π/4 + ψ/2)
sols, s’applique à la nature de l’équation (3), l’auteur
propose de la nommer plutôt équation de dissipation, La figure 2** montre des résultats au pic obtenus par
terme qui correspond mieux à sa nature physique. CHARLES et WATTS (1980) par des essais triaxiaux
drainés 0 230 mm sur un gravier de basalte concassé,
• Par le jeu des valeurs absolues au deuxième pour des contraintes de confinement comprises entre
membre, l’équation de dissipation change de forme 0,027 MPa et 0,7 MPa. La courbe théorique tracée
lorsque changent de signe une ou plusieurs vitesses par l’auteur pour ψ = 45° représente assez bien les
de déformation principale. (Les contraintes principales points expérimentaux, ce qui traduit l’indépendance
effectives sont toujours en compression dans les maté­ approximative de ψ par rapport à la pression de
riaux granulaires, elles ne changent pas de signe). On confinement.
voit donc s’introduire ici, à partir d’une hypothèse
physique sur les flux d’énergie, la notion de zone On trouvera dans Rowe (1969), de nombreuses réfé­
tensorielle développée depuis longtemps par F. DARVE rences expérimentales plus détaillées.
(1983).
Les résultats de BARDEN et KHAYATT (1966) indi­
quent que les relations contraintes-dilatance de Rowe
• Diverses conséquences de l’équation (3) ont déjà restent valides au cours de sollicitations cycliques en
été détaillées en 1983, elle permet bien sûr de symétrie axiale.
retrouver les relations de Rowe ainsi que le critère de
Coulomb à volume constant en symétrie axiale ou
en déformation plane : une écriture au moyen d’in­
variants tensoriels en a été donnée à cette occasion. ** reproduites avec l’autorisation de la revue Géotechnique.
DILATANCE. DISSIPATION D'ÉNERGIE ET CRITÈRE DE RUPTURE TRIDIMENSIONNEL 21

Essai triaxial en compression

Fig. 1 - A Confinement Fig. 2 - A porosité
constant et porosité initiale constante
initiale variable. et confinement variable.
Résultats au pic.

2.2. En déformation plane BARDEN, ISMAIL et TONG (1969) détaillent des

résultats d’essais sur un sable roulé réalisés à l’Uni­
HUGUES. WROTH et WINDLE (1977) rapportent une versité de Manchester avec un appareil de déformation
série d’essais réalisés par STROUD à l’Université de plane sans rotation d’axes, qui indique la validité des
Cambridge, avec des densités et confinements variés, relations de dilatance sur une large plage de pression
sur un appareil de déformation plane avec rotation de confinement (de 0,14 à 6 MPa).
d’axes (cisaillement simple).
Les courbes d’essai sont reproduites figures 3 ’et 4**,
où l’on note que la distorsion cumulée a été 2.3 En chemin tridimensionnel cyclique
poussée jusqu’à 50 % dans certains essais. —Les figures 5 et 6 reproduisent des essais effectués
sur un sable dense à la presse tridimensionnelle de
Les figures 3b et 4b reproduisent la trace de ces l’Institut de Mécanique de Grenoble et communiqués
essais dans le plan de coordonnées : à l’auteur par J. LANIER en 1984.
I x = (σ1 - σ3)/(σ1 + σ3) = t’/s’ Pour tester la validité de la relation de Rowe sur
y = —(ε1 + ε3)/(ε1 - ε3) = v/ÿ chemin tridimensionnel, LANIER utilise le plan dont
les coordonnées sont :
• en abscisses, l’énergie reçue depuis le début de
On peut montrer sans peine à partir de l’équation (3) l’essai :
que la trace des essais de déformation plane dans ce
plan est théoriquement une hyperbole d’équation :
y = (x —sinj ψ)/(x sin ψ —1)
• en ordonnées, le travail des efforts intérieurs ou
La coïncidence avec la courbe théorique tracée par énergie totale de déformation, cumulé depuis le début
l’auteur pour ψ = 36° est nette sur les figures 3b et de l’essai :
4b.
Un écart entre la prévision théorique et la tendance
expérimentale apparaît sur la figure 4b, au tout début
des essais (zone OA approximativement). Cette dis­ L’équation de dissipation (3) prévoit que :
position traduit vraisemblablement l’effet d’un petit
domaine élastique, il en est de même au tout début [P+(t) + P.(t)]/P.(t)
des essais triaxiaux de la figure 1. = 2 sin Tp/(1 + sin ψ) = constante.
22 N° 34 REVUE FRANÇAISE DE GÉOTECHNIQUE

Essai en déformation plane

I7
(a) Courbes d'essai

" (b) Diagrammes de dilatance\

Fig. 3 - Matériau dense. Fig. 4 - Matériau lâche.

La trace de l’essai dans le plan [X(T), Y(T)], que l'on voir la restitution rapide d’une partie de l’énergie
propose de nommer ici diagramme d'énergie, est donc stockée élastiquement dans l’éprouvette, le retour vers
théoriquement une droite, ce qui est étonnamment la compression (repère (2)) est pratiquement indéce­
bien vérifié sur les figures 5b et 6b. lable.
L'usage des coordonnées intégrales [X(T), Y(T)] a —La figure 6 décrit un chemin tridimensionnel dans
bien entendu pour effet de lisser les fluctuations lequel sont maintenus constants :
instantanées, mais par ailleurs une éventuelle dérive
systématique par rapport à la linéarité théorique prévue
I [σx + σy + σz = 1,5 MPa
par (3) n’en serait que mieux soulignée. \ (σz - σy)/(σx - σy) = 0,5

— La figure 5 décrit un chemin axisymétrique en La droite donnée par la trace de l'essai dans le
contrainte qui consiste en une sollicitation équivalente diagramme d'énergie de déformation, figure 6b. pré­
à un essai triaxial en compression suivi d’un décharge sente une pente d’environ 0,68. qui correspond à un
—avec passage aux conditions du triaxial extension — angle d’écoulement ψ = 31°.
puis d’une recharge —avec retour aux conditions du
triaxial compression. —La figure 7 représente les chemins de ces essais
dans l’espace des contraintes principales ordonnées
La droite donnée par la trace de l’essai dans le σ1> σ2 > σ3
diagramme d'énergie, figure 5b, présente une pente On constate que ces chemins constituent un échantil­
valant environ 0,66, ce qui correspond à un angle lonnage réduit mais régulier des directions possibles
d’écoulement = 29,5°. dans le secteur σ1 > σ2 > σ3.
On note enfin sur la figure 5b que si le passage vers La variation de 1,5° trouvée dans l’angle d’écoulement
l’extension (repère (1) sur la figure 5) est marqué par sur ces chemins n’apparaît pas à l’auteur comme
un léger décrochement en y, dans lequel on peut très significative.
DILATANCE, DISSIPATION D'ÉNERGIE ET CRITÈRE DE RUPTURE TRIDIMENSIONNEL 23

(a) Courbes d'essai

(b) Diagramme d' énergie

Fig. 5. - Essai tridimensionnel cyclique (J. Lanier, 1984) σx = σy = 0.5 MPa.

(a) Courbes d' essai

(b) Diagramme d'énergie

Fig. 6. - Essai tridimensionnel cyclique (J. Lanier, 1984) Pm = 0.5 MPa ; b = 0.5.
24 N° 34 REVUE FRANÇAISE DE GÉOTECHNIQUE

(a) dans l'espace fixe des directions principales lié (b) dans l'espace des contraintes
au laboratoire principales ordonnées

Fig. 7. - Représentation des cycles des figures 5 et 6.

3. LE CRITÈRE que F. DARVE (1983) dénomme loi d’écoulement

DE RUPTURE TRIDIMENSIONNEL dérivatoire. Usuellement ce lien est représenté comme
une correspondance scalaire entre les invariants défi­
SOUS GRANDES DÉFORMATIONS nissant les positions angulaires des déviateurs dans le
repère principal.
L’équation de dissipation (3), bien que reposant sur Il peut également être représenté comme une relation
des résultats acquis sur des chemins bidimensionnels, liant les paramètres b = (σ2 — σ3/(σ1 — σ3) et son
permet une description de situations tridimensionnelles, homologue C = (ε2 — ε3)/(σ1 —σ3), que l’on nommera
qui semble confirmée par l’expérience. Elle donne ici « relation déviatoire ».
comme critère de rupture à volume constant celui de
Coulomb, sur ces chemins simples que sont la symétrie —Dans ce qui suit, on admet a priori qu’à l’état
axiale et la déformation plane. On analyse maintenant critique cette relation déviatoire b = f(c) est continue,
ce qui se passe en dehors de ces chemins simples. monotone croissante, et donc inversible ; comme, par
convention 0 < b < 1, on admet par symétrie que
0 < c < 1. Dans un premier temps, on s’attache à
3.1 Le cadre physique et analytique, la forme du critère de rupture, à relation déviatoire
la relation déviatoire donnée ; dans un deuxième temps on relie cette
notion de relation déviatoire aux chemins de moindre
—Le critère de rupture sous grandes déformations est énergie.
la trace, dans l'espace des contraintes, du régime
d’écoulement limite qui peut être observé après de
larges déformations monotones, lorsque le volume 3.2. Détermination du critère de rupture
spécifique devient stationnaire. Ce régime limite, connu pour une relation déviatoire donnée
depuis longtemps (TAYLOR, 1948) est l’origine de la
notion d’état critique (SCHOFIELD et WROTH, 1968). — Le critère de rupture est donné par la résolution
du système formé par:
—La condition de volume constant εv = 0 = ε1
+ ε2 + ε3, et l’équation de dissipation (3), ne suffisent —l’équation de dissipation
pas pour éliminer les vitesses de déformation et en —la condition de volume constant
tirer la relation sur les seules contraintes qui définit —la relation déviatoire
ce critère de rupture, il manque une relation scalaire
entre o et ê.
frit
—Il existe un lien, observé par l’expérience, entre les
déviateurs de contrainte et de vitesses de déformation, (σ2 - σ3)/(σ1 - σ3) = f[(ε2 - ε3)/(ε1 - e3)] (5)
DILATANCE. DISSIPATION D'ÉNERGIE ET CRITÈRE DE RUPTURE TRIDIMENSIONNEL 25

—La présence de valeurs absolues dans l’équation de • Figure 10, la relation déviatoire, définie par
dissipation incite à examiner séparément tous les cas [c = 1/2), sinon b = 0 et c < 1/2, ou bien
possibles, qui heureusement se réduisent à 2 b = 1 et c > 1/2] est un cas limite de fonction f
(Cf Annexe 1), correspondant à deux modes de défor­ envisagée en 3.1.
mations tridimensionnels fondamentaux, séparés par
un mode frontière en déformation plane : Le critère obtenu est précisément la pyramide de
Coulomb. Par rapport aux exemples précédents, le
• le mode de déformation nommé model qui domaine de la déformation plane s’étend au détriment
recouvre les régimes de déformation intermédiaires des modes 1 et II, réduits à des droites.
entre l’essai triaxial en compression et la déformation Le cas de la figure 10 correspond ici à un minimum
plane, pour lesquels, à Volume constant. 0 < c < 1/2, de résistance, en termes de contraintes. On verra plus
soit : loin qu’il s’agit aussi d’une enveloppe de chemins de
ε1, > 0 moindre énergie.
ε2 < 0
ε3< 0 —A volume constant, la notion de relation déviatoire
est équivalente à celle, utilisée en élastoplasticité, de
« potentiel pour les vitesses de déformation » (c’est-à-
dire une fonction des contraintes F(σ1, σ2, σ3), telles
que les vitesses de déformation soient données par e,
= X . ∂F/∂σ1). En effet, la relation déviatoire, ou son
inverse c = f 1 (b) et la condition de volume constant,
peuvent être écrites sous forme d’un système de deux
équations indépendantes, linéaires et homogènes par
rapport aux ε1 suffisantes pour déterminer un tel
potentiel à une constante additive près. Ce potentiel,
MODE I DEFORMATION MODE I qui n’a ici aucune raison d’être « associé », s’obtient
PLANE aisément pour les exemples simples des figures 8 et
10 :
Essai en déformation plane. • avec la colinéarité des déviateurs (fig. 8) on trouve
un potentiel de Von Mises,
• dans le cas du critère de Coulomb (fig. 10), on
• Le mode II qui recouvre les régimes de déformation obtient un potentiel de Tresca.
intermédiaires entre la déformation plane et le triaxial
Enfin, si la relation déviatoire a déjà pu être mesurée
extension, pour lesquels, à volume constant - < c < 1, directement en sollicitation tridimensionnelle —au pic,
il est vrai, cf. DARVE (1983) - la mesure d’un
soit : potentiel des vitesses semble moins évidente.
ε1 > 0
ε2 > 0 3.3. Chemins de moindre énergie
E3 < 0
—On peut se demander si la relation déviatoire est
bien une réalité intrinsèque, c'est-à-dire physiquement
une sorte de résultat statistique lié aux relations entre
—Pour chacun des modes, on peut alors résoudre efforts et mouvement intergranulaires, ou s’il s’agit
analytiquement le système (5), suivant la méthode seulement d’une observation, qui dépend des condi­
exposée en annexe II. Comme on l’a vu au chapitre 1, tions de mesure et de sollicitation.
le critère de rupture obtenu « touche » la pyramide
de Coulomb, en symétrie axiale en contraintes Si la notion de relation déviatoire n’est qu’une obser­
(b = 0 ou b = 1) ou en déformation plane (c = vation qui dépend des conditions de mesure, le
1/2). On montre en annexe II que, pour une relation matériau suit vraisemblablement un chemin de moindre
déviatoire donnée b = f(c), le critère de rupture énergie, qui doit être compatible, à l’état critique,
correspondant est toujours à l'extérieur de la pyramide avec :
de Coulomb. • les conditions imposées par la sollicitation ;
• la relation de volume constant ;
Les figures 8, 9 et 10 présentent trois exemples :
• l'équation de dissipation.
• Figure 8, la relation déviatoire est b = c, les
déviateurs en contraintes et vitesses de déformation Analytiquement, la notion de moindre énergie se
sont colinéaires. Les domaines correspondant aux traduit par un minimum (lié) de la puissance de
modes I et II couvrent des secteurs angulaires égaux, déformation.
dans chacun d’eux le critère est convexe.
—Considérons la famille de sollicitations tridimension­
• Figure 9, la relation déviatoire est b = c” g étant nelles, très utilisée par les expérimentateurs, où les
tel que, en déformation plane, on ait (soit valeurs de o3,
ici: b(l/2) = (1 — tg ψ/2)/2) suivant un résultat de b = (σ2 σ3)/(σ1σ3,), et ε1,
GREEN cité par RAMAMURTHY et TOKHI (1981). sont maintenues constantes.
26 N° 34 REVUE FRANÇAISE DE GÉOTECHNIQUE

(a) Relation déviatoire b. Critère de rupture à l'état critique

Fig. 8 - a. Relation déviatoire.

Fig. 9

Fig. 10
 DILATANCE, DISSIPATION D'ÉNERGIE ET CRITÈRE DE RUPTURE TRIDIMENSIONNEL 27

On peut tirer la puissance de déformation des équa­ —En outre, si l’on impose la symétrie axiale en
tions développées en annexe II, en fonction de contraintes (b = 0), tous les régimes de déformation
σ3, ε1, b et c, puis chercher la relation c(b) qui intermédiaires entre le triaxial compression et la défor­
minimise cette puissance de déformation : mation plane sont également possibles, de manière
• en mode I, on obtient, en notant S = sin ψ : physiquement stable, car la dérivée de P par rapport
à c s’annule pour b = 0, et P est alors égal à sa
P = ∑σ1,ε1, = 2σ3ε1, {c(2b - 1) + (2 - b)} valeur minimale P0
i
• en mode II :
/{c[2b (1 + S) - (1 - S)] - [b(l + S) - 2(1 - S)]}

P prend la même valeur P0 = 2σ3ε1.S/(l — S) en

symétrie axiale en contrainte (b = 0) ou en défor­
mation plane c = 1/2. Cette valeur P0 est bien le Le minimum de cette fonction monotone croissante
minimum de P sur le domaine du mode 1 (0 < c < 1/2) de c a lieu pour la borne inférieure du mode II
et b > 0) car la dérivée de P par rapport à c vaut : c = 1/2 ; et le maximum a lieu pour la borne
∂P/∂c = - 12σ3ε# tb.S2/D2 supérieure c = 1.
où D est le dénominateur dans l’expression précédente Comme en mode I, le chemin de moindre énergie
de P. est encore le régime de déformation plane, indépen­
damment de la valeur de b. A la différence du
Pour b # 0, on voit que ∂P/∂c < 0. mode I, la symétrie axiale en déformation (c = 1) ne
peut être un minimum, il s’agit même d’un maximum,
La fonction P(c) est donc monotone décroissante et physiquement instable car la dérivée de P par rapport
son minimum, sur le domaine du mode 1 (0 < c < 1/2) à c ne s’y annule pas.
est atteint pour la borne supérieure c = 16/2, où P
vaut précisément P0 = 2σ3e1 . S/(l — S). —Dans une sollicitation monotone à b constant, le
matériau évoluerait donc vers un régime limite de
Le chemin de moindre énergie en mode I, pour cette déformation plane à volume constant, quel que soit
famille de sollicitations, est donc le régime de défor­ b, sauf peut être pour b = 0. Le critère de rupture
mation plane c = 1/2 ; indépendamment de la valeur ainsi obtenu, figure 11, est donc le critère de Coulomb,
. donnée à b. (cf. §3.2).

(a) Relation déviatolre (b) Critère de rupture

Fig. 77. - Enveloppe des chemins de moindre énergie, prévue par l'équation de dissipation, à l’état critique,
pour des sollicitations à o3 et b constants.
28 N° 34 REVUE FRANÇAISE DE GÉOTECHNIQUE

Certains expérimentateurs sur presse tridimensionnelle de la pression de confinement, et de la contrainte

ont noté que, dans un essai monotone à b # 0 intermédiaire. Pour l’auteur, ceci suffit à justifier ce
constant, on voit apparaître un mécanisme de défor­ que l’on propose comme un principe, dont on ne
mation plane, localisé, alors que l’emploi de plaques connaît pas pour l’instant la réduction à des causes
rigides aux frontières de l’éprouvette semblait garantir physiques plus élémentaires.
a priori une certaine homogénéité de déformation
(LANIER, 1984-b). Une explication théorique à cette —On peut penser que ce schéma décrit approxima­
tendance à la déformation plane est suggérée par le tivement un mécanisme micro-structural de déformation
développement précédent, sa localisation ou non est correspondant au glissement relatif de grains quasi­
une autre question. (DESRUES, LANIER, STUTZ - rigides en contact frictionnel.
1984).
D’autres mécanismes peuvent devenir prépondérants
L'essai triaxial classique en compression impose la sur certains chemins de sollicitation (cf § 2.2), ils ne
symétrie axiale en contrainte, b = 0 et les expéri­ sont pas représentés ici.
mentateurs savent bien que, malgré les précautions
sévères (embases lubrifiées, homogénéisation soignée Ceci peut se produire lorsque les grains en contact
à la mise en place, etc.) on peut aussi bien voir ne peuvent pas être considérés comme quasi-rigides.
apparaître comme schéma de déformation, sans raison En particulier, pour une sollicitation de type œdomé-
vraiment apparente : trique, qui impose un même signe —au sens large —
—une déformation homogène symétrique ; à toutes les vitesses de déformation principale, l’équa­
—une déformation homogène avec ovalisation de tion de dissipation indique théoriquement un blocage
l’éprouvette ; de la déformation pour ψ ≠ /2 (la seule solution
possible dans (3) est ε = 0).
—la formation d’un genou dans l’éprouvette (que l’on
attribue ici à une zone de déformation plane) qui finit Cet état de fait donne quelques indications sur le
par se localiser en un plan de cisaillement. domaine de validité de l’équation de dissipation : elle
Le développement théorique précédent souligne cette schématiserait approximativement ce qui se passe sur
espèce d’indétermination physique quand b = 0 des chemins de sollicitation drainés, essentiellement
(triaxial compression). déviatoires, à dominante monotone, où la pression
moyenne ne subit que des variations modérées et
Enfin, il est notoire que si les essais triaxiaux en lentes. Sur de tels chemins, l’hypothèse de coaxialité
compression peuvent être stables, ce qui permet de entre les tenseurs σ et ε semble raisonnable.
les poursuivre jusqu’à de larges déformations, les essais
triaxiaux en extension sont toujours instables. —Ce schéma permet une analyse de la notion de
critère de rupture à l’état critique, qui « colle » à bien
L’analyse des chemins de moindre énergie suggère des faits observés, il expliquerait pourquoi le critère
des raisons théoriques à cet état de fait de COULOMB demeure le critère expérimental de
référence (cf. ARTHUR et al. 1977), tout particuliè­
rement en comportement à volume constant de l’étal
critique (VARDOULAlKIS et GOLDSCHEIDER. 1981).
DISCUSSION Le critère de rupture « au pic » n’est pas accessible
par la seule équation de dissipation combinée avec
— ROWE et HORNE, qui décrivaient leurs résultats l’une des notions de relation déviatoire ou de moindre
énergie.
en symétrie axiale en termes de « incrémental input
energy » et de « incrémental output energy » avaient
bien sûr clairement identifié la nature physique du En effet, la notion de « pic » désigne un état extrémal
rapport σ1,ε1,/2σ3è3 (tableau 1). ROWE note de plus en au cours d'une évolution du matériau, dénommée
1969 à propos de la déformation plane : « la relation usuellement « écrouissage ». L’approche des conditions
[contraintes-dilatance]... constitue simplement une « au pic » nécessite donc une schématisation de cette
généralisation de l’expression de COULOMB-RANKINE évolution.
pour inclure des taux de dilatance différents de
l’unité ». Les critères de rupture à l’état critique obtenus par
l’équation de dissipation et la notion de relation
De vives discussions s’étaient alors développées autour déviatoire, coïncident toujours avec celui de COU­
de ces notions en particulier à cause d’un certain LOMB en symétrie axiale et en déformation plane.
manque d’objectivité (cf. DE JOSSELIN DE JONG
(1976)). En dehors de ces chemins simples l’approche en
moindre énergie indique le même critère de COU­
—Avec les définitions données au § 1-2, les termes LOMB.
de « puissance reçue » et de « puissance retransmise »
trouvent une expression objective, que l’on peut Une validation expérimentale directe de l’approche en
étendre aux situations tridimensionnelles. moindre énergie passerait donc par une expérimen­
tation tridimensionnelle à l’état critique, ce qui paraît
Sur les chemins variés, le rapport expérimental de ces difficile compte tenu des effets de localisation de la
puissances paraît constant, indépendant de l’état initial, déformation sous forme de bandes de cisaillement.
DILATANCE. DISSIPATION D'ÉNERGIE ET CRITÈRE DE RUPTURE TRIDIMENSIONNEL 29

Pour l’auteur, ces effets ne remettraient pas en cause ANNEXE II

la notion même d’état critique, mais simplement en
gênent l’observation dans un volume homogène suf­
fisant pour la mesure. A.2.1. Résolution du système (5)
—Il est utile d’exprimer les contraintes en fonction de
σ1, de la pression moyenne p = (σ1, + σ2 + σ3)/3,
ANNEXE 1 et du paramètre b :
i σ2/p = [3 (1 - b) + (2b - 1) 1, /p]/(2 - b)
Cas de signe possibles dans l’équation de σ o3/p = (3 - (1 + b) σ1/p]/(2 - b) (6)
dissipation
De même, pour les vitesses de déformation, en
à volume constant, respectant les conditions sur la fonction de ε, et du paramètre c, compte tenu de la
relation déviatoire : condition de volume constant :
- Les contraintes principales effectives, qui sont des
compressions, sont positives ; / ε2 = - ε1 (1 —2c)/(2 - c)
- Les directions principales étant indicées suivant ε 3 = - ε1, (1 + c)/(2 - c) l '
l’ordre des contraintes principales, les permutations Compte tenu des signes des εi l’équation de dissipation
possibles sur les signes des vitesses de déformation s’écrit :
principales conduisent à distinguer 23 = 8 cas (Tableau • Mode I (1 —sin ψ) σ1,ε1, + (1 + sin ψ)
A-l).
(σ2ε2 + σ3E3) = 0 (8a)
• Mode II (1 - sin xp) (σ1,ε1, + σ2ε2)
Tableau A-1. - Les 8 permutations possibles sur + (1 + sin ψ) σ3ε3 = 0 (8b)
les signes des t
—A partir des systèmes (6), (7) et des équations (8),
1 2 3 4 5 6 7 8 on obtient :
• en mode 1
É, + + + + - - - -
o,/p = 3 (1 + sin xp) (2 - c - b + 2bc)/[6bc sin xp
- - + - + - + + (3 - sin ip) (2 —c - b + 2bc)] (9a)
ε3 - - + + + + - - • en mode II

σ1/P = 3 [ 2 - b - c + 2bc + sin xp.(b + 3c - 2bc)]

/[(3 + sin xp) (2 - b - c + 2bc) - 6 sin xp.
—A volume constant, les vitesses de déformation (1 - b) (1 - c)] (9b)
principales ne peuvent avoir toutes le même signe, ce
qui exclut les permutations 4 et 7. —Etant donné c, la relation déviatoire donne la valeur
de b, d’où o,/p par les équations (9), et ε2/p ainsi
—On a admis que 0 < c < 1, ce qui implique que que ε3/p par le système (6).
ε2 est vitesse de déformation intermédiaire car :

ε2 = cε1, + (1 - c)ε3 ε1 > ε12 > ε13 a)

et ou A.2.2. Position relative par rapport à la
0 < C< 1 ε1, < ε2 < ε3 ß) pyramide de Coulomb
— Pour une relation déviatoire donnée b = f(c), on
obtient un critère de rupture Cf par la méthode
Ceci exclut les permutations 3 et 8 dans le tableau A-l. précédente. En un point courant de Cf, on peut
toujours définir un angle de frottement interne apparent
par :
—Par convention σ1 > σ2 > σ3; par ailleurs les ε,
vérifient la condition de volume constant On vérifie sin (p = (σ1 - σ3)/(σ1 + σ3).
alors sans peine que les permutations correspondant
aux inégalités (3), (5 et 6 dans le tableau A-l) rendent
On sait que ϕ = ψ en déformation plane ou en
négative la puissance des efforts intérieurs symétrie axiale en contraintes.
P= ce qui est exclu par l’équation de dissi­
II Notons Co la pyramide de Coulomb obtenue avec
pation. un frottement interne égal à ψ. Si sin sin ψ xp le
point de Cf considéré est à l’extérieur de la pyramide
Les seuls cas possibles sont donc les permutations 1 de Coulomb Co, si sin ϕ > sin ψ le point considéré
et 2 dans le tableau A-l. sur Cf est à l’intérieur de Co.
30 N° 34 REVUE FRANÇAISE DE GÉOTECHNIQUE

— A partir de l’expression des contraintes donnée en BARDEN L, KHAYATT A.J. (1966). - « Incremental
A-2.1 et de la définition de sin , on trouve : strain rate ratios and strength of sand in the
• en mode 1 (0 < c < 1/2 ) (10a) triaxial test ». Géotechnique 16, n° 4, 338-357.
sin ϕ/sin ψ — 1 = (1 + sin ψ).b.(l — 2c) CAMBOU B. (1985). — « Les micromécanismes de
/{1 + c +(1 - 2c)[l - b(l - sin ψ)]} déformation plastique dans un sol granulaire ».
Revue française de Géotechnique n°31.
• en mode 11 (1/2 < c < 1) (10b) CHARLES J. A., WATTS K.S. (1980). - «The
influence of confining pressure on the shear
sin ϕ/sin ψ - 1 = 2(1 —sin ψ) (1 - b) (2c - 1) strength of compacted rockfill ». Géotechnique 30,
/{3 +(2c - 1) [2b - 1 + 2(1 - b) sin ψ]} n° 4, 143-167.
CORNFORTH D.H. (1964). - «Some experiments on
Compte tenu de la définition des modes et du fait the influence of strain conditions on the strength
que 0 < b < 1, on vérifie que les dénominateurs of sand ». Géotechnique 14, n° 12, 143-167.
et numérateurs des expressions (10) restent positifs,
seuls les numérateurs pouvant s’annuler. DARVE F. (1983). —«Lois de comportement des sols
sous sollicitations tridimensionnelles et cycliques ».
Les cas sin ϕ = sin ne se produisent en mode I, Réunion-débat du 23.03.1983, Annales de
que si b = 0 (symétrie axiale en contraintes) ou si l’I.T.B.T.P., n° 432, février-mars 1985.
c = 1/2 (déformation plane) et en mode 11 que si DE JOSSELIN DE JONG (1976). - « Rowe’s stress
b = 1 (symétrie axiale en contrainte) ou si c = 1/2. dilatancy relation based on friction ». Géotech­
nique 26, n° 3, 527-534 ; discussion by G.
Dans tous les autres cas, les expressions (10) montrent GUDEHUS (1977) 27, n°3, 439.
que sin ϕ > sin ψ, donc : DESRUES J., LANIER J., STUTZ P. - « Localization
of deformation in tests on sand samples ». Engi­
Quelle que soit la relation déviatoire b = f(c), le neering Fracture Mechanics, 1984.
critère de rupture obtenu (Cf) est toujours à l’extérieur
de la pyramide de Coulomb (Co) sauf en déformation FROSSARD E. (1983). - « Une équation d'écoulement
plane, et en symétrie axiale en contraintes, auxquels simple pour les matériaux granulaires ». Géotech­
cas il y a coïncidence. nique 33, n° 1, 21-29.
HORNE M.R. (1965). - « The behaviour of an assem­
bly of rotund, rigid, cohesionless particles ». Part
NOTATIONS 1 and II. Proc. R. Soc. series A-286, 62-97.
b Paramètre de Habib b = (σ2 — σ3)/(σ1 HORNE M.R. (1969). - « The behaviour of an assem­
bly of rotund, rigid, cohesionless particles ». Part
- σ3) III. Proc. R. Soc. series A-310, 21-34.
c Homologue du paramètre b
c = (ε2 - É3)/(ε1, - ε3) HUGUES J.M.O., WROTH C.P., WINDLE D. (1977).
P+, P_, P, Po Puissances développées dans la défor­ — « Pressuremeter tests in sands ». Géotechnique
mation, rapportées à l’unité de volume 27, n° 4, 455-477.
s’ Moyenne des contraintes principales LAMBE T.W., WHITMAN R.V. (1969). - «Soils
extrêmes Mechanics ». New York, Wiley.
t’ Cisaillement maximal
X(T), Y(T) Energies développées dans la défor­ LANIER J. (1984b). - « Développements récents des
mation, rapportées à l’unité de volume essais mécaniques en laboratoire ». Ecole d’hiver
o Contraintes (effectives) de rhéologie des géomatériaux I.M.G., C.N.R.S.
εε Vitesse de déformation RAMAMURTHY T. TOKHI V.K. (1961). - « Relation
| Vitesse de distorsion maximale of triaxial and plane strain strength ». Proc. 10th
v Vitesse de déformation de volume Int. Conf. on Soil Mech. and Found. Eng.
ϕcv Frottement interne à volume constant Stockholm.
à l’état critique ROWE P.W. (1962). - « The stress dilatancy relation
ϕpic Frottement interne au pic for static equilibrium of an assembly of particles
<j) Angle de friction grain sur grain in contact ». Proc. R. Soc. series A-269, 500­
ϕ1, Angle de friction mesuré par la théorie 527.
de Rowe ROWE P.W. (1969). — « The relation between the
Tp Angle d’écoulement
shear strength of sands in triaxial compression,
plane strain and direct shear ». Géotechnique 19,
BIBLIOGRAPHIE n° 1, 75-86.
ARTHUR J.R.F., DUNSTAN T., AL-ANI Q.A.J.L., SCHOFIELD A.N., WROTH C.P. (1968). - « Critical
ASSADI A. (1977). - «Plastic deformation and state soil Mechanics ». Me Graw-Hill.
failure in granular media ». Géotechnique 27, TAYLOR D.W. (1948). - « Fundamentals of soil
n° 1, 53-74. mechanics ». New York, Wiley.
BARDEN L„ ISMAIL H., TONG P. (1969). - «Plane VARDOULAK1S I.. GOLDSCHEIDER M. (1981). -
strain deformation of granular material at low and « Biaxial apparatus for testing shear bands in
high pressures ». Géotechnique 19, n° 4, 441­ soils ». Proc. 10th lnt. Conf. on Soil Mech. and
452. Found. Eng. Stockholm.