Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 1/14
Organization (S): EDF-R & D/AMA
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
Document: R7.01.14
Law of behavior ENDO_ISOT_BETON
Summary:
This documentation presents the theoretical writing and the numerical integration of the law of behavior
ENDO_ISOT_BETON which describes an asymmetrical local damage mechanism of the concretes, with effect of
restoration of rigidity. In addition to the local model, the nonlocal formulation with regularized deformation is also
supported to control the phenomena of localization.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 2/14
Count
matters
1 Introduction Domaine of application .................................................................................................... 3
2 local Law of behavior ................................................................................................................... 3
2.1 Theoretical writing ............................................................................................................................ 3
2.2 Taking into account of the withdrawal and the temperature .............................................................................. 6
2.3 Identification of the parameters .......................................................................................................... 6
2.3.1 Elastic parameters ............................................................................................................ 6
2.3.2 Parameters of damage ............................................................................................. 6
2.3.2.1 Use without dependence with containment ............................................................ 6
2.3.2.2 Use with dependence with containment ............................................................ 6
2.3.2.3 Passage of the values “user” to the values “models” ..................................... 7
2.4 Numerical integration ...................................................................................................................... 7
2.4.1 Evaluation of the damage ........................................................................................... 7
2.4.2 Calculation of the tangent matrix ................................................................................................. 8
2.4.2.1 Stamp tangent with constant damage ......................................................... 8
2.4.2.2 Term of the tangent matrix due to the evolution of the damage ....................... 9
2.4.3 Case of material completely damaged ...................................................................... 10
2.5 Description of the variables intern ................................................................................................. 10
3 Formulation with regularized deformation .......................................................................................... 11
3.1 Formulation .................................................................................................................................... 11
3.2 Integration of the law of behavior ........................................................................................... 11
3.3 Variables intern .......................................................................................................................... 11
4 Control by elastic prediction .......................................................................................................... 12
5 Bibliography ........................................................................................................................................ 12
Appendix 1
Demonstration of the clean reference mark of constraint ........................................................ 13
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 3/14
1
Introduction Domaine of application
The law of behavior ENDO_ISOT_BETON aims at modelling the most simply possible one
behavior of elastic concrete fragile. It can be seen like an extension of the law
ENDO_FRAGILE [R5.03.18] (with which it keeps a proximity of unquestionable formulation) for
applications of Génie Civil.
As for law ENDO_FRAGILE, material is isotropic. Rigidity can decrease, the loss of
rigidity measured by an evolving/moving scalar of 0 (healthy material) to 1 (completely damaged material).
On the other hand, contrary to ENDO_FRAGILE, the loss of rigidity distinguishes traction from
compression, to privilege the damage in traction. Moreover this loss of rigidity can
to disappear by return in compression, it is a question of the phenomenon of restoration of rigidity with
refermeture. It as should be noted as this law of damage aims at describing the rupture of the concrete in
traction; it is not thus adapted at all to the description of the nonlinear behavior of the concrete in
compression. It thus supposes that the concrete remains in a moderate compactness.
Law ENDO_ISOT_BETON present of softening, which generally involves a loss
of ellipticity of the equations of the problem and consequently a localization of the deformations, from where one
pathological dependence with the grid. To mitigate this deficiency of the model, a formulation not
local must be adopted: for law ENDO_ISOT_BETON, modeling GRAD_EPSI [R5.04.02],
based on the regularization of the deformation is usable. In this formulation, it should be noted that
only the relations of behavior are faded compared to a traditional local modeling;
consequently, the constraints preserve their usual direction.
Lastly, that one activates or not the nonlocal formulation, 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 control of the essential loading in statics.
control of the type PRED_ELAS [R5.03.80] then seems the mode of control of the loading it
more adapted.
2
Local law of behavior
2.1 Writing
theoretical
If one seeks to take account of the effect of refermeture, it is necessary to pay a great attention to
continuity of the constraints according to the deformations (what is an essential condition for
a law of behavior in a computation software by finite elements), cf [bib1]. Indeed, if one
model this effect in a too simplistic way, the law of behavior is likely great to present
a discontinuous answer.
To take account of the refermeture (i.e the transition between traction and compression), it is necessary to start
by finely describing what one calls traction and compression, knowing that in traction (resp.
compression) the fissure will be considered “open” (resp. “closed”). A natural solution is of
to place itself in a clean reference mark of deformation. In such a reference mark, the elastic free energy is written (
and µ indicating the coefficients of Lamé):
() =
() 2
tr
+ µ 2
éq 2.1-1
I
2
I
One can then define:
· a traction or voluminal compression, according to the sign of tr,
· a traction or compression in each clean direction, according to the sign of I.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 4/14
According to the rather reasonable principle according to - in a case of traction (“open fissure”), one corrects
the elastic energy of a factor of damage; in a case of compression (“closed fissure”),
one keeps the expression of elastic energy -, the free energy endommageable is written:
(
2
1
2
1
, D) = (
tr)
- D
D
H (-
tr) +
H (
tr)
+
µ H
H
éq 2.1-2
I
(- I) -
+
(I)
2
1+ D
1
I
+ D
It is noticed that the free energy is continuous with each change of mode. It is even
continuously derivable compared to the deformations, since it is nap of derivable functions (
function X 2 H (X) is derivable) and the continuity of the partial derivative at the points tr = 0 and = 0
I
is
immediate. One then clarifies the constraints (by knowing that they will be functions everywhere
continuous of the deformations). As in elasticity, the clean reference mark of the constraints coincides with
clean reference mark of the deformations, result shown in appendix.
One writes the constraints in the clean reference mark:
1 D
1 D
tr
H
tr
H
tr
2µ H
H
éq
2.1-3
II =
()
-
(-
) +
()
+
II
(- II) -
+
(II)
1+ D
1+ D
In this form, the continuity of the constraints with respect to the deformations is clear. The figure opposite
show the constraint
, 1
1 in the plan (
2) with constant damage (case 2D, deformation
plane). The effect of refermeture as well as the continuity of the constraints are quite visible.
Appear 2-a: illustration of continuity
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 5/14
The thermodynamic force F D associated with the variable interns damage is written:
1
D
+
F = -
=
2
tr
tr
éq
2.1-4
2
() H () + µ 2H
I
(I)
D
(1+ D
)
2
I
It remains to define the evolution of the damage. The diagram selected is that of the standard models
generalized. A criterion should be defined, that one takes in the form:
F (Fd) Fd
=
(, D) - K éq
2.1-5
where K defines the threshold of damage. In order to take into account, on the level of the evolution of
the damage, the effect of containment, the threshold K depends on the state of deformation, in the form:
K = K + K
H -
éq
2.1-6
0
(
tr
1
) (
tr)
One compels oneself to remain in the field:
F (Fd) 0 éq 2.1-7
The evolution of the variable of damage is then determined by the conditions of Kuhn-Tucker:
D & = 0 for F < 0
éq 2.1-8
D & 0 for F = 0
Note:
From a formal point of view, the generalized standard materials are characterized by one
potential of dissipation function positively homogeneous of degree 1, transformed
Legendre-Fenchel of the indicating function of the field of reversibility, which is thus worth here:
(&d) = sup Fd &d = K &d + I
éq
2.1-9
IR+
D
D
(&d)
F/F (F) 0
One will note the presence of an indicating function relating to &
D, which ensures that the damage is
growing.
It still remains to take into account the fact that the damage is raised by 1. From a point of view
intuitive, that seems easy. To keep a writing completely compatible with the formalism
generalized standard, it is enough to introduce an indicating function of the acceptable field into
the expression of the free energy:
2
1
(
D
, D)
= (tr) H (- tr) + - H (T
R) +
2
1 + D
éq
2.1-10
µ2
1 D
H (- + -
I)
H (I) + I
D
I
] -; ]
1 ()
1 +
I
D
The introduction of this indicating function prevents the damage from exceeding 1, indeed, for
D = 1, F D = -
= -, and the damage does not evolve/move any more.
D
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 6/14
2.2
Taking into account of the withdrawal and the temperature
The law of behavior takes into account a possible withdrawal of dessication, a possible withdrawal
endogenous and a possible thermal deformation. The deformation in question in it
document being then “elastic strain” ~
HT
rd
Re
= - - -.
On the other hand, the parameters materials in question in the next paragraph are considered
like constants (in particular, they cannot depend on the temperature, in the state of
current development)
2.3
Identification of the parameters
The parameters of the law of behavior are 4 or 5 (see following paragraphs). They
are classically provided in operator DEFI_MATERIAU.
2.3.1 Parameters
rubber bands
They are simplest: it is about the two traditional parameters, Young modulus and coefficient of
Poisson, provided under key word ELAS or ELAS_FO of DEFI_MATERIAU.
2.3.2 Parameters
of damage
According to whether the user wants to use the dependence of the threshold with containment or not, 2 should be provided
or 3 parameters to control the law of damage.
2.3.2.1 Use without dependence with containment
In this case, one considers that the parameter K is null. It should be noted that the compactness of
1
concrete must remain moderate so that the law remains valid (compressive stress about
some times the constraint with the peak of traction, in absolute value).
The user must inform, under key word BETON_ECRO_LINE of DEFI_MATERIAU, the values of:
· SYT: limit of simple tensile stress,
· D_SIGM_EPSI: slope of the curve post-peak in traction.
2.3.2.2 Use with dependence with containment
In this case, the dependence with containment makes it possible the concrete to keep a realistic behavior in
compression until the order of magnitude of appearance of nonthe linearity in compression, given by
SYC, Cf. below (classically, a compressive stress of about ten times the constraint
with the peak of traction, in absolute value)
The user must inform, under key word BETON_ECRO_LINE of DEFI_MATERIAU, the values of:
· SYT: limit of simple tensile stress,
· SYC: limit of the simple compressive stress,
· D_SIGM_EPSI: slope of the curve post-peak in traction.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 7/14
2.3.2.3 Passage of the values “user” to the values “models”
For information, one obtains the values of, K and possibly K (if the user informed
0
1
SYC) by the following formulas:
E
= -
D _ SIGM _ EPSI
2 1 +
1+ -
2 2
K
SYT
0 = (
)
2nd
1+
(1+) 2
E
K = SYC
-
1
(1+) (1 -
2) k0 (1 -
2) (SYC)
2.4 Integration
numerical
Two points are to be regulated before establishing the model: the first relates to the evaluation of
the damage; the second consists in calculating the tangent matrix, calculation made a little more delicate
that usually by the passage in a clean reference mark of deformation.
One places oneself here within the framework of the implicit integration of the laws of behavior. Dependence of
criterion according to containment [éq 2.1-6] is taken into account in explicit form, i.e the threshold K
is entirely determined by the state of deformation of the preceding step, this to simplify integration
model.
2.4.1 Evaluation of the damage
As one will see it, a simple scalar equation makes it possible to obtain the damage, which allows
to avoid a recourse to the iterative methods.
One notes D - the damage with the preceding step and D + the evaluation of the damage to the step
running to the current iteration which will be the damage with the current step when convergence is
attack. Simplest to evaluate the damage of the current iteration is to suppose that one
reached the criterion at the current moment, which results in:
1
2
2
F (Fd) = 0
+
µ
éq
2.4.1-1
2
(tr) H (tr) +
H
I (I) K
(
=
1+ D
) 2
I
what gives:
1
1
2
D test =
+
(tr) H (tr) +
2
µ H
1 éq
2.4.1-2
I
(I) -
K
2
I
3 cases arise:
·
D test D -: that wants to say that at the moment running, the criterion is not reached, one concludes from it that
D +
D -
=
,
·
D - D test 1: the criterion is thus reached, the condition of coherence implies D + = D test,
·
D test 1: the material is then ruined in this point, from where D + = 1.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 8/14
2.4.2 Calculation of the tangent matrix
The tangent matrix is the sum of two terms, the first expressing the relation
constraint/deformation with constant damage, the second being resulting from the condition F = 0. In
effect, one can write:
ij ij
ij D
=
+
éq
2.4.2-1
D
kl
kl
you
=
D C
kl F =0
If the user asks for calculation with tangent matrix (Cf. documentation of STAT_NON_LINE,
[U4.51.03]), the law of behavior provides the expression given by [éq 2.4.2-1]. On the other hand, if
the user asks for calculation with the matrix of discharge, the law of behavior provides the matrix
secant, i.e. the first term of the member of straight line of [éq 2.4.2-1].
2.4.2.1 Stamp tangent with constant damage
As we underlined previously, the calculation of the tangent matrix is a little delicate
fact of the writing of the model in the clean reference mark of deformation. Thus, one knows easily
stamp tangent with constant damage in the clean reference mark of deformation, but what one
seek is this same tangent matrix in the total reference mark.
If the damage does not evolve/move, in the clean reference mark of deformation, the matrix
sought a simple relation of degraded elasticity expresses:
~
1 D
1 D
I
-
~
= H (- tr) +
H (tr)
+
2µ H
H
éq
2.4.2.1-1
ij
(- J) -
+
(J)
1
D
1
D
J
you
+
+
D =c
It is now necessary to express the passage of the total reference mark to the clean reference mark of deformations, at least
if the eigenvalues of deformation are different. The tangent matrix not being
necessary that the algorithms of numerical resolution (diagram of Newton), one will allow oneself, at the time of
calculation of the tangent matrix (and only in this case) to disturb the possible ones numerically
identical eigenvalues (in order to make them distinct). It will be noticed in particular that that
allows, null damage, to find the matrix of elastic rigidity.
One notes with a tilde the tensors in the clean reference mark of deformation (which, one points out it, is too
the clean reference mark of constraints). By definition, by noting U the clean vector associated the i-ème
I
eigenvalue, the matrix basic change Q = (U U U, one a:
1
2
3)
~ T
= Q Q
= Q Q ~ + Q Q ~ + Q Q ~
ij
im
jm
m
im
jm
m
im
jm
m
If the eigenvalues of deformation are distinct, evolution of the clean vectors and
eigenvalues is given by (Cf. previously [§2]):
~&jk
U & U =
for J K éq
2.4.2.1-2
J
K
~ - ~
J
K
&~ = &~ for J K éq
2.4.2.1-3
I
II
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 9/14
One deduces Q easily from it:
~
~
Q =
jk
U
Q éq
2.4.2.1-4
ij
~
~ (K) =
jk
I
~
~ ik
K J
J -
K
K J
J -
K
While using then (the last expression being used only to obtain a clearly symmetrical matrix)
:
=
1
~
Q Q =
Q Q + Q Q
ij
ki
lj kl
(ki lj Li kj) kl
2
One thus obtains:
~
m
= Q Q
~ + Q Q ~ + Q Q ~
ij
im
jm
~
N
im
jm
m
im
jm
m
,
m N
N
m
~
1
Q Q Q Q
m
in
km
ln
jm
= Q Q Q Q
+
ln
~
im
jm
kN
~
kl
~
~
m
kl
-
,
2,
m N K L
N
K L
m
N
Nm
1
Q Q Q Qln
1
Q Q Q Q
1
Q Q Q Q
im
jn
km
in
lm
kN
jm
im
jn
lm
kN
+
~ +
~ +
~
~
~
m
kl
~
~
m
kl
~
~
m
kl
2
-
-
-
,
2,
2,
K L
m
N
K L
m
N
K L
m
N
Nm
Nm
Nm
éq 2.4.2.1-5
The tangent matrix with constant damage is thus written:
~
1
Q Q
Q Q
Q Q
Q Q
ij
m
(km ln + lm kN) (in jm + jn im)
With
Q Q Q Q
~
ijkl =
= im jm kN ln ~
+
~
~
2
kl
m, N
N
m
you
m, N
N -
D =c
m
Nm
éq 2.4.2.1-6
2.4.2.2 Term of the tangent matrix due to the evolution of the damage
The expression to be evaluated is written:
ij D
éq
2.4.2.2-1
D kl F =0
One writes the equation [éq 2.4.1-1] in the form:
1 +
W = K
éq
2.4.2.2-2
2 [()]
(1+d)
with: W () = (tr) 2 H (tr) + µ 2 H.
I
(I)
2
I
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 10/14
While differentiating this expression, it comes:
- 2 (1+) W
D
el
éq
2.4.2.2-3
3
() + 1+
= 0
(1+ D)
(1+ D) 2
W
with: el
=
One uses then the following equality:
1 el
=
+
éq
2.4.2.2-4
D (1+ D) 2
One concludes:
ij
1
=
+
el el
ij
éq
2.4.2.2-5
kl
2 1 + D W
kl
F =0
() ()
2.4.3 Case of completely damaged material
In the case of the completely damaged material, D = 1, the rigidity of the material point can
to cancel itself. That poses problem for the constraint by no means; on the other hand, that can involve
null pivots in the matrix of rigidity. To mitigate this difficulty, one allows oneself to define a rigidity
minimal, for the tangent matrix or the matrix of discharge. This minimal rigidity does not affect
value of the damage (which can reach 1) or the constraint (which can reach 0).
To preserve a reasonable conditioning of the matrix of rigidity, minimal rigidity is taken with
5
10 - rigidity initale. An indicator specifies the behavior during the step of current time:
·
= 0: no evolution of the damage during the step,
· = 1: evolution of the damage during the step,
·
= 2: damage saturated D = 1.
2.5
Description of the internal variables
The model has two internal variables:
· VI (1): damage D,
· VI (2): indicator.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 11/14
3
Formulation with regularized deformation
3.1 Formulation
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 which intervenes in the forced relation deformation and the deformation
regularized which controls the evolution of the damage. This one results from the local deformation
by resolution of the system of partial derivative equations according to:
- 2
B
L = 0
structure
in
éq
3.1-1
N = 0
on
normal
of
edge
N
where the characteristic length L is again indicated under the key word
B
LONG_CARA of
DEFI_MATERIAU. Finally, the relation of behavior is written in the following way, the equation
(2-3) remains identical, while the equation [éq 2.1-5] takes into account the regularized deformation:
1 D
1 D
tr
H
tr
H
tr
2µ H
H
éq
3.1-2
II =
()
-
(-
) +
()
+
II
(- II) -
+
(II)
1+ D
1+ D
F (F D) = F D (, D) - K
éq 3.1-3
3.2
Integration of the law of behavior
One of the advanced advantages for the nonlocal formulation with regularized deformation is the little of
modifications which it involves in the construction of the law of behavior. Indeed, the integration of
internal variables is completely controlled by the regularized deformation.
The method of integration is exactly that described in the paragraph [§2.4.1], in condition of
to replace the deformation by the deformation regularized in the equations.
For the calculation of the tangent matrix, the expressions are the same ones as those given to
paragraph [§2.4.2], certain expressions of the deformation are to be replaced by the deformation
regularized (when the deformation concerns the criterion), while others do not change (when
deformation raises of the forced relation deformation).
3.3 Variables
interns
They are the same internal variables as for the local law:
· VI (1) damage D,
· VI (2) indicating.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 12/14
4
Control by elastic prediction
The control 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. Consequently, only them
points where the damage is not saturated will be taken into account. The algorithm which deals with
this mode of control, 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 4-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:
U = U + U
éq 4-2
0
1
where U and U are given. Thanks to the linearity in small deformations of the operators deformation
0
1
(calculation of the deformations starting from displacements) and regularized deformation, one also obtains
following decompositions:
= 0 + 1
and
= 0 + 1
éq
4-3
The function el
F presenting the good property to be convex, the equation [éq 4-1] presents zero, one or
two solutions, which are required as follows:
· Determination of the number of solutions per study at the boundaries ± and possibly (if
value at the two boundaries is each time positive) determination if el
F presents a minimum
negative;
· Determination of a framing of each solution starting from the preceding study
· Determination of the solution (for a convex function knowing the framing, this
search is simple and fast)
5 Bibliography
[1]
P.B. BADEL: Contributions to the digital simulation of structures out of reinforced concrete. Thesis of
Université Paris VI, 2001.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 13/14
Appendix 1 Démonstration of the clean reference mark of constraint
The term traces some in energy does not pose a problem: it is invariant by any change of reference mark.
2
1 D
Remain the term out of H -
+ - H
I
(I)
(I).
1 + D
I
Notation: one writes with an index (for example) the i-ème eigenvalue of a tensor which is written (while clarifying
I
its two indices).
kl
· If the eigenvalues of the deformation are all distinct, it is shown whereas & =
I
&ii, with
&kl components of & in the fixed reference mark coinciding with the clean reference mark of deformation with
the moment considered (in this reference mark one thus has =).
kl
K
kl
Indeed, let us write the deformations in the form:
= U U
I
I I
I
While differentiating this expression, it comes:
& = & U U
U
U
U
U
I
I
I
+ &
I
I
I + I
I & I
I
By using the fact that the clean vectors are orthonormés:
U U = U & U + U U & = 0
I
J
ij
I
J
I
J
one obtains the variations of the eigenvalues and the clean vectors:
& =
jk
=
I
&ii and &
&
U U
J
K
for J
K
-
J
K
This is obviously valid only if the eigenvalues are distinct (as one can see it
clearly on the expression of the variations of the clean vectors). That comes owing to the fact that them
clean vectors are not continuous functions of the elements of the matrix.
· If two eigenvalues of deformations are equal (and apart from the case very
private individual where they are also null), they are either positive, or negative. Let us take the case
where they are positive (the other case lends itself to a demonstration in any similar point). Energy
3
concerning these two eigenvalues is written then: 2
(the two equal eigenvalues are
I
i=2
considered to have indices 2 and 3). By differentiating this expression, one obtains:
23
3
D
D by noting the common eigenvalue.
I
I = 2
I
i=2
i=2
By invariance of the trace of a matrix, here the restriction of the deformation on the clean plan
considered, one obtains:
3
3
D
D, whatever the evolution which underwent the clean reference mark at this time.
I =
II
i=2
i=2
For the remaining eigenvalue (distinct from both others and index 1 with the selected notations),
one a: D = D.
1
11
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Code_Aster ®
Version
7.4
Titrate:
Law of behavior ENDO_ISOT_BETON
Date:
26/05/05
Author (S):
P. Key BADEL
:
R7.01.04-B Page
: 14/14
By gathering these expressions, one obtains:
D 3 2
H
D H
D
H
H D
I
(I)
3
3
= (2
2
1
()
+
2
1
I
(I)
= II (II)
II
i=1
i=2
i=1
In conclusion, that the eigenvalues are distinct or not, one obtains:
D
2 H () = 2
H
II
(II) D
I
I
II with the adopted notations.
I
I
This reasoning spreads easily with the case of three equal eigenvalues.
The differential of energy with constant damage is written then:
1 D
D (, D)
=
- + -
+
D =cte
(tr) D (tr) H (tr)
H (tr)
1 + D
1 D
2µ
D H - + -
II
II
(II)
H (II)
1 + D
I
On this expression, one observes well that the clean reference mark of deformation is also reference mark
clean of constraint.
Handbook of Référence
R7.01 booklet: Modelings for Génie Civil and the géomatériaux ones
HT-66/05/002/A
Outline document