Acoustic scattering of spherical waves incident on a long fluid-saturated poroelastic cylinder

Acoustic scattering of spherical waves generated by a monopole point source in a perfect (inviscid and ideal) compressible fluid by a fluid-saturated porous cylinder of infinite length is studied theoretically in the present study. The formulation utilizes the Biot theory of dynamic poroelasticity along with the appropriate wave-field expansions, the translational addition theorem for spherical wave functions, and the pertinent boundary conditions to obtain a closed-form solution in the form of infinite series. The analytical results are illustrated with a numerical example in which a monopole point source within water is located near a porous cylinder with a water-saturated Ridgefield sandstone formation. The numerical results reveal the effects of source excitation frequency, the cylinder interface permeability condition, and the location of the point source and the field point on the backscattered pressure magnitudes. Limiting cases are considered, and the obtained numerical results are validated by already well-known solutions.     

Wave propagation passing over a submerged porous breakwater

A linear model of waves propagating over a submerged porous breakwater is derived from two coupled boundary-value problems, each of which represents the governing equation in a different medium. The model is similar to the shallow-water equations (SWE), with a damping term proportional to the character of the porous breakwater. Therefore, waves traveling above the breakwater will be absorbed, and the amplitude decreases. The wave propagation passing over the submerged breakwater for monochromatic and solitary waves is analyzed. For monochromatic waves, the numerical solution agrees with the analytical. The amplitude decreases exponentially with respect to the space variable in the region above the breakwater. The reflected wave is also analyzed when the model is combined with a model using the shallow-water equations.     

Modelling of crack propagation of gravity dams by scaled boundary polygons and cohesive crack model

Crack propagation in concrete gravity dams is investigated using scaled boundary polygons coupled with interface elements. The concrete bulk is assumed to be linear elastic and is modelled by the scaled boundary polygons. The interface elements model the fracture process zone between the crack faces. The cohesive tractions are modelled as side-face tractions in the scaled boundary polygons. The solution of the stress field around the crack tip is expressed semi-analytically as a power series. It reproduces the singular and higher-order terms in an asymptotic solution, such as the William’s eigenfunction expansion when the cohesive tractions vanish. Accurate results can be obtained without asymptotic enrichment or local mesh refinement. The stress intensity factors are obtained directly from their definition and provide a convenient and accurate means to assess the zero-K condition, which determines the stability of a cohesive crack. The direction of crack propagation is determined from the maximum circumferential stress criterion. To accommodate crack propagation, a local remeshing algorithm that is applicable to any polygon mesh is augmented by inserting cohesive interface elements between the crack surfaces as the cracks propagate. Three numerical benchmarks involving crack propagation in concrete gravity dams are modelled. The results are compared to the experimental and other numerical simulations reported in the literature.     

A simple multi‐directional absorbing layer method to simulate elastic wave propagation in unbounded domains

The numerical analysis of elastic wave propagation in unbounded media may be difficult due to spurious waves reflected at the model artificial boundaries. This point is critical for the analysis of wave propagation in heterogeneous or layered solids. Various techniques such as Absorbing Boundary Conditions, infinite elements or Absorbing Boundary Layers (e.g. Perfectly Matched Layers) lead to an important reduction of such spurious reflections. In this paper, a simple absorbing layer method is proposed: it is based on a Rayleigh/Caughey damping formulation which is often already available in existing Finite Element softwares. The principle of the Caughey Absorbing Layer Method is first presented (including a rheological interpretation). The efficiency of the method is then shown through 1D Finite Element simulations considering homogeneous and heterogeneous damping in the absorbing layer. 2D models are considered afterwards to assess the efficiency of the absorbing layer method for various wave types and incidences. A comparison with the PML method is first performed for pure P‐waves and the method is shown to be reliable in a more complex 2D case involving various wave types and incidences. It may thus be used for various types of problems involving elastic waves (e.g. machine vibrations, seismic waves, etc.). Copyright © 2010 John Wiley & Sons, Ltd.     

FREQUENCY-INDEPENDENT INFINITE ELEMENTS FOR ANALYSING SEMI-INFINITE PROBLEMS

Two drawbacks exist with the infinite elements used for simulating the unbounded domains of semi-infinite problems. The first is the lack of an adequate measure for calculating the decay parameter. The second is the frequency-dependent characteristic of the finite/infinite element mesh used for deriving the impedance matrices. Based on the properties of wave propagation, a scheme is proposed in this paper for evaluating the decay parameter. In addition, it is shown that by the method of dynamic condensation, the far-field impedance matrices for waves of lower frequencies can be obtained repetitively from the one for waves of the highest frequency, using exactly the same finite/infinite element mesh. Such an approach ensures that accuracy of the same degree can be maintained for waves of all frequencies within the range of consideration. Effectiveness of the proposed method is demonstrated in the numerical examples through comparison with previous results.     

