Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
1/16
Organization (S): EDF/IMA/MN
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
R5.03.05 document
Viscoplastic relation of behavior of Taheri
Summary
One presents in this document the establishment of the relation of viscoplastic behavior of Taheri,
available for the whole of the isoparametric elements (continuous medium 2D and 3D) except for
plane constraints. After a presentation of the equations of evolution of this law, one describes the system obtained
by implicit discretization; it is shown in particular that it always admits a solution.
This model is well adapted to describe the response of the austenitic steels under cyclic stresses, and
in particular the phenomenon of progressive deformation. On the other hand, because of its complexity (two surfaces of
charge, semi-discrete internal variable), it does not appear desirable to employ it for applications
different (monotonous way of loading, for example).
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
2/16
1
Description of the model
The relation of behavior proposed by Taheri [bib5] makes it possible to describe the response of steels
austenitic under cyclic stresses: it is indeed well adapted to represent it
phenomenon of progressive deformation. Before stating the equations themselves, one can
to specify that this model differs from traditional plasticity (criterion of von Mises with work hardening
kinematics and isotropic) by two characteristics, sources of difficulties in the numerical formulation.
On the one hand, the evolution of the dissipative variables rests on two criteria of load instead of one:
first, traditional, conditions the appearance of plastic deformation, the second makes it possible to keep one
trace “maximum” work hardening reached by material to account for the phenomenon of ratchet.
In addition, satisfactorily to represent the progressive deformation, an internal variable
semi-discrete was introduced. Constant when the behavior is dissipative, it evolves/moves only in
the elastic mode of material. Of original appearance, this model does not rest any less on
physical bases, always exposed in Taheri [bib5]. It is accessible, in a wide version
viscoplastic (necessary to describe the behavior under high temperatures), by
order STAT_NON_LINE under key word RELATION: VISC_TAHERI.
1.1 Behavior
plastic
A detailed description of the law of behavior is given in Taheri and Al [bib6].
Briefly, the state of material is described by its state of deformation, its temperature like
four internal variables:

tensor of total deflection
T
temperature
p
cumulated plastic deformation
p
tensor of plastic deformation
p
constraint of peak, memory of maximum work hardening
p N
plastic tensor deformation due to the last discharge (variable semi-discrete).
The equations of state which express the thermodynamic forces associated according to the variables
of state are written:
= K Tr (- HT) + 2 µ (~ - p)
HT
Id
= (T - Tréf) Id
éq 1.1-1
p

has
has
- B p 1
R
R0 (2
p p
1
éq 1.1-2
3) With (
N)
m E
S
=
+
-





D
D
eq




= -
-

p
B p 1 -
X
[p p p
=
S -



S
C
N]
C =
C + C E
1
éq 1.1-3
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
3/16
~
has
deviatoric part of a tensor has
R
isotropic variable of work hardening
X
kinematic variable of work hardening
K, µ
modules of compressibility and shearing

thermal dilation coefficient
T ref.
temperature of reference
S
constraint of ratchet
B, , A,
has m, C, C

other characteristics of work hardening of material
1
Let us note that the moduli of elasticity and the thermal dilation coefficient are indicated by
the user by the command DEFI_MATERIAU, key word ELAS, while characteristics of
work hardening are fixed by key word TAHERI. These characteristics can depend on
temperature, by employing key words ELAS_FO and TAHERI_FO. Also let us specify that one
example of identification of the characteristics of work hardening on uniaxial situations is given
in Geyer [bib2].
The evolution of the internal variables is defined by two criteria. The first controls plasticity
traditional with work hardenings kinematics and isotropic compounds:
~
0
X
F = (~ - X) - R 0
and
S
=
-
éq 1.1-4
eq
(~ - X) eq
()
1
eq
equivalent standard: has
(a~:a~ 2
eq = 32
)
F
criterion of plasticity
s0
normal external with the criterion F
This criterion is matched traditional condition of load/discharge:
if F < 0 or s0 0
p = 0
(elasticity
! :
!
)

éq 1.1-5
if
F

= 0 and
s0 0
p 0 such as
>
!
!F = 0
(plasticity
! :
)
And the law of flow associated with the criterion F is:
3
2
!p =
! 0
p S
and thus
!p =
!p
éq 1.1-6
2
3 eq
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
4/16
The second criterion controls the evolution of the constraint of peak. Geometrically in the space of
diverters of the constraints, it translates the fact that the first surface of load (F = 0), represented by
a sphere of center X and radius R, remains inside a sphere of center the origin and
p. It is written simply:
G = X + R
p
eq
-
0
éq 1.1-7
G criterion of maximum work hardening
According to the preceding geometrical considerations, the evolution of the constraint of peak is:
if


