background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
1/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Organization (S):
EDF-R & D/AMA















Manual of Reference
R5.03 booklet: Nonlinear mechanics
Document: R5.03.01



Quasi-static nonlinear algorithm
(operator STAT_NON_LINE)




Summary:

The operator
STAT_NON_LINE
[U4.51.03] of Code_Aster allows in the case of a quasi-static stress
to integrate various types of non-linearities coming from the behavior of material or of large
geometrical transformations. One describes the total algorithm of resolution here employed. The integration of the relations
of behavior themselves is described in other documents, like [R5.03.02] for the élasto-
plasticity, to which one will be able to refer for examples. For calculations in great transformations
geometrical, one will be able to consult for example the document [R5.03.20] on nonlinear elasticity into large
displacements, or the document [R5.03.21] on the thermoelastoplasticity with isotropic work hardening.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
2/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Count
matters
1
Presentation ........................................................................................................................................... 3
1.1
General information on the relations of behavior ............................................................................... 3
1.2
Position of the nonlinear quasi-static problem ............................................................................ 4
2
Method of resolution ........................................................................................................................... 6
2.1
Principle of the method of NEWTON .............................................................................................. 6
2.2
Adaptation of the method of NEWTON to the problem arising ........................................................... 7
2.2.1
Resolution without boundary conditions dualized .................................................................. 7
2.2.2
Resolution with boundary conditions dualized ........................................................... 9
2.2.2.1
Phase of prediction ................................................................................................... 9
2.2.2.2
Loop on the iterations of NEWTON ..................................................................... 11
2.2.2.3
RIGI_MECA_TANG
and
FULL_MECA
.......................................................................... 16
2.2.3
Case of the following loadings ............................................................................................ 16
2.3
Linear search ......................................................................................................................... 17
2.3.1
Principle ................................................................................................................................. 17
2.3.2
Minimization of a functional calculus .......................................................................................... 18
2.3.3
Method of minimization ..................................................................................................... 18
2.3.4
Application to the minimization of energy ............................................................................ 20
2.3.5
Determination of the pitch of advance .................................................................................... 21
2.4
Algorithm of
STAT_NON_LINE
................................................................................................... 23
3
Piloting ................................................................................................................................................ 24
4
Great deformations ......................................................................................................................... 24
4.1
Objective ........................................................................................................................................... 24
4.2
Great plastic deformations .................................................................................................. 25
4.3
Functionality
PETIT_REAC
.......................................................................................................... 26
5
Bibliography ........................................................................................................................................ 28
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
3/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
1 Presentation
1.1
General information on the relations of behavior
STAT_NON_LINE
is the operator of Code_Aster allowing to carry out mechanical calculations not
linear when the effects of inertia are neglected.
Calculation relates a priori only to the variables mechanical (displacements, stresses, variables
interns) by excluding any coupling with other physical phenomena (thermics,…). By
consequent, associated fields influencing the mechanical behavior (thermal fields,
metallurgical) are calculated as a preliminary by other operators (
THER_LINEAIRE
[U4.33.01],
THER_NON_LINE
[U4.33.02]), even by other codes (SYRTHES, SATURN,…).
This assumption suffers an exception with regard to modeling
thermo hydro-mechanical (modeling known as '
THM
') for which
STAT_NON_LINE
draft
the whole of the coupled problem of the equations of dissemination of thermics, the pressure of
fluid and of mechanical balance [R7.01.10].
In
STAT_NON_LINE
, two families of behaviors are available:
·
that which corresponds to the key word factor
COMP_ELAS
(
COMP
ortement
ELAS
tick) led to
through equilibrium equation to a nonlinear system depending explicitly on the field
displacements
U
compared to the configuration of reference, and parameterized by the moment of
calculation (through inter alia the thermal evolution). For more details, one will be able
to defer, for example, with the document [R5.03.20] concerning elasticity into large
transformations (hyper-elasticity), or the document [R5.03.21] on the thermoelastoplasticity with
isotropic work hardening,
·
the other family, which corresponds to the key word factor
COMP_INCR
(
COMP
ortement
INCR
émental),
is associated relations of behavior expressed by a differential equation
implicit (for example elastoplasticity, viscoplasticity, hypo-elasticity). In this case,
relation of behavior is integrated as presented for example in [R5.03.02]: in
connecting an increment of displacement
U
calculated starting from a mechanical state given (the state
mechanics being represented by a field of displacements
U
, a stress field
and a field of internal variables
) with the stress field at the moment
T
calculation.
The equilibrium equation thus leads to a nonlinear system in
U
, but which is also
parameterized by the moment of calculation through the facts of the case (variation of the loading
mechanics and thermal evolution for example).
In both cases, one calculates the solution gradually. It is theoretically not
essential in the nonlinear elastic case, but it may be that nonthe linearity of the solution
sought either too strong for the algorithm of resolution used, and that it is essential, for
numerical reasons, to operate step by step.
It is necessary nevertheless to have for the spirit the fundamental difference between the two approaches. The elastic case
suppose the existence of a state of reference per report/ratio to which the elastic strain is written: this
state corresponds in a state without forced deformation nor. It is the “absolute” value of the loading which
create the deformation. The incremental case is based on the state previously calculated and “forgets” all
reference to the former states except that given by the internal variables. In this case, it is
variation of the loading which modifies the state of the system: in particular, one needs a variation of the field of
temperature to create thermal deformations.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
4/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
1.2
Position of the nonlinear quasi-static problem
Consequently [§1.1], it is seen that it is legitimate to consider that the nonlinear problem that one
draft has like unknown factor a displacement and is parameterized by time. That is to say thus the problem not
linear quasi-static according to:
v R U
v L
v
Bv
Drunk
U
T
T
D
T
T
T
(,)
()
()
=
=
=



such as
0
where:
·
T
represent the variable of moment,
·
U
the field of displacement taken starting from a configuration of reference,
·
L
the mechanical loading to which the structure is subjected (pressure, imposed force,…),
·
the relation
Drunk U
=
D
T
()
corresponds to the boundary conditions imposed in displacements
(imposed displacements, connections between degrees of freedom,…) :
B
is a linear operator of
the space of the fields of displacements on a space of functions defined on part of
edge of the structure,
U
D
is a function given on this part,
·
R
represent the internal forces of the problem of quasi-static mechanics nonlinear
(in the linear case, one has
R U
Ku
(,)
T
=
, where
K
is the matrix of rigidity of the structure).
The notation
R U
(,)
T
is a short cut which one will use in the continuation.
In fact, more precisely,
R U
(,)
T
is connected to the stress field
by the operator of
work of the virtual deformations
Q
T
[§ 2.2.1] according to the relation:
R U
Q U
(,)
T
T
=
()
:
éq
1.2-1
In small displacements,
Q
T
is independent of displacements; for the large ones
displacements, to refer to [R5.03.20]. Subsequently, one will place oneself on the assumption
small displacements and small deformations.
The stress field
I
at the moment
T
I
is written
(,
)
U
H
I
I I
I
T T
- 1
, if one notes
T
I
the field of
temperatures and
H
I
- 1
last history of the structure. For the elastic behaviors, the history
does not intervene: the unit
H
I
- 1
is thus empty. For the incrémentaux behaviors, the history is
the whole of the states (fields of displacements, stresses and variables intern) at the moment
precedent:
{
}
H
U
I
I
I
I
I
T
-
-
-
-
-
=
1
1
1
1
1
,
,
,
.
In the general case, the dependence of the operator
R
is, as we saw in [§1.1],
implicit compared to time: it results from the integration of the relation of behavior in
time (for the problems of elastoplasticity for example). The dependence clarifies compared to
time is very rare: it in particular in the case of appears relations of fascinating behavior in
count a phenomenon of work hardening by time known as time-hardening.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
5/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
The dualisation of the boundary conditions of Dirichlet
Drunk U
=
D
T
()
conduit with the problem
according to [R3.03.01]:
R U
B
L
Drunk
U
(,)
()
()
T
T
T
T
D
+
=
=



éq
1.2-2
The unknown factors are now at any moment
T
the couple
()
U,
, where
represent the “multipliers
of Lagrange " of the boundary conditions of Dirichlet [R3.03.01]. The vector
B
T
be interpreted like
opposite of the reactions of support to the corresponding nodes.
The formulation of the quasi-static problem consists in expressing the balance of the structure (forces
external interns = forces) for a succession of moments of calculation
{}
T
I
I I
1
who parameterize it
loading:
R U
B
L
Drunk
U
(,)
()
()
I I
T
I
I
D
I
T
T
T
+
=
=



I
éq
1.2-3
what amounts cancelling in
(,)
U
I
I I
T
the vector
F U
(,)
T
defined by:
F U
R U
B
L
Drunk U
(,)
(,)
()
()
T
T
T
T
T
D
=
+
-
-




The state of the structure in
T
0
is supposed to be known. One carries out
I
increments (or not) of load definite
as follows:
T
0
T
1
T
2
….
T
I-1
T
I
1
2
I
not
charge n°
moment
The unknown factors are calculated in an incremental way by the total algorithm of resolution (even for
elastic behaviors). From
(
)
U
I
-
-
1
1
,
I
, solution satisfying balance in
T
I
- 1
, one
determine
U
I
and
I
, which will make it possible to obtain the solution in
T
I
:
U
U
U
I
I
I
I
I
I
=
+
=
+
-
-


1
1
Increments
U
I
and
I
are initially estimated by linearizing the problem compared to time
around
(
)
U
I
I
T
-
-
-
1
1
1
,
,
I
(phase known as of prediction or Euler [bib 3]). Then a method is used
of NEWTON or one of its alternatives to solve the equation [éq 1.2-3] in an iterative way (one
calculate a continuation
(
U
in
in
,
)
).
For the incremental relations of behavior, one needs to know in
T
I
- 1
stress field
I
- 1
and the field of internal variables
I
- 1
(Cf [R5.03.02] for one
example).
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
6/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
2
Method of resolution
2.1
Principle of the method of NEWTON
The method of NEWTON is a conventional method of resolution of the equations of the type:
F X
()
= 0
where
F
is a vector function (nonlinear) vector
X
.
It consists in building a vector series
{}
X
N
N
converging towards the solution
X
. To find
the new one reiterated
X
N
+1
, one approaches
F X
(
)
N
+1
by a development limited to command 1 around
X
N
, and
one expresses that
F X
(
)
N
+1
must be null:
0
1
1
=
+
-
+
+
F X
F X
F X
X
X
(
)
()
() (
)
N
N
N
N
N
'
,
that is to say:
[
]
F X
X
X
F X
X
X
F X
F X
'
'
() (
)
()
()
()
N
N
N
N
N
N
N
N
+
+
-
-
= -
=
-
1
1
1
or
Recall:
F 'X
()
N
is the tangent linear application associated the function
F
. The derivative at the point
X
in the direction
H
is defined like:
F X H
F X
H
F X
'()
lim
(
)
()
=
+
-
0
.
The matrix of
()
F 'X
in the bases chosen for the vector spaces concerned
be called the matrix jacobienne
F
at the point
X
.
When
F
is related to an Euclidean vector space with actual values,
()
F 'X
is one
linear form, and one can show that there is a vector (single), noted
F X
()
and called it
gradient of
F
, such as:
F 'X H
H
F X
()
()
=
T
(produced scalar of
H
and of the gradient).
When one is close to the solution, the convergence of the method of NEWTON is quadratic (it
numbers of 0 after the comma in the double error with each iteration: 0.19 - 0.036 - 0.0013 -
0.0000017 for example). But this method (using the true tangent) has several disadvantages:
·
it requires the calculation of the tangent to each iteration, which is expensive (especially in
case where a direct solvor is used),
·
if the increment is large, the tangent (known as coherent or consistent) can lead to
divergences of the algorithm,
·
it can not be symmetrical, which obliges to use particular solveurs.
For this reason one can use other matrices in the place of the tangent matrix:
stamp elastic, a tangent matrix obtained before, the symmetrized tangent matrix,…
[§2.2.1].
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
7/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
2.2
Adaptation of the method of NEWTON to the problem arising
2.2.1 Resolution without boundary conditions dualized
If one initially forgets the boundary conditions of Dirichlet, one must solve one
system of the form:
R U
L
(,)
I I
iméca
T
=
,
where
L
iméca
indicate the mechanical loading at the moment
T
I
.
Using the notations of [§ 2.1], that reverts cancelling the vector function
F
defined by:
F U
R U
L
(,)
(,)
I I
I I
iméca
T
T
=
-
Nodal forces
R
can symbolically be noted
Q
T
, where
Q
T
is the matrix associated with
the operator divergence (expression of the agricultural work of virtual deformations), with
(,
)
U
H
I
I
I
I
T T
- 1
.
One a:
(
)
(): (
)
()
Q
U
W
L
F W
W
T
K
I K
I
K
I
K
K
D
D
D
=
=
+
T
where
·
W
K
indicate the basic function associated with
K
ième
degree of freedom of the structure,
·
F
I
indicate the voluminal forces applying to the moment
T
I
on
,
·
T
I
indicate the surface forces applying to the moment
T
I
on the border
of
.
The matrix
Q
depends on displacements
U
in great displacements [R5.03.20].
The application of the method of NEWTON results in solving a linear succession of problems of the type
(
N
is the number of the iteration of NEWTON,
I
that of the variable of moment):
K
U
L
R
in
I
N
iméca
in
+
=
-
1
,
where one notes
U
U
U
in
in
in
+
+
=
-
1
1
the increment of displacement between two iterations of NEWTON
successive. The matrix
K
in
is the matrix of tangent rigidity in
U
in
and the vector
R
in
represent them
forces intern with
N
ième
iteration of NEWTON of
I
ème
no the load. Quantity
L
R
iméca
I
N
-
represent the not balanced forces, which one can also call the “residue”.
The matrix
K
in
is the matrix of the tangent linear application of the function
F
, it is thus worth:
K
F
U
R
U
L
U
U
U
U
in
T
T
iméca
T
I
N I
I
N I
I
N I
=
=
-
(
,)
(
,)
(
,)
.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
8/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
In the absence of following forces [§2.2.3], the second term is null. The matrix
K
in
is thus the derivative
at the point
U
in
nodal forces (or interns) compared to displacements:
K
R
U
U
in
T
I
N I
=
(
,)
A small error in the evaluation of the residue can have serious consequences, because it is its calculation
exact which guarantees, if one converges, that it will be towards the sought solution. On the other hand, it is not
always necessary to use the true tangent matrix, whose calculation and factorization are expensive.
For example, an alternative of the method uses the elastic matrix
K
0
.
Method using the true tangent matrix
K
in
(known as also coherent or consistent matrix)
be called the method of NEWTON; methods using of other matrices (such as for example
stamp elastic
K
0
) are called modified methods of NEWTON. The choice enters a matrix
tangent (the last obtained or a preceding matrix) is carried out in Code_Aster by
the intermediate of the key word
STAMP:
'
TANGENT
“or”
RUBBER BAND
'of the key word factor
NEWTON
. Moreover,
it is possible to use a matrix of discharge, i.e. of a matrix with internal variables
constants (the evolution of nonthe linearities is thus not taken into account in this matrix), in
below of a certain pitch of time, for certain laws of behavior. One will refer to
documentation [U4.51.03] for the use of this functionality.
The method of NEWTON with consistent tangent matrix can be illustrated simply using
diagram of the figure [Figure 2.2.1-a].
R
0
1
R
2
R
U
R
LR
1
LR
2
U
I
1
U
I
0
U
I
2
U
I
L
I
1
U
I
2
U
I
Appear 2.2.1-a
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
9/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
2.2.2 Resolution with boundary conditions dualized
When one takes into account the conditions of imposed displacements, the system to be solved is written:
R U
B
L
Drunk
U
(,)
I I
T
iméca
I
I
D
T
+
=
=



I
The symbol will be used
to note the increments since preceding balance (in
T
I
- 1
) of
displacements, of the parameters of Lagrange, the loading and imposed displacements:


U
U
U
L
L
L
U
U
U
I
I
I
I
I
I
iméca
iméca
iméca
id
id
I
D
=
-
=
-
=
-
=
-


-
-
-
-
1
1
1
1
Knowing
(
)
U
I
- 1
1
,
I
, the couple
(
)
,
U
I
I
will be determined by the resolution of the system:
R U
U
B
L
L
B U
U
U
U
(
,)
)
(
)
I
I
I
T
I
iméca
iméca
I
I
id
id
T
-
-
-
-
-
+
+
+
=
+
+
=
+