A continued‐fraction‐based high‐order transmitting boundary for wave propagation in unbounded domains of arbitrary geometry

A high‐order local transmitting boundary is developed to model the propagation of elastic waves in unbounded domains. This transmitting boundary is applicable to scalar and vector waves, to unbounded domains of arbitrary geometry and to anisotropic materials. The formulation is based on a continued‐fraction solution of the dynamic‐stiffness matrix of an unbounded domain. The coefficient matrices of the continued fraction are determined recursively from the scaled boundary finite element equation in dynamic stiffness. The solution converges rapidly over the whole frequency range as the order of the continued fraction increases. Using the continued‐fraction solution and introducing auxiliary variables, a high‐order local transmitting boundary is formulated as an equation of motion with symmetric and frequency‐independent coefficient matrices. It can be coupled seamlessly with finite elements. Standard procedures in structural dynamics are directly applicable for evaluating the response in the frequency and time domains. Analytical and numerical examples demonstrate the high rate of convergence and efficiency of this high‐order local transmitting boundary. Copyright © 2007 John Wiley & Sons, Ltd.     

Hermite infinite elements and graded quadratic B-spline finite elements

Quadratic B-spline finite elements are defined for a graded mesh. Hermite infinite elements are proposed to extend the applicability of these finite elements to unbounded regions. Test problems used to compare this technique with published procedures show that the quadratic B-spline finite element solution has, as expected, lower error bounds than a linear element solution. These experiments also demonstrate that the Hermite infinite elements used to close the B-spline finite element arrays lead to error norms comparable in size with other infinite element formulations. The generation of solitary waves in a semi-infinite shallow channel by boundary forcing is modelled by the Korteweg-de Vries equation using an array of graded elements closed by a zero pole infinite element. The resulting simulation of solitary wave motion across a non-uniform mesh confirms existing work and illustrates the effectiveness of the present formulation.     

Corner problems and global accuracy in the boundary element solution of nonlinear wave flows

The numerical model for nonlinear wave propagation in the physical space, developed by Grilli, et al.^12,13, uses a higher-order BEM for solving Laplace's equation, and a higher-order Taylor expansion for integrating in time the two nonlinear free surface boundary conditions. The corners of the fluid domain were modelled by double-nodes with imposition of potential continuity. Nonlinear wave generation, propagation and runup on slopes were successfully studied with this model. In some applications, however, the solution was found to be somewhat inaccurate in the corners and this sometimes led to wave instability after propagation in time.

In this paper, global and local accuracy of the model are improved by using a more stable free surface representation based on quasi-spline elements and an improved corner solution combining the enforcement of compatibility relationships in the double-nodes with an adaptive integration which provides almost arbitrary accuracy in the BEM numerical integrations. These improvements of the model are systematically checked on simple examples with analytical solutions. Effects of accuracy of the numerical integrations, convergence with refined discretization, domain aspect ratio in relation with horizontal and vertical grid steps, are separately assessed. Global accuracy of the computations with the new corner solution is studied by solving nonlinear water wave flows in a two-dimensional numerical wavetank. The optimum relationship between space and time discretization in the model is derived from these computations and expressed as an optimum Courant number of 0.5. Applications with both exact constant shape waves (solitary waves) and overturning waves generated by a piston wavemaker are presented in detail.     

Simulation of crack growth in composite material shell structures

This paper implements a domain integral energy method for modelling crack growth in composite material shell structures using the finite element method. Volume integral expressions to evaluate the dynamic energy release rate in a through‐thickness three‐dimensional crack are derived. Using the domain integral, the energy release rate computation is implemented in the DYNA3D explicit non‐linear dynamic finite element analysis program wherein crack propagation is modelled by releasing the constraints between initially constrained node pairs. The implementation enables the program to either determine the energy resistance response for the material (provided experimental data is available) or predict the rate of crack propagation in shell structures. The numerical implementation was verified by simulating mode I and mode III slow crack growth problems in semi‐infinite transversely isotropic media, for which analytic solutions are available. Oscillations of energy following the release of nodal constraints as the crack propagates in discrete increments were suppressed using light mass proportional damping and a moving averaging scheme. Copyright © 2004 John Wiley & Sons, Ltd.     

Harmonic wave propagation in viscoelastic heterogeneous materials Part I: Dispersion and damping relations