G < 0 or
!X
p
eq +!
R 0
! = 0

éq 1.1-8
if
G = 0 and
!X
p
eq +!
R > 0! 0 such as!G = 0


It should be noticed that in the natural state of material, the constraint of peak is not null but is worth
initial elastic limit, namely:
p (init)
ial
= (1 - m) R0
Until now, we did not evoke the evolution of the semi-discrete internal variable p N. In fact,
it evolves/moves only in elastic mode. More exactly, this variable takes account of the state of
plastic deformation during the last discharge; in other words, at the beginning of each discharge, this
variable should take the value of the current plastic deformation instantaneously. However, for
to preserve a continuous behavior, one regularizes the evolution of p N in the following way:
In elastic mode:
p
p

p
p
p
if
elasticity classi
)
N =
! = 0
(
that
! =!
N
(N -)
éq 1.1-9
if p
p
(pseudo - déchar)
N
! 0 tq!F = 0
Ge
In plastic mode:
!pn = 0
The behavior is thus completely given. Before passing to the introduction of viscosity,
the observation of two surfaces of load calls an important remark. One could think that
surface G = 0 is actually activated only in plastic mode. In practice, it of it is nothing. One can
for example to quote the case of a thermal loading: a cooling involves (generally) one
dilation of the surface of load F = 0, so that the constraint of peak is brought to evolve/move for
to preserve G 0, and this same in elastic mode.
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
5/16
1.2
Taking into account of viscosity
To model the behavior of the stainless steels under cyclic loading when
temperature is about 550°C, it is not more possible to neglect the terms of creep. To return
count these effects of viscosity while preserving the properties of the preceding model, a method
simple consists in making viscous the evolution of the plastic deformation. In other words, viscosity
intervenes only in plastic mode: no direct influence on the semi-discrete internal variable nor on
the surface of load G = 0. For that, while following Lemaitre and Chaboche [bib3], one replaces
condition of coherence [éq 1.1-5] by:
NR

F

!p =

éq 1.2-1
K P1 M
F
positive part of F (hooks of Macauley)
K, NR, M
characteristics of viscosity of material
The characteristics of viscosity of material are indicated in command DEFI_MATERIAU,
either by key word LEMAITRE if they do not depend on the temperature, or by the key word
LEMAITRE_FO in the contrary case. In the absence of one of these key words, the behavior is
supposed plastic.
Unchanged all the other equations of the model are left. It will be seen that such an introduction of
viscosity involves only minor modifications of the implicit algorithm of integration of the law of
behavior.
1.3
Description of the internal variables calculated by Code_Aster
The internal variables calculated by Code_Aster are 9. They are arranged in
the following command:
1
p
cumulated plastic deformation
2
p
constraint of peak
3 to 8
p N
plastic tensor of deformation due to the last discharge
(arranged in order xx, yy, zz, xy, xz, yz)
9

loadmeter/discharge (cf [§2.3])
0 elastic discharge
1 traditional plastic load
2 plastic load on two surfaces
3 pseudo-discharge
As for the tensor of the viscoplastic deformations, it is not arranged among the internal variables but
can be calculated in postprocessing via the command CALC_ELEM, options
“EPSP_ELGA” or “EPSP_ELNO”, (cf [U4.61.02]).
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
6/16
2
Numerical formulation of the relation of behavior
In order to be able to treat within the same framework plasticity and viscoplasticity, one chooses to proceed to one
implicit discretization of the relations of behavior, (cf [R5-03-02]). Let us note moreover that one
explicit procedure of integration is delicate for two reasons: on the one hand, processing of the variable
semi-discrete is necessarily implicit and can lead to numerical oscillations (one
pseudo-discharge, therefore one solves F = 0, and of the blow, F can be (very weak but) higher than zero,
from where load with the step following instead of continuing the discharge), and in addition, the equation [éq 1.1-2] is not
not derivable when p = p N.
2.1
Implicit discretization of the equations of behavior
Henceforth, one adopts the convention of following notation. If U indicates a quantity, then:
U
quantity U at the beginning of the step of time
U increment of the quantity U during the step of time
U
quantity U at the end of the step of time (not of exhibitor +)
Let us start by introducing the elastic constraint, i.e. the constraint in the absence of increment
of plastic deformation. One can notice besides that only the term deviatoric cheek a role in
the nonlinear part of the behavior:
E =
(- HT) + µ (~-p-
K tr
2
) and ~=~e - 2µ p
Id
éq 2.1-1
“$ #
$ $ %
$

~ E
By taking account of the equations of state [éq 1.1-1] and [éq 1.1-3] and of the law of flow [éq 1.1-6], one a:
déf
~
~
-
3
S
X
E

C (p
p p
0
=
-
=
-
S
- N) - (2µ + CS) p
S
éq 2.1-2
2
By noting that is not other than S normalized, one deduces some immediately:

3
~
p -
S + (2 µ + C S) p
s0
E

C

-
p

éq 2.1-3


2


=
- (S
p
eq
N)
“#
$
%
$
Consequently, S is entirely determined by:
E
0
0
S
3
S = S S
S =
S
=
with
and
-
2 µ

éq 2.1-4
E
eq
eq
(+ C S) p
eq
S
2
eq
Finally, the functions of load are:
has

has
0
p -
p
3

F = S - D R + (2


0
S
3) With
- N +
p
eq



éq 2.1-5
2


eq



has
p
3

has
p
0
0
2
p -
3

G C S
p


S p

S
D
N
R
(

3) With
p
p 0
p
=
-
+
N
S
-
éq 2.1-6


2

+
+
-
+



2

eq

eq



Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
7/16
2.2
Taking into account of the viscous terms
In the absence of viscous terms, the relation of discretized coherence is:
Elastic mode:
F 0 and
p
= 0
éq 2.2-1
Plastic mode:
F = 0 and
p
0
On the other hand, in the presence of viscosity, the condition of coherence is replaced by the equation
[éq 1.2-1] which, discretized, is written:
1

NR
p

F

1 M p



F
= K p
NR
=
éq 2.2-2
T
K P1 M



T
In other words, while posing:
1
~
1 M
p

F = F -

K p
NR



éq 2.2-3
T

the viscoplastic increment of cumulated deformation is determined by:
~
Elastic mode:
F 0 and
p
= 0
~
éq 2.2-4
Viscoplastic mode:
F = 0 and
p
0
Finally, by adopting an implicit discretization, the only difference between the laws of behavior
plastic and viscoplastic lies in the form of the function of load F: a term there is observed
complementary in the event of viscosity. In fact, incremental plasticity seems the borderline case
(without associated numerical difficulty) of incremental viscoplasticity when viscosity K tends towards
zero. Let us note that this remark was already mentioned by Chaboche and Al [bib1].
2.3
Discretization of the conditions of coherence
Before discretizing the conditions of coherence and describing the various modes of behavior
possible, a remark is essential as for the processing of the semi-discrete variable. Like
only intervenes “to control” pn, one can always bring back itself during a step of time to:
p
p -
p -
N
= N + (1 -)
0 1
éq 2.3-1
The value of is then fixed by the conditions of coherence, which translates the equation of evolution
[éq 1.1-9] on the continuous level. Such a parameter setting with each step of time makes it possible to be freed from
storage of, in condition well-sure of preserving the values of pn.
After this opening remark, one can be interested in the conditions of coherence. For the criterion G
who controls the evolution of the constraint of peak, the discretized form of the condition of coherence is:
G (, p p,) 0
p 0
p G (, p p



,) = 0
éq 2.3-2
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
8/16
The condition of coherence relating to F is more delicate insofar as it controls the evolution
plastic deformation in plastic mode of load and evolution of in mode of discharge.
Once discretized, she is written:
In plastic mode of load (=)
1:
F (, p p, =) 1 = 0
p 0
p F (, p p

, =) 1 = 0
éq 2.3-3
In mode of discharge (p = 0):
F (p, 0 p
,) 0
0 1
F (p, 0 p
=
=

=
,) = 0
éq 2.3-4
To be able to select the mode of behavior of material, and thus equations to be solved,
first question is:
Are we in plastic situation or rubber band?
In fact, there is a solution in elastic mode (pseudo-discharge > 0 or traditional elasticity
= 0) if one can find an increment of constraint of peak such as:
Incremental condition of discharge (scalar equation out of p):
F (p =, 0 p
, =) 1 0
éq 2.3-5
G (p =, 0 p
, =) 1 0
p
0
p
G (p =, 0 p
, =) 1= 0
In the event of plastic load, i.e. when there is not p satisfactory [éq 2.3-5], one then has with
to solve the nonlinear system out of p and p following:
Plastic mode (nonlinear system out of p and p):
F (, p p
, =) 1 = 0
p 0
éq 2.3-6
G (, p p
, =) 1 0
p
0
p
G (, p p
, =) 1= 0
On the other hand, in elastic situation, two choices are still possible: pseudo-discharge (> 0) or
traditional elasticity (= 0). The second case being more favorable, one starts by examining whether it is not
not realizable, i.e. if there is an increment of constraint of peak such as:
Incremental condition of traditional elastic mode (scalar equation out of p):
F (p =, 0 p
, = 0) 0
éq 2.3-7
G (p =, 0 p
, = 0) 0
p
0
p
G (p =, 0 p
, = 0) = 0
Lastly, if it were to be a question of a discharge pseudo-rubber band, it remains to solve the nonlinear system in
and p following:
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
9/16
Discharge pseudo-rubber band (nonlinear system in and p):
F (p =, 0 p
,) = 0
0 1
éq 2.3-8
G (p =, 0 p
,) 0
p
0
p
G (p =, 0 p
,) = 0
Let us note as of now that nonlinear systems [éq 2.3-6] and [éq 2.3-8] can be reduced to
solution of a simple scalar equation if p = 0 makes it possible to obtain a solution.
One can summarize the algorithm of choice of the equations to be solved by the decision tree below.
Charge/discharge
(2.3-5)
Plastic resolution
Elasticity
(2.3-6)
traditional/pseudo
(2.3-7)
Resolution
Resolution
traditional elasticity
pseudo-discharge
(2.3-8)
Appear 2.3-Error! Argument of unknown switch. : Decision tree to choose the mode of
behavior
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
10/16
2.4
Framing of the solutions
With the reading of the preceding paragraph, one could note the need for solving (numerically) one
certain number of scalar equations or nonlinear systems. For that, it is always
interesting to have an interval on which to seek the solution. On the one hand, a framing of
solution shows its existence (what strongly reinforces the chances of success of an algorithm of
resolution!), and in addition, it allows a suitable digital processing, therefore surer.
Concerning p, one undervaluing is of course 0. In addition, the constraint of ratchet S represents
a limit beyond which the model does not have any more a smell. In fact, with the examination of constant materials
obtained by identification, cf Taheri and Al [bib6], G becomes indeed negative when p = S if
the difference between the plastic deformation and the plastic deformation due to the last discharge too is not
important (a few %):
G (, p p
=,
S) = (C
p
p
+ C1) (- N)
S
eq
“#
$ %
$
13
“$ #
$ $ %
$
3%


éq 2.4-1
0

+ (
R
A has
has
1 - m)
+ (2
p p
1
0
3)
(- N) -
“# %
S
S
eq


0,75
&

“# %” #
$
%
$
20%
0,5
0,7

One can also seek one raising for p. By examining the expression of F:
F (, p p
,) S
- D R0
eq

3
0
eq -
(2 µ + CS) p - D R
éq 2.4-2
2

3
- (2 µ + C S
-
0
) p -
(D
max
p) R
2
One deduces one from them raising for p, such as F (, p p
,) 0:

0
max (p
) - (
D p
p
,) R
p
3
éq 2.4-3
(2 µ + C S
)
2
~e -
-
p -
p

(CP p
,)
p

-

N
E
p

(S
)
smax ()
eq
= max
éq 2.4-4
~
E -
p -
C
p



(S
p
-
N) eq

Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
11/16
In particular, one can give one raising (coarse) p independent of p:
max
0
max - (1 - m) R
pmax = 3
éq 2.4-5
(2 µ + C S
)
2

p -
p

-
- p


N
E max
E
[S
]
S
eq
max = eq + (C1 + C) max [
éq 2.4-6
p -
S
- S pn

] eq
One can then notice that the systems [éq 2.3-6] and [éq 2.3-8] always admit a solution. In
effect, if for each system, one writes p respectively ()
p
p and
() solutions of
G = 0, then one a:
· Nonlinear system of plastic load:
F (p =, p
0 (p = 0), =)
1 0 and F (p, p
max (pmax), =)
1 0 éq 2.4-7
· Nonlinear system of pseudo-discharge:
F (p,
p
= 0 (=)
1, =)
1 0
and
F (p =, p



0 (= 0), = 0) 0 éq 2.4-8
3
Methods of numerical resolution
The resolution of the incremental equations confronts us either with a nonlinear scalar equation, or
with a nonlinear system with two unknown factors. The numerical methods below are exposed
employees. One also examines the calculation of the tangent matrix, possibly used by
the total algorithm of STAT_NON_LINE, (cf [R5.03.01]).
3.1
Scalar equation: method of secants
It is a question of solving a nonlinear scalar equation by seeking the solution in an interval of
confidence. For that, one proposes to couple a method of secant with a control of
the interval of search. That is to say the following equation to solve:
F (X) = 0
X [
,
has B] F (A) < 0 F (b) > 0
éq 3.1-1
The method of the secant consists in building a succession of points X N which converges towards the solution.
It is defined by recurrence (linear approximation of the function by its cord):
N
n-1
n+1
n-1
N 1
X - X
X
= X
- F (X -)
éq 3.1-2
F (xn) - F (xn-1)
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
12/16
In addition, if X n+1 were to leave the interval, then one replaces it by the terminal of the interval in
question:
if

X n+1 < has
then
X n+1:= has

éq 3.1-3
if
X n+

1 > B
then
X n+1:= B
On the other hand, if X n+1 is in the interval running, then the interval is reactualized:
if

X n+

1
[
has, B] and F (xn+1) < 0
then
a:= X n+1

éq 3.1-4
if
X n+1 [
has, B] and
1
1


F (xn+) > 0
then
b:= X n+
One considers to have converged when F is sufficiently close to 0 (tolerance to be informed). As
at the first two leader characters, one can choose the terminals of the interval, or, if one lays out
of an estimate of the solution, one can adopt this estimate and one of the terminals of the interval.
3.2
Nonlinear systems: method of Newton and linear search
One presents here a method of Newton which one associated a linear technique of search and
a control of the direction of descent not to leave the field of search (terminals on
unknown factors).
That is to say the following nonlinear system:
F (X, y) = 0
xmin X xmax

with

éq 3.2-1
G
(X, y) = 0
y
min y ymax
If (X, y) is a point of the field of search, then one builds a succession of points (xn yn
,
) which
converge towards a solution (or at least, it is hoped for) by the following process.
· Determination of the direction of descent
A direction of descent (X, y) is given by the resolution of the linear system 2 X 2:
Fn Fn X
N
, X
, y

F




éq 3.2-2
Gn
Gn y = -
N
, X
, y
G




Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
13/16
· Correction of the direction of descent
The direction of descent (X, y) is corrected so that the points considered are
in the field of search (with max the maximum length which one is authorized to describe it
length of the direction of descent):

xmin - X
xmax - X
if X +
X
max < X
X
min
:=
if X +
X
max > X
X
max
:=


max


max
y
éq 3.2-3
min - y
ymax - y
if
y +
y
max < y
y
min
:=
so y +
y
max > y
y
max
:=


max
max
· Seek
linear
It any more but does not remain to minimize the quantity E = (F2 + G2)/2 in the direction of descent.
Let us note that the standard E that one minimizes thus is a measurement of the error made in
resolution of the system: it is null when (X, y) is solution of the system [éq 3.2-1]. For
to minimize E, one simply will seek to cancel his derivative, i.e. to solve
the scalar equation:
[E (X + X, y + y)] = 0 and 0
éq 3.2-4

max
“#
$ %
$
(FF +GG
, X
, X) X + (F F +G G
, y
, y) y
· Criterion of convergence
One considers to have converged when the error E is lower than a prescribed size. By
elsewhere, if the standard of the direction of descent becomes too weak (another size with
to inform), one can think that the algorithm does not manage to converge.
3.3 Criteria
of stop
Until now, values of stop and iteration counts maximum of the methods of resolution
the preceding ones were not specified. Two cases should be distinguished.
· When one seeks to check the conditions of coherence (scalar equation or system not
linear according to the situation), one awaits precise results, whose relative tolerance is
fixed by the user in command STAT_NON_LINE under key word RESI_INTE_RELA,
(cf [U4.32.01]). According to whether one seeks to solve F =,
0 G = 0 or simultaneously
F = G = 0, the criterion of stop is expressed respectively:
F
G
F2 G2
1

or

+

or

R0
R0
R0
2
R0 limit elastic initial, provided by the user, cf [§ 1.1].
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
14/16
In addition, the user always specifies a maximum iteration count in the command
STAT_NON_LINE under key word ITER_INTE_MAXI, (cf [U4.32.01]).
· When one carries out iterations of linear search, one seeks to obtain a convergence
faster (or at least sourer). One should not therefore devoting an excessive time to it.
This is why one fixed once for a a whole iteration count maximum equal to 3, one
limit maximum max equal to 2 and one criterion of relative stop of 1%:
[E (X + X, y + y)]
-

10 2
[E (X + X, y + y)]



=0
3.4 Stamp
tangent
In the optics of a resolution of the equilibrium equations (total) by a method of Newton, it is
essential to determine the consistent matrix of the tangent behavior, (cf Simo and Al [bib4]).
This matrix is composed classically of an elastic contribution and a plastic contribution:

E
p
=
- 2 µ
éq 3.4-1



One immediately deduces from it that in elastic mode (traditional or pseudo-discharge), the matrix
tangent is reduced to the elastic matrix:
Elastic mode:

E
=
éq 3.4-2


On the other hand, in plastic mode, the variation of the plastic deformation is not null any more. Rules
of made up derivation allow to obtain:
p
0
3
p
3
p
p
3

0

S
0


0
0
=
S

+ p


S
Id
S
S éq 3.4-3

=

+
-

~e
2
~e
~e
2
~e

2



eq


tensorial product
~
One can note that one preferred to derive compared to E, knowing that one a:
p
p ~e
p
S S
=
= 2 µ P
with
P
éq 3.4-4

~.
E

~.
:
E
'~
S spaces symmetrical tensors
P projector on the diverters
Finally, it does not remain any more that to calculate the variation of p. Pour that, it is necessary to distinguish if it is about one
traditional mode of plasticity (p = 0) or of plasticity on two surfaces. As follows:
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
15/16
Traditional plasticity: F (, ~
p E) = 0
F E, ~ E
E
E
E
p
p
F
, ~

, p (, ~
p) p

= - F E p



= -
éq 3.4-5
, ~
(, ~) ~
()
~
E
F
E
, p (, ~
p)
Plasticity on two surfaces: F (, p, ~e) 0 and G (, p

, ~
p
p
E
=
) = 0
F p
F
,
, ~
E
F
F p p
F
G
E
p
G
, p
E
,

, ~

E
p
,
, ~








=
G
G

p =
-
éq 3.4-6
E
p
G
. ~
~
, p


E
F
F

,


, ~

, p

, p

G
G
, p
, p

An attentive examination of the expressions [éq 2.1-5] and [éq 2.1-6] makes it possible to note that the variations of
F
~
and G compared to E are not necessarily colinéaires with S °. By taking account of
[éq 3.4-3], one from of deduces whereas the tangent matrix is in general not symmetrical in mode
plastic. Rather than to impose the use of a nonsymmetrical solvor, much more expensive in time
calculation, one prefers to symmetrize this matrix.
3.5 Constraints
plane
The processing of the plane constraints adds a nonlinear equation to solve, coupled with
systems [éq 2.3-6] and [éq 2.3-8], (cf [R5-03-02]). In front of this considerable difficulty and the absence
of expressed need, one preferred not to make it possible to force a state of plane stresses with
level of the law of behavior. In other words, modeling C_PLAN is not available for
law of behavior VISC_TAHERI.
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Code_Aster ®
Version
4.0
Titrate:
Viscoplastic relation of behavior of Taheri
Date:
08/09/97
Author (S):
E. LORENTZ
Key:
R5.03.05-A
Page:
16/16
4 Bibliography
[1]
J.L. Chaboche, G. Cailletaud, 1996. Constitutive integration methods for complex plastic
equations. Comp. Meth. Appl. Mech. Eng., n° 133, pp. 125-155
[2]
PH. Geyer, 1992. Study of the elastoplastic model of behavior for the loading
cyclic develops with EDF/DER/MMN. Note intern EDF/DER, HT-26/92/39/A.
[3]
J. Lemaitre, J.L. Chaboche, 1988. Mechanics of solid materials. Dunod.
[4]
J.C. Simo, R.L. Taylor, 1985. Consist tangent operators for misses-independent
elastoplasticity. Comp. Meth. Appl. Mech. Eng., n° 48, pp. 101-118.
[5]
S. Taheri, 1990. A uniaxial law of behavior in elastoplasticity for the loading
cyclic. Note intern EDF/DER, HI-71/6812.
[6]
S. Taheri, pH. Geyer, J.M. Proix, 1995. Constitutive Three dimensional elastic-plastic law for
the description off ratchetting off 316 stainless steel. Note intern EDF/DER, HI-74/95/012/0.
Handbook of Référence
R5.03 booklet: Nonlinear mechanics
HI-74/97/019/A

Outline document