1
1
1
1
1
(
I
éq
2.2.2-1
One will use a method of NEWTON to solve this system. However, the experiment shows
that the convergence of the method of NEWTON is strongly dependant on a judicious choice on
the initial estimate: “more the initial estimate is close to the solution, plus the algorithm converges quickly”.
To start the iterative process of the method, it is thus useful to determine “a good” increment
initial
(
)
,
U
I
I
0
0
. For that, one linearizes compared to time the continuous problem: it is what one
call the phase of prediction (or of initialization). One connects with the loop of the iterations of
NEWTON which makes it possible, with convergence, to obtain the values of
(
)
,
U
I
I
, and thus those of
(
)
U
I
,
I
.

2.2.2.1 Phase of prediction
One linearizes the system [éq 2.2.2-1] compared to time around
(
)
U
I
- 1
1
,
I
; by taking account of
the balance obtained at the moment
T
I
- 1
, one obtains the equations allowing to calculate values
predictive
(
)
,
U
I
I
0
0
:
R
U
U
B
L
R
B U
U
U
I
I
I
T
I
I
méca
T
I
I
I
D
T
T
-
-
1
1
0
0
+
=
-
=



R
T
T
I
- 1
indicate the total differential of
R
Q U
,
= (): (
())
T T T
compared to T. This notation
particular be to draw the attention to the fact that

T
T
T
T
T
=
+
as one sees it in
continuation.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
10/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Therefore, while using [éq 1.2-1] and dependence of
compared to time and compared to
temperature:
Q
U
Q
U
U
B
L
Q
Q
B U
U
U
U
T
T
I
T
I
iméca
T
T
I
T
T
I
I
id
I
I
I
I
T
T
T
T
:
:
:
:

+
-
-
-
-
1
1
1
1
0
0


+
=
-
-
=


If the stresses depend on other variables, such as for example of the proportions
Z
various metallurgy i.e. phases
(, (), ())
T T T Z T
, one sees appearing to the second member
the term corresponding:
Q
T
T
I
Z
Z
I
:
- 1
. It is supposed thereafter that the effects
metallurgical [R4.04.02] are integrated into the second member representative of the effects
thermics. Moreover, one currently does not hold account for the phase of prediction of
dependence clarifies stresses compared to time. Lastly, dependence of
stamp
Q
T
compared to displacements is neglected on the assumption of small
displacements: the term
Q
U
U
T
I
- 1
, known as term of geometrical rigidity, thus disappears from
the preceding equation.
With obvious notations, the system of equations obtained is also written:
K
U
B
L
L
B U
U
I
I
T
I
iméca
ither
I
id
-
+
=
+
=



1
0
0
0
éq
2.2.2.1-1
Let us note that one can replace the matrix formally
K
I
- 1
, derived from
R
compared to
U
in
U
I
- 1
by the elastic matrix
K
0
.
For the developers, let us specify that the calculation of the tangent matrix at the time of the phase of prediction
is carried out via the option of calculation
RIGI_MECA_TANG
[§ 2.2.2.3].
The effective increment of loading
L
I
appearing to the second member includes/understands two terms:
·
the mechanical increment of loading
L
iméca
, composed of the dead loads
(independent of the geometry, like gravity) and of the following loads (dependant
geometry, like the following pressure [R3.03.04]). Actually, there is cases (it
first increment of load, for example) where
méca
I 1
-
L
is unknown: this fact
L
iméca
is
rather calculated by the relation
1
1
-
-
-
-
=
I
T
I
T
méca
I
méca
I
B
Q
L
L
, which returns exactly to
even, in the case running, by taking account of balance with the increment (i-1). One
will notice whereas the expression utilized the multipliers of Lagrange to the increment
(i-1), which is sometimes unknown (with the first increment of load, for example). But linearity
boundary conditions in imposed displacements, which results in the fact that the matrix
B
is constant, allows in this case an abuse language: one poses
0
=
I
and
truth
I
truth
I
I
0
1
0
+
=
-
(the exponent
truth
corresponds to the variables appearing indeed
in the system of equations above) and one solves:
,
1
0
1
0
0
1



-
=
-
+
=
+
-
-
-
I
D
I
I
I
T
ther
I
méca
I
I
T
I
I
U
B
U
U
B
Q
L
L
B
U
K
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
11/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
·
the increment of loading says “thermal” (and metallurgical [R4.04.02])
L
ither
, resulting from
derivation of the nodal forces compared to the temperature and which is an estimate of the effect
of a variation in temperature. Indeed, if one notes
K
the module of compression
hydrostatic and
the thermal expansion factor, the increment of loading
“thermal” is written (since
R Q
=
T
and
ther
I
K
T
= - 3
I
D
where
I
D
is the matrix
identity):
-
L
Q
Q
I
ther
T
T
I
T
I
T
T
K
T
I
= -
=
:
:

1
3
I
D
,
In the elastic case, they are the nodal forces associated a thermal dilation (it
is not strictly speaking a loading, that is assimilated rather to the effect of a deformation
initial). This estimate is used in the phase of prediction and the criterion of stop. If
thermal dilations make leave the structure of the elastic range (plasticity by
example), this estimate will be corrected at the time of the iterations of NEWTON.

Note:
A particular case relates to the use of an excitation of the type
TYPE_CHARGE:
'
DIDI
'meaning
Differential Dirichlet, i.e. compared to the initial state. That consists, for the conditions with
limits of the blockings type, to impose, not
B U U
=
D
, but
B U U
U
(
)
-
=
0
D
. In this case, it
system to be solved in the phase of prediction for the new increment of load is:



=
-
+
=
+
-
-
D
I
I
I
T
ther
I
méca
I
I
T
I
I
U
U
B
Q
L
L
B
U
K
0
1
0
0
1


2.2.2.2 Loop on the iterations of NEWTON
One must find the values
(
)
,
U
I
I
increments of displacements and parameters of Lagrange
since the values
(
)
U
I
- 1
1
,
I
obtained with preceding balance (urgent
T
I
- 1
). One takes as
initial values
(
)
,
U
I
I
0
0
obtained at the end of the phase of prediction, before beginning them
iterations of the method of NEWTON.
With each iteration, one must solve a system allowing to determine
(
U
I
N
I
N
+
+
1
1
,
)
, increments
displacements and parameters of Lagrange since the result
(
)
U
in
in
,
iteration
the preceding one:
K
U
B
L
R U
B
B U
in
in
T
in
iméca
I
N
T in
in
+
+
+
+
=
-
-
=



1
1
1
0
(
)
éq
2.2.2.2-1
with
R U
Q
(
)
in
T in
=
, stresses
I
N
being calculated starting from displacements
U
in
by
the intermediate of the relation of behavior of the material [§1.1]. In fact, in the case of them
incrémentaux behaviors,
I
N
is calculated from (
I
- 1
,
I
- 1
) and of the increment of
deformation
(U
I
N
+1
)
induced by the increment of displacement
U
U
U
U
in
I
N
in
I
+
+
-
=
+
-
1
1
1
.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
12/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Stamp “tangent”
As in the phase of prediction, one is not obliged to use the true tangent matrix
K
in
. In
private individual, the operator
STAT_NON_LINE
authorize the use of the elastic matrix
K
0
, or
reactualization of the tangent matrix all them
I
0
no time (key word
REAC_INCR
) or all them
N
0
iterations of NEWTON (key word
REAC_ITER
). Thus, the matrix
K
in
perhaps replaced by one
stamp
K
J
, with
J I
, or a matrix
K
im
, with
m N
.
Caution: a “stiff” matrix too does not pose problems of stability but can produce
a very slow convergence; a “flexible” matrix too can lead to divergence, it is
advised in this case to make linear search [§ 2.3].
It is difficult to give a rule making it possible to know when one must reactualize the tangent matrix:
that strongly depends on the degree of nonlinearity of the problem and the size of the increments of load.
In discharge, it is recommended either to use the elastic matrix for the phase of prediction
and of resolution, the elastic matrix for the phase of prediction then the matrix is to use
tangent for the resolution.
The figures hereafter illustrate the various possibilities of reactualization of the tangent matrix:
stamp elastic
K
0
used with each iteration [Figure 2.2.2.2-a], tangent matrix reactualized with
each iteration and with each pitch of time [Figure 2.2.2.2-b], reactualized tangent matrix all them
I
0
no time (
I
0
=1 here) [Figure 2.2.2.2-c], and stamps tangent reactualized all them
N
0
iterations
of NEWTON (
N
0
=2 here) [Figure 2.2.2.2-d].

R
U
L
I
Appear 2.2.2.2-a: use of the elastic matrix
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
13/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
R
U
L
I
Appear 2.2.2.2-b: use of the true revalued tangent matrix
with each iteration
R
U
L
I
L
i+1
Appear 2.2.2.2-c: use of the tangent matrix
revalued with each pitch of time
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
14/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
R
U
L
I
Appear 2.2.2.2-d: use of the revalued tangent matrix
all 2 iterations of NEWTON
Method of NEWTON modified (using another matrix that the consistent tangent matrix)
converge less quickly than the method of NEWTON conventional, but is less expensive. It is
all the more economic as the number of degrees of freedom of the system is high. This is why one
can advise, when the system is of big size, to keep the same tangent matrix during
some iterations. Lastly, let us not forget to announce that in certain cases, it is calculation with
stamp elastic which is fastest in terms of calculating time, even if the iteration count
carried out is much more important (they are cheap iterations since the matrix is not
calculated and factorized that only once); moreover, it is the elastic matrix which it is recommended
to use during the discharges.
As the equation [éq 2.2.2.2-1] shows it, it is necessary to carry out with each iteration of NEWTON calculation
possible of the new tangent matrix
K
in
and of the “nodal forces”
R
B
in
T in
+
: for
developers, let us specify that these operations are carried out by the option of calculation
FULL_MECA
(
RAPH_MECA
if the tangent matrix is not recomputed).
The tangent matrix obtained by the option
RIGI_MECA_TANG
and the tangent matrix obtained
by the option
FULL_MECA
are in general calculated differently [§ 2.2.2.3].
Updating of the unknown factors
The updating of displacements and their increments is done as follows, that of the parameters of
Lagrange is carried out in an identical way:
U
U
U
U
U
U
in
in
in
in
I
N
I
N
+
+
+
+



1
1
1
1
=
+
=
+
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
15/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Criteria of convergence
There are three criteria of convergence.
The criterion
RESI_GLOB_MAXI
consist in checking that the infinite standard of the residue
-
+
méca
I
N
I
T
N
I
T
L
B
Q
, is lower than the value specified by the user. He is not advised
to use this criterion alone, because one cannot easily have an idea of the absolute orders of magnitude
acceptable.

The criterion of convergence chosen by defect amounts checking that the residue is sufficiently small,
as previously, and this relative with a quantity representative of the loading (it is it
criterion
RESI_GLOB_RELA
):
Q
B
L
L
L
B
T
I
N
T
I
N
I
méca
I
méca
I
ther
T
I
N
-
+
+
-
,
éq
2.2.2.2-2
being desired relative precision given by the user (or the default value of 10
- 6
) and
the standard of the maximum.
One can notice that, in the case of use of
RESI_GLOB_RELA
, the criterion can become
singular if the external loading
L
L
B
I
méca
I
ther
T
I
N
+
-
becomes null. This can arrive in
case of total discharge of the structure. If such a case of figure arises (i.e loading
10
- 6
time smaller than the smallest loading observed until now the increment), the code
use the criterion then
RESI_GLOB_MAXI
with like value that observed with convergence
preceding increment. When the loading becomes again not no one, one returns to the criterion
initial.
The third criterion is criterion RESI_REFE_RELA: the idea of this criterion is to build a force
nodal of reference, which will be used to estimate term in the long term, the nullity (approximate) of the residue:
{}
(
)
ref.
J
méca
I
N
I
T
N
I
T
J
F
ddls
J
-
+
L
B
Q
More precisely, the nodal force of reference is built starting from the data of an amplitude of
stress of reference (in mechanics; in the case of the THM, it is necessary to give one
reference to each physical phenomenon entering the generalized stress)
ref.
. From
this amplitude of stress of reference, one defines the tensor
test
J
: it is null for all these
components, except j-ième which is worth
ref.
. One defines then, for each node of each element
following nodal force (the goal being to give an idea of the importance of a component in a point
of Gauss of the stress on the nodal force):
=
=
=
NR
M
J
test
J
J
I
E
I
B
NR
R
1
1
,
1
~
with NR the number of points of Gauss of the element, M the component count of the tensor of
stress; the exponent
being used to note the evalutation of quantity at the point of Gauss
,
are them
weight of the points of Gauss.
The nodal force of reference is then defined by:
E
I
E
ref.
I
R
F
I
~
min
=
where
I
is the whole of the elements connected to the node
I
.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
16/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Convergence is issued carried out when all the criteria specified by the user are checked
at the same time. By defect, one makes a test on the relative total residue (
RESI_GLOB_RELA
) and numbers it
maximum of iterations of NEWTON (
ITER_GLOB_MAXI
).

2.2.2.3
RIGI_MECA_TANG
and
FULL_MECA
It is important to stress that the tangent matrix resulting from the option
RIGI_MECA_TANG
and stamps it
tangent resulting from the option
FULL_MECA
are in general not identical.
Let us suppose that one reached convergence for the moment
T
I
- 1
and that one seeks now with
to obtain balance for the next moment
T
I
. The matrix resulting from IH GI _MECA_TANG comes from one
linearization of the equilibrium equations compared to time around
(
)
U
I
- 1
1
,
I
i.e around
balance at the moment
T
I
- 1
. It is thus the tangent matrix of the system with convergence at the moment
T
I
- 1
.
On the other hand, the matrix resulting from
FULL_MECA
comes from a linearization of the equilibrium equations by
report/ratio with displacement around
(
)
U
I
N
I
N
,
i.e around balance at the moment
T
I
.
Moreover, one can interpret the differences enters
RIGI_MECA_TANG
and
FULL_MECA
in others
terms. One can thus show that the matrix resulting from
RIGI_MECA_TANG
corresponds to the operator
tangent of the problem continuous in time, known as also problem of speed (and connects the speed of stress
at the speed of deformation), whereas the matrix resulting from
FULL_MECA
corresponds to the operator
tangent of the problem discretized in time [bib1]. The document [R5.03.02] gives the expression in
each of the two cases for the relation of elastoplasticity of Von Mises with isotropic work hardening.
It is reminded the meeting that the processing of a relation of behavior [R5.03.02 § 5] consists with:
·
to calculate the stresses
I
N
and internal variables
in
starting from the initial state
(
,
)
I
I
-
-
1
1
and of the increment of displacement
U
I
N
,
·
to calculate the nodal forces
R
Q
in
T in
=
and reactions of support
B
T in
,
·
to calculate (possibly) the tangent matrix (option
RIGI_MECA_TANG
for the phase of
prediction, option
FULL_MECA
for the iterations of NEWTON
)
.

2.2.3 Case of the following loadings
A following loading (in mechanics) is a loading which depends on the geometry of the structure,
as for example the pressure which is exerted in the direction opposed to the normal (or the forces
of inertia in a reference mark not galiléen). Thus, when the structure becomes deformed with the evolution of
charge, the loading, expressed in an absolute reference mark, is transformed. The loads which do not depend
no the geometry of the structure are called dead or fixed loads (for example,
gravity). To indicate that a load must be treated like a following load in
STAT_NON_LINE
, one indicates
TYPE_CHARGE:
“SUIV”
under the key word
EXCIT
.
A mechanical loading
L
méca
T
()
comprising following loads is written
L
L
U
fixed
suiv
T
T
()
(,)
+
(the fixed exponent indicates here the died loads, and suiv the following loads). The system of equations
to solve becomes then:
R U
B
L
L
U
Drunk
U
(,)
()
(,)
I I
T
fix I
suiv
I I
I
id
T
T
T
+
=
+
=



I
Operations of derivation allowing to write the phase of prediction and the iterations of the method
of NEWTON thus utilize the derivative of
L
suiv
compared to displacements
U
.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
17/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
The phase of prediction becomes:
(K
L
U
U
B
L
L
L
B U
U
U
I
suiv
I
T
I
ifixe
isuiv
ither
I
id
I
-
-
+
=
+
+
=




1
0
0
0
1
-
)
and the iterations of NEWTON consist in solving the system:
(
(
)
(
)
K
L
U
U
B
L
L
U
R U
B
B U
U
I
N
suiv
I
N
T
in
ifixe
suiv
I
N
in
T in
in
I
N
-
+
=
+
-
-
=




+
+
+
)
1
1
1
0
Thus, with the beginning of each pitch of load (prediction) and with each iteration of NEWTON, one must
to calculate a matrix of “rigidity”
(
)
-
L
U
U
suiv
I
N
/
and a vector
L
U
suiv
I
N
(
)
dependant on the loadings
follower.
The only loads which can be treated like following loads in the current state of
the operator
STAT_NON_LINE
are:
·
pressure for modelings
3D
,
3d_SI
,
D_PLAN
,
D_PLAN_SI
,
AXIS
,
AXIS_SI
,
C_PLAN
,
C_PLAN_SI
[R3.03.04],
·
the loading of gravity for the elements
CABLE_POULIE
[R3.08.05], elements with three
nodes comprising a pulley and two strands of cables: the force of gravity being exerted on
the element depends on the respective lengths of the two strands,
·
the centrifugal force in great displacements, which for a rotational speed
is given
by:
[
]
[
]
+
=
OM
OM
U
D
D
(
)
0
;
·
the loading of gravity for all modelings
THM
porous environments not
saturated [R7.01.10]: indeed, the density depends on the nodal variables
U,
p
T
and
to take account of the relations of behavior of the géomatériaux one.
2.3 Seek
linear
Linear search here exposed relates to linear search in the absence of piloting. For
description of linear search in the presence of piloting, one will refer to documentation
[R5.03.80].
2.3.1 Principle
The introduction of linear search into the operator
STAT_NON_LINE
result from an official report:
method of NEWTON with consistent matrix does not converge in all the cases of figure,
in particular when one leaves too much far from the solution. In addition, the use of matrices other than
stamp tangent consistent can, when they are too “flexible”, lead to divergence.
linear search makes it possible to be guarded against such divergences.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
18/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
It consists in considering
(
U
in
in
+
+
1
1
,
)
, either like the increment of displacements and of
parameters of Lagrange, but as a direction of search in which one will seek with
to minimize a functional calculus (the energy of the structure). A pitch of advance will be found
in this
direction, and the updating of the unknown factors will consist in making:
=
+
= +
U
U
U
in
I
N
in
in
in
in
+
+
+
+



1
1
1
1
In the absence of linear search (by defect) the scalar
is of course equal to 1.
2.3.2 Minimization of a functional calculus
In order to be better convinced of the founded good of linear search, one can interpret the method of
NEWTON like a method of minimization of a functional calculus (if matrices
tangents are symmetrical). We insist on the fact that the equations obtained are
rigorously those of the method of NEWTON exposed in [§2.2] and that only the way of y
to arrive is different.
“The talk the dualisation of the boundary conditions of Dirichlet and we place Forget” to simplify
on the assumption of the small deformations. The functional calculus is considered:
J: V
J
=
-
-
IR
()
(())
.
.
U
U
U
F U
U
W
D
D
D
T
where density of free energy
W
allows to connect the tensor of the stresses
with the tensor of
linearized deformations
by the relation
=
W
.
The functional calculus
J
being convex, to find the minimum of
J
is equivalent to cancel its gradient, that is to say:
=
J
V
().
U v
v
0
,
what is exactly the Principle of Virtual Work since:
=
-
-
J ().
(): ()
.
.
U v
U
F
v
v
T v
D
D
D
Thus, to solve the equations resulting from the Principle of Virtual Work (bases problem formulated in
[§1.2]) is equivalent to minimize the functional calculus
J
who represents the energy of the structure (energy
decreased intern of the work of the external forces
F
and
T
).
2.3.3 Method of minimization
Minimization is made in an iterative way, classically in two times with each iteration:
·
calculation of a direction of search
along which one will seek reiterated according to,
·
calculation of the “best” pitch of advance
in this direction:
U
U
N
N
+
=
+
1
.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
19/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
In a problem of minimization, the natural idea is to advance in the direction opposed to the gradient
functional calculus, which is locally the best direction of descent since this direction is
normal with the lines of isovaleurs and directed in the direction of the decreasing values [Figure 2.3.3-a].
U
N
-
J
(U
N
)
J
(U) =
cste
Appear 2.3.3-a
However, it is possible to improve the choice of the direction of descent by using this method
of gradient in a metric news. It is what will enable us to find the equations
conventional of the method of NEWTON [§2.3.4].
Let us take the simple example of a problem with two variables
X
and
y
for which the functional calculus with the form
of an ellipse whose minimum is in (
/
has
,
/
b):
J (,)
X y
ax
by
X
y
=
+
-
-
1
2
1
2
2
2
By choosing like direction of descent the reverse of the gradient of
J
, one passes from one reiterated to the following
(let us reason on
X
only) by:
X
X
ax
N
N
N
+
=
-
-
1
(
)
,
who does not point towards the solution since the normal in a point of an ellipse does not pass in general
not by the center of the ellipse [Figure 2.3.3-b].
has
B
y
X
Appear 2.3.3-b
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
20/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
On the other hand, if one carries out a change of variables so that the isovaleurs of
J
become
circles [Figure 2.3.3-c]:
X
X has
y
B y
X y
X
y
X has
B y
=
=
=
+
-
-




J (,)
(
)
1
2
2
2
B
has
X
y
Appear 2.3.3-c
The use of the direction reverses gradient of
J
then allows to obtain
the solution in one
iteration:
X
X
X
has
has
X
has
N
N
N
N
+
+
=
-
-
=
=
1
1
(
)
Thus, the use of the method of gradient in the metric news allows a convergence
immediate. In a more intricate case (functional calculus convex but different from an ellipse),
convergence is not instantaneous but the change of variables makes it possible to reduce appreciably
the iteration count necessary.
2.3.4 Application to the minimization of energy
To simplify, one will place oneself in the discretized linear case where the functional calculus
J
is worth:
J ()
U
the U.K.U
U L
=
-
1
2
T
T
O one notes
K
the matrix of rigidity of the structure, and
L
the vector of the imposed loadings.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
21/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
To minimize
J
, we will use the same method of descent as previously while making with
precondition a completely similar change of variables. The matrix
K
being symmetrical defined
positive, its eigenvalues are real positive: one can thus define the “square root” of
K
that
one will note
K
(also symmetrical). One poses
U
K U
=
, the minimization of
J
is then
equivalent to that of:
J ()
U
U U
U
K L
=
-
-
1
2
1
T
T
While taking as direction of descent the direction reverses gradient of
J
, one obtains:
U
U
U
K L
N
N
N
+
-
=
-
-
1
1
(
)
Maybe, while returning to the initial variables:
U
U
K
K U
L
N
N
N
+
-
=
-
-
1
1
(
)
Or:
K U
U
L
K U
(
)
N
N
N
+
-
=
-
1
One finds the equations of the method of NEWTON: the matrix
K
is Hessienne of
functional calculus
J
(matrix of the derived second) and the second member is the difference of the loading
and of the nodal forces, also called residue.
Thus method of NEWTON perhaps interpreted like resulting from the minimization of energy
structure via a method of gradient applied after a change of metric.
2.3.5 Determination of the pitch of advance
Let us return to the real problem, that of the resolution of
R U
L
(,)
I I
I
T
=
. This problem can be interpreted
like the minimization of:
W
I
iT I
()
U
U L
-
,
where
W
I
()
U
corresponds to the discretization, on the basis of function of form, of energy interns
structure
W
D
(())
U
.
One calculated by the method of NEWTON an increment of displacement
U
in+1
who, in the problem
of minimization, is interpreted like a direction of search. One will calculate the pitch of advance
in this direction allowing to minimize the value of the functional calculus:
{
}
Min
+
-
+
+
+
R
W
in
in
I
I
N
I
N
(
)
(
)
U
U
L U
U
1
1
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
22/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
To find the minimum of this scalar function of
that one will note
()
F
, one seeks the point where
its derivative is cancelled (that amounts making orthogonal the final residue and the direction of search):
[
] [
]
F
I
N
T
T
in
in
I
'
()
(
)
=
+
-
=
+
+
U
Q
U
U
L
1
1
0
(
F
'
()
is the projection of the residue on the direction of search)
With the notations of [§2.2.2] and by taking of account reactions of support, the scalar equation with
to solve to determine the pitch of advance
, is written:
[
] [
]
U
Q
U
U
B
L
in
T
T
I
N
I
N
T
in
iméca
+
+
+
+
+
+
-
=
1
1
1
0
(
)
(
)
I
N
So that determination of
is not too expensive, one uses a method of secant of which it
numbers maximum iterations is fixed by the user. The method of secant can be interpreted like
a method of NEWTON where the derivative at the point running is approached by the direction joining it
not running and the preceding point:
p
p
p
p
p
p
p
p
p
p p
p
p
G
G
G
G
G
G
G
+
-
-
-
-
-
=
-
-
-
=
-
-
1
1
1
1
1
1
,
where one noted
G
F
p
p
=
'
(
)
.
One leaves
0
0
=
and
1
1
=
. The method of secant has a command of convergence of about 1.6
[bib2]. It is represented schematically on the figure [Figure 2.3.5-a].
0
1
0
1
2
3
4
G
= F '(
)
Appear 2.3.5-a
At the end of linear search, one brings up to date displacements and parameters of Lagrange by:
=
+
= +
U
U
U
in
I
N
in
in
in
in
+
+
+
+



