background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
1/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Organization (S):
EDF/AMA, CS IF















Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
Document: R7.01.03



Law of behavior to double Drücker criterion
Prager for cracking and compression
concrete


Summary

The model presented in this document is a nonlinear law of behavior for the concrete. It rests on
theory of plasticity, it is valid for the three-dimensional states of stress. Assumptions of
modeling retained are as follows:
·
a field of reversibility of the stresses delimited by two criteria of the type Drücker Prager,
·
a work hardening of each criterion,
·
in compression, a positive work hardening to a peak, then a negative work hardening,
·
in traction, a negative work hardening exclusively,
·
a dependence of the shape of the curves post-peak in both cases (traction/compression) with the size
finite element (the shape of this curve is related on negative work hardening and the energy of cracking),
·
normal plastic rules of flow (associated plasticity) and a formulation of work hardening
isotropic,
·
the taking into account of the dependence of the thresholds of elasticity compared to the temperature,
·
the taking into account of the dependence of the Young modulus compared to the temperature.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
2/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Count
matters
1
Notations ............................................................................................................................................... 4
2
Introduction ............................................................................................................................................ 5
2.1
Outstanding features of the model .......................................................................................... 5
2.2
Why two criteria of Drücker Prager ...................................................................................... 5
3
Field of reversibility and functions thresholds .......................................................................................... 6
3.1
Pace of the field and the thresholds of reversibility .............................................................................. 6
3.2
Mathematical expression of the field of reversibility .................................................................. 8
3.3
Criterion of rupture. choice of the coefficients has, B, C and D ........................................................................ 9
3.4
Analyze field of reversibility retained .................................................................................. 12
3.5
Work hardening .................................................................................................................................. 19
3.5.1
Functions of work hardening ...................................................................................................... 19
3.5.2
Curves of work hardening and modules post peak ........................................................................ 21
3.5.2.1
Model of cracking distributed ................................................................................. 21
3.5.2.2
Behavior of the concrete in traction and linear curve post-peak ............................. 24
3.5.2.3
Behavior of the concrete in traction and exponential curve post-peak .................... 24
3.5.2.4
Behavior of the concrete in compression and linear curve post-peak ..................... 25
3.5.2.5
Behavior of the concrete in compression and nonlinear curve post-peak .............. 26
4
Plastic flow .......................................................................................................................... 27
4.1
General form of the rule of normality ...................................................................................... 27
4.2
Expression of the plastic flow partly current ............................................................ 28
4.3
Expression of the plastic flow at the top of a cone ....................................................... 29
4.3.1
Demonstration by the general theory of standard materials ........................................ 29
4.3.2
Demonstration by plastic work ................................................................................. 32
4.4
Together equations of behavior (summarized) .................................................................. 33
5
Numerical integration of the law of behavior ............................................................................... 35
5.1
The total problem and the local problem: recalls ........................................................................ 35
5.2
Digital processing of the regular case. ......................................................................................... 37
5.3
Existence of a solution and condition of applicability ..................................................................... 40
5.4
Processing of the nonregular cases .................................................................................................. 41
5.4.1
Calculation of the stresses and plastic deformations .............................................................. 41
5.4.2
Acceptability ......................................................................................................................... 42
5.4.2.1
Acceptability a priori and a posteriori ......................................................................... 42
5.4.3
Existence of a regular solution and a singular solution. .......................................... 44
5.4.4
Inversion of the nodes of the cones of traction and compression ....................................... 45
5.4.5
Projection at the top of the two cones ................................................................................ 46
5.5
Determination of the tangent operator ............................................................................................ 47
5.5.1
Tangent operator of speed with only one active criterion ...................................................... 47
5.5.2
Tangent operator of speed with two active criteria ....................................................... 48
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
3/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
5.5.3
Derivative successive of the criteria in traction and compression ...................................... 49
5.5.3.1
Successive drifts of the criteria compared to the stress .................................. 49
5.5.3.2
Successive drifts of the criteria compared to the plastic multipliers ......... 49
5.6
Variables intern model ....................................................................................................... 50
5.7
Top-level flowchart of resolution ........................................................................................... 50
Appendix 1
snap-back with the initial values of the coefficients C and D .................................... 56
6
Bibliography ....................................................................................................................................... 60
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
4/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
1 Notations
indicate the tensor of stress, arranged in the form of vector according to convention:














23
13
12
33
22
11
One notes:
()
=
Trace
I
1
()
=
tr
H
3
1
the hydrostatic stress
()
I
tr
S
-
=
3
1
the diverter of the stresses
()
=
tr
H
3
1
voluminal deformation
()
I
tr
-
=
3
1
~
the diverter of the deformations
()
&
&
&
~
~
.
trace
2
3
~
=
eq
the rate of deformation is equivalent
()
2
2
2
1
S
trace
J
=
the second invariant of the stresses
()
2
2
2
3
3
S
trace
J
eq
=
=

the equivalent stress
()
3
3
2
2
2
S
trace
J
Oct.
=
=
()
3
3
1
=
=
=
trace
I
H
Oct.
C
F
initial limit of rupture in simple compression
DC
F
initial limit of rupture out of Bi compression
C
F
elastic limit in compression
T
F
initial limit of rupture in traction
C
T
F
F
=
relationship between ultimate strength in traction and compression
C
DC
F
F
=
relationship between ultimate strength in Bi-compression and simple compression
p
T
plastic deformation in traction
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
5/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
T
plastic multiplier in traction
p
C
plastic deformation in compression
C
plastic multiplier in compression
()
p
C
C
F
curve of work hardening in compression
()
p
T
T
F
curve of work hardening in traction
U
T
ultimate plastic deformation in traction
U
C
ultimate plastic deformation in compression
F
C
G
energy of rupture in compression (characteristic of material)
F
T
G
energy of rupture in traction (characteristic of material)
the maximum of the temperature during the history of loading
2 Introduction
2.1
Outstanding features of the model
The model presented in this document is a nonlinear law of behavior for the concrete. It
be based on the theory of plasticity, it is valid for the three-dimensional states of stress.
assumptions of modeling selected partly take again the models developed per G. Heinfling
[bib2] and J.F. Georgin [bib1] and are as follows:
·
there is a field of reversibility of the stresses delimited by two criteria of the Drücker type
Prager,
·
each criterion is hammer-hardened, the field of rupture corresponds to the maximum of the field of
reversibility,
·
in compression, work hardening is positive to a peak, then it becomes negative,
·
in traction, work hardening is negative exclusively,
·
the curves post-peak in both cases (traction/compression) vary with the size of
the finite element (the shape of this curve is related on negative work hardening and the energy of
cracking),
·
the plastic flow is governed by a rule of normality (associated plasticity) the formulation
work hardenings is of isotropic type,
·
the modulus of elasticity and the thresholds of reversibility vary with the temperature.
Note:
The terminology of criterion of traction and criterion of compression is debatable. Us
will use by practice, while being quite conscious that a state of tensile stresses
can lead to the activation of the criterion known as of compression.
2.2
Why two criteria of Drücker Prager
The authors of the theses referred to [bib1] and [bib2] use a criterion of Drücker Prager in
compression and a criterion of Rankine in traction. They justify these choices by considerations
physics by showing that the field of reversibility thus obtained is close to reality
experimental. On the other hand they limit their modelings in states of two-dimensional stresses.
We preferred to also replace the criterion of traction by a surface of the type Drücker Prager.
By this choice, one frees oneself from certain difficulties particularly in the formulations
three-dimensional. Surface 3D defining the working states of stresses with respect to
traction is not any more one pyramid (Rankine 3D) but a conical surface whose node is located on
the hydrostatic axis. The trace of the criterion “known as of traction” on the plan deviatoric is not any more one triangle,
but a circle. The formulation obtained is simpler. The difference between the two criteria is
tiny for states of stress close to states of plane stress. On the other hand, for the states of
stress strongly confined, the two approaches (of Rankine and Drücker Prager) are different,
what is a limit of the model suggested.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
6/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
3
Field of reversibility and functions thresholds
3.1
Pace of the field and the thresholds of reversibility
The field of reversibility is the field of the space of the stresses inside of which ways of
stress are reversible. In the space of the main stresses
(
)
3
2
1
,
,
, they are two
cones whose axis is the trisecting one of equation
3
2
1
=
=
. [Figure 3.1-a] one gives some
chart.
1
2
3
Appear 3.1-a


