Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
1/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
Organization (S):
EDF-R & D/AMA
Instruction manual
U4.5- booklet: Method of resolution
Document: U4.51.03
Operator
STAT_NON_LINE
1 Goal
To calculate the mechanical or thermo evolution hydro-mechanical coupled, into quasi-static, of a structure
into nonlinear.
Nonthe linearity is related either to the behavior of material (for example plastic), or with the geometry
(for example in great displacements). To have details on the method of resolution employed,
one will refer to the reference material [R5.03.01].
The evolution can be studied in several successive work (réentrant concept), that is to say in continuation (it
last calculated moment is the initial moment of following calculation), that is to say in recovery on the basis of one moment
former.
If time necessary to carry out calculation is not sufficient, the program stops, but them
already calculated results are backed up if a data base were defined in the profile of study of
the user. Product a structure of data of the type
evol_noli
.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
2/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
Count
matters
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
3/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
4/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
2 Syntax
statnl [evol_noli] = STAT_NON_LINE
(
reuse
= statnl, [evol_noli]
MODEL
= Mo,
[model]
CHAM_MATER
=
chmat,
[cham_mater]
CARA_ELEM
=
carac,
[cara_elem]
EXCIT =
_F (
CHARGE
=
chi,
[char_meca]
FONC_MULT
= fi,
[function/formula]
TYPE_CHARGE
=
/
“FIXE_CSTE”
[DEFECT]
/
“FIXE_PILO”
/
“SUIV”
/
“DIDI”
),
SOUS_STRUC = _F (
CAS_CHARGE
=
chi,
[char_meca]
/
ALL =
“YES”,
[DEFECT]
/
NET
=
lma, [l_maille]
),
| COMP_INCR = _F (see the document [U4.51.11]
),
| COMP_ELAS
=_F
(see the document [U4.51.11]
),
VARI_COMM
=_F
(
/
IRRA
=
will irra [evol_varc]
/
CORROSION
=
corro
[evol_varc]
),
ETAT_INIT
=_F
(
/
|
SIGM
=
sig, [cham_elem_SIEF_R]
[carte_SIEF_R]
|
VARI
=
vain,
[cham_elem_VARI_R]
|
DEPL
=
depl,
[cham_no_DEPL_R]
| VARI_NON_LOCAL
=
vanolo, [cham_no_VANL_R]
/
EVOL_NOLI
=
evol,
[evol_noli]
/
NUME_ORDRE=
nuini,
[I]
/
INST
=
instini,
[R]
PRECISION
=
/
1.0E-3,
[DEFECT]
/
prec, [R]
CRITERION =
/
“RELATIVE”,
[DEFECT]
/
“ABSOLUTE”,
NUME_DIDI
=
nudidi, [I]
INST_ETAT_INIT
=
istetaini
[R]
),
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
5/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
INCREMENT
=_F
(
LIST_INST
=
litps,
[listr8]
EVOLUTION
=
/
“CHRONOLOGICAL”,
[DEFECT]
/
“RETROGRESSES”,
/
“WITHOUT”,
/NUME_INST_INIT
=
nuini,
[I]
/
INST_INIT
=
instini,
[R]
/
NUME_INST_FIN
=
nufin,
[I]
/
INST_FIN
= instfin,
[R]
PRECISION
=
/
1.0E-3,
[DEFECT]
/
prec, [R]
SUBD_PAS
=/
1, [DEFECT]
/
subpas,
[I]
SUBD_PAS_MINI
=
submini,
[R]
COEF_SUBD_PAS_1
=/
1.,
[DEFECT]
/
coefsub, [R]
OPTI_LIST_INST
=
/
“INCR_MAXI”,
[DEFECT]
NOM_CHAM
= nomch,
[KN]
NOM_CMP
=
nomcmp, [kN]
VALE
=
valley
[R]
),
NEWTON
=_F (
PREDICTION
=
/
“TANGENT”,
[DEFECT]
/
“ELASTIC”,
/
“EXTRAPOL”,
/
“DEPL_CALCULE”,
EVOL_NOLI
=
evol_noli, [evol_noli]
STAMP
=
/
“TANGENT”,
[DEFECT]
/
“ELASTIC”
REAC_INCR
=
/
1, [DEFECT]
/
MF,
[I]
REAC_ITER
=
/
0, [DEFECT]
/
it,
[I]
REAC_ITER_ELAS
=
/
0, [DEFECT]
/
it,
[I]
PAS_MINI_ELAS
=
/
0, [DEFECT]
/
pasmini, [R]
),
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
6/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
RECH_LINEAIRE =_F (
RESI_LINE_RELA
=
/
1.E-1,
[DEFECT]
/
reslin,
[R]
ITER_LINE_MAXI
=
/
3
[DEFECT]
/
itelin
[I]
PAS_MINI_CRIT
=
/
0. [DEFECT]
/
pmicri
[R]
ITER_LINE_CRIT
=
/
20 [DEFECT]
/
itelic
[I]
RHO_MIN
=
/
1.E-2 [DEFECT]
/
rmin
[R]
RHO_MAX
=
/
1.E+1 [DEFECT]
/
rmax
[R]
RHO_EXCL
=/
9.E-3 [DEFECT]
/
rexc
[R]
),
PARM_THETA
=
/
1.,
[DEFECT]
/
theta,
[R]
PILOTING =_F (
TYPE
=
/
“DDL_IMPO”,
/
“LONG_ARC”,
/“ANA_LIM”,
/
“DEFORMATION”,
/
“PRED_ELAS”,
/
ALL =
“YES”,
[DEFECT]
/
GROUP_MA
= lgrma,
[l_gr_maille]
/
NET
=
lma, [l_maille]
/
NODE
=
No,
[node]
/
GROUP_NO
= grno,
[gr_noeud]
NOM_CMP =
nomcmp, [kN]
COEF_MULT
=
/
1.,
[DEFECT]
/
cmult,
[R]
ETA_PILO_R_MAX
=
etarmax,
[R]
ETA_PILO_R_MIN
=
etarmin,
[R]
ETA_PILO_MAX
=
etamax, [R]
ETA_PILO_MIN
=
etamin
[R]
PROJ_BORNES
=
/
“YES” [DEFECT]
/“NOT”
SELECTION
=
/
“NORM_INCR_DEPL”,
[DEFECT]
/
“ANGL_INCR_DEPL”,
/
“RESIDUE”,
),
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
7/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
SOLVEUR =_F (see the document [U4.50.01]
),
CONVERGENCE =_F (
/RESI_GLOB_RELA
=
1.E-6,
[DEFECT]
/
|
RESI_GLOB_MAXI
=
resmax,
[R]
|
RESI_GLOB_RELA
=
resrel,
[R]
|
RESI_REFE_RELA
=
resref,
[R]
SIGM_REFE
=
sigref, [R]
EPSI_REFE
=
sigref, [R]
FLUX_THER_REFE
=
sigref, [R]
FLUX_HYD1_REFE
=
sigref, [R]
FLUX_HYD2_REFE
=
sigref, [R]
ITER_GLOB_ELAS
=
/
25,
[DEFECT]
/
maxelas, [I]
ITER_GLOB_MAXI
=
/
10,
[DEFECT]
/
maglob,
[I]
STOP
=/
“YES”,
[DEFECT]
/
“NOT”,
RESI_INTE_RELA
=
/
1.E-6,
[DEFECT]
/
resint,
[R]
ITER_INTE_MAXI
=
/
10,
[DEFECT]
/
iteint,
[I]
ITER_INTE_PAS
=
/
0, [DEFECT]
/
itepas,
[I]
RESO_INTE
=
/
“IMPLICIT”,
[DEFECT]
/
“RUNGE_KUTTA_2”,
/
“RUNGE_KUTTA_4”,
),
CRIT_FLAMB
=_F (
NB_FREQ =
/
3, [DEFECT]
/
nbfreq,
[I]
CHAR_CRIT
=
/
(- 10,10),
[DEFECT]
/
intcc,
),
SENSITIVITY (see the document [U4.50.02]
),
FILING
=_F
(
/
LIST_INST
=
list_r8,
[listr8]
/
INST
=
l_r8,
[R]
/
PAS_ARCH
= npas,
[I]
PRECISION
=
/
1.E-3,
[DEFECT]
/
prec, [R]
/ARCH_ETAT_INIT
=
“YES”,
/
NUME_INIT
=
nuinit, [I]
DETR_NUME_SUIV
=
“YES”,
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
8/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
CHAM_EXCLU
=
|
“DEPL”,
|
“SIEF_ELGA”,
|
“VARI_ELGA”,
|
“VARI_NON_LOCAL”,
|
“LANL_ELGA”,
),
DISPLAY
=_F
(
/
LIST_INST
=
list_r8,
[listr8]
/
INST
=
l_r8,
[R]
/
PAS_ARCH
= npas,
[I]
UNIT =
/unit
[I]
LONG_R
=/
12 [DEFECT]
/
long_r
[I]
PREC_R =
/5 [DEFECT]
/
prec_r
[I]
LONG_I =
/6 [DEFECT]
/
long_i
[I]
NOM_COLONNE
=
|
“STANDARD”,
|
“MINIMUM”,
|
“ITER_NEWT”,
|
“INCR_TPS”,
|
“RESI_RELA”,
|
“RELA_NOEU”,
|
“RESI_MAXI”,
|
“MAXI_NOEU”,
|
“RESI_REFE”,
|
“REFE_NOEU”,
|
“RELI_ITER”,
|
“RELI_COEF”,
|
“PILO_PARA”,
|
“LAGR_ECAR”,
|
“LAGR_INCR”,
|
“LAGR_ITER”,
|
“MATR_ASSE”,
|
“ITER_DEBO”,
|
“CTCD_ITER”,
|
“CTCD_INFO”,
|
“CTCD_GEOM”,
|
“CTCD_NOEU”,
|
“CTCC_CONT”,
|
“CTCC_FROT”,
|
“CTCC_GEOM”,
INFO_RESIDU
=
“YES”,
[DEFECT]
“NOT”
),
OBSERVATION =_F (see the document [U4.53.01]
),
LAGR_NON_LOCAL =_F (
ITER_PRIM_MAXI
=
/
10,
[DEFECT]
/
iterprimmax,
[I]
RESI_PRIM_ABSO = resiprimab,
[R]
ITER_DUAL_MAXI
=
/
50,
[DEFECT]
/
iterdmax,
[I]
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
9/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
RESI_DUAL_ABSO
=
residabso, [R]
R
=
/
1000. [DEFECT]
/
rho
[R]
),
SOLV_NON_LOCAL =_F (
to see the document [U4.50.01]
),
INFORMATION
=
/
1, [DEFECT]
/2,
TITRATE
=
tx [kN]
);
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
10/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3 Operands
3.1 Operands
MODEL
/
CHAM_MATER
/
CARA_ELEM
MODEL = Mo
CHAM_MATER = chmat
CARA_ELEM = carac
These key words make it possible to inform:
· the name of the model (
Mo
) whose elements are the subject of mechanical calculation,
· the name of the material field (
chmat
) affected on the mesh. Attention, all meshs
model must be associated a material (if not fatal error with message little
clarify),
· the name of the characteristics (
carac
) of the elements of hull, beam, pipe, bars, cable,
and discrete elements affected on the model Mo. Obviously, this key word is optional: if
the model does not contain such elements, it is not useful; on the other hand, if the model
contains such elements, it is obligatory.
3.2 Word
key
EXCIT
EXCIT:
This key word factor makes it possible to describe with each occurrence a load (stresses and conditions
with the limits), and possibly a multiplying coefficient and/or a type of load.
3.2.1 Operands
CHARGE
CHARGE: CH
I
CH
I
is the mechanical loading (possibly comprising the evolution of a field of
temperature) specified with I
ème
occurrence of
EXCIT
.
One and only one load can comprise the evolution of a field of temperature, which will have
previously be defined thanks to the key word
TEMP_CALCULEE
control
AFFE_CHAR_MECA
.
Caution:
In a thermomechanical calculation, if the initial temperature is different from the temperature
of reference (given in the operator
AFFE_MATERIAU
), the field of deformation associated with
the initial moment can be incompatible and thus lead to a state of stresses and variables
interns associated not no one. If one uses a relation of behavior incremental (key word
factor
COMP_INCR
) and if one explicitly does not define a state of stresses and variables
interns initial (associate with a field of initial temperature different from the temperature from
reference), the internal variable and stress field calculated to the first increment
account will hold that only variation in temperature enters the initial moment and the first
moment, and not of the possible stresses of compatibility associated with the initial temperature.
To take this initial state hopes some, it should be given explicitly, for example thanks to
key words
SIGM
,
DEPL
,
VARI
and
VARI_NON_LOCAL
in
ETAT_INIT
.
To avoid such situations which can lead to computational errors, it is worth
to better begin a calculation while considering than it is necessary to start from a virgin state.
Caution:
If one carries out a calculation into axisymmetric and that one imposes nodal forces, these efforts
must be divided by 2 * pi (one works on a sector of 1 radian) compared to
real loadings. In the same way, if one wishes to calculate the resultant of the efforts, the result is
to multiply by 2 * pi to have the total resultant on the complete structure. In the same way in
plane stresses or in plane deformation, one works on a thickness unit: efforts
(on the thickness) applied must be divided by the thickness, the real efforts are obtained
by multiplying by the thickness the efforts “of calculation”.
Caution:
Loadings resulting from
AFFE_CHAR_CINE
are not usable with
STAT_NON_LINE
.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
11/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.2.2 Operand
FONC_MULT
FONC_MULT: F
I
F
I
is the multiplying function of the time of the loading specified with I
ème
occurrence of
EXCIT
.
The loading and boundary conditions for
N
occurrences of the key word factor
EXCIT
are:
=
=
N
I
I
I
CH
F
CH
1
For the conditions of Dirichlet, of course, only the specified value is multiplied by
F
I
.
By defect:
F
I
= 1.
Note:
The field of temperature is not multiplied by
F
I
.
3.2.3 Operand
TYPE_CHARGE
TYPE_CHARGE: tchi
By defect,
tchi
is worth
“FIXE_CSTE”
: that corresponds to a loading applied to
initial geometry and not controlled. It can however be a function, and, in particular, to depend on
time.
If
tch
I
is worth
“FIXE_PILO”
, the loading is always fixed (independent of the geometry) but
will be controlled thanks to the key word
PILOTING
[§3.11]. The loads controllable must result
of
AFFE_CHAR_MECA
or of
AFFE_CHAR_MECA_F
and not to be affected key word
FONC_MULT
.
One cannot control the loadings of gravity, the centrifugal force, the forces of Laplace,
thermal loadings or of initial or anelastic deformations, and conditions of
connection.
If
tch
I
is worth
“SUIV”
, the loading is known as “follower”, i.e. it depends on the value on
unknown factors: for example, pressure, being a loading applying in the normal direction
with a structure, depends on the geometry brought up to date of the aforementioned, and thus on displacements. One
following loading is revalued with each iteration of the algorithm of resolution. A loading
fix is revalued only at each new moment, and only if
CH
I
depends on time (defined in
AFFE_CHAR_MECA_F
and parameterized by the moment).
Currently the loadings which can be described as
“SUIV”
are the loading of
gravity for the element of
CABLE_POULIE
, pressure for modelings
3D
,
3d_SI
,
D_PLAN
,
D_PLAN_SI
,
AXIS
,
AXIS_SI
,
C_PLAN
,
C_PLAN_SI
and for all modelings
THM
(
3d_HHM
,
3d_HM
,
3d_JOINT_CT
,
3d_THH
,
3d_THHM
,
3d_THM
,
AXIS_HHM
,
AXIS_HM
,
AXIS_THH
,
AXIS_THHM
,
AXIS_THM
,
D_PLAN_HHM
,
D_PLAN_HM
,
D_PLAN_THH
,
D_PLAN_THHM
,
D_PLAN_THM)
and the centrifugal force in great displacements (key word
ROTATION
in
AFFE_CHAR_MECA
).
If
tchi
is worth
“DIDI”
then conditions of Dirichlet (imposed displacements, conditions
linear) will apply to the increment of displacement as from the moment given under
ETAT_INIT/NUME_DIDI
(by defect the moment of resumption of calculation) and not on displacement
total. For example for a displacement imposed (key word
DDL_IMPO
of
AFFE_CHAR_MECA
)
condition will be form:
U
U
D
-
=
0
where
U
0
is the displacement defined by
NUME_DIDI
and
not:
U
D
=
.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
12/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.3 Key word
SOUS_STRUC
For more precision concerning the use of substructures (elastic linear) in one
nonlinear structure, one will refer to documentation [U2.07.02]
SOUS_STRUC
This key word factor makes it possible to specify which are the loadings to be used for
substructures. In its absence, the loadings on under structures are null.
These loadings are added to the loadings “finite elements” which can be applied to
remain model.
CAS_CHARGE = nocas
nocas
is the name of the loading case to be used. See operator
MACR_ELEM_STAT
[U4.62.01].
/ALL = “YES”
This key word makes it possible to affect the loading
nocas
with all under structures of
model.
/MESH = l_mail
This key word factor makes it possible not to affect the loading
nocas
that with some
substructures.
3.4 Key words
COMP_INCR and COMP_ELAS
The syntax of these key words common to several controls is described in the document
[U4.51.11].
3.5 Word
key
VARI_COMM
VARI_COMM
:
Variables of controls which control the laws of behavior (as well as
temperature).
3.5.1 Operand
IRRA & CORROSION
/
IRRA
:
irr
Exposure fields.
/CORROSION
: corro
Fields of corrosion.
3.6 Word
key
ETAT_INIT
ETAT_INIT:
Initial State of reference selected. By defect, all the fields are identically null. This initial state
can be defined either by specifying each field of the initial state, or in extraction since one
concept of the type
evol_noli
preexistent.
The data of an initial state does not have a direction (and is not thus taken into account) only for the part of
field treated in incremental behavior (
COMP_INCR
); if the behavior is elastic
(
COMP_ELAS
) that does not have any angle of attack.
If one wants to take into account an initial state in elasticity, it is the key word
ELAS
located under
COMP_INCR
that it is necessary to use.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
13/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
Note:
If the user specified that the concept result is réentrant (by the reserved word
reuse
), key word ETAT_INIT is obligatory.
3.6.1 Operands
SIGM
/
VARI
/
DEPL
/
VARI_NON_LOCAL
/
|
SIGM
=
sig
| VARI
= vain
| DEPL
= depl
| VARI_NON_LOCAL = vanolo
Respectively, stress fields at the points of Gauss, variables intern at the points of
Gauss, of displacements to the nodes and nonlocal variables to the nodes (for models
not buildings) taken in an initial state. If one of these fields is not specified, it is taken null by defect. They
can for example be resulting from the control
CREA_CHAMP
, or to be read in one
file with format I-DEAS by the control
LIRE_RESU
(attention format MED only reads
fields with the nodes).
3.6.2 Operands
EVOL_NOLI
/
EVOL_NOLI
:
evol
Name of the concept of the evol_noli type from where will be extracted the initial state.
3.6.3 Operand
NUME_ORDRE
/
INST
/
NUME_DIDI
/NUME_ORDRE
= nuini
/
INST =
instini
Extraction of the initial mechanical state in evol starting from the number of filing
NUME_ORDRE
or of the moment of filing
INST
to carry out the continuation of calculation.
If
NUME_ORDRE
or
INST
are not filled, one takes the last filed number
existing in evol.
NUME_DIDI: nudidi
In the case of loadings of the differential type DIRICHLET (“DIDI”), one gives under NUME_DIDI
the number of filing of the mechanical state (displacement) which is used as reference for the application
of these boundary conditions (cf [§3.2.2]). By defect one takes the definite mechanical state under
NUME_ORDRE or INST.
3.6.4 Operand
INST_ETAT_INIT
INST_ETAT_INIT: istetaini
One can associate a value of moment
istetaini
in this initial state.
By defect:
· when the initial state is defined by the data of the fields, an associated moment ago.
· when the state is given by a concept
evol_noli
, it is about the moment in the precedent
calculation (
istetaini = instini
).
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
14/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
With - Simple example (by defect)
LIST1 = DEFI_LIST_REEL (
=0 BEGINNING.,
INTERVAL =_F (UNTIL = 4., A =4 NUMBER)),
U = STAT_NON_LINE (INCREMENT =_F (LIST_INST =LIST1)) ,
LIST2 = DEFI_LIST_REEL (
=4 BEGINNING.,
INTERVAL =_F (UNTIL = 10., A =6 NUMBER)),
U = STAT_NON_LINE
(reuse=U,
INCREMENT
=_F (LIST_INST
=LIST2),
ETAT_INIT
=_F (EVOL_NOLI
=U))
,
1
er
STAT_NON_LINE
: carry out calculation for the moments 1, 2, 3 and 4s.
2
Nd
STAT_NON_LINE
: carry out calculation for the moments 5, 6, 7, 8, 9 and 10s, the initial state
corresponding to time 4s.
B - Example to show the interest of
INST_ETAT_INIT
(two different lists of moments)
LIST1 = DEFI_LIST_REEL (
=0 BEGINNING.,
INTERVAL =_F (UNTIL = 10., A =10 NUMBER)),
U = STAT_NON_LINE (INCREMENT =_F (LIST_INST =LIST1)) ,
LIST2 = DEFI_LIST_REEL (
=20 BEGINNING.,
INTERVAL =_F (UNTIL = 30., A =10 NUMBER)),
U = STAT_NON_LINE
(reuse=U
INCREMENT
=_F (LIST_INST
=LIST2),
ETAT_INIT
=_F (EVOL_NOLI
=U,
INST_ETAT_INIT
=
20.))
,
1
er
STAT_NON_LINE
: carry out the calculation of moments 1 with 10s.
2
Nd
STAT_NON_LINE
: carry out the calculation of moments 21 with 30s, the initial state corresponding to the moment
t=10s of the 1
er
STAT_NON_LINE
(by defect
INST=10
.). This initial state corresponds for this 2
Nd
STAT_NON_LINE
at the moment t=20s. (
INST_ETAT_INIT=20.)
.
C - Example to show the interest of
INST_ETAT_INIT
(practical when the cyclic one is made)
LIST1 = DEFI_LIST_REEL (
=0 BEGINNING.,
INTERVAL =_F (UNTIL = 10., A =10 NUMBER)),
U1 = STAT_NON_LINE (INCREMENT =_F (
LIST_INST =LIST1)) ,
U2 = STAT_NON_LINE (INCREMENT =_F (
LIST_INST =LIST1),
ETAT_INIT
=_F (
EVOL_NOLI
=U1,
INST_ETAT_INIT
=
0.))
,
1
er
STAT_NON_LINE:
carry out the calculation of moments 1 with 10s.
2
Nd
STAT_NON_LINE:
carry out the calculation of moments 1 with 10s, the initial state corresponding to the moment
t=10s of the 1
er
STAT_NON_LINE
(by defect
INST=10
.). This initial state corresponds for this 2
Nd
STAT_NON_LINE
at the moment t=0s. (
INST_ETAT_INIT: 0.)
.
3.6.5 Operand
PRECISION/CRITERION
Cf [U4.71.00].
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
15/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.7 Word
key
INCREMENT
INCREMENT:
Defines the intervals of time taken in the incremental method.
The moments thus defined have physical direction only for relations of behavior where
time intervenes explicitly (viscoelastic or viscoplastic for example). In the others
cases, they allow only indicer the increments of load and to parameterize the evolution of one
possible field of temperature.
3.7.1 Operands
LIST_INST/EVOLUTION
LIST_INST: litps
The moments of calculation are those defined in the concept litps by operator DEFI_LIST_REEL
[U4.34.01].
EVOLUTION
: /“CHRONOLOGICAL” [DEFECT]
/
“RETROGRESSES”
/
“WITHOUT”
The “CHRONOLOGICAL” key word makes it possible to check if the list of moments given by the user is
strictly increasing (so not an error message is transmitted).
The “RETROGRADE” key word makes it possible to reverse the list of moments given by the user and of
to check that after this operation, it is well strictly decreasing.
There is no checking when one specifies an evolution “WITHOUT”.
3.7.2 Operands
NUME_INST_INIT/INST_INIT/NUME_INST_FIN/INST_FIN
/
NUME_INST_INIT = nuini
/
INST_INIT
= instini
The initial moment of the calculation (which thus (Re) is not calculated) is indicated either by its value
(INST_INIT), that is to say by its sequence number in the list of moments litps (NUME_INST_INIT).
To be able to reach by value, it is necessary that the list is ordered (EVOLUTION:
“CHRONOLOGICAL” or “RETROGRESSES”).
In the absence of key words INST_INIT or NUME_INST_INIT, the defect is calculated
following manner:
· if an initial state is specified (operand ETAT_INIT) and if it definite one moment
corresponding (by EVOL_NOLI or INST_ETAT_INIT) then the initial moment is that
defined by the initial state,
· if there is no initial state (operand ETAT_INIT) or that it does not define a moment
corresponding (the fields are given in ETAT_INIT without specifying
INST_ETAT_INIT), then one takes the first moment of the list of moments litps
(NUME_INST_INIT:0), or the last when the evolution is retrograde.
· In the event of filing (see keyword FILING), the initial moment in continuation is it
last pitch filed and not that defined in INST_INIT.
/
NUME_INST_FIN =
nufin
/
INST_FIN
=
instfin
The final moment (last calculated pitch) is indicated same manner that the initial moment (either
NUME_INST_FIN, is INST_FIN), except that it is not possible to refer to the moment
initial state.
Caution: with an evolution RETROGRAGE, INST_INIT > INST_FIN.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
16/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
With - Simple example (by defect)
LIST = DEFI_LIST_REEL (=0 BEGINNING.,
INTERVAL =_F (UNTIL = 10., A =10 NUMBER)),
U = STAT_NON_LINE (
INCREMENT =_F (LIST_INST =LIST,
INST_FIN
=4.))
,
U = STAT_NON_LINE (
reuse=U,
INCREMENT
=_F
(
LIST_INST
=LIST),
ETAT_INIT
=_F
(EVOL_NOLI
:U))
,
1
er
STAT_NON_LINE:
carry out calculation for the moments 1, 2, 3 and 4s.
2
Nd
STAT_NON_LINE:
carry out calculation for the moments 5, 6, 7, 8, 9 and 10s, the initial state
corresponding to time 4s. (by defect
INST_INIT=INST_ETAT_INIT=INST=4.
).
B - Example to show the interest of
INST_INIT
LIST = DEFI_LIST_REEL (=0 BEGINNING.,
INTERVAL = _F (UNTIL = 10., A =10 NUMBER)),
U = STAT_NON_LINE
(INCREMENT = _F (LIST_INST = LIST,
INST_FIN
=
4.))
,
U = STAT_NON_LINE
(reuse = U,
INCREMENT
=_F
(
LIST_INST
=LIST,
INST_INIT =8.),
ETAT_INIT
=_F
(
EVOL_NOLI
=U))
,
1
er
STAT_NON_LINE:
carry out the calculation of moments 1 with 4s.
2
Nd
STAT_NON_LINE:
carry out calculation for moments 9 and 10s (does not do anything for t=5, 6, 7 and 8s),
the initial state corresponding to time t=4s (by defect
INST=4
.).
3.7.3 Operand
PRECISION
PRECISION: prec
Cf [U4.71.00]
3.7.4 Operand
SUBD_PAS/SUBD_PAS_MINI/COEF_SUBD_PAS_1
SUBD_PAS
=
subpas
SUBD_PAS_MINI = submini
COEF_SUBD_PAS_1
= coefsub
Allows to carry out an automatic recutting of the pitch of time when the algorithm of Newton
do not converge.
The pitch of time is redécoupé in
subpas
under pitch. By defect there is no recutting
(
subd_pas
: 1). The automatic subdivision stops when the new pitches created are more
small that
SUBD_PAS_MINI
. The new pitches created are of identical size, except the first
who is equal to this size multiplied by
COEF_SUBD_PAS_1
(by defect 1). This allows best
to take into account the problems of discharge of the structure (change of tangent matrix)
without using the elastic matrix (
PREDICTION:“ELASTIC”
or
STAMP: “ELASTIC”
under the operand
NEWTON
).
When a pitch of time was redécoupé several times (let us call N the number of times where one has
proceeded to a subdivision of the same pitch), the following pitch is automatically subdivided (n-1) time,
this to avoid, in the event of convergence difficult to try a pitch of too important time.
Notice concerning the key word
CUT OUT
under
SOLVEUR
:
During elastoplastic calculation of buckling, it can happen that the tangent matrix of the system
that is to say singular during iterations of Newton. By redécoupant the pitch of time, one can
to pass these hard points. Under operand SOLVEUR, the key word CUTS OUT under
STOP_SINGULIER is used to manage these hard points. It is then necessary to inform the words
keys relating to recutting so that the method CUTTING is activated.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
17/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.7.5 Operand
OPTI_LIST_INST/NOM_CHAM/NOM_CMP/VALE
OPTI_LIST_INST =
“INCR_MAXI”
[DEFECT]
NOM_CHAM
=
“TEMP”
[DEFECT]
NOM_CMP =
“TEMP”
[DEFECT]
VALE =
vale
These operands have interest only when a thermomechanical calculation is carried out. Allows to create
if need be a new list of mechanical pitch of time so that, between each increment of
time, the increment of temperature is lower than a value given by the user and indicated
by the key word
VALE
.
The creation of this new list is done in the following way:
· List moments initial in mechanics: Ti
· List moments thermics:
· New final list of moments mechanical to create if need be: Tf
· One inserts between each interval of the basic list mechanical Ti, the thermal moments
include in this interval. One then recovers for each interval a list of moments
=
[
0
,
1
,
2
,
NR
]
· Construction of the list final Tf
- Initialization:
F
=
0
- 1
er
Test:
If T (
J
) - T (
F
) > value with T (T) the temperature at time T and
F
the last moment
inserted in the new Tf list, then one keeps in the new list Tf, the moment
j-1
- 2
ème
Test:
If T (
J
) - T (
j-1
) > value then one redécoupe uniformly this interval in way
to satisfy the condition on the increment of temperature.
Example: IF T (
) = [T (
1
) = 20°C, T (
2
) = 30°C, T (
3
) = 55°C, T (
4
) = 65°C] with
VALE
= 15°C
Initialization:
F
=
1
Interval 1:
1
er
Test = 2
ème
Test: T (
2
) - T (
1
) = 10°C < 15 thus one Tf = [
1
]
Interval 2:
1
er
Test: T (
3
) - T (
F
) = 35°C > 15 thus one has Tf = [
1
,
2
] and
F
=
2
2
ème
Test: T (
3
) - T (
2
) = 25°C > 15 thus one Tf = [
1
,
2
, T
3
such as T (T
3
) = 42.5°C,
3
] and
F
=
3
Interval 3:
1
er
Test = 2
ème
Test: T (
4
) - T (
3
) = 10°C < 15°C
from where the following final list:
Tf = [
1
,
2
, T
3
such as T (T
3
) = 42.5°C,
3
,
4
]
3.8 Word
key
NEWTON
NEWTON
:
Specify the characteristics of the method of resolution of the nonlinear incremental problem
(method of NEWTON-RAPHSON).
3.8.1 Operand
PREDICTION
PREDICTION
=
/“TANGENT”
/“ELASTIC”
/“EXTRAPOL”
/“DEPL_CALCULE”
The purpose of the phase of prediction (cf [R5.03.01]) is to calculate an estimate of the field of
displacements in order to allow the method of NEWTON more quickly to converge.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
18/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
When the key word misses, it is the tangent matrix of speed (option
RIGI_MECA_TANG
in the file .mess) which is used if one chose for the method of NEWTON one
STAMP:“TANGENT”
, and it is the elastic matrix (option
RIGI_MECA
in the file .mess)
who is used if one chose
STAMP:“ELASTIC”
.
/“TANGENT”
One uses the tangent matrix of the problem of speed (option
RIGI_MECA_TANG
in the file
.mess).
/“ELASTIC”
One uses the elastic matrix (option
RIGI_MECA
in the file .mess).
/“EXTRAPOL”
One calculates the estimate of the increment of displacement starting from the total increment obtained like
solution with the pitch of previous time (balanced by the report/ratio of the pitches of time). One projects this
estimate on the whole of the fields kinematically acceptable (i.e satisfying them
boundary conditions of DIRICHLET) according to the standard given by the elastic matrix, which must
thus to be calculated. This functionality is interesting in the case of the use of diagrams
of explicit integration local of type RUNGE-KUTTA which does not provide a tangent matrix:
in this case the method of NEWTON uses an elastic matrix, but the iteration count
necessary can be high. The use of extrapolation can improve the performances.
/“DEPL_CALCULE”
Allows to propose like displacement for the prediction with each pitch of time, it
displacement given by a mechanical history specified under the key word
EVOL_NOLI
([§3.8.3]).
Utility:
· let us suppose that one carries out the first calculation with a coarse mesh. One wishes to carry it out
even calculation but on a finer mesh. One can suppose that the solution in displacement
for this second calculation is not distant from that of the first calculation and thus only one good
prediction of displacement for this second calculation is the projection of displacements of calculation
1 on the nodes of the new mesh (the projection of displacements on the new mesh
must be realized beforehand with the operator
PROJ_CHAMP
[U4.72.05]). This key word allows
to carry out this mode of prediction.
· that makes it possible to reduce the place memory and to preserve these results for a continuation
later. For a large calculation, one can store only displacements with all them
moments with formats IDEAS or MED in
IMPR_RESU.
If one wants to recompute the stresses and
internal variables, one is done
LIRE_RESU
with the adequate format then one uses
DEPL_CALCULE
with
ITER_GLOB_MAXI: 0
(only one iteration is carried out) and
STOP:NOT
(there is not
convergence, one does not check balance). It is however necessary for reasons of
syntax to give a loading (to avoid the loadings dirichlet which impose one
linear resolution) as well as a criterion of convergence, even if this information is not
takings into account.
3.8.2 Operand
STAMP
STAMP =
/
“TANGENT”
REAC_INCR
=
/
1
[DEFECT]
/MF
REAC_ITER
=
/
0
[DEFECT]
/it
The matrix used for the total iterations of the method is the tangent matrix [R5.03.01]
who is revalued all them
MF
increments of time (
MF
positive or null) and all them
it
iterations of
NEWTON for an increment of time given (precisely to the iterations of number
it
,
2it
,
3it
…). Thus with the first iteration of NEWTON, one reassembles the tangent matrix only if
it
1 is worth: if not one keeps the matrix used in the phase of prediction. By convention if
it
0 are worth
the matrix is not revalued during all the pitch of time.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
19/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
PAS_MINI_ELAS
=
/
0.
[DEFECT]
/
pasmini [R]
REAC_ITER_ELAS =/
0
[DEFECT]
/
it
[I]
Allows to pass from the tangent matrix to the matrix of discharge (i.e by considering that not
linearities do not evolve/move) when the pitch of time is or becomes (by recutting) lower than
pasmini. This matrix of discharge is the elastic matrix for the models of behavior
of plastic type; for the models of damage it is identified with the secant matrix.
As convergence with the elastic matrix is slower than that with the matrix
tangent, the key word
ITER_GLOB_ELAS
under the key word factor
CONVERGENCE
allows to define
an iteration count maximum specific to the use of the matrix elastic and different from that
associated the use of the tangent matrix.
One can define a frequency of reactualization of the matrix of discharge with the key word
REAC_ITER_ELAS
(analog of
REAC_ITER
). If the matrix of discharge does not depend on the state
of deformation, to take
REAC_ITER_ELAS
= 0 (since it will be the same one during
iterations).
Utility:
This option can be useful when the automatic recutting of the pitch of time (cf [§ 3.7.4])
converge a calculation is not enough to make. For example, in the case of lenitive laws, the matrix
tangent can become singular and it is thus to better use the elastic matrix to converge.
/
“ELASTIC”
The matrix used corresponds to the elastic design: it is evaluated only once at the moment
initial, at the beginning of algorithm.
This “elastic” matrix is calculated by using the YOUNG modulus given under the key word
ELAS
of the operator
DEFI_MATERIAU
, and not the slope at the origin of the traction diagram
data under the key word
TRACTION
(and which is useful, it, in the expression of the relation of
behavior).
3.8.3 Operand
EVOL_NOLI
EVOL_NOLI: evol_noli
Name of the concept of the evol_noli type which will be useful in the prediction by DEPL_CALCULE.
3.9 Word
key
RECH_LINEAIRE
RECH_LINEAIRE:
Linear search can make it possible to improve convergence of the method of Newton
(Cf [R5.03.01] for more details).
Caution:
It is disadvised using linear search with deformations GREEN_GR for
modelings COQUE_3D and in the presence of contact.
3.9.1 Operand
RESI_LINE_RELA
/
ITER_LINE_MAXI
RESI_LINE_RELA =/
1.E-1 [DEFECT]
/
reslin
ITER_LINE_MAXI
=
/
3
[DEFECT]
/
itelin
They are the parameters of linear search. The maximum iteration count is given
itelin
to carry out and precision
reslin
to reach to carry out the convergence of
linear search. It is advised not to use linear search with contact.
It is not necessary to specify a precision nor an iteration count very high, the practice
showing that 2 or 3 iterations of linear search are sufficient. One can thus be satisfied
to ask 3 iterations with the precision by defect.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
20/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.9.2 Operand
PAS_MINI_CRIT
/
ITER_LINE_CRIT
PAS_MINI_CRIT
=
/
0.
[DEFECT]
/
pmicri
[R]
ITER_LINE_CRIT =
/
20
[DEFECT]
/
itelic
[I]
At the time of pitch of time when convergence is delicate, one can want to increase the number
maximum of iterations of linear search. It is what the key words allow
PAS_MINI_CRIT and ITER_LINE_CRIT. When the pitch of time (directly fixed by the user
or consequence of cuttings of pitch of time) becomes lower than the value pmicri, the number
iterations of linear search for search passes from itelin (well informed by
ITER_LINE_MAXI) with itelic (informed by ITER_LINE_MAXI)
3.9.3 Operands
RHO_MIN
/
RHO_MAX
/
RHO_EXCL
RHO_MIN =
/
1.E-2
[DEFECT]
/
rmin [R]
RHO_MAX =
/
1.E+1
[DEFECT]
/
rmax [R]
RHO_EXCL
=
/
9.E-3
[DEFECT]
/
rexc [R]
These key words fix interval I of linear search, in the form
:
[
] [
]
rexc
rexc
R
R
I
,
max
min,
-
-
=
.
3.10 Operand
PARM_THETA
PARM_THETA
=
/
1.
[DEFECT]
/
theta
For modelings
THM
, the argument
theta
is the parameter of the theta-method used for
to solve the evolutionary equations of thermics and hydraulics (cf [R5.03.60] for more
details). Its value must lie between 0 (explicit method) and 1 (method completely
implicit).
For the laws of behavior
ROUSS_VISC
,
ASSE_COMBU
,
ZIRC_CYRA2
and
ZIRC_EPRI
,
the argument
theta
is used for integration of the law of behavior (for the model
ASSE_COMBU
, it
is used to integrate the law of Lemaitre in 1D). It can take values 0.5 (semi-implicit) or 1
(implicit).
3.11 Word
key
PILOTING
PILOTING:
When intensity
of part of the loading is not known a priori (loading known as of
reference defined in
AFFE_CHAR_MECA
or
AFFE_CHAR_MECA_F
with load of the type
FIXE_PILO
), the key word
PILOTING
allows to control this loading via a node
(or node groups) on which one can impose various modes of piloting (key word
TYPE
).
Caution:
With FIXE_PILO, one cannot use for the loading of reference the key word
FONCT_MULT.
Caution:
When the loading of reference is defined by AFFE_CHAR_MECA_F, this loading can
to be a function of the variables of space but not of time.
Caution:
The key word PILOTING is interdict with the contact.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
21/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.11.1 Operand
TYPE
TYPE:
/“DDL_IMPO”
/“LONG_ARC”
/“ANA_LIM”
/“DEFORMATION”
/“PRED_ELAS”
It is the type of piloting carried out. Five modes of piloting are available (cf [R5.03.80] for
more details):
/“DDL_IMPO”
Allows to impose a given value of increment of displacement (only one component
I
possible) in a single node
No
(or of a group of nodes comprising one node).
With each increment of time, one seeks the amplitude
loading of reference which
will allow to satisfy the following incremental relation:
C
U No
T
mult
I
()
=
/“LONG_ARC”
Allows to control the intensity
loading of reference by the length (X-coordinate
curvilinear) of the response in displacement of a group of nodes (to be used for example
when one wants to control the buckling of a test-tube). The following relation is checked:
C
U
T
mult
=
with
U
U
N C
C
N
=
,
2
where N are the nodes of piloting and C the components of the displacement of the nodes
considered. Even if the group of node of piloting is tiny room to only one node, it is necessary when
to even use
GROUP_NO
.
/“ANA_LIM”
This mode of piloting is specific to the calculation of load limits (law
NORTON_HOFF
) by
approach kinematic (cf [R7.07.01] for more detail). If F indicates the loading
assembled controlled,
TYPE
_
CHARGE
= `
FIXED
_
PILO
'
, then the function of piloting is written
simply:
1
.
)
P (
=
= U
F
U
Except for the calculation of limiting load, this functionality is not of interest a priori.
For this mode of piloting, no other key word is to be specified.
Note:
The use of lenitive laws of behavior can lead to snap backs
brutal which makes delicate the course of calculation. Two modes of piloting
following cures it (cf [R5.03.80] for more detail).
/“DEFORMATION”
DEFORMATION
guarantees that at least a point of Gauss of the structure sees its deformation
to evolve/move in a monotonous way. The relation is checked:
T
C
mult
=
-
-
)
(
max
Gauss
of
not
This mode of piloting is valid for all the laws of behavior including into large
deformations
SIMO_MIEHE
.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
22/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
/“PRED_ELAS”
PRED_ELAS
ensure that at least a point of Gauss of the structure left the threshold of elasticity
linearized
F
préd élas
-
of a quantity
T/C
mult
. The relation is checked:
T
F
C
élas
préd
mult
=
-
)
(
max
Gauss
of
not
This mode of piloting is valid only for the laws
ENDO_FRAGILE
(
with
version
local and two nonlocal versions),
ENDO_ISOT_BETON
and ENDO_ORTH_BETON (with
local version and the nonlocal version),
BARENBLATT
and
BETON_DOUBLE_DP
.
Use Attention:
When one wants to use these the last two modes of piloting, it is essential to make a first
STAT_NON_LINE
without the key word
PILOTING
to start the problem and to obtain an initial state
-
different from zero (if not divide check for piloting by increment of deformation). One carries out
after a recovery starting from this initial state not no one and one use piloting.
Moreover, the resolution of the two preceding equations makes it possible to obtain the intensity of the loading
unknown factor. In certain cases, the solution of these equations can lead to several solutions for
intensity. One then chooses always the solution which is closest to
-
. This is why, when one
wants to impose an alternated loading, one is obliged with each change of sign of the loading of
to carry out a first
STAT_NON_LINE
without the key word
PILOTING
in order to obtain an initial state
-
of
traction or of compression. One carries out then a second
STAT_NON_LINE
in continuation from
the preceding initial state with the key word
PILOTING
.
Note:
DEFORMATION
and
PRED_ELAS
are not available for the elements of structures.
3.11.2 Operands
NODE/GROUP_NO
/
NODE = No
/
GROUP_NO
=
grno
One gives the name of the node or the name of group of nodes on which one will impose it
piloting. To use only with `
DDL_IMPO
'or `
LONG_ARC
'.
For `
DDL_IMPO'
, if the operand is used
GROUP_NO
, the group of nodes in question
must contain that only one node. For `
LONG_ARC'
, one only uses
GROUP_NO
(which can
if required to contain one node).
3.11.3 Operands ALL/MESH/GROUP_MA
/
ALL =
“YES”
[DEFECT]
/
GROUP_MA
=
lgrma
/
NET
=
lma
One gives the meshs or groups of meshs being used to control calculation. To use only with
DEFORMATION
or
PRED_ELAS
. Interesting to reduce the resolution of the equations of these
three modes of pilotings.
3.11.4 Operand
NOM_CMP
NOM_CMP
:
nomcmp
It is the name of the component (corresponding to the degree of freedom
I
) used for piloting
(“DX” for example). To use only with “DDL_IMPO” or “LONG_ARC”.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
23/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.11.5 Operand
COEF_MULT
COEF_MULT
: cmult
It is the value (noted
C
mult
in the formula of definition) by which one multiplies the degree of
freedom used for piloting. By defect, this value is worth 1. With not using with
ANA_LIM
.
Example with
DDL_IMPO:
Let us suppose that one wants to know the limiting load of a structure.
The loading imposed on the structure is the pressure of unknown intensity (P=
* value of
Px reference) on the group of mesh A. to find the load limits P
limit
, one will control it
displacement of node NO1. It is wanted that final displacement according to X of this node is equal to 2.
(either according to the list of moments of the pitches of 0.2, or a coefficient cmult=1/0.2=5.)
PRESSURE = AFFE_CHAR_MECA (CLOSE = (GROUP_MA =A, PX = 1.0)) ,
LIST =
DEFI_LIST_REEL (=0 BEGINNING.,
INTERVAL
=_F (UNTIL
=
10,
NUMBERS
=10),
RESU =
STAT_NON_LINE (
EXCIT =
_F (
CHARGE = PRESSURE,
TYPE_CHARGE
=
“FIXE_PILO”),
PILOTING
=_F (
TYPE = “DDL_IMPO”,
NODE = NO1,
NOM_CMP
=
“DX”,
COEF_MULT
=
5.))
,
In the fichier.resu, the value of
calculation will be displayed at every moment. To know
charge limit, it is enough to make P
limit
=
* Px. (Here Px is worth 1 thus one has the limiting load directly).
If one imposes on the structure a pressure P close to the limiting load without using piloting, it
calculation will not converge if one is close to the limiting load.
3.11.6 Operand
ETA_PILO_R_MAX/ETA_PILO_R_MIN
ETA_PILO_R_MAX =
etarmax,
[R]
ETA_PILO_R_MIN =
etarmin,
[R]
These two key words make it possible to specify the interval of awaited values of piloting.
principle of operation is as follows: with each iteration of Newton, if values are found
of piloting in the interval
[
]
max
min, etar
etar
, all values of piloting apart from this
interval are not considered. On the other hand, if no value of piloting is found in
this interval, all the values of piloting are preserved.
If one does not specify values, it is
-
for etarmin and
+
for etarmax.
A possible use of this interval is as follows. one wishes for example, to control a pressure
some share on the structure and one expects to keep this positive pressure. By fixing etarmin at
0, that make it possible to preserve only the positive values of piloting, if at least one is found
positive value of piloting at the time of the resolution of piloting.
3.11.7 Operand
ETA_PILO_MAX/ETA_PILO_MIN
ETA_PILO_MAX:
etamax
Stop of calculation when the parameter of piloting reaches the value given etamax.
ETA_PILO_MIN:
etamin
Allows to stop calculation when the parameter
ETA_PILOTAGE
reached this minimal value
etamin (for lenitive models, makes it possible to stop calculation when the structure is
sufficient softened).
Caution:
With law ENDO_ISOT_BETON, these two words key are obligatory.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
24/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.11.8 Operand
PROJ_BORNES
PROJ_BORNES
=
/
“YES” [DEFECT]
/“NOT”
In the event of going beyond of the interval (etamin, etamax), the user can indicate if he wants
to project the value of piloting on (etamin, etamax).
With PROJ_BORNE=' OUI', projection will be carried out (if eta>etamax - > eta=etamax; if
eta<etamin - > eta=etamin), which allows, in the event of convergence to stop calculation
precisely on etamin or etamax.
With PROJ_BORNE=' NON', one does nothing, therefore calculation will stop, in the event of convergence, with
a value higher than etamax or lower than etamin.
3.11.9 Operand
SELECTION
/SELECTION
=/“NORM_INCR_DEPL”, [DEFECT]
/
“ANGL_INCR_DEPL”,
/
“RESIDUE”,
This operand makes it possible to select the method allowing for choice of the value of piloting
if several solutions are provided by the resolution of piloting.
“NORM_INCR_DEPL” makes it possible to select the value of piloting by the smallest standard of
the increment of displacement on the pitch of time considered.
“ANGL_INCR_DEPL” makes it possible to select the value of piloting by the smallest angle enters
the displacement obtained for the pitch of current time and the displacement obtained for the pitch of
previous time.
“RESIDUE” makes it possible to select the value of piloting leading to the smallest residue.
3.12 Word
key
SOLVEUR
The syntax of this key word common to several controls is described in the document [U4.50.01].
3.13 Word
key
CONVERGENCE
CONVERGENCE:
If none of the two operands following is present, then all occurs like if:
RESI_GLOB_RELA = 1.E-6
.
3.13.1 Operand
RESI_GLOB_RELA/
RESI_GLOB_MAXI
| RESI_GLOB_RELA
=
resrel
The algorithm continues the total iterations as long as:
Max
,…,
_
Max
F
I
Nb ddl
in
=
>
1
resrel
L
where
F
N
is the residue of the iteration
N
and
L
the vector of the imposed loading and the reactions
supports (cf [R5.03.01] for more details).
When the loading and the reactions of support become null, i.e. when
L
is
no one (for example in the case of a total discharge), one passes from the criterion of convergence
relating to the absolute criterion of convergence RESI_GLOB_MAXI. This operation is transparent
for the user (message of alarm emitted in the file .mess). When the vector
L
becomes again different from zero, one passes by again automatically with the relative criterion of convergence
RESI_GLOB_RELA.
If this operand misses, the test is carried out with the default value, except if
RESI_GLOB_MAXI is present.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
25/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
| RESI_GLOB_MAXI
=
resmax
The algorithm continues the total iterations as long as:
Max
,…,
_
F
I
Nb ddl
in
=
>
1
resmax
where
F
N
is the residue of the iteration
N
(Cf [R5.03.01] for more details).
If this operand misses, the test is not carried out.
If
RESI_GLOB_RELA
and
RESI_GLOB_MAXI
both are present, the two tests are
carried out.
3.13.2 Operand
RESI_REFE_RELA
|
RESI_REFE_RELA =
resref,
[R]
SIGM_REFE
=
sigref, [R]
EPSI_REFE
=
epsref, [R]
FLUX_THER_REFE =
fthref, [R]
FLUX_HYD1_REFE =
fh1ref, [R]
FLUX_HYD2_REFE =
fh2ref, [R]
This operand results in estimating the convergence of the algorithm of Newton in the manner
following (cf [R5.03.01] for more details). From the stress of reference
sigref
(and/or a deformation of reference
epsref
if one uses nonlocal laws with gradient of
deformation, and/or a heat flux of reference
fthref
in a case THM, and/or two
hydrous references of flow
fh1ref
and
fh2ref
in a case HHM), one calculates a reference
of residue
Fref
(a of the same vector length than the vector residue). Convergence will be
carried out if and only if:
[
]
ref.
I
F
ddl
Nb
I
N
I
resref
_
,…,
1
F
<
3.13.3 Operand
ITER_GLOB_MAXI
ITER_GLOB_MAXI =/10
[DEFECT]
/
maglob
Maximum iteration count carried out to solve the total problem at every moment
(10 per defect). This test is always carried out.
3.13.4 Operand
ITER_GLOB_ELAS
ITER_GLOB_ELAS =/25
[DEFECT]
/
maxelas
Maximum iteration count carried out with the elastic matrix when the word is used
key
PAS_MINI_ELAS
key word factor
NEWTON
(see [§3.8.2]) .pour to solve the problem
total at every moment (25 per defect).
It is reminded the meeting that
PAS_MINI_ELAS
allows to pass from the tangent matrix to the matrix
rubber band when the pitch of time is or becomes (by recutting) lower than one
certain value specified under
PAS_MINI_ELAS
.
3.13.5 Operand
STOP
STOP
=
/
“YES”
[DEFECT]
If one of the criteria of total convergence chosen is not checked afterwards
maglob
iterations, then the program stops (the preceding results are backed up).
/
“NOT”
If
maglob
is insufficient to check the criteria of convergence given by
the user, one passes nevertheless at the next moment. Use to be avoided.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
26/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.13.6 Operands
RESI_INTE_RELA
/
ITER_INTE_MAXI
RESI_INTE_RELA =/
1.E-6
[DEFECT]
/
resint
ITER_INTE_MAXI =/
10
[DEFECT]
/
iteint
In the majority of the relations of behavior, a nonlinear equation or a system
nonlinear must be solved locally (in each point of GAUSS). These
operands (residue and a maximum number of iterations known as intern) are used to test
the convergence of this iterative algorithm of resolution. For more details, to refer to
reference material, for example with the document [R5.03.02]. These operands
are useless with the behaviors
ELAS
,
VMIS_CINE_LINE
,
VMIS_ECMI, LINE
,
VMIS_ECMI_TRAC
,
VMIS_ISOT_LINE
,
VMIS_ISOT_TRAC
,
VISC_ISOT_LINE,
VISC_ISOT_TRAC, BARENBLATT,
NORTON_HOFF
,
DIS_CONTACT
,
DIS_CHOC
,
ARM
,
ASSE_CORN
,
DIS_GOUJ 2e_PLAS
,
DIS_GOUJ 2e_ELAS
,
VMIS_ASYM_LINE
,
GRILL_ISOT_LINE
,
GRILL_CINE_LINE
,
GRILL_PINTO_MEN
,
PINTO_MENEGOTTO
,
GRANGER_FP
and
GRANGER_FP_V
(except stress planes),
BAZANT_FD
and all them
relations META_XXX.
3.13.7 Operand
ITER_INTE_PAS
ITER_INTE_PAS
=
0
[DEFECT]
itepas
Redécouper locally the pitch of time allows to facilitate the integration of the relation of
behavior at the points of GAUSS (for the relations of
CHABOCHE
,
VISC_TAHERI
,
LMARC
,
LAIGLE
,
MONOCRYSTAL
,
ROUSS_PR
,
ROUSS_VISC
,
CJS
and
BETON_DOUBLE_DP
). If
itepas
is worth
0, 1 or - 1 it does not have there recutting. If
itepas
is positive, one redécoupe systematically it
no time locally in
itepas
small pitches of time before carrying out the integration of
relation of behavior. If
itepas
is negative, recutting in |
itepas
| small pitches of
time is carried out only in the event of nonlocal convergence.
3.13.8 Operand
RESO_INTE
RESO_INTE
=/“IMPLICIT”
[DEFECT]
/
“RUNGE_KUTTA_2”
/
“RUNGE_KUTTA_4”
Allows to specify the type of diagram of integration to solve the system of equations not
linear formed by the equations constitutive of the models of behavior to variables
interns:
· models
POLY_CFC
and POLYCRYSTAL are treated only by the explicit diagram
RUNGE-KUTTA of command 2,
· two models
VMIS_POU_LINE
and
VMIS_POU_FLEJOU
can be treated by
two diagrams
IMPLICIT
and
RUNGE_KUTTA_4
,
· two models MONOCRYSTAL and VENDOCHAB can be treated by both
implicit schemes and RUNGE_KUTTA_2,
· the other models use the diagram
IMPLICIT
.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
27/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.14 Key word
CRIT_FLAMB
CRIT_FLAMB
=_F (
NB_FREQ
=
/
3, [DEFECT]
/
nbfreq,
[I]
CHAR_CRIT =/
(- 10,10),
[DEFECT]
/
intcc,
),
This key word makes it possible to start calculation, at the end of each increment of time, of a criterion of
stability.
This criterion is useful to detect, during the loading, the point from which one loses
stability (by buckling for example).
This criterion is calculated in the following way: at the end of a pitch of time, in small disturbances,
one solves
(
)
0
det
=
-
G
T
K
K
. K
T
is the coherent tangent matrix at this moment. K
G
is
stamp geometrical rigidity, calculated starting from the stress field at this moment.
In practice, the loading is unstable if
1
<
(in fact
0
1
<
<
-
). The values are calculated
clean by the method of Sorensen (Cd
MODE_ITER_SIMULT
). This can be rather expensive
for the problems of big size.
The key word
CHAR_CRIT
allows to save time by making only one test of Sturm in
provided frequency band. If at least a frequency is found, then one calculates really them
values of the critical loads in this interval.
For great displacements and the great deformations
GREEN (_GR)
or
SIMO_MIEHE
, one
solves
(
)
0
det
=
- Id
K
T
because K
T
K contains then
G
(and possibly K
p
).
The criterion is then a criterion of instability: when
change sign (thus passes by 0) it
loading is unstable.
The key word
NB_FREQ
(3 per defect) the number of critical loads indicates to calculate. In fact
only the first is enough but there can be multiple modes
One stores the clean mode the corresponding to smallest critical load (in absolute value) in
S.D.
RESULT
, under the name
MODE_FLAMB
. This clean mode can be extracted and visualized
(like a field of displacements or a conventional clean mode). It is standardized to 1 on more
large component of displacement.
3.15 Key word
SENSITIVITY
The syntax of this key word common to several controls is described in the document [U4.50.02].
3.16 Word
key
FILING
FILING =
Allows to file or certain results with all or certain moments of calculation.
In the absence of this key word all the pitches of time are filed, including the moments of calculations
lately created by automatic recutting of the pitch of time. Filing allows
to appreciably reduce the size of the bases by selecting the backed up moments.
Note:
In the presence of contact, one cannot file more than 99.999 moments of calculations.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
28/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.16.1 Operand
LIST_INST
/
INST
/
PAS_ARCH
/“LIST_INST”
= list_r8
/
“INST”
=
l_r8
/
“PAS_ARCH”
=
npas
The designation of the moments to be stored is carried out either by a list of moments (
list_r8
or
l_r8
) provided that the evolution is ordered (
EVOLUTION: CHRONOLOGICAL
or
RETROGADE
, cf [§3.6.1]) or then by a frequency of filing (all them
npas
time).
In the absence of these key words all the pitches of time are filed.
Two note:
· the last pitch of calculation is always stored to be able to carry out a recovery,
· if one employs an access by list of moments, then the moments of calculations lately
created by automatic recutting of the pitch of time are not filed
3.16.2 Operand
PRECISION
PRECISION = prec
Cf [U4.71.00]
3.16.3 Operand
ARCH_ETAT_INIT/NUME_INIT/DETR_NUME_SUIV
/`ARCH_ETAT_INIT `
= “NOT”
[DEFECT]
“YES”
Only for one concept not réentrant if not error message. Allows to impose
the filing of the initial state in the sequence number 0 (interesting when the initial state comes
of another
STAT_NON_LINE
. Allows to have the 1
er
not on a curve).
/`NUME_INIT `
= nuinit
Only for one réentrant concept if not error message. Allows to specify from
which sequence number one files.
By defect:
· if the initial state is not fixed by the calculated concept, it is about the last sequence number
+1 (example A),
· if the calculated concept coincides with the concept which fixes L `initial state, it is about the number
of command +1 pennies
ETAT_INIT
(example B and C).
DETR_NUME_SUIV
=
“NOT”
[DEFECT]
“YES”
This operation can result in crushing preexistent sequence numbers: the key word
DETR_NUME_SUIV
confirm this destruction, while its absence puts an end to calculation.
With - Simple example
LIST =
DEFI_LIST_REEL (=0 BEGINNING.,
INTERVAL
=_F (UNTIL
=5.,
NUMBERS
=5)),
U1 = STAT_NON_LINE (INCREMENT =_F (
LIST_INST
=LIST,
INST_FIN
=3.))
,
U2 = STAT_NON_LINE (INCREMENT =_F (LIST_INST =LIST)) ,
U2 = STAT_NON_LINE (reuse=U2,
ETAT_INIT
=_F (EVOL_NOLI
=U1),
INCREMENT
=_F (LIST_INST
=LIST),
FILING
=_F (LIST_INST
=
LIST))
,
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
29/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
The final result for the filing of U2 is as follows:
number of filing
: 1
2
3
4
5
6
7
corresponding moments
: 1. 2. 3. 4. 5. 4. 5.
B - Simple example
LIST = DEFI_LIST_REEL (=0 BEGINNING.,
INTERVAL
=_F (UNTIL
=10.,
NUMBERS
=5)),
U2 = STAT_NON_LINE (
INCREMENT =_F (LIST_INST =LIST)) ,
&U2 = STAT_NON_LINE (
reuse=U2,
ETAT_INIT
=_F (
EVOL_NOLI
=U2,
INST
=4.),
INCREMENT
=_F (
LIST_INST
=LIST),
FILING
=_F (
LIST_INST
=LIST,
DETR_NUME_SUIV
= ' OUI'))
,
The result of filing for 1st U2 is as follows:
number of filing
: 1
2
3
4
5
corresponding moments
: 2.
4.
6.
8.
10.
The final result of filing for U2 is as follows (by defect nuinit = 3):
number of filing
: 1
2
3
4
5
corresponding moments
: 2.
4.
6.
8.
10.
C - Example with NUME_INIT
LIST = DEFI_LIST_REEL (=0 BEGINNING.,
INTERVAL
=_F (UNTIL
=10.,
NUMBERS
=5)),
U2 = STAT_NON_LINE (
INCREMENT =_F (LIST_INST =LIST)) ,
U2 = STAT_NON_LINE (
reuse=U2,
ETAT_INIT
=_F (
EVOL_NOLI
=U2,
INST
=4.),
INCREMENT
=_F (
LIST_INST
=LIST),
FILING
=_F (
LIST_INST
=LIST,
NUME_INIT
=2
,
DETR_NUME_SUIV
= ' OUI'))
,
The result of filing for 1st U2 is as follows:
number of filing
: 1
2
3
4
5
corresponding moments
: 2.
4.
6.
8.
10.
The final result of filing for U2 is as follows:
number of filing
: 1
2
3
4
corresponding moments:
2. 6. 8.
10.
3.16.4 Operand
CHAM_EXCLU
CHAM_EXCLU = |
“DEPL”
|
“SIEF_ELGA”
|
“VARI_ELGA”
|
“VARI_NON_LOCAL”
|
“LANL_ELGA”
Allows to specify the fields which will not be filed, except with the last pitch of time.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
30/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.17 Word
key
DISPLAY
This keyword factor makes it possible to personalize the display of the table of convergence in
STAT_NON_LINE.
DISPLAY:
If this keyword is not indicated, the table is displayed in “STANDARD” mode and with
INFO_RESIDU=' NON'.
Each occurrence of DISPLAY relates to the display of a column and its format. The command of
columns given by the succession of the NOM_COLONNE is respected.
3.17.1 Operand
UNIT
UNIT =
links
The table of convergence will be duplicated in the file of unit links.
Note:
The unit can be repeated with each occurrence of the keyword factor but only first is
taking into account (with display of an alarm).
3.17.2 Operand NOM_COLONNE
NOM_COLONNE
=
|
“STANDARD”,
|
“MINIMUM”,
|
“ITER_NEWT”,
|
“INCR_TPS”,
|
“RESI_RELA”,
|
“RELA_NOEU”,
|
“RESI_MAXI”,
|
“MAXI_NOEU”,
|
“RESI_REFE”,
|
“REFE_NOEU”,
|
“RELI_ITER”,
|
“RELI_COEF”,
|
“PILO_PARA”,
|
“LAGR_ECAR”,
|
“LAGR_INCR”,
|
“LAGR_ITER”,
|
“MATR_ASSE”,
|
“ITER_DEBO”,
|
“CTCD_ITER”,
|
“CTCD_INFO”,
|
“CTCD_GEOM”,
|
“CTCD_NOEU”,
|
“CTCC_CONT”,
|
“CTCC_FROT”,
|
“CTCC_GEOM”,
Type of the column to be displayed (each value corresponds to a displayed column):
ITER_NEWT: number of the iteration of Newton in progress. The column is marked by “X” as long as it
convergence there on all the criteria did not have.
INCR_TPS: moment of current calculation.
RESI_RELA and RELA_NOEU: value of RESI_GLOB_RELA and display of the node where it is maximum.
column is marked by X as long as the residue is larger than that specified by the user
(operand RESI_GLOB_RELA).
RESI_MAXI and MAXI_NOEU: value of RESI_GLOB_MAXI and display of the node where it is maximum.
column is marked by X as long as the residue is larger than that specified by the user
(operand RESI_GLOB_MAXI).
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
31/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
RESI_REFE and REFE_NOEU: value of RESI_REFE_RELA and display of the node where it is maximum.
column is marked by X as long as the residue is larger than that specified by the user
(operand RESI_REFE_RELA).
RELI_ITER and RELI_COEF: iteration count and linear coefficient of search.
PILO_PARA: value of the parameter of piloting.
LAGR_ECAR,
LAGR_INCR and LAGR_ITER: parameters of Lagrangian increased
(see LAGR_NON_LOCAL)
MATR_ASSE: option of assembly for the matrix (elastic, tangent, secant)/
ITER_DEBO: indicate an iteration of Borst for the plane stresses or the behaviors
unidimensional (see COMP_INC)
CTCD_ITER: iteration count intern contact/friction, methods discrete. The column is
marked by X as long as the contact did not converge on the geometry.
CTCD_INFO: information on the state of contact for the discrete methods:
· ALGO: resolution of the problem of contact (iterations intern)
· ALGO/REAC_GEOM: resolution of the problem of contact (internal iterations) and updated of
geometry for reactualization
· INIT_GEOM/ALGO: initialization of the geometry for the contact and resolution of the problem of
contact.
· ATT_PT_FIXE: do not make an attempt fixes for the contact discrete methods
CTCD_GEOM: value of maximum displacement for the geometrical reactualization of the contact,
discrete methods.
CTCD_NOEU: node where the value of displacement is maximum during the geometrical reactualization
contact, discrete methods.
CTCC_GEOM: number of the iteration of contact continuous method at the time of the loop on the geometry.
column is marked by X as long as one did not converge.
CTCC_FROT: number of the iteration of contact continuous method at the time of the loop on the threshold of
friction. The column is marked by X as long as one did not converge.
CTCC_CONT: number of the iteration of contact continuous method at the time of the loop on the state of contact
(active stresses). The column is marked by X as long as one did not converge.
Composite types (displays several columns):
STANDARD: standard display (by defect) of the table of convergence. Contains:
· The number of the iteration of Newton (ITE_NEWT)
· All columns necessary according to functionalities' activated (linear search, contact,
piloting,…)
· The value of residues (RESI_MAXI and RESI_RELA)
MINIMUM: minimum display of the table of convergence. Contains:
· The number of the iteration of Newton (ITER_NEWT)
· The value of residues (RESI_MAXI and RESI_RELA)
Note:
· One cannot ask more than sixteen columns (16 columns of 16 characters, that is to say a width
total of 256)
· The columns are cumulable: one can ask for the MINIMUM display and add one
unspecified column
· One can have several times the same column
· As long as “X” is displayed in column ITER_NEWT, calculation did not converge. This
depends of course on the value of the residues but also of the convergence of the contact or on
De Borst.
· For the method of contact continues, the iterations of Newton constitutes an internal loop
with three other loops (CTCC_GEOM, CTCC_FROT and CTCC_CONT). ITER_NEWT is not thus
not in first position in “STANDARD” mode and it is the marking of columns CTCC_ *
who exploits the part of final Justice of the Peace convergence.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
32/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
3.17.3 Operand INFO_RESIDU
INFO_RESIDU
=
“NOT”,
[DEFECT]
“YES”
This operand makes it possible to add a column for each residue evaluated (RESI_RELA, RESI_MAXI
and RESI_REFE). This column will indicate the node where the residue is maximum, which can help
the user when there are difficulties of convergence. For example, to see whether the material were
badly definite with an incorrect value on an element.
This option is strictly equivalent to the addition of columns RELA_NOEU, RELA_MAXI or
RELA_REFE when one completely describes the display of the table of convergence but allows
to display information on the nodes when one is in STANDARD” or “MINIMUM” mode “,
without needing to describe all the other columns.
3.17.4 Operands LONG_R, PREC_R and LONG_I
LONG_R
=/
12
[DEFECT]
/
long_r
[I]
PREC_R =
/5
[DEFECT]
/
prec_r
[I]
LONG_I =
/6
[DEFECT]
/
long_i
[I]
These operands make it possible to modify the display of information in the table of
convergence. All the columns have a fixed width of 16 characters. When information is
a reality, one can require a personalized display: the length long_r of displayed reality
(maximum 16) and numbers it significant digits.
When it is an entirety, one can regulate the length by long_i. For a it, character string
format is always of 16 characters.
3.18 Operand
OBSERVATION
The syntax of this key word common to the control
DYNA_NON_LINE
is described in the document
[U4.53.01].
3.19 Operand
SOLV_NON_LOCAL
The syntax of this key word is identical to key word SOLVEUR describes in the document [U4.50.01]. With
to use for a nonlocal model.
3.20 Operand
LAGR_NON_LOCAL
The integration of nonlocal laws of behavior imposes the resolution of a total problem (on all
the structure): the minimization of a functional calculus energy (the expression of Lagrangian increased) by
report/ratio with a scalar nodal variable.
The resolution of this problem is carried out by means of an algorithm primal newton and dual BFGS
compound, which consists of two phases:
· Resolution of the primal problem:
- Minimization compared to the variable interns nonlocal and its gradient (
cham_elem
)
- Minimization compared to the variable interns with the nodes (
cham_no
)
- Primal Test of convergence: the largest component of the assembled residue
· Resolution of the dual problem: (Maximization compared to the multipliers of Lagrange)
- Calculation of a direction of descent BFGS
- Linear Search by method of Wolfe
- Dual Test of convergence: the largest component of the gradient
- Reactualization of the multipliers of Lagrange
ITER_PRIM_MAXI: iterprimmax
(10 per defect)
Iteration count maximum for the resolution of the primal problem.
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
33/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
RESI_PRIM_ABSO: resiprimab
Precision for the test of convergence for the primal problem.
ITER_DUAL_MAXI: iterdmax
(50 per defect)
Iteration count maximum for the resolution of the dual problem.
RESI_DUAL_ABSO: residabso
Precision for the test of convergence for the dual problem.
R: rho
(1000 per defect)
Coefficient of penalization of Lagrangian increased.
Note:
As the precision of the dual problem strongly depends on that of the primal problem, one
advise to choose a better precision for the primal problem, for example 100 or
1000 times more than for the dual problem.
3.21 Operand
INFORMATION
INFORMATION
:
inf
Allows to carry out in the file message various intermediate impressions in the presence of
unilateral contact treaty by the method of the active stresses.
inf =
1
impression of the list of the nodes in contact after convergence with each
iteration of Newton.
= 2
idem
1
more impression of associations/dissociations of nodes enters
iterations of the method of the active stresses.
Other impressions are made systematically during nonlinear calculation, independently
value assigned to the key word
INFORMATION
: they are the impressions of the residues and the increments
relative of displacement during iterations of Newton.
3.22 Operand
TITRATE
TITRATE: tx
tx
is the title of calculation. It will be printed at the head results. See [U4.03.01].
Code_Aster
®
Version
8.2
Titrate:
Operator
STAT_NON_LINE
Date
:
31/01/06
Author (S):
Mr. ABBAS
Key
:
U4.51.03-H1
Page:
34/34
Instruction manual
U4.5- booklet: Method of resolution
HT-62/06/004/A
Intentionally white left page.