Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
1/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
Organization (S):
EDF/IMA/MN
Manual of Reference
R7.02 booklet: Breaking process
R7.02.07 document
Rate of refund of energy in thermo elastoplasticity
Summary:
One presents the calculation of the rate of refund of total mechanical energy by the method théta in 2D or 3D
for an elastoplastic thermo problem. The relations of elastoplastic thermo behavior are described
in detail in the document [R5.03.02].
This rate of refund of total mechanical energy makes it possible to analyze the situations of loadings not
monotonous of the defect, for irreversible material behaviors.
Let us note that the problem of the elastoplastic thermo rupture is a delicate problem. It is advised of
to consult the references before a first use.
Caution:
The defect must be modelized by a notch and not by a fissure [§5].
Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
2/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
Contents
1 Choice of the formulation of the rate of refund of energy in thermo elastoplasticity .............................. 3
2 Relation of behavior ..................................................................................................................... 4
3 Lagrangian Expression of the rate of refund of energy in thermo elastoplasticity ......................... 6
4 Establishment in Code_Aster ........................................................................................................... 8
5 Restrictions ............................................................................................................................................ 8
6 Bibliography .......................................................................................................................................... 8
Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
3/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
1
Choice of the formulation of the rate of refund of energy in
thermo elastoplasticity
One considers a notched elastic solid occupying the field
space
R
2
or
R
3
. That is to say:
·
U
the field of displacement,
·
T
the field of temperature,
·
the tensor of the stresses,
·
F
the field of the voluminal forces applied to
,
·
G
the field of the surface forces applied to a part
S
of
,
·
U
the field of displacements imposed on a part
S
D
of
.
F
G
S
S
D
In thermo or not-linear linear elasticity, the rate of refund of energy
G
is defined by the opposite
derivative of the potential energy compared to the field
[bib1]:
G
W
= -
Total potential energy with the balance of the system is:
()
W U
F U
G U
=
-
-
D
D
D
I
I
I
S
I
where
is the density of free energy. In elasticity,
is equal to the density of free energy elastic
[R7.02.01]
One extends this definition for the elastoplastic thermo problem, while choosing to replace
by
total mechanical energy
~
. This choice is justified in the document [bib2].
~
is a function of the following variables of state:
·
the tensor of the total deflections,
·
p
the tensor of the plastic deformations,
·
T
the field of temperature,
·
p
the variable interns scalar isotropic work hardening (cumulated plastic deformation),
·
one or more tensorial or scalar variables of kinematic work hardening.
Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
4/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
(
)
(
)
()
()
~,
()
,
,
D
S
D
S
p
p
p
T
p
O
T
O
T
T p
D
D
D
p T
D
D
D T
D
T
=
=
+
+
0
where
is the density of free energy
is the density power density of plastic dissipation
is the density of entropy
()
It is noted that
~
is the density of increased free energy of dissipated voluminal energy
plastically during all the evolution, and to which energy is added
()
S D
D T
D
O
T
(contribution of the temperature to the variation of free energy).
Caution:
One limits oneself to a solid
notched (cf [§5]).
2
Relation of behavior
The behavior of the solid is supposed to be thermo elastoplastic associated a criterion of Von Mises
with isotropic or kinematic work hardening. This type of behavior is currently treated in
the operator
STAT_NON_LINE
[U4.32.01] under the key word factor
COMP_INCR
. Treated relations
in this document are:
VMIS_ISOT_LINE
: Von Mises with linear isotropic work hardening,
VMIS_ISOT_TRAC
: Von Mises with isotropic work hardening given by a traction diagram,
VMIS_CINE_LINE
: Von Mises with linear kinematic work hardening.
For more details, to consult the documents [R5.03.02] and [U4.32.01].
is connected to the field of displacement
U
by:
()
(
)
U
U
U
=
+
1
2
I J
J I
,
,
The density of free energy is written:
(
)
(
)
,
,
,
,
H (,) Z ()
E
p
p
T p
T
p
T
T
=
+
+
where
H
is the density of energy of work hardening
Z
an arbitrary function of the temperature
E
density of thermo energy elastic defined by:
(
)
(
)
(
)
(
)
(
)
E
,
,
p
ij ije
ijkl
ij
ijp
ref.
ij
kl
p
ref.
T
T T
kl
T T
kl
=
=
-
-
-
-
-
-
1
2
1
2 A
,
with
the thermal expansion factor, and (
With
ijkl
) the tensor of elasticity.
In the particular case where there no was plastic evolution, one finds the expression of the density
of elastic energy for a thermo problem elastic linear [R7.02.01 §1.1].
Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
5/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
Free energy of work hardening
H
is deduced from:
(
)
(
)
H
,
R,
p p T
p T
=
, for isotropic work hardening where
(
)
R,
p T
is the radius of the surface of load
and
(
)
(
)
H,
X,
T
T
=
, for kinematic work hardening where
(
)
X,
T
is the translation of surface
of load in the space of the stresses (in the case of linear kinematic work hardening
=
p
)
For the relation of behavior of Von Mises with linear isotropic work hardening:
(
)
R,
p T
E E
E E p
T
T
= -
E
T
E
y
Traction diagram
Characteristics of the material (Young Modulus
E
and
D_SIGM_EPSI
:
E
T
) can depend on
temperature [R5.03.02 §3.2.1].
Plastic dissipation for a law of behavior of Von Mises checks:
()
D
p
D
p
O
T
y
=
where
y
is the initial linear elastic limit.
Finally density of total mechanical energy
~
is written:
(
)
(
)
()
()
(
)
(
)
(
)
(
)
~,
,
,
D
S
!
R,
X,
p
p
p
O
T
O
T
ij
ij
ijp
ref.
ij
O
p
y
O
T p
p
T
D
T
D
T T
p T dp
T D
p
=
+
+
=
-
-
-
+
+
+
1
2
For a linear isotropic work hardening:
(
)
(
)
(
)
~
~,
=
=
-
-
-
+
-
+
p
ij
ij
ijp
ref.
ij
T
T
y
T p
T T
E E
E E p
p
1
2
1
2
2
Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
6/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
3
Lagrangian expression of the rate of refund of energy
in thermo elastoplasticity
The rate of refund is calculated in Code_Aster by the method théta [R7.02.01 §1.3]. One notes by
!Q
the Lagrangian derivative of the quantity
Q
in a virtual propagation of the notch of
,
being a small real parameter and
a field of vector representing the direction of propagation of
the notch (one thus has
!Q
Q
Q.
((),)
X
=
+
).
The rate of refund of total mechanical energy in this propagation
is:
()
(
)
(
)
-
=
-
+
-
-
+
-
=
+
+
+
+
+
=
= -
= -
G
~
~
~!
,
,
~. !
~. ! ~! ~! ~. ! S!
~
~
~
S
~
,
,
F U
F U
G U
G U
N N
I
I
I
I
K K
I
I
I
I
K K
K
K
S
p
I
ij
p
ijp
I
I
ij
ij
ijp
ij
D
D
T p
T T
p p
T
T
However
with
ij
ij
ij
ij
(
)
(
)
p
p T
p
X
T
Y
I
I
I
=
+
=
R,
~
,
Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
7/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
()
()
(
)
that is to say
conventional terms
with
-
=
-
+
+
+
+
+
=
+
-
+
=
+
=
+
=
+
=
+
G
ij
ij
ij
ijp
I
I
K K
ij
I J
J I
I K K J
J K K I
ijp
ijp
ij K
p
K
K K
K K
ij
ij
ij K K
T T
p
D
F G
T
T
p
p
!
!
~! R! ~! ~
,
!
!
!
!
!
,
,
,
,
,
,
,
,
,
,
,
1
2
1
2
U
U
U
U
T
p
One can eliminate
!U
expression of
()
G
by noticing that
!U
is kinematically acceptable and
by using the equilibrium equation [R7.02.01 §1.3]. In the same way, terms
(
.
R
X.
)
-
+
+
+
p
Y
p
p
are eliminated as well as the terms
~.
T
T
S T
-
.
The following expression then is obtained:
()
(
)
(
)
()
()
(
)
()
-
=
-
-
+
+
+
+
+
=
-
+
-
+ +
+
-
+
G
G
~
~
R,
~
~
~
R
~
,
,
,
,
,
,
,
,
,
,
,
,
,
,
K K
ij
I K
K J
ij
ij K
p
K
K
K
y
K
K
ij
ij K
K
K K
ij
I K
K J
K
y
K
ij
ij K
ij
ij K
p
K
U
T T
p T
p
D
U
T T
p
D
conventional terms
and finally
conventional terms
F G
F G
,
,
:
For a radial and monotonous loading:
(
)
ij
ij K
p
y
K
ij
ij K
p
,
,
,
R
~
=
+
+
and one finds
the expression of
()
G
in nonlinear thermoelasticity [R7.02.03].
Code_Aster
®
Version
4.0
Titrate:
Rate of refund of energy in thermo elastoplasticity
Date:
25/04/97
Author (S):
G. DEBRUYNE, E. SCREWS
Key:
R7.02.07-A
Page:
8/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/97/006/A
4
Establishment in Code_Aster
The comparison of the formulations of
()
G
in thermo linear elasticity and thermo elastoplasticity
show that the terms of the two formulations differ only by terms from transport from
internal variables.
The presence of the key word factor
COMP_INCR
, and of the key word factor
RELATION
:
“VMIS_ISOT_LINE”
(or
“VMIS_ISOT_TRAC”
or
“VMIS_CINE_LINE”
) indicates that it is necessary to recover the field
displacements
U
, stresses
, and characteristics of elastoplastic material. It is
also necessary to recover the fields of the tensors of plastic deformation by the operator
CALC_ELEM
[U4.61.02].
The types of finite elements which support these options are the same ones as in elasticity [R7.02.01 §2.4].
They are the isoparametric elements 2D and 3D.
The supported loadings are the same ones as in the elastic case.
5 Restrictions
Caution:
This formulation of G for an elastoplastic thermo relation is valid only for one
solid notched and not for a fissured solid. One will choose for example (but the user will be able
to choose its own regular notch):
OK
OK
NOT
Indeed, the main difficulty in the establishment of this formulation is impossibility of showing
the existence of derived from total mechanical energy for a field comprising a fissure, and this
mainly by the absence of knowledge of the singularities of the fields in plasticity. For
to circumvent the problem, one regularizes the field by representing the defect in the form of notch. For
more details, it is advised to consult [bib2].
The validation of this formulation is carried out in test SSNP102 [V6.03.102] - Calculation of the rate of
restitution of energy for an elastoplastic problem.
6 Bibliography
[1]
BUI H.D.: Fragile breaking process, Masson, 1977.
[2]
DEBRUYNE G.: Proposal for an energy parameter of ductile rupture in
thermo plasticity, HI-74/95/027/0, 23/02/96