An analytico-numerical method is presented to study the propagation of plane harmonic waves in infinite periodic linear viscoelastic media. Part I considers only the dispersion and attenuation of acoustical longitudinal and shear waves. To show the accuracy of the method, examples of plane harmonic wave propagation in an infinite homogeneous medium and in a periodic layered viscoelastic medium are presented. The method is then used to calculate the damping and dispersion relations for a fibre-reinforced viscoelastic composite material. The results show clearly the influence of materials' viscoelastic properties and heterogeneities on the propagation of plane harmonic waves through the media.     

Transmission of anti-plane shear waves past an interface crack in dissimilar media

The scattered field generated by horizontally polarized incident shear waves in a system of layers having a flaw at the interface is studied. The physical model is idealized to the case of a crack sandwiched between two dissimilar media of infinite height. The local intensification of the dynamic stresses due to the crack is analyzed for incident waves directed at an arbitrary angle. Results of numerical computations are obtained by solving pairs of coupled integral equations and reveal the variations of the stress and displacement fields with ratios of densities and shear moduli of the two adjoining materials.     

Hydroelastic Behavior of a Floating Ring-Shaped Plate

The interaction between incident surface water waves and floating elastic plate is studied. This paper considers the diffraction of plane incident waves on a floating flexible ring-shaped plate and its response to the incident waves. An analytic and numerical study of the hydroelastic behavior of the plate is presented. An integro-differential equation is derived for the problem and an algorithm of its numerical solution is proposed. The representation of the solution as a series of Hankel functions is the key ingredient of the approach. The problem is first formulated. The main integro-differential equation is derived on the basis of the Laplace equation and thin-plate theory. The free-surface elevation, plate deflection and Green’s function are expressed in polar coordinates as superpositions of Hankel and Bessel functions, respectively. These expressions are used in a further analysis of the integro-differential equation. The problem is solved for two cases of water depth infinite and finite. For the coefficients in the case of infinite depth a set of algebraic equations is obtained, yielding an approximate solution. Then a solution is obtained for the general and most interesting case of finite water depth analogously in the seventh section. The exact solution might be approximated by taking into account a finite number of the roots of the plate dispersion relation. Also, the influence of the plate’s motion on wave propagation in the open water field and within the gap of the ring is studied. Numerical results are presented for illustrative purposes.     

Application of 3D time domain boundary element formulation to wave propagation in poroelastic solids

The dynamic responses of fluid-saturated semi-infinite porous continua to transient excitations such as seismic waves or ground vibrations are important in the design of soil-structure systems. Biot's theory of porous media governs the wave propagation in a porous elastic solid infiltrated with fluid. The significant difference to an elastic solid is the appearance of the so-called slow compressional wave. The most powerful methodology to tackle wave propagation in a semi-infinite homogeneous poroelastic domain is the boundary element method (BEM). To model the dynamic behavior of a poroelastic material in the time domain, the time domain fundamental solution is needed. Such solution however does not exist in closed form. The recently developed ‘convolution quadrature method’, proposed by Lubich, utilizes the existing Laplace transformed fundamental solution and makes it possible to work in the time domain. Hence, applying this quadrature formula to the time dependent boundary integral equation, a time-stepping procedure is obtained based only on the Laplace domain fundamental solution and a linear multistep method. Finally, two examples show both the accuracy of the proposed time-stepping procedure and the appearance of the slow compressional wave, additionally to the other waves known from elastodynamics.     

Rayleigh surface waves problem in linear thermoviscoelasticity with voids

In this paper, we study the propagation of the Rayleigh surface waves in a half-space filled by an exponentially functionally graded thermoviscoelastic material with voids. We take into consideration the dissipative character of the porous thermoviscoelastic models upon the propagation waves and study the damped in time wave solutions. The propagation condition is established in the form of an algebraic equation of tenth degree whose coefficients are complex numbers. The eigensolutions of the dynamical system are explicitly obtained in terms of the characteristic solutions. The concerned solution of the Rayleigh surface wave problem is expressed as a linear combination of the five analytical solutions, while the secular equation is established in an implicit form. The explicit secular equation is obtained for an isotropic and homogeneous thermoviscoelastic porous half-space, and some numerical simulations are given for a specific material.     

Wave propagation in micropolar mixture of porous media

