3.7.3 Transient acoustic wave propagation

Products: ABAQUS/Standard  ABAQUS/Explicit  

Elements tested

AC1D2    AC1D3   

AC2D3    AC2D4    AC2D4R    AC2D6    AC2D8   

AC3D4    AC3D6    AC3D8    AC3D8R    AC3D10    AC3D15    AC3D20   

ACAX3    ACAX4    ACAX4R    ACAX6    ACAX8   

ASI1    ASI2    ASI3   

ASI4    ASI8   

ASI2A    ASI3A   

Feature tested

Transient wave propagation in an acoustic medium with nonreflective end condition.

Problem description

The model consists of a column of fluid 1 m long with a constant cross-sectional area. The tube lies horizontally (along the -axis), and the acoustic medium has a prescribed constant inward particle acceleration of 1 m/s2 at 0 m. A nonreflective boundary is specified at 1 m using the nonreflective feature of the *IMPEDANCE and *SIMPEDANCE options.

The acoustic material in the column is air with a bulk modulus 1.424 × 105 N/m2 and a density 1.21 kg/m3. The speed of sound is calculated as = 343.05 m/s. The analytical result for the pressure is

The column is modeled using either 100 first-order or 50 second-order acoustic elements. For each acoustic element tested, the acceleration is specified in each of two ways:

  1. There is no ASI element or *TIE option, and an inward volume acceleration is specified on degree of freedom 8 as a *CLOAD (“afav” files).

  2. In ABAQUS/Standard an ASI element is placed at 0 with its normal pointing into the fluid (this activates the displacement degree of freedom on the node at 0), and in ABAQUS/Explicit a structural element with the *TIE option is used to define the interaction between the fluid and structure. An acceleration is prescribed directly with *BOUNDARY, TYPE=ACCELERATION (“afas” files). In these cases the first time interval in the analysis is performed using the *IMPEDANCE option; the analysis continues in time using the *SIMPEDANCE option.

A transient dynamic analysis is performed for a period long enough to allow the wave to propagate past the nonreflective boundary.

Results and discussion

For both ABAQUS/Standard and ABAQUS/Explicit the results for pressure are within 0.4% of the analytical result for all tests, except for linear tetrahedra, which are within 3% of analytical results.

Input files

ABAQUS/Standard input files

ec12afav.inp

AC1D2 elements.

ec13afav.inp

AC1D3 elements.

ec23afav.inp

AC2D3 elements.

ec24afav.inp

AC2D4 elements.

ec26afav.inp

AC2D6 elements.

ec28afav.inp

AC2D8 elements.

ec34afav.inp

AC3D4 elements.

ec36afav.inp

AC3D6 elements.

ec38afav.inp

AC3D8 elements.

ec3aafav.inp

AC3D10 elements.

ec3fafav.inp

AC3D15 elements.

ec3kafav.inp

AC3D20 elements.

eca3afav.inp

ACAX3 elements.

eca4afav.inp

ACAX4 elements.

eca6afav.inp

ACAX6 elements.

eca8afav.inp

ACAX8 elements.

ec12afas.inp

ASI1/AC1D2 elements.

ec13afas.inp

ASI1/AC1D3 elements.

ec23afas.inp

ASI2/AC2D3 elements.

ec24afas.inp

ASI2/AC2D4 elements.

ec26afas.inp

ASI3/AC2D6 elements.

ec28afas.inp

ASI3/AC2D8 elements.

ec34afas.inp

ASI4/AC3D4 elements.

ec36afas.inp

ASI4/AC3D6 elements.

ec38afas.inp

ASI4/AC3D8 elements.

ec3aafas.inp

ASI4/AC3D10 elements.

ec3fafas.inp

ASI4/AC3D15 elements.

ec3kafas.inp

ASI8/AC3D20 elements.

eca3afas.inp

ASI2A/ACAX3 elements.

eca4afas.inp

ASI2A/ACAX4 elements.

eca6afas.inp

ASI3A/ACAX6 elements.

eca8afas.inp

ASI3A/ACAX8 elements.

ABAQUS/Explicit input files

eca3afav_xpl.inp

ACAX3 elements.

eca4arav_xpl.inp

ACAX4R elements.

ec23afav_xpl.inp

AC2D3 elements.

ec24arav_xpl.inp

AC2D4R elements.

ec34afav_xpl.inp

AC3D4 elements.

ec36afav_xpl.inp

AC3D6 elements.

ec38arav_xpl.inp

AC3D8R elements.

eca3afas_xpl.inp

ACAX3 elements.

eca4aras_xpl.inp

ACAX4R elements.

ec23afas_xpl.inp

AC2D3 elements.

ec24aras_xpl.inp

AC2D4R elements.

ec34afas_xpl.inp

AC3D4 elements.

ec36afas_xpl.inp

AC3D6 elements.

ec38aras_xpl.inp

AC3D8R elements.