background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
1/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Organization (S):
EDF-R & D/AMA















Manual of Reference
R5.03 booklet: Nonlinear mechanics
Document: R5.03.04



Relations of behavior élasto-visco-plastic
of Chaboche




Summary:

This document describes the integration of the model of behavior élasto-visco-plastic of J.L. Chaboche with
nonlinear and isotropic kinematic work hardening. Two versions of the model are available in
Code_Aster:
·
a version with one or two variables kinematics, introduced recently, takes into account all them
variations of the coefficients with the temperature, and has an effect of work hardening on the variables
tensorial of recall. This version also makes it possible to modelize (in an optional way) the character
viscous of the material (viscosity of Norton). It is integrated by the solution of only one equation
nonlinear scalar.
This model is available in 3D, plane deformation, axisymetry. Modeling in plane stress
use a method of condensation static (of Borst).
·
a version with two variables kinematics which exists in Code_Aster since the version 2, which
does not take into account all the variations of the coefficients compared to the temperature, but which was
used for several studies, and for which one has sets of identified parameters. This
version is integrated in environment PLASTI. It does not make it possible to modelize viscosity.
It does not take into account of effect of work hardening on the tensorial variable of recall.
This model is available in 3D, plane strain, plane stress and axisymetry.
One gives also elements to identify the coefficients of the relation of behavior.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
2/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
1 Models élasto-visco-plastics of J.L. Chaboche
available in Code_Aster
For the structural analysis subjected to cyclic loadings, work hardenings isotropic (linear
or not) and linear kinematics conventional [R5.03.02] and [R5.03.16] are not sufficient any more. In
private individual, one cannot correctly describe the stabilized cycles obtained in experiments on
a tensile specimen subjected to an alternated imposed deformation or a traction and compression.
If one seeks to precisely describe the effects of a cyclic loading, it is desirable to adopt
modelings more sophisticated (but easy to use) such as the model of Said Taheri, by
example, cf [R5.03.05], or if the number of cycles is limited the model of Jean-Louis Chaboche
who is introduced here.
Actually, the model of J.L. Chaboche can be more or less sophisticated. Developed models
in Code_Aster comprise is a kinematic variable (
VISC_CIN1_CHAB
) is two
(
VISC_CIN2_CHAB
and
CHABOCHE),
and of isotropic work hardening.
The choice to use two variables kinematics complicates certainly the model, but makes it possible to identify
correctly uniaxial tests in a broader range of deformations [bib2], [bib7]. One
certain number of indentifications of the parameters of this model were carried out mainly for
the stainless steels A316 and A304 ([bib7], [bib8]).
The models comprise 8 parameters (a kinematic variable) or 10 (two variables
kinematics), introduced into the control
DEFI_MATERIAU
:
CIN1_CHAB (CIN1_CHAB_FO)
= _F (
R_0 =
R_0,
R_I
=
R_I, (useless
if
B=0)
B
=
B,
(defect: 0.)
C_I =
C_I,
K
=
K,
(defect: 1.)
W
=
W,
(defect: 0.)
G_0
=
G_0,
A_I
=
A_I, (defect
:
0.)
)
CIN2_CHAB (CIN2_CHAB_FO)
= _F (
R_0 =
R_0,
R_I
=
R_I, (useless
if
B=0)
B
=
B,
(defect: 0.)
C1_I = C1_I,
C2_I = C2_I,
K
=
K, (defect: 1.)
W
=
W, (defect: 0.)
G1_0 = G1_0,
G2_0 = G2_0,
A_I
=
A_I
,
(defect
:
0.)
)
The 8 or 10 parameters are real constants. All these parameters can depend on
temperature (key words
CIN1_CHAB_FO
or
CIN2_CHAB_FO
) and the awaited values are of type
function
.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
3/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
If one wants to introduce in addition to viscosity (models
VISC_CIN1_CHAB
and
VISC_CIN2_CHAB
), it is also necessary to provide in the control
DEFI_MATERIAU
, under the key word
LEMAITRE
(or
LEMAITRE_FO
) parameters
NR
and
UN_SUR_K
, which can depend on
temperature.
The parameter
UN_SUR_M
key word
LEMAITRE
(respectively
LEMAITRE_FO
) must obligatorily
to be put at zero (respectively with the identically null function). In the absence of one of the key words
LEMAITRE
or
LEMATIRE_FO
, the behavior is supposed plastic.
Model CHABOCHE is a model with two variables kinematics with isotropic work hardening, but
without the effect of work hardening on the term of recall and without taking into account of the variation of C
1
and C
2
with the temperature. The characteristics of work hardening are given by 9 constants, introduced
in the control
DEFI_MATERIAU
:
CHABOCHE = _F
(
R_0
=
R_0,
R_I
=
R_I,
B = B,
K = K,
W = W,
A1
=
A1,
A2
=
A2,
C1
=
C1,
C2
=
C2,
)
In this case the characteristics do not depend any more a temperature since version 5, because these
variations were badly taken into account by this model.
The use of these laws of behavior is accessible in the controls
STAT_NON_LINE
or
DYNA_NON_LINE
by the key words
VISC_CIN1_CHAB
,
VISC_CIN2_CHAB
or
CHABOCHE
of
COMP_INCR
.
In the continuation of this document, one describes the models precisely
VISC_CIN1_CHAB
and
VISC_CIN2_CHAB
. One presents then the detail of his numerical integration in link with
construction of the coherent tangent matrix. One also briefly describes the integration of the model with
two variables kinematics CHABOCHE. Lastly, one also gives some elements for
identification of the characteristics of material.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
4/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
2 Models
VISC_CIN1_CHAB
and
VISC_CIN2_CHAB
2.1
Description of the models
At any moment, the state of material is described by the deformation
, the temperature
T
, deformation
plastic
p
, cumulated plastic deformation
p
and the tensor of recall
X
. Equations of state
then define according to these variables of state the stress
=
H
Id
+
~
(broken up into
parts hydrostatic and deviatoric), the isotropic share of work hardening
R
and the kinematic share
X
:
()
(
)
(
)
H
ref.
tr
K tr
T-T
=
=
-
1
3
HT
HT
Id
with
=
éq 2.1-1
(
)
~
~
= -
=
-
µ
H
Id
p
2
éq 2.1-2
()
R
p
=
R
éq 2.1-3
(
)
(
)
(
)
X X
X
X
=
=
+
p
p
p
p
p
p
,
,
,
1
2
éq
2.1-4
where
K,
µ
and coefficients of
()
()
X p
R p
and
are characteristics of the material which can
to depend on the temperature. More precisely, they are respectively the modules of compressibility
and of shearing, the thermal expansion factor, functions of isotropic work hardening and
kinematics. As for
T
ref.
, it is about the temperature of reference, for which one considers
thermal deformation as being null.
Note:
For the model
VISC_CIN1_CHAB
only the only tensorial variable is considered
()
X
1
p
thus
()
X
2
0
p
=
. This remains valid for all the continuation: one will describe the two models formally of
the same way, the model
VISC_CIN1_CHAB
resulting from
VISC_CIN2_CHAB
while supposing
()
X
2
0
p
=
.
The evolution of the plastic deformation is controlled by a normal law of flow to a criterion of
plasticity of von Mises:
(
)
(
)
()
F,
~
~: ~
R
R p
eq
eq
X
X
X
With
WITH A
=
-
-
-
=
1
2
3
2
with
éq
2.1-5
(
)
&
& F
& ~
~
p
X
X
X
X
=
=
-
-
-
-
3
2
1
2
1
2 eq
éq 2.1-6
&
&
&
&
p
= =
2
3
p
p
:
éq 2.1-7
As for the plastic multiplier
&
, it is obtained by the condition of coherence:
if
or
if
and
F
&F
&
F
&F
&
<0
0
0
0
0
0
<
=
=
=