This paper is concerned with the possible propagation of waves in an infinite porous continuum consisting of a micropolar elastic solid and a micropolar viscous fluid. Micropolar mixture theory of porous media developed by Eringen [A.C. Eringen, Micropolar mixture theory of porous media, J. Appl. Phys. 94 (2003) 4184–4190] is employed. It is found that there exist four coupled longitudinal waves (two coupled longitudinal displacement waves and two coupled longitudinal microrotational waves) and six coupled transverse waves in a continuum of this micropolar mixture. All the waves are found to attenuate and dispersive in nature. A problem of reflection of coupled longitudinal waves from a free boundary surface of a half-space consisting the mixture of a micropolar elastic solid and Newtonian liquid, is investigated. The expressions of various amplitude ratios and surface responses are derived. Numerical computations are performed to find out the phase velocity and attenuation of the waves. The variation of amplitude ratios, energy ratios and surface responses are also computed for a specific model. All the numerical results are depicted graphically. Some limiting cases have also been discussed.     

A dynamic infinite element for three-dimensional infinite-domain wave problems

A three-dimensional dynamic infinite element which satisfies the following requirements: (1) displacement compatibility on the interface between finite and infinite elements; (2) definition of the wave propagation and amplitude attenuation behaviours in the infinite element using wave propagation functions; (3) convergence of the generalized integrals related to mass and stiffness matrices of the infinite element: and (4) displacement continuity along the common boundary of neighbouring infinite elements in the case of simulating multiple material layers or multiple wave numbers within the foundation, is presented in this paper. Since P-waves, S-waves and R-waves in the foundation can be simulated Simultaneously in the present infinite element, the seismic response of an arch-dam-foundation system, especially a thin double-curvature arch-dam-foundation system where the boundary element loses its competitive capacity with the finite element, can be economically calculated by coupling this infinite element with conventional finite elements. The good accuracy obtained using the present infinite element and finite element coupling model to simulate foundation wave problems has been proven by comparing the current numerical results with previous analytical results.     

Existence of longitudinal waves in anisotropic poroelastic solids

In anisotropic fluid-saturated porous solids, four waves can propagate along a general phase direction. However, solid particles in different waves may not vibrate in mutually orthogonal directions. In the propagation of each of these waves, the displacement of pore–fluid particles may not be parallel to that of solid particles. The polarization for a wave is the direction of aggregate displacement of the particles of the two constituents of a porous aggregate. These polarizations, for different waves, are not mutually orthogonal. Out of the four waves in anisotropic poroelastic medium, two are termed as quasi-longitudinal waves. The prefix ‘quasi’ refers to their polarization being nearly, but not exactly, parallel to the direction of propagation. The existence of purely longitudinal waves in an anisotropic poroelastic medium is ensured by the stationary characters of two expressions. These expressions involve the elastic (stiffness and coupling) coefficients of a porous aggregate and the components of phase direction. Necessary and sufficient conditions for the existence of longitudinal waves are discussed for different anisotropic symmetries. Conditions are also discussed for the existence of the apparent longitudinal waves, i.e., the propagation of wave motion with the particle displacement parallel to the ray direction instead of the phase direction. A graphical solution of a numerical example is shown to check the existence of these apparent longitudinal waves for general directions of phase propagation.     

Elastic Wave Propagation in a Class of Cracked,Functionally Graded Materials by BIEM

Elastic wave propagation in cracked, functionally graded materials (FGM) with elastic parameters that are exponential functions of a single spatial co-ordinate is studied in this work. Conditions of plane strain are assumed to hold as the material is swept by time-harmonic, incident waves. The FGM has a fixed Poisson’s ratio of 0.25, while both shear modulus and density profiles vary proportionally to each other. More specifically, the shear modulus of the FGM is given as μ (x)=μ ₀ exp (2ax ₂), where μ ₀ is a reference value for what is considered to be the isotropic, homogeneous material background. The method of solution is the boundary integral equation method (BIEM), an essential component of which is the Green’s function for the infinite inhomogeneous plane. This solution is derived here in closed-form, along with its spatial derivatives and the asymptotic form for small argument, using functional transformation methods. Finally, a non-hypersingular, traction-type BIEM is developed employing quadratic boundary elements, supplemented with special edge-type elements for handling crack tips. The proposed methodology is first validated against benchmark problems and then used to study wave scattering around a crack in an infinitely extending FGM under incident, time-harmonic pressure (P) and vertically polarized shear (SV) waves. The parametric study demonstrates that both far field displacements and near field stress intensity factors at the crack-tips are sensitive to this type of inhomogeneity, as gauged against results obtained for the reference homogeneous material case     

Diffraction of torsional waves by a penny-shaped crack in an infinitely long cylinder bonded to an infinite medium

This paper contains an analysis of the interaction of torsional waves with penny-shaped crack located in an infinitely long cylinder which is bonded to an infinite medium. Both the cylinder and infinite medium are of homogeneous and elastic but dissimilar materials. The solution of the problem is reduced to a Fredholm integral equation of the second kind which is solved numerically. The numerical solution is used to calculate the stress intensity factor at the rim of the penny-shaped crack.