1
1
1
1
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
23/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
The test of convergence carries:
·
on the maximum number of iterations of linear search indicated by the user under
key word
ITER_LINE_MAXI
key word factor
NEWTON
(the default value 0 inhibits
linear search, and
1 is worth then),
·
on the criterion
RESI_LINE_RELA
given by
F
F
()
. ()
0
, where
is worth by defect 0.1.
Linear search is to some extent a “policy insurance” allowing of
to secure against serious divergences of the method of NEWTON. When direction of
seek
(
U
in
in
+
+
1
1
,
)
is “bad” (if the tangent matrix is too flexible, by
example), the linear algorithm of search leads to a low value of
, which avoids
of going “in the decoration”. It is not necessary to do many iterations in the method
of secant (2 or 3 are enough to avoid the catastrophes) because each one is rather expensive (it
is necessary to recompute the internal forces) and there is not the ambition to find with each iteration of
NEWTON it
really optimal.
2.4 Algorithm
of
STAT_NON_LINE
One will use like previously the index
I
(like “moment”) to note the number of an increment of
charge and the exponent
N
(like “newton”) to note the number of the iteration of NEWTON in progress.
The algorithm used in the operator
STAT_NON_LINE
can then be written way schematically
following:
(
)
U
0
0
,
and
0
known
Loop over moments
T
I
(or increments of load): loading
()
L
L
I
I
T
=
·
(
)
U
I
-
-
1
1
,
I
known
·
Prediction: calculation of
U
I
0
and
i0
·
Loop on iterations of NEWTON: calculation of a continuation
(
)
U
in
I
N
,
-
(
)
U
in
I
N
,
and
(
)
U
in
I
N
,
known
-
Calculation of the matrices and vectors associated with the following loads
-
Expression of the relation of behavior
-
calculation of the stresses
I
N
and of the internal variables
in
starting from the values
I
- 1
and
I
- 1
with preceding balance (
T
I
- 1
) and of the increment of displacement
U
U
U
I
N
in
I
=
-
- 1
since this balance
-
calculation of the “nodal forces”:
Q
B
in
T in
+
-
possible calculation of the matrix of tangent stiffness:
K
K U
in
in
= ()
-
Calculation of the direction of search
(
U
in
in
+
+
1
1
,
)
by resolution of a system
linear
-
Iterations of linear search:
-
Updating of the variables and their increments:
U
U
U
U
U
U
in
I
N
in
I
N
I
N
in
I
N
in
in
I
N
in
I
N
+
+
+
+
+
+
+
+