éq 2.1-8
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
5/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Note:
Evolution of the variables
X
1
and
X
2
is given by:
()
()
()
()
X
X
p
p
1
1
2
2
1
2
2
3
2
3
,
,
&
&
&
&
&
&
=
=
=
-
=
-
C p
C p
p
p
p
p
1
2
1
1
2
2
éq 2.1-9
Functions
() ()
C p
p
and
()
R p
are defined, in accordance with [bib2] by:
()
(
)
()
(
)
(
)
()
(
)
(
)
R p
R
R
R E
C p
C
K
E
C p
C
K
E
LP
wp
wp
=
+
-
=
+
-
=
+
-
-
-
-
0
1
1
2
2
1
1
1
1
()
(
)
(
)
()
(
)
(
)
1
1
0
2
2
0
1
1
p
has
E has
p
has
E has
LP
LP
=
+ -
=
+ -
-
-
The presence of viscosity can modelize in a simple way (cf Lemaitre and Chaboche [bib2]) in
replacing the condition of coherence [éq 2.1-8] by:
NR
K
F


=
&
éq 2.1-10
F
positive part of
F
(hooks of Macauley)
NR
K,
characteristics of viscosity (Norton) of material
Unchanged all the other equations of the model are left. It will be seen that such an introduction of
viscosity involves only minor amendments of the implicit algorithm of integration of the law of
behavior.

