13.38. FLUID38 - Dynamic Fluid Coupling

Matrix or VectorShape Functions Integration Points
Mass Matrix

None
Damping MatrixNot definedNone

Reference: Fritz([12])

13.38.1. Description

This element is used to represent a dynamic coupling between two points of a structure. The coupling is based on the dynamic response of two points connected by a constrained mass of fluid. The points represent the centerlines of concentric cylinders. The fluid is contained in the annular space between the cylinders. The cylinders may be circular or have an arbitrary cross-section. The element has two DOFs per node: translations in the nodal x and z directions. The axes of the cylinders are assumed to be in the nodal y directions. These orientations may be changed with KEYOPT(6).

13.38.2. Assumptions and Restrictions

  1. The motions are assumed to be small with respect to the fluid channel thickness.

  2. The fluid is assumed to be incompressible.

  3. Fluid velocities should be less than 10% of the speed of sound in the fluid.

  4. The flow channel length should be small compared to the wave length for propagating vibratory disturbances (less than about 10%), in order to avoid the possibility of standing wave effects.

13.38.3. Mass Matrix Formulation

The mass matrix formulation used in the element is of the following form:

(13–50)

The m values are dependent upon the KEYOPT(3) value selected. For KEYOPT(3) = 0 (concentric cylinder case):

(13–51)

(13–52)

(13–53)

where:

ρ = fluid mass density (input as DENS on MP command)
R1 = radius of inner cylinder (input as R1 on R command)
R2 = radius of outer cylinder (input as R2 on R command)
L = length of cylinders (input as L on R command)

Note that the shape functions are similar to that for PLANE25 with MODE = 1. The element mass used in the evaluation of the total structure mass is .

For KEYOPT(3) = 2, which is a generalization of the above cylindrical values but for different geometries, the m values are as follows:

(13–54)

(13–55)

(13–56)

(13–57)

(13–58)

(13–59)

where:

M1 = mass of fluid displaced by the inner boundary (Boundary 1) (input as M1 on R command)
M2 = mass of fluid that could be contained within the outer boundary (Boundary 2) in absence of the inner boundary (input as M2 on R command)
Mhx, Mhz = hydrodynamic mass for motion in the x and z directions, respectively (input as MHX and MHZ on R command)

The element mass used in the evaluation of the total structure mass is M2 - M1.

The lumped mass option (LUMPM,ON) is not available.

13.38.4. Damping Matrix Formulation

The damping matrix formulation used in the element is of the following form:

(13–60)

The c values are dependent upon the KEYOPT(3) value selected. For KEYOPT(3) = 0:

(13–61)

(13–62)

(13–63)

(13–64)

where:

Wx, Wz = estimate of resonant frequencies in the x and z response directions, respectively (input as WX, WZ on RMORE command)
f = Darcy friction factor for turbulent flow (input as F on R command)
Δx, Δz = estimate of peak relative amplitudes between inner and outer boundaries for the x and z motions, respectively (input as DX, DZ on R command)

For KEYOPT(3) = 2, the c values are as follows:

(13–65)

(13–66)

(13–67)

(13–68)

where:

Cx, Cz = flow and geometry constants for the x and z motions, respectively (input as CX, CZ on RMORE command)