Products: ABAQUS/Standard ABAQUS/Explicit
Problems in unbounded domains, or problems in which the region of interest is small compared to the surrounding medium, are common in acoustic analysis. For example, a submarine deep underwater may experience loads due to the fluid and radiate sound into the fluid, as if the ocean were infinitely large. The extent to which an exterior fluid may be considered “unbounded” or “infinite” depends on the number of wavelengths between the body or region of interest and the nearest boundary: the higher this number, the more likely that the influence of these boundaries is small enough to be neglected. For example, the effect of the surrounding fluid on a submarine in a relatively shallow harbor may be like an infinite medium at high frequencies, but the harbor bottom and free surface may exert an effect at lower frequencies. Similarly, a loudspeaker in air may radiate sound from the high-frequency tweeter as if the surrounding air were infinite, but the effects of the walls of the room may affect the radiation pattern of the low-frequency woofer.
ABAQUS provides a range of features to model exterior fluid effects. All of the surrounding fluid may be modeled with finite elements; clearly, this is practical only if the extent of the surrounding medium is small. At the next level of sophistication, the user may model a small region of fluid and apply a simple radiation boundary condition to the terminating surface. These radiation boundary conditions are derived using simple models of waves passing through a boundary. ABAQUS provides several alternative formulations of these boundary conditions (see “Coupled acoustic-structural medium analysis,” Section 2.9.1, and “Acoustic, shock, and coupled acoustic-structural analysis,” Section 6.9.1 of the ABAQUS Analysis User's Manual). These radiation conditions do not add degrees of freedom to the system and do not affect the symmetry of the matrix.
Finally, acoustic infinite elements are provided, which allow the retained finite element fluid region to be even smaller, with similar accuracy. The acoustic infinite element formulation differs from the radiation boundary condition formulation in several key respects. In the infinite elements the infinite exterior is subdivided into elements, and a method of weighted residuals statement is enforced on the elements in a manner entirely analogous to the usual finite element method. Degrees of freedom, corresponding to interpolation functions in the infinite direction, are added to the overall matrix system. The method of weighted residuals used in ABAQUS results in nonsymmetric infinite element matrices, so the relative cost of these elements is higher than that of a simple radiation boundary condition. However, the accuracy of the infinite elements is sufficiently high that the finite element region can be reduced considerably, offsetting the cost in many applications.
The solution in the unbounded acoustic medium is assumed to be linear and governed by the same equations as the finite acoustic region:
is used here to denote the volumetric drag parameter. Consider the infinite exterior of a region of acoustic fluid bounded by a convex surface and a conventional finite element mesh defined on this surface. Each facet of this surface mesh, together with the normal vectors at the nodes, defines a subdivision of the infinite exterior that will be referred to as the “infinite element.” Application of the method of weighted residuals results in a weak form of this equation over the infinite element volume:
, and for the solution field, 
, will be different. Selection of these functions has been the subject of considerable experimentation and analysis (for example, see Allik, 1991, and Burnett, 1994). In ABAQUS considerations of accuracy, numerical stability, time-domain well-posedness, and economy have led to the selection of a form proposed by Astley (Astley, 1994). Here, the Fourier transform is applied to the equilibrium equation, and the functions for the solution field are chosen as tensor products of conventional finite element shape functions in the directions tangential to the terminating surface and polynomials in 
, where r is a coordinate in the infinite direction, and a spatially oscillatory factor. The weights are chosen as complex conjugates of the solution field functions times a 
factor. This combination has shown to result in an element with several desirable properties. First, if the material properties are constant with respect to frequency in steady-state analysis, the mass, damping, and stiffness element matrices are constant as well; equivalently, in transient analysis the element results in a second-order differential equation for the pressure. Moreover, the damping matrix has positive eigenvalues, a necessary condition for well-posedness in transient analysis. Analytical investigations of the formulation demonstrate that the element captures the exact solutions for radiation impedances for modes of a sphere; the order of the spherical mode modeled exactly is equal to the order of the polynomial used in the infinite direction. The element integrals do not contain oscillatory kernels and can be evaluated using standard Gauss quadrature methods. Finally, the element can be formulated using the usual coordinate map due to Bettess (1984) so that arbitrary convex terminating surfaces can be used.To continue with the derivation, we transform the weighted residual statement into the frequency domain:
Now the weight and solution interpolation functions are defined as
, 
, and 
are the parent element coordinates; 
; 
; 
; and 
are element shape functions with indices varying over the number of degrees of freedom of the element. The mapping and shape functions are described below.Inserting the shape functions into Equation 3.3.2–1 and integrating by parts, we obtain the following element equation:
This equation is clearly nonsymmetric, due to the fact that the weight and trial functions are not identical. Nevertheless, if the material properties are constant as a function of frequency, each element matrix corresponding to the terms above is constant as well. Gradients of density inside the element volume have been ignored in the formulation.
The element shape functions are defined as follows:
ranges over the n nodes of the infinite element on the terminating surface, while the subindex 
ranges over the number of functions used in the infinite direction. The function index i is equal to the subindex for the first n functions, 
for the second n functions, and so on.The functions 
are conventional two-dimensional shape functions (in three dimensions) or, with 
, one-dimensional shape functions for axisymmetric or two-dimensional elements. The role of these functions is to specify the variation of the acoustic field in directions tangent to the terminating surface. The variation of the acoustic field in the infinite direction is given by the functions 
, which are members of a set of ten ninth-order polynomials in 
. The members of this set are constructed to correspond to the Legendre modes of a sphere. That is, if infinite elements are placed on a sphere and if tangential refinement is adequate, an ith order acoustic infinite element will absorb waves associated with the (
)th Legendre mode. The first member of this set corresponds to the value of the acoustic pressures on the terminating surface; the other functions are generalized degrees of freedom. Specifically,
are equal to zero at the terminating surface, except for 
. The coordinate map is described in part by the element shape functions, in the usual isoparametric manner, and in part by a singular function. Together they map the true, semi-infinite domain onto the parent element square or cube. To specify the map for a given element, we first define distances 
between each infinite element node 
on the terminating surface and the element's reference node, located at 
:
is the “node ray” at the node. The “node ray” at a node is the unit vector in the direction of the line between the reference point and the node. We use the interpolated reference distance on the terminating surface,
, corresponding to the infinite direction:
. With these definitions, the geometric map can be specified as
The infinite elements are not isoparametric, since the map uses a lower-order function of the parent coordinates than the interpolation scheme does. However, this singular mapping is convenient and invertible.
The phase factor 
appears in the element formulation. This factor models the oscillatory nature of the solution inside the infinite element volume but must also satisfy some additional properties. First, it must be zero at the face of the infinite element that connects to the finite element mesh. Second, it must be 
continuous across infinite element lateral boundaries. Third, it should be such that the mass-like term in the infinite element equation,
, the definition 
is the only departure from the definition used in Astley (1994).In a transient analysis the element matrix equation is transformed back to the time domain, using constant material properties taken from the value at zero frequency. The resulting second-order ordinary differential equation for the pressure degrees of freedom is added to the overall system in the model and integrated in the usual manner.
In steady-state analysis the formulation for the acoustic infinite elements is consistent with that used for acoustic finite elements (see “Coupled acoustic-structural medium analysis,” Section 2.9.1, and “Acoustic, shock, and coupled acoustic-structural analysis,” Section 6.9.1 of the ABAQUS Analysis User's Manual). First, the complex density
The expression for the basis and weight functions is formally identical to the preceding transient development, but with the complex wave number
is the same as for the transient case. It is real-valued. The remaining terms are complex; and, as in the case for acoustic finite elements, they are manipulated into real and imaginary parts. We use
An infinite element volume may have an impedance boundary condition defined on a face extending to infinity. Such a definition can be used to model, for example, an infinite half-space of water above a lossy seabed or under a layer of ice. The derivation of the contribution of the impedance condition to the infinite element matrix is straightforward and follows the derivation for finite elements (see “Coupled acoustic-structural medium analysis,” Section 2.9.1). The contributing term in the transient case is
Lateral faces of acoustic and solid infinite elements (that is, those that extend to infinity) can be coupled to each other. For example, two infinite half-spaces may interact at their interface. As in the case of lateral-face impedance terms, the derivation of the contribution of the coupling condition to the infinite element matrix follows the derivation for finite elements to the greatest degree possible. We restrict the consideration to solid medium infinite elements that share an edge with the acoustic medium infinite elements and where the acoustic infinite element is the slave in a tied contact condition. The contributing terms in the transient case are
In tied contact between fluid and solid medium elements, a tributary area of the slave node is defined, and the conservation of momentum and continuity of displacement is enforced between the slave node and the master surface. On a lateral face of an infinite element, the “area” is unbounded, so an area measure based on the acoustic infinite element formulation,
, is evaluated at the shared edge to compute this coupling condition. For consistency the boundary traction on the solid medium infinite element due to the acoustic pressure is defined as