Note:
The definition of
X
1
and
X
2
in the form [éq 2.1-9]:
·
allows to keep a formulation which takes into account the variations of the parameters with
temperature without introducing term in
&T
as in [bib.4], in the same way that it
viscoplastic model of Chaboche. These terms necessary because their are not taken in
account would lead to inaccurate results [bib4].
·
allows to have a coherent writing with the thermodynamic expression of the potential
plastic [bib2] (p.221).
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
6/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Note:
Significance of the functions
() () () () ()
C p
PC p
p R p
1
1
2
2
,
,
,
,
:
it is noted that the functions
() () () () ()
C p
PC p
p R p
1
1
2
2
,
,
,
,
intervening in the equations
the preceding ones allow all the three to modelize various nonlinear effects of work hardening.
The introduction of work hardening, is on the level of the kinematic part, by
()
C p
, that is to say on the level
term of recall, by the function
()
p
, the same effect on the classification tests does not have
[bib2]. The use of a model with
()
p
allows in particular to identify more easily of
strong cyclic work hardenings. Several work of identification of the coefficients of the models of
Chaboche were carried out besides on the basis of model with a work hardening represented by
()
p
([bib5], [bib6]), in particular for stainless steel A316L.
2.2
Integration of the relations
VMIS_CIN1_CHAB
and
VMIS_CIN2_CHAB
To numerically carry out the integration of the law of behavior, one carries out a discretization in
time and one adopts a diagram of implicit, famous Euler adapted for relations of behavior
elastoplastic. Henceforth, the following notations will be employed:
With
With
With
-
,
and
represent
respectively values of a quantity at the beginning and the end of the pitch of time considered thus that sound
increment during the pitch. The problem is then the following: knowing the state at time
T
-
like
increments of deformation
(resulting from the phase of prediction (cf Doc. R
STAT_NON_LINE
[R5.03.01])) and of temperature
T
, to determine the state of the variables intern at time
T
as well as
stresses
.
One takes into account the variations of the characteristics compared to the temperature while noticing
that:
(
)
H
H
K
K
K
=
+
-
-
-
tr
HT
éq
2.2-1
(
)
~
~
~
~
-
=
+
=
-
-
-
µ
µ
µ
µ
2
2
p
p
E
éq
2.2-2
with
~
~
~
E
=
+
-
-
µ
µ
µ
2
Within sight of the equation [éq 2.2-1], one notes that the hydrostatic behavior is purely elastic
if
K
is constant. Only the processing of the deviatoric component is delicate.
In the absence of viscous term, the relation of discretized coherence is:
Elastic mode:
0
0
=
p
F
and
Plastic mode:
0
0
=
p
F
and
On the other hand, in the presence of viscosity, the condition of coherence is replaced by the equation
[éq 2.1-10] which, discretized, is written:
NR
NR
T
p
K
F
K
F
T
p
/
1


=




=
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
7/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
In other words, while posing:
NR
T
p
K
F
F
/
1
~


-
=
the viscoplastic increment of cumulated deformation is determined by:
0
0
~
0
0
~
=
=
p
F
p
F
and
:
ic
viscoplast
Mode
and
:
rubber band
Mode
éq 2.2-3
Finally, by adopting an implicit discretization, the only difference between the laws of behavior
plastic and viscoplastic lies in the form of the function of load
F
: a term there is observed
complementary in the event of viscosity. In fact, incremental plasticity seems the borderline case
incremental viscoplasticity when
K
tends towards zero. This convergence was already described by
J.L. Chaboche and G. Cailletaud in [bib3].
In the continuation of this paragraph, one will thus detail the integration of the viscoplastic law. To find
the case of the plastic behavior, it is enough to take
0
=
K
in the equations below (one
recall that the user to place itself in this case must obligatorily remove the key word
LEMAITRE
or
LEMAITRE_FO
control
DEFI_MATERIAU
).
(
)
~
~
+
1
2
-
-
=
-
-
-
-
-
-
X
X
p
1
2
1
1
2
2
1
2
2
3
2
3
2
2
3
E
C
C
C
C
µ
Equations of flow [éq 2.1-6] and [éq 2.1-7], once discretized, and the condition of coherence
[éq 2.2-3] are written (by noticing that
p
=
):
1
2
-
1
2
1
2
-
1
2
p
E
p
E
p
=
-
-
-
-
-
-
-
-
-
-




3
2
2
3
2
3
2
2
3
2
3
2
3
2
3
2
2
3
2
3
1
2
1
2
1
2
1
2
p
C
C
C
C
C
C
C
C
eq
~
~
µ
µ
éq 2.2-4
0
F
~
0
0
F
~
=
p
p
éq
2.2-5
The processing of the condition of coherence (preceding equation) is conventional. One starts with one
elastic test (
p
=
0
) which is well the solution if the criterion of plasticity is not exceeded, it be-with
to say if:
()
()
()
~
-
-
E
eq
C p
C p
R p
-
-




-
<
-
-
-
2
3
2
3
0
1
1
2
2
éq
2.2-6
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
8/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
In the contrary case, the solution is plastic (
p
>
0
) and the condition of coherence is reduced to
0
F
~
=
. To solve it, it is shown that one can bring back oneself to a scalar problem while expressing
p
and
1
2
,
according to
p
. By gathering the equations of the problem resulting from the discretization
implicit, the system of equations is obtained:
()
NR
eq
T
p
K
p
R
C
C
C
C
/
1
2
1
2
1
3
2
3
2
2
3
2
3
2
~


+
=


-
-
-
-
-
-
-
2
1
2
1
p
E
µ
éq 2.2-7
()
NR
T
p
K
p
R
C
C
C
C
p
/
1
2
1
2
1
3
2
3
2
2
3
2
3
2
~
2
3


+
-
-
-
-
-
=
-
-
2
1
2
1
p
p
E
µ
éq
2.2-8
1
1
2
2
=
-
=
-
p
p
1
2
p
p
éq
2.2-9
In this writing, it should well be noted that
p
p
p
=
+
-
and
I
I
I
=
+
-
and that
C
I
I
,
are
functions of p. By considering the three last equations, this linear system in
p
and
I
can
to be solved to express these quantities according to
p
. Indeed, it is equivalent to:


