Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
1/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
Organization (S):
EDF-R & D/AMA
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
Document: V5.03.103
SDNV103 - Impact of a bar of Taylor
elastoplastic
Summary:
One studies the impact rubbing of an elastoplastic bar on a rigid solid mass in nonlinear dynamics.
modeling includes/understands: contact, friction, elastoplasticity, great deformations.
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
2/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
1
Problem of reference
1.1 Geometry
1.2
Properties of material
E = 117. E3 MPa
naked = 0.35
= 8.93 E-9 G/mm
3
Y
= 100. MPa
E
T
= 400. MPa
Coefficient of friction of Coulomb:
µ
= 0.25
E
E
T
Y
1.3
Boundary conditions and loadings
The rigid foundation is completely locked throughout calculation.
The bar is free of any blocking.
There is a relation between unilateral contact and friction of Coulomb between the face inferior of the bar
and the higher face of the rigid foundation.
1.4 Conditions
initial
The bar is subjected at initial 227. E3 m/s.
B
With
32.4 mm
6.4 mm
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
3/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
2
Reference solution
2.1
Method of calculation
The reference solution comes from [bib1]. They are explicit axisymmetric calculations carried out with
various codes. One is unaware of almost all modeling: presence of contact, presence of
friction, coefficient of friction? In this measurement, one uses this reference in manner
indicative. The other tests will be of nonregression.
2.2
Sizes and results of reference
The sizes tested are:
·
radial displacement of point a:
mm
87
.
3
5
96
.
3
88
.
3
72
.
3
86
.
3
93
.
3
=
+
+
+
+
·
vertical displacement of point b:
mm
46
.
13
5
24
.
13
57
.
13
62
.
13
63
.
13
24
.
13
-
=
-
-
-
-
-
2.3
Uncertainties on the solution
Uncertainties on the reference solution are very important (see [§2.1]).
2.4 References
bibliographical
[1]
L. STAINIER, P. pH. PONTHOT: “Year improved broad one-point integration method for strain
elastoplastic analysis “, Comput. Methods Appl. Mech. Engrg. 118 (1994).
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
4/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
3 Modeling
With
3.1
Characteristics of modeling
3.2
Characteristics of the mesh
A number of nodes: 2850
Numbers and types of elements: 480 HEXA20, 200 PENTA15, 224 QUAD8, 6 TRIA6, 280 SEG3,
1 HEXA8, 6 QUAD4, 8 SEG2,101 POI1
3.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
METHOD
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
ALL
“YES”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR
RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
QUICKLY
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
5/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
3.4
Sizes tested and results
Being given the heaviness of modeling, one takes only some pitches of time and one carries out
only tests of not-regression.
3.5 Parameters
of execution
Version: 6.3.07
Machine: CLASTER
System: IRIX64
Overall dimension memory:
300 megawords
Time CPU To use: 720 seconds
4 Modeling
B
4.1
Characteristics of modeling
4.2
Characteristics of the mesh
A number of nodes: 359
Numbers and types of elements: 101 QUAD8, 55 SEG3
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
6/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
4.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
METHOD
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
ALL
“YES”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
QUICKLY
4.4
Sizes tested and results
Size tested
Reference:
Code_Aster
Difference (%)
radial displacement of
not A
3.87 2.92
24.55
vertical displacement of
not B
13.46 12.69
5.72
4.5 Parameters
of execution
Version: 6.3.07
Machine: CLASTER
System: IRIX64
Overall dimension memory:
16 megawords
Time CPU To use: 434. seconds
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
7/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
5 Modeling
C
5.1
Characteristics of modeling
5.2
Characteristics of the mesh
A number of nodes: 74
Numbers and types of elements: 4 HEXA20, 4 PENTA15, 10 QUAD8, 2 TRIA6, 14 SEG3,
1 HEXA8, 6 QUAD4, 8 SEG2,17 POI1
5.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION “COULOMB”
METHOD
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
ALL
“YES”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR
RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
QUICKLY
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
8/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
5.4
Sizes tested and results
For this very reduced modeling which has only one role of algorithmic control, one only tests
values of not-regression.
5.5 Parameters
of execution
Version: 6.3.07
Machine: CLASTER
System: IRIX64
Overall dimension memory:
300 megawords
Time CPU To use: 185 seconds
6 Modeling
D
6.1
Characteristics of modeling
6.2
Characteristics of the mesh
A number of nodes: 29
Numbers and types of elements: 5 QUAD8, 12 SEG3
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
9/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
6.3 Functionalities
tested
Controls
AFFE_CHAR_MECA CONTACT
FRICTION
“COULOMB”
METHOD
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
ALL
“YES”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
QUICKLY
6.4
Sizes tested and results
For this very reduced modeling which has only one role of algorithmic control, one only tests
values of not-regression.
6.5 Parameters
of execution
Version: 6.3.07
Machine: CLASTER
System: IRIX64
Overall dimension memory:
16 megawords
Time CPU To use: 37. seconds
Code_Aster
®
Version
6.2
Titrate:
SDNV103 - Impact of an elastoplastic bar of Taylor
Date:
15/04/03
Author (S):
NR. TARDIEU
Key
:
V5.03.103-A
Page:
10/10
Manual of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
7
Summary of the results
Present modelings differ from the reference of the literature by the taking into account from not
additional linearities (contact, friction), which explains the differences between their results
respective.
It is also noted that calculation 3D presents a overcost of enormous time CPU compared to
the axisymetry, which is explained at the same time by the greatest number of degrees of freedom but also by
processing of the friction which is much more complex in 3D than in 2D.