Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
1/8
Organization (S): EDF/IMA/MN, SAMTECH
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
Document: V3.04.104
SSLV104 - Poutre in rotation
Summary:
This test makes it possible to validate the linear elastic design of a slim beam subjected to a rotation of one of its
ends. Three modelings are tested: elements 3D (HEXA20) and elements of COQUE_3D (QUAD9 and
TRIA7). That tests the inertias of rotation, without taking account of the elementary terms of stiffening
centrifuge (cf [V3.04.105]).
The reference solution is analytical (1D). The results coincide perfectly with the reference solution.
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A
Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
2/8
1
Problem of reference
1.1 Geometry
Directed slim beam carried in space by the axis of directing vector (1, 1, 1).
0.5 m
Z
0.02 m
Y
Embedded face
X
Square section of surface: 4.0 104 m2
Length of the beam: 0.5 m
1.2
Material properties
E = 2. 1011 Pa
= 0
= 7800 kg/m3
A_CIS = 0.8333 (factor of correction of transverse shearing equal to 5/6 for a theory of the type
Thin Reissner hull)
1.3
Boundary conditions and loadings
Free beam fixed in rotation around an axis perpendicular to its greater dimension and
passing by the center of the embedded face.
Component of the vector rotation: (1, 0, 1).
Disk speed: = 3000 rd/S.
The important value disk speed does not have anything physics.
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A
Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
3/8
2
Reference solution
2.1
Method of calculation used for the reference solution
In the local reference mark of the beam:
2 Ux
+ 2x = 0 with U () =
0
0
2
X
X
E
Ux (L) = (L) =
0
X
xx
By integrating the preceding differential equation one obtains, in the reference mark of the beam:
2
x3
U (X) =
X L2 -
U = U2 =
X
0
2nd
y
3
Displacements of all points of the beam are thus written in the total reference mark:
2
R 3
U
=
2 -
X (X, Y, Z)
R L
2 3rd
3
2
R 3
U
=
2 -
y (X, Y, Z)
R L
2 3rd
3
2
R 3
U
=
2 -
Z (X, Y, Z)
R L
2 3rd
3
front
R
EC. =
X 2 + Y 2 + Z 2
2.2
Results of reference
Values of three displacements in the center of the section furthest away from the axis of rotation.
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A
Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
4/8
3 Modeling
With
3.1
Characteristics of modeling
Elements 3D (HEXA20)
Regulated grid including:
4 elements on the section
50 elements over the length
200 elements on the whole
0.5 m
Z
0.02 m
Y
Embedded face
X
3.2
Characteristics of the grid
A number of nodes: 1521
A number of meshs and types: 200 HEXA20
3.3 Functionalities
tested
Commands
Keys
AFFE_MODELE
AFFE
PHENOMENE
“MECANIQUE”
[U4.22.01]
MODELISATION “3D”
AFFE_CHAR_MECA
DDL_IMPO
DX, DY, DZ
[U4.25.01]
ROTATION
CALC_MATR_ELEM
OPTION
“RIGI_MECA”
[U4.41.01]
CALC_VECT_ELEM
OPTION
“CHAR_MECA”
[U4.41.02]
RESO_LDLT
MATR_FACT
[U4.51.02]
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A
Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
5/8
4
Results of modeling A
4.1 Values
tested
Identification
Reference
Aster
% difference
DX in L
8.44 103
8.44 103
0.05
DY in L
8.44 103
8.44 103
0.04
DZ in L
8.44 103
8.44 103
0.04
4.2 Parameters
of execution
Version: 4.00.02
Machine: CRAY C90
Obstruction memory:
8 MW
Time CPU To use:
6 seconds
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A
Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
6/8
5 Modeling
B
5.1
Characteristics of modeling
Elements of hull MEC3QU9H
0.02
Grid including:
2 elements according to the dispatcher
4 elements according to the length
8
elements
with
total
0.5 m
Z
y
X
Embedded face
5.2
Characteristics of the grid
A number of nodes: 45
A number of meshs and types: 8 QUAD9
5.3 Functionalities
tested
Commands
Keys
AFFE_MODELE
AFFE
PHENOMENE
“MECANIQUE”
[U4.22.01]
MODELISATION “3D”
AFFE_CHAR_MECA
DDL_IMPO
DX, DY, DZ
[U4.25.01]
ROTATION
CALC_MATR_ELEM
OPTION
“RIGI_MECA”
[U4.41.01]
CALC_VECT_ELEM
OPTION
“CHAR_MECA”
[U4.41.02]
RESO_LDLT
MATR_FACT
[U4.51.02]
6
Results of modeling B
6.1 Values
tested
Identification
Reference
Aster
% difference
DX in L
8.44 103
8.44 E3
0.04
DY in L
8.44 103
8.44 E3
0.04
DZ in L
8.44 103
8.44 E3
0.04
6.2 Parameters
of execution
Version: 4.00.14
Machine: CRAY C90
Obstruction memory:
16 MW
Time CPU To use:
4.3 seconds
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A
Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
7/8
7 Modeling
C
7.1
Characteristics of modeling
Elements of hull MEC3TR7H
0.02
Grid including:
4 elements according to the dispatcher
8 elements according to the length
64
elements
with
total
0.5 m
Z
y
X
Embedded face
7.2
Characteristics of the grid
A number of nodes: 217
A number of meshs and types: 64 TRIA7
7.3 Functionalities
tested
Commands
Keys
AFFE_MODELE
AFFE
PHENOMENE
“MECANIQUE”
[U4.22.01]
MODELISATION “3D”
AFFE_CHAR_MECA
DDL_IMPO
DX, DY, DZ
[U4.25.01]
ROTATION
CALC_MATR_ELEM
OPTION
“RIGI_MECA”
[U4.41.01]
CALC_VECT_ELEM
OPTION
“CHAR_MECA”
[U4.41.02]
RESO_LDLT
MATR_FACT
[U4.51.02]
8
Results of modeling C
8.1 Values
tested
Identification
Reference
Aster
% difference
DX in L
8.44 103
8.44 E3
0.08
DY in L
8.44 103
8.44 E3
0.02
DZ in L
8.44 103
8.44 E3
0.03
8.2 Parameters
of execution
Version: 4.00.14
Machine: CRAY C90
Obstruction memory:
16 MW
Time CPU To use:
7.6 seconds
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A
Code_Aster ®
Version
4.0
Titrate:
SSLV104 Poutre in rotation
Date:
07/12/98
Author (S):
J.M. PROIX, P. MASSIN, A. LAULUSA
Key:
V3.04.104-B Page:
8/8
9
Summary of the results
The coincidence of the results with the analytical solution makes it possible to validate the loading due to the forces
of inertia of rotation.
Modeling COQUE_3D with MEC3QU9H gives the solution with very few elements.
One will refer to test SSLV105 [V3.04.105] to evaluate the effect of the stiffening centrifuges for
the element 3D, HEXA20.
Handbook of Validation
V3.04 booklet: Linear statics of the voluminal structures
HI-75/98/040 - Ind A