-
-
-
-
=






+
+
-
-
2
1
2
1
2
1
2
1
/
1
~
2
3
3
)
(
C
C
C
C
p
T
p
K
p
p
R
NR
E
p
µ
éq
2.2-10
(
)
(
)
1
2
1
1
1
1
1
2
2
2
+
=
-
+
=
-
-
-
p
p
p
p
p
p
éq
2.2-11
While calculating
C
1
1
and
C
2
2
and by replacing them in the expression of
p
one is obtained
expression of
p
according to
p
only:
()
()
()
()
()
()
()
C
C
p
C
p
p
M p
M p
p
C
C
p
C
p
p
M p
M p
p
M p
C p
p p
I
I
I
1
1
1
1
1
1
1
1
1
2
2
2
2
2
2
2
2
2
1
1
1
1
1
1
1
1
2
2
-
2
-
=
+




-
+






=
-
=
+




-
+






=
-
= +
p
p
p
p
with
éq 2.2-12
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
9/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
By deferring this expression in the expression of
p
one finds:
(
)
(
)
(
)
(
)


-
+
-
-






+
+
+
+
=
-
-
2
1
p
M
C
p
M
C
p
p
T
p
K
p
M
M
p
R
NR
2
2
2
1
1
1
/
1
2
1
~
2
3
3
)
(
1
µ
E
p
what is simplified in:
(
)
1
2
-
p
E
=
-
+




1
3
2
1
2
D p
p
p M
M
()
~
éq
2.2-13
with:
() ()
()
()
(
)
NR
T
p
K
p
p
M
p
M
p
R
p
D
/
1
2
1
3


+
+
+
+
=
µ
It now only remains to replace
p
in the expressions of
C
1
1
and
C
2
2
for
to express this term according to
p
by:
(
)
(
)
C
M
D
p
p M
M
M
p
C
M
D
p
p M
M
M
p
1
1
1
1
2
1 1
2
2
2
1
2
2 2
3
2
3
2
1
2
-
1
1
2
-
2
-
=
-
+




-
=
-
+




-
~
~
E
E
then to substitute the expression obtained thus that
p
according to
p
in the equation
0
F
~
=
, and
one obtains a scalar equation in
p
to solve, namely:
()
()
0
3
2
3
2
2
3
2
3
2
~
~
/
1
2
1
2
1
=


-
-


-
-
-
-
-
=
-
-
NR
eq
T
p
K
p
R
C
C
C
C
p
F
2
1
2
1
p
E
µ
what is simplified in:
()
()
()
()
0
3
2
3
2
~
~
/
1
2
1
/
1
=


-
-


-
-


+
=
-
-
NR
eq
NR
T
p
K
p
R
M
M
p
D
T
p
K
p
R
p
F
2
1
E
éq 2.2-14
This scalar equation in
p
is solved numerically, by a method of search for zero of
function (method of secants which one briefly describes in appendix 2).
Once determined
p
, one can calculate
p
using the equation [éq 2.2-13] then
1
and
2
with
the aid of the equations [éq 2.2-11]. It any more but does not remain to calculate the tensor of the stresses, by
equations [éq 2.2-1] and [éq 2.2-2], and to bring up to date the internal variables
1
and
2
.
Note:
·
an interesting borderline case (for the validation of this model) arises while posing
I
=
0
. One
finds itself then exactly in the situation of linear kinematic work hardening (if
R p
y
()
=
, [R5.03.02]) or of mixed work hardening for
R p
()
unspecified (cf [R5.03.16]),
·
these models are also available in plane stresses, by a total method
(static condensation due to R. of Borst) [R5.03.03].
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
10/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
2.3
Calculation of tangent rigidity
In order to allow a resolution of the total problem (equilibrium equations) by a method of
Newton [R5.03.01], it is necessary to determine the coherent tangent matrix of the problem
incremental.
This matrix is composed classically of an elastic contribution and a plastic contribution:
µ
=
-
E
p
2
éq
2.3-1
with
E
= +
2
µ
p
, which gives again in particular
~
~
~
E
=
+
-
-
µ
µ
µ
2
One immediately deduces from it that in elastic mode (conventional or pseudo-discharge), the matrix
tangent is reduced to the elastic matrix:
=
E
éq 2.3-2
For that, one once more adopts the convention of writing of the symmetrical tensors of command 2 pennies
form vectors with 6 components. Thus, for a tensor
has
:
[
]
has
=
T
xx
yy
zz
xy
xz
yz
has
has
has
has
has
has
2
2
2
éq
2.3-3
If moreover the hydrostatic vector is introduced
1
and stamps it deviatoric projection
P
:
[
]
1
=
T
1 1 1 0 0 0
éq 2.3-4
P
Id
1 1
=
-
1
3
éq 2.3-5
where
is the tensorial product
Then the matrix of coherent tangent rigidity is written for an elastic behavior:
µ
E
1 1
P
=
+
K
2
éq 2.3-6
On the other hand, in plastic mode, the variation of the plastic deformation is not null any more.
One derives compared to
~
E
, knowing that one a:
µ
p
p
E
E
p
E
P
=
=
~.
~
~.
2
éq
2.3-7
S
space symmetrical tensors
P
projector on the diverters
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
11/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
To calculate
p
E
~
, one uses the expression of
p
according to
~
E
and p:
(
)
1
2
-
p
E
=
-
+




1
3
2
1
2
D p
p
p M
M
()
~
what is written in the form:
()
()
()
1
2
-
p
E
=
+
+
With p
B p
B p
~
1
2
Thus:
()
()
()
()
1
2
-
p
E
E
E
E
E
Id
~
~
~
~
~
=
+
+
+
With p
With p
B p
B p
1
2
Quantities of the type
()
With p
~
E
are calculated using:
()
()
With p
With p
p
p
~
~
E
E
=
Finally, it any more but does not remain to calculate the variation of
p
:
p
~
E
One uses for that:
()
0
~
,
F
~
=
E
p
(
)
()
()
()
0
3
2
3
2
~
~
,
~
/
1
2
1
/
1
=


-
-


-
-


+
=
-
-
NR
eq
NR
T
p
K
p
R
M
M
p
D
T
p
K
p
R
p
F
2
1
E
E
()
()
()
()
E
E
E
E
E
E
E
E
~
,
F
~
~
,
F
~
~
~
~
,
F
~
~
,
F
~
,
~
,
~
,
,
p
p
p
p
p
p
p
p
-
=
-
=
éq
2.3-8
The detail of calculations is given in appendix 1.
The initial tangent matrix, used by the option
RIGI_MECA_TANG
is obtained by adopting it
behavior of the preceding pitch (elastic or plastic, meant by internal variable
being worth 0 or
1) and while making tend
p
towards zero in the preceding equations.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
12/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
2.4
Significance of the internal variables
The variables intern of the two models at the points of Gauss (
VELGA
) are:
·
V1
= p: cumulated plastic deformation (positive or null)
·
V2
=
: being worth 1 if the point of Gauss plasticized during the increment or 0 if not.
The following internal variables are, for modeling
3D:
·
For the model
VISC_CIN1_CHAB
-
V3 =
1xx
-
V4 =
1
yy
-
V5 =
1
zz
-
V6 =
1
xy
-
V7 =
1
xz
-
V8 =
1
yz
·
For the model
VISC_CIN2_CHAB
-
V3 =
1
xx
-
V4 =
1
yy
-
V5 =
1
zz
-
V6 =
1
xy
-
V7 =
1
xz
-
V8 =
1
yz
-
V9 =
2
xx
-
V10 =
2
yy
-
V11 =
2
zz
-
V12 =
2
xy
-
V13 =
2
xz
-
V14 =
2
yz
For modelings
C_PLAN
,
D_PLAN
, and
AXIS
:
-
V7 =
0
-
V8 =
0
-
V13 =
0
-
V14 =
0
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
13/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
3
Model with two variables kinematics:
CHABOCHE
3.1
Description of the model
This model comprises two tensorial variables which describe the kinematic share of work hardening:
X
X
X
=
+
1
2
. The equations of the behavior are then:
()
(
)
(
)
H
ref.
tr
K tr
T-T
=
=
-
1
3
HT
HT
Id
with
=
éq 3.1-1
(
)
~
~
= -
=
-
µ
H
Id
p
2
éq 3.1-2
()
R
p
=
R
éq 3.1-3
() ()
X
X
X
=
+
1
2
p
p
p
p
,
,
éq 3.1.4
where
K,
µ
and coefficients of the functions
() ()
R p
p
,
X
1
and
()
X
2
p
are characteristics of
material which can depend on the temperature. More precisely, they are respectively them
modules of compressibility and shearing, the thermal expansion factor, functions
of isotropic and kinematic work hardening. As for
T
ref.
, it is about the temperature of reference, for
which one regards the thermal deformation as being null.
The evolution of the internal variables is controlled by a normal law of flow to a criterion of
plasticity:
(
) (
)
()
F,
,
~
~: ~
R
R p
eq
eq
X X
X
X
With
WITH A
1
2
1
2
3
2
=
-
-
-
=
with
éq
3.1-5
(
)
&
& F
& ~
~
p
X
X
X
X
=
=
-
-
-
-
3
2
1
2
1
2 eq
éq 3.1-6
&
&
&: &
p
= =
2
3
p
p
éq 3.1-7
As for the plastic multiplier
&
, it is obtained by the condition of coherence:
if
or
if
and
F
&F
&
F
&F
&
<0
0
0
0
0
0
<
=
=
=



éq 3.1-8
Evolution of the variables
X
I
is given in the model
CHABOCHE
by:
()
&
&
&,
,
X
X
I
I
I
p
I
C
has
p
p
I
=
-




=
2
3
1 2
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
14/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Functions
()
p
and
()
R p
are defined by:
()
(
)
(
)
p
K
E
wp
= +
-
-
1
1
()
(
)
R p
R
R
R E
LP
=
+
-
-
0
3.2
Integration of the relation of behavior
CHABOCHE
As for the relations of behavior
VISC_CIN1_CHAB
and
VISC_CIN2_CHAB
, one is adopted
diagram of implicit Euler. One takes into account the variations of the elastic characteristics of
even way that previously [éq.2.1-1], [éq.2.1-2].
One starts with an elastic test, while taking for stress field:
~
~
~
E
=
+
-
-
µ
µ
µ
2
who is well the solution of the problem if:
()
()
()
~
-
-
E
C p
C p
R p
-
-
-
<
-
-
-
2
3
2
3
0
1
1
2
2
In the contrary case, the solution is elastoplastic. It is then necessary to solve the system of equations not
linear according to:
(
)
~
~
-
-
-
µ
µ
µ
-
-




-
=
1
2
0
p
()
()
X
X
X
X
p
p
1
1
2
2
-
+
=
-
+
=
2
3
0
2
3
0
1 1
1
2 2
2
A.c.
p
C
p
A.c.
p
C
p
(
)
()
~
-
-
1
2
X
X
eq
R p
-
=
0
that one can write in a way more contracted in the following form:
()
()
()
()
()
F
y
G
y
H
y
I
y
J
y
y
p
L
= =








=






0
1
2
with
X
X
One solves this system by the method of Newton proposed in environment PLASTI, (described in
detail in [R5.03.10]), is:
()
()
()
F
y
D
y
F
y
y
y
D
y
L
K
K
L
K
K
K
K
= -
=
+


+
1
While reiterating in
K
until convergence. This resolution takes place for each point of integration.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
15/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
The resolution requires the calculation of Jacobien of the local system
F
L
. His general expression is given
hereafter; analytical calculations are not detailed in this document.
J
G
G
G
G
p
H
H
H
H
p
I
I
I
I
p
J
J
J
J
p
=


















X
X
X
X
X
X
X
X
1
2
1
2
1
2
1
2
3.3
Operator of tangent behavior
After resolution of the preceding discretized system, the solution obtained is such as
the equation
()
(
)
F
y
L
=
0
is checked at the end of the increment. One seeks to evaluate the tangent operator in
this point, i.e.




+
T
T
.
For a small variation of
F
L
, by considering this time
like variable and not like
parameter, the system remains with balance and one checks
dF
L
=
0
, i.e.:


F
F
F
F
F
p
p
L
L
L
L
L
+
+
+
+
=
X
X
X
X
1
1
2
2
0
This system can be still written:
()
F
y
y
X
X
H
L
=
=










,
with
0
0
0
0
This system of equations can be put in the form:
K
H
=
from where the required tangent operator:




=
+
-
T
T
K H
1
One is led to re-use the same matrix jacobienne
J
that previously to evaluate the operator
tangent. The calculation of
K
-
1
is carried out numerically by a method of decomposition of Gauss.
Its expression is detailed in appendix 1.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
16/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
3.4
Significance of the internal variables
The variables intern model at the points of Gauss (
VARI_ELGA
) are, for modeling
3D
:
·
V1 =
X
xx
1
·
V2 =
X
yy
1
·
V3 =
X
zz
1
·
V4 =
X
xy
1
·
V5 =
X
xz
1
·
V6 =
X
yz
1
·
V7 =
X
xx
2
·
V8 =
X
yy
2
·
V9 =
X
zz
2
·
V10 =
X
xy
2
·
V11 =
X
xz
2
·
V12 =
X
yz
2
·
V13
= p: cumulated plastic deformation (positive or null)
For modelings D_PLAN, and AXIS,
·
V5 =
0
·
V6 =
0
·
V11 =
0
·
V12 =
0

4
Comparison of the models
VMIS_CIN2_CHAB
and
CHABOCHE
The main difference between the models
VISC_CIN2_CHAB
and
CHABOCHE
relate to the evolution of
variables kinematics.
In the case
VISC_CIN2_CHAB
, one a:
(
)
()
()
&
&
&
&
&
&
X
X
I
I
I
I
I
I
I
I
p
I
I
C
C
C
C p
p
p
=
+
=
+
-




2
3
2
3
2
3
In the case
CHABOCHE
, one a:
()
&
&
&,
,
X
X
I
I
I
p
I
C
has
p
p
I
=
-




=
2
3
1 2
These two models are not equivalent: in particular, variation of the coefficients with
temperature is not well taken into account in the model
CHABOCHE.
In the particular case where coefficients
C
I
are constant, the two models are equivalent. For
that, it is necessary to choose:
·
()
I
I
p
C
=
what implies:
I
I
A.c.
0
1
=
=
,
·
()
C p
C
A.c.
I
I
I I
=
=
and
W
=
0
, because of the additional term:
2
3 &
C
I
who allows
to take into account the variation of
C
I
with the temperature and the plastic deformation
cumulated p.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
17/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
5
Principle of the identification of the parameters of the model.
In the simplest case (only one kinematic variable,
1
1
=
=
cste C
cste
,
, ()
R p
y
=
) them
coefficients of the model
1
1
, C
can be identified on a simple tensile test uniaxial, or
on a cyclic curve of work hardening.
Indeed in the uniaxial case, the model is reduced in 1D to [bib2]:
dX
C D
D
X
p
p
y
1
1
1 1
1
1
,
=
-
= ±
-
=
X
that one can integrate (in monotonous loading) in the following way:
(
)
(
)
X
C
X
C
X
p
p
y
1
1
1
10
1
1
1
0
1
1
exp
,
=
+
-




-
-
= ±
=
+
whose asymptote of the traction diagram makes it possible to obtain
C
1
1
by:
p
X
C
1
1
1
thus
+




y
C
1
1
and whose slope in the beginning provides
(
)
C
X
1
10
0
If
=
:
p
X
C
y X
X
C
y X
-
=
-
0
1
1
1 10
10
1
1 1
&
For a model has two variables kinematics, without isotropic work hardening, a traction diagram
still allows to find these relations:

p
y
C
C
+
+








1
1
2
2
and the slope in the beginning is worth
C
C
1
2
+

But apart from these simple cases a numerical identification is necessary to obtain them
parameters. One will be able to make this identification for example on tensile tests compression with
imposed deformation.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
18/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
6 Bibliography
[1]
P. MIALON, Elements of analysis and numerical resolution of the relations of elastoplasticity.
EDF - Bulletin of the Management of the Studies and Search - Series C - N° 3 1986, p. 57 - 89.
[2]
J.LEMAITRE, J.L.CHABOCHE, Mechanics of solid materials. Dunod 1996
[3] J.L.CHABOCHE,
G.CAILLETAUD, constitutive Integration methods for complex equations,
Methods computer in Applied Machanics Engineering, N°133 (1996), p 125-155
[4]
J.L.CHABOCHE, Cyclic viscoplastic constitutive equations, Newspaper off Applied Mechanics,
Vol.60, December 1993, pp. 813-828
[5]
R. FORTUNIER, Law of behavior of Chaboche
: identification of the parameters
elastoplastic and élasto-visco-plastics of steel EDF-SPH between 20°C and 600°C. Note
FRAMATOME/Novatome, NOVTUDD90011, October 1990
[6]
C.MIGNE, Retiming of the parameters of the model of kinematic plasticity nonlinear of
SYSTUS. Modeling of the phenomenon of progressive deformation with consolidation
cyclic of material. Note FRAMATOME EE/R. DC.0286. September 1992.
[7]
J.J.ENGEL, G.ROUSSELIER, Behavior in uniaxial stress under loading
cyclic of the austenitic stainless steel 17-12 Mo with very low carbon and nitrogenizes control.
Identification of 20) C with 600°C of a model of elastoplastic behavior to work hardening
nonlinear kinematics. Note EDF/DER/EMMA N°D599 CHECHMATE/T43 (1985)
[8]
P. GEYER, C.COUTEROT, characterization of steel 304L used during the tests
“deformation, progressive” on CUMULUS and identification of the parameters of the model of
Chaboche, Note EDF/DER/HT-26/93/040/A
[9]
R of Borst “the zero normal stress condition in plane stress and Shell elastoplasticity”
Communications in applied numerical methods, Flight 7, 29-33 (1991)

