Generalized debye series for acoustic scattering from objects of separable geometric shape

Summary In a classical paper of 1908, Debye has resolved the electromagnetic field scattered by a dielectric cylinder into a series of waves multiply internally reflected in the cylinder. For acoustic scattering by elastic cylinders, a corresponding series was derived from the conventional solution (obtained by satisfying the overall or global continuity conditions) by Brill and Überall, taking into account mode conversions of longitudinal (L) into transverse (T, shear) waves, or vice versa, upon internal scattering in a some-what involved fashion. In a series of papers, Gérard has shown that this approach could be greatly simplified by introducing local reflection and transmission coefficients at each interface, which is suitable for generalizing the Debye series to the case of elastic waves coupled by the continuity conditions at the external and each of any possible (multiple) internal interfaces of the scattering object. The approach is then applicable to all elastic objects for which surface and interfaces form coordinate surfaces of any separable geometry; the corresponding derivation is given here in the most general fashion, and is concretely illustrated by the examples of an elastic plate, infinite cylinder and sphere.     

Improvement of PUFEM for the numerical solution of high‐frequency elastic wave scattering on unstructured triangular mesh grids

The purpose of this paper, which builds on previous work (Int. J. Numer. Meth. Engng 2009; 77 :1646–1669), is to improve a numerical scheme based on the partition of unity finite element method (PUFEM) for the solution of the time harmonic elastic wave equations. The approach consists to approximate the displacement field by the standard finite element shape functions, enriched locally by superimposing pressure (P) and shear (S) plane waves. The aim is to accurately model two‐dimensional elastic wave problems on relatively coarse mesh grids, capable of containing many wavelengths per nodal spacing, for wide ranges of frequencies. This allows us to relax the traditional requirement of about 10 nodal points per S wavelength. In this work, an exact integration scheme for the linear triangular finite element is developed to evaluate the oscillatory integrals arising from the use of the PUFEM. The main contribution here consists in developing an explicit closed‐form solution for two‐dimensional wave‐based integrals, when the phase variation is linear in the local coordinate element system. The evaluation of the element mass matrix is performed from appropriate edge integrals. All other element matrices, obtained by adequate splitting of the element stress tensor matrix, are simply deduced from the element mass matrix entries. The results show clearly that the proposed integration scheme evaluates accurately the entries of the global matrix with drastic reduction of the computational time. Numerical tests dealing with the scattering of S elastic plane waves by a circular rigid body show that, for the same discretization level, it is possible to improve the accuracy by using large elements associated with high numbers of approximating plane waves rather than using small elements with less plane waves. However, this increases the conditioning and the fill‐in of the global matrix. At high frequency, it is even possible to push the number of degrees of freedom per S wavelength under 2 and still achieve good accuracy. Finally, some remarks on the choice of the numbers of P and S plane waves leading to better accuracy and conditioning are discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

Guided waves in functionally graded viscoelastic plates

In this paper, a dynamic solution for the propagating viscoelastic waves in functionally graded material (FGM) plates subjected to stress-free conditions is presented in the context of the Kelvin–Voigt viscoelastic theory. The FGM plate is composed of two orthotropic materials. The material properties are assumed to vary in the thickness direction according to a known variation law. The three obtained wave equations are divided into two groups, which control viscoelastic Lamb-like wave and viscoelastic SH wave, respectively. They are solved respectively by the Legendre orthogonal polynomial series approach. The validity of the method is confirmed through a comparison with the Lamb wave solution of a pure elastic FGM plate and a comparison with the SH wave solution of a viscoelastic homogeneous plate. The dispersion curves and attenuation curves for the graded and homogeneous viscoelastic plates are calculated to highlight their differences. The viscous effect on dispersion curves is shown. The influences of gradient variations are illustrated.     

Generalized thermoelastic waves in functionally graded plates without energy dissipation

