Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 1/60

Organisme(s) : EDF/AMA, CS SI
















Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
Document : R7.01.03




Loi de comportement à double critère Drücker
Prager pour la fissuration et la compression
du béton



Résumé

Le modèle présenté dans ce document est une loi de comportement non linéaire pour le béton. Il s'appuie sur la
théorie de la plasticité, il est valable pour les états de contrainte tridimensionnels. Les hypothèses de
modélisation retenues sont les suivantes :

· un domaine de réversibilité des contraintes délimité par deux critères de type Drücker Prager,
· un écrouissage de chacun des critères,
· en compression, un écrouissage positif jusqu'à un pic, puis un écrouissage négatif,
· en traction, un écrouissage négatif exclusivement,
· une dépendance de l'allure des courbes post-pic dans les deux cas (traction/compression) avec la taille
de l'élément fini (l'allure de cette courbe est liée à l'écrouissage négatif et à l'énergie de fissuration),
· des règles d'écoulement plastiques normales (plasticité associée) et une formulation d'écrouissage
isotrope,
· la prise en compte de la dépendance des seuils d'élasticité par rapport à la température,
· la prise en compte de la dépendance du module d'Young par rapport à la température.
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 2/60


Table
des
matières

1 Notations ............................................................................................................................................... 4
2 Introduction............................................................................................................................................ 5
2.1 Caractéristiques principales du modèle .......................................................................................... 5
2.2 Pourquoi deux critères de Drücker Prager...................................................................................... 5

3 Domaine de réversibilité et fonctions seuils.......................................................................................... 6
3.1 Allure du domaine et des seuils de réversibilité.............................................................................. 6
3.2 Expression mathématique du domaine de réversibilité .................................................................. 8
3.3 Critère de rupture. choix des coefficients a,b,c et d........................................................................ 9
3.4 Analyse du domaine de réversibilité retenu.................................................................................. 12
3.5 Ecrouissage .................................................................................................................................. 19
3.5.1 Fonctions d'écrouissage ...................................................................................................... 19
3.5.2 Courbes d'écrouissage et modules post pic ........................................................................ 21
3.5.2.1 Modèle de fissuration répartie ................................................................................. 21
3.5.2.2 Comportement du béton en traction et courbe post-pic linéaire ............................. 24
3.5.2.3 Comportement du béton en traction et courbe post-pic exponentielle.................... 24
3.5.2.4 Comportement du béton en compression et courbe post-pic linéaire..................... 25
3.5.2.5 Comportement du béton en compression et courbe post-pic non linéaire.............. 26
4 Ecoulement plastique.......................................................................................................................... 27
4.1 Forme générale de la règle de normalité ...................................................................................... 27
4.2 Expression de l'écoulement plastique en partie courante ............................................................ 28
4.3 Expression de l'écoulement plastique au sommet d'un cône ....................................................... 29
4.3.1 Démonstration par la théorie générale des matériaux standards........................................ 29
4.3.2 Démonstration par le travail plastique ................................................................................. 32
4.4 Ensemble des équations de comportement (résumé) .................................................................. 33
5 Intégration numérique de la loi de comportement............................................................................... 35
5.1 Le problème global et le problème local : rappels ........................................................................ 35
5.2 Traitement numérique du cas régulier. ......................................................................................... 37
5.3 Existence d'une solution et condition d'applicabilité ..................................................................... 40
5.4 Traitement des cas non réguliers.................................................................................................. 41
5.4.1 Calcul des contraintes et déformations plastiques .............................................................. 41
5.4.2 Acceptabilité......................................................................................................................... 42
5.4.2.1 Acceptabilité a priori et a posteriori ......................................................................... 42
5.4.3 Existence d'une solution régulière et d'une solution singulière. .......................................... 44
5.4.4 Interversion des sommets des cônes de traction et compression....................................... 45
5.4.5 Projection au sommet des deux cônes................................................................................ 46
5.5 Détermination de l'opérateur tangent............................................................................................ 47
5.5.1 Opérateur tangent en vitesse avec un seul critère actif ...................................................... 47
5.5.2 Opérateur tangent en vitesse avec deux critères actifs....................................................... 48
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 3/60

5.5.3 Dérivées successives des critères en traction et en compression...................................... 49
5.5.3.1 Dérives successives des critères par rapport à la contrainte.................................. 49
5.5.3.2 Dérives successives des critères par rapport aux multiplicateurs plastiques......... 49

5.6 Variables internes du modèle ....................................................................................................... 50
5.7 Organigramme général de résolution ........................................................................................... 50

Annexe 1
snap-back avec les valeurs initiales des coefficients c et d .................................... 56
6 Bibliographie ....................................................................................................................................... 60

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 4/60


1 Notations

désigne le tenseur de contrainte, rangé sous forme de vecteur selon la convention :

11


22


33

12
13


23

On note :

I

1 = Trace()

= 1 tr

la contrainte hydrostatique
H
()
3
1

s = - tr()I
le déviateur des contraintes
3

= 1 tr
la déformation volumique
H
()
3

~ = - 1 tr()I
le déviateur des déformations
3

&
~
3


le taux de déformation équivalent
eq =
trace(~.
2
& &~)
1

J = trace( 2
s
le second invariant des contraintes
2
)
2

eq
3
= 3J =
trace( 2
s
la contrainte équivalente
2
)
2
2
trace( 2
s )

=
J =

oct
3 2
3
I
trace

1
( )


oct = H =
=
3
3

f


c
limite initiale de rupture en compression simple
f
limite initiale de rupture en bi compression
cc
f



c
limite d'élasticité en compression
f


t
limite initiale de rupture en traction
f

t
=

rapport entre limite de rupture en traction et en compression
f c
f

cc
=

rapport entre limite de rupture en bi-compression et compression simple
f c
p



déformation plastique en traction
t
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 5/60




t
multiplicateur plastique en traction
p



déformation plastique en compression
c

multiplicateur plastique en compression
c

f ( p

c
courbe d'écrouissage en compression
c )

f ( p

t
courbe d'écrouissage en traction
t )
u



t
déformation plastique ultime en traction
u



c
déformation plastique ultime en compression
f

G
c
énergie de rupture en compression (caractéristique du matériau)
f

G
t
énergie de rupture en traction (caractéristique du matériau)

le maximum de la température au cours de l'historique de chargement

2 Introduction

2.1
Caractéristiques principales du modèle

Le modèle présenté dans ce document est une loi de comportement non linéaire pour le béton. Il
s'appuie sur la théorie de la plasticité, il est valable pour les états de contrainte tridimensionnels. Les
hypothèses de modélisation retenues reprennent en partie les modèles développés par G. Heinfling
[bib2] et J. F. Georgin [bib1] et sont les suivantes :

· il existe un domaine de réversibilité des contraintes délimité par deux critères de type Drücker
Prager,
· chacun des critères s'écrouit, le domaine de rupture correspond au maximum du domaine de
réversibilité,
· en compression, l'écrouissage est positif jusqu'à un pic, puis il devient négatif,
· en traction, l'écrouissage est négatif exclusivement,
· les courbes post-pic dans les deux cas (traction/compression) varient avec la taille de
l'élément fini (l'allure de cette courbe est liée à l'écrouissage négatif et à l'énergie de
fissuration),
· l'écoulement plastique est régi par une règle de normalité (plasticité associée) la formulation
des écrouissages est de type isotrope,
· le module d'élasticité et les seuils de réversibilité varient avec la température.

Remarque :

La terminologie de critère de traction et critère de compression est discutable. Nous
l'utiliserons par habitude, en étant bien conscient qu'un état de contraintes de traction
peut conduire à l'activation du critère dit de compression.


2.2
Pourquoi deux critères de Drücker Prager

Les auteurs des thèses citées en référence [bib1] et [bib2] utilisent un critère de Drücker Prager en
compression et un critère de Rankine en traction. Ils justifient ces choix par des considérations
physiques en montrant que le domaine de réversibilité ainsi obtenu est proche de la réalité
expérimentale. Par contre ils limitent leurs modélisations à des états de contraintes bidimensionnelles.
Nous avons préféré remplacer le critère de traction par une surface également du type Drücker Prager.
Par ce choix, on s'affranchit de certaines difficultés particulièrement dans les formulations
tridimensionnelles. La surface 3D définissant les états de contraintes admissibles vis à vis de la
traction n'est plus une pyramide (Rankine 3D) mais une surface conique dont le sommet est situé sur
l'axe hydrostatique. La trace du critère « dit de traction » sur le plan déviatorique n'est plus un triangle,
mais un cercle. La formulation obtenue est plus simple. La différence entre les deux critères est
minime pour des états de contrainte proches d'états de contrainte plane. Par contre, pour les états de
contrainte fortement confinés, les deux approches (de Rankine et de Drücker Prager) sont différentes,
ce qui est une limite du modèle proposé.
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 6/60


3
Domaine de réversibilité et fonctions seuils

3.1
Allure du domaine et des seuils de réversibilité

Le domaine de réversibilité est le domaine de l'espace des contraintes à l'intérieur duquel les trajets de
contrainte sont réversibles. Dans l'espace des contraintes principales ( , , , il s'agit de deux
1
2
3 )
cônes dont l'axe est la trisectrice d'équation = = . La [Figure 3.1-a] en donne une
1
2
3
représentation graphique.


2
1
3
Figure 3.1-a




25
oct
C
20
15
10
Domaine de
réversibilité
5
B
oct
0
-
-
-
0
10
P
P
t
c

Figure 3.1-b
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 7/60


Dans un plan ( ,
le domaine de réversibilité est déterminé par deux droites comme indiqué sur
oct
oct )
la [Figure 3.1-b].
Pour un état de contrainte plane, le domaine de réversibilité est la coupe du domaine tridimensionnel
par un plan d'équation = cste , comme indiqué sur la [Figure 3.1-c] , le résultat dans un plan
3
( , étant représenté sur la figure [Figure 3.1-d].
1
2 )


2
1
3
Figure 3.1-c


2
1

Figure 3.1-d

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 8/60


3.2
Expression mathématique du domaine de réversibilité

Il est défini par l'inéquation :

f (, A) 0