background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
19/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Appendix 1 Stamps behavior tangent for the models
VMIS_CIN1_CHAB
and
VMIS_CIN2_CHAB
To obtain the tangent behavior in the elastoplastic case, it is necessary to calculate
p
E
~
[éq 2.3-7].
One uses for that the expression of
p
according to
~
E
and p, which is written in the form:
()
()
()
-
-
+
+
=
2
1
p
B
p
B
p
D
p
* 2
*
1
~
2
3
E
p
with
()
()
()
p
D
p
M
p
p
B
I
I
-
=
*
()
()
()
M p
C p
p p
I
I
I
= +
1
() ()
()
()
(
)
NR
T
p
K
p
p
M
p
M
p
R
p
D
/
1
2
1
3


+
+
+
+
=
µ
The following definitions are pointed out:
()
(
)
()
(
)
(
)
()
(
)
(
)
R p
R
R
R E
C p
C
K
E
p
has
has
E
LP
I
I
wp
I
I
LP
=
+
-
=
+
-
=
+ -
-
-
-
0
0
1
1
1
thus:
()
()
()
()
-
-
+
+




+
=
2
1
E
E
E
E
E
p
Id
~
~
~
~
2
3
2
3
~
* 2
*
1
p
B
p
B
p
D
p
p
D
p
Quantities of the type
()
With p
~
E
are calculated using:
()
()
With p
With p
p
p
~
~
E
E
=
These various terms are expressed by:
·
()
()
3
2
3
2
p
D p
p
I p




