Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
1/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
Organization (S):
EDF-R & D/AMA
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
Document: V6.03.121
SSNP121 Intégration of the terms of contact in 2D
and 3D
Summary:
This problem corresponds to a quasi-static analysis of a problem of mechanics with contact without
friction. One is interested particularly here in integration of the terms of contact. It is a question of studying two
identical blocks subjected to symmetrical imposed displacements.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
2/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
This test comprises three modelings in 2D (linear elements SEG2):
· modeling a: METHODE= `CONTINUE'. Method of integration by subelements
also proposed by Bathe [bib1] is used with three subelements,
· modeling b: METHODE= `CONTRAINTE',
· modeling C: METHODE= `PENALIZATION',
a modeling in 2D (linear elements SEG2 in with respect to quadratic elements SEG3):
· modeling H
: METHODE= `CONTINUE'. Surfaces of contact are made up
elements SEG2 in with respect to elements SEG3,
three modelings in 3D (quadratic elements):
· modeling D
: METHODE= `CONTINUE'. Surfaces of contact are made up
elements QUAD8 in with respect to elements QUAD8,
· modeling E: METHODE= `CONTINUE'. Surfaces of contact are made up of elements
TRIA6 in with respect to elements TRIA6,
· modeling F: METHODE= `CONTINUE'. Surfaces of contact are made up of elements
TRIA6 in with respect to elements QUAD8.
· modeling G: METHODE= `CONTINUE'. An alternative of modeling A where one tests
the initial activation of the statute of contact.
and three modelings in 3D (linear elements in with respect to quadratic elements):
· modeling I: METHODE= `CONTINUE'. Surfaces of contact are made up of elements
QUAD4 in with respect to elements QUAD8,
· modeling J: METHODE= `CONTINUE'. Surfaces of contact are made up of elements
TRIA3 in with respect to elements TRIA6,
· modeling K: METHODE= `CONTINUE'. Surfaces of contact are made up of elements
TRIA3 in with respect to elements QUAD8.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
3/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
1
Problem of reference
1.1 Geometry
Length has = 2 Mr.
Width B = 1 Mr.
O not medium of segment AB.
1.2
Material properties
Plates 1 and 2:
Poisson's ratio: 0.0
Young modulus: 2. 10
6
NR/m
2
1.3
Boundary conditions and loadings
Plate 1 is locked:
· On HG DX = 0 and DY = 0.
Plate 2 is subjected to an imposed displacement:
· On CD: DY = U
0
= 0.1 m and DX = 0.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
4/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution, analytical, can be deduced from a very simple calculation. Deformation
= U
0
/HD = 0.1/2. The pressure is worth E then
= 1.10
5
AP.
2.2
Results of reference
The contact pressure is constant and equal to 1.10
5
AP on all the surface of contact. In the same way it
vertical displacement (according to y) is constant on the surface of contact and equal to U
0
/2 = 0.05 Misters.
2.3 Reference
bibliographical
[1]
NR. EL-ABBASI and K.J. BATHE: “Stability and Patch Test Performance off Contact
Discretizations and has New Solution Algorithm ", Computers & Structures, 79, 1473-1486, 2001
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
5/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
3 Modeling
With
3.1
Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with the method CONTINUES for
processing of CONTACT with integration of the type “SIMPSON2”.
12 finite elements SEG2 are laid out on the initial surface of contact of plate 1 and only 11
on the surface of contact of other surface. By activating key word “SIMPSON2” three subelements
are used for the integration of the terms of contact. The mesh comprises in all 265 elements
QUAD4 for the two plates.
3.2
Characteristics of the mesh
A number of nodes:
313
A number of meshs and types:
265 QUAD4 and 132 SEG2
3.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“SIMPSON2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
4
Results of modeling A
4.1 Values
tested
Identification Reference
Aster
% difference
DY at point A
0.05
0.0499
0.007%
LAGR_C at point A
1.E+5
1.0008 E+5
0.08%
DY at the point O
0.05
0.0499
0.043%
LAGR_C at the point O
1.E+5
1.0035 E+5
0.351%
DY at the point B
0.05
0.0499
0.007%
LAGR_C at the point B
1.E+5
1.0008 E+5
0.08%
One as checks by a IMPR_RESU as the vertical pressure and displacements are identical on
all segment AB.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
6/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
5 Modeling
B
5.1
Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with the method FORCED for
processing of CONTACT.
5.2
Characteristics of the mesh
One uses the same mesh as for preceding modeling.
A number of nodes:
313
A number of meshs and types:
265 QUAD4 and 132 SEG2
5.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“FORCED”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
6
Results of modeling B
6.1 Values
tested
Identification Reference
Aster
% difference
DY at point A
0.05
0.05017
0.347%
SIYY at point A
1.E+5
9.9624 E+5
0.376%
DY at the point O
0.05
0.04916
1.666%
SIYY at the point O
1.E+5
1.0600 E+5
6.004%
DY at the point B
0.05
0.04916
1.666%
SIYY at the point B
1.E+5
1.0600 E+5
6.004%
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
7/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
7 Modeling
C
7.1
Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with the method PENALIZATION with
a coefficient of penalization of 1.E7 for the processing of the contact.
7.2
Characteristics of the mesh
One uses the same mesh as for preceding modeling.
A number of nodes:
313
A number of meshs and types:
265 QUAD4 and 132 SEG2
7.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
AFFE_CHAR_MECA CONTACT
E_N
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
8
Results of modeling C
8.1 Values
tested
Identification Reference
Aster
% difference
DY at point A
0.05
0.05060
2.121%
SIYY at point A
1.E+5
9.7875 E+5
2.125%
DY at the point O
0.05
0.04965
0.699%
SIYY at the point O
1.E+5
1.0741 E+5
7.413%
DY at the point B
0.05
0.04965
0.699%
SIYY at the point B
1.E+5
1.0741 E+5
7.413%
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
8/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
9 Modeling
D
9.1
Characteristics of modeling
One uses a modeling 3D for the solid elements with the method CONTINUES for
processing of CONTACT. Surfaces of contact in opposite consist of elements QUAD8.
9.2
Characteristics of the mesh
One uses the same mesh as for preceding modeling.
A number of nodes:
850
A number of meshs and types:
128 HEXA20 and 64 QUAD8
9.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
10 Results of modeling D
10.1 Values
tested
Identification Reference
Aster
% difference
DZ at point A
0.05
0.05
0.0%
SIYY at point A
1.E+5
1. E+5
0.0%
Max SIYY on the surface of contact
1.E+5
1. E+5
0.0%
Min SIYY on the surface of contact
1.E+5
1. E+5
0.0%
DZ at the point B
0.05
0.05
0.0%
SIYY at the point B
1.E+5
1. E+5
0.0%
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
9/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
11 Modeling
E
11.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with the method CONTINUES for
processing of CONTACT. Surfaces of contact in opposite consist of elements TRIA6.
11.2 Characteristics of the mesh
One uses the same mesh as for preceding modeling.
A number of nodes:
1010
A number of meshs and types:
256 PENTA15 and 128 TRIA6
11.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
12 Results of modeling E
12.1 Values
tested
Identification Reference
Aster
% difference
DZ at point A
0.05
0.050002
0.004%
SIYY at point A
1.E+5
1.001 E+5
0.01%
Max SIYY on the surface of contact
1.E+5
9.9992 E+4
0.008%
Min SIYY on the surface of contact
1.E+5
1.00015 E+5
0.015%
DZ at the point B
0.05
0.05
0.0%
SIYY at the point B
1.E+5
9.9996 E+4
0.003%
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
10/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
13 Modeling
F
13.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with the method CONTINUES for
processing of CONTACT. Surfaces of contact in opposite consist of elements TRIA6 in
with respect to meshs QUAD8.
13.2 Characteristics of the mesh
One uses the same mesh as for preceding modeling.
A number of nodes:
930
A number of meshs and types:
64 HEXA20, 128 PENTA15, 32 QUAD8 and 64 TRIA6
13.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
14 Results of modeling F
14.1 Values
tested
Identification Reference
Aster
% difference
DZ at point A
0.05
0.05
0.0%
SIYY at point A
1.E+5
1. E+5
0.0%
Max SIYY on the surface of contact
1.E+5
1. E+5
0.0%
Min SIYY on the surface of contact
1.E+5
1. E+5
0.0%
DZ at the point B
0.05
0.05
0.0%
SIYY at the point B
1.E+5
1. E+5
0.0%
14.2 Notice
The surface of contact slave is with a grid in TRIA6 and the surface of main contact is with a grid in
QUAD8. A use with elements QUAD8 for surface slave and of the TRIA6 for
surface main does not satisfy the conditions of compatibility necessary to the good integration of
terms of contact [R5.03.52]. In a general way, if one does not apprehend correctly this
concept of compatibility, one advises with the user not informed to use the same elements for
main mesh of surfaces and slave.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
11/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
15 Modeling
G
15.1 Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with the method CONTINUES for
processing of CONTACT with integration of the type “SIMPSON2”.
Plate 2 is subjected to a pressure E
= 1.10
5
AP on all surface CD. This pressure
is equivalent to a displacement imposed U
0
= 0.1 m on CD according to following formulas':
= U
0
/HD = - 0.1/2 and P
0
= E
= 1.10
5
AP
Note:
It should be taken care of the orientation of the group of mesh of plate 2 in such way that the normal
that is to say outgoing. One had to use the key word ORIENT_PEAU_2D.
Note:
In this modeling, one imposes a pressure on plate 2. So that the problem is
soluble, it is necessary CONTACT_INIT = “YES”. In the contrary case, with the first iteration of
first pitch of time, plate 2 would have a rigid movement of body.
15.2 Characteristics of the mesh
A number of nodes:
313
A number of meshs and types:
265 QUAD4 and 132 SEG2
15.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“SIMPSON2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
CONTACT_INIT
“YES”
AFFE_CHAR_MECA DDL_IMPO
AFFE_CHAR_MECA PRES_REP
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
Boundary conditions:
To avoid the rigid movements of body, the node D has a horizontal displacement no one.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
12/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
16 Results of modeling G
16.1 Values
tested
Identification Reference
Aster
% difference
DY at point A
0.05
0.0499
0.007%
LAGR_C at point A
1.E+5
1.0008 E+5
0.08%
DY at the point O
0.05
0.0499
0.043%
LAGR_C at the point O
1.E+5
1.0035 E+5
0.351%
DY at the point B
0.05
0.0499
0.007%
LAGR_C at the point B
1.E+5
1.0008 E+5
0.08%
One as checks by a IMPR_RESU as the vertical pressure and displacements are identical on
all segment AB.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
13/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
17 Modeling
H
17.1 Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with the method CONTINUES for
processing of the CONTACT between linear mixed elements and quadratic elements with presence
incompatibilities of mesh.
12 finite elements SEG2 are laid out on the initial surface of contact of the plate slave and
only 11 finite elements SEG3 on the surface of main contact. By activating the key word
“NCOTES2”, a diagram of the Newton-Dimensions type coupled to a technique of subdivision in under
elements was used for the integration of the terms of contact. The mesh comprises in all
144 elements QUAD4 for the plate slave opposite 121 elements QUAD8 for the plate Master.
17.2 Characteristics of the mesh
A number of nodes:
721
A number of meshs and types:
144 QUAD4, 121 QUAD8, 48 SEG2 and 44 SEG3
17.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
18 Results of modeling H
18.1 Values
tested
Identification Reference
Aster
% difference
DY at point A
0.05
0.0499
0.012%
LAGR_C at point A
1.E+5
1.00034 E+5
0.034%
DY at the point B
0.05
0.0499
0.012%
LAGR_C at the point B
1.E+5
1.00034 E+5
0.034%
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
14/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
19 Modeling
I
19.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with the method CONTINUES for
processing of the CONTACT between linear mixed elements and quadratic elements with presence
incompatibilities of mesh. Surfaces of contact consist of 5 elements QUAD4 in
with respect to 4 elements QUAD8. The diagram of integration used is of Newton-Dimensions type coupled to
a technique of subdivision in subelements.
19.2 Characteristics of the mesh
A number of nodes:
227
A number of meshs and types:
16 HEXA20, 25 HEXA8, 32 QUAD8 and 50 QUAD8
19.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
20 Results of modeling I
20.1 Values
tested
Identification Reference
Aster
% difference
DZ at point A
0.05
0.0499996
0.000061%
SIYY at point A
1.E+5
1.000045 E+5
0.005%
Max SIYY on the surface of contact
1.E+5
9.9988 E+4
- 0.012%
Min SIYY on the surface of contact
1.E+5
1.000093 E+5
0.009%
DZ at the point B
0.05
0.0499996
0.000061%
SIYY at the point B
1.E+5
1.000045 E+5
0.005%
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
15/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
21 Modeling
J
21.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with the method CONTINUES for
processing of the CONTACT between linear mixed elements and quadratic elements. Surfaces of
contact consist of elements TRIA3 in with respect to elements TRIA6. The mesh is
compatible. The diagram of integration used is of Newton-Dimensions type.
21.2 Characteristics of the mesh
A number of nodes:
630
A number of meshs and types:
128 PENTA6, 128 PENTA15, 64 TRIA3 and 64 TRIA6
21.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
22 Results of modeling J
22.1 Values
tested
Identification Reference
Aster
% difference
DZ at point A
0.05
0.05
0.0%
SIYY at point A
1.E+5
1.E+5
0.0%
Max SIYY on the surface of contact
1.E+5
1.E+5
0.0%
Min SIYY on the surface of contact
1.E+5
1.E+5
0.0%
DZ at the point B
0.05
0.05
0.0%
SIYY at the point B
1.E+5
1.E+5
0.0%
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
16/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
23 Modeling
K
23.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with the method CONTINUES for
processing of the CONTACT for the linear/quadratic mixed elements. Surfaces of contact in
opposite consist of elements TRIA3 in with respect to meshs QUAD8. The mesh is
compatible. The diagram of integration used is of Newton-Dimensions type.
23.2 Characteristics of the mesh
A number of nodes:
550
A number of meshs and types:
64 HEXA20, 128 PENTA6, 32 QUAD8 and 64 TRIA3
23.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT
STAT_NON_LINE COMP_ELAS
RELATION
“ELAS”
24 Results of modeling K
24.1 Values
tested
Identification Reference
Aster
% difference
DZ at point A
0.05
0.05
0.0%
SIYY at point A
1.E+5
1. E+5
0.0%
Max SIYY on the surface of contact
1.E+5
1. E+5
0.0%
Min SIYY on the surface of contact
1.E+5
1. E+5
0.0%
DZ at the point B
0.05
0.05
0.0%
SIYY at the point B
1.E+5
1. E+5
0.0%
24.2 Notice
The surface of contact slave is with a grid in TRIA3 and the surface of main contact is with a grid in
QUAD8. A use with elements QUAD4 for surface slave and of the TRIA6 for
surface main does not satisfy the conditions of compatibility necessary to the good integration of
terms of contact [R5.03.52].
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
17/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
25 Summary of the results
One seeks on this example very simple to test a new technique of integration of the terms of
contact based on the subdivision by subelements available only for the method
“CONTINUOUS”. This technique aims to attenuate the amplitude of oscillation of the pressure of
contact. In the case studied here, the pressure is constant on all the surface of contact (notice
that the Poisson's ratio is null).
One also implemented a new diagram of integration of the Newton-Dimensions type allowing of
to solve the problems of contact for Linear/Quadratic mixed elements with
incompatibility of mesh.
One notes thus that with the methods “FORCED” and “PENALIZATION” the solution present
nonphysical oscillations of about 6 to 7%. By using the “CONTINUOUS” method, them
oscillations disappear almost completely and the results obtained are very close to the solution
of reference (<0,5%). Moreover, one notes that when this technique is used, the interpenetration
(due to the use of master-slave modeling) decreases also significantly.
The “CONTINUOUS” method allows moreover a nonapproximate processing of surfaces of contact
with a grid with quadratic elements, whereas the “FORCED” method is used with one
linearization of the elements on the surface of contact. One advises with the users not informed on
compatibility of the elements to be used for the main mesh of surfaces and slave to use the same one
type of meshs for two surfaces.
Code_Aster
®
Version
8.2
Titrate:
SSNP121 Intégration of the terms of contact in 2D and 3D
Date
:
20/12/05
Author (S):
S. LAMARCHE, Mr. TORKHANI, NR. TARDIEU
Key
:
V6.03.121-C
Page:
18/18
Manual of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
Intentionally white left page.