Coupled Frequency Response Analysis of Fluid-Structure Models

Coupled Frequency Response Analysis of fluid-structure models, commonly termed Acoustic Analysis, is generally performed to model sound propagation within a structural cavity, such as the interior of a vehicle or a musical instrument.

OptiStruct allows both direct and modal frequency response analysis for fluid-structure models. The responses of both the structural and fluid domains are computed; for the structural domain, these responses are the displacements and rotations of the structural grids, and for the fluid domain, these responses are the pressures at the fluid grid points.

The accelerations of structural grids at the fluid-structure interface excite the fluid domain and conversely, the pressures on the fluid grids at the fluid-structure interface excite the structural domain. Hence the problem is coupled and the motions of structural and fluid degrees-of-freedom are solved simultaneously.

Loading is sinusoidal with excitation frequency $Ω$ , and can be in the form of forces, acoustic sources enforced displacements, enforced velocities, and/or enforced accelerations. Frequency response loads and boundary conditions are defined in the Bulk Data section of the input deck. These are then referenced in a Subcase definition through an SPC or DLOAD data selector.

Damping may be defined for both the structural and fluid domains. For the structural domain damping may be defined through structural damping elements, material damping, structural damping (PARAM,G), or modal damping (SDAMPING referenced by a SDAMPING(STRUCT) subcase data selector). For the fluid domain damping may be defined through material damping, fluid damping (PARAM,GFL), or modal damping (SDAMPING referenced by a SDAMPING(FLUID) subcase data selector). In addition, the normalized admittance coefficient for porous materials can be specified by ALPHA on the MAT10 data.

Frequency dependent fluid acoustic absorber elements can be specified on the fluid faces of the fluid-structure boundary using the CAABSF elements. The absorber elements can be point, line, triangular, or quadrilateral in shape. The CAABSF data references the PAABSF data which is used to specify the frequency dependent resistance (real part of the impendence) and reactance (imaginary part of the impendence) as well as the area factors for point and line elements.

PARAM, LFREQFL and PARAM, HFREQFL can be used to exclude modes from a Coupled Modal Frequency Response Analysis (Acoustic Analysis).

Wave Equation

The acoustic analysis is based on inviscid flow with linear pressure-density relation as:(1)

\frac{1}{ρ} \nabla p + \ddot{u} = 0

and the continuity equation is:(2)

p + β (\nabla \cdot u) = 0

Where,

$p$: Pressure of the fluid domain
$u$: Displacement of the fluid domain
$β$: Compressibility of the fluid domain
$ρ$: Density of the fluid domain

Combining the above equations, the governing equation of the fluid domain is:(3)

\frac{\ddot{p}}{β} - \frac{1}{ρ} \nabla^{2} p = 0

Effect of the Structure on the Fluid Domain at the Interface

After finite element discretization, the assembly of equations for the fluid domain is:(4)

M_{F} \ddot{p} + C_{F} \dot{p} + K_{F} p - A_{i n t} \ddot{u} = s_{F}

Where, $M_{F}$ , $C_{F}$ , $K_{F}$ , and $s_{F}$ are the mass matrix, damping matrix, stiffness matrix and source vector respectively, of the fluid domain.

The matrix $A_{i n t}$ represents the interface matrix and $\ddot{u}$ is the acceleration of the structural grids at the fluid-structure interface. (The pressure gradient at the interface will be influenced by the acceleration of the structural grids).

Effect of Fluid on the Structural Domain at the Interface

The structural equation assembly can be written as:(5)

M_{S} \ddot{u} + C_{S} \dot{u} + K_{S} u - A_{i n t}^{T} p = s_{S}

Where, $M_{S}$ , $C_{S}$ , $K_{S}$ , and $s_{S}$ are the mass matrix, damping matrix, stiffness matrix and source vector, respectively of the structural domain.

The matrix $A_{i n t}^{T}$ represents the transpose of the interface matrix and $p$ is the pressure at the interface fluid grids at the fluid-structure interface (The displacement, velocity and acceleration of the structural grids at the interface will be influenced by the pressure at the interface fluid grids).

Coupled Fluid-Structure Interface Equation

Therefore the combined fluid-structure interface equation is:(6)