In this paper, a dynamic solution of the propagating thermoelastic waves in functionally graded material (FGM) plate subjected to stress-free, isothermal boundary conditions is presented in the context of the Green–Naghdi (GN) generalized thermoelastic theory. The FGM plate is composed of two orthotropic materials. The materials properties are assumed to vary in the direction of the thickness according to a known variation law. The coupled wave equation and heat conduction equation are solved by the Legendre orthogonal polynomial series expansion approach. The convergency of the method is discussed through a numerical example. The dispersion curves of the inhomogeneous thermoelastic plate and the corresponding pure elastic plate are compared to show the characteristics of thermal modes and the influence of the thermoelasticity on elastic modes. The displacement, temperature and stress distributions of elastic modes and thermal modes are shown to discuss their differences. A plate with a different gradient variation is calculated to illustrate the influence of the gradient field on the wave characteristics.     

Effects of inhomogeneity on surface waves in anisotropic media

This paper investigates the effects of anisotropy and inhomogeneity on surface waves in elastic media. Exponential variation in properties are assumed for the elastic parameters and material density. The classical equations of motion for propagation of waves in an inhomogeneous transversely isotropic elastic solid are deduced. The equations of motion for surface waves are derived and general surface waves are investigated. This general theory is then utilized to investigate Rayleigh, Love and Stoneley waves. Results obtained in the above cases reduce to the corresponding well-known classical results when inhomogeneity and anisotropy are not present. It is seen that inhomogeneity has significant effects on dispersion characteristics. Numerical calculations are included for Love waves and some conclusions have been drawn from the above calculations.     

Longitudinal and transverse waves in anisotropic elastic materials

Summary It is shown that a necessary and sufficient condition for a longitudinal wave to propagate in the direction n in an anisotropic elastic material is that the elastic stiffness C ₁₁ (n) is a stationary value (maximum, minimum or saddle point) at n. Explicit expressions of all n and the corresponding elastic stiffness C ₁₁ (n) for which a longitudinal wave can propagate are presented for orthotropic, tetragonal, trigonal, hexagonal and cubic materials. As to longitudinal waves in triclinic and monoclinic materials, only few explicit expressions are possible. We also present necessary and sufficient conditions for a transverse wave to propagate in the direction n. As an illustration, explicit expressions of all n, the polarization vector a and the wave speed c for which a transverse wave can propagate in cubic and hexagonal materials are given. The search for n in hexagonal materials confirms the known fact that a transverse wave can propagate in any direction. A longitudinal wave is necessarily accompanied by two transverse waves. However, a transverse wave can propagate without being accompanied by a longitudinal wave.     

Circular crested waves in anisotropic thermoelastic plates bordered with inviscid liquid

The propagation of circularly crested waves in a homogeneous, transversely isotropic, thermally conducting plate bordered with layers of inviscid liquid or half space of inviscid liquid on both sides is investigated in the context of conventional coupled thermoelasticity, Lord-Shulman and Green-Lindsay theories of thermoelasticity. Secular equations for circular homogeneous transversely isotropic plate in closed form and isolated mathematical conditions for symmetric and antisymmetric wave modes in completely separate terms are derived. The results for isotropic materials and uncoupled theories of thermoelasticity have been obtained as particular cases. The special cases such as short wavelength waves, thin plate waves and leaky Lamb waves of the secular equation are also deduced and discussed. The amplitudes of displacement components and temperature change have also been computed and studied. Finally, the numerical solution is carried out for transversely isotropic circular plate of cobalt material bordered with water. The dispersion curves for symmetric and antisymmetric wave modes, attenuation coefficient and amplitudes of displacement and temperature change in case of fundamental symmetric (S₀) and skew symmetric (A₀) modes are presented in order to illustrate and compare the theoretical results. The analytical and numerical results are found to be in close agreement.     

材料性质对纤维增强复合材料动力学特性的影响

基于弹性动力学理论,对纤维增强复合材料结构中弹性波散射及动应力集中问题进行了分析研究。将非均匀界层区处理为具有相同厚度的多层弹性介质。给出了结构各区域弹性波分析解的表达式。根据位移与应力在界面处的连续条件,确定了弹性波模式系数。给出了界面处动应力集中系数的表达式。作为算例,计算了两种纤维增强复合材料结构中各界面动应力集中系数的数值结果,并分析了界层材料性质以及结构尺寸对各界面动应力集中系数的影响。     

瑞利波作用下径向非匀质地基中的单桩竖向响应研究

