Transient acoustic wave propagation

This problem contains basic test cases for one or more Abaqus elements and features.

This page discusses:

ProductsAbaqus/StandardAbaqus/Explicit

Elements tested

  • AC1D2
  • AC1D3
  • AC2D3
  • AC2D4
  • AC2D4R
  • AC2D6
  • AC2D8
  • AC3D4
  • AC3D5
  • AC3D6
  • AC3D8
  • AC3D8R
  • AC3D10
  • AC3D15
  • AC3D20
  • ACAX3
  • ACAX4
  • ACAX4R
  • ACAX6
  • ACAX8
  • ASI1
  • ASI2D2
  • ASI2D3
  • ASI3D3
  • ASI3D4
  • ASI3D6
  • ASI3D8
  • ASIAX2
  • ASIAX3

Features 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 x-axis), and the acoustic medium has a prescribed constant inward particle acceleration of a0= 1 m/s2 at x= 0 m. A nonreflective boundary is specified at x= 1 m using nonreflective impedance conditions.

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

p(x,t)=Kfa0c2(ct-x),        x<ct,
p(x,t)=0,        xct.

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 constraint, and an inward volume acceleration is specified on degree of freedom 8 as a concentrated load (“afav” files).

  2. In Abaqus/Standard an ASI element is placed at x= 0 with its normal pointing into the fluid (this activates the displacement degree of freedom on the node at x= 0), and in Abaqus/Explicit a structural element with a tie constraint is used to define the interaction between the fluid and structure. An acceleration is prescribed directly with a boundary condition (“afas” files). In these cases the first time interval in the analysis is performed using boundary impedances; the analysis continues in time using surface impedances

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

Results and discussion

Numerical results for the pressure distribution along the length of the acoustic medium are compared with the analytical results at time t= 0.0044 s.

In Abaqus/Standard for cases where the volume acceleration is specified (“afav” input files), the maximum errors are below 0.7%, except for linear tetrahedra, which are within 3% of the analytical results. For cases where a boundary acceleration is specified on the nodes of the ASI elements (“afas” input files), the maximum errors are below 4.2%, except for linear tetrahedra, which are within 5.9% of the analytical results.

In Abaqus/Explicit the results for pressure are within 0.2% of the analytical results for all tests, except for linear tetrahedra, which are within 3.2% of the 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.

ec35afav.inp

AC3D5 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

ASI2D2/AC2D3 elements.

ec24afas.inp

ASI2D3/AC2D4 elements.

ec26afas.inp

ASI2D3/AC2D6 elements.

ec28afas.inp

ASI2D3/AC2D8 elements.

ec34afas.inp

ASI3D4/AC3D4 elements.

ec34afas_po.inp

POST OUTPUT analysis.

ec36afas.inp

ASI3D3/ASI3D4/AC3D6 elements.

ec35afas.inp

ASI3D4/AC3D5 elements.

ec38afas.inp

ASI3D4/AC3D8 elements.

ec3aafas.inp

ASI3D6/AC3D10 elements.

ec3fafas.inp

ASI3D6/AC3D15 elements.

ec3kafas.inp

ASI3D8/AC3D20 elements.

eca3afas.inp

ASIAX2/ACAX3 elements.

eca4afas.inp

ASIAX2/ACAX4 elements.

eca6afas.inp

ASIAX3/ACAX6 elements.

eca8afas.inp

ASIAX3/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.