0
5
10
15
20
25
-
-
-
0
10
Oct.
Oct.
T
P
B
C
C
P
Field of
reversibility
Appear 3.1-b
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
7/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
In a plan
(
)
Oct.
Oct.
,
the field of reversibility is determined by two lines as indicated on
[Figure 3.1-b].
For a state of stress planes, the field of reversibility is the cut of the three-dimensional field
by a plan of equation
cste
3
=
, as indicated on [Figure 3.1-c], the result in a plan
(
)
2
1
,
being represented on the figure [Figure 3.1-d].
1
2
3
Appear 3.1-c
1
2
Appear 3.1-d
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
8/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
3.2
Mathematical expression of the field of reversibility
It is defined by the inequation:
()
0
,
With
F
éq 3.2-1
in which
With
represent the thermodynamic forces associated with the variables intern (us
let us note
the whole of the internal variables).
For the model concrete that we present here, the equation [éq 3.2-1] takes the particular form
(
)
0
.
,
+
-
+
=
C
C
Oct.
Oct.
C
comp
With
F
B
has
With
F
éq
3.2-2
(
)
0
.
,
+
-
+
=
T
T
Oct.
Oct.
T
trac
With
F
D
C
With
F
éq
3.2-3
(
)
0
.
,
+
-
=
C
C
Oct.
C
H
With
F
B
has
With
F
comp
éq
3.2-4
(
)
0
.
,
+
-
=
T
T
Oct.
T
H
With
F
D
C
With
F
trac
éq
3.2-5
The equations [éq 3.2-2] and [éq 3.2-3] correspond respectively to the thresholds of “compression” and
of “traction”. The equations [éq 3.2-4] and [éq 3.2-5] limit the threshold of reversibility in the field
isotropic traction, they amount excluding the X-axis on [Figure 3.1-b] beyond
points
T
P
or
C
P
. It is clear that only one of these two last conditions is enough. For material not
hammer-hardened, the choice of the coefficients is such as
C
T
COp
COp
<
and the condition [éq 3.2-5] involves [éq 3.2-4].
We will see later that work hardening can reverse the command of the points
T
P
and
C
P
, returning the condition
[éq 3.2-4] more constraining than [éq 3.2-5].
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
9/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
3.3
Criterion of rupture. choice of the coefficients has, B, C and D
When the state of stress reaches the edge of the field of reversibility, of the plastic deformations
develop and the thresholds move: they are hammer-hardened. The threshold of compression “increases” in
the first time, then decreases, whereas the threshold of traction can only decrease. The threshold of rupture
corresponds to the maximum field being able to be reached, it is represented on [Figure 3.3-a] in one
diagram of plane stress:
Initial threshold of reversibility
in compression
Threshold of rupture
in compression
1
2
Appear 3.3-a
The work hardening of the thresholds results mathematically in the evolution of the quantities
C
With
and
T
With
, them
thresholds of rupture corresponding to the maximum of the functions
C
C
C
With
F
F
-
=
and
T
T
T
With
F
F
-
=
. In
models selected, these functions are such as:
C
C
F
F
=
Max
and
T
T
F
F
=
Max
;
The coefficients has, B, C, and D are thus defined from:
·
F
T
': resistance in axial traction plain of the concrete,
·
F
C
': resistance in axial compression plain of the concrete,
·
F
DC
': resistance in axial compression Bi of the concrete,
One defines moreover coefficients:
C
T
F
F
=
and
C
DC
F
F
=
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
10/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
To determine the coefficients has, B, C and D it is necessary to give itself 4 equations which express in fact that them
criteria are reached for states of stresses particular and judiciously selected.
A first possibility consists in writing that the two criteria are cut on the axes compression
simple (points C of [Figure 3.3-b]).
With
C
D
Appear 3.3-b
By recalling that:
In simple compression:
-
=
=
<
3
2
;
3
;
0
Oct.
Oct.
Out of Bi compression
-
=
=
<
3
2
;
3
2
;
0
Oct.
Oct.
In simple traction
=
=
>
3
2
;
3
;
0
Oct.
Oct.
The following relations then are obtained:
Number
of condition
State of stress
Criterion reached
relation obtained
1 Compression
simple
Compression
2
3
=
+
B
has

2

Bi compression

Compression
2
3
2
=
+
B
has
3 Traction
simple Traction
2
3
=
+
-
D
C
4 Compression
simple
Traction
2
3
=
+
D
C
Table 3.3-a
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
11/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Who gives, while posing:
C
T
F
F
=
and
C
DC
F
F
=
1
2
3
2
1
2
1
2
-
=
-
-
=
B
has
éq
3.3-1
+
=
+
-
=
1
1
3
2
2
1
1
2
D
C
éq
3.3-2
But this choice is problematic.
Indeed, after work hardening of the criterion of traction, and for a limit of traction become null it
field of admissibility takes the form indicated on [Figure 3.3-c], making nonacceptable of
Bi compressions states.
Appear 3.3-c
Moreover, with this choice of the coefficients, certain ways of simple traction compression presented
snap-back as indicated in appendix.
We then preferred to replace the condition number 4 of [Table 3.3-a] by a condition
expressing that, after the limit of traction fell down to zero, the field of reversibility is that
represented on [Figure 3.3-d].
Appear 3.3-d
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
12/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
This resulted in replacing the relation
2
3
=
+
D
C
by
2
=
C
The choice of the coefficients has, B, C and D is finally:
1
2
3
2
1
2
1
2
-
=
-
-
=
B
has
éq
3.3-3
3
2
2
2
=
=
D
C
éq
3.3-4
model selected
Appear 3.3-e

[Figure 3.3-e] shows the difference between the two models for a state of plane stress.
3.4
Analyze field of reversibility retained
In this chapter, we give indications on the order of magnitude of working stresses to
feel criterion selected. We endeavor to give indications on tensile stresses,
in particular for three-dimensional states of stress.
[Figure 3.4-a] shows the initial fields (i.e. before work hardening) for the values
following of the parameters materials:
C
F
initial limit of rupture in simple compression:
C
F
= 40 Mpa
DC
F
initial limit of rupture out of Bi compression
DC
F
= 44 Mpa
C
DC
F
F
=
relationship between ultimate strength in Bi-compression and simple compression
= 1.1
C
F
elastic limit in compression;
33
,
0
=
T
F
initial limit of rupture in traction
T
F
= 4 Mpa
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
13/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
1
2
3
Appear 3.4-a

