8.4. Acoustic Fluid-Structure Interaction (FSI)

8.4.1. Coupled Acoustic Fluid-Structural System with an Unsymmetric Matrix Equation

The coupling conditions on the interface between the acoustic fluid and the structure are given by:

(8–150)

(8–151)

where:

p = acoustic pressure

Equation 8–150 is a kinetic condition relating the solid stress to the pressure imposed on the interface by sound. Equation 8–151 is a kinematic condition that assumes that there is no friction between the solid and acoustic fluid on the interface.

In order to completely describe the FSI problem, the fluid pressure load acting at the interface is added to Equation 15–6. This effect is included in FLUID29, FLUID30, FLUID220, and FLUID221 only if KEYOPT(2) ≠ 1. Hence, the structural equation is rewritten as:

(8–152)

The fluid pressure load vector at the interface S is obtained by integrating the pressure over the area of the surface as follows:

(8–153)

where:

{N'} = shape functions employed to discretize the displacement components u, v, and w (obtained from the structural element).

Substituting the finite element approximating function for pressure given by Equation 8–34 into Equation 8–153 leads to:

(8–154)

By comparing the integral in Equation 8–154 with the matrix definition of [R]T in Equation 8–37, the following relation becomes clear:

(8–155)

Substituting Equation 8–155 into Equation 8–152 results in the dynamic elemental equation of the structure, expressed as:

(8–156)

Equation 8–37 and Equation 8–156 describe the complete finite element discretized equations for the FSI problem. These equations are written in assembled form as:

(8–157)

The acoustic fluid element in an FSI problem will generate all the submatrices with a superscript F in addition to the coupling submatrices , [R]T, and [R]. Submatrices with a superscript S will be generated by the compatible structural element used in the model.

Assuming that the actual surface is at an elevation η relative to the mean surface in z-direction, the pressure for a sloshing (free) surface is given by:

(8–158)

By utilizing the definition of velocity and the momentum conservation equation in addition to Equation 8–158, pressure can be expressed as:

(8–159)

The surface integration of the "weak" form (Equation 8–37) on the sloshing surface is given by:

(8–160)

The acoustic fluid matrix equation with sloshing effect is expressed as:

(8–161)

where:

Substituting Equation 8–161 into Equation 8–157 yields:

(8–162)

If the impedance boundary is exerted on the FSI interface (input as IMPD on the SF command), the coupling condition expressed in Equation 8–151 is rewritten as:

(8–163)

Substituting Equation 8–163 into Equation 8–37 yields:

(8–164)

where:

Damping matrix [CFSI] on the impedance FSI interface has been shown to be the same as the damping matrix in Equation 8–40. Therefore, the coupling matrix expressed in Equation 8–162 can still be the final matrix equation. In an incompressible fluid the fluid density is independent of the pressure. This implies the speed of sound equivalently tends toward infinity. The matrix in Equation 8–162 is set to zero.

8.4.2. Coupled Acoustic Fluid-Structural System with Symmetric Matrix Equation for Full Harmonic Analysis

The matrix in Equation 8–162 has been shown to be unsymmetric. Solving Equation 8–162 may consume more computer resources and time than solving the symmetric matrix equation. For the frequency domain, assume that:

(8–165)

Substituting Equation 8–165 into Equation 8–156 and Equation 8–37 yields:

(8–166)

(8–167)

Dividing coupled Equation 8–167 by yields the acoustic matrix equation, written as:

(8–168)

The coupled matrix equation is given by:

(8–169)

After solving Equation 8–169, the pressure is obtained using Equation 8–165.