Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
1/8
Organization (S): EDF/EP/AMV
Handbook of Validation
V8.01 booklet: Fluid
V8.01.101 document
FDLV101 - Deux cylinders separated by a fluid
incompressible
Summary:
This test of the field of the fluids (fluid coupling/structure) validates the calculation of matrix of mass added in
case where one has several structures immersed in the same fluid.
By a modal analysis, one thus determines the coupled modes of the two structures because of the mass of
fluid which separates them. One adopts a modeling planes (thermal for the fluid, and deformation planes for
cylinders).
One finds the modes coupled of the system with less than 0.1% of the analytical result.
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
2/8
1
Problem of reference
1.1 Geometry
N
R2
M
R R
X
R1
()
1
()
2
Two cylinders separated by incompressible fluid:
interior radius R1 = 1.0 m external radius R2 = 1.1 m
1.2
Material properties
Fluid:
Water: O = 1000.0 Kg.m3
Solid:
Steel: S = 7800.0 Kg.m-3; E = 2.E11 Pa; = 0.3
Arises connecting the piston to the solid mass:
One places a discrete element on mesh POI1 in the center of the cylinder 1 of K1 stiffness and two
discrete elements on mesh POI1 on cylinder 2 on the level of the axis OX whose stiffness is worth
K2.
Discrete elements of type K_T_D_L:
K1= (1.E7, 1.E7, 1.E7) NR/m
K2 = (5.E6, 5.E6, 5.E6) NR/m
1.3
Boundary conditions and loading
One imposes a pressure (IE by analogy thermal a null temperature [R4.07.03]) in a node
unspecified of the fluid.
One imposes a null displacement of the cylinders according to OY.
1.4 Conditions
initial
Without object for the calculation of added mass and the modal analysis.
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
3/8
2
Reference solution
2.1
Method of calculation used for the reference solution
Analytical calculation:
One will suppose that the movements of the cylinders and the fluid are primarily plane. Effects
longitudinal will be neglected in front of the transverse effects. The problem is two-dimensional. Count
held of symmetry, the reference mark used is a cylindrical reference mark (R,) related to the central cylinder (see figure
above). In this frame of reference and with this particular geometry, the normal derivative
is equal to the derivative
compared to R.
N
R
In all this part, the variable p indicates the hydrodynamic field of pressure in the fluid
created by the natural vibrations from the structures, X1 or 2 indicates the clean modes of cylinder 1 or 2
respectively.
Clean modes of the hulls of border (
1) and (2) in the absence of fluid is form (N
indicate the command of the mode):
cos
N
or
sin N
0
1n
X
(R)
and
X2 N (R)
0
cos N
or
sin N
is the azimuth angle. These modes are uncoupled of course. The component 1ère corresponds to
normal displacement of the interior hull, 2nd with that of the external hull. In volume
fluid, one has thus two problems to solve:
cos N
p
p
1n
1n
p
= 0
= - or
= 0
éq 2.1-1
1n
F
N
N
sin N
1
2
and:
cos N
p
p
2n
2
p
= 0
= 0
N
= or
éq 2.1-2
2n
F
N
N
1
2
sin N
The field p N
1 corresponds to the field of pressure generated in the fluid if the central hull 1 vibrates
only, the field p2n is that created by the external hull 2 if it only vibrates. The linearity of
the equation of Laplace makes it possible to solve independently each problem and then of
to superimpose to find the field of pressure total.
The solution of the problem [éq 2.1-1] is, in polar co-ordinates, of the type [bib1]:
cos
N
N
1
p (R,)
= Arn +
B
or
1n
R
sin N
One must have N 0, bus if not one with the not-conservation of the volume of the fluid.
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
4/8
Constant A and B are determined by the boundary conditions:
cos
N
p
p
1n
1n
= - or
and
= 0
F
N
N
sin N
1
R
2
R
It is found whereas the field of pressure for each of the two problems is written:
N
N
N
N
R1
1
+
F
(R R) (R2 R) 1 (R R
2
) cos
p (R,)
=
or
2
1n
N
(
N
R2 R)
1
- 1
sin N
and:
N
N
N
N
R2
2
1
1
+
F
(R R) (R R) (R R
2
) cos
p
(R,)
=
or
2
2 N
N
(
N
R2 R)
1
- 1
sin N
With
The modal coefficients of added mass mijnm are calculated starting from the following formula
[R4.07.03] if I = 1 or 2, J = 1 or 2, (N, m) belongs to ¦2.
m With
=
p
X
im (R)
N
jn
(J) D
ijnm
J
J
The indexing is a little more complex here than in the formula presented in [R4.07.03]: indices I and
J refer to hulls 1 and 2, and indices m and N are associated the modes of hull. One
notice that there is coupling of the modes of the various hulls, external and intern.
It is noticed, on the one hand, that the fluid does not couple the modes of different indices N bus them
integrals
cos N cos m D
cancel themselves; in addition, the fluid does not couple either the modes
cos N and sin N bus
cos N sin N D =
0. The only existing coupling is a coupling between
()
two hulls for the modes of comparable nature.
With each mode N, one associates a matrix of command 4 symmetrical. A submatrix corresponding to
projection on mode N is written:
With
With
m
m
With
11nn
12nn
1
M
= A
With
m21nn m22nn
With
1
M
0
The total matrix is written:
With
With
with M = M
With
1
2
0
M
2
cos N
2
with
With
=
p
or
11
m N
L 1
R
,
1 (1
R
)
D
0
sin N
That is to say:
2n
R R
+ 1
With
2
(2) 1
=
11
m N
R L
éq 2.1-3
F
1
N
(R R) 2N -
2
1
1
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
5/8
one will obtain:
2n
R R
+ 1
With
2
(2) 1
m
=
22nn
R L
éq 2.1-4
F
2
N
(R R) 2N -
2
1
1
and:
N
R R
With
With
(2 2) 1
m
=
= -
21nn
12
m N
éq 2.1-5
F
1
R
2
R L
N
(R R) 2N -
2
1
1
L indicates here the height of the hulls cylinders in the longitudinal direction.
In our case, one considers only the modes of command N = 1 of the hulls: they correspond
respectively with the modes of translation of each hull along an axis passing by the center
central tube: one takes those arbitrarily corresponding to axis OX: coefficients of mass
added linear are written:
2
R R
+ 1
With
2 (2
) 1
=
11
m
R
F
1 (R R) 2 -
2
1
1
2
R R
+ 1
With
2 (2
) 1
m
=
22
R
F
2 (R R) 2 -
2
1
1
R R
With
With
(2 2) 1
m
=
= -
21
12
m
R R
F
1
2 (R R) 2 -
2
1
1
The equation of the generalized movement of the two coupled hulls is written:
m
0 X K
0 X
m
m X
1
1
1
1
11
12
1
+
= -
0 m X 0 K X
m
m X
2
2
2
2
12
22
2
The own pulsations of the coupled system are given by the equation of degree 4:
m + m
m
K
1
11
12
0
det
2
1
-
= 0
m
m + m
0 K
12
2
22
2
Numerical application:
K1 = 107 NR/m K2 = 107 NR/m
m11 = 33.060 kg/m
m22 = 40.004 kg/m
m12 = 36 200 kg/m
One obtains two Eigen frequencies:
f1 = 1.696 Hz f2 = 4.128 Hz
2.2
Results of reference
Analytical
2.3 References
bibliographical
[1]
R.J GIBERT. Vibrations of Structures. Interactions with fluids. Eyrolles (1988).
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
6/8
3 Modeling
With
3.1
Characteristics of modeling
Thermal formulation planes for fluid (QUAD4 and SEG2)
Plane and discrete deformation formulation for solid (TRIA3, QUAD4 and POI1)
This modeling is designed to determine the modes of command N = 1 of the cylinders. Modes of
hulls of a higher nature cannot be simulated by this type of model, but by one
modeling of the type COQUE_CYL [U4.22.01].
Cutting =
64 meshs QUAD4 on the circumference of the cylinders
64 meshs TRIA3 on the interior of the interior cylinder
64 meshs SEG2 on the fluid interface/cylinders
2 meshs QUAD4 following the thickness of the fluid
2 meshs QUAD4 following the thickness of the external cylinder
Boundary conditions:
DDL_IMPO: (Group_no: HANG DY: 0. )
DDL_IMPO: (Group_no: ACCREXT DY:0. )
TEMP_IMPO: (Group_no: TEMPIMPO TEMP: 0.)
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
7/8
3.2
Characteristics of the grid
A number of nodes: 356 QUAD4
A number of meshs and types: 64 TRIA3, 128 SEG2, 3 POI1
3.3 Functionalities
tested
Commands
Keys
AFFE_MODELE
“THERMIQUE”
“PLAN”
[U4.22.01]
CALC_MASS_AJOU
MODE_MECA
[U4.??.??]
NUME_DDL_GENE
NUME_DDL_GENE
“PLEIN”
STOCKAGE
[U4.55.07]
MODE_MECA
MODE_ITER_SIMULT
“BANDE”
FREQ
[U4.52.01]
concept
“matr_asse_gene_R”
COMB_MATR_ASSE
COMB_R
[U4.53.01]
concept
“matr_asse_gene_R”
4
Results of modeling A
4.1 Values
tested
Identification
Reference (Hz)
Aster (Hz)
% difference
Command of the clean mode I: 1
1.696
1.695808
0.011
Command of the clean mode I: 2
4.128
4.12473
0.079
4.2 Remarks
Calculations of modes carried out by:
MODE_ITER_SIMULT OPTION: “PLUS_PETITE” NMAX_FREQ: 2.
4.3 Parameters
of execution
Version: 3.05.24
Machine: CRAY C98
System:
UNICOS 8.0
Obstruction memory:
8 megawords
Time CPU To use:
10.28 seconds
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A
Code_Aster ®
Version
3
Titrate:
FDLV101 Deux cylinders separated by an incompressible fluid
Date:
24/08/99
Author (S):
G. ROUSSEAU
Key:
V8.01.101-A Page:
8/8
Intentionally white left page.
Handbook of Validation
V8.01 booklet: Fluid
HP-51/96/031 - Ind A