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:
To continue with the derivation, we transform the weighted residual statement into the frequency domain:
Now the weight and solution interpolation functions are defined asInserting 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:
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,
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 :
We use the interpolated reference distance on the terminating surface,
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,
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
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,