Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
1/20
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
Organization (S):
EDF-R & D/AMA, DeltaCAD
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
Document: V6.04.504
SSNV504 - Extrusion of a piece
Summary:
This test simulates the extrusion (without friction) of a piece, case frequently encountered in the industrial studies
of working. The interest of this test is mainly to validate the “passage” of geometrical singularities
(acute angles and blunt) and to test the taking into account of two areas of contact (contact multi-areas).
Modelings selected are as follows:
·
Modeling A (AXIS):
CONTACT node nets, associated with meshs SEG2,
·
Modeling B (3D):
CONTACT node nets, associated with meshs QUAD4,
·
Modeling C (AXIS): CONTACT node nets SLIDE, associated with meshs SEG2,
·
Modeling D (AXIS): CONTACT, method CONTINUES associated with meshs SEG2,
·
Modeling E (AXIS):
CONTACT node nets SLIDE, associated with meshs SEG3,
·
Modeling F (3D):
CONTACT node nets SLIDE, associated with meshs QUAD4,
·
Modeling G (3D):
CONTACT node nets SLIDE, associated with meshs QUAD8,
·
Modeling H (2D):
CONTACT, method CONTINUES associated with SEG3
Modelings with SLIDE make it possible to maintain the contact throughout extrusion, in way
mathematics. The result is physically close to modeling without slide, since it is about one
extrusion without friction.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
2/20
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
1
Problem of reference
1.1 Geometry
With
B
C
D
E
G
F
Die
Piece 50.
65.
R
45.
48.
70. 15.
30. 20.
70.
15.10.
H
I
J
K
L
NR
M
Z
Axis of symmetry
Dimension in mm
1.2
Properties of material
Piece:
·
E = 5000. MPa
Young modulus
·
= 0.45
Poisson's ratio
Die:
·
E = 200.000. MPa
Young modulus
·
= 0.3
Poisson's ratio
Piece/die
·
µ = 0
Coefficient of friction
1.3
Boundary conditions and loadings
- C.L. : lines HI, IJ and JA embedded
- Loading: Piloting in displacement imposed on the rear face of the piece
1.4 Conditions
initial
Without object.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
3/20
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
2
Reference solution
2.1
Method of calculation used for the reference solution
The objective of this case-test is to analyze the feasibility of calculation into axisymmetric and 3D.
2.2
Results of reference
No value of reference is available. The values of reference which will be retained to test
the future versions of Code_Aster are those obtained during the first execution with
modeling A and which will be considered to be acceptable. The analyzed values of reference are:
·
displacements of the front face of the piece,
·
the maximum iteration count to converge.
Piece
(localization)
Déplacemen
T (mm)
Comments
K
5.000
Piece in the vicinity of the point B of the die
K
20.825
Piece in the vicinity of the point C of the die
K
55.880
Piece in the vicinity of the point D of the die
K
78.690
Piece in the vicinity of the point E of the die
K
144.895
Piece in the vicinity of the point F of the die
K
155.096
Piece at the point G of the die
2.3
Uncertainties on the solution
< 0.1%
2.4 References
bibliographical
None
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
4/20
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
X
7
3
4
3
5
5
5
4
Solid: Modeling AXIS (TRIA3)
Contact: CONTACT (SEG2)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Area of contact n°2: LK/FG
Surface slave
Surface main
Y
M
NR
Limiting conditions:
- dimensioned JA, IJ, HI: DX=0., DY=0.
- dimensioned MN: DY=1.
H
I
J
3.2
Characteristics of the mesh
A number of nodes: 333
A number of meshs: 510 TRIA3 and 152 SEG2
3.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_2D
AFFE_MODELE
AFFE
MODELING = “AXIS”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
NOEUD_BOUCLE
METHOD = “FORCED”
FRICTION = “WITHOUT”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
5/20
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
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.000 0.
·
Numbers iterations
1
1
0.
Piece at the point C
·
Displacement at the point K (mm) 20.825
20.825 0.
·
Numbers iterations
1
1
0.
Piece at the point D
·
Displacement at the point K (mm) 55.880
55.880 0.
·
Numbers iterations
1
1
0.
Piece at the point E
·
Displacement at the point K (mm) 78.690
78.690 0.
·
Numbers iterations
1
1
0.
Piece at the point F
·
Displacement at the point K (mm) 144.895
144.895 0.
·
Numbers iterations
1
1
0.
Piece at the point G
·
Displacement at the point K (mm) 155.096
155.096 0.
·
Numbers iterations
1
1
0.
4.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 5 pitches
-
of
70.mm with
140.mm in 5 pitches
-
of
140.mm with
155.mm in 5 pitches
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
6/20
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
X
7
3
4
3
5
5
5
4
Solid: Modeling 3D (HEXA8, PENTA6)
Contact: CONTACT (QUAD4)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Area of contact n°2: LK/FG
Surface slave
Surface main
Y
Discretization: mesh ¼ of the structure
Z
M
NR
Boundary conditions on the groups of following nodes:
- “FIL_EXT”: group nodes located on surface external of die (HI, IJ, JA)
=> DX=0., DY =0., DZ=0.
- “FIL_SYM1”: Group nodes of the die located in plan XOY: DZ=0.
- “FIL_SYM2”: Group nodes of the die located in plan YOZ: DX=0.
- “LOP_SYM1”: Group nodes of the piece located in plan XOY: DZ=0.
- “LOP_SYM2”: Group nodes of the piece located in plan YOZ: DX=0
- “LOP_DDL”: Group nodes located on the rear face of the piece: DY =1
H
I
J
Z
X
5.2
Characteristics of the mesh
A number of nodes: 3292
A number of meshs: 2150 HEXA8, 260 PENTA6, 1814 QUAD4 and 68 TRIA3
5.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_3D
AFFE_MODELE
AFFE
MODELING = “3D”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
METHOD = “FORCED”
FRICTION = “WITHOUT”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
7/20
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
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.000
0.
·
Numbers iterations
1
1
0.
Piece at the point C
·
Displacement at the point K (mm) 20.825
20.796
0.14
·
Numbers iterations
1
1
0.
Piece at the point D
·
Displacement at the point K (mm) 55.880
55.642
0.425
·
Numbers iterations
1
1
0.
Piece at the point E
·
Displacement at the point K (mm) 78.690
78.627
0.081
·
Numbers iterations
1
1
0.
Piece at the point F
·
Displacement at the point K (mm) 144.895
144.983 0.061
·
Numbers iterations
1
1
0.
Piece at the point G
·
Displacement at the point K (mm) 155.096
155.065 0.02
·
Numbers iterations
1
1
0.
6.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
15.mm in 10 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 10 pitches
-
of
70.mm with
140.mm in 35 pitches
-
of
140. mm with
155.mm in 15 pitch
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
8/20
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
X
7
3
4
3
5
5
5
4
Solid: Modeling AXIS (TRIA3)
Contact: CONTACT (SEG2)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Surface slave
Surface main
Y
M
NR
Limiting conditions:
- dimensioned JA, IJ, HI: DX=0., DY=0.
- dimensioned MN: DY=1.
H
I
J
One uses here the function SLIDE which makes it possible to maintain the contact throughout extrusion, of
mathematical way. The result is physically close to modeling without slide, since it
acts of an extrusion without friction.
7.2
Characteristics of the mesh
A number of nodes: 333
A number of meshs: 510 TRIA3 and 152 SEG2
7.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_2D
AFFE_MODELE
AFFE
MODELING = “AXIS”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
NOEUD_BOUCLE
METHOD = “FORCED”
FRICTION = “WITHOUT”
SLIDE = “YES”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
9/20
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
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.000 0.
·
Numbers iterations
1
1
-
Piece at the point C
·
Displacement at the point K (mm) 20.825
20.585
- 1.152
·
Numbers iterations
1
1
-
Piece at the point D
·
Displacement at the point K (mm) 55.880
55.541
- 0.847
·
Numbers iterations
1
3
-
Piece at the point E
·
Displacement at the point K (mm) 78.690
78.372
- 0.404
·
Numbers iterations
1
3
-
Piece at the point F
·
Displacement at the point K (mm) 144.895
144.224 - 0.463
·
Numbers iterations
1
3
-
8.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 5 pitches
-
of
70.mm with
140.mm in 5 pitches
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
10/20
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
X
7
3
4
3
5
5
5
4
Solid: Modeling AXIS (TRIA3)
Contact: CONTACT (SEG2)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Area of contact n°2: LK/FG
Surface slave
Surface main
Y
M
NR
Limiting conditions:
- dimensioned JA, IJ, HI: DX=0., DY=0.
- dimensioned MN: DY=1.
H
I
J
9.2
Characteristics of the mesh
A number of nodes: 333
A number of meshs: 510 TRIA3 and 152 SEG2
9.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_2D
AFFE_MODELE
AFFE
MODELING = “AXIS”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
NOEUD_BOUCLE
METHOD = “CONTINUES”
FRICTION = “WITHOUT”
SLIDE = “NOT”
ITER_GEOM_MAXI = 1
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
11/20
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
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.000 0.
·
Numbers iterations
1
1
0.
Piece at the point C
·
Displacement at the point K (mm) 20.825
20.790
- 0.167
·
Numbers iterations
1
1
0.
Piece at the point D
·
Displacement at the point K (mm) 55.880
55.831
- 0.088
·
Numbers iterations
1
1
0.
Piece at the point E
·
Displacement at the point K (mm) 78.690
78.684
- 0.008
·
Numbers iterations
1
1
0.
Piece at the point F
·
Displacement at the point K (mm) 144.895
144.893 - 0.002
·
Numbers iterations
1
1
0.
Piece at the point G
·
Displacement at the point K (mm) 155.096
155.071 0.016
·
Numbers iterations
1
1
0.
10.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 5 pitches
-
of
70.mm with
140.mm in 5 pitches
-
of
140.mm with
155.mm in 5 pitches
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
12/20
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
X
7
3
4
3
5
5
5
4
Solid: Modeling AXIS (TRIA6)
Contact: CONTACT (SEG3)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Surface slave
Surface main
Y
M
NR
Limiting conditions:
- dimensioned JA, IJ, HI: DX=0., DY=0.
- dimensioned MN: DY=1.
H
I
J
11.2 Characteristics of the mesh
A number of nodes: 333
A number of meshs: 510 TRIA6 and 152 SEG3
11.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_2D
AFFE_MODELE
AFFE
MODELING = “AXIS”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
NOEUD_BOUCLE
METHOD = “FORCED”
FRICTION = “WITHOUT”
SLIDE = “YES”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
13/20
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
One uses here the function SLIDE which makes it possible to maintain the contact throughout extrusion, of
mathematical way. The result is physically close to modeling without slide, since it
acts of an extrusion without friction. The values of reference are regarded as being those
modeling C.
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.000 0.
·
Numbers iterations
1
1
0.
Piece at the point C
·
Displacement at the point K (mm) 20.585
20.612
- 0.001
·
Numbers iterations
1
2
-
Piece at the point D
·
Displacement at the point K (mm) 55.407
55.412 0.
·
Numbers iterations
3
2
-
Piece at the point E
·
Displacement at the point K (mm) 78.372
78.379 0.
·
Numbers iterations
3
2
-
Piece at the point F
·
Displacement at the point K (mm) 144.224
144.224 0.
·
Numbers iterations
3
2
-
12.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 5 pitches
-
of
70.mm with
140.mm in 5 pitches
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
14/20
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
X
7
3
4
3
5
5
5
4
Solid: Modeling 3D (HEXA8, PENTA6)
Contact: CONTACT (QUAD4)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Surface slave
Surface main
Y
Discretization: mesh ¼ of the structure
Z
M
NR
Boundary conditions on the groups of following nodes:
- “FIL_EXT”: group nodes located on surface external of die (HI, IJ, JA)
=> DX=0., DY =0., DZ=0.
- “FIL_SYM1”: Group nodes of the die located in plan XOY: DZ=0.
- “FIL_SYM2”: Group nodes of the die located in plan YOZ: DX=0.
- “LOP_SYM1”: Group nodes of the piece located in plan XOY: DZ=0.
- “LOP_SYM2”: Group nodes of the piece located in plan YOZ: DX=0
- “LOP_DDL”: Group nodes located on the rear face of the piece: DY =1
H
I
J
Z
X
13.2 Characteristics of the mesh
A number of nodes: 3292
A number of meshs: 2150 HEXA8, 260 PENTA6, 1814 QUAD4 and 68 TRIA3
13.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_3D
AFFE_MODELE
AFFE
MODELING = “3D”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
METHOD = “FORCED”
FRICTION = “WITHOUT”
SLIDE = “YES”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
15/20
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
One uses here the function SLIDE which makes it possible to maintain the contact throughout extrusion, of
mathematical way. The result is physically close to modeling without slide, since it
acts of an extrusion without friction. The values of reference are regarded as being those
modeling C.
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.000
0.
·
Numbers iterations
1
1
0.
Piece at the point C
·
Displacement at the point K (mm) 20.585
20.612
0.001
·
Numbers iterations
1
2
-
Piece at the point D
·
Displacement at the point K (mm) 55.407
55.412 0.
·
Numbers iterations
3
2
-
Piece at the point E
·
Displacement at the point K (mm) 78.372
78.379 0.
·
Numbers iterations
3
2
-
Piece at the point F
·
Displacement at the point K (mm) 144.224
144.224 0.
·
Numbers iterations
3
2
-
14.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
15.mm in 10 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 10 pitches
-
of
70.mm with
140.mm in 35 pitches
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
16/20
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
X
7
3
4
3
5
5
5
4
Solid: Modeling 3D (HEXA20, PENTA15)
Contact: CONTACT (QUAD8)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Surface slave
Surface main
Y
Discretization: mesh ¼ of the structure
Z
M
NR
Boundary conditions on the groups of following nodes:
- “FIL_EXT”: group nodes located on surface external of die (HI, IJ, JA)
=> DX=0., DY =0., DZ=0.
- “FIL_SYM1”: Group nodes of the die located in plan XOY: DZ=0.
- “FIL_SYM2”: Group nodes of the die located in plan YOZ: DX=0.
- “LOP_SYM1”: Group nodes of the piece located in plan XOY: DZ=0.
- “LOP_SYM2”: Group nodes of the piece located in plan YOZ: DX=0
- “LOP_DDL”: Group nodes located on the rear face of the piece: DY =1
H
I
J
Z
X
15.2 Characteristics of the mesh
A number of nodes: 12213
A number of meshs: 2150 HEXA20, 260 PENTA15, 1814 QUAD8 and 68 TRIA6
15.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_3D
AFFE_MODELE
AFFE
MODELING = “3D”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
METHOD = “FORCED”
FRICTION = “WITHOUT”
SLIDE = “YES”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
17/20
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
One uses here the function SLIDE which makes it possible to maintain the contact throughout extrusion, of
mathematical way. The result is physically close to modeling without slide, since it
acts of an extrusion without friction. The values of reference are regarded as being those
modeling C.
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.221
0.044
·
Numbers iterations
1
1
-
Piece at the point C
·
Displacement at the point K (mm) 20.585
20.822
0.012
·
Numbers iterations
1
1
-
Piece at the point D
·
Displacement at the point K (mm) 55.407
55.532
0.002
·
Numbers iterations
3
1
-
Piece at the point E
·
Displacement at the point K (mm) 78.372
78.476
0.001
·
Numbers iterations
3
1
-
Piece at the point F
·
Displacement at the point K (mm) 144.224
144.364 0.001
·
Numbers iterations
3
1
-
16.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
15.mm in 10 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 10 pitches
-
of
70.mm with
140.mm in 35 pitches
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
18/20
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 with quadratic elements of contact (SEG3) by the method CONTINUES. It
case-test presents phenomena of flip-flop (oscillations of the statute of the contact related to difficulties
of zero numerical), solved here by an internal test on the maximum number of flip-flop (transparent with
the user, and fixed at 15), i.e. the maximum number of changes of statute of the contact. One
message of alarm indicates to the user that the maximum number of changes of statute (15) was
reached.
X
7
3
4
3
5
5
5
4
Solid: Modeling AXIS (TRIA6)
Contact: CONTACT (SEG3)
With
B
C
D
F
E
K
G
L
Area of contact n°1: NL/ABCDEF
Surface slave
Surface main
Y
M
NR
Limiting conditions:
- dimensioned JA, IJ, HI: DX=0., DY=0.
- dimensioned MN: DY=1.
H
I
J
17.2 Characteristics of the mesh
A number of nodes: 1174
A number of meshs: 510 TRIA6 and 152 SEG3
17.3 Functionalities
tested
Controls Key word
factor
Key word
MODI_MAILLAGE ORIE_PEAU_2D
AFFE_MODELE
AFFE
MODELING = “AXIS”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
SEEK = “NOEUD_VOISIN”
METHOD = “CONTINUES”
FRICTION = “WITHOUT”
SLIDE = “NOT”
ITER_GEOM_MAXI = 3
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
19/20
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
Identification
Reference
Aster %
difference
Piece at the point B
·
Displacement at the point K (mm) 5.000
5.000 0.
·
Numbers iterations
1
0
100.
Piece at the point C
·
Displacement at the point K (mm) 20.825
20.801
- 0.115
·
Numbers iterations
1
1
0.
Piece at the point D
·
Displacement at the point K (mm) 55.880
55.678
- 0.361
·
Numbers iterations
1
1
0.
Piece at the point E
·
Displacement at the point K (mm) 78.690
78.691
- 0.002
·
Numbers iterations
1
1
0.
Piece at the point F
·
Displacement at the point K (mm) 144.895
144.949 - 0.037
·
Numbers iterations
1
1
0.
Piece at the point G
·
Displacement at the point K (mm) 155.096
154.674 0.272
·
Numbers iterations
1
1
0.
18.2 Remarks
·
Calculation is carried out by imposing a displacement on the rear face of the piece (MN).
displacement is imposed in the following way:
-
of
0.mm with
5.mm in 5 pitches
-
of
5.mm with
20.mm in 5 pitches
-
of
20.mm with
50.mm in 5 pitches
-
of
50.mm with
70.mm in 5 pitches
-
of
70.mm with
140.mm in 5 pitches
-
of
140.mm with
155.mm in 5 pitches
·
Although the data file passes without problem in linear statics (MECA_STATIQUE), them
calculations do not converge with the key words by defect used for convergence in
STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E
- 6
is too constraining
(forces to which the piece is subjected being at the beginning relatively weak). To mitigate it
problem, it is necessary to use the key word RESI_GLOB_MAXI = 1.E
- 6
.
Code_Aster
®
Version
8.3
Titrate:
SSNV504 - Extrusion of a piece
Date:
04/05/06
Author (S):
Mr. KHAM, P. MASSIN
Key
:
V6.04.504-B
Page:
20/20
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
19 Summary of the results
Two modelings (Axisymmetric and 3D) pass without problem. Convergence is fast.
The results between modeling AXIS and modeling 3D are similar.
In the same way, the methods of the STRESSES and CONTINUOUS also give results
satisfactory. On this example, the method CONTINUES seems to converge more quickly than that
STRESSES, in particular in 3D with SLIDE.
The particular case of the contact slide is attached to this case-test: modeling C is then taken
like reference. By the method CONTINUES, one obtains solutions close to these values of
reference, with generally a lower iteration count.
On the figure above we present the position of the piece for a displacement of the face postpones
piece of 20, 70, 155 mm
20mm 70mm 155mm
This test made it possible to validate:
·
the “passage” of geometrical singularities (acute angles and blunt),
·
the taking into account of two areas of contact (contact multi-areas),
·
the modeling of a contact-slide.