éq 3.2-1
dans laquelle A représente les forces thermodynamiques associées aux variables internes (nous
notons l'ensemble des variables internes).
Pour le modèle béton que nous présentons ici, l'équation [éq 3.2-1] prend la forme particulière

a
f
A
f
A

éq
3.2-2
comp (
c )
oct + .
,
=
oct - c + c 0
b
c
f
A
f
A
éq
3.2-3
trac (
t )
oct + .
,
=
oct - t+ t 0
d
a
H
f
( A ) .
,
=
oct
f
A

éq
3.2-4
c
- c + c 0
comp
b
c
H
f ( A )
.
,
=
oct
f
A

éq
3.2-5
t
- t+ t 0
trac
d

Les équations [éq 3.2-2] et [éq 3.2-3] correspondent respectivement aux seuils de « compression » et
de « traction ». Les équations [éq 3.2-4] et [éq 3.2-5] limitent le seuil de réversibilité dans le domaine
de la traction isotrope, elles reviennent à exclure l'axe des abscisses sur la [Figure 3.1-b] au delà des
points P ou P . Il est clair qu'une seule de ces deux dernières conditions suffit. Pour le matériau non
t
c
écroui, le choix des coefficients est tel que OP < OP et la condition [éq 3.2-5] entraîne [éq 3.2-4].
t
c
Nous verrons plus tard que l'écrouissage peut inverser l'ordre des points P et P , rendant la condition
t
c
[éq 3.2-4] plus contraignante que [éq 3.2-5].
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 9/60


3.3
Critère de rupture. choix des coefficients a,b,c et d

Quand l'état de contrainte atteint le bord du domaine de réversibilité, des déformations plastiques se
développent et les seuils se déplacent : ils s'écrouissent. Le seuil de compression « s'agrandit » dans
un premier temps, puis diminue, alors que le seuil de traction ne peut que diminuer. Le seuil de rupture
correspond au domaine maximal pouvant être atteint, il est représenté sur la [Figure 3.3-a] dans un
diagramme de contrainte plane :

Seuil initial de réversibilité
en compression
2
1
Seuil de rupture
en compression

Figure 3.3-a

L'écrouissage des seuils se traduit mathématiquement par l'évolution des quantités A et A , les
c
t
seuils de rupture correspondant aux maximum des fonctions f = f - A et f = f - A . Dans les
c
c
c
t
t
t
modèles retenus, ces fonctions sont telles que : Max f = f et Max f = f ;
c
c
t
t

Les coefficients a, b, c, et d sont donc définis à partir de :

· ft' : la résistance en traction uni axiale du béton,
· fc' : la résistance en compression uni axiale du béton,
· fcc' : la résistance en compression bi axiale du béton,
f
f
On définit de plus les coefficients :
t
=
et
cc
=

f
f
c
c
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 10/60


Pour déterminer les coefficients a, b, c et d il faut se donner 4 équations qui expriment en fait que les
critères sont atteints pour des états de contraintes particuliers et judicieusement choisis.
Une première possibilité consiste à écrire que les deux critères se coupent sur les axes compression
simple (points C de la [Figure 3.3-b]).


A
C
D

Figure 3.3-b

En rappelant que :


2
En compression simple :
<
0 ;

oct =
;oct = -

3
3

2
En bi compression <
0 ;

oct = 2
;oct = -

3
3


2
En traction simple >
0 ;

oct =
;oct =

3
3

On obtient alors les relations suivantes :

Numéro
Etat de contrainte
Critère atteint
relation obtenue
de condition
1 Compression
simple
Compression
a + 3b = 2



3
2
Bi compression
Compression
2a + b = 2

3 Traction
simple Traction - c + 3d = 2
4 Compression
simple
Traction c + 3d = 2
Tableau 3.3-a
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 11/60


f
f
Qui donne, en posant :
t
=
et
cc
=

f
f
c
c
-1
2
a =

2
b =

éq
3.3-1
2 -1
3 2 -1
1-
c = 2
d = 2 2 1 éq
3.3-2
1+
3 1+
Mais ce choix est problématique.

En effet, après écrouissage du critère de traction, et pour une limite de traction devenue nulle le
domaine d'admissibilité prend la forme indiquée sur la [Figure 3.3-c], rendant non admissibles des
états de bi compressions.


Figure 3.3-c

De plus, avec ce choix des coefficients, certains trajets de traction compression simple présentaient
des snap-back comme indiqué en annexe.

Nous avons alors préféré remplacer la condition numéro 4 du [Tableau 3.3-a] par une condition
exprimant que, après que la limite de traction soit retombée à zéro, le domaine de réversibilité est celui
représenté sur la [Figure 3.3-d].


Figure 3.3-d
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 12/60


Ceci conduit à remplacer la relation c + 3d = 2 par c = 2
Le choix des coefficients a,b,c et d est finalement :

-1
2
a =

2
b =

éq
3.3-3
2 -1
3 2 -1
2 2
c = 2 d =

éq
3.3-4
3


modèle retenu

Figure 3.3-e


La [Figure 3.3-e] montre la différence entre les deux modèles pour un état de contrainte plane.

3.4
Analyse du domaine de réversibilité retenu

Dans ce chapitre, nous donnons des indications sur l'ordre de grandeur des contraintes admissibles au
sens du critère retenu. Nous nous attachons à donner des indications sur les contraintes de traction,
notamment pour des états de contrainte tridimensionnels.
La [Figure 3.4-a] montre les domaines initiaux (c'est-à-dire avant écrouissage) pour les valeurs
suivantes des paramètres matériaux :

f


f = 40 Mpa
c
limite initiale de rupture en compression simple :
c
f


f = 44 Mpa
cc
limite initiale de rupture en bi compression
cc
f


cc
=

rapport entre limite de rupture en bi-compression et compression simple
= 1.1
f c
f



=

c
limite d'élasticité en compression ;
33
,
0
f


f = 4 Mpa
t
limite initiale de rupture en traction
t
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 13/60



2
1
3
Figure 3.4-a


Les figures [Figure 3.4-b], [Figure 3.4-c] et [Figure 3.4-d] montrent les coupes du domaine
tridimensionnel par des planes = 0 et = -
Mpa

25
.
3
3


Plan = Mpa

0
3
2
1
3

Figure 3.4-b
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 14/60



Plan = - Mpa

25
3
2
1
3

Figure 3.4-c


Plan = Mpa

0
3
Plan = - Mpa

25
3
2
1
3

Figure 3.4-d

La [Figure 3.4-e] montre les domaines de réversibilité dans un plan ( , pour des états de
1
2 )
contrainte constants, domaines paramétrés par la valeur de . Nous représentons les domaines
3
3
pour = -25 Mpa, = 0 Mpa, = 4 Mpa, = 10 Mpa, = 15 Mpa. On y voit que pour un
3
3
3
3
3
confinement de 25 Mpa de compression, les contraintes de traction peuvent atteindre 15 Mpa, et que,
parallèlement, le domaine de réversibilité pour = 15 Mpa n `est pas vide et correspond à des
3
contraintes de compression et de l'ordre de -25 Mpa. On voit également, que, pour une valeur
1
2
donnée de , la valeur maximale de traction
3
Obtenue pour et est atteinte à l'intersection des critères de traction et de compression.
1
2
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 15/60



Domaine de rupture = - Mpa

25
3
Domaine de rupture = Mpa

0
3
Domaine de rupture = Mpa

4
3
Domaine de rupture =
Mpa

10
3
Domaine de rupture
= Mpa

15
3
Domaine de rupture = - Mpa

25
3
Figure 3.4-e

Nous étudions donc le lieu d'intersection des critères de traction et de compression. Nous notons
( 0 eq
, le point d'intersection des deux critères dans le plan (
eq
, (point C de la [Figure 3.4-f]).
H
)
H
0 )


eq

C
eq
0
Domaine de
réversibilité
H
0

P
P
H
t
c

Figure 3.4-f
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 16/60


Le lieu d'intersection des deux critères dans l'espace des contraintes est donné par :


2 eq



sin
1 =

+ + 0
0
H

3

6

2 eq



sin

2 =

- + + 0
0
H

3

6
= 0
3
3
H - 1 -

2

Où est un paramètre.



Figure 3.4-g

La [Figure 3.4-g] montre les projections de ce lieu dans les plans ( , et ( , .
2
3 )
1
2 )
On peut calculer facilement la valeur maximum de la contrainte le long de cette courbe :

f
2
c


éq
3.4-1
max =
+
f t
3

3

Cette équation montre que, quelle que soit la valeur choisie pour la limite de rupture en traction, la
contrainte maximale atteinte en traction est supérieure au tiers de la limite de rupture en compression.
La [Figure 3.4-h] montre les trois contraintes principales en fonction du paramètre .
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 17/60



Figure 3.4-h

On voit que l'on peut atteindre un niveau de traction de 15 Mpa, mais pour un confinement de
= - Mpa
25
et = - Mpa
25
.
2
3

Pour essayer d'éviter cet inconvénient, qui est important, on peut essayer de jouer sur les valeurs de la
résistance en compression et le paramètre .
A titre d'exemple, nous avons choisi le jeu suivant de paramètres :

f = 20 Mpa
c
f = 40 Mpa
cc
= 2
f = 4 Mpa
t

La [Figure 3.4-i] montre les critères avec ce choix de paramètres. La [Figure 3.4-j] montre la valeur des
contraintes principales à l'intersection des deux critères pour ce nouveau choix de paramètres. La
traction maximum obtenue est plus faible (8 Mpa), mais elle est atteinte pour un niveau de confinement
également plus bas (-7Mpa).
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 18/60



Critère de traction initial à = Mpa

0
3
Critère de traction initial à = - Mpa

8
3
Limite du domaine élastique initial
en compression à
= Mpa

0
3
Critère du pic initial de compression à
= Mpa

0
3
Critère du pic initial de compression à
= - Mpa

8
3
Figure 3.4-i



Figure 3.4-j
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 19/60


3.5 Ecrouissage

Comme nous l'avons déjà mentionné au paragraphe [§3.3] , quand l'état de contrainte atteint le bord
du domaine de réversibilité, les déformations plastiques et les variables internes se développent, les
seuils se déplacent : ils s'écrouissent. Pour notre modèle, les variables internes sont au nombre de
deux, elles sont notées p
pour la variable interne « dite de compression » et p
pour celle « dite de
c
t
traction
». Ces variables déterminent l'évolution des seuils de compression et de traction
respectivement, les forces thermodynamiques leurs sont reliées par les relations :

A = f - f
éq
3.5-1
c
c
c ( p
c )
et
A = f - f








éq 3.5-2
t
t
t ( p
t )
f et f représentent les valeurs des résistances en compression et traction
t ( p
t )
c ( p
c )
respectivement.

3.5.1 Fonctions
d'écrouissage

La fonction f est d'abord croissante puis décroissante, la partie décroissante étant soit linéaire
c ( p
c )
[Figure 3.5.1-a], soit quadratique [Figure 3.5.1-b],


f
p

c ( c )
f'c
f c
pc

u
e
c

Figure 3.5.1-a

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 20/60


f
p

c ( c )
f'c
f c
p

c
e
uc

Figure 3.5.1-b

est une donnée du modèle. L'allure de la courbe entre e
et u
(écrouissage négatif) dépend de
c
c
l'élément, et plus précisément de ses dimensions, suivant un critère analogue à celui choisi par
G. Heinfling, [Erreur ! Source du renvoi introuvable.] pour la prise en compte de la localisation des
déformations.
En traction, l'allure de la courbe donnant la valeur de la limite d'élasticité f en fonction de la
t ( p
t )
déformation plastique cumulée p
ne comporte pas de partie "pré-pic", la partie « post-pic » étant soit
t
linéaire [Figure 3.5.1-c], soit exponentielle [Figure 3.5.1-d].



f
p

t ( t )
ft
g ft
pt
ut

Figure 3.5.1-c
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 21/60



f
p

t ( t )
ft
u
p
t
t

Figure 3.5.1-d

3.5.2 Courbes
d'écrouissage et modules post pic

3.5.2.1 Modèle de fissuration répartie

L'introduction d'un comportement adoucissant post-pic dans les relations contraintes-déformations
pose un problème majeur. Sous sollicitation statique, au delà d'un certain niveau de contrainte,
correspondant à l'amorce du comportement adoucissant, les équations régissant l'équilibre de la
structure perdent leur nature elliptique. Ces équations du problème mécanique forment alors un
système d'équations aux dérivées partielles mal posées dont le nombre de solutions est multiple. Ce
problème se traduit par une non-objectivité par rapport au maillage. Il en découle une sensibilité
pathologique de la solution numérique à la finesse et à l'orientation du maillage.
Afin de résoudre ce problème, ou au moins, d'en limiter les conséquences sur la fiabilité de la solution
prédite, il est nécessaire d'utiliser des techniques dites de régularisation. L'objet de ces techniques est
d'enrichir la description mécanique du milieu, pour pouvoir décrire des états non homogènes de
déformation, et pour préserver la nature mathématique du problème. On opère cette régularisation en
introduisant, dans la loi de comportement, une longueur caractéristique ou longueur interne, reliée à la
largeur de la zone de localisation. Plusieurs techniques sont possibles pour améliorer la description
mécanique du milieu adoucissant. Elles constituent des limiteurs de localisation. La mise en oeuvre de
ces techniques nécessite en général, des développements numériques délicats. Une approche
intermédiaire entre l'utilisation des modèles classiques et la mise en oeuvre de ces limiteurs de
localisation consiste à faire dépendre la pente post-pic de la relation contrainte-déformation, de la taille
de l'élément, de manière à dissiper à la rupture une énergie constante. Cette approche constitue un
pas vers une description non locale du milieu continu.
Considérons d'abord une fissure réelle de surface S dont la mesure est A [Figure 3.5.2.1-a]. S est une
surface de discontinuité du champ de déplacement u . On suppose que pour créer cette discontinuité,
il faut dépenser une énergie W dont l'expression est : W = G (x)dS , G étant une propriété du
S
f
f
matériau.
Considérons maintenant que l'on veuille représenter le même phénomène, en représentant non pas
une discontinuité de déplacement mais une déformation plastique uniformément répartie dans un
volume V .
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 22/60


tr
d
ij
p
L'énergie dissipée sera : W = dV
ij dt , où on a noté t le « temps à rupture ».
V

dt
r
0
S
V




Lver

Figure 3.5.2.1-a
En faisant la série d'hypothèses suivantes :

· la fissure est plane,
·
G est constant le long de la fissure et donc W = A G
.
,
f
f
·
V est un cylindre de base S et d'une hauteur L ,
ver
tt
d p
·
ij
ij
g
dt est constant dans V .
f =

dt
0

On aboutit finalement à la relation :
tt
p

d
ij
ij
W = Vg = V

dt

= AG
.
éq
3.5.2.1-1
f
f

dt
0
Ou encore :
tt
p

d
G
f
ij
ij
f
g =

dt

=
éq
3.5.2.1-2

dt
L
0
ver
On voit facilement que : g f = u f
, écriture dans laquelle les quantités ( f
g , f , ,
u )
0
()d
représentent respectivement ( f
u
p
g , f , , en traction et ( f
u
p
g , f , , en compression. La
c
c
c
c )
t
t
t
t )
donnée de f
g détermine donc , ceci en traction comme en compression :
u
u
c

G
f
g =
f



=

c
( p d
c
c )
f
p
c

L
0
ver
u
t

G
f
g =
f



=

t
( p d
t
t )
f
p
t

L
0
ver
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 23/60


La quantité f
g est donc liée à la pente de la courbe post pic dans un diagramme contrainte-variable
d'écrouissage, laquelle est liée à la pente post pic dans un diagramme contrainte déformation.
Supposons par exemple que la relation contrainte déformation soit linéaire en régime post pic.
Appelons E
la pente post pic dans le diagramme (,) et h < 0 la pente correspondante dans
T < 0
EE
hE
de diagramme ( f , ) [Figure 3.5.2.1-b]. On a les relation h
T
=
E =
qui montrent
E - E
T
E + h
T
que on doit avoir :
- h < E , faute de quoi le diagramme (,) présente un snap back.

ft

ft
EE
E
h
T
=
T
g f
E
E - E
t
T
t


Figure 3.5.2.1-b
La condition - h < E est dite condition d'applicabilité, elle se traduira par une inégalité sur f
g et donc
sur L .
ver
Dans le cadre d'une résolution par la méthode des éléments finis, le volume élémentaire représentatif
du milieu fissuré peut être assimilé à un élément du maillage. La longueur caractéristique (notée par la
suite lc) introduite dans la méthode de l'énergie de rupture équivalente correspond à la longueur Lver.
Lors d'un calcul correspondant à une structure quelconque, la détermination de cette longueur
caractéristique est délicate. Elle dépend de la position du plan de fissure, des dimensions et du type
des éléments ...
Une estimation simple pour les cas bidimensionnels peut s'exprimer sous la forme :
l = r A A
c
e
e est l'aire de l'élément considéré, et r, un facteur correcteur, valant 1 pour les
éléments quadratiques, et 2 pour les éléments linéaires.
On peut étendre cette formulation au cas 3D :
3
l = r V V
c
e
e désigne le volume de l'élément.
Concernant l'évolution de l'écrouissage avec la température, nous considérons comme dans [bib2] que
les énergies de rupture et les résistances à rupture dépendent non pas de la température courante T
du point matériel considéré au temps t , mais de la température maximale atteinte en ce point depuis
le début du chargement jusqu'au temps t . Quand nous aurons besoin de montrer la dépendance des
quantités par rapport à la température, nous noterons :

· le maximum de la température depuis le début de chargement,
·
f
la résistance en compression,
c ()
·
f ( désigne la résistance en traction,
t )
·
f , la courbe d'écrouissage en compression,
c (
p
c )
·
f , la courbe d'écrouissage en traction.
t (
p
t )
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 24/60


3.5.2.2 Comportement du béton en traction et courbe post-pic linéaire

Dans cette modélisation, le béton est supposé élastique jusqu'à sa résistance en traction f . La
t
courbe f en traction est représentée sur la [Figure 3.5.1-c] et est entièrement définie par la
t ( p
t )
résistance en traction du matériau, l'énergie de fissuration
f
G , et la longueur caractéristique lc.
t
L'expression mathématique de cette courbe est :
f
,
f
1

éq
3.5.2.2-1
t (
p
pt )= ()


-
t
t

ut ()

L'équivalence de l'énergie dissipée permet d'écrire :

ut
u
f
G
l
f
,
l f
1
d

t (
t
p
) = c t( pt )

= c t ()

t
p

-
0
0
ut ()


t
d'où


l f
f
G

éq
3.5.2.2-2
t (
u
)
.
c
t ().t ()
=
2
et
f

G
.
2
u
éq
3.5.2.2-3
t ()
t ()
=

l . f
c
t ()
La condition d'applicabilité s'écrit :
.
2 E(). f
Gt ()
l

éq
3.5.2.2-4
c
2
ft ()

3.5.2.3 Comportement du béton en traction et courbe post-pic exponentielle

Dans cette modélisation, le béton est supposé élastique jusqu'à sa résistance en traction f . La
t
courbe f en traction est représentée sur la [Figure 3.5.1-d] et est entièrement définie par la
t ( p
t )
résistance en traction du matériau, l'énergie de fissuration
f
G , et la longueur caractéristique lc.
t
L'expression mathématique de cette courbe est :

f
,
f
.exp
a

éq
3.5.2.3-1
t (
p
pt )= ()



-
t
t

ut ()
L'équivalence de l'énergie dissipée permet d'écrire :



f
G
l
f
l f
a
d

t (
p
) = c t(,pt )

= c t ()

t
p


exp -
u
t
0
0

t ()



Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 25/60


D'où :

u
.
.

G ft ( ) l f
c
t ( ) t ( )
=
éq
3.5.2.3-2
a
et
ut ( )
f
Gt ( )
=

éq
3.5.2.3-3
a

lc. ft ( )


p
p
Soit encore : f
,
f
.exp
l . f

t ( t ) =
t (
)

- c t () t

f
Gt ()

2
l . f
c
t

La pente maximale de la courbe est alors h

max ()
( )
= -
f
Gt ()
et la condition d'applicabilité s'écrit :

E(). f
Gt ()
l
éq
3.5.2.3-4
c


2
ft ()

3.5.2.4 Comportement du béton en compression et courbe post-pic linéaire

Dans cette modélisation, le comportement du béton est supposé élastique jusqu'à la limite d'élasticité,
donnée par un coefficient de proportionnalité (noté ) en pourcentage de la résistance au pic f
.
c ()
Pour les bétons standard est de l'ordre de 30%. La courbe f en compression est représentée
c ( p
c )
sur la [Figure 3.5.1-a] et est entièrement définie par la résistance en traction du matériau, l'énergie de
fissuration
f
G , et la longueur caractéristique lc.
c
L'expression mathématique de cette courbe est :



u
p 2
fc(, pc )
c

= f 'c ( ) + (2 -
2 )
c

+ -1
si p
c e
éq 3.5.2.4 - 1





2
e ( )
(
)
( )
e ( )






p
u
f
,
f '
( c - c )
c ( c ) = c ( )
( )
si e pc u

u

c
éq 3.5.2.4 - 2

( e ( )- c ( ))
( )
( )


Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 26/60



fc
La résistance en compression maximum est atteinte lorsque :
2 2

e () = (
- ) ( )
E()
L'équivalence de l'énergie dissipée permet d'écrire :
u
c
f
G = l
f , d
c ( )
c c (
p
c )
p
c
0
d'où
f
2
1
1
G
l . f

éq
3.5.2.4-3
c ()

= c c () + e()+ uc()


6
2

et
f
u
G
.
2
2
1

éq
3.5.2.4-4
c ()
c ()
+
=
-
e()

l . f
3
c
c ()


fc
La pente de la courbe est alors h()
( )
= -

uc()- e()
et la condition d'applicabilité s'écrit :

E(). f
Gc ()
l

éq
3.5.2.4-5
c
6


2
fc (
2
) 11-
4 -
4

3.5.2.5 Comportement du béton en compression et courbe post-pic non linéaire

Dans cette modélisation, le comportement du béton est supposé élastique jusqu'à la limite d'élasticité,
donnée par un coefficient de proportionnalité (noté ) en pourcentage de la résistance au pic f
.
c ()
Pour les bétons standard est de l'ordre de 30%. La courbe f en compression est représentée
c ( p
c )
sur la [Figure 3.5.1-b] et est entièrement définie par la résistance en traction du matériau, l'énergie de
fissuration
f
G , et la longueur caractéristique lc.
c
L'expression mathématique de cette courbe est :




p
p 2
fc(, pc )
c

= f c( ) + (2 -
2 )
c

+ -1
si p
c
e
éq 3.5.2.5 - 1



2
e ( )
(
)

( )


e ( )






2

( pc -e )
fc ( , c ) = f c( )
( )
1-
si e pc uc
éq 3.5.2.5 - 2


u
2


(c ()-e())
( )
( )


Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 27/60



fc
La résistance en compression maximum est atteinte lorsque :
2 2

e () = (
- ) ( )
E()
L'équivalence de l'énergie dissipée permet d'écrire :

u
c
f
G = l
f , d
c ( )
c c (
p
c )
p
c
0
d'où :
f
2
G
l . f

éq
3.5.2.5-3
c ()

= c c () uc()
+ e()


3
3

et :
f
u
3 G

éq
3.5.2.5-4
c ()
c ()

=
- e()
2

l f
2
.
c
c ()

.
2 fc
La pente maximale de la courbe post-pic est alors h

max ()
( )
= - uc()-e()
et la condition d'applicabilité s'écrit :
f
3 E( ).G ( )
l
c
c
1


éq
3.5.2.5-5
2
2
fc (
2
) 4 - -



4 Ecoulement
plastique

Dans ce paragraphe, nous donnons l'expression des vitesses de déformation plastique, en distinguant
le cas dit général où l'état de contrainte est situé sur une zone « régulière » du bord du domaine de
réversibilité et le cas où il est au sommet d'un des cônes.

4.1
Forme générale de la règle de normalité

Dans l'espace (, A) , les inégalités [éq 3.2-2], [éq 3.2-3], [éq 3.2-4], [éq 3.2-5], définissent un
domaine convexe que nous noterons C
la fonction indicatrice de ce convexe :
(
. Nous noterons
, A)
c

(
si

0
, A) C
(, A)

éq
4.1-1
c
=
(,A)
sinon

Quand la frontière du domaine de réversibilité est atteinte, des déformations plastiques irréversibles se
développent, selon la théorie classique de la plasticité.
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 28/60


Pour un matériau standard, [bib4] la loi d'écoulement vérifie le principe du travail plastique maximum,
ce qui se traduit par l'équation : ( p ) éq 4.1-2
c
& ,&

note le sous différentiel de la fonction . Nous rappelons [bib3] que le sous différentiel d'une
c
c
fonction convexe en un point x est l'ensemble des vecteurs z tels que
:
f ( *
x ) f (x)
*
*
+ z, x - x x

On voit alors facilement que [éq 4.1-2] entraîne :


éq
4.1-3
c ( *
*
, A ) (, A)+ ( *
- )+ ( *
A - A)
*
*


A
et

&p
c
&

Compte tenu de la définition de la fonction caractéristique, on voir facilement que [ éq
4.1-3]
est
équivalent à :

p
p
+ A
*
+
*
A
*
et

*
A C(
éq
4.1-4
, A)
& & &
&

En d'autres termes l'écoulement plastique est tel que le couple (, A) réalise le maximum de la
dissipation plastique parmi les forces thermodynamiques admissibles.


4.2
Expression de l'écoulement plastique en partie courante

Quand la fonction f est différentiable au point considéré (,) la règle de normalité s'écrit
simplement
p
f
=
& & éq
4.2-1


f
=
& & éq
4.2-2
A


& et f vérifiant les conditions de Kuhn-Tucker :
0
&
f 0 éq
4.2-3

. f = 0
&

La variable d'écrouissage est liée au multiplicateur plastique par la loi d'écrouissage. En utilisant le
travail plastique, on peut écrire :
p
f


& = & .
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 29/60


f

Si f est une fonction homogène d'ordre 1 par rapport à la variable tensorielle , on a
= f , ce

qui conduit à l'égalité : =
& & et donc finalement aux équations :
f
p

éq
4.2-4
c = p
comp
c
& &


f
p
p


éq
4.2-5
t =
trac
t
& &


4.3
Expression de l'écoulement plastique au sommet d'un cône

Nous donnons deux présentations du même résultat. La première présentation utilise la théorie des
matériaux standards généralisés et les sous différentiels, la seconde part d'une égalité posée a priori
sur le travail plastique.

4.3.1 Démonstration par la théorie générale des matériaux standards

Le domaine C(


est constitué de deux cônes. La fonction
n'est pas différentiable soit à
, A)
c
l'intersection de ces deux cônes, soit au sommet de chacun de ces cônes. Quand le point (, A)
appartient à l'intersection des deux cônes, les équations précédentes restent valables, avec la
précision que les déformations plastique de compression et de traction se développent en même
temps. Ce cas dit « multi critère » est du reste traité dans [bib4]. Nous nous contenterons ici de traiter
le cas où (, A) est au sommet d'un cône, et nous choisirons le cas le plus fréquent du sommet du
cône de traction, sachant que le cas du sommet du cône de compression se traite exactement de la
même façon.
Les critères sont réécrits en utilisant les variables eq
et , plus pratiques dans les développement
H
analytiques.
2
c
f
(, A )
f
A

4.3.1-1
trac
=
eq
t
+ H - t+ t 0
3d
d
H
f (
c
, A
f
A

4.3.1-2
t ) =
H - t+ t 0
trac
d
Nous considérons donc un cas où :
eq = 0


c









4.3.1-3

f
A
H -
t+ t = 0
d

En partant de [éq 4.1-4], nous allons calculer la dissipation plastique comme étant le maximum de
p *
*

A
& +& pour tous les couples * , *A C(
, A)

p
D = Max ( p *
*

A )
4.3.1-4
,
C ,A &
+ &
* A*
(
)
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 30/60


En écrivant alors que ce maximum est fini et atteint quand * = et A* = A , nous trouverons des
conditions sur p
& et &. En fait, le caractère fini suffira.
En utilisant la décomposition des tenseurs en partie isotrope et déviatoire, et la forme particulière des
variables d'écrouissage, on trouve facilement :
p *
*
p *
*
p
p
*
+ A = ~
s + 3 + A
&
&
&

4.3.1-5
H H
t
t
& &
Considérons alors l'ensemble des vecteurs contraintes de trace nulle et dont la contrainte
1
équivalente de Von Mises vaut 1 : = eq =
=
1
{ ,
,
1
(
trace )
}0
=
eqs
I
1 +

H
s1
1
( , A)


2
c
C(

4.3.1-6
,A)

eq +
f
A
H -
t+ t 0
3d
d

c
f
A
H -
t+ t 0
d
En d'autres termes, la « direction » du déviateur des contraintes est quelconque pour un couple
( , A)C( .
, A)
On peut donc écrire :
p

*
eq
p *
*
p
p *
D =

Max ~
Max

4.3.1-7
s
& s1 +3H + A
H
t
t
& &
*


1 1

eq
*
* *

,
,
,s

t
A
H
1
1
2
*
eq

+ c *
- f + *0

3d
d H
t
A
t

c *
*

- f + 0
d H
t
A
t
~
Il est clair que le maximum de
*
ps
& est atteint quand *sest « parallèle » à p
&~ et que l'on a alors :
1
1
p *
2 p
Max ~
s
.
1 =

1
s
1
&
~eq
3 &



[éq 4.3.1-7] peut donc s'écrire :
p
2
*
~ p eq
D =
eq +
Max
&
Max
( *
3
pH
H
)

+
&
*
3

*
eq
*
*
eq
*
*



,
,


,
,
t
A
H
t
A
H

2
*

eq
c *
*
2
*

f
0


+ - +


eq + c * - f +
*0

3d
d H
t
A
t
3d
d H
t
A
t



c

*
*
c *
*

- f +

0

- f +

0
d H
t
A
t
d H
t
A
t




4.3.1-8
+
Max
(p *A
t
t )


&
*

eq *
*
,
, t
A
H

2
*
eq
c *
*



+ - f +

0
3d
d H
t
A
t


c *
*


- f +

0
d H
t
A
t



Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 31/60


Comme ~
peq 0
&
, on a pour le premier terme :
2
*
p
eq
2
~
p
3
~
d * 3d
3c *
Max

A
f
4.3.1-9
eq
= eq -
+
t-

&
&
*
3
3

2 t
2
2 H
eq *
*
,
, t
A
H
2
*

eq + c * - f +
*0
3d
d H
t
A
t

c *
*

- f +

0
d H
t
A
t
[éq 4.3.1-9] reporté dans [éq 4.3.1-8] donne :

p
~ p
D = 2d f +
*
3
-
c ~
2
+
A* -
d ~
2
4.3.1-10
eq
t
&
Max ( ( p
p
H
eq
*
*
&
&
Max
*
*
&
&
H
)
( t ( p
p
t
eq )
, tA
,
H
t
A
H
c *
*
c
- f +

0
* - f +
*0
d H
t
A
t
d H
t
A
t
Posons alors :
p
~ p
m =
3
- 2c ,
p
~ p
n = - 2d et q =
d p

eq f
H
eq
&
&
t
eq
&
&
&
~
2


Avec ces notations, [éq 4.3.1-10] devient :

p
D = q +
Max ( *
*
m + nA
4.3.1-11
t
H
)
*
*
, tA
H
c *
*
- f +

0
d H
t
A
t

Il s'agit d'un problème de type « simplexe ». Le domaine de *
*
, A est représenté sur la
H
t
[Figure 4.3.1-a].


*
At
f t
*
*
Domaine , A
H
t
*
d
H
f t
c

Figure 4.3.1-a
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 32/60


Comme le domaine *
*
, A s'étend vers - à la fois pour *
et *
A , pour que p
D soit fini, il faut
H
t
H
t
que m et n soient positifs. Le maximum de
*
*
m + nA est atteint pour un couple *
*
, A situé sur le
t
H
H
t
bord du domaine de *
*
, A .
H
t
p
*
c
On a alors : D = q + nf

t + Max m - n

H


*


d
H
Pour que
p
D soit fini, il faut que :
c
m = n
d
Reprenant les définitions de m et n , cette relation donne :
p
c p
=
H
t
&3










éq 4.3.1-12
d &
Par ailleurs, les contraintes m 0 et n 0 donnent :
p
~ p

3
2c








éq 4.3.1-13
H
eq
&
&
et
p
~ p
2d








éq 4.3.1-14
t
eq
&
&
ces deux dernières inégalités étant équivalentes du fait de [éq 4.3.1-11].
Les équations [éq 4.3.1-11] et [éq 4.3.1-12] définissent l'écoulement plastique au sommet d'un des
cônes du domaine de réversibilité.

4.3.2 Démonstration par le travail plastique

Le point de départ est de considérer que par rapport aux développements faits en des points réguliers,
ce sont essentiellement les relations [éq 4.2-1] et [éq 4.2-2], dites règles de normalité, qui ne peuvent
plus être écrites. Or la relation [éq 4.2-1] implique l'égalité p
f
, qui peut, elle, être
t
t ( p
t )
p
&
= &
maintenue.
Nous partirons donc de l'équation :
p
f

éq
4.3.2-1
t
t ( p
t )
p
&
= &
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 33/60


Nous utilisons la décomposition en partie isotrope et déviatoire des tenseurs et trouvons :
p
f = ~
s +
3

&
&
éq
4.3.2-2
t
t ( p
t )
p
H H
&
Au sommet du cône de traction, on a les relations [éq 4.3.1-3], qui, portées dans [éq 4.3.2-2] donnent,
en utilisant également [éq 3.5-2] :
p
d
f =
f
t
t ( p
t )
t ( p
t ) H
&
3

éq
4.3.2-3
c
&
c
Et on retrouve donc la relation [éq 4.3.1-12] :
p
=
H
t
&3
.
d &

4.4
Ensemble des équations de comportement (résumé)

On note H la matrice d'élasticité :

+ µ
2


0
0
0



+ µ
2

0
0
0


+ µ
2
0
0
0
H =

0
0
0
µ
2
0
0
0
0
0
0
µ
2
0



0
0
0
0
0
µ
2
Avec :
=
E
E
3 + 2µ
(
et µ =
, et K =

1 + )(1-
2 )
(21+ )
3

Les relations contraintes déformations s'écrivent finalement :

= H(
p
p
- -
éq
4.4-1
c
t )

Pour un point régulier du cône de compression :
a
f
A
f
A

éq
4.4-2
comp (
) 2
,
=
eq
c
+ H - c+ c 0
3b
b
f
p f
0
éq
4.4-3
c
comp =
pc = p comp
c
&
;
& &


Pour un point régulier du cône de traction :
2
c
f
(, A )
f
A

éq
4.4-4
trac
=
eq
t
+ H - t+ t 0
3d
d
f
p
p
p
f
0

éq
4.4-5
t
trac =
t =
trac
t
&
;
& &

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 34/60


Pour un point au sommet du cône de compression :
s = 0











éq 4.4-6
H
f
(
a
, A
f
A
éq
4.4-7
c ) =
H - c+ c = 0
comp
b
p
a p

=
c H
c
&3











éq 4.4-8
b &
p
~ p

3
2a








éq 4.4-9
c
&
&
H
c eq
Pour un point au sommet du cône de traction :

s = 0











éq 4.4-10
H
f (
c
, A
f
A

éq
4.4-11
t ) =
H - t+ t = 0
trac
d
p
c p

=
t H
t
&3











éq 4.4-12
d &
p
~ p

3
2c









éq 4.4-13
t
&
&
H
t eq

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 35/60


5
Intégration numérique de la loi de comportement

5.1
Le problème global et le problème local : rappels

Pour une structure donnée (géométrie et matériau), et pour un chargement donné, les champs de
déplacement, contrainte et variables internes se trouvent en résolvant un ensemble d'équations aux
dérivées partielles non linéaire formées à partir des équations d'équilibre et des lois de comportement.
Le document [bib5] présente l'algorithme dont nous donnons ici un résumé :

u et connu
0
0
Boucle instants t : chargement L = L t
i
( i)
i
u connu ; calcul de la prédiction
0
u

i 1
-
i
Itérations d'équilibre de Newton n
n
u connu ; n
u
u
u
i =
n
i
-
i
i 1
-
Boucle éléments el
Boucle points de gauss g
n
calcul
el
el


=
u


g
g
i
( ni)
loi de comportement :
n
n
n
el
calcul de : el
et el
à partir de el

, el

et


g i
g i
g i 1
-
g i 1
-
g i
n
el


calcul de
g i (selon option)
n
el


g i
Accumulation dans vecteurs et matrices assemblés :
el
n
Accumulation de T
el
Q
T
g
.
dans
n
Q
.
g i
i
n
el
el


Accumulation de calcul de T
g i
el
Q g
Q dans n
K (selon option)
n
g
el
i


g i
Calcul de
n
u
par :
i
n
n
T
n
K . u
= Q
-

. + L
i
i
i
i
itération de recherche linéaire pour déterminer
Actualisation :
n
n
n
u
+1 = u
+ u

i
i
i
SI test convergence OK
fin Newton : pas de temps suivant i = i+1
Sinon
n = n+1
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 36/60


n
Le calcul des contraintes et variables internes el

n
el
à l'itération de Newton n et au temps t à
g i
g i
i
n
el
partir des contraintes et variables internes el

, el

au temps t et de la valeur et

de
g i 1
-
g i 1
-
i 1
-
g i
l'accroissement de déformation dans l'intervalle de temps estimé à l'itération de Newton n consiste à
intégrer les équations [éq 4.4-1], [éq 4.4-2] à [éq 4.4-5] ou [éq 4.4-6] à [éq 4.4-9] ou [éq 4.4-10] à
[éq 4.4-13] selon les cas avec les conditions initiales :
(t

éq
5.1-1
i
= el
1
- )
g i 1
-
(t

éq
5.1-2
i
= el
1
- )
g i 1
-
p
(t
= el
p


éq
5.1-3
i 1
- )
g i 1
-

Avec la condition de chargement en déformation imposée :
(t =
éq
5.1-4
i )
n
el
g i
Le résultat de cette intégration fournira :
n
el

= t
g
( i)
i
n
el
p
p
g = t
i
(
i )
n
el

= t
g
( i)
i

L'objet de ce chapitre est de présenter l'intégration numérique de ces équations. Il s'agit d'un système
d'équations différentielles non linéaires que nous résolvons par une méthode d'Euler implicite. A partir
de maintenant, les quantité au début du pas de temps (connues) seront notés avec un indice -, alors
n
que les inconnues en fin de pas de temps (toutes inconnues sauf
el
= ) seront notées sans indice.
g i
Pour une quantité quelconque a on note
-
a = a - a .
On commence toujours par calculer une solution élastique e
, en supposant qu'il n'y a pas d'évolution
des déformations plastiques et des variables internes. Si au moins un des critères est violé par cette
solution élastique, il y a lieu de calculer des écoulements plastiques. Il faut alors distinguer les cas
réguliers pour lesquels la solution est sur la partie régulière d'un des cônes ou à leur intersection
des cas dits singuliers où la solution est au sommet d'un des deux cônes. La logique permettant
d'examiner et de choisir ces différents cas, et l'algorithme qui en découle sont relativement complexes.
Nous présentons donc en premier le traitement de chacun des cas et expliquons leur enchaînement
postérieurement, dans le chapitre [§5.7].
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 37/60


5.2
Traitement numérique du cas régulier.

On ne présente en détail que le cas où se développent à la fois des déformations plastiques en
traction et en compression et où donc la solution appartient à l'intersection des deux cônes. Notons
toutefois que, même si e
viole à la fois les deux critères, pour autant la solution finale peut très
bien n'appartenir finalement qu'à un des cônes écroui. On est donc amené à rechercher des soldions
dont on postule qu'elles appartiennent à un des deux cônes ou aux deux. Le cas ou elle appartient à
un seul des deux cônes se déduit facilement du cas plus général présenté ici. Les équations que nous
avons à résoudre sont finalement :

µ
se =
-
s +

µ ~
2

éq
5.2-1
µ-
e
K
=
- + 3K
éq
5.2-2
H
-
H
H
K
e
s = s - 2µ( ~ p
~ p
+
éq
5.2-3
c
t )
e
= - 3K

+

éq
5.2-4
H
H
( p
p
H c
H t )
a
f

f

éq
5.2-5
comp (
pc )
2
,
=
eq + H - c(pc )= 0
3b
b
f

p

éq
5.2-6
c =

p comp
c


p
2
c
f
(,
f

éq
5.2-7
trac
)
t
=
eq + H - (tp)
t
= 0
3d
d
f
p
p


éq
5.2-8
t =


trac
t


Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 38/60


En prenant les parties isotropes et déviatoriques des déformations plastiques, les équations [éq 5.2-6]
et [éq 5.2-8] donnent :
p
~


p
s
c
=

éq
5.2-9
c
eq
b
2
p
p a


=









éq 5.2-10
H c
c
b
3
p
~


p
s
t
=








éq 5.2-11
t
eq
2d
p
p c


=









éq 5.2-12
H t
t
d
3
En reportant [éq 5.2-9] et[éq 5.2-11] dans [éq 5.2-3], on trouve :
p
p




e
s
c
t
s = s - µ
2
+


éq
5.2-13
eq

b
2
2d
qui montre que s est parallèle à e
s d'où l'on déduit :
e
s = s éq 5.2-14
eq
eq
e


En reportant [éq 5.2-14] dans [éq 5.2-9] et [éq 5.2-11] on trouve :

p
e
~


p
s
c
=








éq 5.2-15
c
eq
e
b
2
p
e
~


p
s
t
=








éq 5.2-16
t
eq
e
2d
On reporte alors [éq 5.2-15] et [éq 5.2-16] dans [éq 5.2-3] et [éq 5.2-4], et on exprime les critères
[éq 5.2-5] et [éq 5.2-7] avec ces nouveaux résultats. Cela conduit à deux équations ayant comme
inconnues
p

et
p

:
c
t

2
eq
a

e
e
p
2
2
µ Ka
p
Kac
+
f
éq
5.2-17
H -

c
+
-
t
+
- c -
p
c
+
pc =

2
2


(
) 0
3b
b
3b
b
3bd
bd
2
eq
e
c e
p
Kac

p
2
2
µ
Kc
+
f
éq
5.2-18
H -

c
+
-
t
+
- t -
p
t
+
pt =


2
2
(
) 0
3d
d
3bd
bd
3d
d
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 39/60


C'est ce système de deux équations à deux inconnues qu'il faut finalement résoudre. Dans le cas où
les fonctions f et f sont linéaires, c'est à dire dans le cas où on est en régime post-pic linéaire en
c
t
compression comme en traction, il s'agit d'un système linéaire qui sera donc résolu en une itération.
Dans le cas, soit du régime pré-pic en compression (qui est toujours non linéaire), soit de
modélisations à régimes post pic non linéaires, le système [éq 5.2-17] et [éq 5.2-18] est résolu par une
méthode de Newton :
On note f *

,

le critère de compression considéré comme fonction des seules variables
comp (
p
p
c
t )
p

et
p

, de même pour la traction :
c
t

*
2
a
µ
2
Ka2
µ
2
Kac
f

,

=
+ -

+
-

+
- f - +


comp (
p
p
c
t )
eq
e
e
p
p
H
c
2
2


t
c ( p
p
c
c )
b
3
b
b
3
b
bd
3
bd
*
2
c
µ
2
Kac
µ
2
Kc2
f

,

=

+ -

+
-

+
- f - +


trac (
p
p
c
t )
eq
e
e
p
p
H
c
t
2
2
t ( p
p
t
t )
d
3
d
bd
3
bd
d
3
d

La ième itération de Newton pour système [éq 5.2-17] - [éq 5.2-18] est :

i+1
i
p
p




,
-1
f






comp

c
c
( p p
c
t ) i

=
- J

p
p
i





t
t
f

,


trac ( p
p
c
t )
Le jacobien J vaut :
i
i
f


f


comp
comp


p
p







c
t
J

=

i
f

f

trac
trac


p
p







c
t
Avec :

f


µ
2
2

-
comp
Ka
f
+

c ( p
p
c
c )
= -
+
-

p

2
2
p



b
3
b





c
c


fcomp
µ
2
Kac
= -
+



p
3
t
bd
bd

f
2
trac
µ
Kac
=
-
+



p
3
c
bd
bd
f


µ
2
Kc2
f
- +

trac
t ( p
p
t
t )
= -
+
-

p


2
2
p



d
3
d



t
t

Le Jacobien initial du système se déduit des valeurs des dérivées en
p
et
p
, ce qui
t = 0
c = 0
revient à résoudre le système non linéaire en partant de la solution nulle. Les non linéarités sont
introduites par les courbes d'adoucissement. Dans la partie post-pic, lorsqu'elles sont linéaires, la
convergence se fait en une itération. Lorsqu'elles sont non linéaires, la convergence ne nécessite que
quelques itérations. Partir de la solution nulle ne pose donc pas de problème de convergence. Cela
revient à partir de la linéarisation des critères au voisinage de la prédiction élastique.
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 40/60


5.3
Existence d'une solution et condition d'applicabilité

Nous rappelons que la solution du problème [éq 5.2-17] et [éq 5.2-18] doit vérifier les conditions
[éq 4.2-3] et donc entre autres la positivité des accroissements des multiplicateurs plastiques.

p
éq
5.3-1
c 0

p
éq
5.3-2
t 0
Supposons que nous soyons dans un cas de comportement post pic linéaire en traction comme en
compression et appelons respectivement h et h les pentes des parties post pic. Les accroissements
c
t
des multiplicateurs plastiques sont obtenus en résolvant le système linéaire :

µ
2
Ka2
µ
2
Kac
p


+
+
2
a
h
f
2
2
c
+


c

eq
-


e + eH -
p


c (c )
b
3
b
bd
3
bd
b
3
b


éq
5.3-3
2

=
µ
2
Kac
µ
2
Kc



2
c
p
eq
-
e
e
p

+
+
+ h
f

2
2
t
bd
3
bd
d
3
d

t
+ H - t (t )

d
3
d


Puisque les critères de traction et de compression ont été activés en traction comme en compression,
le second membre de ce système est positif. Mais rien n'assure pour autant que la solution de
[éq 5.3-3] sera positive.
Si l'on pose :
µ
2
Ka2
µ
2
Kac

+
+ h
2
2
c
+

HTC = b
3
b
bd
3
bd

éq
5.3-4
2


µ
2
Kac
µ
2
Kc


+
+
+ h

2
2
t
bd
3
bd
d
3
d


On a :
p



2
a
c

eq
e
+ e
f
H -
-
p


c (c )
-1 b
3
b


= HTC

éq
5.3-5

2
c
p
eq
-


f
t

e + eH - t (pt )
d
3
d


Avec :
h
9

t
HTC = 3
+ µ -
+ µ +
+
+
+
+
+
éq
5.3-6
2 (3
) h
K
t
6
(3K 2 ) h h
h K 9h K
16
9
c t
c
t
(h
h
K
c
t
)


4
4

3
9K
9K
µ
3
ht + (3K + µ)

-
+
- µ
3

-
HTC 1 =
1

4

2
2

2
éq
5.3-7
HTC 9K
9K
µ
3
K
(- )

1
3
h
9K 18
12

-
+
- µ
c +
-
+ µ


2
2

2


Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 41/60


On voit que les conditions de positivité [éq 5.3-1] et [éq 5.3-2] conduisent à des relations relativement
compliquées. Si la solution du problème [éq 5.2-17] et [éq 5.2-18] ne vérifie pas les conditions de
positivité [éq 5.3-1] et [éq 5.3-2], cela peut correspondre soit au fait que les coefficients h et h sont
c
t
tels qu'il n'y a pas de solution positive (cela correspondrait à un snap-back dans un diagramme
(,)), soit au fait que la solution n'active finalement qu'un des deux critères.
Examinons le cas plus simple d'un seul critère activé. Supposons pour fixer les notations que le seul
critère activé soit le critère de traction.
On doit avoir :

2
p
µ
2
Kc



2
c
h
f
éq
5.3-8
t
+
+
=
eq
e
2
2
t
+ e
H
-
p




t (
-
t )
d
3
d

d
3
d
On voit réapparaître la condition dite d'applicabilité :
2
2
µ
Kc

+
+ h

éq
5.3-9
t
> 0

3 2
2


d
d

Cette condition est la généralisation de la condition - h < E présentée au paragraphe [§3.5.2.1] dans
un cas particulier de sollicitation uni axiale.
On retiendra donc la stratégie suivante :
2
eq
a
-
2
eq
a
-
Si
e + e
f
et
e + e
f

H -
c ( pc )> 0
H -
c ( pc )> 0
3b
b
3b
b
Activation a priori des deux critères : résolution problème [éq 5.2-17] et [éq 5.2-18]
Si non convergence ou si non vérification conditions de positivité [éq 5.3-1] et [éq 5.3-2]
Recherche solution avec un seul critère activé
Si non convergence ou non vérification condition positivité
Arrêt sur diagnostique de non vérification de condition d'applicabilité du type [éq 5.3-9]

5.4
Traitement des cas non réguliers

Dans ce paragraphe, nous décrivons le traitement discret des équations correspondant à la projection
au sommet du cône de traction, [éq 4.4-10] à [éq 4.4-13], sachant que la projection au sommet du
cône de compression se fait de la même façon. Les équations [éq 4.4-10] et [éq 4.4-12] définissent
l'écoulement plastique dans ce cas, alors que l'équation [éq 4.4-13] est une condition d'acceptabilité
de la projection au sommet du cône.

5.4.1 Calcul des contraintes et déformations plastiques

Les formes discrètes de [éq 4.4-10] à [éq 4.4-13], sont :
s = 0 éq
5.4.1-1
d
=
f - +


éq
5.4.1-2
H
t ( p
p
t
t )
c
p
c
p


3
=


éq
5.4.1-3
t H
t
d
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 42/60


La relation [éq 5.2-4] établie dans le cas régulier est toujours valable, on l'utilise conjointement à
[éq 5.4.1-3] dans [éq 5.4.1-1] et on obtient :

e
c
p
d
- K

=
f - +

éq
5.4.1-4
t
t
H
( p
p
t
t )
d
c
La relation [éq 5.4.1-4] est une équation non linéaire par rapport à la variable
p

que l'on résout par
t
un algorithme de Newton, ce qui permet de calculer
p


par [éq 5.4.1-3] et par [éq 5.4.1-2].
t H
H
Compte tenu de [5.4.1-1], les contraintes sont donc complètement connues. [éq 5.4.1-1] donne
encore :
e
~
s = s - 2
µ
p


éq
5.4.1-5
t
= 0
Cette dernière équation permet de calculer
p
~
et les déformations plastiques sont complètement
t
connues.

5.4.2 Acceptabilité

La forme discrète de la relation [éq 4.4-13] est :
p
p
3


c~
2

éq
5.4.2-1
t H
t eq
[éq 5.4.1-5] donne :
eq
~ p
e

éq
5.4.2-2
t
=
eq
µ
2
En utilisant [éq 5.2-4] et [éq 5.4.2-2], [éq 5.4.2-1] s'écrit :
eq
e
cK
e

- éq
5.4.2-3
H
H
µ
2
5.4.2.1 Acceptabilité a priori et a posteriori
d
-
Pour le critère de traction et pour la partie post pic du critère de compression, =
f +


H
t ( p
p
t
t )
c
est une fonction décroissante de la variable d'écrouissage
p

. On en déduit que
t
e
-
e
- - et donc que :
H
H
H
H
eq cK
eq
e
e
-
e
cK
e

-


-
H
H
H
H
µ
2
µ
2
eq
La condition e
cK
e
-

- est dite condition d'acceptabilité a priori car elle peut être calculée
µ
H
H
2
eq cK
dès la prédiction élastique. La condition e
e

- est dite condition d'acceptabilité a
H
H
µ
2
posteriori.
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 43/60



Direction de projection
eq

( eq
e
e
,
H
)
zone de projection
au sommet du cône
De traction
H
Pt

Figure 5-1

Ces conditions ont une interprétation graphique simple. On peut voir facilement que, dans le cas d'une
solution régulière, on a :
eq
e
eq
-
µ
= 2
e
-
cK
H
H
eq
Cela montre que la solution en contrainte s'obtient en projetant le point ( e
e
,
parallèlement à une
H
)
direction (cK,
µ
2 ) dans un diagramme (
eq
, , comme indiqué sur la [Figure 5.4.2-a]. Les zones
H
)
d'acceptabilité de la projection au sommet sont des cônes dont le sommet et celui du cône de
réversibilité et délimités d'une part par l'axe > OP et par une demi-droite issue du même point et
H
t
de direction (cK,
µ
2 ).
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 44/60


5.4.3 Existence d'une solution régulière et d'une solution singulière.

Dans le cas où la projection au sommet du cône est acceptable a posteriori, il se peut qu'existe aussi
une solution régulière comme on peut le voir sur la [Figure 5.4.2-b].


zone de projection
au sommet a posteriori
eq

zone de projection
au sommet a priori
( eq
e
e
,
H
)
H
2
1
-
P
P
P
t
t
t
Direction de projection pour le
cas régulier

Figure 5.4.2-b
Le sommet du cône de traction avant écrouissage est noté -
P , celui du cône écroui avec un
t
1
accroissement de variable d'écrouissage
p


est noté 1
P , celui du cône écroui avec un
t
t
2
1
accroissement de variable d'écrouissage
p
p


>

est noté 2
P . On voit qu'il existe une solution
t
t
t
1
2
régulière avec
p


et une solution avec projection au sommet du cône pour
p


. Etant donné que
t
t
la solution régulière correspond à un moindre écrouissage, dans le processus d'évolution, elle sera
rencontrée avant la solution avec projection au sommet : c'est donc la solution régulière qu'il faut
retenir.
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 45/60


Voilà pourquoi l'enchaînement des recherches de solution régulières et avec projection au sommet est
le suivant :
eq
i projection au sommet acceptable a priori : e
cK
e
-

-
µ
H
H
2
Calcul de la solution avec projection au sommet :
p

par [éq 5.4-4]
t
Si non
Recherche solution régulière
Si non convergence ou non vérification condition positivité
Calcul de la solution avec projection au sommet :
p

par [éq 5.4-4]
t
eq cK
Vérification d'acceptabilité a posteriori : e
e

-
H
H
µ
2
eq cK
Si pas acceptable : e
e

> -
H
H
µ
2
Arrêt sur diagnostique de non vérification de condition d'applicabilité

5.4.4 Interversion des sommets des cônes de traction et compression

A priori, le sommet du cône de compression correspond à une pression hydrostatique de traction
beaucoup plus grande que celle du sommet du cône de traction. Mais, comme on peut le voir sur la
[Figure 5.4.4-a], on peut trouver une histoire de chargement qui n'active jamais le critère de traction,
qui active et écrouit fortement le critère de compression jusqu'à le rendre strictement inclus dans le
domaine de réversibilité du critère de traction.

Trajet de contrainte
eq

Critère de compression
H
Pt
Critère de compression écroui

Figure 5.4.4-a
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 46/60


Quand les deux critères ont été ainsi inversés, ne devrait plus intervenir a priori que le critère de
compression. Il se peut alors que la solution soit une projection au sommet du cône de compression,
qui se traite exactement comme la projection au sommet du cône de traction.

5.4.5 Projection au sommet des deux cônes

Dans le cas où les deux cônes ont été intervertis et dans le cas où la prédiction élastique viole les deux
critères, il se peut que soient finalement acceptables à la fois la solution de projection au sommet du
cône de compression et au sommet du cône de traction. Dans ces situations, aucun critère ne permet
de sélectionner une solution plutôt que l'autre et on recherchera donc une projection simultanée au
sommet des deux cônes, qui devront donc partager le même sommet, comme indiqué sur la
[Figure 5.4.5-a].

eq
Direction de projection
Critère de
Critère de
le critère de compression
compression
traction
Prédiction élastique
initial
initial
( e
h , eq e) à projeter
Solution avec projection
sur le critère de traction
Critères de
compression et
traction écrouis
h
Solution avec projection
au sommet des deux cônes

Figure 5.4.5-a
La solution avec projection aux deux sommets est obtenue en résolvant le système :
e
a
p
c
p
b
- K
- K

= f - +


éq
5.4.5-1
c
t
c
H
( p
p
c
c )
b
d
a
e
a
p
c
p
d
- K
- K

=
f - +


éq
5.4.5-2
c
t
t
H
( p
p
t
t )
b
d
c
L'état de contrainte est donné par :
s = 0
d
-
b
=
f +

= f - +


éq
5.4.5-3
H
t ( p
p
t
t )
c ( p
p
c
c )
c
a
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 47/60


5.5
Détermination de l'opérateur tangent

Au cours des itérations de l'algorithme de Newton-Raphson, il est nécessaire de calculer la matrice de
raideur tangente. La construction de celle-ci joue un rôle important dans la stabilité, la rapidité et la
précision de la méthode de résolution. Pour conserver ces propriétés, la matrice de raideur tangente
doit être construite à partir d'un opérateur liant l'incrément de contrainte à l'incrément de déformation
de façon précise à la fin du processus de retour sur les surfaces de charge. La matrice de Hooke, ainsi
que les déformations thermiques interviennent comme des constantes lors de la détermination de
l'opérateur tangent cohérent, construit à la fin de l'itération dans l'incrément concerné.
Le calcul de l'opérateur de comportement tangent cohérent prend en compte les déformations
plastiques. Pour des raisons de simplicité, nous avons choisi de calculer l'opérateur de comportement
tangent en vitesse.

5.5.1 Opérateur
tangent
en
vitesse avec un seul critère actif

Dans le cas d'un seul critère actif, par exemple, le critère en compression, le calcul de l'opérateur de
comportement tangent en vitesse est le suivant :
On utilise donc les équations en vitesse, en charge élastoplastique :

f
p
comp
- H - c
= 0


&

éq
5.5.1-1
& &



T
f
f
comp

+ comp p

éq
5.5.1-2
p
c = 0

&


&
c
L'opérateur tangent en vitesse est défini par :


D
& = &
éq
5.5.1-3
En identifiant [éq 5.5.1-3] avec [éq 5.5.1-1] et [éq 5.5.1-2], on trouve classiquement :
1

T
f
f
comp

D = H - H
comp H éq
5.5.1-4





avec :
T
f

f

f

comp
comp
comp
=
H
-

éq
5.5.1-5
p


c
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 48/60


5.5.2 Opérateur
tangent
en
vitesse avec deux critères actifs

Dans le cas où les deux critères sont activés, le critère en compression et le critère en traction, le
calcul de l'opérateur de comportement tangent en vitesse est le suivant :
On part de :

f
f
p
comp
p

- H - c
-
trac
t
= 0


&

éq
5.5.2-1
& &


&


T
f
f
comp

+ comp p

éq
5.5.2-2
p
c = 0

&


&
c

T
f
f
trac

+ trac p

éq
5.5.2-3
p
t = 0

&


&
t
On aboutit à :
f

f
f
f
f
comp

T
T

comp

f
trac


T
T
D = H - H


trac
comp
trac

éq
5.5.2-4
cc
+ ct
+
tc
+





tt














avec :

T
f
f
f
trac
trac


H
- trac

p


t


cc =
T
f
f
f
f
f
f
f
f
f
f
comp
comp


comp

T
trac
trac

T
trac
comp



trac

T



H
-


H
-

-
H


trac
H comp

p




p
c





t











éq 5.5.2-5


T
f

comp
ftrac
-
H






ct =
T
f
f
f
f
f
f
f
f
f
f
comp
comp


comp

T
trac
trac

T
trac
comp



trac

T






H
-
H
-
-
trac
comp

p





p



H




H




c



t






éq 5.5.2-6


T
f
f

- trac H comp






tc =
T
f
f
f
f
f
f
f
f
f
f
comp
comp


comp

T
trac
trac

T
trac
comp



trac

T



H
-


H
-

-
H


trac
H comp

p




p
c





t











éq 5.5.2-7
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 49/60



T
f
f
f
comp
comp


comp

H
-


p

c



tt =
T
f
f
f
f
f
f
f
f
f
f
comp
comp


comp

T
trac
trac

T
trac
comp



trac

T



H
-


H
-

-
H


trac
H comp

p




p
c





t











éq 5.5.2-8

L'expression semble coûteuse à exprimer en terme de produits de matrice et de calcul. Mais, lorsque
les opérations sont faites dans l'ordre qui convient, il suffit de calculer peu de termes. De plus, ce sont
les mêmes termes qui interviennent à plusieurs reprises. Il faut calculer les dérivées des critères par
rapport à la contrainte, et par rapport aux multiplicateurs plastiques, puis les sommes et les produits
avec les valeurs réelles, pour terminer par la constitution des matrices et leurs sommes.
Enfin, la matrice résultante à l'avantage d'être symétrique, ce qui convient à la résolution standard
avec le Code_Aster.

5.5.3 Dérivées successives des critères en traction et en compression

5.5.3.1 Dérives successives des critères par rapport à la contrainte

Les dérivées des composantes isotropes et déviatoriques de la contraintes par rapport au tenseur de
contraintes s'expriment de la façon suivante :
1

1
1
En définissant le vecteur

0 =
0
0


0
les dérivées des critères par rapport au tenseur de contraintes s'expriment de la façon suivante :
f
comp
s
a
=
+

0


b eq
2
b
3
f

s
c
trac =
+

0


2d eq

d
3

5.5.3.2 Dérives successives des critères par rapport aux multiplicateurs plastiques

Dérivée du critère de compression dans le cas d'une courbe post-pic linéaire :
p
f

4

1

comp
c
p
= - . f .
-
si

p
c
2
c
e





3


c
e
e

f


1

comp = f .
si
p
c



-
c
( ut
e )
p
u
e
c
c

Dérivée du critère de compression dans le cas d'une courbe post-pic non linéaire :
p
f

4

1

comp
c
p
= - . f .
-
si

p
c
2
c
e





3


c
e
e

p
f



comp
-
= .
2 f

.
si
p
c
2
c
( c
e
u
-
c
e )
p
u
e
c
c

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 50/60


Dérivée du critère de traction dans le cas d'une courbe post-pic linéaire :
f

f
trac
t
p
u
=
si
p
u
t
t


t
t
Dérivée du critère de traction dans le cas d'une courbe post-pic exponentielle :
a. p
-
t
f
a
trac

u

t
= f
e

p
t u
t
t

5.6
Variables internes du modèle

Nous montons ici les variables internes stockées en chaque point de Gauss dans l'implémentation du
modèle

Numéro de variable interne
Sens physique
1
p
: déformation plastique cumulée en compression
c
2
p
: déformation plastique cumulée en traction
t
3
: température maximum atteinte au point de gauss
4 Indicateur
de
plasticité

5.7
Organigramme général de résolution

L'organigramme comprend les différentes étapes de la résolution, avec le traitement des projections
aux sommets des cônes de compression et de traction de la façon suivante :
en début d'algorithme,
on effectue une projection au sommet du cône de traction :
· lorsque la prédiction élastique vérifie la condition de projection a priori en traction,
· lorsque la prédiction élastique vérifie la condition de projection a priori en compression et
que les sommets des cônes de traction et de compression se sont intervertis sur l'axe
hydrostatique,
on effectue une projection simultanée aux sommets des cônes de traction et de compression :
· lorsque la prédiction élastique vérifie la condition de projection a priori en compression et
que les sommets des cônes de traction et de compression se sont intervertis sur l'axe
hydrostatique, et que la projection au sommet du cône de traction n'a pas donné de solution
valide,
on effectue une projection au sommet du cône de compression :
· lorsque la prédiction élastique vérifie la condition de projection a priori en compression et
que les sommets des cônes de traction et de compression se sont intervertis sur l'axe
hydrostatique, et que la projection au sommet du cône de traction n'a pas donné de solution
valide, et que la projection simultanée aux sommets des deux cônes n'a pas donné de
solution valide,
en milieu d'algorithme,
on effectue une, deux ou trois résolutions standards avec projection sur le critère de compression ou
sur le critère de traction ou sur les deux critères à la fois,
et en fin d'algorithme,
lorsque que les résolutions standards avec l'activation d'un critère (traction ou compression) ou des
deux critères à la fois n'ont pas donné de solution,
on effectue une projection au sommet du cône de traction :
· lorsque la prédiction élastique vérifie la condition de projection a posteriori en traction ,
· lorsque la prédiction élastique vérifie la condition de projection a posteriori en compression
(et que les sommets des cônes de traction et de compression se sont intervertis sur l'axe
hydrostatique,
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 51/60


on effectue une projection simultanée aux sommets des cônes de traction et de compression :
· lorsque la prédiction élastique vérifie la condition de projection a posteriori en compression
et que les sommets des cônes de traction et de compression se sont intervertis sur l'axe
hydrostatique, et que la projection au sommet du cône de traction n'a pas donné de solution
valide,
on effectue une projection au sommet du cône de compression :
· lorsque la prédiction élastique vérifie la condition de projection a posteriori en compression
et que les sommets des cônes de traction et de compression se sont intervertis sur l'axe
hydrostatique, et que la projection au sommet du cône de traction n'a pas donné de solution
valide, et que la projection simultanée aux sommets des deux cônes n'a pas donné de
solution valide.

A l'issue de chaque résolution ayant convergé, on effectue les vérifications de conformité de la solution
suivante :

· validité de la solution par rapport au second critère, lorsque la résolution a été faite avec un
seul critère. Dans tous les cas, il suffit de vérifier que les deux critères calculés avec la
contrainte finale, sont négatifs ou nuls,
· conformité de la solution : on calcule en cours de résolution la contrainte équivalente finale. Il
arrive parfois que la solution se situe au delà du sommet du cône qui s'écrouit, ce qui conduit à
une contrainte équivalente « négative ». Numériquement, cela se traduit par un critère final
strictement positif. Pour vérifier la conformité de la solution, il suffit de vérifier que les deux
critères calculés avec la contrainte finale, sont négatifs ou nuls,
· validité des projections aux sommets des cônes. Il faut vérifier qu'après résolution, lorsqu'on
connaît l'écrouissage du critère, la pente de la droite reliant la prédiction élastique à la
projection est inférieure à la pente de la direction de projection. (condition de projection a
posteriori au sommet des cônes). Dans le cas contraire, cela signifie qu'il existe une solution
avec résolution standard.

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 52/60






Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 53/60



















Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 54/60



















































Remarque :

Dans le cas des projections aux sommets des cônes, on commence systématiquement par la
projection au sommet du cône de traction. Si cette solution est valide, on conserve celle-ci. Dans
le cas contraire, si le critère de traction est activé, on effectue une résolution avec projection aux
sommets des deux cônes. Si la nouvelle solution est valide, on conserve celle-là. Sinon, on
effectue une résolution avec projection au sommet du cône de compression seule.

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 55/60


Remarque :

Dans le cas où les conditions de projection au sommets des cônes sont activées, l'une des trois
solutions doit être valide, mais pour des sauts élastiques particulièrement importants, il se peut
que la résolution n'aboutisse pas. La solution est alors d'effectuer un redécoupage du pas de
temps.



Remarque :

La projection aux sommets des deux cônes simultanément suppose que le critère de traction soit
activé. Il se peut très bien qu'en cas de permutation des sommets, seul le critère de compression
se retrouve activé. Il faut faire alors une projection au sommet du cône de compression seul.

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 56/60


Annexe 1 snap-back avec les valeurs initiales des coefficients c
et d

Nous montrons dans cet annexe le problème de snap-back rencontré dans la simulation d'un essai de traction
suivi compression simples, dans le cas où le choix des coefficients c et d du critère de traction correspond à
une situation où le critère de traction coupe les axes dans un diagramme de contrainte planes, c'est à dire le
choix des coefficients [éq 3.3-1] et [éq 3.3-2] conduisant à un domaine de réversibilité représenté sur la
[Figure 3.3-b] :
Les hypothèses sont les suivantes :

· on ne prend en compte que le critère de traction,
· la courge d'écrouissage est du type : f = f + h ,
t ( p
t )
p
t
t
· on note simplement :
p
= si bien que la courbe d'écrouissage s'écrit :
f
,
t = f
t +

h
t
· le coefficient de poisson nul,
· l'écrouissage est négatif ,
· la condition d'applicabilité est remplie : - E < h < 0 .

On considère un essai uni axial piloté en déformation selon x, comme indiqué sur la Figure 5-



xx
xx
1
xx
P0
0

P1
xx
Temps t
2
xx
0
xx

1

xx
xx
2
xx
P2

Figure 5-a

Dans les autres directions y et z, les conditions imposées sont des conditions de contraintes nulles :

.
yy = zz = 0
On commence par imposer une déformation de traction 1
telle qu'il y ait une plastification en traction, mais
xx
sans que la limite de traction retombe à 0. C'est le point P dans le diagramme contrainte-déformation. On note
1
0
la déformation pour laquelle la limite de rupture en traction apparaît pour la première fois. Au delà du point
xx
P , on impose une décroissance de la déformation, qui entraîne un déchargement élastique du matériau,
1
jusqu'à la valeur 2

de la déformation pour laquelle on a à nouveau une plastification, mais cette fois sous
xx < 0
contrainte de compression. L'objet de cette annexe est essentiellement d'étudier le comportement du modèle
non retenu (celui correspondant aux formules [éq 3.3-1] et [éq 3.3-2]) au delà du point P .
2
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 57/60


A1.1 Calcul des contraintes et des déformations au cours du chargement

Compte tenu des hypothèses rappelées plus haut, les invariants de contrainte valent :

I










éq A1.1-1
1 = xx
xx
J =









éq A1.1-2
2
3
eq
=








éq A1.1-3
xx
L'écoulement plastique se calcule par :
p
c +
=
2
&








éq A1.1-4
xx
&
d
3
A1.1.1 Chemins 0P0 et P0P1

Compte tenu des relations [éq A1.1-1] à [éq A1.1-4], et compte tenu du fait que le long de ce trajet les
contraintes sont positives, on trouve facilement
xx



c + 2 '

h xx +
f
1
3d
t








éq A1.1.1-1
xx = E
1

(c+ 2)2
E

h +
2



9d

A1.1.2 Chemin P1P2

Par définition, P1P2 est un chemin de décharge élastique, le point P2 étant tel que le critère y est à
nouveau atteint
3d
=







éq A1.1.2-1
2
xx
1
c - 2 xx
3d 1
xx
=

+
éq
A1.1.2-2
2
xx
1
c - 2
xx
E

A1.1.3 Au delà du point P2

On s'intéresse maintenant à la pente de la courbe au point P2 dans le repère ( , . Plus
xx
xx )
précisément on s'intéresse à la pente de cette courbe pour une solution dissipative.
En écrivant que le matériau reste plastique au delà du point P2, c'est à dire que l'état de contrainte
reste sur le critère-lequel s'écrouit-, on trouve :
=
d
3







éq A1.1.3-1
xx

&
h&
c - 2
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 58/60


Par ailleurs, en calculant l'accroissement d'écoulement plastique, et en le reportant dans le calcul de
l'accroissement de contrainte, on trouve :
xx


&
&
xx =
+
(c- 2)
&






éq A1.1.3-2
E
3d
On peut alors éliminer & entre [éq A1.1.3-1] et [éq A1.1.3-2] et on obtient :


hE
xx
& =

éq
A1.1.3-3
&
(
= ET
xx
c - 2)2
2
h + E
2
9d
Cette formule donne la pente de la réponse dans le plan ( ,

xx
xx )
Le numérateur est toujours négatif puisque h est négatif .

Le signe de xx
& = E dépend du signe du dénominateur, ainsi deux cas se posent :
2
T
xx 2
&
(c- 2)2

Si h < -E
alors xx
& = E est positif on a une configuration de snap-back.
2
9d
2
T
xx 2
&


(c- 2)2

Si h > -E
alors xx
& = E est négatif et il n'y a pas de snap-back.
2
9d
2
T
xx 2
&

Une nouvelle condition apparaît pour éviter les snap back, condition que nous comparons à celle déjà
évoquée, mais qui concernait en fait les éventuels snap back au point P1.
hE
Pente au point P1 dans le repère ( , :
E =
,
xx
xx )
T1
h + E
hE
Pente au point P2 dans le repère ( , :
E =
.
xx
xx )
2
T
(c- 2)2
h + E
2
9d
(c- 2)2
2
'
f
Or
2
= 3 = 3 t
.
2

'
9


d
fc
On peut par exemple exprimer E en fonction de E en éliminant h :
2
T
1
T

E E
1
E
.
T =
T
2
2
3 E + ET -
1 (
2
1 3
)
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 59/60


Par exemple,
3 2
E
E =
pour E
.
T = -
2
T

1
1
( - 3 2
)
A titre d'exemple, pour E=32000 Mpa, f ' = Mpa
3
et f ' =
Mpa
3
,
38
on trouve E
. Ainsi
T
601
-
t
c
1
E est très faible par rapport à E. comme illustré sur la [Figure A1.1.3-a].
1
T

xx
xx

Figure A1.1.3-a

Ainsi donc, une condition impliquant qu'il n'y ait pas de snap back au point P2 serait trop restrictive et
aboutirait à choisir un matériau pratiquement pas fragile en traction.
C'est pour cette raison que nous avons préféré modifier l'expression des coefficients c et d comme
indiqué au paragraphe [§3.3].

Ceci dit, et même si la solution adoptée, consistant à modifier les coefficients c et d semble
raisonnable, l'exemple traité dans cet annexe montre qu'un problème très simple peut finalement être
un problème de structure : il y a dans cet exemple des conditions d'équilibre, ce sont les conditions :

.
yy = zz = 0
Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Code_Aster ®
Version
5.0

Titre :

Loi de comportement double Drücker Prager

Date
:

19/03/02
Auteur(s) :
C. CHAVANT, B. CIREE Clé
:
R7.01.03-A Page
: 60/60


6 Bibliographie

[1]
J. F. GEORGIN « Contribution à la modélisation du béton sous sollicitation de dynamique
rapide. La prise en compte de l'effet de vitesse par la viscoplasticité » - Thèse (15/01/98).
[2]
G. HEINFLING « Contribution à la modélisation numérique du comportement du béton et des
structures en béton armé sous sollicitations thermo-mécaniques à haute température » -
Thèse (14/01/98)
[3]
R. T. ROCKAFELLAR « Convex analysis » Princeton University Press 1972
[4]
B. HALPHEN et NGUYEN QUOC SON, « Sur les matériaux standards généralisés » Journal
de Mécanique Vol 14 n° 1, 1975
[5]
N. TARDIEU, I. VAUTIER, E. LORENTZ « Algorithme non linéaire quasi statique » document
de référence Aster [R5.03.01].

Manuel de Référence
Fascicule R7.01 : Modélisation pour le Génie Civil
HT-66/02/004/A

Document Outline