Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
1/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
Organization (S):
EDF-R & D/AMA
Manual of Reference
R7.02 booklet: Breaking process
Document: R7.02.11
Law of behavior of Barenblatt
Summary:
The law of behavior of Barenblatt is a law of interface making it possible to modelize the opening of a fissure
while taking account of a force of cohesion enters the lips of the aforementioned. An energy of surface allows
to take into account the energy cost of the opening of the fissure. The latter will be represented by
finite elements of gasket type. An energy of penalization will make it possible to take into account the condition of not
interpenetration of the lips of the fissure.
We present here the form of energy and the stress which in drift according to the jump of displacement
as well as the internal variables. The existence of instability requires a piloting by the elastic prediction of which one
will detail the elements specific to this law.
Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
2/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
Count
matters
Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
3/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
1
Law of Barenblatt behavior
To raise the problem of infinite stresses on the face of fissure Barenblatt introduces forces of
cohesion enters the lips of the fissure. They are gravitational attractions being exerted between the particles of
leaves and other of the plan of separation of the fissure.
It is considered that the opening of the fissure costs an energy defined on the surface of discontinuity of
the fissure. It is called energy of surface. The problem is solved by minimizing the sum of this
energy and of elastic energy.
The surface of discontinuity is modelized by a finite element of gasket (Doc. [R3.06.09]) only element of
code supporting this law of behavior.
The condition of noninterpenetration of the lips of the fissure is taken into account by a method of
penalization.
To preserve a local processing of the conditions of opening, one carried out a regularization of
the energy of surface for jumps close to zero. With this intention a small parameter was introduced
(key word: SAUT_C) defining the width of these “small” jumps.
1.1
Energy of surface
On the surface of discontinuity
()
modelized by elements of the gasket type (
QUAD4
)
(see Doc. [R3.06.09]) an energy of surface is defined
S
E
depending on the standard of the jump on
displacement.
[]
[]
()
N
R
S
U
I
D
E
+
+
=
)
(U
The first term corresponds to energy necessary for the opening of the fissure, the second is one
energy of penalization which will make it possible to take into account the condition of noninterpenetration
lips of the fissure.
The function is defined:
[]
(
)
[]
[]
<
+
=
+
0
0
0
N
N
N
R
U
if
U
if
U
I
In other words, to want to impose the interpenetration of the lips will cost an infinite energy. In
practical, one will approach this function by a continuous function which tends quickly towards the infinite one
when the normal jump becomes negative:
[]
()
[]
-
+
D
U
C
U
I
N
N
R
2
2
1
~
with
C
coefficient of penalization and
-
.
negative part.
Note:
C
is adjusted starting from critical jump SAUT_C entered by the user, it is the same coefficient
that that taken for the regularization: cf slope of the stress in 0
+
.
Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
4/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
Moreover one poses
[]
[]
[]
2
2
T
N
U
U
U
+
=
+
so that a negative normal jump does not influence
the energy of surface and does not make evolve/move the threshold of the criterion.
One defines the energy of surface according to the cases of figure:
·
If
[]
U
[SAUT_C, +
[and
[]
0
N
U
:
[]
=
D
E
S
)
(
U
with
[]
[]
-
=
-
U
U
C
C
G
C
E
G
1
)
(
C
G
is the rate of critical refund of energy and
C
the stress criticizes in the beginning.
K
having the following properties (choice carried out starting from article [F&M]):
[]
()
=
concave
0
0
C
G
U
·
If
[]
U
[0, SAUT_C [and
[]
0
N
U
:
[]
=
D
E
S
)
(U
with
K
continuous and derivable in SAUT_C and such as
()
0
0
K'
=
:
[]
[]
constant
with
C
C
U
C
~
~
2
1
)
(
2
+
=
U
SAUT_C
G
C
K
[]
U
SAUT_C
K
G
C
[]
U
Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
5/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
·
If
[]
0
<
N
U
and
[]
0
=
T
U
:
Energy of penalization:
[]
-
=
D
U
C
E
N
S
2
2
1
1.2
Stress in the element of gasket
The stress in the element of gasket derives from the energy of surface. It is given by:
[]
[]
=
=
T
N
T
N
U
U
·
If
[]
U
[SAUT_C, +
[and
[]
0
N
U
:
[]
[]
[]
[]
[]
[]
=
-
-
U
U
U
U
C
C
C
C
G
T
C
G
N
C
E
U
E
U
·
If
[]
U
[0, SAUT_C [and
[]
0
N
U
:
[]
[]
=
-
-
C
JUMP
G
T
C
C
JUMP
G
N
C
C
C
C
C
E
U
E
U
_
_
SAUT_C
C
[]
U
SAUT_C
C
[]
U
K
G
C
[U]
N
Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
6/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
·
If
[]
0
<
N
U
and
[]
0
=
T
U
:
[]
=
0
N
U
C
Note:
The regularization of energy in zero makes it possible to define a stress in the element for one
null jump. Physically in other words, as soon as the stress in the element of gasket
will increase, a small jump will appear. The behavior of the Barenblatt type
will be carried out that when the standard of the jump in the element exceeds SAUT_C, before the gasket
comprise as a spring.
1.3 Variables
interns
The law of Barenblatt behavior at summer implemented with three internal variables.
The first
1
v
is used for piloting and for the discharges, it is a threshold which corresponds to
greater jump (in standard) ever reached (see: [§2] Piloting).
The second
2
v
allows to know if the element is in elastic mode (
[]
U
< SAUT_C) or polishing substance
(
[]
U
> SAUT_C).
The third corresponds to the percentage of energy of surface dissipated during the loading:
-
-
=
1
3
exp
1
v
G
v
C
C
SAUT_C
C
N
[]
N
U
Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
7/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
2 Piloting
The piloting recommended with the law of Barenblatt behavior is of type
PRED_ELAS
[R5.03.80], it
allows to follow dissipative solutions presenting of instabilities.
The implementation of this technique of piloting requires the resolution of the following local equation:
[] []
)
_
(
)
_
(
1
0
C
JUMP
L
C
JUMP
F
N
C
N
el
+
+
=
+
-
+
U
U
éq
2-1
with
unknown factor which one obtains by the solution of a quadratic equation,
C
L
length
characteristic of the model, and
N
threshold which will make it possible to keep in memory the standard of the jump with
the previous moment (stored in the variable interns No 1 of the law of Barenblatt behavior).
The threshold evolves/moves as follows:
If
[]
C
JUMP
N
_
1
+
+
U
then
-
=
+
If not
[]
C
JUMP
N
_
1
-
=
+
+
U
Piloting imposes ultimately that the jump of displacement continues to increase some share length
potential fissure.
The implementation of this piloting makes it possible to follow unstable branches of balance of the curve
total force/imposed displacements (see [Figure 2-a]).
Appear total Curved 2-a: Force/Déplacements imposed with unstable branch
Code_Aster
®
Version
7.2
Titrate:
Law of behavior of Barenblatt and piloting
Date:
12/03/04
Author (S):
J. LAVERNE
Key
:
R7.02.11-B
Page
:
8/8
Manual of Reference
R7.02 booklet: Breaking process
HT-66/04/002/A
3 Words
keys
The law of Barenblatt behavior is used in STAT_NON_LINE and DYNA_NON_LINE with the word
key BARENBLATT. This law of behavior is usable on elements of gasket with
modeling PLAN_FISSURE and AXIS_FISSURE.
Two parameters are to be seized in DEFI_MATERIAU:
SIGM_C
:
critical stress
SAUT_C
:
small parameter of regularization
Controls
STAT_NON_LINE COMP_INCR
RELATION
BARENBLATT
DYNA_NON_LINE COMP_INCR
RELATION
BARENBLATT
AFFE_MODELE MODELING PLAN_FISSURE
AXIS_FISSURE
DEFI_MATERIAU RUPT_FRAG
SIGM_C
SAUT_C
This law of behavior was validated by the case test SSNP118 (see Doc. [V6.03.118]).
4 Bibliography
[1]
Mr. CHARLOTTE, J.J. MARIGO, G. FRANKFURT, L. TRUSKINOVSKY: “Revisiting brittle
fracture ace year energy minimization problem: Comparisons off Griffith and Barenblatt surfaces
energy models ", Symposium one continuous ramming and fracture.
[2]
J. LAVERNE, E. LORENTZ, J.J. MARIGO: “A model of rupture with energy of
Barenblatt: theoretical considerations and numerical establishments. 16
ème
French congress
from Mechanics, Nice September 1-5, 2003.