在单桩动力响应简化算法研究的基础上，提出了一个用于计算瑞利波激励下的单桩竖向动力响应的简化解析方法，得到了桩底固支时单桩竖向动力位移分布的闭合解。采用考虑桩周土的非匀质性的Novak薄层法计算地基土动力阻抗的方法，获得桩底固支时单桩竖向动力响应的计算公式。分析结果表明非匀质区刚度与匀质区刚度比、非匀质区大小对单桩竖向动力位移响应有一定影响。     

Dynamic Effective Properties of Particle-Reinforced Composites with Viscoelastic Matrix

The frequency-dependent dynamic effective properties of the particle-reinforced composites with the viscoelastic matrix are studied. Several equations to predict the effective wavenumber of the coherent plane waves propagating through particle-reinforced composites are discussed and the equation given by Gubernatis, J.E., [‘Effects of microstructure on speed and attenuation of elastic waves in porous materials’, Wave Motion, 6, 1984, 579–589] based on the independent scattering approximation is used in this paper. The effective phase velocity, the effective attenuation and the effective elastic moduli are evaluated. Numerical calculations are carried out for two kinds of composites, namely, Lead-Epoxy and Glass-Epoxy and the numerical results show that the frequency-dependent dynamic effective properties are related to both the multiple scattering effects among the distributed particles and the viscous dissipative effects of the viscoelastic matrix. However, these effects in the composites with distributed heavy particles (lead) and light particles (glass) are of evidently different features.     

Guided thermoelastic waves in functionally graded plates with two relaxation times

Functionally graded material (FGM) is a promising heat insulation material. Wave propagation in FGM structures has received much attention for the purpose of non-destructive testing and evaluation. Few literatures dealt with the thermoelastic wave in FGM structures although the thermal effect would cause attenuations of elastic waves. In this paper, guided thermoelastic waves in FGM plates subjected to stress-free, isothermal boundary conditions are investigated in the context of the Green–Lindsay (GL) generalized thermoelastic theories (with two relaxation times). Coupled wave equations and heat conduction equation are solved by the Legendre polynomial approach. Dispersion curves for a pure elastic graded plate are calculated to make a comparison with the published data. For the thermoelastic graded plate, dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Attenuation curves for graded plates with different relaxation times are compared. The influences of different material gradient shapes are discussed. Two homogeneous thermoelastic plates with different volume fractions are obtained to show their differences from graded plates. Finally, thermoelastic wave dispersion curves for a homogeneous plate and a graded plate are calculated in the context of the classical coupled thermoelastic theory (CT) to show its differences and similarities to the generalized theory.     

Axisymmetric elasticity solutions for a uniformly loaded annular plate of transversely isotropic functionally graded materials

Summary The axisymmetric problem of a functionally graded, transversely isotropic, annular plate subject to a uniform transverse load is considered. A direct displacement method is developed that the non-zero displacement components are expressed in terms of suitable combinations of power and logarithmic functions of r, the radial coordinate, with coefficients being undetermined functions of z, the axial coordinate. The governing equations as well as the corresponding boundary conditions for the undetermined functions are deduced from the equilibrium equations and the boundary conditions of the annular plate, respectively. Through a step-by-step integration scheme along with the consideration of boundary conditions at the upper and lower surfaces, the z-dependent functions are determined in explicit form, and certain integral constants are then determined completely from the remaining boundary conditions. Thus, analytical elasticity solutions for the plate with different cylindrical boundary conditions are presented. As a promising feature, the developed method is applicable when the five material constants of a transversely isotropic material vary along the thickness arbitrarily and independently. A numerical example is finally given to show the effect of the material inhomogeneity on the elastic field in the annular plate.     

Elastic stress transmission in cellular systems—Analysis of wave propagation

A closely packed array of thin-walled rings constitutes an idealisation of a cellular structure. Elastic waves propagating through such structures must do so via the ring (cell) walls. A theoretical investigation into the propagation of elastic stresses in thin-walled circular rings is undertaken to examine the nature of wave transmission. Three modes of motion, corresponding to shear, extensional and flexural waves, are established and their respective velocities defined by a cubic characteristic equation. The results show that all three waves are dispersive. By neglecting extension of the centroidal axis and rotary inertia, explicit approximate solutions can be obtained for flexural waves. Employment of Love's approach for extensional waves [Love AEH. A treatise on the mathematical theory of elasticity, 4th ed. New York: Dover Publications; 1944. p. 452–3] enables approximate solutions for shear waves to be derived. The three resulting approximate solutions exhibit good agreement with the exact solutions of the characteristic equation over a wide range of wavelengths. The effects of material property, ring wall thickness and ring diameter on the three wave modes are discussed, and the results point to flexural waves as the dominant means of elastic energy transmission in such cellular structures. Wave velocities corresponding to different frequency components determined from experimental results are compared with theoretical predictions of group velocity for flexural waves and good correlation between experimental data and theory affirms this conclusion.     

