COHF-8406 47-9 DE84 009185: H. Chung and S. S. Chen
COHF-8406 47-9 DE84 009185: H. Chung and S. S. Chen
COHF-8406 47-9 DE84 009185: H. Chung and S. S. Chen
COHF-8406 47—9
DE84 009185
HYDRODYNAMIC MASS
",. « a'S?
by
DISCLAIMER
This report was prepared as an account of work sponsored by an agency of the United States
Government. Neither the United States Government nor any agency thereof, nor any of their
employees, makes any warranty, express or implied, or assumes any legal liability or responsi-
bility for the accuracy, completeness, or usefulness of any information, apparatus, product, or
process disclosed, or represents that its use would not infringe privately owned rights. Refer-
ence herein to any specific commercial product, process, or service by trade name, trademark,
manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recom-
mendation, or favoring by the United States Government or any agency thereof. The views
and opinions of authors expressed herein do not necessarily state or reflect those of the
United States Government or any agency thereof.
effect of the fluid on natural frequencies can be accounted for using the
c
v
Viscous damping coefficient kg/m-sec
e
i' g
i
Hydrodynamic forces per unit length on cylinder in N/m
x and y directions, respectively
(r. , 6., Z.) Local cylindrical coordinates associated with cylinder i (m, dimensionlesi
X X X *
m)
V Volume of a body m
13, *±y
Hydrodynamic masses per unit length kg/m
O". . , T . .
fluid pressures are generated. In turn, the integrated effect of these pressures
The fluid moving with the structure has an important effect on the dynamics
cases. In many practical situations, the effects of fluid viscosity and com-
mass.
The objective of this paper is to provide formulas and graphs for computing
to obtain. Therefore, fluid damping effect will be presented only for a few simple
the body behaves as though an added mass of fluid were rigidly attached to and
moving with it. When the body is subjected to an excitation, not only must the
mass of the body be accelerated, but also that of the added fluid mass. The
F - (M + m) ^ f + cv £ . (1)
Where M is the body mass, m is referred to as hydrodynamic mass (or added mass),
and c is the viscous damping coefficient. Note that two force components are
8u
(1) in — r , in phase with the acceleration, arises
3t
The hydrodynamic mass is proportional to the fluid density p and the body
m = p V cm , (2)
where c is the hydrodynamic mass coefficient. The sum of the hydrodynamic mass
and the body mass is called the virtual mass or the apparent mass.
Equation (1) corresponds to a body moving in a particular direction. For
hydrodynamic masses for multiple bodies with N degrees of freedom are represented
for multiple bodies with N degress of freedom, the number of terms needed to
motion and composite motion. In the former one, each body oscillates as an in-
dividual unit, while in the latter, all bodies oscillate as a rigid body. For
considered.
1. Individual Motion
fluid, as shown in Fig. 1. The axes of the cylinders are parallel and per-
1
f. = - I a. .
I
(3)
where a.., a.., T.., and f3.. are hydrodynamic masses. It also has been shown that
^3 i3 iJ ij
a.. = a.. ; 3.. = 3.. ;
Equations (3) show the hydrodynamic forces in the x and y direction due
u. = ulcosi - vising ,
v. = ujsir.v + v!cos'4> ,
1 1 1
(5)
f.. = f_!co3v - gjsinif; ,
8 =
i ,
where u' and v! are displacement components in the x' and y1 directions for
cylinder i. Substituting Eqs. (5) into (3), one can solve for f! and g! ;
k
f! = - : / a!. — ^ + a]
lj
at at
(6)
2 2 '
3 u! 3 v!
e' = - Z f T' 1 +
where
2 2
a1. . = a. . cos "> + S. . sin ^ + (a. . + x. . )sin^ COSIJJ ,
2 2
a! . = a. . cos ^ - x . .sin 41 + (6 . . - a. .) sini) cosij; ,
2 2
x! . = -a . .sin~!!> + x . .cos ty + ( 3 . . - a. ^) sin6 costjj , (7)
? 2
3 ! . = a. . sin~0 + 3 . . cos ^ - (a. . + x . .) sini> cosij;
Equations (7) define the transformation of the hydrodynamic mass. Naturally,
the symmetric properties are still preserved in the new coordinate system.
2. Composite Motion
A. Translational Motion
U
i " Ul •
v± = U 2 , (i = 1, 2, 3, ... k ) . (3)
The resultant force acting on the whole cylinder array can be resolved
k
r1 — i J. , v. --1
f
Substituting Eqs. (3) into (9) and using 8) yields
F
2 =-
where
a
12
" i = l 3=1
k k (11)
a = : i
21
i=l j=l
k k
a = Z I
22
i=l j=l
Since a.. =x
a12 - a n _ (12)
The symmetry of the hydrodynamic mass for composite motion is similar to that of
2 2
a' = a^ cos \}i + a. s i n ty + (a „ +• =, )sini> cosii/ ,
2 2
a
12 = a 2 i c o s ^ ~ a 2 1 s i n \p + ( a 2 2 - = i ; L )sin^ cosif- ,
2 2
a = - a s i n vp + a c o s
2
a = a s i n
22
B. Rotational Motion
Assume that the cylinder array rotates around the z-axis with an
90
angular velocity — . The cylinder accelerations associated with the
01
3 u. .2.
1 - . . 3 0
- rSln9
3
^4
where r. is the distance from the center of the cylinder to the z-axis, and r.
is the angle between the x-axis and the position vector r.. The torque is
given by
k
T = Z [-f.r.sine. + g.r.cosS.] . (15)
1 1 1 X 1 1
i=l
Substituting Eqs. (3) into (15) and using (14) yields
where
k k _
Y = E E (r r )
3
i=l j=l
+ £..cosO.cosS.)
The hydrodynainic masses for multiple bodies based on Eq. (3) can be
each body possesses N T degrees of freedom, the order of the matrix will be
N(=N' x k). Since [m..] is symmetric, one can find the N eigenvalues of the
hydrodynamic mass matrix [m..]. These eigenvalues are called effective hydro-
masses for multiple bodies play the same role as the hydrodynamic mass for a
presented in this section are also valid for other structural shapes, such as a
A generic word "hydrodynainic mass" will be used for referring to both hydrod;—.arsic
mass and hydrodynamic moment of inertia.
III. NATURAL FREQUENCIES
1. Single 3ody
equation:
f
f " fv (
M-hn > (17)
2. Multiple Bodies
The natural frequency in this case can be obtained in the same way as a
single body.
When all bodies can move freely these bodies will oscillate in a co-
tained for the natural frequencies of ccupled modes. However, when all bodies
liquid is available:
H
f - f
f (
<
f v M+p
frequencies.
Similar techniques can be employed for composite motion, provided the
given by [1]:
c = Re(H) , (19)
in
and
2
c 7 = -OTTR wlm(H) ,
where
4K (a)
H=1 i -*-
a = \/±
2
6jR
S (20)
V
, 2>/2~
c - 1
and
2/TT 2
-I
rz— S
- IQC-^^K.CCJ) + ^ ( 0 ) ^ ( 8 ) - I ^ C O K Q C S ) ] } - 1 ,
where
a = I'TS ,
3 =^>1S ,
r = ?v. ,'R, ,
and -,
c == —
C :
m" " " S±nh\\ l
ilShCto) (22)
}C—_L
where
a = In
11
except for a shell with simply supported end conditions. For a simply
(23)
c =
~ [I (a ))
s da n s
(24)
2-R
a =
Note that -_(:0 is the modified Bessel function of the first kind.
different r..
The hydr^dynaE: c mass matrices are as follows [6] for motion in the x
'R + R
1 2
"12
(25)
[m..] = :- R +
l V R
9
2 J 22
12
R + R
R l 2
lvll
R + (26)
i V
V'22
are given by
0
/ O O O 0 0 «»
2R
K < •o -2p (R. + R~) + & ; - R,) „ . „, ,„ , , „. 1
l"2 ll + c
1 2 22 2_ 1 _ £=1
Z £ coth(xa)exp(-2Za)
'12
1/2
0 9?
2 2 2
p - R - R (27)
2R R
1 2
2 ,"
1/2,
2
* < -
- R
-1
2pR,
2 1/2.
2 2 „? , T.2
1 3 R + R
2 ~ 1 2
a 2 = 2 Sin -1
A
2
where p = R. + ?..-, + G.
R /R =2.0. From Fig. 7, v1;L = 1.058. Then let R = 2.0, and R 2 = 1.0;
these give G/R, = 0.5 and Rj/R.. = 0.5. Again, from Fig. 7, v,, = 1.036.
=
Therefore, using the original subscripts of 3 = 1.0 and R ? = 2.0, v « 2
1.036.
F. A Group of Cylinders
ideal fluid (Fig. 1 ) . Each cylinder can nove in the x and y directions.
given by [2,3]
2k 3 2 u.
F, = 2 ™ i , i = i, 2, 3, .., 2k (28)
1 J
j=l 9t
where m.. is the element of the hydrodynanic rsass matrix. Equation (28)
2k/R. + R.\ 2 3 2 u-
F. ^iJ^T^j P'Vy ^ , (29)
where R. is the radius of cylinder i, and v.. is the element of the hydro-
dynamic mass coefficient matrix. When all cylinders are of the same radius,
m,. is equal to v.. multiplied by the fluid nass displaced by the cylinder.
Two coaxial shells are shown in Fig. 12. The hydrodynamic forces
acting on the two shells per unit area are 7, and F 2 , these are given by
a
n
ll
(30)
m« m.
'22
The hydrodynariic mass matrix [m..] per ur.it shell surface area is given by
R
1 V 11
- p
. .
"R2V21 R
2
(31)
(5 K
'21 > n
a)]
a3 [r (a) 10 (3) - i; (3) K^ (
where the prine denotes differentiation -..-i-h respect to the argument of the
(32)
V = V
12 21 (Contd.)
Table 1 shows the hydrodynamic masses taken from Refs. [7-12]. Note
C T = I /pi • (33)
X a.
presented in Table 2. These data are taken from Refs. [7-9,11,14-22] and
3. Computer Programs
Detailed descriptions of the program will not be presented here. However, the
Fortran listings of these programs are available upon request from the authors.
A. AMASS
diameters and arranged in any pattern [2,3]. The numerical technique used
in the program is based on the series expansion of the solution for the
potential flow.
16
2. AMASS-FEM
and continuity equations. This problem can be solved numerically via a finite-
velocity and pressure components, and used the Galerkin method to develop a
finite element computer program. This program calculates the hydrodynamic mass
V. CONCLUDING REMARKS
The methods for calculating hydrodynamic mass presented in this paper are
structure displacement becomes large, nonlinear effects of the fluid field become
important. In fact, it has been shown that the hydrodynamic mass is a function
of the vibration amplitude [24], and that the damping is a function of the first
The results presented for cylinders are based on two-dimensional flow theory;
i.e., the axial motion of the fluid is neglected. Two-dimensional theory tends to
is much larger than the cylinder diameter. Many practical situations satisfy
17
this requirement. However, when the wavelength is comparable with the diamecer,
important.
for large wavelengths. In the case of shell vibration, since the wavelength is
comparable with the diameter, the hydrodynamic mass becomes dependent on the
modes of vibration.
Very limited data are given for irregular shapes. The hydrodynamic masses
for two-dimensional bodies can be obtained using the finite element program in-
program is needed.
18
References
2. S. S. Chen, Vibration of Nuclear Fuel Bundles, Nucl. Eng. Des. ^5.. 399-422
(1976).
Confined Fluid, Trans. ASME, J. Appl. Mech. 99_, 213-217 (June 1977).
6. V. Yu. Mazurj Motion of Two Circular Cylinders in an Ideal Fluid, Izv. Akad.
REFERENCES (Contd.)
22. T . Sarpkaya, Added Masses of Lenses and Parallel Plates, J . Eng. Mech.
Div. , ASCE _86_5 141-151 (1960).
ACKNOWLEDGMENT
This work was performed under the sponsorship of the United States
spheres in an infinite
-*-f;
/»?•'
10
1.0
0.8
0.6
50
0.4
100
—•*
E 0.2
500
0 . 10
0 .08
0 .06
5000
0.04
0.0 2
50000
0.01 i I I I II II
20 40 60 80 100
or
CD
HYuRODYNAMiC MASS COEFFICIENT (C m )
HYDRODYNAMIC MASS COEFFICIENT {v\\)
C7>
0.
10
6/Rj
HYDRODYNAMIC MASS COEFFICIENT (?/,,)
o
0.2 0.4 0.6 .0
ECCENTRICITY, €
0
0 0.4 0.6 0.8 1.0
ECCENTRICITY, €
SHELL
SHELL
z
/ / / / / / / / / / / / / / I / / / / / / / / 7
FLUID
o >- z
2b
t
2b
£ !
o
»•-
tu 2a
o
Lu
U_
UJ
8
CO
CO
0.5
O
0.2-
— RECTANGULAR CYLINDER: WENDEL [ l l ]
ELLIPTIC CYLINDER : LAMB [ 7 ] AMD WENDEL [ l l ] i
0.1 I __J 1 1 . I
0. 0.2 0.5 I 2 10
ASPECT RATIO,a/b
8
T
16 2o
J_
2b
14
O
UJ
O
O
12 t
2a
UJ
10 2b-
O
8
CO
CO
ELLIPTIC
CYLINDER
x —RECTANGULAR
CYLINDER
I 4
IN AN INFINITE
FLUID MEDIUM
CONTAINING FLUID
0
0.01 0.02 0.5
/o.
p
o
Q
CO
p
to
p
CO
o
<O
O
o
CO
o
c>
CO
p
o
CM
O
O
O
O O o P pco* p p o
O CO o o
CO — CO
'1N3I31JJ300 SSWl OiWVNAGOaOAH
1.2
i.O
e
o
ma=
p.8
o
o
o
co 0.6
CO O
o
O
3 0.4
O
cc
THEORY
CD 6 WAUGH AND ELLIS [16]
0.2 EXPERIMENTS
A WAUGH AND ELLIS [16]
O PATTON [ 8 ] I
0 I , 1 ,' I
0 6 8 10
d/R
i—AT
PERPENDICULAR MOTION
.PARALLEL MOTION
0
3
n
CO.
3
li
1
1N3IOUJ300 SSVW
HORIZONTAL
MOTION, c m x
I i I 7-l-T-M
8 10 CD
a/b
HYDRODYNAMIC MASS COEFFICIENTS .
p p p p p
" l en b> '
0.80
0.75 <r
- 0.55 V =
IN-PLANE MOTION
OUT-OF-PLANE MOTION
0.50
0.25 0.125
CJ
ZSL
o
d
CD
- 0.025 5c
.2
£
.0
o
a
8 0.8 = /° a 2 b c m
CO
CO
<3C
o
0.6
5
o
0.4 EXPERIMENTS
SARPKAYA [22]
lul STELSON AND MAVIS [23]
0.2
0
0.2 0.5 10
b/a
Table 1. Simple two-direensional bodies
Rotational 0
b
Vertical
a
'[7,11],
Horizontal Trpa
a.
b
(See Fig. 13)
[l-(a/b)2]2
Rotational 3 P(b -a )
2(a/b)[l+(a/b) Z ]
(See Fig. 14)
Serii-circle
T i Rotational
(4 l\ 16
- 1
R \-2 " 4 ''70R
Thi--strip
v / / x /'/A h
Vertical trpa
77 4 3
Kf)
(a\
[8,10],
Rotational
8P a 16 ' (hi
(h « a)
Regular Triangle
—I
Table 1. Simple tvo-dimensional bodies (Contd.)
Sauare
yy/y/y
f
2a
Rotational 0.234 pna 0.2757
I
-*— 2a—*J
Rectangle
Hexagon
(with fluid gap)
E l l i p t i c Cylinder Containing
Fluid
(b 2 -a 2 )
Rotational 7- pab 2 2
b +a
(See Fig. 14)
Table 1. Simple two-dimensional bodies (Contd.)
I - -•
T~ /
/
/
Perpendicu-
lar to
2
I p
a3
T
1 (a/h)
3 (G/a)
[9]
2o /
/
/
/
a Wall
_L $
-—G -*-
Circular Disk
8 ,,3 m 8
Vertical p R
3 "m ,,3 3ff
p ii i\
16 5 .i 16
Hot at ion.i 1
4') ~3 " 4!iii
(h « a)
Elliptic Disk
4 2,, m 4k
— -rrpa b k c =
m 2, 3
]C I Vertical /2
pna b
[7]
r r~i—rm— 0 dO
2b V s i i i G + f —1 c o s
(h << a)
Rectangular Disk
m = pa be
a in
(b/a) (b/a)
Vertical c - -r 1-0.425 [14].
vn 4 l+(b/a)
(h « a)
Table 2. Simple three-dimensional bodies (Contd.)
h
/tan30 3
Vertical c
111
. J LJ - . J
J'M
- rsi
(in
(h « a)
2. Bodies o£ Revolution'^
Sphere (in an Infinite Fluid)
Any 2 D3
Direction - irpR
0.5 HI.
Table 2. Simple three-dimensional bodies (Contd.)
Any
Direction 17],
2 3 I8,16] e
Vertical 0.5 k
-j irpR k
(See Fig. 16)
- 2
T7.9].
v » + I]
,) ~ 2
Any
pUc cm • Re
«„ [17,18],
Direction
3 «u
(t:, - tipcud of
(See Fig. 18)
[7]t
pic,
Rotation
(about
(See F i g . 20)
minor axis)
3 > b , 1) = C I = ~ - 1 ,pab 2 (a ? +b 2 )
Ctrctil.ir Toroid
Vertical. pVc
Rotational,
(about Pic
x-axis)
Vertical
or "12
Horizontal
sym
Veitical
sym
ZV22
(See Figs. 22 and 23)
Table 2. Simple three-dimensional bodies (Contd.)
Vertica 1. pa be c [22,23]^
in in
(Sec F i g . 24)
b
J
Hydrodynamic mass coefficient is defined noc as the ratio of hydrodynamic mass to displaced fluid mass by
structure but as a ratio independent of thickness, h. Its definition is shown in the fourth column.
ftHycirodynamic mass coefficient is defined as the ratio of hydrodynamic mass to displaced fluid mass,
tt+t refers to theoretical results and e refers to experimental values.