[\begin{matrix} M_{S} & 0 \\ - A_{i n t} & M_{F} \end{matrix}] [\begin{matrix} \ddot{u} \\ \ddot{p} \end{matrix}] + [\begin{matrix} C_{S} & 0 \\ 0 & C_{F} \end{matrix}] [\begin{matrix} \dot{u} \\ \dot{p} \end{matrix}] + [\begin{matrix} K_{S} & A_{i n t}^{T} \\ 0 & K_{F} \end{matrix}] [\begin{matrix} u \\ p \end{matrix}] = [\begin{matrix} s_{S} \\ s_{F} \end{matrix}]

The above equations are solved simultaneously for unknowns in the structural and the fluid domains, either by direct frequency response or modal frequency response. For modal frequency response, OptiStruct will calculate the eigenspace for both structure and fluid domain automatically.

Loads on the Fluid Mesh

The fluid grid points can be loaded by specifying the load magnitude and GRID using the SLOAD data.

The SLOAD data is referenced by the ACSRCE data which defines the dynamic characteristics of the load (DELAY, DPHASE, and a tabular listing of the load scale factor versus frequency). In addition, the density and bulk modulus of the loaded fluid are specified on the ASRCE data. The material characteristics of the fluid must be specified in the ASRCE in case the same fluid GRID is shared by two different fluid meshes.

Fluid-Structure Interface Visualization and Refinement

OptiStruct has support for both grid-to-grid matching and non-matching interfaces. The interface is specified through the ACMODL card. If an ACMODL card is not specified in the input deck, the fluid-structure interface is automatically defined by OptiStruct based on default values for the ACMODL parameters.

Based on a search box specified on the ACMODL card, OptiStruct outputs an *.interface file, containing information about the fluid-structure interface. With the model loaded in HyperMesh, you can import the *.interface file to visualize the fluid-structure interface (ensure that the FE overwrite option is activated on import). The "^Fluid Faces at Interface" component is created, which allows you to view the interface found between the structural and fluid domains. If a component "^Acoustically Rigid Fluid Faces" is created, that means at those fluid surfaces, there are no structural grids found. A structural grid set "^Structural grids at Interface" is also created to display the structural grids found at the interface.

There are several steps to take to improve the interface:

Perform an OptiStruct check run. This will create the *.interface file, allowing you to visualize the interface.
If the interface is not as desired, you may create a new SET containing those fluid grids that describe the fluid boundary.
You can then specify the newly created set on the ACMODL card under FSET.
Perform another OptiStruct check run, and review the new interface.

External Fluid-Structure Coupling Definition File

OptiStruct can use the binary ftn.70 coupling file generated by AKUSMOD instead of internally calculating the coupling.

To use this option, add PARAM,AKUSMOD,YES to the input deck.

Modal and Panel Participation

Modal participation is a measure of how much each mode participates at a given frequency in a modal frequency response calculation. Output of modal participation may be requested for structural degrees of freedom as well as for fluid degrees of freedom.

Calculation and output of modal participation can be requested, for any number of degrees of freedom, using the PFMODE I/O Options Entry.

Panel participation is a measure of the influence of sets of specified structural grids, defined by PANEL Bulk Data Entries. The response of a fluid grid is influenced through each panel or each grid at the acoustic interface. Calculation and output of the contribution from each panel at specific loading frequencies can be requested through the PFPANEL I/O Options Entry, for modal frequency response. Also, the calculation and output of the contribution from each grid at the interface can be requested through the PFGRID I/O Options Entry.

A file *.pfmode.pch is generated based on the definition of the PFMODE and PFPANEL I/O Options Entry. The output for PFGRID would be in a H3D file. The results for PFMODE and PFPANEL are best plotted in HyperGraph, whereas the contour results for PFGRID are best visualized in HyperView.

Non-Reflecting Boundary

Allows the acoustic wave to propagate normally through the boundary, without reflection. This condition is called the Sommerfeld boundary condition.

To create a non-reflecting boundary, set the values of the TABLEDi entry referenced by the TZREID field (Resistance-real part of Impedance) in the PAABSF Bulk Data Entry to be equal to $(ρ_{f l u i d}) * (c_{f l u i d})$ for all frequencies.

Where,

$ρ_{f l u i d}$: Density of the fluid
$c_{f l u i d}$: Speed of sound in the fluid