Shear band systems in plane strain extension: analytical solution and comparison with experimental results

Shear banding represents a local failure mechanism of a soil structure as a response to shear loading. In soil structures of different spatial scales systems of regularly spaced shear bands can be observed as a consequence of extensional loading. The phenomenon of single shear bands, defined as thin zones of localized deformation with a discontinuity of the strain field at its boundaries, is well understood. Inside the shear band the material undergoes inelastic strain softening accompanied by shearing and dilation, whereas the material outside the shear band unloads accompanied by elastic contraction in extension tests. Despite numerous experimental and numerical investigations, the physical mechanisms and parameters determining the spacing of parallel shear bands remained unknown. The paper in hand presents an analytical solution for the spacing of the shear bands and a comparison with a large base of experimental data gained from 1g and ng (geotechnical centrifuge) model experiments. The analytical solution is based on the assumption that the elastic energy rate in the unloaded zone between the shear bands tends to a minimum value. The spacing was calculated as the energetically preferred solution for a broad range of cohesive-frictional granular materials. The dependency of the calculated spacing on initial and boundary conditions as well as on material parameters was found to be in good agreement with the experimental results.     

Perfectly matched layers for transient elastodynamics of unbounded domains

One approach to the numerical solution of a wave equation on an unbounded domain uses a bounded domain surrounded by an absorbing boundary or layer that absorbs waves propagating outward from the bounded domain. A perfectly matched layer (PML) is an unphysical absorbing layer model for linear wave equations that absorbs, almost perfectly, outgoing waves of all non‐tangential angles‐of‐incidence and of all non‐zero frequencies. In a recent work [Computer Methods in Applied Mechanics and Engineering 2003; 192: 1337–1375], the authors presented, inter alia, time‐harmonic governing equations of PMLs for anti‐plane and for plane‐strain motion of (visco‐) elastic media. This paper presents (a) corresponding time‐domain, displacement‐based governing equations of these PMLs and (b) displacement‐based finite element implementations of these equations, suitable for direct transient analysis. The finite element implementation of the anti‐plane PML is found to be symmetric, whereas that of the plane‐strain PML is not. Numerical results are presented for the anti‐plane motion of a semi‐infinite layer on a rigid base, and for the classical soil–structure interaction problems of a rigid strip‐footing on (i) a half‐plane, (ii) a layer on a half‐plane, and (iii) a layer on a rigid base. These results demonstrate the high accuracy achievable by PML models even with small bounded domains. Copyright © 2004 John Wiley & Sons, Ltd.     

Elastic Property Measurement Using Rayleigh-Lamb Waves

Abstract

A nondestructive technique is described for the measurement of elastic constants of isotropic plates using ultrasonic Rayleigh-Lamb waves. The experimental method employs continuous harmonic waves and a pair of variable-angle contact transducers in pitch-catch mode. The phase velocity of the R-L waves at a particular frequency is determined from the phase shift over a measured path length. This simple experimental technique can measure phase velocity over the range 1–10 mm/μs with an error of less than 0.5% over a frequency range of 50 kHz-2 MHz. Individual symmetric and antisymmetric modes can be generated through the selection of transducer angle and frequency. Young's modulus and Poisson's ratio for the material are calculated from measurements of frequency and phase velocity by a nonlinear least squares solution to the dispersion equations. The sensitivity of the nonlinear least squares function to the measurement region of the dispersion curve is investigated. It was found that estimations of material properties are more accurate and less sensitive to small experimental errors when only selected frequencies and R-L modes are used in the least squares calculation. This technique is demonstrated with several isotropic materials and with both thick (6 mm) and thin (0.8 mm) plates. Values for elastic constants determined by the contact transducer Lamb wave technique compare favorably with values measured using the pulse-echo-overlap method. The uncertainty in measurements of Young's modulus and Poisson's ratio was less than 1% and 2%, respectively. The technique has advantages over more traditional methods for measuring elastic properties when it is desirable to use wavelengths greater than the plate thickness, when properties may vary with frequency, or when it is necessary to measure in-plane elastic properties of thin plate structures.     