Figures [Figure 3.4-b], [Figure 3.4-c] and [Figure 3.4-d] the cuts of the field show
three-dimensional by the plane ones
0
3
=
and
Mpa
25
3
-
=
.
1
2
3
Plan
Mpa
0
3
=
Appear 3.4-b
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
14/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
1
2
3
Plan
Mpa
25
3
-
=
Appear 3.4-c
1
2
3
Plan
Mpa
0
3
=
Plan
Mpa
25
3
-
=
Appear 3.4-d
[Figure 3.4-e] shows the fields of reversibility in a plan
(
)
2
1
,
for states of
stress
3
constant, fields parameterized by the value of
3
. We represent the fields
for
25
3
-
=
Mpa,
0
3
=
Mpa,
4
3
=
Mpa,
10
3
=
Mpa,
15
3
=
Mpa. One sees there that for one
containment of 25 Mpa of compression, stresses tensile can reach 15 Mpa, and that,
in parallel, the field of reversibility for
15
3
=
Mpa N `is not empty and corresponds to
compressive stresses
1
and
2
about
25
-
Mpa. It is also seen, that, for a value
data of
3
, the maximum value of traction
Obtained for
1
and
2
is reached with the intersection of the criteria of traction and compression.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
15/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Field of rupture
Field of rupture
Field of rupture
Field of rupture
Field of rupture
Field of rupture
Mpa
25
3
-
=
Mpa
0
3
=
Mpa
4
3
=
Mpa
10
3
=
Mpa
15
3
=
Mpa
25
3
-
=
Appear 3.4-e
We thus study the place of intersection of the criteria of traction and compression. We note
(
)
eq
H
0
0
,
the point of intersection of the two criteria in the plan
(
)
eq
H
,
(not C of [Figure 3.4-f]).
H
eq
C
T
P
C
P
Field of
reversibility
eq
0
0
H
Appear 3.4-f
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
16/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
The place of intersection of the two criteria in the space of the stresses is given by:




-
-
=
+


+
-
=
+


+
=
2
1
0
3
0
0
2
0
0
1
3
6
sin
3
2
6
sin
3
2
H
H
eq
H
eq
Where
is a parameter.
Appear 3.4-g
[Figure 3.4-g] shows projections of this place in the plans
(
)
2
1
,
and
(
)
3
2
,
.
One can easily calculate the maximum value of the stress along this curve:
T
C
F
F
+
=
3
2
3
max
éq
3.4-1
This equation shows that, whatever the value chosen for the ultimate strength in traction,
maximum stress attack in traction is higher than the third of the ultimate strength in compression.
[Figure 3.4-h] shows the three main stresses according to the parameter
.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
17/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Appear 3.4-h
It is seen that one can reach a level of traction of 15 Mpa, but for a containment of
Mpa
25
2
-
=
and
Mpa
25
3
-
=
.
To try to avoid this disadvantage, which is important, one can try to exploit the values of
resistance in compression and the parameter
.
As example, we chose the following play of parameters:
C
F
= 20 Mpa
DC
F
= 40 Mpa
= 2
T
F
= 4 Mpa
[Figure 3.4-i] shows the criteria with this choice of parameters. [Figure 3.4-j] the value shows of
main stresses with the intersection of the two criteria for this new choice of parameters.
maximum traction obtained is weaker (8 Mpa), but it is reached for a level of containment
also low (- 7Mpa).
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
18/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Initial criterion of traction with
Limit of the initial elastic range
in compression with
Criterion of the initial peak of compression with
Initial criterion of traction with
Criterion of the initial peak of compression with
Mpa
0
3
=
Mpa
0
3
=
Mpa
0
3
=
Mpa
8
3
-
=
Mpa
8
3
-
=
Appear 3.4-i
Appear 3.4-j
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
19/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
3.5 Work hardening
As we already mentioned in the paragraph [§3.3], when the state of stress reaches the edge
field of reversibility, the plastic deformations and the variables intern develop, them
thresholds move: they are hammer-hardened. For our model, the variables intern are
two, they are noted
p
C
for the internal variable “known as of compression” and
p
T
for that “known as of
traction
”. These variables determine the evolution of the thresholds of compression and traction
respectively, the thermodynamic forces theirs are connected by the relations:
()
p
C
C
C
C
F
F
With
-
=
éq
3.5-1
and
()
p
T
T
T
T
F
F
With
-
=
éq 3.5-2
where
()
p
C
C
F
and
()
p
T
T
F
the values of resistances in compression and traction represent
respectively.
3.5.1 Functions
of work hardening
The function
()
p
C
C
F
is initially increasing then decreasing, the decreasing part being is linear
[Figure 3.5.1-a], that is to say quadratic [Figure 3.5.1-b],
f'
C
E
F
C
C
U
()
F
C
C
p
C
p
Appear 3.5.1-a
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
20/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
E
F
C
C
U
f'
C
()
F
C
C
p
C
p
Appear 3.5.1-b
is a data of the model. The shape of the curve enters
E
C
and
U
C
(negative work hardening) depends on
the element, and more precisely of its dimensions, according to a criterion similar to that chosen by
G. Heinfling, [Error! Source of the untraceable reference.] for the taking into account of the localization of
deformations.
In traction, shape of the curve giving the value of the elastic limit
()
p
T
T
F
according to
cumulated plastic deformation
p
T
do not comprise a part “pre-peak”, the part “post-peak” being is
linear [Figure 3.5.1-c], that is to say exponential [Figure 3.5.1-d].

F
T
G
T
F
T
U
T
p
()
F
T
T
p
Appear 3.5.1-c
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
21/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
T
U
()
F
T
T
p
F
T
T
p
Appear 3.5.1-d
3.5.2 Curves
of work hardening and modules post peak

3.5.2.1 Model of cracking distributed
The introduction of a behavior softening post-peak into the relations stress-strains
pose a major problem. Under statical stress, beyond a certain level of stress,
corresponding to the starter of the lenitive behavior, equations governing the balance of
structure lose their elliptic nature. These equations of the mechanical problem form one then
system of partial derivative equations evil posed of which the number of solutions is multiple. It
problem results in an not-objectivity compared to the mesh. It results from this a sensitivity
pathological of the numerical solution to the smoothness and the orientation of the mesh.
In order to solve this problem, or at least, to limit the consequences of them on the reliability of the solution
predicted, it is necessary to use techniques known as of regularization. The object of these techniques is
to enrich the mechanical description of the medium, to be able to describe nonhomogeneous states of
deformation, and to preserve the mathematical nature of the problem. One operates this regularization in
introducing, in the law behavior, a characteristic length or internal length, connected to
width of the area of localization. Several techniques are possible to improve description
mechanics of the lenitive medium. They constitute limitings device of localization. The implementation of
these techniques requires in general, of the delicate numerical developments. An approach
intermediate enters the use of the conventional models and the placement of these limitings device of
localization consists in making depend the slope post-peak on the relation stress-strain, of the size
element, so as to dissipate with the rupture a constant energy. This approach constitutes one
not towards a nonlocal description of the continuous medium.
Let us consider initially a real fissure of surface
S
whose measurement is A [Figure 3.5.2.1-a].
S
is one
surface discontinuity of the field of displacement
U
. It is supposed that to create this discontinuity,
it is necessary to spend an energy
W
whose expression is:
S
)
(
=
W
dS
G
F
X
,
F
G
being a property of
material.
Let us consider now that one wants to represent the same phenomenon, while representing not
a discontinuity of displacement but a plastic deformation uniformly distributed in one
volume
V
.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
22/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Dissipated energy will be:
V
0
=
W
R
T
p
ij
ij
dt
dt
D
FD
, where one noted
R
T
“time-to-failure”.
S
V
L
worm
Appear 3.5.2.1-a
By making the series of following assumptions:
·
the fissure is plane,
·
F
G
is constant along the fissure and thus
F
G
With
.
=
W
,
·
V
is a basic cylinder
S
and a height
worm
L
,
·
=
T
T
p
ij
ij
F
dt
dt
D
G
0
is constant in
V
.
One leads finally to the relation:
F
T
p
ij
ij
F
G
With
dt
dt
D
V
Vg
W
T
.
0
=
=
=
éq
3.5.2.1-1
Or:
worm
F
T
p
ij
ij
F
L
G
dt
dt
D
G
T
=
=
0
éq
3.5.2.1-2
It is seen easily that:
()
U
D
F
G
F
=
0
, writing in which quantities
(
)
,
,
,
U
F
F
G
represent respectively
(
)
p
T
U
T
T
F
T
F
G
,
,
,
in traction and
(
)
p
C
U
C
C
F
C
F
G
,
,
,
in compression.
data of
F
G
thus determine
U
, this in traction as in compression:
()
worm
F
p
C
p
C
F
L
G
D
F
G
C
U
C
C
=
=
0
()
worm
F
p
T
p
T
F
L
G
D
F
G
T
U
T
T
=
=
0
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
23/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Quantity
F
G
is thus related to the slope of the curve post peak in a diagram stress-variable
of work hardening, which is related to the forced slope post peak in a diagram deformation.
Let us suppose for example that the forced relation deformation is linear in mode post peak.
Let us call
0
<
T
E
the slope post peak in the diagram
()
,
and
0
<
H
the corresponding slope in
of diagram
()
,
F
[Figure 3.5.2.1-b]. There is the relation
H
E
He
E
E
E
EE
H
T
T
T
+
=
-
=
who show
that one must have:
E
H
<
-
, or else the diagram
()
,
present a snap back.
F
T
T
F
T
G
T
F
E
E
T
H
EE
E
E
T
T
= -
Appear 3.5.2.1-b
The condition
E
H
<
-
is known as condition of applicability, it will result in an inequality on
F
G
and thus
on
worm
L
.
Within the framework of a resolution by the finite element method, representative elementary volume
fissured medium can be compared to an element of the mesh. The characteristic length (noted by
continuation L
C
) introduced into the method of the energy of equivalent rupture corresponds to the Lver length.
During a calculation corresponding to an unspecified structure, determination this length
characteristic is delicate. It depends on the position of the plan of fissure, dimensions and the type
elements…
A simple estimate for the two-dimensional cases can be expressed in the form:
E
C
With
R
L
=
where A
E
is the surface of the element considered, and R, a correct factor, being worth 1 for
quadratic elements, and
2
for the linear elements.
One can extend this formulation to the case 3D:
3
E
C
V
R
L
=
where V
E
indicate the volume of the element.
Concerning the evolution of work hardening with the temperature, we regard as in [bib2] that
energies of rupture and resistances to rupture depend not on the current temperature
T
material point considered at time
T
, but of the maximum temperature reached in this point since
the beginning of the loading until time
T
. When we need to show the dependence of
quantities compared to the temperature, we will note:
·
the maximum of the temperature since the beginning of loading,
·
()
C
F
resistance in compression,
·
()
T
F
indicate resistance in traction,
·
()
p
C
C
F
,
the curve of work hardening in compression,
·
()
p
T
T
F
,
the curve of work hardening in traction.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
24/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A

3.5.2.2 Behavior of the concrete in traction and linear curve post-peak
In this modeling, the concrete is supposed to be elastic until its resistance in traction
T
F
.
curve
()
p
T
T
F
in traction is represented on [Figure 3.5.1-c] and is entirely defined by
resistance in traction of material, the energy of cracking
F
T
G
, and the characteristic length L
C
.
The mathematical expression of this curve is:
()
()
()




-
=
U
T
p
T
T
p
T
T
F
F
1
,
éq
3.5.2.2-1
The equivalence of dissipated energy makes it possible to write:
()
()
()
()




-
=
=
U
T
U
T
p
T
U
T
p
T
T
C
p
T
T
C
F
T
D
F
L
F
L
G
0
0
1
,
from where
()
() ()
2
.
.
=
U
T
T
C
F
T
F
L
G
éq
3.5.2.2-2
and
()
()
()
=
T
C
F
T
U
T
F
L
G
.
.
2
éq
3.5.2.2-3
The condition of applicability is written:
() ()
()
2
.
.
2
T
F
T
C
F
G
E
L
éq
3.5.2.2-4

3.5.2.3 Behavior of the concrete in traction and exponential curve post-peak
In this modeling, the concrete is supposed to be elastic until its resistance in traction
T
F
.
curve
()
p
T
T
F
in traction is represented on [Figure 3.5.1-d] and is entirely defined by
resistance in traction of material, the energy of cracking
F
T
G
, and the characteristic length L
C
.
The mathematical expression of this curve is:
()
()
()




-
=
U
T
p
T
T
p
T
T
has
F
F
exp
.
,
éq
3.5.2.3-1
The equivalence of dissipated energy makes it possible to write:
()
()
()
()




-
=
=
0
0
exp
,
p
T
U
T
p
T
T
C
p
T
T
C
F
T
D
has
F
L
F
L
G
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
25/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
From where:
()
() ()
has
F
L
G
U
T
T
C
F
T
.
.
=
éq
3.5.2.3-2
and
()
()
()
=
T
C
F
T
U
T
F
L
G
has
.
éq
3.5.2.3-3
That is to say still:
()
()
() ()




-
=
F
T
p
T
T
C
T
p
T
T
G
F
L
F
F
.
exp
.
,
The maximum slope of the curve is then
()
()
()
-
=
F
T
T
C
G
F
L
H
2
max
.
and the condition of applicability is written:
() ()
()
2
.
T
F
T
C
F
G
E
L
éq
3.5.2.3-4

3.5.2.4 Behavior of the concrete in compression and linear curve post-peak
In this modeling, the behavior of the concrete is supposed to be elastic until the elastic limit,
data by a proportionality factor (noted
)
expressed as a percentage of resistance to the peak
()
C
F
.
For the standard concretes
is about 30%. The curve
()
p
C
C
F
in compression is represented
on [Figure 3.5.1-a] and is entirely defined by resistance in traction of material, the energy of
cracking
F
C
G
, and the characteristic length L
C
.
The mathematical expression of this curve is:
()
()
(
) () () ()
()
(
)
()
()
(
)
()
()
(
)
()
()








-
-
=




-
+
-
+
=
-
-
2
3.5.2.4
éq
1
3.5.2.4
éq
U
C
p
C
E
U
C
E
U
C
p
C
C
C
C
E
p
C
E
p
C
E
U
C
C
p
C
C
F
F
F
F
if
if
'
,
1
2
2
'
,
2
2
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
26/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Resistance in maximum compression is reached when:
() (
) ()
()
-
=
E
F
C
E
2
2
The equivalence of dissipated energy makes it possible to write:
()
()
p
C
p
C
C
C
F
C
D
F
L
G
U
C
=
0
,
from where
()
()
()
()


+
+
=
U
C
E
C
C
F
C
F
L
G
2
1
6
1
2
.
éq
3.5.2.4-3
and
()
()
()
()
+
-
=
E
C
C
F
C
U
C
F
L
G
3
1
2
.
.
2
éq
3.5.2.4-4
The slope of the curve is then
()
()
()
()
-
-
=
E
U
C
C
F
H
and the condition of applicability is written:
() ()
()
-
-
2
2
4
4
11
6
.
C
F
C
C
F
G
E
L
éq
3.5.2.4-5

3.5.2.5 Behavior of the concrete in compression and nonlinear curve post-peak
In this modeling, the behavior of the concrete is supposed to be elastic until the elastic limit,
data by a proportionality factor (noted
)
expressed as a percentage of resistance to the peak
()
C
F
.
For the standard concretes
is about 30%. The curve
()
p
C
C
F
in compression is represented
on [Figure 3.5.1-b] and is entirely defined by resistance in traction of material, the energy of
cracking
F
C
G
, and the characteristic length L
C
.
The mathematical expression of this curve is:

()
()
(
) () () ()
()
(
)
()
()
(
)
()
()
(
)
()
()








-
-
-
=




-
+
-
+
=
-
-
2
3.5.2.5
éq
1
3.5.2.5
éq
U
C
p
C
E
E
U
C
E
p
C
C
C
C
E
p
C
E
p
C
E
p
C
C
p
C
C
F
F
F
F
if
if
2
2
2
2
1
,
1
2
2
,
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
27/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Resistance in maximum compression is reached when:
() (
) ()
()
-
=
E
F
C
E
2
2
The equivalence of dissipated energy makes it possible to write:
()
()
p
C
p
C
C
C
F
C
D
F
L
G
U
C
=
0
,
from where:
()
()
()
()


+
=
E
U
C
C
C
F
C
F
L
G
3
3
2
.
éq
3.5.2.5-3
and:
()
() () ()
-
=
E
C
C
F
C
U
C
F
L
G
2
.
2
3
éq
3.5.2.5-4
The maximum slope of the curve post-peak is then
()
()
()
()
-
-
=
E
U
C
C
F
H
.
2
max
and the condition of applicability is written:
() ()
()
-
-
2
2
4
1
.
2
3
C
F
C
F
G
E
L
C
éq
3.5.2.5-5


4 Flow
plastic
In this paragraph, we give the expression speeds of plastic deformation, while distinguishing
the case says general where the state of stress is located on a “regular” area of the edge of the field of
reversibility and the case where it is at the top of one of the cones.
4.1
General form of the rule of normality
In space
)
(
With
,
, the inequalities [éq 3.2-2], [éq 3.2-3], [éq 3.2-4], [éq 3.2-5], define one
convex field which we will note
)
,
(A
C
. We will note
C
the indicating function of this convex:


=
if not
)
(
if
0
)
(
)
,
(A
C
C
With
,
With
,
éq
4.1-1
When the border of the field of reversibility is reached, of the irreversible plastic deformations
develop, according to the classical theory of plasticity.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
28/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
For a standard material, [bib4] the law of flow checks the principle of maximum plastic work,
what results in the equation:
()
C
p
&
&
,
éq 4.1-2
where
C
note under differential of the function
C
. We remind the meeting [bib3] that under differential of one
convex function in a point
X
is the whole of the vectors
Z
such as
:
()
()
*
*
*
,
X
X
X
Z
X
F
X
F
-
+
It is then seen easily that [éq 4.1-2] involves:
(
)
()
(
) (
)
*
*
*
*
*
*
With
With
With
With
,
With
,
and
-
+
-
+
&
&
p
C
C
éq
4.1-3
Taking into account the definition of the characteristic function, one to see easily that [éq
4.1-3]
is
equivalent with:
(
)
With
p
p
C
,
and
+
+
*
*
*
*
With
With
With
&
&
&
&
éq
4.1-4
In other words the plastic flow is such as the couple
()
With
,
carry out the maximum of
plastic dissipation among the acceptable thermodynamic forces.

4.2
Expression of the plastic flow partly current
When the function
F
is differentiable at the point considered
()
,
the rule of normality is written
simply
=
F
p
&
&
éq
4.2-1
With
=
F
&
&
éq
4.2-2
&
and
F
checking the conditions of Kuhn-Tucker:


=
0
.
0
0
F
F
&
&
éq
4.2-3
The variable of work hardening is related to the plastic multiplier by the law of work hardening. By using it
plastic work, one can write:
p
F
&
&
=
.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
29/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
If
F
is a homogeneous function of command 1 compared to the tensorial variable
, one has
F
F
=
, it
who leads to the equality:
=
&
&
and thus finally with the equations:
=
comp
p
C
p
C
F
&
&
éq
4.2-4
=
trac
p
T
p
T
F
&
&
éq
4.2-5
4.3
Expression of the plastic flow at the top of a cone
We give two presentations of the same result. The first presentation uses the theory of
standard materials generalized and under differentials, the second share of an equality posed a priori
on plastic work.
4.3.1 Demonstration by the general theory of standard materials
The field
(
)
With
C
,
consists of two cones. The function
C
is not differentiable is with
the intersection of these two cones, is at the top of each one of these cones. When the point
()
With
,
belongs to the intersection of the two cones, the preceding equations remain valid, with
precision that the deformations figure of compression and traction develop into same
time. This case known as “multi criterion” is remainder treated in [bib4]. We will be satisfied here to treat
the case where
()
With
,
is at the top of a cone, and we will choose the most frequent case of the node of
cone of traction, knowing that the case of the node of the cone of compression is treated exactly
even way.
The criteria are rewritten by using the variables
eq
and
H
, more practical in the development
analytical.
0
3
2
)
,
(
+
-
+
=
T
T
H
eq
T
trac
With
F
D
C
D
With
F
4.3.1-1
(
)
0
,
+
-
=
T
T
H
T
H
With
F
D
C
With
F
trac
4.3.1-2
We thus consider a case where:


=
+
-
=
0
0
T
T
H
eq
With
F
D
C
4.3.1-3
On the basis of [éq 4.1-4], we will calculate plastic dissipation as being the maximum of
*
*
With
&
&
+
p
for all the couples
(
)
With
C
,
,
*
*
With
(
)
(
)
*
*
With
With
*
*
&
&
+
=
p
C
p
Max
With
D
,
,
4.3.1-4
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
30/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
By writing whereas this maximum is finished and reached when
*
=
and
With
With
*
=
, we will find
conditions on
p
&
and
&
. In fact, the finished character will be enough.
By using the partly isotropic decomposition of the tensors and déviatoire, and the particular form of
variables of work hardening, one finds easily:
*
*
*
3
~
T
p
T
p
H
H
p
p
With
+
+
=
+
&
&
&
&
&
S
With
*
*
4.3.1-5
Let us consider the unit then
1
vectors forced of null trace and of which the stress
equivalent of Von Mises 1 is worth:
()
{
}
0
trace
,
1
,
1
=
=
=
eq
(
)
(
)




+
-
+
-
+
+
=
0
0
3
2
1
1
1
,
T
T
H
T
T
H
eq
H
eq
With
With
F
D
C
With
F
D
C
D
C
S
I
S
With
,
4.3.1-6
In other words, the “direction” of the diverter of the stresses is unspecified for a couple
(
)
(
)
With
C
,
,
With
.
One can thus write:


+
+
=




+
-
+
-
+
*
*
* 1
*
3
~
1
1
0
*
*
0
*
*
*
3
2
1
* 1
,
*
,
*
,
*
T
p
T
p
H
H
p
eq
p
With
Max
D
Max
T
With
T
F
H
D
C
T
With
T
F
H
D
C
eq
D
T
With
H
eq
&
&
&
S
S
S
4.3.1-7
It is clear that the maximum of
*
1
~
S
p
&
is reached when
*
1
S
is “parallel” with
p
&
~
and that one has then:
p
p
eq
Max
=
&
&
~
3
2
*
1
~
1
1
S
S
.
[éq 4.3.1-7] can thus be written:
(
)
(
)












+
+
+


=












+
-
+
-
+
+
-
+
-
+
+
-
+
-
+
*
*
*
0
*
*
0
*
*
*
3 2
*
,
*
,
*
0
*
*
0
*
*
*
3 2
*
,
*
,
*
0
*
*
0
*
*
*
3 2
*
,
*
,
*
3
~
3
2
T
p
T
p
H
eq
p
p
With
D
Max
Max
Max
T
With
T
F
H
cd.
T
With
T
F
H
cd.
eq
D
T
With
H
eq
T
With
T
F
H
cd.
T
With
T
F
H
cd.
eq
D
T
With
H
eq
T
With
T
F
H
cd.
T
With
T
F
H
cd.
eq
D
T
With
H
eq
H
eq
&
&
&
4.3.1-8
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
31/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Like
0
~
p
eq
&
, one has for the first term:


-
+
-
=






+
-
+
-
+
*
*
*
2
3
2
3
2
3
~
3
2
~
3
2
0
*
*
0
*
*
*
3 2
*
,
*
,
*
H
T
eq
eq
C
F
D
With
D
T
p
eq
p
Max
T
With
T
F
H
cd.
T
With
T
F
H
cd.
eq
D
T
With
H
eq
&
&
4.3.1-9
[éq 4.3.1-9] deferred in [éq 4.3.1-8] gives:
(
)
(
)
(
)
(
)
p
p
T
T
p
p
H
T
p
p
eq
eq
H
eq
D
With
C
F
D
D
Max
Max
T
With
T
F
H
cd.
T
With
H
T
With
T
F
H
cd.
T
With
H
-
+
-
+
=
+
-
+
-
&
&
&
&
&
~
2
~
2
3
~
2
*
*
0
*
*
*
,
*
0
*
*
*
,
*
4.3.1-10
Let us pose then:
p
p
H
eq
C
m
-
=
&
&
~
2
3
,
p
p
T
eq
D
N
-
=
&
&
~
2
and
F
D
Q
p
eq
=
&
~
2
With these notations, [éq 4.3.1-10] becomes:
(
)
*
*
0
*
*
*
,
*
T
p
Na
m
Q
D
H
Max
T
With
T
F
H
cd.
T
With
H
+
+
=
+
-
4.3.1-11
It is about a problem of the type “simplex”. The field of
*
*
,
T
H
With
is represented on
[Figure 4.3.1-a].
*
H
*
T
With
T
F
C
D
T
F
*
*
,
T
H
With
Field
Appear 4.3.1-a
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
32/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Like the field
*
*
,
T
H
With
extends towards
-
at the same time for
*
H
and
*
T
With
, so that
p
D
that is to say finished, it is necessary
that
m
and
N
are positive. The maximum of
*
*
T
Na
m
H
+
is reached for a couple
*
*
,
T
H
With
located on
edge of the field of
*
*
,
T
H
With
.
One has then:






-
+
+
=
D
C
N
m
F
N
Q
D
H
H
Max
T
p
*
*
So that
p
D
that is to say finished, it is necessary that:
D
C
N
m
=
Taking again the definitions of
m
and
N
, this relation gives:
p
T
p
H
D
C
=
&
&
3
éq 4.3.1-12
In addition, stresses
0
m
and
0
N
give:
p
p
H
eq
C
&
&
~
2
3
éq 4.3.1-13
and
p
p
T
eq
D
&
&
~
2
éq 4.3.1-14
these two last inequalities being equivalent because of [éq 4.3.1-11].
The equations [éq 4.3.1-11] and [éq 4.3.1-12] define the plastic flow in the node of one of
cones of the field of reversibility.
4.3.2 Demonstration by plastic work
The starting point is to consider that compared to the developments made in regular points,
they are primarily the relations [éq 4.2-1] and [éq 4.2-2], known as rules of normality, which cannot
to be written more. However the relation [éq 4.2-1] implies the equality
()
p
p
T
T
p
T
F
&
&
=
, which can, it, being
maintained.
We will thus leave the equation:
()
p
p
T
T
p
T
F
&
&
=
éq
4.3.2-1
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
33/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
We use the partly isotropic decomposition and déviatoire tensors and find:
()
H
H
p
p
T
T
p
T
F
+
=
&
&
&
3
~
S
éq
4.3.2-2
At the top of the cone of traction, there are the relations [éq 4.3.1-3], which, carried in [éq 4.3.2-2] give,
while also using [éq 3.5-2]:
()
()
H
p
T
T
p
T
T
p
T
F
C
D
F
=
&
&
3
éq
4.3.2-3
And one thus finds the relation [éq 4.3.1-12]:
p
T
H
D
C
=
&
&
3
.
4.4
Together equations of behavior (summarized)
One notes
H
the matrix of elasticity:














µ
µ
µ
µ
+
µ
+
µ
+
=
2
0
0
0
0
0
0
2
0
0
0
0
0
0
2
0
0
0
0
0
0
2
0
0
0
2
0
0
0
2
H
With:
(
) (
)
-
+
=
2
1
1
E
and
(
)
+
=
µ
1
2
E
, and
3
2
3
µ
+
=
K
The forced relations deformations are written finally:
(
)
p
T
p
C
H
-
-
=
éq
4.4-1
For a regular point of the cone of compression:
(
)
0
3
2
,
+
-
+
=
C
C
H
eq
C
comp
With
F
B
has
B
With
F
éq
4.4-2
=
=
comp
p
C
p
C
comp
p
C
F
F
&
&
&
;
0
éq
4.4-3
For a regular point of the cone of traction:
0
3
2
)
,
(
+
-
+
=
T
T
H
eq
T
trac
With
F
D
C
D
With
F
éq
4.4-4
=
=
trac
p
T
p
T
trac
p
T
F
F
&
&
&
;
0
éq
4.4-5
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
34/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
For a point at the top of the cone of compression:
0
=
S
éq 4.4-6
(
)
0
,
=
+
-
=
C
C
H
C
H
With
F
B
has
With
F
comp
éq
4.4-7
p
C
H
p
C
B
has
=
&
&
3
éq 4.4-8
eq
p
C
H
p
C
has
&
&
~
2
3
éq 4.4-9
For a point at the top of the cone of traction:
0
=
S
éq 4.4-10
(
)
0
,
=
+
-
=
T
T
H
T
H
With
F
D
C
With
F
trac
éq
4.4-11
p
T
H
p
T
D
C
=
&
&
3
éq 4.4-12
eq
p
T
H
p
T
C
&
&
~
2
3
éq 4.4-13
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
35/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
5
Numerical integration of the law of behavior
5.1
The total problem and the local problem: recalls
For a given structure (geometry and material), and for a given loading, fields of
displacement, stress and variables intern are by solving a whole of equations with
nonlinear derivative partial formed starting from and the law equilibrium equations of behavior.
The document [bib5] presents the algorithm of which we give a summary here:
0
U
and
0
known
Loop urgent
I
T
: loading
()
I
I
T
L
L
=
1
-
I
U
known; calculation of the prediction
0
I
U
Iterations of balance of Newton N
N
I
U
known;
1
-
-
=
I
N
I
N
I
U
U
U
Loop elements el
Loop points of gauss G
calculation
()
N
I
el
G
N
I
el
G
U
=
law of behavior:
calculation of:
N
I
el
G
and
N
I
el
G
from
1
-
I
el
G
,
1
-
I
el
G
and
N
I
el
G
calculation of
N
I
el
G
N
I
el
G
(
according to option
)
Accumulation in vectors and matrices assembled:
Accumulation of
N
I
el
G
el
G
T
Q
.
in
N
I
T
Q
.
Accumulation of calculation of
el
G
N
I
el
G
N
I
el
G
el
G
T
Q
Q
in
N
I
K
(
according to option
)
Calculation of
N
I
U
by:
I
N
I
T
N
I
N
I
L
Q
U
K
+
-
=
.
.
linear iteration of search to determine
Updating:
N
I
N
I
N
I
U
U
U
+
=
+
1
IF test convergence OK
fine Newton: no time following I = i+1
If not
N = n+1
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
36/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
The calculation of the stresses and variables internal
N
I
el
G
N
I
el
G
with the iteration of Newton
N
and at time
I
T
with
to leave the internal stresses and variables
1
-
I
el
G
,
1
-
I
el
G
at time
1
-
I
T
and of the value and
N
I
el
G
of
the increase in deformation in the interval of time estimated with the iteration of Newton
N
consist with
to integrate the equations [éq 4.4-1], [éq 4.4-2] with [éq 4.4-5] or [éq 4.4-6] with [éq 4.4-9] or [éq 4.4-10] with
[éq 4.4-13] according to the cases with the initial conditions:
()
1
1
-
-
=
I
el
G
I
T
éq
5.1-1
()
1
1
-
-
=
I
el
G
I
T
éq
5.1-2
()
1
1
-
-
=
I
el
G
p
I
p
T
éq
5.1-3
With the condition of loading in imposed deformation:
()
N
I
el
G
I
T
=
éq
5.1-4
The result of this integration will provide:
()
()
()
I
N
I
el
G
I
p
N
I
el
G
p
I
N
I
el
G
T
T
T
=
=
=
The object of this chapter is to present the numerical integration of these equations. It is about a system
nonlinear differential equations which we solve by a method of implicit Euler. To leave
of now, the quantity at the beginning of the pitch of time (known) will be noted with an index
-
, then
that unknown factors at the end of the pitch of time (all unknown factors except
N
I
el
G
=
) will be noted without index.
For an unspecified quantity
has
one notes
-
-
=
has
has
has
.
One always starts by calculating an elastic solution
E
, by supposing that there is no evolution
plastic deformations and internal variables. So at least one of the criteria is violated by this
elastic solution, it is necessary to calculate plastic flows. The cases should then be distinguished
regular for which the solution
is on the regular part of one of the cones or with their intersection
cases known as singular where the solution
is at the top of one of the two cones. Logic allowing
to examine and choose these various cases, and the algorithm which results from this are relatively complex.
We thus present in first the processing of each case and explain their sequence
subsequently, in the chapter [§5.7].
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
37/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
5.2
Digital processing of the regular case.
One presents in detail only the case where at the same time plastic deformations develop in
traction and in compression and where thus the solution
belongs to the intersection of the two cones. Let us note
however that, even if
E
violate at the same time the two criteria, for as much the final solution
can very
to belong well finally only to one of the cones hammer-hardened. One is thus brought to seek balanced
which one postulates that they belong to one of the two cones or both. The case or it belongs to
only one of the two cones results easily from the more general case presented here. Equations that us
let us have to solve are finally:
µ
+
µ
µ
=
-
-
~
2
S
S
E
éq
5.2-1
H
H
E
H
K
K
K
+
=
-
-
3

éq
5.2-2
(
)
p
T
p
C
E
S
S
~
~
2
+
µ
-
=
éq
5.2-3
(
)
T
p
H
C
p
H
E
H
H
K
+
-
=
3
éq
5.2-4
()
()
0
3
2
,
=
-
+
=
p
C
C
H
eq
p
C
comp
F
B
has
B
F
éq
5.2-5
=
comp
p
C
p
C
F
éq
5.2-6
0
)
(
3
2
)
,
(
=
-
+
=
p
T
T
H
eq
p
T
trac
F
D
C
D
F
éq
5.2-7
=
trac
p
T
p
T
F
éq
5.2-8
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
38/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
By taking the isotropic and deviatoric parts plastic deformations, the equations [éq 5.2-6]
and [éq 5.2-8] give:
eq
p
C
p
C
B
=
S
2
~
éq
5.2-9
B
has
p
C
C
p
H
3
=
éq 5.2-10
eq
p
T
p
T
D
=
S
2
~
éq 5.2-11
D
C
p
T
T
p
H
3
=
éq 5.2-12
While deferring [éq 5.2-9] and [éq 5.2-11] in [éq 5.2-3], one finds:
eq
p
T
p
C
E
D
B




+
µ
-
=
S
S
S
2
2
2
éq
5.2-13
who shows that
S
is parallel to
E
S
from where one deduces:
eq
E
E
eq
=
S
S
éq 5.2-14
While deferring [éq 5.2-14] in [éq 5.2-9] and [éq 5.2-11] one finds:
eq
E
E
p
C
p
C
B
=
S
2
~
éq 5.2-15
eq
E
E
p
T
p
T
D
=
S
2
~
éq 5.2-16
One defers then [éq 5.2-15] and [éq 5.2-16] in [éq 5.2-3] and [éq 5.2-4], and one expresses the criteria
[éq 5.2-5] and [éq 5.2-7] with these new results. That led to two equations having like
unknown factors
p
C
and
p
T
:
(
)
0
3
2
3
2
3
2
2
2
2
=
+
-


+
µ
-




+
µ
-
+
-
p
C
p
C
C
p
T
p
C
E
H
eq
E
F
data base
Kac
data base
B
Ka
B
B
has
B
éq
5.2-17
(
)
0
3
2
3
2
3
2
2
2
2
=
+
-




+
µ
-


+
µ
-
+
-
p
T
p
T
T
p
T
p
C
E
H
eq
E
F
D
Kc
D
data base
Kac
data base
D
C
D
éq
5.2-18
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
39/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
It is this system of two equations to two unknown factors which should finally be solved. If
functions
C
F
and
T
F
are linear, i.e. if one is in linear mode post-peak in
compression as in traction, it acts of a linear system which will thus be solved in an iteration.
In the case, is mode pre-peak in compression (which is always nonlinear), that is to say of
modelings with nonlinear modes post peak, the system [éq 5.2-17] and [éq 5.2-18] is solved by one
method of Newton:
One notes
(
)
p
T
p
C
comp
F
,
*
the criterion of compression regarded as function of the only variables
p
C
and
p
T
, in the same way for traction:
(
)
(
)
p
C
p
C
C
p
T
p
C
E
H
eq
E
p
T
p
C
comp
F
data base
Kac
data base
B
Ka
B
B
has
B
F
+
-


+
µ
-




+
µ
-
+
=
-
3
2
3
2
3
2
,
2
2
2
*
(
)
(
)
p
T
p
T
T
p
T
p
C
E
H
eq
E
p
T
p
C
trac
F
D
Kc
D
data base
Kac
data base
D
C
D
F
+
-




+
µ
-


+
µ
-
+
=
-
2
2
2
*
3
2
3
2
3
2
,
The ième iteration of Newton for system [éq 5.2-17] - [éq 5.2-18] is:
(
)
(
)
I
p
T
p
C
trac
p
T
p
C
comp
I
I
p
T
p
C
I
p
T
p
C
F
F
J






-




=




-
+
,
,
1
1
The jacobien
I
J
is worth:
I
p
T
trac
p
C
trac
p
T
comp
p
C
comp
I
F
F
F
F
J






=
With:
(
)
p
C
p
C
p
C
C
p
C
comp
F
B
Ka
B
F
+
-




+
µ
-
=
-
2
2
2
3
2


+
µ
-
=
data base
Kac
data base
F
p
T
comp
3
2


+
µ
-
=
data base
Kac
data base
F
p
C
trac
3
2
(
)
p
T
p
T
p
T
T
p
T
trac
F
D
Kc
D
F
+
-




+
µ
-
=
-
2
2
2
3
2
Initial Jacobien of the system results from the values of derived in
0
=
p
C
and
0
=
p
T
, which
amounts solving the nonlinear system on the basis of the null solution. Nonthe linearities are
introduced by the curves of softening. In the post-peak part, when they are linear,
convergence is done in an iteration. When they are nonlinear, convergence only requires
some iterations. To leave the null solution thus does not pose a problem of convergence. That
returns starting from the linearization of the criteria in the vicinity of the elastic prediction.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
40/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
5.3
Existence of a solution and condition of applicability
We remind the meeting that the solution of the problem [éq 5.2-17] and [éq 5.2-18] must check the conditions
[éq 4.2-3] and thus inter alia the positivity of the increases in the plastic multipliers.
0
p
C
éq
5.3-1
0
p
T
éq
5.3-2
Let us suppose that we are in a case of linear behavior post peak in traction as in
compression and let us call respectively
C
H
and
T
H
slopes of the parts post peak. Increases
plastic multipliers are obtained by solving the linear system:
()
()






-
+
-
+
=










+
+
µ
+
µ
+
µ
+
+
µ
-
-
p
T
T
E
H
eq
E
p
C
C
E
H
eq
E
p
T
p
C
T
C
F
D
C
D
F
B
has
B
H
D
Kc
D
data base
Kac
data base
data base
Kac
data base
H
B
Ka
B
3
2
3
2
3
2
3
2
3
2
3
2
2
2
2
2
2
2
éq
5.3-3
Since the criteria of traction and compression were activated in traction as in compression,
the second member of this system is positive. But nothing ensures in so far as the solution of
[éq 5.3-3] will be positive.
If one poses:






+
+
µ
+
µ
+
µ
+
+
µ
=
T
C
H
D
Kc
D
data base
Kac
data base
data base
Kac
data base
H
B
Ka
B
2
2
2
2
2
2
3
2
3
2
3
2
3
2
HTC
éq
5.3-4
One a:
()
()






-
+
-
+
=




-
-
-
p
T
T
E
H
eq
E
p
C
C
E
H
eq
E
p
T
p
C
F
D
C
D
F
B
has
B
3
2
3
2
1
HTC
éq
5.3-5
With:
(
)
(
)
(
)
K
H
H
K
H
K
H
H
H
K
H
K
H
T
C
T
C
T
C
T
T
9
16
4
3
9
4
9
2
3
6
3
3
2
+
+
µ
+
+
+
+
µ
+
-
µ
+
=
HTC
éq
5.3-6
(
)
()






-
µ
+
-
+
µ
-
µ
+
-
µ
-
µ
+
-
µ
+
+
=
-
1
12
18
9
3
2
3
2
9
2
9
3
2
3
2
9
2
9
3
4
3
1
K
K
H
K
K
K
K
K
H
C
T
HTC
HTC
1
éq
5.3-7
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
41/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
It is seen that conditions of positivity [éq 5.3-1] and [éq 5.3-2] lead to relations relatively
intricate. If the solution of the problem [éq 5.2-17] and [éq 5.2-18] does not check the conditions of
positivity [éq 5.3-1] and [éq 5.3-2], that can correspond is with the fact that the coefficients
C
H
and
T
H
are
such as there is no positive solution (that would correspond to a snap-back in a diagram
()
,
), that is to say with the fact that the solution activates finally only one of the two criteria.
Let us examine the simpler case of only one activated criterion. Let us suppose to fix the notations that only
activated criterion is the criterion of traction.
One must have:
()
-
-
+
=




+
+
µ
p
T
T
E
H
eq
E
T
p
T
F
D
C
D
H
D
Kc
D
3
2
3
2
2
2
2
éq
5.3-8
One sees reappearing the condition known as of applicability:
0
3
2
2
2
2
>




+
+
µ
T
H
D
Kc
D
éq
5.3-9
This condition is the generalization of the condition
E
H
<
-
presented at the paragraph [§3.5.2.1] in
a particular case of axial stress plain.
The following strategy will thus be retained:
If
()
0
3
2
>
-
+
-
p
C
C
E
H
eq
E
F
B
has
B
and
()
0
3
2
>
-
+
-
p
C
C
E
H
eq
E
F
B
has
B
Activation a priori of the two criteria: resolution problem [éq 5.2-17] and [éq 5.2-18]
So not convergence or so not checking conditions of positivity [éq 5.3-1] and [éq 5.3-2]
Seek solution with only one activated criterion
So not convergence or not checking condition positivity
Stop on diagnostic of nonchecking of condition of applicability of the type [éq 5.3-9]
5.4
Processing of the nonregular cases
In this paragraph, we describe the discrete processing of the equations corresponding to projection
at the node of the cone of traction, [éq 4.4-10] with [éq 4.4-13], knowing that projection at the top of
cone of compression is done in the same way. The equations [éq 4.4-10] and [éq 4.4-12] define
the plastic flow in this case, whereas the equation [éq 4.4-13] is a condition of acceptability
projection at the top of the cone.
5.4.1 Calculation of the stresses and plastic deformations
Discrete forms of [éq 4.4-10] with [éq 4.4-13], are:
0
=
S
éq
5.4.1-1
(
)
p
T
p
T
T
H
F
C
D
+
=
-
éq
5.4.1-2
p
T
H
p
T
D
C
=
3
éq
5.4.1-3
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
42/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
The relation [éq 5.2-4] established in the regular case is always valid, one jointly uses it with
[éq 5.4.1-3] in [éq 5.4.1-1] and one obtains:
(
)
p
T
p
T
T
p
T
E
F
C
D
K
D
C
H
+
=
-
-
éq
5.4.1-4
The relation [éq 5.4.1-4] is a nonlinear equation compared to the variable
p
T
that one solves by
an algorithm of Newton, which makes it possible to calculate
H
p
T
by [éq 5.4.1-3] and
H
by [éq 5.4.1-2].
Taking into account [5.4.1-1], the stresses are thus completely known. [éq 5.4.1-1] gives
still:
0
~
2
=
µ
-
=
p
T
E
S
S
éq
5.4.1-5
This last equation makes it possible to calculate
p
T
~
and the plastic deformations are completely
known.
5.4.2 Acceptability
The discrete form of the relation [éq 4.4-13] is:
eq
p
T
H
p
T
C ~
2
3
éq
5.4.2-1
[éq 5.4.1-5] gives:
µ
=
2
~
eq
E
eq
p
T
éq
5.4.2-2
While using [éq 5.2-4] and [éq 5.4.2-2], [éq 5.4.2-1] is written:
H
E
eq
E
H
cK
-
µ
2

éq
5.4.2-3
5.4.2.1 Acceptability a priori and a posteriori
For the criterion of traction and the part post peak of the criterion of compression,
(
)
p
T
p
T
T
H
F
C
D
+
=
-
is a decreasing function of the variable of work hardening
p
T
. One deduces from it that
H
E
E
H
H
H
-
-
-
and thus that:
H
E
eq
E
E
eq
E
H
H
H
cK
cK
-
µ
-
µ
-
2
2
The condition
-
-
µ
H
H
E
eq
E
cK
2
is known as condition of acceptability a priori because it can be calculated
as of the elastic prediction. The condition
H
E
eq
E
H
cK
-
µ
2
is known as condition of acceptability has
posteriori.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
43/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Direction of projection
area of projection
at the node of the cone
Of traction
H
eq
(
)
eq
E
E
H
,
T
P
Figure 5-1
These conditions have a simple graphic interpretation. One can see easily that, in the case of one
regular solution, one a:
cK
H
E
eq
eq
E
H
µ
=
-
-
2
That shows that the solution in stress is obtained by projecting the point
(
)
eq
E
E
H
,
parallel to one
direction
(
)
µ
2
,
cK
in a diagram
(
)
eq
H
,
, as indicated on [Figure 5.4.2-a]. Areas
of acceptability of projection at the top are cones of which the node and that of the cone of
reversibility and delimited on the one hand by the axis
T
H
COp
>
and by a half-line resulting from the same point and
of direction
(
)
µ
2
,
cK
.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
44/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
5.4.3 Existence of a regular solution and a singular solution.
If projection at the top of the cone is acceptable a posteriori, it may be which exists too
a regular solution as one can see it on [Figure 5.4.2-b].
Direction of projection for
regular case
area of projection
at the node a posteriori
H
eq
(
)
eq
E
E
H
,
area of projection
at the node a priori
-
T
P
1
T
P
2
T
P
Appear 5.4.2-b
The node of the cone of traction before work hardening is noted
-
T
P
, that of the cone hammer-hardened with one
increase in variable of work hardening
1
p
T
is noted
1
T
P
, that of the cone hammer-hardened with one
increase in variable of work hardening
1
2
p
T
p
T
>
is noted
2
T
P
. It is seen that there is a solution
regular with
1
p
T
and a solution with projection at the top of the cone for
2
p
T
. Since
the regular solution corresponds to a less work hardening, in the process of evolution, it will be
met before the solution with projection at the top: it is thus the regular solution which it is necessary
to retain.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
45/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
For this reason the sequence of regular search of solution and with projection at the top is
the following:
I projection at the acceptable top a priori:
-
-
µ
H
H
E
eq
E
cK
2
Calculation of the solution with projection at the top:
p
T
by [éq 5.4-4]
So not
Seek regular solution
So not convergence or not checking condition positivity
Calculation of the solution with projection at the top:
p
T
by [éq 5.4-4]
Checking of acceptability a posteriori:
H
E
eq
E
H
cK
-
µ
2
If not acceptable:
H
E
eq
E
H
cK
-
>
µ
2
Stop on diagnostic of nonchecking of condition of applicability
5.4.4 Inversion of the nodes of the cones of traction and compression
A priori, the node of the cone of compression corresponds to a hydrostatic pressure of traction
much larger than that the node of the cone of traction. But, as one can see it on
[Figure 5.4.4-a], one can find a history of loading which never activates the criterion of traction,
who activates and strongly hammer-hardens the criterion of compression until it to return strictly included in
field of reversibility of the criterion of traction.
H
eq
T
P
Way of stress
Criterion of compression
Criterion of compression hammer-hardened
Appear 5.4.4-a
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
46/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
When the two criteria were thus reversed, but the criterion should not intervene a priori any more of
compression. It may be whereas the solution is a projection at the top of the cone of compression,
who treats himself exactly like projection with the node of the cone of traction.
5.4.5 Projection at the top of the two cones
If the two cones were inverted and if the elastic prediction violates both
criteria, it may be which are finally acceptable at the same time the solution of projection at the top of
cone of compression and at the top of the cone of traction. In these situations, no criterion allows
to select a solution rather than the other and one will thus seek a simultaneous projection with
node of the two cones, which will have to thus share the same node, as indicated on
[Figure 5.4.5-a].
eq
H
Criterion of
traction
initial
Criterion of
compression
initial
Direction of projection
the criterion of compression
Elastic prediction
(
He
,
eq
E
) to project
Solution with projection
on the criterion of traction
Solution with projection
at the node of the two cones
Criteria of
compression and
traction hammer-hardened
Appear 5.4.5-a
The solution with projection at the two tops is obtained by solving the system:
(
)
p
C
p
C
C
p
T
p
C
E
F
has
B
K
D
C
K
B
has
H
+
=
-
-
-
éq
5.4.5-1
(
)
p
T
p
T
T
p
T
p
C
E
F
C
D
K
D
C
K
B
has
H
+
=
-
-
-
éq
5.4.5-2
The state of stress is given by:
0
=
S
(
) (
)
p
C
p
C
C
p
T
p
T
T
H
F
has
B
F
C
D
+
=
+
=
-
-
éq
5.4.5-3
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
47/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
5.5
Determination of the tangent operator
During iterations of the algorithm of Newton-Raphson, it is necessary to calculate the matrix of
tangent stiffness. The construction of the aforementioned plays an important part in stability, the speed and
precision of the method of resolution. To preserve these properties, the matrix of tangent stiffness
must be built starting from an operator binding the increment of stress to the increment of deformation
in a precise way at the end of the process of return on surfaces of load. The matrix of Hooke, thus
that the thermal deformations intervene like constants at the time of the determination of
the coherent tangent operator, built at the end of the iteration in the increment concerned.
The calculation of the operator of coherent tangent behavior takes into account the deformations
plastics. For reasons of simplicity, we chose to calculate the operator of behavior
tangent of speed.
5.5.1 Operator
tangent
in
speed with only one active criterion
In the case of only one active criterion, for example, the criterion in compression, the calculation of the operator of
tangent behavior speed is as follows:
One thus uses the equations of speed, in elastoplastic load:
0
H
=




-
-
comp
p
C
F
&
&
&
éq
5.5.1-1
0
=
+
p
C
p
C
comp
T
comp
F
F
&
&
éq
5.5.1-2
The tangent operator of speed is defined by:
D
&
&
=
éq
5.5.1-3
While identifying [éq 5.5.1-3] with [éq 5.5.1-1] and [éq 5.5.1-2], one finds classically:
H
H
H
D
-
=
comp
T
comp
F
F
1
éq
5.5.1-4
with:
p
C
comp
comp
T
comp
F
F
F
-
=
H
éq
5.5.1-5
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
48/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
5.5.2 Operator
tangent
in
speed with two active criteria
If the two criteria are activated, the criterion in compression and the criterion in traction, it
calculation of the operator of tangent behavior speed is as follows:
One leaves:
0
H
=




-
-
-
trac
p
T
comp
p
C
F
F
&
&
&
&
éq
5.5.2-1
0
=
+
p
C
p
C
comp
T
comp
F
F
&
&
éq
5.5.2-2
0
=
+
p
T
p
T
trac
T
trac
F
F
&
&
éq
5.5.2-3
One leads to:








+
+




+
-
=
T
trac
tt
T
comp
Tc
trac
T
trac
ct
T
comp
DC
comp
F
F
F
F
F
F
H
H
D
éq
5.5.2-4
with:








-




-




-




-
=
H
H
H
H
H
comp
T
trac
trac
T
comp
p
T
trac
trac
T
trac
p
C
comp
comp
T
comp
p
T
trac
trac
T
trac
DC
F
F
F
F
F
F
F
F
F
F
F
F
F
éq 5.5.2-5








-




-




-




-
=
H
H
H
H
H
comp
T
trac
trac
T
comp
p
T
trac
trac
T
trac
p
C
comp
comp
T
comp
trac
T
comp
ct
F
F
F
F
F
F
F
F
F
F
F
F
éq 5.5.2-6








-




-




-




-
=
H
H
H
H
H
comp
T
trac
trac
T
comp
p
T
trac
trac
T
trac
p
C
comp
comp
T
comp
comp
T
trac
Tc
F
F
F
F
F
F
F
F
F
F
F
F
éq 5.5.2-7
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
49/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A








-




-




-




-
=
H
H
H
H
H
comp
T
trac
trac
T
comp
p
T
trac
trac
T
trac
p
C
comp
comp
T
comp
p
C
comp
comp
T
comp
tt
F
F
F
F
F
F
F
F
F
F
F
F
F
éq 5.5.2-8
The expression seems expensive to express in term of products of matrix and calculation. But, when
the operations are made in the order which is appropriate, it is enough to calculate terms little. Moreover, it is
the same terms which intervene on several occasions. It is necessary to calculate the derivative of the criteria by
report/ratio with the stress, and the plastic multipliers, then sums and products
with the actual values, to finish by the constitution of the matrices and theirs let us be.
Lastly, the resulting matrix with the advantage of being symmetrical, which is appropriate for the standard resolution
with Code_Aster.
5.5.3 Derivative successive of the criteria in traction and compression

5.5.3.1 Successive drifts of the criteria compared to the stress
The derivative of the isotropic and deviatoric components of the stresses compared to the tensor of
stresses are expressed in the following way:
By defining the vector














=
0
0
0
1
1
1
0
the derivative of the criteria compared to the tensor of stresses are expressed in the following way:
0
S
B
has
B
F
eq
comp
3
2
+
=
0
D
C
D
S
F
eq
trac
3
2
+
=

5.5.3.2 Successive drifts of the criteria compared to the plastic multipliers
Derived from the criterion of compression in the case of a linear curve post-peak:
(
)
U
C
p
C
E
E
U
T
C
p
C
comp
E
p
C
E
p
C
E
C
p
C
comp
if
F
F
if
F
F




-
=




-
-
=
1
.
1
.
.
3
4
2
Derived from the criterion of compression in the case of a nonlinear curve post-peak:
(
)
U
C
p
C
E
E
U
C
E
p
C
C
p
C
comp
E
p
C
E
p
C
E
C
p
C
comp
if
F
F
if
F
F




-
-
=




-
-
=
2
2
.
.
2
1
.
.
3
4
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
50/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Derived from the criterion of traction in the case of a linear curve post-peak:
U
T
p
T
U
T
T
p
T
trac
if
F
F
=
Derived from the criterion of traction in the case of an exponential curve post-peak:




-




=
U
T
p
T
has
U
T
T
p
T
trac
E
has
F
F
.
5.6
Variables intern model
We install here the internal variables stored in each point of Gauss in the implementation of
model
Internal number of variable
Feel physical
1
p
C
: plastic deformation cumulated in compression
2
p
T
: plastic deformation cumulated in traction
3
: maximum temperature attack at the point of gauss
4 Indicator
of
plasticity
5.7
Top-level flowchart of resolution
The flow chart includes/understands the various stages of the resolution, with the processing of projections
at the nodes of the cones of compression and traction in the following way:
at the beginning of algorithm,
one carries out a projection at the top of the cone of traction:
·
when the elastic prediction checks the condition of projection a priori in traction,
·
when the elastic prediction checks the condition of projection a priori in compression and
that the nodes of the cones of traction and compression were inverted on the axis
hydrostatic,
one carries out a simultaneous projection with the nodes of the cones of traction and compression:
·
when the elastic prediction checks the condition of projection a priori in compression and
that the nodes of the cones of traction and compression were inverted on the axis
hydrostatic, and that projection at the top of the cone of traction did not give a solution
validate,
one carries out a projection at the top of the cone of compression:
·
when the elastic prediction checks the condition of projection a priori in compression and
that the nodes of the cones of traction and compression were inverted on the axis
hydrostatic, and that projection at the top of the cone of traction did not give a solution
validate, and that simultaneous projection with the nodes of the two cones did not give
valid solution,
in medium of algorithm,
one carries out one, two or three standards resolutions with projection on the criterion of compression or
on the criterion of traction or the two criteria at the same time,
and at the end of the algorithm,
when that the standards resolutions with activation of a criterion (traction or compression) or
two criteria at the same time did not give a solution,
one carries out a projection at the top of the cone of traction:
·
when the elastic prediction checks the condition of projection a posteriori in traction,
·
when the elastic prediction checks the condition of projection a posteriori in compression
(and that the nodes of the cones of traction and compression were inverted on the axis
hydrostatic,
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
51/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
one carries out a simultaneous projection with the nodes of the cones of traction and compression:
·
when the elastic prediction checks the condition of projection a posteriori in compression
and that the nodes of the cones of traction and compression were inverted on the axis
hydrostatic, and that projection at the top of the cone of traction did not give a solution
validate,
one carries out a projection at the top of the cone of compression:
·
when the elastic prediction checks the condition of projection a posteriori in compression
and that the nodes of the cones of traction and compression were inverted on the axis
hydrostatic, and that projection at the top of the cone of traction did not give a solution
validate, and that simultaneous projection with the nodes of the two cones did not give
valid solution.
With the exit of each resolution having converged, one carries out the checks of conformity of the solution
following:
·
validity of the solution compared to the second criterion, when the resolution was made with one
only criterion. In all the cases, it is enough to check that the two computed criterions with
final stress, are negative or null,
·
conformity of the solution: one calculates in the course of resolution the final equivalent stress. It
arrive sometimes that the solution is beyond the node of the cone which is hammer-hardened, which leads to
an equivalent stress “negative”. Numerically, that results in a final criterion
strictly positive. To check the conformity of the solution, it is enough to check that both
computed criterions with the final stress, are negative or null,
·
validity of projections at the tops of the cones. It should be checked that after resolution, when one
knows the work hardening of the criterion, the slope of the straight line connecting the elastic prediction to
projection is lower than the slope of the direction of projection. (condition of projection has
posteriori at the node of the cones). In the contrary case, that means that there is a solution
with standard resolution.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
52/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A




background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
53/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A

















background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
54/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A

















































Note:
In the case of projections at the tops of the cones, one starts systematically with
projection at the top of the cone of traction. If this solution is valid, the aforementioned is preserved. In
the contrary case, if the criterion of traction is activated, one carries out a resolution with projection with
nodes of the two cones. If the new solution is valid, that one is preserved. If not, one
carry out a resolution with projection at the top of the cone of compression alone.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
55/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Note:
If the conditions of projection at the tops of the cones are activated, one of the three
solutions must be valid, but for particularly important elastic jumps, it may be
that the resolution does not succeed. The solution is then to carry out a recutting of the pitch of
time.

Note:
Projection at the tops of the two cones simultaneously supposes that the criterion of traction is
activated. It may be very well that in the event of permutation of the nodes, only the criterion of compression
finds itself activated. It is necessary to then make a projection at the top of the cone of compression alone.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
56/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
Appendix 1 snap-back with the initial values of the coefficients C
and D

We show in this appendix the problem of snap-back met in the simulation of a tensile test
follow-up compression simple, if the choice of the coefficients
C
and
D
criterion of traction corresponds to
a situation where the criterion of traction cuts the axes in a diagram of stress plane, i.e. it
choice of the coefficients [éq 3.3-1] and [éq 3.3-2] leading to a field of reversibility represented on
[Figure 3.3-b]:
The assumptions are as follows:
·
one takes into account only the criterion of traction,
·
the marrow of work hardening is of the type:
()
p
T
T
p
T
T
H
F
F
+
=
,
·
one notes simply:
p
T
=
so that the curve of work hardening is written:
+
=
H
F
F
T
T
,
·
the null Poisson's ratio,
·
work hardening is negative,
·
the condition of applicability is met:
0
<
<
-
H
E
.

One considers an axial plain test controlled in deformation according to X, as indicated on Figure 5 -
xx
0
xx
1
xx
2
xx
Time T
xx
0
xx
1
xx
2
xx
xx
0
P
1
P
2
P
Appear 5-a

Under the other the imposed directions y and Z, conditions are conditions of null stresses:
0
=
=
zz
yy
.
One starts by imposing a deformation of traction
1
xx
such that there is a plasticization in traction, but
without the limit of traction falling down to 0. It is the point
1
P
in the diagram stress-strain. One notes
0
xx
the deformation for which the ultimate strength in traction appears for the first time. Beyond the point
1
P
, one imposes a decrease of the deformation, which involves an elastic unloading of material,
up to the value
0
2
<
xx
deformation for which one has a plasticization again, but this time under
compressive stress. The object of this appendix is primarily to study the behavior of the model
not retained (that corresponding to the formulas [éq 3.3-1] and [éq 3.3-2]) beyond the point
2
P
.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
57/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
A1.1 Calculation of the stresses and the deformations during the loading
Taking into account the assumptions pointed out higher, the invariants of stress are worth:
xx
I
=
1
éq A1.1-1
3
2
xx
J
=
éq A1.1-2
xx
eq
=
éq A1.1-3
The plastic flow is calculated by:
+
=
&
&
D
C
p
xx
3
2
éq A1.1-4
A1.1.1 Paths 0P
0
and P
0
P
1
Taking into account the relations [éq A1.1-1] with [éq A1.1-4], and the fact that along this way them
stresses
xx
are positive, one finds easily
(
)








+
+
+
+
=
2
2
'
9
2
3
2
1
1
D
E
C
H
F
D
C
H
E
T
xx
xx
éq A1.1.1-1
A1.1.2 Path P1P2
By definition, P1P2 is an elastic path of discharge, the P2 point being such as the criterion is there with
new reached
1
2
2
3
xx
xx
C
D
-
=
éq A1.1.2-1
1
1
2
2
3
xx
xx
xx
E
C
D
+


-
=
éq
A1.1.2-2
A1.1.3 Beyond the P2 point
One is interested now in the slope of the curve at the P2 point in the reference mark
(
)
xx
xx
,
. More
precisely one is interested in the slope of this curve for a dissipative solution.
By writing that the material remains plastic beyond the P2 point, i.e. that the state of stress
remain on criterion-which is hammer-hardened, one finds:
-
=
&
&
H
C
D
xx
2
3
éq A1.1.3-1
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
58/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
In addition, by calculating the increase in flow plastic, and by deferring it in the calculation of
the increase in stress, one finds:
(
)
2
3
-
+
=
C
D
E
xx
xx
&
&
&
éq A1.1.3-2
One can then eliminate
&
between [éq A1.1.3-1] and [éq A1.1.3-2] and one obtains:
(
)
2
2
2
9
2
T
xx
xx
E
D
C
E
H
He
=
-
+
=
&
&
éq
A1.1.3-3
This formula gives the slope of the response in the plan
(
)
xx
xx
,
The numerator is always negative since H is negative.
The sign of
2
2
T
xx
xx
E
=
&
&
depends on the sign of the denominator, thus two cases are posed:
If
(
)
2
2
9
2
D
C
E
H
-
-
<
then
2
2
T
xx
xx
E
=
&
&
one is positive has a configuration of snap-back.

If
(
)
2
2
9
2
D
C
E
H
-
-
>
then
2
2
T
xx
xx
E
=
&
&
is negative and there is no snap-back.
A new condition appears to avoid the snap back, condition which we already compare with that
evoked, but which related to in fact possible the snap back at the P1 point.
Slope at the P1 point in the reference mark
(
)
xx
xx
,
:
E
H
He
E
T
+
=
1
,
Slope at the P2 point in the reference mark
(
)
xx
xx
,
:
(
)
2
2
9
2
2
D
C
E
H
He
E
T
-
+
=
.
However
(
)
2
'
'
2
2
2
3
3
9
2




=
=
-
C
T
F
F
D
C
.
One can for example express
2
T
E
according to
1
T
E
while eliminating
H
:
(
)
2
2
3
1
3
1
1
2
-
+
=
T
T
T
E
E
E
E
E
.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
59/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
For example,
=
2
T
E
for
)
3
1
(
3
2
2
1
-
-
=
E
E
T
.
As example, for E=32000 Mpa,
Mpa
F
T
3
'
=
and
Mpa
F
C
3
,
38
'
=
one finds
601
1
-
T
E
. Thus
1
T
E
is very weak compared to E. as illustrated on [A1.1.3-a Figure].
xx
xx
Appear A1.1.3-a
Thus, a condition implying that there is no snap back at the P2 point would be too restrictive and
would lead to practically choose a material not fragile in traction.
For this reason we preferred to modify the expression of the coefficients C and D like
indicated at the paragraph [§3.3].
This said, and even if the adopted solution, consisting in modifying the coefficients
C
and
D
seem
reasonable, the example treated in this appendix shows that a very simple problem can finally be
a problem of structure: there is in this example of the equilibrium conditions, they are the conditions:
0
=
=
zz
yy
.
background image
Code_Aster
®
Version
5.0
Titrate:
Law of behavior doubles Drücker Prager
Date
:
19/03/02
Author (S):
C. CHAVANT, B. CIREE
Key
:
R7.01.03-A
Page
:
60/60
Manual of Reference
R7.01 booklet: Modeling for the Civil Engineering
HT-66/02/004/A
6 Bibliography
[1]
J.F. GEORGIN “Contribution to the modeling of the concrete under stress of dynamics
rapid. The taking into account of the effect speed by viscoplasticity “- Thesis (15/01/98).
[2]
G. HEINFLING “Contribution to the numerical modeling of the behavior of the concrete and of
concrete structures reinforced under mechanical thermo stresses at high temperature “-
Thesis (14/01/98)
[3]
R.T. ROCKAFELLAR “Convex analysis” Princeton University Close 1972
[4]
B. HALPHEN and NGUYEN QUOC HIS, “On generalized standard materials” Newspaper
of Mechanics Flight 14 n° 1, 1975
[5]
NR. TARDIEU, I. VAUTIER, E. LORENTZ “quasi static nonlinear Algorithm” document
of reference Aster [R5.03.01].