1
1
1
1
1
1
1
1
=
+
= +
=
+
= +
and
-
Test of convergence
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
24/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
·
Filing of the results at the moment
T
I
U
U
U
I
I
I
I
I
I
I
I
=
+
=
+

-
-


1
1
It is noticed that there are three overlapping levels of iterations: a loop external on the pitches of time,
a loop of iterations (known as total) of NEWTON and subloops possible for
linear search (if it is asked by the user) and certain relations of behavior
requiring iterations (known as interns), for example for elastoplasticity in plane stresses.

3 Piloting
One will refer to documentation [R5.03.80].

4 Large
deformations
4.1 Objective
Until now, we made the assumption that displacements and deformations remained moderate
so as to respect the assumption of the small disturbances. This assumption becomes null and void in
many cases and the method of resolution of the problem must then integrate the evolution of
geometry of the problem, to handle a particular kinematics and to use an adequate formulation of
the law of behavior.
In practice, the assumption of the small deformations can be applied as long as the square of
modulus of deformation remains lower than the precision of calculations considered. In the same way,
the assumption of small rotations can be applied as long as the product enters the square of
the swing angle and the modulus of deformation remain lower than the precision of calculations
considered.
Various alternatives exist within Code_Aster; our objective is not here to do one of them
detailed presentation and we return to the various documents treating specifically each
problems:
·
Relation of nonlinear elastic behavior in great displacements: [R5.03.20],
·
Beams in great displacements: [R5.03.40],
·
Voluminal elements of hulls into nonlinear geometrical: [R5.07.05],
·
Elasto (visco) plasticity, metallurgy and great deformations: [R4.04.03],
·
Elastoplastic modeling with isotropic work hardening in great deformations: [R5.03.21].
The objective is here to present a general functionality of the code which allows a simple processing of
problems great deformations: the argument
PETIT_REAC
key word
DEFORMATION
under
key word factor
COMP_INCR
.
In the continuation, we will limit ourselves to the case of the great plastic deformations, which allows of good
to include/understand the problems
PETIT_REAC
. We base ourselves for that on [bib4] and on
documentation [R5.03.21].
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
25/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
4.2
Great plastic deformations
Initially let us point out the requirements of a modeling great plastic deformations
in mechanics of the continuous mediums in term of equilibrium equations, of kinematic description and of
relation of behavior so as to determine the limitations of modeling well
PETIT_REAC
.
Equilibrium equations
If one makes the choice write the equilibrium equations on the current configuration
and to use it
tensor of stresses of Cauchy
, they are summarized with:
:
.
.
=
+
v
F v
T v
v V
D
D
D
AD
éq
4.2-1
The preceding equation can be written on the configuration of reference; it is besides what is made
for nonlinear elasticity in great deformations [R5.03.20]. In the case of laws of behavior
of incremental nature i.e whose evolution is controlled by the current state, this writing loses of sound
interest. One prefers then the writing [éq 4.2-1], although current configuration
that is to say an unknown factor of
problem.
Kinematic description
For the laws of plastic behavior, one uses commonly like measures deformation it
rate of deformation
D
:
D
U
U
= 12 (&
&)
+
T
Moreover, one defines also the rate of rotation
W
:
W
U
U
= 12 (&
&)
-
T
Let us stress that in the preceding expressions, the operator gradient is defined on the configuration
current.
Elastoplastic relation of behavior
Before writing the relations of behavior, we will make three simplifying assumptions:
·
the elastic strain are small in front of the plastic deformations (what is checked
generally well in the case of metals). This makes it possible to break up the rate of
deformation in an additive way in a plastic part and an elastic part:
D D
D
=
+
E
p
·
the plastic deformation is isochoric. This makes it possible to write that the tensors of stresses of
Kirchhoff
(stress measurement adequate to the great deformations) and of Cauchy
(intervening in the equilibrium equations) are identical to the first command.
·
the behavior is isotropic.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
26/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Under these assumptions, the laws of behavior in great deformations can be written under one
form near to the small transformations. Thus for the plasticity of Von-Settings to isotropic work hardening
one obtains:
J
eq
eq
=
=




