background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
1/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
Organization (S):
EDF-R & D/AMA














Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
Document: D4.06.14



Structures of data related to contact-friction




Summary:

This document describes the structures of data necessary to definition (SD '
DEFI_CONT
') and with the resolution
(SD '
RESO_CONT
') of the problems of contact-friction defined by the key word
CONTACT
of the operator
AFFE_CHAR_MECA
. One endeavors to give the detailed instructions of the majority of the tables used in
corresponding routines. Description is purely data-processing, and it is advised to read them before
reference materials [R5.03.50], [R5.03.51], of use [U4.25.01], implementation practical
[U2.04.04] and of maintenance of the contact [D9.05.02].
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
2/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
1 Philosophy
general
Potential areas of contact are defined in the operator
AFFE_CHAR_MECA
for each
occurrence of the key word
“CONTACT”
. Each area is defined by an occurrence of the key word. One
area of contact comprises several surfaces (two in general) which one seeks to prevent
the interpenetration two to two. In the case of methods
“NODAL”
,
“MAIT_ESCL”
and
“MAIT_ESCL_SYME”
, there are two surfaces whose composition is given under the key words
GROUP_MA_MAIT
(or
MAILLE_MAIT
) and
GROUP_MA_ESCL
(or
MAILLE_ESCL
). For the methods
“TERRITORY”
and
“HIERARCHICAL”
(not implemented to date), the key word will be used
GROUP_MA
(or
NET
): in this case, each group of meshs (each mesh) of the list will define
a potential surface of contact (there could thus be more than two surfaces per area).
data relating to the various areas and surfaces of contact are stored in a structure of
data of the type
“DEFI_CONT”
whose name is
TANK (1:8)//“.CONTACT”
.
In the operators
STAT_NON_LINE
and
DYNA_NON_LINE
, it is supposed that only one load contains
contact (one checks it in the routine
nmdome
). During pitch of time and iterations of
Newton, one fills of the tables of size fixes (dimensioned using the maximum of nodes slaves)
who contain the data necessary to the processing of the contact (structure of data of the type
“RESO_CONT”
). They are sometimes under-tables resulting from the tables created in
AFFE_CHAR_MECA
:
they relate to only the nodes or meshs in the course of processing with the algorithm (areas of contact
effective current). In these tables, information is followed sequentially without concept of
area or of surface of contact: very coarsely, one stores the couples node slave - mesh (or
node) main and characteristics associated (ddls concerned, coefficients, components of
normal, play running).
Note:
·
The system of contact is composed of several areas, themselves consistent in
surfaces, made up of meshs, container of the nodes,
·
Surfaces of contact are identified by their absolute number
I
in the list of all
surfaces of contact, all areas confused,
·
Only tables
CONTMA
,
CONTNO
and
SANSNO
index the nodes and the meshs by
their absolute number in the code; all the other tables use the index in
CONTMA
and
CONTNO
to indicate a mesh or a node,
·
One always takes three components for the normal, in 2D as in 3D. On the other hand,
the table of the degrees of freedom contains of them two per node in 2D, and three in 3D.
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
3/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
2
Structure of data
DEFI_CONT

The structure of data
DEFI_CONT
contains the tables defining the potential areas of contact (created
in
AFFE_CHAR_MECA
, except
DDLCO
and
PDDLCO
created in
STAT_NON_LINE
).
2.1
List variables

It is supposed that only one load contains the key word factor
“CONTACT”
.
All the tables start with
TANK (1:8)//“.CONTACT”
(
variable
DEFICO (1:16)
), with the suffix:
Variable suffix
Type
Dimension Subscripted
by
Contents
.NDIMCO NDIMCO I 9+NZOCO
list useful entireties, and numbers effective
nodes slaves for each area
(the max at the beginning)
.METHCO METHCO I 1+8 * NZOCO
number of
area of contact
a number of areas and characteristic of
method of pairing for each one
.TOLECO TOLECO R 2 * NZOCO
number of
area of contact
parameters of geometrical tolerance
for pairing
.CONVCO CONVCO I 3 * NZOCO
number of
area of contact
parameters of convergence
.SYMECO SYMECO I NZOCO+1
information on the symmetrical areas
of pairing (
MAIT_ESCL_SYME
)
.PZONECO PZONE
I
NZOCO+1
number of
area of contact
number of the last surface of
each area
.MAILCO CONTMA I
NMACO
pointer
PSURMA
list numbers of the meshs of
contact of various surfaces
potential
.PSUMACO PSURMA
I
1+NSUCO
number of
surface
index of the last mesh of each
surface in the vector
CONTMA
.NOEUCO CONTNO I
NO
pointer
PSURNO
list numbers of the nodes of
contact of various surfaces
potential
.PSUNOCO PSURNO
I
1+NSUCO
number of
surface
index of the last node of each
surface in the vector
CONTNO
.NOEUQU CONOQU I 3 * NO/2
pointer
PNOQUA
list numbers of the nodes of
contact “quadratic” of different
potential surfaces
.PNOEUQU PNOQUA
I
1+NSUCO
number of
surface
index of the last node of each
surface in the vector
CONOQU
.MANOCO MANOCO I
NMANO
pointer
PMANO
indices of the meshs of
CONTMA
containing a node given of
CONTNO
.PMANOCO PMANO
I
1+NNOCO
index of node
in
CONTNO
index in
MANOCO
of the last
net containing a given node
.NOMACO NOMACO I
NO
pointer
PNOMA
indices of the nodes of
CONTNO
belonging to a mesh given of
CONTMA
.PNOMACO PNOMA
I
1+NMACO
index of mesh
in
CONTMA
index in
NOMACO
last node
belonging to a given mesh
.MAMACO MAMACO I
NMAMA
pointer
PMAMA
indices of the meshs of
CONTMA
adjacent with a given mesh
.PMAMACO PMAMA
I
1+NMACO
index of mesh
in
CONTMA
index in
MAMACO
of the last
adjacent mesh with a given mesh
.SSNOCO SANSNO I
NO
pointer
PSANS
absolute numbers of the nodes to be excluded
nodes slaves
.PSSNOCO PSANS
I
1+NZOCO
number of
area of contact
index of the last node to be excluded in
SANSNO
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
4/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
.JSUPCO JEUSUP R8
NZOCO
number of
area of contact
value of the fictitious play
.NOZOCO NOZOCO I
NO
index of node
in
CONTNO
number of the area of contact to which
the node belongs
.CHAMCO CHAMCO I
NZOCO
number of
area of contact
code field to which apply
unilateral conditions
.COEFCO COEFCO R8
NZOCO
number of
area of contact
coefficient of the unilateral relation for
pressure or the temperature
.DDLCO DDLCO I
NDDL
pointer
PDDL
numbers of the degrees of freedom
potentially implied in the writing
unilateral relations
.PDDLCO PDDL
I
NO
index of node
in
CONTNO
index in
DDLCO
the last ddl of one
node of
CONTNO
given
.JEUFO1 JJFO1 K8
NZOCO
Number of
surface
contact
Fictitious play when it is given by one
function of space in
AFFE_CHAR_MECA_F
.JEUFO2 JJFO2 K8
NZOCO
Number of
surface
contact
Fictitious play when it is given by one
function of space in
AFFE_CHAR_MECA_F
DIRCO JDIR R8
3 * NZOCO
Number of
surface
contact
Direction fixes nodal pairing
data by
VECT_Y
RUB IFRO R8 NESMAX
Number of
node slave
Coefficient of friction of Coulomb
PENAL IPENA
R8
2 * NESMAX
Number of
node slave
Coefficient of regularization of the contact
and of friction
COMAFO ICOMA R8 NESMAX
Number of
node slave
Value of
COEF_MATR_FROT
TANDEF JTGDEF
R8 3 * NZOCO
Number of
surface
contact
Value of
VECT_Y
NORLIS JNORLI I NZOCO+1
Number of
surface
contact
Indicate the presence of smoothing of
normals


This part gathers the objects suitable for the method
CONTINUOUS
:
.CARACF JCMCF
R 6 * NZOCO+1
number of
area of contact
integration and coefficients of
regularization
.ECPDON JECPD
I 5 * NZOCO+1
number of
area of contact
parameters of the loops of
method
CONTINUOUS
.MAESCL JMAESC I 3 * NTMA+1
number of
net slave
for each mesh one gives the number
of its area numbers points of
contact
.NOESCL JNOESC R 10 * NO
number of node
of contact
vectors tangent and normal of
each point.
.TABFIN JTABF
R 16 * NTPC+1
number of point
of contact
characteristics of pairing.
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
5/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
2.2 Table
NDIMCO
(addresses
JDIM
)
ZI (JDIM) =
NDIM
dimension of space (two or three)
ZI (JDIM+1) =
NZOCO
a number of areas of contact
ZI (JDIM+2) =
NSUCO
a number of surfaces of contact
ZI (JDIM+3) =
NMACO
a number of meshs of contact
ZI (JDIM+4) =
NO
a number of nodes of contact
ZI (JDIM+5) =
NMANO
dimension of
MANOCO
ZI (JDIM+6) =
NO
dimension of
NOMACO
ZI (JDIM+7) =
NMAMA
dimension of
MAMACO
ZI (JDIM+8)
= NESMAX
a maximum number of nodes slaves
ZI (JDIM+8+IOC) =
numbers effective nodes slaves in the area
IOC
(a number
maximum at the time of initialization),
IOC=1, NZOCO

2.3 Table
METHCO
(addresses
JMETH
)
ZI (JMETH) = NZOCO
: a number of areas of contact
For the area
N
:
ZI (JMETH+8 * (n-1) +1) =
- 1 if
APPARIEMENT= “NOT”
0 if
“NODAL” APPARIEMENT=
1 if
APPARIEMENT= “MAIT_ESCL”
or
“MAIT_ESCL_SYME”
2 if
APPARIEMENT= “TERRITORY”
3 if
“HIERARCHICAL” APPARIEMENT=
4 if
VECT_NORM_2
is defined
ZI (JMETH+8 * (n-1) +2) =
1
VECT_Y
is informed and 0 if not
ZI (JMETH+8 * (n-1) +3) =
not used
ZI (JMETH+8 * (n-1) +4) =
1 if linear projection (rectilinear segment or plane triangle)
2 if quadratic projection
ZI (JMETH+8 * (n-1) +5) =
+1
if
RECHERCHE= “NOEUD_BOUCLE”
+
/-
2
if
RECHERCHE=
“NOEUD_VOISIN”/“MAILLE_VOISIN”
+/
-
3
if
RECHERCHE= “NOEUD_BOITE”/“MAILLE_BOITE”
ZI (JMETH+8 * (n-1) +6) =
- 1 if method of
CONTACT
used is `
PENALISA'
0 if method of
CONTACT
used is `
CONTRAIN'
1 if method of
CONTACT
used is `
LAGRANGI'
2 if method of
FRICTION
2D used is `
LAGRANGI'
3 if method of
FRICTION
2D or 3D used is `
PENALISA'
4 if method of
FRICTION
3D used is `
LAGRANGI'
5 if method of
CONTACT
and of
FRICTION
used is
`
PENALISA'
6 if the method used is `
CONTINUE'
ZI (JMETH+8 * (n-1) +7) =
0 if
REAC_GEOM= “WITHOUT”
- 1 if
“AUTOMATIC” REAC_GEOM=
if not
NB_REAC_GEOM
(geometrical frequency of reactualization)
ZI (JMETH+8 * (n-1) +8) =
0 if
NORMALE= “MAIT”
1 if
NORMALE= “MAIT_ESC”
ZI (JMETH+9 * (n-1) +9) =
0 if
STOP_SINGULIER= “YES”
1 if
STOP_SINGULIER= “NOT”
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
6/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
2.4 Table
CONVCO
(addresses
JCONV
)
For the area
N
:
ZI (JCONV+3 * (n-1)) =
0 if
STOP_SINGULIER= “YES”
1 if
STOP_SINGULIER= “NOT”
ZI (JCONV+3 * (n-1) +1) =
NB_RESOL
ZI (JCONV+3 * (n-1) +2) =
ITER_MULT_MAXI
2.5 Table
TOLECO
(addresses
JTOLE
)
For the area
N
:
ZI (JTOLE+2 * (n-1)) =
TOLE_PROJ_EXT
ZI (JTOLE+2 * (n-1) +1) =
TOLE_PROJ_INT

2.6 Table
SYMECO
(addresses
JSYME
)
ZI (JSYME)
: A number of symmetrical areas of contact
For the symmetrical area
N
:
ZI (JSYME+n) =
Number of the area main partner at symmetrical area N

2.7 Table
CARACF
(addresses
JCMCF
)
ZI (JCMCF) = NZOCO
: numbers total areas of contact.
In this table some parameters for the methods are stored
CONTINUOUS
,
LAGRANGIAN
and
PENALIZATION
. For the method
CONTINUOUS
, one specifies, inter alia, the diagram of integration with
to use for the terms of contact and friction and the coefficients of increase. Let us recall that
integration with the nodes is taken by defect
INTEGRATION=' NOEUD'
and that diagrams
SIMPSON
,
SIMPSON1
and
SIMPSON2
are available only in 2D.
For the area
N
:
CARACF (1+6 * (n-1) +1) =
1
if
INTEGRATION=' NOEUD'
2
if
INTEGRATION=' GAUSS'
3
if
INTEGRATION=' SIMPSON'
4
if
INTEGRATION=' SIMPSON1'
5
if
INTEGRATION=' SIMPSON2'
CARACF (1+6 * (n-1) +2) =
Coefficient of increase for the contact
COEF_REGU_CONT
CARACF (1+6 * (n-1) +3) =
Coefficient of increase for friction
COEF_REGU_FROT
CARACF (1+6 * (n-1) +4) =
Coefficient of Coulomb for friction.
CARACF (1+6 * (n-1) +5) =
1
if
FROTTEMENT=' SANS'
3
if
FROTTEMENT=' COULOMB'
CARACF (1+6 * (n-1) +6) =
the value of
COEF_MATR_FROT

2.8
System of contact: pointer
PZONE
(addresses
JZONE
)
absolute number (
I
) of the first surface of the area
N
:
ZI (JZONE+n-1) +1
absolute number (
I
) of the last surface of the area
N
:
ZI (JZONE+n)
a number of surfaces of the area
N
:
ZI (JZONE+n) - ZI (JZONE+n-1)
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
7/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
2.9 Area
N
: pointers
PSURMA
,
PSURNO
and
PNOQUA
(addresses
JSUMA
,
JSUNO
and
JNOQUA
)
number of the first
surface
:
i1 = ZI (JZONE+n-1) +1
number of the last surface:
i2 = ZI (JZONE+n)
a number of meshs of the area:
ZI (JSUMA+i2) - ZI (JSUMA+i1-1)
= ZI (JSUMA+ZI (JZONE+n)) - ZI (JSUMA+ZI (JZONE+n-1))
a number of nodes of the area:
ZI (JSUNO+i2) - ZI (JSUNO+i1-1)
= ZI (JSUNO+ZI (JZONE+n)) - ZI (JSUNO+ZI (JZONE+n-1))
a number of nodes of the area
being node mediums of one
quadratic element
:
ZI (JNOQUA+i2) - ZI (JNOQUA +I1-1)
= ZI (JNOQUA +ZI (JZONE+n)) - ZI (JNOQUA +ZI (JZONE+n-1))

2.10 Surface
I
: tables
CONTMA
,
CONTNO
and
CONOQU
(addresses
JMACO
,
JNOCO and JNOQU
), pointers
PSURMA
,
PSURNO
and
PNOQUA
(addresses
JSUMA
,
JSUNO
and
JNOQUA
)
The number of meshs of surface
I
is:
nbma = ZI (JSUMA+i) - ZI (JSUMA+i-1)
The index in
CONTMA
first mesh of surface
I
is:
ZI (JSUMA+i-1) +1
The index in
CONTMA
of the last mesh of surface
I
is:
ZI (JSUMA+i)
The list of the numbers of the meshs of surface
I
is
ZI (JMACO+jdecma+ima-1)
for
ima=1, nbma
,
with
jdecma = ZI (JSUMA+i-1)
.
ima
ième
net surface
I
has as an absolute number:
ZI (JMACO+jdecma+ima-1)
; its index
in
CONTMA
is:
posma = jdecma+ima.
Net index
posma
in
CONTMA
: its absolute number is
ZI (JMACO+posma-1).
The number of nodes of surface
I
is:
nbno = ZI (JSUNO+i) - ZI (JSUNO+i-1)
The index in
CONTNO
first node of surface
I
is:
ZI (JSUNO+i-1) +1
The index in
CONTNO
last node of surface
I
is:
ZI (JSUNO+i)
The list of the numbers of the nodes of surface
I
is
ZI (JNOCO+jdecno+ino-1)
for
ino=1, nbno
,
with
jdecno = ZI (JSUNO+i-1)
.
ino
ième
node of surface
I
has as an absolute number:
ZI (JNOCO+jdecno+ino-1)
; its index
in
CONTNO
is:
posno = jdecno+ino.
Node of index
posno
in
CONTNO
: its absolute number is
ZI (JNOCO+posno-1).
The number of nodes of surface
I
being node medium of a quadratic mesh is:
No =
ZI (JNOQUA+i) - ZI (JNOQUA +i-1)
The index in
CONOQU
first “quadratic” node medium of surface
I
is: 3 *
ZI (JNOQUA
+i-1) +1
The index in
CONOQU
last node “quadratic” medium of surface
I
is: 3 *
ZI (JNOQUA
+i) - 2
The list of the numbers of the nodes medium “
quadratic
” of surface
I
is
ZI (JNOQU
+jdecqu+3 * (inq-1))
for
inq=1, No
, with
jdecqu = 3 * ZI (JNOQUA +i-1) +1
.
The list of the numbers of the associated nodes node for surface
I
is
ZI (JNOQU
+jdecqu+3 * (inq-1) +1)
and
ZI (JNOQU +JDECQU+3 * (INQ-1) +2)
for
inq=1, No
, with
jdecqu = 3 * ZI (JNOQUA +i-1) +1
.
inq
ième
“quadratic” node medium of surface
I
has as an absolute number:
ZI (JNOQU
+jdecqu+3 * (inq-1) - 1)
; its index in
CONOQU
is:
posqu = jdecqu+3 * (inq-1).
“Quadratic” node of index
posqu
in
CONOQU
: its absolute number is
ZI (JNOQU+posqu-1).
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
8/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
2.11 Tables
MANOCO
,
NOMACO
and
MAMACO
(addresses
JMANO
,
JNOMA
and
JMAMA
) pointer
PMANO
,
PNOMA
and
PMAMA
(addresses
JPOMA
,
JPONO
and
JPOIN
)
For the node of index
posno
in
CONTNO
, the index in
CONTMA
meshs of contact of same
surface containing the node is
ZI (JMANO+jdec+ima-1)
with
jdec = ZI (JPOMA+posno-1)
and
ima
varying
1
with
nbma
, for
nbma = ZI (JPOMA+posno) - ZI (JPOMA+posno-1)
, or even
ZI (JMANO+k-1)
for
K
varying
ZI (JPOMA+posno-1) +1
with
ZI (JPOMA+posno)
.
For the mesh of index
posma
in
CONTMA
, the index in
CONTNO
nodes of contact of this
mesh is
ZI (JNOMA+jdec+ino-1)
with
jdec = ZI (JPONO+posma-1)
and
ino
varying
1
with
nbno
, for
nbno = ZI (JPONO+posma) - ZI (JPONO+posma-1)
, or even
ZI (JNOMA+k-1)
for
K
varying
ZI (JPONO+posma-1) +1
with
ZI (JPONO+posma)
.
The list of the indices in
CONTMA
meshs close to the mesh of index
posma
is:
ZI (JMAMA+jdec+ima-1)
, with
jdec = ZI (JPOIN+posma-1)
and
ima
varying
1
with
nbma
,
for
nbma = ZI (JPOIN+posma) - ZI (JPOIN+posma-1)
.

2.12 Table
SANSNO
and pointer
PSANS
(addresses
JSANS
and
JPSANS
)
For the area
N
:
a number of nodes to be excluded from the nodes slaves:
nsans = ZI (JPSANS+n) - ZI (JPSANS+n-1)
absolute numbers of the nodes to be excluded:
ZI (JSANS+jdec+ino-1)
, for
ino = 1
,
nsans
, with
jdec = ZI (JPSANS+n-1)
.

2.13 Table
JEUSUP
(addresses
JJSUP
)
For the area
N
:
ZR (JJSUP+n-1)
= value for the area of the fictitious play
(DIST_1+DIST_2
, or
COEF_IMPO)
given by the user.

2.14 Tables
NOZOCO
(addresses
JZOCO
),
CHAMCO
(addresses
JCHAM
) and
COEFCO
(addresses
JCOEF
)
For the node of index
posno
in
CONTNO
, the number of the area of contact to which the aforementioned
belongs is:
N = ZI (JZOCO+posno-1)
.
For the area
N
:
Code field to which applies the unilateral relation:
icode = ZI (JCHAM+n-1)
icode = +1
: relation on displacements (with pairing: “deformable” contact)
icode = - 1
: relation on displacements (without pairing: “rigid” contact)
icode = - 2
: relation on the pressure (without pairing: only for one modeling
“THM”)
icode = - 3
: relation on the temperature (without pairing “: only for one modeling
“THM”)
icode = - 4
: relation on pressure 1 (without pairing “: only for one modeling
“THM”)
icode = - 5
: relation on pressure 2 (without pairing “: only for one modeling
“THM”)
Multiplying coefficient of the unilateral relation:
ZR (JCOEF+n-1)
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
9/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
2.15 Table
DDLCO
and pointer
PDDLCO
(addresses
JDDL
and
JPDDL
)
These two tables are in the structure of data
DEFI_CONT
but are defined in the routine
crsdco
called by
STAT_NON_LINE
.
For the node of index
posno
in
CONTNO
, the degrees of freedom are
ZI (JDDL+jdecdl+iddl-1)
for
iddl
varying
1
with
nbddl with
nbddl = ZI (JPDDL+posno) - ZI (JPDDL+posno-1)
and
jdecdl = ZI (JPDDL+posno-1)
.
One gathered in the continuation the structures of data suitable for the method
CONTINUOUS
.

2.16 Table
ECPDON
(addresses
JECPD
)
This table gives some parameters total necessary for the algorithmic one used by
method
CONTINUOUS
. Let us note that the parameters
ITER_CONT_MAX
,
ITER_FROT_MAX
and
ITER_GEOM_MAX
fix the maximum number of the loops of contact, threshold of friction and of
geometry. These numbers can not be reached if the test of stop of each ball is satisfied. It
test, to see routine
mmmcri.f
, relates to the value of the relative increment of displacement:
U
U
with
10
=
- 2
.
NR being the number of the area of contact.
ECPDON (1+5 * (N-1) +1) =
1
if
MODL_AXIS=' OUI'
0
if
MODL_AXIS=' NON'
ECPDON (1+5 * (N-1) +2) =
The value of
ITER_CONT_MAX
ECPDON (1+5 * (N-1) +3) =
The value of
ITER_FROT_MAX
ECPDON (1+5 * (N-1) +4) =
The value of
ITER_GEOM_MAX
ECPDON (1+5 * (N-1) +5) =
The initial threshold value
SEUIL_INIT
MAESCL (1):
numbers total areas of contact.

2.17 Table
MAESCL
(addresses
JMAESC
)
NR being the number of the mesh slave.
MAESCL (1+3 * (N-1) +1) =
index of the mesh NR in the table
CONTMA
MAESCL (1+3 * (N-1) +2) =
number of the area of contact of NR
MAESCL (1+3 * (N-1) +3) =
numbers points of contact in NR
MAESCL (1):
numbers total meshs slave.
2.18 Table
NOESCL
(addresses
JNOESC
)
NR is the number of node of contact. The whole parameter I varies from 1 to 3.
NOESCL (1+3 * (N-1) +1) =
0 if NR is slave 1 if not.
NOESCL (1+3 * (N-1) +1+I) =
components of the first tangent vector
NOESCL (1+3 * (N-1) +4+I) =
components of the first tangent vector
NOESCL (1+3 * (N-1) +7+I) =
components of the first tangent vector
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
10/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
NOESCL (1):
numbers total nodes of contact. Let us recall that in 2D the 3rd component of
vectors tangent or normal is taken equalizes to zero.
This table is used for smoothing (routine `
lissag
') which makes it possible to smooth the normals on the surfaces
of contact intervening in the calculation of the matrix of contact. Let us recall that smoothing is made in
two stage. The first consists in carrying out an average of the normals to the meshs which contain
the node of contact. The second consists in interpolating, with the functions of form associated with
the element, a field of normals in any point of the element.

2.19 Table
TABFIN
(addresses
JTABF
)
In this table is classified all information necessary for the resolution concretes
problem of the rubbing contact. This information is described in the routine
mappar
. Let us recall that
pairing is made in an exact way using a method of Newton for the resolution of one
problem of optimization with stresses (cf routine
mprojp
) and which allows us, at the same time
to recover the values of the tangent vectors.
NR being the number of the point of contact.
TABFIN (1):
numbers total points of contact.
whole parameter I varies from 1 to 3.
TABFIN (1+16 * (N-1) +1) =
absolute number of the mesh Master
TABFIN (1+16 * (N-1) +2) =
absolute number of the mesh slave
TABFIN (1+16 * (N-1) +3) =
first barycentric parameter of the point NR
TABFIN (1+16 * (N-1) +4) =
first barycentric parameter of the point in screw-with
live
TABFIN (1+16 * (N-1) +5) =
second barycentric parameter of the point in opposite
TABFIN (1+16 * (N-1) +5+I) =
3 components of the 1st tangent vector
TABFIN (1+16 * (N-1) +8+I) =
3 components of the 2nd tangent vector
TABFIN (1+16 * (N-1) +12) =
second barycentric parameter of the point NR
TABFIN (1+16 * (N-1) +13) =
statute of contact
TABFIN (1+16 * (N-1) +14) =
initial value of contact pressure
TABFIN (1+16 * (N-1) +15) =
number of the area of contact of the point NR
TABFIN (1+16 * (N-1) +16) =
weight of the point of contact
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
11/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
3
Structure of data
RESO_CONT

The structure of data
RESO_CONT
contains the tables defining the effective couples of contact (created
in
STAT_NON_LINE
or
DYNA_NON_LINE
) and variables used in the method of resolution by
methods
STRESSES
,
LAGRANGIAN
and
PENALIZATION
.
3.1
List variables

All the tables start with
RESOCO (1:14) (“&&RESOCO”)
, followed suffix:
Variable suffix Standard Dimension
Subscripted
by
Contents
.APREAC APREAC
I
4 * NZOCO
number of
area of contact
indicator of reactualization for
pairing, fixed counter of pairing,
type of projection and reactualization of
normals, number of surface slave
.APPARI APPARI
I 1+3 * NESMAX
number of
node slave
a number of nodes slaves, and for each one
: index of the node slave, index of the mesh
Master paired, indicating of reactualization
.APMEMO APMEMO
I
4 * NO
index of the node
in
CONTNO
data relating to pairing of the last
blow where this node was slave
.APPOIN APPOIN
I
1+NESMAX
number of
node slave
pointer of navigation in
APCOEF,
APCOFR
and
APDDL
.APCOEF APCOEF R8 30 * NESMAX
pointer
APPOIN
multiplying coefficients of the ddls (1 by
ddl) for imposition of nonthe penetration
(+1 for the node slave, and the opposite of
value of the function of form for each
main node)
.APCOFR APCOFR R8 60 * NESMAX
pointer
APPOIN
multiplying coefficients of the ddls (1 by
ddl) for imposition friction (+1 for
node slave, and opposite of the value of
function of form for each main node)
.APDDL APDDL
I 30 * NESMAX
pointer
APPOIN
numbers of the degrees of freedom of the node
slave and of the nodes of the mesh Master
paired
.NORINI NORINI R8
3 * NO
index of the node
in
CONTNO
direction of evaluation of the normal play on
the whole of the potential nodes of contact
.NORMCO NORMCO R8
3 * NESMAX
number of
node slave
direction of evaluation of the normal play on
connection of contact
.TANGCO TANGCO R8
6 * NESMAX
number of
node slave
direction of evaluation of the tangent play on
connection of contact
.APJEU APJEU R8
NESMAX
number of
node slave
value of the normal play running between the node
slave and the mesh Master paired
.APJEFX APJEFX R8
NESMAX
number of
node slave
value of the tangent play in direction 1
between the node slave and the mesh Master
paired
.APJEFY APJEFY R8
NESMAX
number of
node slave
value of the tangent play in direction 2
between the node slave and the mesh Master
paired
.JEUINI JEUINI R8
NESMAX
number of
node slave
value of the initial play when
REAC_GEOM=' SANS'
.COCO COCO
I
8
to remember the state of preceding contact
.LIAC LIAC
I
3 * NESMAX+1
number of
node slave
absolute numbers of the active connections of
contact-friction
.CONVEC CONVEC K8 3 * NESMAX+1
number of
node slave
Type of active connections
: contact or
adherent friction
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
12/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
.LIOT LIOT
I
4 * NESMAX+4
number of
node slave
removed connections of contact-friction
.MU DRIVEN R8
6 * NESMAX
number of
node slave
multipliers of Lagrange related to
contact-friction
.COEFMU COEFMU R8
NESMAX
number of
node slave
coefficient by which it is necessary to multiply it
multiplier of Lagrange of the front contact
to test its sign
.ATMU ATMU R8
NEQ
number of ddl
nodal forces of contact
.AFMU AFMU R8
NEQ
number of ddl
nodal forces of friction
.DEL0 DELT0 R8
NEQ
number of ddl
vector used in the algorithm of resolution
.DELT R8 DELTA
NEQ
number of ddl
vector used in the algorithm of resolution
.CM1A CM1A
second members used in the algorithm
of resolution of the contact
.CM2 A CM2 A
second members used in the algorithm
of resolution of friction
.CM3 A CM3 A
second members used in the algorithm
of resolution of friction
.MATR MATR
SD of the type
MATR_ASSE
[D4.06.10]: stamp
used in the algorithm of resolution
.SLCS STOC
SD of the type
STOC_LCIEL
[D4.06.07]:
description of the storage of the matrix
MATR
For the link of the variables described with the resolution of the problem of contact by the method of the stresses
active, one will refer to the document [R5.03.50] and for the resolution of the problem of contact-friction to
document [R5.03.51].
3.2 Table
APREAC
(addresses
JREAC
)
For the area
N
:
ZI (JREAC+4 * (n-1))
: reactualization of pairing
0:
not
- 1:
no pairing but initial passage in
rechno
1:
by double loop on the nodes ('
NOEUD_BOUCLE
')
+/­ 2: by vicinity of “last” ('
NOEUD_VOISIN
“/”
MAILLE_VOISIN
')
+/­ 3: by boxes of position ('
NOEUD_BOITE
“/”
MAILLE_BOITE
')
ZI (JREAC+4 * (n-1) +1)
: a number of times where pairing was kept fixed
ZI (JREAC+4 * (n-1) +2)
: type of projection and normal reactualization geometry/
+/­ 1: linear projection
+/­ 2: quadratic projection
> 0
so normal and recomputed co-ordinates
< 0
if not
ZI (JREAC+4 * (n-1) +3)
: not used
3.3 Table
APPARI
(addresses
JAPPAR
)
ZI (JAPPAR) = NESCL
: numbers effective nodes slaves
For
iescl
ième
node slave
ZI (JAPPAR+3 * (iescl-1) +1)
: index in
CONTNO
node slave
ZI (JAPPAR+3 * (iescl-1) +2)
: index in
CONTMA
mesh Master paired
(negative if nodal pairing: opposed index in
CONTNO
paired main node)
(0 if not of pairing)
ZI (JAPPAR+3 * (iescl-1) +3)
: indicator of reactualization
0
: no the reactualization of projection
+1
: reactualized linear projection + normal
+2
: reactualized quadratic projection + normal
­ 1
: reactualized linear projection, not normals (not used)
­ 2
: reactualized quadratic projection, not normals (not used)
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
13/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
3.4 Table
APMEMO
(addresses
JAPMEM
)
For
posno
ième
node of
CONTNO
(slave or not)
ZI (JAPMEM+4 * (posno-1))
:
1 if the node is slave
0 if the node is a Master
- 1 if the node is excluded (belongs to
SNAS_GROUP_NO
)
- 2 if the node is excluded (null pivot during symmetrical pairing)
ZI (JAPMEM+4 * (posno-1) +1)
:
index in
CONTNO
main node nearest the last time that
this node was slave
ZI (JAPMEM+4 * (posno-1) +2)
:
index in
CONTMA
mesh Master paired the last time that it
node was slave
ZI (JAPMEM+4 * (posno-1) +3)
:
number of the box of current position

3.5 Pointer
APPOIN
and tables
APCOEF
,
APCOFR
and
APDDL
(addresses
JAPPTR
,
JAPCOE
,
JAPCOF
and
JAPDDL
)
Tables
APCOEF
and
APDDL
have same dimension (a coefficient by ddl implied). They are subscripted by
even pointer
APPOIN
.

For
iescl
ième
node slave:
ZI (JAPPTR+iescl-1) + 1
: beginning of the arrangement in
APCOEF
and
APDDL
ZI (JAPPTR+iescl)
: end of the arrangement in
APCOEF
and
APDDL
ZI (JAPPTR+iescl) - ZI (JAPPTR+iescl-1) = nbddl1 +
summon
nbddl2
main nodes
nbddl1 =
numbers ddls node slave of index
posno1
in
CONTNO
:
NBDDL1 = ZI (JPDDL+posno1) - ZI (JPDDL+posno1-1),
with posno1 = ZI (JAPPAR+3 * (iescl-1) +1)
nbddl2 =
ddls of each main node of index numbers
posno2
in
CONTNO
:
NBDDL2 = ZI (JPDDL+posno2) - ZI (JPDDL+posno2-1)
jdec1 = ZI (JAPPTR+iescl-1)
ZI (JAPDDL+jdec1+k-1),
k=1, nbddl1
:
number of K
ième
ddl of the node slave
ZR (JAPCOE+jdec1+k-1),
k=1, nbddl1
:
coefficient associated with K
ième
ddl of the node slave
jdec2 = ZI (JAPPTR+iescl-1) + nbddl1
for
m = 1, nmaitr
(
nmaitr
main nodes with each one
nbddl2
ddls) and
K = 1, nbddl2
ZI (JAPDDL+jdec2+ (M-1) * nbddl2+k-1)
: number of K
ième
ddl of the m
ième
main node
ZR (JAPCOE+jdec2+ (M-1) * nbddl2+k-1)
: coefficient associated with K
ième
ddl of the m
ième
main node


Tables
APCOFR
and
APDDL
, subscripted by the same pointer
APPOIN
, are used in the case of presence of
friction. The arrangement is exactly the same one as in what precedes with the details close following:
·
APCOFR
connect the ddl nodes Master and slave concerning displacements in
tangent plan on the surface of contact
·
APCOFR
contains
60 * NESMAX
terms is twice as much as
APCOEF
. Indeed, of 1 with
30 * NESMAX
the relations in a direction of the tangent plan are stored, of
30 * NESMAX+1
with
60 * NESMAX
are
stored relations in the orthogonal direction with the preceding one in the tangent plan (useful
only in 3D).

NB:
In the case of the contact between two surfaces, the coefficients of the ddls are then multiplied by
components of the entering normal of the mesh Master paired. In the case of the rigid contact, them
coefficients are then multiplied by the components of the normal outgoing slave.
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
14/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
3.6 Table
NORINI
(addresses
JNRINI
)
Components of the current normal with
posno
I
ème
node:
ZI (JNRINI+3 * (posno-1) +k-1)
,
k=1,3
3.7 Table
NORMCO
(addresses
JNORMO
)
For
iescl
ième
node slave
ZI (JNORMO+3 * (iescl-1) +k-1)
,
k=1,3
: components of the direction (normalized) of evaluation of the play
normal.
3.8 Table
TANGCO
(addresses
JTANGO
)
For
iescl
ième
node slave
ZI (JTANGO+6 * (iescl-1) +k-1)
,
k=1,3
: components of the first direction (normalized)
of evaluation of the tangent play.
ZI (JTANGO+6 * (iescl-1) +k-1)
,
k=4,6
: components of the second direction (normalized)
of evaluation of the tangent play.
3.9 Table
APJEU
(addresses
JAPJEU
)
For
iescl
ième
node slave
ZI (JAPJEU+iescl-1):
play enters the node slave and the mesh (or the node) main
or: specified value of the second member (case without pairing)
3.10 Table
COCO
(addresses
JCOCO
)
It contains the memories of the state of preceding contact.
ZI (JCOCO)
= NDIM
: dimension of space (2 or 3)
ZI (JCOCO+1) = INDIC
: 0 if initialization
+1 if one added a connection
­ 1 if a connection were removed
ZI (JCOCO+2) = NBLIAC
: a number of active connections in the preceding state
ZI (JCOCO+3) = AJLIAI
: index in the list of the active connections of the last connection having been
calculated for the vector
CM1A
ZI (JCOCO+4) = SPLIAI
: index in the list of the active connections of the last correct line of
calculation of the matrix
T
1
.A
A.C
-
ZI (JCOCO+5) = LLF
: a number of connections of adherent friction in the preceding state
ZI (JCOCO+6) = LLF1
: a number of connections of adherent friction following the first direction
in the preceding state
ZI (JCOCO+7) = LLF
: a number of connections of adherent friction following the second direction
in the preceding state
3.11 Table
LIAC
(addresses
JLIAC
)
It contains the absolute numbers of the active connections of contact and adherent friction.
The list is not ordered
3.12 Table
CONVEC
(addresses
JVECC
)
It is fixed on table LIAC. It contains the type of the connection:
C0 if it is about a connection in contact
F0 if it is about a connection in adherent friction following the two directions of slips
F1 if it is about a connection in adherent friction following the first direction of slips
F2 if it is about a connection in adherent friction following the second direction of slips
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
15/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
3.13 Table
LIOT
(addresses
JLIOT
)
It contains the absolute numbers of the connections of contact-friction causing the appearance of a null pivot
at the time of the resolution. This connection will thus be removed system.
ZI (JLIOT) =
A number of connections of contact to null pivot
ZI (JLIOT+1)
with
ZI (JLIOT+NBLIAC)
=
Connections of contact to null pivot
ZI (JLIOT+NBLIAI+1) =
A number of connections of friction to null pivot in
directions 1 and 2
ZI (JLIOT+NBLIAI+2)
with
ZI (JLIOT+2 * NBLIAC+1)
=
Connections of friction to null pivot in
directions 1 and 2
ZI (JLIOT+2 * NBLIAI+2) =
A number of connections of friction to null pivot in
direction 1
ZI (JLIOT+2 * NBLIAI+3)
with
ZI (JLIOT+3 * NBLIAC+2)
=
Connections of friction to null pivot in the direction
1
ZI (JLIOT+3 * NBLIAI+3) =
A number of connections of friction to null pivot in
direction 2
ZI (JLIOT+3 * NBLIAI+4)
with
ZI (JLIOT+4 * NBLIAC+3)
=
Connections of friction to null pivot in the direction
2
Caution:
Each under-vector of
LIOT
is length
NBLIAI
, this is why one stores at the beginning of these
the last their working length.
3.14 Table
DRIVEN
(addresses
JMU
)
It contains the multipliers of Lagrange associated with contact-friction. Its maximum length is
6 * NESMAX
, but its effective length with a given iteration is based on the number of connections
active
NBLIAC
. It is organized as follows:
ZR (JMU)
with
ZR (JMU+NBLIAC-1)
=
Lagrange of the contact
ZR (JMU+NBLIAC)
with
ZR (JMU+2 * NBLIAC-1)
=
Lagrange of the adherent connections in the direction
1
ZR (JMU+2 * NBLIAC)
with
ZR (JMU+3 * NBLIAC-1) =
Lagrange of the adherent connections in the direction
2
ZR (JMU+3 * NBLIAC)
with
ZR (JMU+4 * NBLIAC-1) =
Lagrange of the sliding joints
ZR (JMU+6 * NBLIAC-1) =
Useful size for the resolution
3.15 Table
COEFMU
(addresses
JCMU
)
It contains the coefficient by which it is necessary to multiply the multiplier of Lagrange DRIVEN in the routine
algoco
before testing its sign. This coefficient is worth +1 in the case of a unilateral relation on
displacement, - 1 in the case of a unilateral relation on the pressure or the temperature of the elements
THM (this in order to be coherent with the fact that the hydraulic equation and the thermal equation of
problem coupled THM are multiplied by - 1).
3.16 Table
ATMU
(addresses
JATMU
)
It contains the nodal reactions of contact, i.e.
DRIVEN
With
T
.
, where
With
is the matrix of contact. Its
dimension is the total number of degrees of freedom of the problem, that is to say
NEQ
.
3.17 Tables
DELT0
and
DELTA
(addresses
JDELT0
and
JDELTA
)
They are auxiliary vectors, dimensioned with the total number of degrees of freedom
NEQ
, used in
the algorithm of active stresses.
background image
Code_Aster
®
Version
7.4
Titrate:
Structures of data related to contact-friction
Date:
06/12/04
Author (S):
Mr. ABBAS, NR. TARDIEU
Key
:
D4.06.14-D
Page
:
16/16
Data-processing manual of Description
D4.06 booklet: Structures related to the finite elements
HT-66/04/003/A
3.18 Variables
CM1A
,
MATR
and
STOC
CM1A
is a collection of
NBLIAI
objects length
NEQ
: each one of these objects contains one
column of the matrix
T
1
.A
C
-
, where
C
is the matrix of tangent rigidity (including/understanding the terms of
Lagrange) and
With
the matrix of contact. These vectors are used in the calculation of the matrix
T
1
.A
A.C
-
-
, stored in the matrix
MATR
, with a line storage of sky describes by the variable
STOC
. In these vectors and matrices, matrix A is reduced to the only active connections.

3.19 Variables
CM2 With
and
CM3 With
CM2 With
and
CM3 With
are collections of
NBLIAI
objects length
NEQ
: each one of these objects contains
a column of the tangent matrices of friction. For more precise details, to refer to the document
[R5.03.51].