Code_Aster ®
Version
4.0
Titrate:
FDLV104 Masse added on model generalized 3D
Date:
12/01/98
Author (S):
G. ROUSSEAU
Key:
V8.01.104-A Page:
1/6
Organization (S): EDF/EP/AMV
Handbook of Validation
V8.01 booklet: Fluid
V8.01.104 document
FDLV104 - Calcul of mass added on model
generalized 3D
Summary:
This test of the field of the fluids (fluid coupling/structure) validates the calculation of mass added on a model
generalized of type 3D, namely two spheres coupled by an incompressible fluid.
By a modal analysis carried out by dynamic under-structuring, one obtains the coupled modes of both
spheres because of the presence of the fluid. The Eigen frequencies are compared with those determined by one
direct calculation.
One obtains the Eigen frequencies with a margin of less than 3%.
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
4.0
Titrate:
FDLV104 Masse added on model generalized 3D
Date:
12/01/98
Author (S):
G. ROUSSEAU
Key:
V8.01.104-A Page:
2/6
1
Problem of reference
1.1 Geometry
Sphere 1
Sphere 2
K1
K3
K2
K2
K3
K1
K3
K1
H
0
Y
K3
X
R
K1
K2
L
Z
L
Length:
L = 4 m
width:
L = 2 m
height:
H = 2 m
Thickness of the hulls:
E = 103 m
Radius of the spheres:
R = 0.5 m
1.2
Properties of materials
Structure: steel - material elastic
E = 2.1011 Pa
= 0.3
S = 7800 kg/m3
The springs have as respective stiffnesses:
K1 = 106 NR/m
K2 = 105 NR/m
K3 = 107 NR/m
Fluid: water - thermal material are equivalent
= 1.
the specific heat plays the part of density of the fluid: F = 1000 kg/m3
1.3
Boundary conditions and loadings
Side structure: the ddl of rotation DRX, DRY, DRZ of sphere 1 are blocked:
DRX: 0.0
DRY: 0.0
DRZ: 0.0
Fluid side: one imposes a pressure (i.e temperature) null in a node of the fluid grid.
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
4.0
Titrate:
FDLV104 Masse added on model generalized 3D
Date:
12/01/98
Author (S):
G. ROUSSEAU
Key:
V8.01.104-A Page:
3/6
2
Reference solution
2.1
Method of calculation used for the reference solution
One uses a direct modal calculation to obtain the Eigen frequencies modified by the fluid around
spheres.
The two spheres result by a translation from vector from components (2, 0., 0.) and for one
rotation of nautical analyzes (90., 0., 90.) what makes it possible to plan to build a model
generalized on the 1ère sphere. One will be able to carry out same modal calculation with coupling
fluid/structure by using the dynamic under-structuring. Direct calculation with taking into account of
mass added leads to the Eigen frequencies of the immersed system following:
f1 =
3.1172 Hz
f2 =
3.1183 Hz
f3 =
9.2727 Hz
f4 =
9.8267 Hz
f5 =
22.4400 Hz
f6 =
30.3295 Hz
2.2
Results of reference
Direct modal calculation by Code_Aster.
2.3
Uncertainty on the solution
Uncertainties on the solution are related to the discretization of the fluid interface/structure and on
the invariance of the grid by rotation of nautical angles (90°, 0°, 90°).
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
4.0
Titrate:
FDLV104 Masse added on model generalized 3D
Date:
12/01/98
Author (S):
G. ROUSSEAU
Key:
V8.01.104-A Page:
4/6
3 Modeling
With
3.1
Characteristics of modeling
Modeling includes/understands:
· side structure:
336 elements of hulls DKT modelling sphere 1,
6 discrete elements of type K_TR_L (modeling DIS_TR) modelling the six springs connecting
sphere 1 with the solid mass,
· fluid side:
672 THER_FACE3 thermal elements modelling the fluid interface/structure,
39.445 Thermal THER_TETRA4 elements modelling the fluid.
3.2
Characteristics of the grid
A number of meshs and types: 672 meshs TRIA3, 39.445 TETRA4, 6 SEG2
3.3 Functionalities
tested
Commands
Keys
CALC_MATR_AJOU
MODELE_GENE
[U4.55.10]
OPTION
“MASS_AJOU”
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
4.0
Titrate:
FDLV104 Masse added on model generalized 3D
Date:
12/01/98
Author (S):
G. ROUSSEAU
Key:
V8.01.104-A Page:
5/6
4
Results of modeling A
4.1 Values
tested
Identification
Hz
Hz
% difference
Reference
Aster
N°1 mode
3.1172
3.1189
+0.056%
N°2 mode
3.1183
3.1321
+0.44%
N°3 mode
9.2727
9.3144
+0.45%
N°4 mode
9.8267
9.8533
+0.27%
N°5 mode
22.4400
23.1291
+3.1%
N°6 mode
30.3295
31.1779
+2.8%
4.2 Remarks
The variation on the two last Eigen frequencies is explained by a not-symmetry of the grid of
the fluid interface/structure (thermal) compared to that of the first. Grid of the interface
fluid/structure (thermal elements) of the second sphere does not result exactly from that from
the first by rotation of nautical angles of 90°, 0°, 90°.
4.3 Parameters
of execution
Version: NEW 3.5.24
Machine: CRAY C90
Obstruction memory:
30 MW
Time CPU to use:
161.3 seconds
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
4.0
Titrate:
FDLV104 Masse added on model generalized 3D
Date:
12/01/98
Author (S):
G. ROUSSEAU
Key:
V8.01.104-A Page:
6/6
5
Summary of the results
The results obtained by calculation by dynamic under-structuring show perfect
agreement with direct calculation. Nevertheless, to have this agreement, it should be taken care that them
grids of the interfaces fluid/structure (thermal elements) of the substructures which are repeated with
the interior of the fluid grid result perfectly by rotations of nautical angle defined in
DEFI_MODELE_GENE.
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A