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
Organization (S): EDF-R & D/AMA
Handbook 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.
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
1
Problem of reference
1.1 Geometry
6.4 mm
B
32.4 mm
With
1.2
Properties of material
E = 117. E3 MPa
naked = 0.35
And
= 8.93 E-9 G/mm3
Y
Y = 100. MPa
AND = 400. MPa
Coefficient of friction of Coulomb:
E
µ = 0.25
1.3
Boundary conditions and loadings
The rigid foundation is completely blocked 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.
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
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:
93
.
3
+ 86
.
3
+ 72
.
3
+ 88
.
3
+ 96
.
3
· radial displacement of point a:
= 87
.
3
mm
5
-
24
.
13
-
63
.
13
-
62
.
13
-
57
.
13
-
24
.
13
· vertical displacement of point b:
= - 46
.
13
mm
5
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).
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
3 Modeling
With
3.1
Characteristics of modeling
3.2
Characteristics of the grid
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
Commands
AFFE_CHAR_MECA CONTACT
FROTTEMENT
“COULOMB”
METHODE
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
TOUT
“OUI”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR
RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
VITE
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
3.4
Sizes tested and results
Being given the heaviness of modeling, one takes only some steps of time and one carries out
only tests of not-regression.
3.5 Parameters
of execution
Version: 6.3.07
Machine: CLASTER
System: IRIX64
Obstruction memory:
300 megawords
Time CPU To use: 720 seconds
4 Modeling
B
4.1
Characteristics of modeling
4.2
Characteristics of the grid
A number of nodes: 359
Numbers and types of elements: 101 QUAD8, 55 SEG3
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
4.3 Functionalities
tested
Commands
AFFE_CHAR_MECA CONTACT
FROTTEMENT
“COULOMB”
METHODE
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
TOUT
“OUI”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
VITE
4.4
Sizes tested and results
Size tested
Reference:
Code_Aster
Difference (%)
radial displacement of
3.87 2.92
24.55
not A
vertical displacement of
13.46 12.69
5.72
not B
4.5 Parameters
of execution
Version: 6.3.07
Machine: CLASTER
System: IRIX64
Obstruction memory:
16 megawords
Time CPU To use: 434. seconds
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
5 Modeling
C
5.1
Characteristics of modeling
5.2
Characteristics of the grid
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
Commands
AFFE_CHAR_MECA CONTACT
FRICTION “COULOMB”
METHODE
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
TOUT
“OUI”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR
RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
VITE
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
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
Obstruction memory:
300 megawords
Time CPU To use: 185 seconds
6 Modeling
D
6.1
Characteristics of modeling
6.2
Characteristics of the grid
A number of nodes: 29
Numbers and types of elements: 5 QUAD8, 12 SEG3
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
6.3 Functionalities
tested
Commands
AFFE_CHAR_MECA CONTACT
FROTTEMENT
“COULOMB”
METHODE
“PENALIZATION”
CREA_CHAMP OPERATION
“AFFE”
TYPE_CHAM
“NOEU_DEPL_R”
AFFE
TOUT
“OUI”
NOM_CMP
(“DX”, “DY”, “DZ”)
DYNA_NON_LINE COMP_INCR RELATION
“VMIS_ISOT_LINE”
DEFORMATION
“SIMO_MIEHE”
INCREMENT
SUBD_PAS
4
ETAT_INIT
VITE
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
Obstruction memory:
16 megawords
Time CPU To use: 37. seconds
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
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
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.
Handbook of Validation
V5.03 booklet: Nonlinear dynamics of the voluminal structures
HT-66/02/001/A
Outline document