=
with
() () () ()
I p
D p
D p
D p
p
=
-
1
2
·
()
()
()
() ()
()
p
H
p
I
p
M
p
p
D
p
M
p
p
B
I
I
I
I
=
-
-
=
.
*
It remains to calculate:
p
~
E
One thus uses, following [éq 2.3-8]:
()
()
E
E
E
E
~
,
F
~
~
,
F
~
~
,
~
,
p
p
p
p
-
=
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
20/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
() ()
()
()
()
NR
NR
eq
T
p
K
p
R
p
G
T
p
K
p
R
p
S
p
/
1
/
1
~
,
~
,
~
,
F
~


-
-
=


-
-
=
E
E
E
with
()
()
() ()
()
p
D
T
p
K
p
R
p
M
B
p
D
T
p
K
p
R
With
B
B
With
NR
I
I
NR
E






+
-
=


+
=
+
+
=
-
-
/
1
/
1
2
2
1
1
3
2
~
S
Then, while posing
() ()
NR
v
T
p
K
p
R
p
R
/
1


+
=
:
()
()
()
()
()
()
(
)
()
()
()
()
()
()
()
()
()
()
()
() ()
()
()
() ()
()
()
()
()
(
)
()
1
1
2
'2
1
'1
3
2
22
1
21
2
2
2
1
3
2
22
1
21
1
'
,
2
2
1
1
'
,
'
,
~
,
~
:
2
3
1
1
*
~
.
2
3
:
2
3
~
2
3
:
2
3
2
3
~
,
~
,
~
-
-
-
-
-
-
-