WITH D D
D
:(
(
(
&
~
E
p
p
-
)
F
p)
p)
, =
-
R
p
3
2
éq 4.2-2
where
·
~
is the diverter of
,
·
With
is the tensor of elasticity,
·
J
indicate it derived from Jaumann of
, which takes account of terms of transport defined in
the aid of the rate of rotation
W
and which is given by:
J
= &
-
.
W
W
.
+
4.3 Functionality
PETIT_REAC
The principle of modeling
PETIT_REAC
simply consist in reactualizing the geometry of
problem during iterations of Newton (and not at the end of each pitch of time). This means
that all the quantities intervening in the equations of the problem are evaluated on the configuration
current. Anything else is not modified compared to the case small disturbances.
We have just seen which products are necessary to a “clean” integration of a law of
elastoplastic behavior in great deformations, with the help of three simplifying assumptions.
Now let us detail the differences which the resolution implies by
PETIT_REAC
.
Equilibrium equations
Taking into account the geometrical reactualization, the interior efforts, first term of [éq 4.2-1], are
estimated in an exact way, with the help of the approximation by finite elements. The calculation of the efforts
outsides east him independent of the resolution by
PETIT_REAC
: the dead loads are calculated
on the configuration of reference and the loads following on the current configuration. Equations
with balance are thus dealt in an exact way.
Kinematic description
Kinematic description is the same one as that of the small disturbances. This means that one
increment of deformation is calculated by:
=
+
1
2 (
(
(
)
U
U
)
)
T
The total deflection is then the sum of each one of these increments of linearized deformation,
calculated on different configurations. It is thus delicate to give him a physical direction and better
is worth to use it like an indicator of the level of deformation reached.
Elastoplastic relation of behavior
The difference between the formulation [eq 4.2-2] and the formulation small disturbances lies in
replacement of derived from Jaumann by the simple derivative in time. The disadvantage of this
last is due to the fact that it is not incrémentalement objective i.e not invariant by rigid rotation
structure (for more details, to see documentation [R5.03.21]).
The employment of
PETIT_REAC
is thus not appropriate to great rotations but it is it with large
deformations.
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
27/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
Integration within Code_Aster
In term of finite elements, the resolution by
PETIT_REAC
imply with each pitch of load
resolution of the same nonlinear system as previously:
R U
B
L
Drunk
U
(,)
()
()
I I
T
I
I
D
I
T
T
T
+
=
=



