Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
1/12
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.18
Law of damage of an elastic material
fragile
Summary:
This document describes the elastic model of behavior fragile ENDO_FRAGILE available in statics and in
dynamics. The damage is modelized in a scalar way; loadings in compression and in
traction are not distinguished. In addition to the local model, the nonlocal formulations with gradient
of damage and with regularized deformation are also supported to control the phenomena of
localization.
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
2/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
1 Field
of application
Law ENDO_FRAGILE aims at modelizing in the manner a simplest possible behavior
fragile rubber band. The material is elastic isotropic. Its rigidity can decrease in an irreversible way
when the deformation energy becomes important, without distinguishing traction from compression.
This loss of rigidity is measured by a variable interns scalar called damage which
evolve/move of 0 (healthy material) to 1 (completely damaged material, i.e. without rigidity). Moreover,
the stress cannot exceed a threshold which also decrease him with the level of damage for
to reach 0 when the material is completely damaged. One will refer to [bib1] for one
description of this type of phenomenology.
The property of decrease of the threshold in stress with the level of damage is called
softening and generally involves a loss of ellipticity of the equations of the problem. It results from it
a localization of the deformations and damage in tapes of which the thickness is
directly controlled by the size of the finite elements. To mitigate this deficiency of the model, two
nonlocal formulations are proposed, one founded on the introduction of the gradient of
the damage and activated by modeling
*
_GRAD_VARI
[R5.04.01], the other resting on one
regularization of the deformations and activated by modeling
*
_GRAD_EPSI
[R5.04.02]. In a case
as in the other, the width of the tapes of localization is henceforth controlled by a parameter
material, well informed in operator DEFI_MATERIAU under key word LONG_CARA of the key word factor
NON_LOCAL [U4.43.01]. However, obtaining a physical problem posed again well is not
obtained that at the price of an important overcost in time calculation. In addition, it should well be noticed that
only the relations of behavior are deteriorated and not the equilibrium equations. Consequently,
the stresses preserve their usual direction.
Lastly, that one activates or not these nonlocal formulations, softening character of the behavior
also involve the appearance of instabilities, physics or parasites, which result in
snap-backs on the total answer and returns the piloting of the essential loading in statics.
piloting of the type PRED_ELAS [R5.03.80] then seems the mode of control of the level of
the most suitable loading.
2
Local law of behavior
2.1
Relations of behavior
The state of material is characterized by the deformation
and the damage
D
ranging between 0 and 1.
The forced relation deformation is elastic, rigidity is affected in a linear way by
the damage:
(
)
E
-
=
D
1
éq
2.1-1
with
E
the tensor of Hooke. In addition, the evolution of the damage, always increasing, is
controlled by the following function threshold:
()
()
()
2
1
1
K
K
2
1
,
F
-
+
+
=
-
=
D
W
D
D
D
y
where
E
éq
2.1-2
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
3/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Coefficients
y
W
and
, both positive, is parameters of the model. The condition of
coherence then determines completely the rate of damage
d&
:
()
()
0
,
F
0
0
,
F
=
D
D
D
D
&
&
éq
2.1-3
The equations [éq 2.1-1] with [éq 2.1-3] are enough to entirely describe the law of behavior
ENDO_FRAGILE, indeed very simple. One can also notice that it fits in
formalism suggested by Marigo [bib2].
2.2
Identification of the parameters of the model
The parameters of this law of behavior are four. On the one hand, the module of
Young
E
and the Poisson's ratio
who determine the tensor of Hooke by:
()
Id
E
tr
1
1
E
E
-
+
=
-
éq
2.2-1
In addition,
y
W
and
who define the lenitive behavior. They are determined by a test
of simple traction, cf [Figure 2.2-a]. To simplify the input of the data of the model, one informs not
not
y
W
and
but directly the tangent module
T
E
and the stress with the peak
y
under the key word
factor ECRO_LINE or ECRO_LINE_FO of operator DEFI_MATERIAU. As for
E
and
, they are
given classically under the key word factor ELAS or ELAS_FO.
For whatever purpose it may serve, here also expressions of the deformation with rupture
R
in this test of
simple traction, as well as voluminal energy
0
K
consumed to damage a point completely
hardware, this last expression being valid whatever the history of loading:
+
=
=
-
=
-
=
1
2
1
1
1
2
1
1
1
2
0
y
y
R
y
T
y
T
R
W
E
E
K
E
E
éq 2.2-2
E
W
y
y
2
2
=
E
E
T
-
=
Appear 2.2-a: Simulation of a simple tensile test
y
R
0
<
T
E
E
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
4/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
2.3
Integration of the law of behavior
Temporal discretization of the equations [éq 2.1-1] with [éq 2.1-3] on a pitch of time
]
[
T
T
-
is
realized by a diagram of implicit Euler. For any function of time
Q
, one notes
()
-
-
= T
Q
Q
and
()
T
Q Q
=
. To integrate in time the law of behavior then means to determine the state of stress and
of damage solution of the following nonlinear system, where deformation
and the state of material
at the beginning of the pitch of time
(
)
-
-
D
,
are given:
(
)
E
-
=
D
1
éq 2.3-1
()
(
)
()
0
,
F
0
0
,
F
=
-
-
-
-
D
D
D
D
D
D
éq
2.3-2
A method of resolution was proposed by [bib3]. It starts by examining the solution without
evolution of the damage (also called elastic test) then, if necessary, proceeds to one
correction to check the condition of coherence. In this case, the existence and the unicity of
solution guarantee the correct operation of the method. Let us consider the elastic test:
()
()
0
,
F
F
el
=
=
-
-
D
D
D
if
solution
éq
2.3-3
In the contrary case, the damage is obtained while solving
()
0
,
F
=
D
:
(
)
E
=
-
+
=
2
1
1
1
W
W
W
D
y
where
éq
2.3-4
As for the stress, it is given by [éq 2.3-1] in all the cases.
It still remains to be made sure that the damage does not exceed value 1. In fact, when
1
=
D
,
the rigidity of the material point considered is cancelled. Insofar as no technique of suppression
finite elements “broken” is not implemented (technical possibly delicate when them
finite elements have several points of Gauss), of the null pivots can appear in
stamp rigidity. This is why a numerical threshold is introduced
C
D
beyond which one is considered
elastic residual rigidity for the tangent matrix, the equations of remaining behavior
unchanged.
To preserve a reasonable conditioning of the matrix of rigidity, one chooses
5
10
1
-
-
=
C
D
. One
indicator
, arranged in the second internal variable, the behavior specifies then during the pitch
time running:
·
0
=
elastic behavior (deformation energy lower than the threshold)
·
1
=
evolution of the damage
·
2
=
(saturated damage) (
1
=
D
).
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
5/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
2.4
Description of the internal variables
The variables intern are two:
·
VI (1) damage
D
·
VI (2) indicating
3
Formulation with gradient of damage
3.1
Tally standard generalized
To mitigate the pathological localization of the deformations and the damage inherent in
softening character of material, a first alternative consists in extending the law of behavior
by taking of account the gradient of the damage on a material point scale [bib4]. This
approach is restricted with generalized standard materials. Model ENDO_FRAGILE enters
indeed within this framework; microscopic free energy
µ
and potential of dissipation
microscopic
µ
are written:
() (
)
()
()
[]
()
()
()
(
)
-
+
-
-
=
=
+
+
-
=
D
D
D
K
D
D
D
D
D
µ
1
1
2
1
I
1
,
0
bl
0
el
1
0
bl
0
el
E
where
éq
3.1-1
()
[
[
()
D
D
K
D
&
&
&
+
µ
+
=
0
0
I
éq
3.1-2
where
K
I
indicate the indicatrix of the convex unit
K
, null in
K
and being worth
+
elsewhere. One
the convexity of these potentials compared to the damage, strict convexity will notice for
free energy: that confirms the existence and the unicity of the solution for the integration of the law.
To preserve a simple model, one will be satisfied to introduce a quadratic term in gradient
of damage in the free energy. The potential of dissipation remains unchanged. Consequently, them
macroscopic potentials are written:
(
) (
) ()
()
[]
()
·
·
·
·
+
+
+
-
=
D
D
D
C
D
D
D
D
2
1
I
1
,
,
1
0
bl
0
el
éq
3.1-3
(
)
[
[
()
·
+
·
·
+
=
D
D
K
D
&
&
&
&
0
0
I
, D
éq
3.1-4
where the factor
C
depends on the characteristic length
B
L
material (distance characteristic
of interaction enters the microscopic cracks given under key word LONG_CARA of DEFI_MATERIAU) of
following manner:
()
(
)
2
0
2
2
2
2
1
13
8
13
4
0
,
B
B
L
K
L
D
C
+
=
=
µ
éq
3.1-5
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
6/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
The macroscopic potentials thus built then make it possible to define the total potentials in
the origin of the variational formulation of the law of behavior. In this relatively simple case,
one can interpret this variational formulation in term of laws of local state and evolution of which them
expressions are close to those for the local law:
(
)
E
-
=
D
1
éq
3.1-6
()
()
D
C
D
D
+
-
=
K
2
1
,
F
G
E
éq
3.1-7
()
()
0
,
F
0
0
,
F
G
G
=
D
D
D
D
&
&
éq
3.1-8
The essential difference with the local model lies in the definition of the function threshold which depends
henceforth Laplacian of the damage. This last confers on the condition coherence one
differential character. The interpretation of the variational formulation also provides the conditions
with the limits for the field of damage as well as the conditions of interface (or jump):
N
N
normal
of
edge
on
0
=
D
éq
3.1-9
[
]
N
N
normal
of
interface
one
of
lies
crosses
with
0
=
D
C
éq
3.1-10
3.2
Integration of the law of behavior
The integration of the laws of behavior to gradients of internal variables is described in the booklet
[R5.04.01]. It is based inter alia on a local stage which depends explicitly on the relation on
behavior treated. It is the aforementioned which one proposes to describe here.
First of all, it is necessary to define the standard used to build the term of penalization of
Lagrangian increased, i.e. the matrix
Nrr
. As recommended, one adopts an estimate
diagonal of the matrix hessienne of energy compared to the variables of damage, which
conduit with:
(
)
(
)
+
=
+
=
=
·
·
·
·
13
4
1
2
1
2
2
0
2
0
2
B
L
K
NR
K
NR
NR
NR
NR
NR
NR
with
R
R
éq
3.2-1
In accordance with the general theory, it is now necessary to solve the following nonlinear system, in
which the state at the beginning of the pitch of time as well as the deformation
, multipliers of Lagrange
µr
,
the coefficient of penalization
R
and the nodal damage evaluated at the point of current Gauss
D
Br
R =
are given:
()
()
(
)
-
-
-
+
µ
+
-
D
D
D
NR
R
D
D
D
R
R
R
R
R
R
R
R
R
R
,
éq
3.2-2
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
7/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
where, once more, the compact notation
(
)
·
=
D
,
D
Dr.
is used. Being given decoupling enters
the value and the gradient of the damage in the free energy and the potential of dissipation, it
system [éq 3.2-2] makes it possible to separately treat the gradient of the damage and its value:
()
(
)
0
=
-
+
+
-
D
µ
D
D
NR
R
éq
3.2-3
(
)
(
)
(
)
-
·
·
µ
·
·
·
·
·
·
·
-
-
+
µ
+
-
D
D
D
NR
R
D
D
,
éq
3.2-4
On the one hand, the equation [éq 3.2-3] makes it possible to determine the gradient of damage immediately:
+
+
=
NR
R
C
NR
R
µ
D
éq
3.2-5
As for the equation [éq 3.2-4], it is interpreted like an equation of coherence:
()
(
) ()
()
0
F
0
F
0
0
F
K
With
With
With
D
D
D
D
With
R
R
R
R
-
=
=
-
-
·
·
·
-
·
·
-
·
·
·
with
éq
3.2-6
in which the thermodynamic force
R
With
·
has as an expression:
(
)
(
)
·
·
·
·
·
·
·
·
-
+
µ
+
-
=
D
NR
R
D
D
With
R
,
éq
3.2-7
The solution of this equation of coherence is similar to that for the local model
[éq 2.3-1] - [éq 2.3-4]. Initially, an elastic test then a correction so necessary. Being given
strict convexity of the potentials, the solution is single. When the elastic test is not solution,
resolution of the equation
()
0
F
=
·R
With
conduit with:
()
E
=
+
=
+
+
·
·
·
·
·
·
·
·
2
1
K
W
D
D
NR
R
NR
R
W
where
µ
éq
3.2-8
After some handling, this equation is reduced to search roots of a polynomial of
degree 3 in
·
D
, which one knows that only one (that which one seeks) is higher than
-
·
D
.
3.3
Estimate of the Laplacian of the damage
Like the local law, the law with gradients can lead to structural instabilities. Piloting by
elastic prediction PRED_ELAS can then appear essential. But this last requires
given of a threshold of elasticity point by point, whose existence is not necessarily acquired in
presence of nonlocal models. In fact, in this case, such a threshold exists thanks to
the interpretation [éq 3.1-6] - [éq 3.1-8] of the variational formulation of the behavior: it is the function
G
F
.
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
8/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
However, the evaluation of this threshold requires to know the Laplacian of the field of damage
who is not calculated by the algorithm of integration of the nonlocal laws. To know the nodal field
to damage is not enough besides because the functions to form, polynomials of degree 1,
do not allow to estimate a derivative directly second. Fortunately, the recourse to an artifice
allows despite everything to estimate the Laplacian with a precision sufficient for the algorithm of piloting.
Indeed, at the end of a pitch of time, the equation of coherence [éq 3.1-8] is checked, and this of as much
better than the precision requested by the user (key word RESI_DUAL_ABSO under the key word factor
LAGR_NON_LOCAL) is large. In each point of Gauss, one can thus be confronted with two cases of
appear. That is to say the threshold is not reached and, in this case, the damage in this point does not vary. One
suppose whereas it is the same for the Laplacian:
-
-
=
=
D
D
D
D
éq
3.3-1
In the contrary case, the damage evolves/moves and the equation
()
0
,
F
G
=
D
is checked, so that one
in deduced:
()
-
=
E
2
1
K
1
D
C
D
éq
3.3-2
In any rigor, this does not make it possible to determine the Laplacian when the damage is saturated.
However, that does not have importance since these points are not taken into account to control it
loading (one seeks points which dissipates).
3.4
Description of the internal variables
In the case of the formulation with gradient of damage, the variables intern are now with
numbers of six:
·
VI (1) damage
D
·
VI (2) gradient of the damage
D
, component
X
(
R
into axisymmetric)
·
VI (3) gradient of the damage
D
, component
y
(
Z
into axisymmetric)
·
VI (4) gradient of the damage
D
, component
Z
(0 in 2D)
·
VI (5) estimate of the Laplacian of the damage
D
(as long as
1
<
D
)
·
VI (6) indicating
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
9/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
4
Formulation with regularized deformation
4.1
Formulation continues in time
The approach with regularized deformation [R5.04.02] also makes it possible it to control the phenomena of
localization and for this reason seems an alternative to the formulation with gradient
of damage. But at the difference in the latter, this formulation has the advantage of
to resort to the standard algorithms for the nonlinear problems. Indeed, the only difference by
report/ratio with the local law of behavior lies in the data of two deformations instead of one,
local deformation
who intervenes in the forced relation deformation and the deformation
regularized
who controls the evolution of the damage. The aforementioned results from the local deformation
by resolution of the system of partial derivative equations according to:
=
=
-
N
N
0
0
normal
of
edge
on
structure
in
2
B
L
éq
4.1-1
where the characteristic length
B
L
is again well informed under key word LONG_CARA of
DEFI_MATERIAU. Finally, the relation of behavior is written in the following way, where
function threshold F was already defined in [éq 2.1-2]:
(
)
E
-
=
D
1
éq 4.1-2
()
()
0
,
F
0
0
,
F
=
D
D
D
D
&
&
éq
4.1-3
4.2
Integration of the law of behavior
One of the advanced advantages for the nonlocal formulation with regularized deformation is the little of
amendments which it involves in the construction of the law of behavior. Indeed, the integration of
internal variables is completely controlled by the regularized deformation
. They thus are found
expressions of the local law:
()
()
()
()
(
)
E
=
-
+
=
>
=
=
=
-
-
-
2
1
1
1
0
,
F
F
0
,
F
F
el
el
W
W
W
D
D
D
D
D
y
with
if
if
éq 4.2-1
The stress is then obtained directly by the relation [éq 4.1-2]. Moreover, one preserves well-sure
the introduction of a damage criticizes [éq 2.3-5] to preserve a residual rigidity.
4.3 Variables
interns
They are the same internal variables as for the local law:
·
VI (1) damage
D
·
VI (2) indicating
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
10/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
5
Piloting by elastic prediction
The piloting of the type PRED_ELAS controls the intensity of the loading to satisfy some
equation related to the value of the function threshold
el
F
during the elastic test [bib5]. Consequently, only
the points where the damage is not saturated will be taken into account. The algorithm which takes in
charge this mode of piloting, cf [R5.03.80], requires the resolution of each one of these points of Gauss of
the following scalar equation in which
is a data and
the unknown factor:
()
=
el
f~
éq 5-1
The function
el
f~
provides the value of the function threshold during an elastic test when the field of
displacement breaks up in the following way according to the scalar parameter
:
1
0
U
U
U
+
=
éq 5-2
where
0
U
and
1
U
are given. Thanks to the linearity in small deformations of the operators deformation
(calculation of the deformations starting from displacements) and regularized deformation, one also obtains
following decompositions:
1
0
1
0
+
=
+
=
and
éq
5-3
Consequently, whatever the adopted, local, nonlocal modeling with gradient
of damage or not local with regularized deformation, the function
el
f~
can put itself under
following form:
()
(
)
(
)
()
()
()
=
=
=
-
=
=
=
=
=
=
-
+
+
=
-
-
-
-
D
S
D
C
D
S
D
S
S
K
K
K
2
1
f~
1
1
0
0
1
1
0
0
1
1
0
0
1
0
1
0
el
E
E
E
E
E
E
E
E
E
E
E
E
regularized
N
déformatio
ent
endommagem
of
gradient
local
law
éq
5-4
The equation [éq 5-1] is reduced thus to search roots of the polynomial of degree 2 following:
()
=
=
-
=
+
+
=
1
1
0
1
0
1
0
0
0
2
2
1
0
2
1
2
12
1
2
P
E
E
E
E
E
E
E
E
E
P
P
P
P
P
P
where
éq
5-5
One provides to the algorithm piloting the linear approximation of P in the vicinity of the roots (real)
when they exist or if not the value of
for which P reaches its minimum.
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
11/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
6 Bibliography
[1]
LEMAITRE J., CHABOCHE J.L. : Mechanics of solid materials. Dunod: Paris, 1988.
[2]
MARIGO J.J.
: Formulation of a law of damage of an elastic material.
Report of the Academy of Science, Paris 1981; series II, 292 (19): 1309-1312.
[3]
SIMO J.C., TAYLOR R.L.: Consist tangent operators for misses-independent elastoplasticity.
Methods computer in Applied Mechanics and Engineering 1985; 48: 101-118.
[4]
LORENTZ E., ANDRIEUX S.: With variational formulation for nonlocal ramming models.
International Newspaper off Plasticity 1999; 15: 119-138.
[5]
LORENTZ E., BADEL P.: With load control method for ramming finite element simulations.
International Newspaper for Numerical Methods in Engineering, submitted 2002.
Code_Aster
®
Version
6.4
Titrate:
Law of damage of a fragile elastic material
Date
:
18/11/03
Author (S):
E. LORENTZ
Key
:
R5.03.18-A
Page
:
12/12
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/03/005/A
Intentionally white left page.