-
-
+
+
=
=
=
=
=
+
+
-
=
-
+
+
-
=
-
-
=
-
-
=
NR
E
eq
eq
v
eq
v
eq
eq
v
E
v
p
eq
E
eq
v
v
p
eq
eq
v
v
p
T
p
T
NR
K
p
R
p
B
p
B
p
With
S
p
L
S
p
B
p
D
p
R
p
L
S
p
B
p
D
p
R
p
L
S
p
With
S
p
D
p
R
p
L
p
L
p
L
p
L
p
L
p
R
S
B
B
With
S
D
R
p
R
S
S
p
D
p
R
p
R
p
G
p
G
p
S
S
S
S
S
S
E
E
E
E
with
Finally,
p
E
~
puts itself in the form:
()
()
()
()
(
)
(
)
(
)
1
2
-
p
E
E
E
E
E
Id
~
~
~
~
~
=
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
3
2
3
2
1
2
1
1
1
1
2
1
2
1
2
1
2
1
2
2
2
2
p
D p
I
p
I
p
I
p
H
H
H
H
H
H
S
has
has
S
has
has
S
has
has
with:
()
() ()
()
()
() ()
()
()
() ()
()
()
() ()
()
()
() ()
()
()
() ()
()
()
() ()
()
()
() ()
()
()
() ()
I
p
I p L p
L p
I
p
I p L
p
L p
I
p
I p L
p
L p
H p
H p L p
L p
H
p
H p L
p
L p
H
p
H p L
p
L p
H p
H p L p
L p
H
p
H p L
p
L p
H
p
H p L
p
S
has
has
S
has
has
S
has
has
= -
= -
= -
= -
= -
= -
= -
= -
= -
3
2
3
2
3
2
3
2
3
2
3
2
3
2
3
2
3
2
1
3
1
21
3
2
22
3
1
1
1
3
1
1
1
21
3
2
1
1
22
3
2
2
1
3
1
2
2
21
3
2
2
2
22
.
.
.
()
L p
3
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
21/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Appendix 2 Resolution of the equation F (
p) = 0