Thermoelastic Lamb waves in a transversely isotropic plate bordered with layers of inviscid liquid

Analysis for the propagation of thermoelastic waves in a homogeneous, transversely isotropic, thermally conducting plate bordered with layers of inviscid liquid or half space of inviscid liquid on both sides, is investigated in the context of coupled theory of thermoelasticity. Secular equations for homogeneous transversely isotropic plate in closed form and isolated mathematical conditions for symmetric and anti-symmetric wave modes in completely separate terms are derived. The results for isotropic materials and uncoupled theories of thermoelasticity have been obtained as particular cases. It is shown that the purely transverse motion (SH mode), which is not affected by thermal variations, gets decoupled from rest of the motion of wave propagation and occurs along an in-plane axis of symmetry. The special cases, such as short wavelength waves and thin plate waves of the secular equations are also discussed. The secular equations for leaky Lamb waves are also obtained and deduced. The amplitudes of displacement components and temperature change have also been computed and studied. Finally, the numerical solution is carried out for transversely isotropic plate of zinc material bordered with water. The dispersion curves for symmetric and anti-symmetric wave modes, attenuation coefficient and amplitudes of displacement and temperature change in case of fundamental symmetric (S₀) and skew symmetric (A₀) modes are presented in order to illustrate and compare the theoretical results. The theory and numerical computations are found to be in close agreement.     

Analytical solutions describing the consolidation of a multi-layered soil under circular loading

Analytical solutions describing the consolidation of a multi-layered soil under circular loading are presented. From the governing equations of saturated poroelastic soil in a cylindrical coordinate system, the eighth-order state-space equation of consolidation is obtained by eliminating the variation of time t using the Laplace transform together with the technique of Fourier expansions with respect to the coordinate θ and the Hankel transform with respect to coordinate r. The solution of the eighth-order state-space equation is derived directly by using the Laplace transform and its inversion of the z-domain. Based on the continuity between layers and the boundary conditions, the transfer-matrix method is utilized to derive the solutions for the consolidation of a multi-layered soil under circular loading in the transformed domain. By the inversion of the Laplace transform and the Hankel transform, the analytical solutions in the physical domain are obtained. A numerical analysis based on the solutions is carried out by a corresponding program.     

Scaled boundary finite‐element analysis of a non‐homogeneous elastic half‐space

As a result of stresses experienced during and after the deposition phase, a soil strata of uniform material generally exhibits an increase in elastic stiffness with depth. The immediate settlement of foundations on deep soil deposits and the resultant stress state within the soil mass may be most accurately calculated if this increase in stiffness with depth is taken into account. This paper presents an axisymmetric formulation of the scaled boundary finite‐element method and incorporates non‐homogeneous elasticity into the method. The variation of Young's modulus (E) with depth (z) is assumed to take the form E=m_Ez^α, where m_E is a constant and αis the non‐homogeneity parameter. Results are presented and compared to analytical solutions for the settlement profiles of rigid and flexible circular footings on an elastic half‐space, under pure vertical load with αvarying between zero and one, and an example demonstrating the versatility and practicality of the method is also presented. Known analytical solutions are accurately represented and new insight regarding displacement fields in a non‐homogeneous elastic half‐space is gained. Copyright © 2003 John Wiley & Sons, Ltd.     

非均匀复合材料板中剪切波传播的研究

基于弹性界层中弹性波干涉理论,采用有效介质法,研究了剪切波在非均匀、纤维随机分布复合材料板中的传播,得到了非均匀弹性介质内的有效波数。通过满足弹性有界层的上、下边界条件,得到了非均匀界层中的频散方程。作为特例,绘出了不同参数下板中的前四阶频散曲线。可以看出,非均匀弹性有界层中的频散曲线和均匀界层中的有很大不同。最后分析了纤维和基体特性比、纤维的体积份数以及板厚与纤维半径比对频散曲线的影响。