Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
1/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
Organization (S):
EDF-R & D/AMA, SAMTECH
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
Document: V6.04.140
SSNV140 - Embedded cylindrical panel
Summary:
One presents in this test a quasi static calculation of embedded elastic cylindrical panel subjected to one
surface force data either in the total reference mark, or in the local reference mark. This force is constant in
first case and following in the second case. Modeling thus is tested
COQUE_3D
nonlinear
geometrical by using the algorithm of update of great rotations 3D
GREEN_GR
of
STAT_NON_LINE
,
as well as the processing of the following pressures for modeling
COQUE_3D
. The reference solution which
is numerical is obtained with software the SAMCEF software.
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
2/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
1
Problem of reference
1.1 Geometry
The lengths are expressed in meters.
X
y
Z
125
.
0
=
H
10
2 =
L
100
=
R
1
P
2
P
3
P
4
P
rad
2
.
0
2
=
Z
F
The cylindrical panel is embedded along its 4 sides and subjected to a surface force.
·
data in the total reference mark for modeling a:
F
E
= -
>
F
F
Z
Z
Z
;
0
·
data in the local reference mark for modeling b:
F
N
= -
>
F
F
Z
Z
;
0
who leads to the membrane compression of the panel, accompanied by localized bendings. In
case of modeling B, N is the normal with the reactualized geometry of the panel.
Because of the geometrical and physical symmetry of the problem, only the quarter
P P P P
1 2 3 4
panel
is modelized, by taking account of the conditions of symmetry.
1.2
Material properties
Elastic behavior:
E
=
=
450000
0 3
AP;
.
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
3/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
1.3
Boundary conditions and loadings
C.L. :
P2P3: DX=0.
DY=0.
DZ=0.
DRX=0
DRY=0
DRZ=0
P3P4: DX=0.
DY=0.
DZ=0.
DRX=0
DRY=0
DRZ=0
Symmetry: P1P2:
DY=0.
DRX=0.
DRZ=0.
P4P1: DX=0.
DRY=0.
DRZ=0.
One seeks the successive states of balance under the loading made up of the surface force
data in the total reference mark:
()
F T
T
Z
=
T
being pseudo time.
One is only interested in the component vertical of displacement in
P
1
.
2
Reference solution
2.1
Method of calculation used for the reference solution
This solution [bib4] is that which is obtained with software the SAMCEF software [bib1]. Modeling is
based on a theory of hull in resulting efforts with a Co-rotational formulation [bib3] and
a discretization DSQ [bib2].
elements
16
elements
16
The mesh considered in the calculation of reference is a regular mesh of
16 16
×
elements
quadrilateral with 4 nodes each one.
2.2
Results of reference
History of vertical displacement DZ in meters with the node charged P1
Moment
Surface force
fz (NR)
DZ in P1
(m)
0.2 0.2
1.043E01
0.4 0.4
3.935E01
0.7 0.7
5.13E01
2.3 References
bibliographical
[1]
The SAMCEF software, Manual of reference V7.1 Elements Volume, 1998
[2]
J-L. Batoz, G.Dhatt, “Modeling of the Structures by Finite elements:Beams and Plates”,
Hermès, Paris, 1992
[3]
Crisfield Mr. A., “Non-linear Finite Element Analysis off Solids and Structures”, Volume 1:
Essentials, John Wiley, Chichester, 1994
[4]
PH. JETTEUR, Nonlinear Kinematics of the Hulls. Report/ratio SAMTECH, Contract
PP/GC-134/96, 1998
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
4/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
3 Modeling
With
3.1
Characteristics of modeling
Modeling
COQUE_3D
Element MEC3QU9H (voluminal hull 3D) - regular Mesh
elements
8
elements
8
3.2
Characteristics of the mesh
A number of nodes: 320
A number of meshs and type: 64 QUAD9
The mesh is twice less fine than that of the reference solution.
3.3 Functionalities
tested
·
Modeling COQUE_3D into nonlinear geometrical.
·
The static algorithm of update of great rotations
GREEN_GR
of
STAT_NON_LINE
.
Order Mots
keys
STAT_NON_LINE
COMP_ELAS
RELATION: “elas”
COQUE_NCOU: 1
DEFORMATION: “green_gr”
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
5/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
4
Results of modeling A
4.1 Values
tested
History of vertical displacement DZ to the node charged
Moment
Force
surface fz
Aster
Reference
% difference
0.2 0.2
1.040E01
1.043E01
0.351
0.4 0.4
3.920E01
3.935E01
0.362
4.2 Remarks
The default value of COEF_RIGI_DRZ = 0.00001 was selected (value generally used in
analyze linear).
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
6/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
5 Modeling
B
It is identical to modeling A, except for loading which is following: it is reactualized
according to the geometry. One goes until moment 0.7 per pitch of time of 0.1. There are thus 7 pitches of
charge.
5.1
Characteristics of modeling
Modeling
COQUE_3D
Element MEC3QU9H (voluminal hull 3D) - regular Mesh
elements
8
elements
8
5.2
Characteristics of the mesh
A number of nodes: 320
A number of meshs and type: 64 QUAD9
The mesh is twice less fine than that of the reference solution.
5.3 Functionalities
tested
·
Modeling COQUE_3D into nonlinear geometrical.
·
The static algorithm of update of great rotations
GREEN_GR
of
STAT_NON_LINE
.
The following loading is obtained in the following way:
Order Mots
keys
STAT_NON_LINE
COMP_ELAS
RELATION: “elas”
COQUE_NCOU: 1
DEFORMATION: “green_gr”
EXCIT
TYPE_CHARGE: SUIV
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
7/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
6
Results of modeling B
6.1 Values
tested
The values of reference come from software the SAMCEF software. They are given as an indication.
History of vertical displacement DZ to the node charged
Moment
Force
surface fz
Aster
Reference
(The SAMCEF software)
% difference
0.4 0.4
3.92E01
3.93E01
0.25%
0.7 0.7
5.11E01
5.13E01
0.28%
6.2 Remarks
Rotations not being very large, it is enough to use the default value of COEF_RIGI_DRZ =
0.00001 (value generally used in linear analysis).
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
8/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
7 Modeling
C
It is identical to modeling A, except for the functionality tested, COMP_INCR instead of
COMP_ELAS. As in both cases, one remains in the elastic range, one must thus find them
same results.
7.1
Characteristics of modeling
Modeling
COQUE_3D
Element MEC3QU9H (voluminal hull 3D) - regular Mesh
elements
8
elements
8
7.2
Characteristics of the mesh
A number of nodes: 320
A number of meshs and type: 64 QUAD9
The mesh is twice less fine than that of the reference solution.
7.3 Functionalities
tested
·
Modeling COQUE_3D into nonlinear geometrical.
·
The static algorithm of update of great rotations
GREEN_GR
of
STAT_NON_LINE
.
The following loading is obtained in the following way:
Order Mots
keys
STAT_NON_LINE
COMP_INCR
RELATION: “elas”
COQUE_NCOU: 1
DEFORMATION: “green_gr”
EXCIT
TYPE_CHARGE: SUIV
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
9/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
8
Results of modeling C
8.1 Values
tested
History of vertical displacement DZ to the node charged
Moment
Force
surface fz
Aster
Reference
% difference
0.2 0.2
1.040E01
1.043E01
0.352
0.4 0.4
3.921E01
3.935E01
0.360
8.2 Remarks
The default value of COEF_RIGI_DRZ = 0.00001 was selected (value generally used in
analyze linear).
Code_Aster
®
Version
7.4
Titrate:
SSNV140 - Embedded cylindrical panel
Date:
02/11/05
Author (S):
P. MASSIN, Mr. Al MIKDAD
Key
:
V6.04.140-B
Page:
10/10
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-66/05/005/A
9
Summary of the results
The facts of the case correspond to a thin hull
H
L = 0 625%
.
what is severe for
the finite element triangle MEC3TR7H (not presented here bus presenting a case of blocking at shearing
transverse). The deformation is primarily membranous and rotations remain moderate.
The reference solution obtained by software the SAMCEF software is numerical and its mesh twice more
end that that used here. Its model finite elements rests on an approach in resulting efforts with
a rotational formulation [bib3]. The approach chosen in Code_Aster is 3D forced plane
with a formulation Lagrangian total [R3.07.04]. The quality of the results obtained is good. The variation
compared to this solution is lower than 0.5% for the vertical displacement of the central point.