background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
1/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
Organization (S):
EDF-R & D/MMC















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



Model of Rousselier for the ductile rupture



Summary

The model of Rousselier describes the damage due to the plastic growth of cavities in a metal. It
allows to modelize cracking and the ductile rupture. The relation of behavior is elastoplastic or
viscoplastic with isotropic work hardening. It allows the changes of plastic volume and is written in
small deformations. The writing in great deformations with a formulation of Simo and Miehe modified, in
the elastoplastic case only, is described in [R5.03.06].
This model is available in the control
STAT_NON_LINE
via the key word
RELATION =
“ROUSS_PR”
or
“ROUSS_VISC”
under the key word factor
COMP_INCR
and with the key word
DEFORMATION =
“PETIT_REAC”
.
This model is established for modelings three-dimensional (
3D
), axisymmetric (
AXIS
), in stresses
plane and in plane deformations (
C_PLAN, D_PLAN
).

One presents the writing and the digital processing of this model.
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
2/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
Count
matters
1
Introduction ............................................................................................................................................ 3
2
Notations ................................................................................................................................................ 4
3
Model of Rousselier ............................................................................................................................ 5
3.1
Derivation of the equations of the model ............................................................................................... 5
3.2
Equations of the model ....................................................................................................................... 7
4
Numerical formulation .......................................................................................................................... 8
4.1
Key words, given material and variables intern ......................................................................... 8
4.2
Expression of the model discretized ..................................................................................................... 8
4.3
Resolution of the nonlinear scalar equation ............................................................................... 10
4.4
Form of the tangent matrix of the behavior ................................................................... 10
5
Bibliography ........................................................................................................................................ 12
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
3/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
1 Introduction
The mechanisms at the origin of the ductile rupture of metals are associated the development of
cavities within material. Three phases are generally distinguished:
· germination: it is about the priming or nucleation of the cavities, in sites which correspond
preferentially with the particles of second phase present in material,
· growth: it is the phase which corresponds to the development itself of the cavities,
controlled primarily by the plastic flow of the metal matrix which surrounds these
cavities,
· coalescence: it is the phase which corresponds to the localization of the deformation between the cavities
to create macroscopic fissures.
The model of Rousselier [bib1], [bib2], [bib3] presented here is based on assumptions
microstructural which introduces a microstructure made up of cavities and of a matrix of which them
elastic strain negligible are compared with the plastic deformations. In this case, and in
the absence of nucleation of new cavities, porosity
F
, definite like the relationship between volume
cavity
C
V
and total volume
V
representative elementary volume, is directly connected to
macroscopic plastic deformation by:
(
)
p
C
F
F
V
V
F
F
F
&
&
tr
1
1
1
with
0
0
-
=
=
-
-
=
éq
1-1
where
0
F
indicate initial porosity,
O
and
are respectively the density in
configurations initial and current (one takes in the continuation
1
=
O
) and
p
&
the rate of deformation
plastic of total volume
V
.
The construction of the model rests on a thermodynamic and phenomenologic analysis which brings
to write the plastic potential
F
in the following form:
(
)
)
(
exp
,
,
F
1
1
1
p
R
F
D
F
p
m
eq
-




+
=
éq
1-2
where
/
=
is the stress of Kirchhoff,
is the stress of Cauchy,
R
isotropic work hardening
function of the cumulated plastic deformation
p
,
1
and
1
D
parameters of material. The presence
in the plastic potential of the hydrostatic stress
m
authorize the changes of volume
plastic.
In the event of nucleation of new cavities, one considers that the voluminal fraction created is
proportional to the cumulated plastic deformation. It is thus enough to replace
F
by
p
With
F
N
+
in
equations of the model.
N
With
is a parameter of material. The equation [éq 1-1] is not modified.
In the viscoplastic case, one writes the viscoplastic potential
vp
F
like a function of the potential
plastic
F
:
(
)
F
p
vp
,
,
F
F
=
éq 1-3
One will consider only the particular case such as:
m
HS
p






=
=
0
0
F
F
&
&

éq 1-4
who is reduced to a function power (law of the Norton type) when two parameters of material
0
&
and
0
are very large.
Thereafter, one presents the relations of behavior of the model of Rousselier and his integration
numerical.
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
4/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
2 Notations
One will note by:
Id
tensor second-order identity
II
tensor identity of the fourth command
With
tr
trace second-order tensor
With
~
With
deviatoric part of the tensor
With
defined by
Id
With
With
With
)
tr
3
1
(
~
-
=
m
With
hydrostatic part of the tensor
With
defined by
3
tr A
=
m
With
With
eq
equivalent value of von Mises defined by
With
eq
= 32 ~: ~
WITH A
:

doubly contracted product:
)
tr (
:
,
T
J
I
ij
ij
B
With
AB
B
With
=
=
tensorial product:
(
)
WITH B
=
ijkl
ij kl
WITH B
,
µ
, E,
, K
moduli of the isotropic elasticity
p&
speed of equivalent plastic deformation
p
p
p
&
&
&
~
:
~
3
2
=

In addition, within the framework of a discretization in time, all the quantities
Q
evaluated at the moment
precedent are subscripted by
-
, quantities evaluated at the moment
T
T
T
+
=
-
are not subscripted and
the increments are indicated by
. One has as follows:
Q
Q
Q
+
=
-
The numerical resolution is carried out by one
- method, with
1
0
. For all the quantities,
one defines:
Q
Q
Q
+
=
-
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
5/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
3
Model of Rousselier
We now describe the derivation of the equations of the model of Rousselier presented in
introduction.
3.1
Derivation of the equations of the model
It is supposed that the specific free energy breaks up into three parts: a hyperelastic part
who depends only on the elastic strain, a part related to the mechanism of work hardening and one
part related to the damage:
(
)
()
()
)
(
,
,
F
p
F
p
F
p
E
E
E
+
+
=
éq
3.1-1
The inequality of Clausius-Duhem is written (one does not consider the thermal part):
0
:
-
&
&
éq 3.1-2
expression in which
p
E
&
&
&
+
=
represent the rate of deformation.
Dissipation is still written:
0
:
-
-
+


-
F
F
p
p
p
E
E
&
&
&
&
:
éq
3.1-3
The second principle of thermodynamics then requires the following expression for the relation
elastic stress-strain:
E
=
éq 3.1-4
One defines the thermodynamic forces associated with the elastic strain, with the deformation
figure cumulated and with porosity in accordance with the framework of generalized standard materials:
()
E
E
=
éq 3.1-5
()
p
p
=
With
éq 3.1-6
()
F
F
=
B
éq 3.1-7
It remains then for dissipation:
0
B
With
:
-
-
F
p
p
&
&
&
éq
3.1-8
The principle of maximum dissipation applied starting from the viscoplastic potential
B)
With,
,
(
F
vp
allows
to deduce the laws of evolution from them from the plastic deformation, the plastic deformation cumulated and of
porosity, is:
=
vp
p
F
&
éq 3.1-9
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
6/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
With
F
-
=
vp
p&
éq 3.1-10
B
F
-
=
vp
f&
éq 3.1-11
It is supposed that
B)
With,
,
(
F
vp
is a function of the plastic potential
B)
With,
,
F (
and that this last
break up into two terms depending respectively on the second invariant on
coupled to A and of
first invariant of
coupled to b:
()
(
)
(
)
(
)
B
,
F
With
,
F
F
F
m
m
eq
VM
vp
+
=
=
éq
3.1-12
By assumption, the first term breaks up in an additive way like the potential of von Mises:
(
)
()
)
(
With
With
,
F
0
p
R
R
p
eq
eq
eq
VM
-
=
-
-
=
éq
3.1-13
Not to obtain a commonplace result, the decomposition of the second term must be multiplicative:
(
)
() ()
B
B
,
F
H
G
m
m
m
=
éq
3.1-14
Taking into account the equation [éq 1-1], the laws of evolution for
p
&
tr
and
f&
lead to the equality:
()
()
()
(
)
()
(
)
F
H
F
H
F
G
G
m
m
B
B
1
1
'




-
-
=
éq
3.1-15
The two members of this equation are functions of the two independent variables
m
and
F
,
thus they is equal to a constant of dimension the reverse of a stress, it is the parameter of
material
1
/
1
. The parameter without dimension
1
D
appears in the integration of
G
G/
:
()




=
1
1
1
exp
m
m
D
G
éq
3.1-16
The function
()
F
B
and the opposite function
()
B
1
H
F
=
are unknown. The simplest choice and more
naturalness is to take
H
H
1
, which gives:
()
()
F
H
H
=
B
B
1
éq 3.1-17
()
(
)
F
F
F
H
-
-
=
=
1
1
dB
D
B
1
éq
3.1-18
The plastic potential is written finally:
)
(
exp
F
1
1
1
p
R
F
D
m
eq
-




+
=
éq
3.1-19
The law of evolution for
p&
give:
()
()
F
dF
F
D
V
p
=
=
&
éq 3.1-20
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
7/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
The function
()
F
V
the viscosity of material defines. One will consider only the particular case such as:
()
m
HS
V






=
0
0
F
F
&
éq
3.1-21
who is reduced to a function power (law of the Norton type) when two parameters of material
0
&
and
0
are very large. Conversely one a:
0
)
(
F
=
-
p
S &
éq 3.1-22








=
-
m
p
HS
p
S
1
0
1
0
)
(
&
&
&
éq
3.1-23
In the case of plasticity independent of time, the preceding equation becomes
0
F
=
(criterion or
threshold of plasticity) and
p&
is given by the equation of consistency
0
F
=
&
if
0
F
=
and
0
=
p&
if
0
F
<
.
The equations of the model are now completely defined, in the case without nucleation of
new cavities. In the event of nucleation of new cavities, one considers that the voluminal fraction
created is proportional to the cumulated plastic deformation. It is thus enough to replace
F
by
p
With
F
N
+
in the equations of the model.
N
With
is a parameter of material. The equation [éq 1-1] is not
not modified.
3.2
Equations of the model
One summarizes the equations of the model deduced from the thermodynamic and phenomenologic analysis which
precede:
(
)
0
)
(
exp
1
0
1
0
1
1
1
=








-
-




+
+
=
-
m
m
N
eq
vp
p
HS
p
R
p
With
F
D
&
&
éq 3.2-1
(
)
[
]
E
II
Id
Id
:
2
µ
+
=
=
éq
3.2-2
0
1
1
F
p
With
F
N
-
-
-
=
éq 3.2-3
eq
eq
p
p
p
2
~
3
2
~
3
~
&
&
&
=
=
éq 3.2-4
(
)




+
=
1
1
exp
tr
m
N
p
p
With
F
D
p&
&
éq
3.2-5
(
)
p
F
With
F
&
&
tr
1
1
-
=
éq 3.2-6
with
1
1
=
With
, this parameter being introduced only for numerical reasons.
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
8/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
4 Formulation
numerical
4.1
Key words, given internal material and variables
For the foreseeable applications, the model was established under two distinct key words:
`
ROUSSELIER_PR
`for the plastic model with nucleation of cavities or `
ROUSSELIER_VISC
`for
the viscoplastic model without nucleation. That makes it possible to avoid useless numerical calculations.
corresponding simplified equations are obtained starting from the general equations while posing
respectively
0
0
=
or
0
=
N
With
.
The whole of the parameters of the model is provided under the key words factors `
ROUSSELIER
`or
`
ROUSSELIER_FO
`and `
TRACTION
`(to define the traction diagram) control
DEFI_MATERIAU
([U4.43.01]). Parameters of the viscoplastic model (
0
,
0
&
and
m
) are provided
by the key word `
ROUSSELIER_VISC
`.
The internal variables produced in Code_Aster are:
·
V1, cumulated plastic deformation
p
,
·
V2, porosity
F
,
·
V3 with V8, the tensor of elastic strain
E
,
·
V9, the indicator of plasticity (0 if the last calculated increment is elastic, 1 if solution
figure regular, 2 if singular plastic solution).
We now present the numerical integration of the law of behavior and give
the form of the tangent matrix (options
FULL_MECA
and
RIGI_MECA_TANG
).
4.2
Expression of the discretized model
The numerical resolution is carried out by one
- method, with
1
0
. For all the quantities
Q
, one defines:
Q
Q
Q
+
=
-
Q
Q
Q
+
=
-
The system of equations discretized is:
(
)
p
E
~
~
2
~
2
~
-
=
=
µ
µ
éq
4.2-1
(
)
p
E
m
K
K
-
=
=
tr
tr
tr
éq
4.2-2
eq
p
p 2 ~
3
~
=
éq 4.2-3
(
)




+
=
1
1
exp
tr
m
N
p
p
With
F
Pd
éq
4.2-4
(
)
p
F
With
F
-
=
tr
1
1
éq 4.2-5
(
)
0
)
(
exp
1
0
1
0
1
1
1
=








-
-




+
+
=
-
m
m
N
eq
vp
T
p
HS
p
R
p
With
F
D
&
éq 4.2-6
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
9/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
This system is reduced to the solution of only one equation scalar for the unknown factor
F
, knowing
,
T
and quantities
-
Q
. It is noted that
does not intervene in the algorithm, on the other hand it
will intervene in the calculation of the coherent tangent matrix. One calculates successively:
(
)




-
-
+
=
-
F
With
F
K
m
m
1
tr
1
éq
4.2-7
p
is the positive root of the quadratic equation:
()
(
)
(
)
(
)
0
/
exp
1
1
1
1
1
2
=
-
-
+
+
-
m
N
N
D
F
With
F
p
p
With
F
p
With
éq
4.2-8
[
]
(
)
~
2
~
~
2
~ 3
1
~
+




+
-
=
-
-
µ
µ
µ
eq
p
éq
4.2-9
[
]
p
eq
eq
-
+
=
-
µ
µ
3
~
2
~
éq
4.2-10
The scalar equation for
F
is the equation [éq 4.2-6]
0
=
vp
.
Notice 1:
Like
F
is very weak in most of the structure, it would be preferable
to use
p
like main unknown factor. But in this case it is not possible
to bring back to a scalar equation, which makes more difficult the use of a method of the type
Newton. It is also one of the reasons why the equations [éq 1-1], [éq 3.2-6] and
[éq 4.2-5] were not modified by the introduction of the nucleation of the cavities.
Notice 2:
The equation [éq 3.2-6] can be integrated exactly:




-
-
=
F
F
With
p
1
1
ln
1
tr
0
1
from where:




-
-
=
-
F
F
With
p
1
1
ln
1
tr
1
Like the numerical parameter
1
With
can be modified in a discontinuous way, the derived form
[éq 4.2-5] was preserved, including in the calculation of the coherent tangent matrix. If
the use of the parameter
1
With
was to be abandoned in a later version, it would be necessary
to consider the use of the integrated form.
Notice 3:
The integrated form
0
=
vp
is used, including in plasticity instead of the relation of
consistency
0
F
=
&
who gives
p&
. The coherent tangent matrix is calculated with this
form integrated.
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
10/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
4.3
Resolution of the nonlinear scalar equation
The resolution of the equation
()
0
=
F
vp
be carried out by an algorithm of Newton on controlled terminals
in the routine
LCROUS
.
()
F
vp
and its derivative compared to
F
are calculated in the routine
RSLPHI
called by
LCROUS
. The initial values of the terminals are:
· limit lower:
0
1
=
F
since
()
0
0
<
vp
(it was checked as a preliminary that the branch
rubber band (negative threshold) is not solution),
· limit higher:
2
F
such as
()
0
0
>
vp
sought by dichotomy between 0 and
-
- F
1
(first value for this search:
2
1
-
- F
).
The algorithm of Newton begins with the value
0
=
F
. Whatever the value found for
F
one
thus note for the continuation that the function
()
F
vp
and its derivative compared to
F
are at least
calculated for
0
=
F
and
2
1
-
- F
.
The developments carried out to improve convergence and the robustness of the algorithm are
described in [bib5].
4.4
Form of the tangent matrix of the behavior
One gives the form of the tangent matrix here (option FULL_MECA during iterations of
Newton, option RIGI_MECA_TANG for the first iteration).
For option RIGI_MECA_TANG, the tangent operator is the same one as that which connects
E
with
in
[éq 3.2-2].
For option FULL_MECA, the tangent matrix is obtained by linearizing the system of equations which
governs the law of behavior: [éq 4.2-1] with [éq 4.2-6]. It is thus about a coherent tangent matrix.
To simplify the expressions, one notes in this paragraph [§4.5]:
Q
for
Q
, quantities being
all expressed at the moment
T
T
T
+
=
-
. The coherent tangent matrix is:






+


-
+


+
+


+
-
+
=
Id
Id
Id
Id
II
~
1
3
3
~
~
~
3
5
1
4
5
4
2
3
1
3
K
y
K
has
y
has
has
has
has
has
has
éq 4.4-1
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
11/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
This operator is calculated in the routine
RSLJPL
. The coefficients are calculated as follows:
)
(
3
2
7
1
1
p
Z
Z
K
y
K
has
eq
+
+
=
éq
4.4-2
1
3
1
2
)
(
µ
y
y
has
+
=
éq 4.4-3
5
3
2
Z
has
eq
µ
=
éq 4.4-4
2
2
4
3
X
y
has
µ
=
éq 4.4-5
1
1
5
3
µ
y
has
=
éq 4.4-6
1
2
6
3
µ
eq
has
p
K
has
-
=
éq 4.4-7
eq
N
X
p
With
F
Z
Kz
y
1
1
6
1
)
(
3
+
-
=
éq 4.4-8
2
5
1
2
3
eq
Z
X
y
µ
-
=
éq 4.4-9
eq
N
X
p
With
Z
Kz
y
1
1
6
3
3
-
=
éq 4.4-10
p
Z
Z
Z
Z
Z
With
y
+
+
=
2
7
1
9
1
8
1
4
éq 4.4-11
)
(
2
7
6
9
1
8
2
1
5
p
Z
Z
has
Z
Z
Z
has
With
y
eq
+
-
=
éq
4.4-12




+
+
=
1
1
1
1
exp
)
(
1
m
N
p
With
F
Pd
With
Z
éq
4.4-13
vp
R
Z
+
=
µ
3
2
éq 4.4-14
)
1
(
)
(
1
1
1
3
F
With
Z
p
With
F
K
Z
N
-
-
+
=
éq
4.4-15
eq
vp
p
R
Z
-
=
4
éq 4.4-16
p
Z
eq
+
=
µ
3
5
éq 4.4-17




=
1
1
6
exp
m
D
Z
éq 4.4-18
background image
Code_Aster
®
Version
6.4
Titrate:
Model of Rousselier for the ductile rupture
Date:
10/09/04
Author (S):
G. ROUSSELIER, R. MASSON, G. BARBER
Key
:
R5.03.07-A
Page
:
12/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/04/002/A
)
(
1
6
7
p
With
F
Z
Z
N
+
=
éq 4.4-19
p
With
F
F
Z
N
-
-
-
= 1 1
8
éq 4.4-20
p
With
F
With
Z
N
N
-
-
= 1
9
éq 4.4-21
3
1
2
1
2
7
6
3
1
)
(
X
Z
Z
p
Z
Z
Z
Z
X
-
+
+
=
éq
4.4-22
p
X
Z
Z
Z
Z
p
Z
Z
X
+
-
+
-
=
3
1
4
1
7
4
6
3
2
)
(
éq
4.4-23
2
1
6
1
3
Z
Z
With
X
N
=
éq 4.4-24
p
D
T
p
dS
T
dp
p
Dr.
R
vp
&
)
/
(
1
)
(
+
=
éq
4.4-25
For the plastic model with nucleation of cavities `
ROUSSELIER_PR
`and for the model
viscoplastic without nucleation `
ROUSSELIER_VISC
`, the corresponding simplified equations are
obtained starting from the equations above while posing respectively
dp
p
Dr.
R
vp
/
)
(
=
and
0
=
N
With
.


5 Bibliography
[1]
ROUSSELIER G.: “Finite constitutive deformation relations including ductile fracture
ramming ", in Three-Dimensional Constitutive Relations and Ductile Fracture, ED. Nemat-
Nasser, North Holland publishing company, pp. 331-355, 1981.
[2]
ROUSSELIER G.: “The Rousselier model for porous metal plasticity and ductile fracture”, in
Handbook off Materials Behavior Models, ED. J. Lemaitre, Academic Close, pp. 436-445,
2001.
[3]
ROUSSELIER G.: “Dissipation in porous metal plasticity and ductile fracture”, J. Mech. Phys.
Solids, vol. 49, pp. 1727-1746, 2001.
[4]
BARBER G.: “Model of Rousselier in Code_Aster: new implementation”, notes
EDF R & D HT-2C/98/007/A, 1998.
[5]
MASSON R., ROUSSELIER G., BONNAMY Mr.: “Improvement of convergence and evolution
implementation of the model of Rousselier in Code_Aster ", note EDF R & D
HT-26/01/037/A, 2002.
[6]
MIALON P.
: “Elements of analysis and numerical resolution of the relations of
elastoplasticity ", EDF, bulletin of DER, data-processing series C mathematical, 3, pp. 57-89,
1986.