Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
1/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
Organization (S):
EDF-R & D/AMA
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
V6.04.128 document
SSNV128 - Plate with contact and friction
on a rigid level
Summary:
This problem corresponds to a quasi-static analysis of a problem of mechanics with contact and friction.
A square plate is subjected to compressive forces and is compressed on 1 indeformable plan where it undergoes
forces of contact and friction. This test is based on results 2D suggested by the GRECO in
great deformations.
This test comprises 18 modelings making it possible to test:
·
modelings 2D (QUAD4, QUAD8, TRIA3 and TRIA6) and modelings 3D (HEXA8, HEXA20,
HEXA27, PENTA6, PENTA15, TETRA4 and TETRA10),
·
methods of processing of the contact with friction
`
LAGRANGIAN
'
,
`
PENALIZATION
'
and
`
CONTINUOUS
'
for the control CONTACT,
·
the geometrical reactualization:
`
AUTOMATIC
'
,
`
Control
'
or
`
WITHOUT
'
,
·
pairing
`
NODAL
'
or
`
MAIT_ESCL
'
,
·
the search of the nodes in contact with
`
NOEUD_BOUCLE
'
or
`
NOEUD_VOISIN
'
,
·
the option
`
MAIT
'
or
`
MAIT_ESCL
'
for the key word NORMAL,
·
the use of key word VECT_Y,
·
the use of key word COEF_MATR_FROT,
·
coefficients of penalization for the elements of contact: E_T and E_N,
·
key words of exclusion of a group of node: SANS_ (GROUP_) NO and VECT_Y.
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
2/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
1
Problem of reference
1.1 Geometry
before deformation
locked face
in X = 40 mm
With
40 mm
40 mm
F
Surface contact
with friction
X
y
Z
F
Thickness of the plate ep = 1 Misters.
Side of the plate has = 40 Misters.
Position of the points of reference under the surface of contact (mm)
X y Z
To 0 0 0
B 1.25
0 0
C 5. 0 0
D 7.5 0 0
E 11.25 0 0
1.2
Material properties
Plate:
Poisson's ratio: 0.2
Young modulus: 1.3 10
11
NR/m
2
Jig (only if it is modelized by of the same elements dimension than the plate):
Poisson's ratio: 0.2
Young modulus:
E
(10
16
NR/m
2
).
The coefficient of friction under the rigid plan is
µ = 1.
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
3/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
1.3
Boundary conditions and loadings
The jig, when it is of N-1 dimension compared to the dimension of the plate, is locked:
·
by a complete embedding.
The jig, when it is of the same dimension than the plate, is locked:
·
in the plan
X
= 40 mm for displacements according to
X
(symmetry of the problem),
·
by an embedding of its lower face,
The plate is locked:
·
in the plan
X
= 40 mm for displacements according to
X
(symmetry of the problem),
·
according to Y, with the node located at the intersection of, the face symmetry plane of contact and plan
Z = 0 to prevent the movements of rigid body,
In 3D, to bring back itself to a problem 2D:
·
displacement following Z is locked for all the nodes,
The plate is subjected to two pressures distributed:
·
a vertical acting under the face of the top: F = 5 daN/mm
2
,
·
horizontal acting under the face initially in
X
= 0, F = 15.daN/mm
2
.
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution comes from results obtained by an average on other computer codes
[bib1].
2.2
Results of reference
Tangential displacements (according to
X
) at points A B C D E of the surface of contact.
2.3 Reference
bibliographical
[1]
R.A. FEIJOO H.J.C. BARBOSA and NR. ZOUAIN “Numerical formulations for contact problems
with friction “Newspaper off Theorical and Applied Mechanics GAUTHIER-VILLARS
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
4/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
3 Modeling
With
3.1
Characteristics of modeling
Modeling:
D_PLAN
to test the operand “COULOMB” of the key word FRICTION under
key word factor CONTACT for elements SEG2.
The plate is with a grid with elements QUAD4 and the jig with elements SEG2.
3.2
Characteristics of the mesh
A number of nodes: 396
A number of meshs and types: 352 QUAD4 for the plate and 32 SEG2 for the jig.
3.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
STAT_NON_LINE NEWTON
“MATRIX”
“TANGENT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
STAT_NON_LINE NEWTON
“MATRIX”
“TANGENT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
5/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
4
Results of modeling A
4.1
Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.487
DX at the point B
2.72E5
2.708 E5
- 0.439
DX at the point C
2.28E5
2.274 E5
- 0.259
DX at the point D
1.98E5
1.972 E5
- 0.366
DX at the point E
1.50E5
1.536 E5
2.43
4.2
Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846E5
- 0.487
DX at the point B
2.72E5
2.708E5
- 0.439
DX at the point C
2.28E5
2.274E5
- 0.259
DX at the point D
1.98E5
1.973E5
- 0.366
DX at the point E
1.50E5
1.536E5
2.43
4.3
Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.871E5
0.398
DX at the point B
2.72E5
2.733E5
0.492
DX at the point C
2.28E5
2.298E5
0.819
DX at the point D
1.98E5
1.998E5
0.936
DX at the point E
1.50E5
1.566E5
4.38
4.4
Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.871E5
0.398
DX at the point B
2.72E5
2.733E5
0.492
DX at the point C
2.28E5
2.298E5
0.819
DX at the point D
1.98E5
1.998E5
0.936
DX at the point E
1.50E5
1.566E5
4.38
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
6/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
5 Modeling
B
5.1
Characteristics of modeling
Modeling:
D_PLAN
to test the operand “COULOMB” of the key word FRICTION under
key word factor CONTACT for elements SEG3.
The plate is with a grid with elements QUAD8 and the jig with elements SEG3.
5.2
Characteristics of the mesh
A number of nodes: 558
A number of meshs and types: 160 QUAD8 for the plate and 1 SEG3 for the jig.
5.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“CONTROL”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
STAT_NON_LINE NEWTON
“MATRIX”
“TANGENT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
and
E_N
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“WITHOUT”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
STAT_NON_LINE NEWTON
“MATRIX”
“TANGENT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
7/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
6
Results of modeling B
6.1
Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.859 E5
- 0.005
DX at the point B
2.72E5
2.720 E5
0.028
DX at the point C
2.28E5
2.283 E5
0.118
DX at the point D
1.98E5
1.979 E5
- 0.042
DX at the point E
1.50E5
1.541 E5
2.79
6.2
Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.859 E5
- 0.005
DX at the point B
2.72E5
2.729 E5
0.028
DX at the point C
2.28E5
2.283 E5
0.118
DX at the point D
1.98E5
1.979 E5
- 0.042
DX at the point E
1.50E5
1.541 E5
2.709
6.3
Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.874 E5
0.486
DX at the point B
2.72E5
2.735 E5
0.545
DX at the point C
2.28E5
2.296 E5
0.708
DX at the point D
1.98E5
1.993 E5
0.667
DX at the point E
1.50E5
1.556 E5
3.767
6.4
Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.874 E5
0.486
DX at the point B
2.72E5
2.735 E5
0.545
DX at the point C
2.28E5
2.296 E5
0.708
DX at the point D
1.98E5
1.993 E5
0.667
DX at the point E
1.50E5
1.556 E5
3.767
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
8/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
7 Modeling
C
7.1
Characteristics of modeling
Modeling:
D_PLAN
to test the operand “COULOMB” of the key word FRICTION under
key word factor CONTACT for elements SEG2.
The plate is with a grid with elements TRIA3 and the jig with elements SEG2.
7.2
Characteristics of the mesh
A number of nodes: 431
A number of meshs and types: 732 TRIA3 for the plate and 32 SEG2 for the jig.
7.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“WITHOUT”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
STAT_NON_LINE NEWTON
“MATRIX”
“TANGENT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
and
E_N
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
STAT_NON_LINE NEWTON
“MATRIX”
“TANGENT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
9/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
8
Results of modeling C
8.1
Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.848 E5
- 0.420
DX at the point B
2.72E5
2.712 E5
- 0.300
DX at the point C
2.28E5
2.279 E5
- 0.042
DX at the point D
1.98E5
1.978 E5
- 0.095
DX at the point E
1.50E5
1.538 E5
2.58
8.2
Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.848 E5
- 0.420
DX at the point B
2.72E5
2.712 E5
- 0.300
DX at the point C
2.28E5
2.279 E5
- 0.042
DX at the point D
1.98E5
1.978 E5
- 0.095
DX at the point E
1.50E5
1.538 E5
2.528
8.3
Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.873 E5
0.462
DX at the point B
2.72E5
2.737 E5
0.630
DX at the point C
2.28E5
2.303 E5
1.023
DX at the point D
1.98E5
2.004 E5
1.193
DX at the point E
1.50E5
1.567 E5
4.453
8.4
Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.873 E5
0.462
DX at the point B
2.72E5
2.737 E5
0.630
DX at the point C
2.28E5
2.303 E5
1.023
DX at the point D
1.98E5
2.004 E5
1.193
DX at the point E
1.50E5
1.567 E5
4.453
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
10/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
9 Modeling
D
9.1
Characteristics of modeling
Modeling:
D_PLAN
to test the operand “COULOMB” of the key word FRICTION under
key word factor CONTACT for elements SEG3.
The plate is with a grid with elements TRIA6 and the jig with elements SEG3.
9.2
Characteristics of the mesh
A number of nodes: 763
A number of meshs and types: 325 TRIA6 for the plate and 32 SEG3 for the jig.
9.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
11/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
10 Results of modeling D
10.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.858 E5
- 0.065
DX at the point B
2.72E5
2.721 E5
0.050
DX at the point C
2.28E5
2.286 E5
0.279
DX at the point D
1.98E5
1.980 E5
0.027
DX at the point E
1.50E5
1.541 E5
2.751
10.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.858 E5
- 0.065
DX at the point B
2.72E5
2.721 E5
0.050
DX at the point C
2.28E5
2.286 E5
0.279
DX at the point D
1.98E5
1.980 E5
0.027
DX at the point E
1.50E5
1.541 E5
2.751
10.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.872 E5
0.428
DX at the point B
2.72E5
2.736 E5
0.573
DX at the point C
2.28E5
2.299 E5
0.860
DX at the point D
1.98E5
1.994 E5
0.724
DX at the point E
1.50E5
1.557 E5
3.792
10.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.872 E5
0.428
DX at the point B
2.72E5
2.735 E5
0.573
DX at the point C
2.28E5
2.299 E5
0.860
DX at the point D
1.98E5
1.994 E5
0.724
DX at the point E
1.50E5
1.557 E5
3.792
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
12/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
11 Modeling
E
11.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements QUAD4.
The plate is with a grid with elements HEXA8 and the jig with elements QUAD4.
11.2 Characteristics of the mesh
A number of nodes: 792
A number of meshs and types: 320 HEXA8 for the plate and 32 QUAD4 for the jig.
11.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
and
E_N
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
13/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
12 Results of modeling E
12.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.487
DX at the point B
2.72E5
2.708 E5
- 0.439
DX at the point C
2.28E5
2.274 E5
- 0.259
DX at the point D
1.98E5
1.973 E5
- 0.366
DX at the point E
1.50E5
1.536 E5
2.429
12.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.487
DX at the point B
2.72E5
2.708 E5
- 0.439
DX at the point C
2.28E5
2.274 E5
- 0.259
DX at the point D
1.98E5
1.973 E5
- 0.366
DX at the point E
1.50E5
1.536 E5
2.429
12.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.491
DX at the point B
2.72E5
2.708 E5
- 0.444
DX at the point C
2.28E5
2.274 E5
- 0.262
DX at the point D
1.98E5
1.973 E5
- 0.368
DX at the point E
1.50E5
1.536 E5
2.428
12.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.491
DX at the point B
2.72E5
2.708 E5
- 0.444
DX at the point C
2.28E5
2.274 E5
- 0.262
DX at the point D
1.98E5
1.973 E5
- 0.368
DX at the point E
1.50E5
1.536 E5
2.428
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
14/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
13 Modeling
F
13.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements QUAD8.
The plate is with a grid with elements HEXA20 and the jig with elements QUAD8.
13.2 Characteristics of the mesh
A number of nodes: 1316
A number of meshs and types: 160 HEXA20 for the plate and 1 QUAD8 for the jig.
13.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
VECT_Y
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
15/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
14 Results of modeling F
14.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.912 E5
1.817
DX at the point B
2.72E5
2.766 E5
1.680
DX at the point C
2.28E5
2.310 E5
1.330
DX at the point D
1.98E5
2.002 E5
1.135
DX at the point E
1.50E5
1.559 E5
3.996
14.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.912 E5
1.817
DX at the point B
2.72E5
2.766 E5
1.680
DX at the point C
2.28E5
2.310 E5
1.330
DX at the point D
1.98E5
2.002 E5
1.135
DX at the point E
1.50E5
1.539 E5
3.996
14.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.911 E5
1.804
DX at the point B
2.72E5
2.765 E5
1.666
DX at the point C
2.28E5
2.309 E5
1.312
DX at the point D
1.98E5
2.002 E5
1.112
DX at the point E
1.50E5
1.559 E5
3.963
14.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.911 E5
1.804
DX at the point B
2.72E5
2.765 E5
1.666
DX at the point C
2.28E5
2.309 E5
1.312
DX at the point D
1.98E5
2.002 E5
1.112
DX at the point E
1.50E5
1.559 E5
3.963
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
16/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
15 Modeling
G
15.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements QUAD8.
The plate is with a grid with elements HEXA27 and the jig with elements QUAD8.
15.2 Characteristics of the mesh
A number of nodes: 408
A number of meshs and types: 24 HEXA27 for the plate and 8 QUAD8 for the jig.
15.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
17/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
16 Results of modeling G
16.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.330
DX at the point C
2.28E5
2.275 E5
- 0.213
16.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.330
DX at the point C
2.28E5
2.275 E5
- 0.213
16.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.324
DX at the point C
2.28E5
2.275 E5
- 0.205
16.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.324
DX at the point C
2.28E5
2.275 E5
- 0.205
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
18/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
17 Modeling
H
17.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements QUAD4.
The plate is with a grid with elements PENTA6 and the jig with elements QUAD4.
17.2 Characteristics of the mesh
A number of nodes: 670
A number of meshs and types: 323 PENTA6 for the plate and 1 QUAD4 for the jig.
17.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
19/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
18 Results of modeling H
18.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.320
DX at the point B
2.72E5
2.716 E5
- 0.140
DX at the point C
2.28E5
2.284 E5
0.159
DX at the point D
1.98E5
1.978 E5
- 0.082
DX at the point E
1.50E5
1.540 E5
2.678
18.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.320
DX at the point B
2.72E5
2.716 E5
- 0.140
DX at the point C
2.28E5
2.284 E5
0.159
DX at the point D
1.98E5
1.978 E5
- 0.082
DX at the point E
1.50E5
1.540 E5
2.678
18.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.320
DX at the point B
2.72E5
2.716 E5
- 0.140
DX at the point C
2.28E5
2.284 E5
0.159
DX at the point D
1.98E5
1.978 E5
- 0.082
DX at the point E
1.50E5
1.540 E5
2.678
18.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.851 E5
- 0.320
DX at the point B
2.72E5
2.716 E5
- 0.140
DX at the point C
2.28E5
2.284 E5
0.159
DX at the point D
1.98E5
1.978 E5
- 0.082
DX at the point E
1.50E5
1.540 E5
2.678
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
20/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
19 Modeling
I
19.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements QUAD8.
The plate is with a grid with elements PENTA15 and the jig with elements QUAD8.
19.2 Characteristics of the mesh
A number of nodes: 1297
A number of meshs and types: 158 PENTA15 for the plate and 32 QUAD8 for the jig.
19.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
21/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
20 Results of modeling I
20.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.865 E5
0.179
DX at the point B
2.72E5
2.727 E5
0.277
DX at the point C
2.28E5
2.290 E5
0.453
DX at the point D
1.98E5
1.983 E5
0.164
DX at the point E
1.50E5
1.543 E5
2.871
20.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.865 E5
0.179
DX at the point B
2.72E5
2.727 E5
0.277
DX at the point C
2.28E5
2.290 E5
0.453
DX at the point D
1.98E5
1.983 E5
0.165
DX at the point E
1.50E5
1.543 E5
2.871
20.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.866 E5
0.207
DX at the point B
2.72E5
2.728 E5
0.306
DX at the point C
2.28E5
2.291 E5
0.486
DX at the point D
1.98E5
1.984 E5
0.205
DX at the point E
1.50E5
1.544 E5
2.932
20.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.866 E5
0.207
DX at the point B
2.72E5
2.728 E5
0.306
DX at the point C
2.28E5
2.291 E5
0.486
DX at the point D
1.98E5
1.984 E5
0.205
DX at the point E
1.50E5
1.544 E5
2.932
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
22/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
21 Modeling
J
21.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements TRIA3.
The plate is with a grid with elements TETRA4 and the jig with elements TRIA3.
21.2 Characteristics of the mesh
A number of nodes: 441
A number of meshs and types: 980 TETRA4 for the plate and 64 TRIA3 for the jig.
21.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
23/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
22 Results of modeling J
22.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.831 E5
- 1.005
DX at the point B
2.72E5
2.698 E5
- 0.794
DX at the point C
2.28E5
2.277 E5
- 0.096
DX at the point D
1.98E5
1.977 E5
- 0.125
DX at the point E
1.50E5
1.539 E5
2.624
22.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.831 E5
- 1.005
DX at the point B
2.72E5
2.698 E5
- 0.794
DX at the point C
2.28E5
2.270 E5
- 0.096
DX at the point D
1.98E5
1.977 E5
- 0.125
DX at the point E
1.50E5
1.539 E5
2.624
22.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.831 E5
- 1.007
DX at the point B
2.72E5
2.698 E5
- 0.796
DX at the point C
2.28E5
2.278 E5
- 0.098
DX at the point D
1.98E5
1.977 E5
- 0.126
DX at the point E
1.50E5
1.539 E5
2.624
22.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.831 E5
- 1.007
DX at the point B
2.72E5
2.698 E5
- 0.796
DX at the point C
2.28E5
2.278 E5
- 0.098
DX at the point D
1.98E5
1.977 E5
- 0.126
DX at the point E
1.50E5
1.539 E5
2.624
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
24/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
23 Modeling
K
23.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements TRIA6.
The plate is with a grid with elements TETRA10 and the jig with elements TRIA6.
23.2 Characteristics of the mesh
A number of nodes: 1236
A number of meshs and types: 526 TETRA10 for the plate and 32 TRIA6 for the jig.
23.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“PENALIZATION”
E_T
and
E_N
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
“NODAL” PAIRING
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“Control”
NB_REAC_GEOM
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
25/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
24 Results of modeling K
24.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.876 E5
0.551
DX at the point C
2.28E5
2.289 E5
0.435
24.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.875 E5
0.551
DX at the point C
2.28E5
2.289 E5
0.435
24.3 Values tested by penalized method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.876 E5
0.581
DX at the point C
2.28E5
2.291 E5
0.469
24.4 Values tested by penalized method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.876 E5
0.581
DX at the point C
2.28E5
2.291 E5
0.469
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
26/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
25 Modeling
L
25.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements QUAD4.
The plate is with a grid with elements HEXA8 and the jig with elements QUAD4. 32 finite elements
QUAD4 are laid out on the faces of contact (2D) of the plate and the jig. The plate is tilted with
45° (rotation around the axis Y).
25.2 Characteristics of the mesh
A number of nodes: 792
A number of meshs and types: 320 HEXA8 for the plate and 32 QUAD4 for the jig.
25.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
VECT_Y
AFFE_CHAR_MECA CONTACT
COEF_MATR_FROT
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“AUTOMATIC”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
27/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
26 Results of modeling L
26.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.02E5
2.012 E5
- 0.377
DX at the point B
1.92E5
1.915 E5
- 0.271
DX at the point C
1.61E5
1.608 E5
- 0.126
DX at the point D
1.40E5
1.395 E5
- 0.363
DX at the point E
1.06E5
1.086 E5
2.492
26.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.02E5
2.012 E5
- 0.377
DX at the point B
1.92E5
1.915 E5
- 0.271
DX at the point C
1.61E5
1.608 E5
- 0.126
DX at the point D
1.40E5
1.395 E5
- 0.363
DX at the point E
1.06E5
1.086 E5
2.492
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
28/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
27 Modeling
M
27.1 Characteristics of modeling
Modeling:
D_PLAN
to test the operand
“COULOMB”
key word
FRICTION
key word
factor
CONTACT
for elements SEG2.
The plate is with a grid with elements QUAD4 and the jig with elements SEG2.
27.2 Characteristics of the mesh
A number of nodes: 396
A number of meshs and types: 352 QUAD4 for the plate and 32 SEG2 for the jig.
27.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION “NODE”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
SEUIL_INIT
0
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI 1
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI 4
AFFE_CHAR_MECA CONTACT
ITER_FROT_MAXI 10
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT 1
AFFE_CHAR_MECA CONTACT
COEF_REGU_FROT 1
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
STAT_NON_LINE COMP_ELAS RELATION
“ELAS”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
TEST_RESU
VALE_CONT
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
29/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
28 Results of modeling M
28.1 Values tested by the method continues “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.484
DX at the point B
2.72E5
2.708 E5
- 0.436
DX at the point C
2.28E5
2.274 E5
- 0.255
DX at the point D
1.98E5
1.973 E5
- 0.365
DX at the point E
1.50E5
1.536 E5
2.429
Identification Reference
Aster %
difference
PLAY in N365
0.0
- 1.6156E-27
CONT in N365
1
1
RN in N365
1.0486E5
1.0486E5
- 6.04E-4
28.2 Values tested by the method continues “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.484
DX at the point B
2.72E5
2.708 E5
- 0.436
DX at the point C
2.28E5
2.274 E5
- 0.255
DX at the point D
1.98E5
1.973 E5
- 0.365
DX at the point E
1.50E5
1.536 E5
2.429
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
30/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
29 Modeling
NR
29.1 Characteristics of modeling
Modeling:
3D
to test the operand
“COULOMB”
key word
FRICTION
key word factor
CONTACT
for elements QUAD4.
The plate is with a grid with elements HEXA8 and the jig with elements QUAD4. 32 finite elements
QUAD4 are laid out on the faces of contact (2D) of the plate and the jig.
29.2 Characteristics of the mesh
A number of nodes: 792
A number of meshs and types: 320 HEXA8 for the plate and 32 QUAD4 for the jig.
29.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION “NODE”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
SEUIL_INIT
0
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI 1
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI 4
AFFE_CHAR_MECA CONTACT
ITER_FROT_MAXI 8
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT 100
AFFE_CHAR_MECA CONTACT
COEF_REGU_FROT 100
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
STAT_NON_LINE COMP_ELAS RELATION
“ELAS”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
31/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
30 Results of modeling NR
30.1 Values tested by the method continues “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.922 E5
2.174
DX at the point B
2.72E5
2.775 E5
2.024
DX at the point C
2.28E5
2.319 E5
1.727
DX at the point D
1.98E5
2.012 E5
1.617
DX at the point E
1.50E5
1.569 E5
4.652
30.2 Values tested by the method continues “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.922 E5
2.174
DX at the point B
2.72E5
2.775 E5
2.024
DX at the point C
2.28E5
2.319 E5
1.727
DX at the point D
1.98E5
2.012 E5
1.617
DX at the point E
1.50E5
1.569 E5
4.652
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
32/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
31 Modeling
O
31.1 Characteristics of modeling
Modeling:
D_PLAN
to test the operand
“COULOMB”
key word
FRICTION
key word
factor
CONTACT
for elements SEG2.
The plate is with a grid with elements QUAD4 and the jig with elements SEG2.
31.2 Characteristics of the mesh
A number of nodes: 396
A number of meshs and types: 352 QUAD4 for the plate and 32 SEG2 for the jig.
31.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“LAGRANGIAN”
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
AFFE_CHAR_MECA CONTACT
REAC_GEOM
“WITHOUT”
STAT_NON_LINE COMP_ELAS RELATION
“ELAS”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
“LDLT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
33/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
32 Results of modeling O
32.1 Values tested by Lagrangian method “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.492
DX at the point B
2.72E5
2.708 E5
- 0.444
DX at the point C
2.28E5
2.274 E5
- 0.262
DX at the point D
1.98E5
1.973 E5
- 0.369
DX at the point E
1.50E5
1.536 E5
2.427
32.2 Values tested by Lagrangian method “LDLT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.846 E5
- 0.492
DX at the point B
2.72E5
2.708 E5
- 0.444
DX at the point C
2.28E5
2.274 E5
- 0.262
DX at the point D
1.98E5
1.973 E5
- 0.369
DX at the point E
1.50E5
1.536 E5
- 0.427
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
34/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
33 Modeling
P
33.1 Characteristics of modeling
Modeling:
D_PLAN
to test the operand
“COULOMB”
key word
FRICTION
key word
factor
CONTACT
for elements SEG3.
The plate is with a grid with elements QUAD8 and the jig with elements SEG3.
33.2 Characteristics of the mesh
A number of nodes: 620
A number of meshs and types: 160 QUAD8 for the plate and 106 SEG3 for the jig.
33.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION “NODE”
AFFE_CHAR_MECA CONTACT
MODL_AXIS
“NOT”
AFFE_CHAR_MECA CONTACT
SEUIL_INIT
100.
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI 3
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI 1
AFFE_CHAR_MECA CONTACT
ITER_FROT_MAXI 10
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT 100.
AFFE_CHAR_MECA CONTACT
COEF_REGU_FROT 100.
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
STAT_NON_LINE COMP_ELAS RELATION
“ELAS”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
35/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
34 Results of modeling P
34.1 Values tested by the method continues “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.866 E5
- 0.013
DX at the point B
2.72E5
2.720 E5
- 0.020
DX at the point C
2.28E5
2.282 E5
- 0.109
DX at the point D
1.98E5
1.979 E5
- 0.054
DX at the point E
1.50E5
1.54 E5
2.692
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
36/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
35 Modeling
Q
35.1 Characteristics of modeling
Modeling:
D_PLAN
to test the key words
SANS_GROUP_NO (or SANS_NO)
and VECT_Y of
key word factor
CONTACT
for elements SEG2.
The plate is with a grid with elements QUAD4 and the jig with elements SEG2.
Vector VECT_Y gives the direction of the slip, whose direction is given compared to
normal vector
NR
(directed towards the interior of the main field) and so that doublet (VECT_Y,
nv
) form a direct base.
35.2 Characteristics of the mesh
A number of nodes: 396
A number of meshs and types: 352 QUAD4 for the plate and 32 SEG2 for the jig.
35.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION “NODE”
AFFE_CHAR_MECA CONTACT
CONTACT_INIT “YES”
AFFE_CHAR_MECA CONTACT
SANS_GROUP_NO
“PS”
AFFE_CHAR_MECA CONTACT
VECT_Y
(- 1., 0., 0.)
AFFE_CHAR_MECA CONTACT
SEUIL_INIT
0.1
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI 2
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI 4
AFFE_CHAR_MECA CONTACT
ITER_FROT_MAXI 10
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT 1.
AFFE_CHAR_MECA CONTACT
COEF_REGU_FROT 1.
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_BOUCLE”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
STAT_NON_LINE COMP_ELAS RELATION
“ELAS”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
37/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
36 Results of modeling Q
36.1 Values tested by the method continues “MULT_FRONT”
Identification Reference
Aster %
difference
PLAY at the point F
0
2.832 E27
2.83 E27
CONTACT at the point F
1.
1.
0
RN at the point F
1.486 E+5
1.48 E+5
- 2.16 E4
DX at point A
2.86E5
2.845 E5
- 0.492
DX at the point B
2.72E5
2.707 E5
- 0.444
DX at the point C
2.28E5
2.274 E5
- 0.262
DX at the point D
1.98E5
1.972 E5
- 0.369
DX at the point E
1.50E5
1.536 E5
2.427
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
38/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
37 Modeling
R
37.1 Characteristics of modeling
Modeling:
D_PLAN
to test the key words
SANS_GROUP_NO (or SANS_NO)
and VECT_Y of
key word factor
CONTACT
for elements QUAD4.
The plate and the jig are with a grid with elements HEXA8. A line of 32 finite elements QUAD4
is laid out on each face of contact (2D) of the plate and the jig.
Vector VECT_Y gives the direction of the slip, corresponding to one of the vectors
tangent
2
,
1
=
I
I
R
, where the trihedron
(
)
nv
v
v
,
,
2
1
form a direct base. It is given here
automatically.
37.2 Characteristics of the mesh
A number of nodes: 792
A number of meshs and types: 320 HEXA8 for the plate and the jig and 64 QUAD4 for the contact between
plate and the jig.
37.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
COULOMB
AFFE_CHAR_MECA CONTACT
METHOD
“CONTINUES”
AFFE_CHAR_MECA CONTACT
INTEGRATION “NODE”
AFFE_CHAR_MECA CONTACT
SANS_GROUP_NO
“PP”
AFFE_CHAR_MECA CONTACT
ITER_GEOM_MAXI 1
AFFE_CHAR_MECA CONTACT
ITER_CONT_MAXI 4
AFFE_CHAR_MECA CONTACT
ITER_FROT_MAXI 12
AFFE_CHAR_MECA CONTACT
COEF_REGU_CONT 100
AFFE_CHAR_MECA CONTACT
COEF_REGU_FROT 100
AFFE_CHAR_MECA CONTACT
PAIRING “MAIT_ESCL”
AFFE_CHAR_MECA CONTACT
SEEK
“NOEUD_VOISIN”
AFFE_CHAR_MECA CONTACT
NORMAL
“MAIT”
STAT_NON_LINE COMP_ELAS RELATION
“ELAS”
STAT_NON_LINE SOLVEUR
“MULT_FRONT”
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
39/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
38 Results of modeling R
38.1 Values tested by the method continues “MULT_FRONT”
Identification Reference
Aster %
difference
DX at point A
2.86E5
2.97E5
3.846
DX at the point B
2.72E5
2.82E5
3.676
DX at the point C
2.28E5
2.36E5
2.330
DX at the point D
1.98E5
2.06E5
3.340
DX at the point E
1.50E5
1.62E5
8.000
The relative error increases slightly in comparison with modeling NR, in particular for
not E.
Code_Aster
®
Version
8.2
Titrate:
SSNV128 - Plate with contact and friction on a rigid level
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. TORKHANI
Key
:
V6.04.128-D
Page:
40/40
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
39 Summary of the results
The results obtained on the whole of modelings of this case test are satisfactory, as well in
2D that in 3D. On points A, B, C and D, one obtains a relative error lower than 1% compared to
results of the GRECO. On the other hand, for the point E, the relative error, about 2,5%, remains
acceptable.
From a time point of view calculating, one notes that the method by penalization is in general more
rapid that Lagrangian method. One notes, on the other hand, an ease of use of the method
Lagrangian (cf Doc. [R5.03.01]) compared to the penalized method, since convergence and
quality of the results obtained by this last method are conditioned by the coefficient of
penalization which is automatically given in the Lagrangian method.
From a qualitative point of view:
·
in 2D, one notes a positive effect on the results of the quadratic mesh compared to
linear mesh (cf modelings A and B),
·
in 3D, one notes that calculations with quadratic mesh does not improve the results.
This is explained by the fact why one decreased by 50% the number of meshs compared to
linear mesh to have a number of nodes are equivalent.
As regards modeling L, inclined plate, one notes a good convergence and results
satisfactory.
Lastly, for the method continues (modelings M, NR, P and Q), one obtains in 3D results with
relative errors slightly more important than those obtained by the methods Lagrangian and
penalized. These errors can come from the choice of the coefficients (cf the documentation of use
controls AFFE_CHAR_MECA and AFFE_CHAR_MECA_F for more information on these
coefficients). In 2D, one notes a convergence slower than the methods Lagrangian and penalized.
Lastly, the exclusion of the nodes rubbing in redundancy with the boundary conditions imposed (pivot
no one) using SANS_GROUP_NO and VECT_Y functions correctly in 2D (modeling Q) but
increase the level of error in 3D (modeling R).