It is a question of solving a nonlinear scalar equation by seeking the solution in a confidence interval.
For that, one proposes to couple a method of secant with a control of the interval of search. That is to say
the following equation to solve:
()
[]
()
()
F
,
F
F
X
X
B has
has
B
=
<
>
0
0
0
éq A2-1
The method of the secant consists in building a succession of points
X
N
who converges towards the solution. It is
defined by recursion (linear approximation of the function by its cord):
() () ()
X
X
X
X
X
X
X
N
N
N
N
N
N
N
+
-
-
-
-
=
-
-
-
1
1
1
1
1
F
F
F
éq A2-2
In addition, if
X
N
+
1
was to leave the interval, then one replaces it by the terminal of the interval in question:
if
then
if
then
X
has
X
has
X
B
X
B
N
N
N
N
+
+
+
+
<
=
>
=




1
1
1
1
:
:
éq A2-3

On the other hand, if
X
N
+
1
is in the interval running, then one reactualizes the interval:
[]
()
[]
()
if
and
then
if
and
then
X
B has
X
has
X
X
B has
X
B
X
N
N
N
N
N
N
+
+
+
+
+
+
<
=
>
=



1
1
1
1
1
1
0
0
,
F
,
F
éq A2-4

One considers to have converged when
F
is sufficiently close to 0 (tolerance to be informed). As for both
first leader characters, one can choose the terminals of the interval, or, if one has an estimate of
the solution, one can adopt this estimate and one of the terminals of the interval.

Note:
This method functions well if there is only one solution in the interval
[]
B has
,
. Without that being
formally shown, one can note that
()
F 0
0
>
.
One seeks then
B
such as
()
F B
<
0
.
One leaves for that
()
B
C
C
R p
E
eq
=
-




-
-
-
-
~
µ
2
3
2
3
3
1 1
2 2
If
()
F B
is
>
0
, one multiplies
B
by 10 and one tests if
()
F B
>
0
, and so on, until finding a value
B
such as
()
F B
<
0
.
One is sure that there is then at least a solution on
[]
B has
,
.
background image
Code_Aster
®
Version
6.4
Titrate:
Behavior élasto-visco-plastic of J.L.Chaboche
Date:
08/12/03
Author (S):
P. of BONNIERES,
J.M. PROIX, E.LORENTZ
Key
:
R5.03.04-B
Page
:
22/22
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A

























Intentionally white left page.