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
Organization (S): EDF-R & D/AMA
Handbook 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.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
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.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
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. 106 NR/m2
1.3
Boundary conditions and loadings
Plate 1 is blocked:
· On HG DX = 0 and DY = 0.
Plate 2 is subjected to an imposed displacement:
· On CD: DY = U0 = 0.1 m and DX = 0.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
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
= U0/HD = 0.1/2. The pressure is worth then E = 1.105 Pa.
2.2
Results of reference
The contact pressure is constant and equal to 1.105 Pa 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 U0/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
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
3 Modeling
With
3.1
Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with method CONTINUE 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 grid comprises in all 265 elements
QUAD4 for the two plates.
3.2
Characteristics of the grid
A number of nodes:
313
A number of meshs and types:
265 QUAD4 and 132 SEG2
3.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“SIMPSON2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
5 Modeling
B
5.1
Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with method CONTRAINTE for
processing of CONTACT.
5.2
Characteristics of the grid
One uses the same grid 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
METHODE
“CONTRAINTE”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIQUE”
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%
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
7 Modeling
C
7.1
Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with method PENALIZATION with
a coefficient of penalization of 1.E7 for the processing of the contact.
7.2
Characteristics of the grid
One uses the same grid 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
METHODE
“PENALIZATION”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIQUE”
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%
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
9 Modeling
D
9.1
Characteristics of modeling
One uses a modeling 3D for the solid elements with method CONTINUE for
processing of CONTACT. Surfaces of contact in opposite consist of elements QUAD8.
9.2
Characteristics of the grid
One uses the same grid 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
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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%
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
11 Modeling
E
11.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with method CONTINUE for
processing of CONTACT. Surfaces of contact in opposite consist of elements TRIA6.
11.2 Characteristics of the grid
One uses the same grid 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
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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%
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
13 Modeling
F
13.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with method CONTINUE for
processing of CONTACT. Surfaces of contact in opposite consist of elements TRIA6 in
with respect to meshs QUAD8.
13.2 Characteristics of the grid
One uses the same grid 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
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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 grid of surfaces and slave.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
15 Modeling
G
15.1 Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with method CONTINUE for
processing of CONTACT with integration of the type “SIMPSON2”.
Plate 2 is subjected to a pressure E = 1.105 Pa on all surface CD. This pressure
is equivalent to a displacement imposed U0 = 0.1 m on CD according to following formulas':
= U0/HD = - 0.1/2 and P0 = E = 1.105 Pa
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 = “OUI”. In the contrary case, with the first iteration of
first step of time, plate 2 would have a rigid movement of body.
15.2 Characteristics of the grid
A number of nodes:
313
A number of meshs and types:
265 QUAD4 and 132 SEG2
15.3 Functionalities
tested
AFFE_CHAR_MECA CONTACT
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“SIMPSON2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
AFFE_CHAR_MECA CONTACT
CONTACT_INIT
“OUI”
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.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
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.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
17 Modeling
H
17.1 Characteristics of modeling
One uses a modeling 2d_PLAN for the solid elements with method CONTINUE for
processing of the CONTACT between linear mixed elements and quadratic elements with presence
incompatibilities of grids.
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-Cotes type coupled to a technique of subdivision in under
elements was used for the integration of the terms of contact. The grid comprises in all
144 elements QUAD4 for the plate slave opposite 121 elements QUAD8 for the plate Master.
17.2 Characteristics of the grid
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
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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%
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
19 Modeling
I
19.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with method CONTINUE for
processing of the CONTACT between linear mixed elements and quadratic elements with presence
incompatibilities of grids. Surfaces of contact consist of 5 elements QUAD4 in
with respect to 4 elements QUAD8. The diagram of integration used is of Newton-Cotes type coupled to
a technique of subdivision in subelements.
19.2 Characteristics of the grid
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
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES2”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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%
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
21 Modeling
J
21.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with method CONTINUE 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 grids are
compatible. The diagram of integration used is of Newton-Cotes type.
21.2 Characteristics of the grid
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
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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%
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
23 Modeling
K
23.1 Characteristics of modeling
One uses a modeling 3D for the solid elements with method CONTINUE 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 grids are
compatible. The diagram of integration used is of Newton-Cotes type.
23.2 Characteristics of the grid
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
METHODE
“CONTINUE”
AFFE_CHAR_MECA CONTACT
INTEGRATION
“NCOTES”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NON”
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].
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
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
“CONTINUE”. 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-Cotes type allowing of
to solve the problems of contact for mixed elements Linéaires/Quadratiques with
incompatibility of grids.
One notes thus that with methods “CONTRAINTE” and “PENALIZATION” the solution present
nonphysical oscillations of about 6 to 7%. By using method “CONTINUE”, 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.
Method “CONTINUE” allows moreover a nonapproximate processing of surfaces of contact
with a grid with quadratic elements, whereas method “CONTRAINTE” 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 grid of surfaces and slave to use the same one
type of meshs for two surfaces.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
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
Intentionally white left page.
Handbook of Validation
V6.03 booklet: Nonlinear statics of the plane systems
HT-66/05/005/A
Outline document