I
With the difference close the forces intern are formally calculated by:
R U
U
(,)
():
I
I
T
I
T
= Q
where the operator
Q
depends on displacements.
Within this framework, the calculation of the tangent matrix carries out formally to:
K
R
U
U
U
U
U
U
U
U
I
N
T
T
T
I
N
I
I
N
I
I
N
I
=
=
+
(
,)
(
,)
(
,)
Q
Q
():
()
:
The first term is the contribution of the behavior, similar to what was presented into small
transformations, with the difference which this contribution is evaluated here in current configuration.
second term is the contribution of the geometry which is not present in small transformations.
Within the framework of the resolution
PETIT_REAC
, this term is not present in the calculation of the matrix
tangent. One thus has:
K
U
U
U
I
N
T
I
N
I
= Q ():
(
,)
One can announce that the absence of the geometrical contribution in the tangent matrix can sometimes
to make convergence difficult. Moreover, the presence of important plastic deformations can
to make the problem quasi-incompressible, because of isochoric nature of these deformations.
problems of numerical blocking can then appear which can be circumvented by employment
elements under integrated
background image
Code_Aster
®
Version
7.4
Titrate:
Quasi-static nonlinear algorithm
Date:
06/07/05
Author (S):
P. BADEL, J. LAVERNE, NR. TARDIEU
Key
:
R5.03.01-D
Page
:
28/28
Manual of Reference
R5.03 booklet: Nonlinear mechanics
HT-66/05/002/A
5 Bibliography
[1]
P. MIALON: Elements of analysis and numerical resolution of the relations of elastoplasticity.
EDF - Bulletin of the Management of the Studies and Search - Series C - N° 3 1986, p. 57 - 89.
[2]
J.F. MASTER: Numerical analysis, duplicated course of the ENTPE.
[3]
J. SHI & Mr. A. CRISFIELD: “Combining arc-length control and line searches in path
following ", Communications in Numerical Methods in Engineering, vol. 11, 793-803 (1995).
[4]
E. LORENTZ: Great plastic deformations, modeling in Aster by
PETIT_REAC
,
EDF - CR MN 1536/07.