A coupled electromechanical analysis of a piezoelectric layer bonded to an elastic substrate: Part I,development of governing equations

This two-part contribution presents a novel and efficient method to analyze the two-dimensional (2-D) electromechanical fields of a piezoelectric layer bonded to an elastic substrate, which takes into account the fully coupled electromechanical behavior. In Part I, Hellinger–Reissner variational principle for elasticity is extended to electromechanical problems of the bimaterial, and is utilized to obtain the governing equations for the problems concerned. The 2-D electromechanical field quantities in the piezoelectric layer are expanded in the thickness-coordinate with seven one-dimensional (1-D) unknown functions. Such an expansion satisfies exactly the mechanical equilibrium equations, Gauss law, the constitutive equations, two of the three displacement–strain relations as well as one of the two electric field-electric potential relations. For the substrate the fundamental solutions of a half-plane subjected to a vertical or horizontal concentrated force on the surface are used. Two differential equations and two singular integro-differential equations of four unknown functions, the axial force, N, the moment, M, the average and the first moment of electric displacement, D₀ and D₁, as well as the associated boundary conditions have been derived rigorously from the stationary conditions of Hellinger–Reissner variational functional. In contrast to the thin film/substrate theory that ignores the interfacial normal stress the present one can predict both the interfacial shear and normal stresses, the latter one is believed to control the delamination initiation.     

Radiation from a sphere in compressible plasma

Summary The Hansen's vector wave functions have been modified so it could apply directly for the solution of general exterior boundary value problems in compressible plasma for a spherical geometry. The vector wave functionL has been included to represent an acoustic wave and the three angular orthogonal functions have been defined using complex Fourier series in the azymuthal direction. Using this modified method of vector wave functions, the exterior spherical boundary value problem in compressible isotropic plasma has been solved. The boundary conditions prescribed over the surface of the sphere have been the tangential electric field (or the tangential magnetic field) and the radial component of the velocity vector (or the pressure). From those four basic boundary value problems the coefficients have been derived and several particular cases has been discussed.The research reported in this paper was supported in part by the National Science Foundation, U.S.A.     

On solitary waves. Part 2 A unified perturbation theory for higher-order waves

A unified perturbation theory is developed here for calculating solitary waves of all heights by series expansion of base flow variables in powers of a small base parameter to eighteenth order for the one-parameter family of solutions in exact form, with all the coefficients determined in rational numbers. Comparative studies are pursued to investigate the effects due to changes of base parameters on (i) the accuracy of the theoretically predicted wave properties and (ii) the rate of convergence of perturbation expansion. Two important results are found by comparisons between the theoretical predictions based on a set of parameters separately adopted for expansion in turn. First, the accuracy and the convergence of the perturbation expansions, appraised versus the exact solution provided by an earlier paper [1] as the standard reference, are found to depend, quite sensitively, on changes in base parameter. The resulting variations in the solution are physically displayed in various wave properties with differences found dependent on which property (e.g. the wave amplitude, speed, its profile, excess mass, momentum, and energy), on what range in value of the base, and on the rank of the order n in the expansion being addressed. Secondly, regarding convergence, the present perturbation series is found definitely asymptotic in nature, with the relative error δ(n) (the relative mean-square difference between successive orders n of wave elevations) reaching a minimum, δ _m , at a specific order, n=n _m , both depending on the base adopted, e.g. n _{m , α} =11-12 based on parameter α (wave amplitude), n _{m , β} =15 on β (amplitude-speed square ratio), and n _{m , ∈} =17 on ∈ ( wave number squared). The asymptotic range is brought to completion by the highest order of n=18 reached in this work.     

SH波绕界面孔的散射

用波函数展开方法研究了SH波绕界面孔的散射问题。由入射、反射和透射波组成的自由波场与孔的散射场叠加成总波场。按照一定方式将两个半平面散射波场延拓于全平面，通过Hankel-Fourier展开方法求得了任意形状孔散射场的级数解。以椭圆形孔为例计算了孔边缘的动应力集中系数。     

The mathematical model of reflection and refraction of longitudinal waves in thermo-piezoelectric materials

In this paper, the basic equations of motion, of Gauss and of heat conduction, together with constitutive relations for pyro- and piezoelectric media, are presented. Three thermoelastic theories are considered: classical dynamical coupled theory, the Lord–Shulman theory with one relaxation time and Green and Lindsay theory with two relaxation times. For incident elastic longitudinal, potential electric and thermal waves, referred to as qP, φ-mode and T-mode waves, which impinge upon the interface between two different transversal isotropic media, reflection and refraction coefficients are obtained by solving a set of linear algebraic equations. A case study is investigated: a system formed by two semi-infinite, hexagonal symmetric, pyroelectric–piezoelectric media, namely Cadmium Selenide (CdSe) and Barium Titanate (BaTiO₃). Numerical results for the reflection and refraction coefficients are obtained, and their behavior versus the incidence angle is analyzed. The interaction with the interface give rises to different kinds of reflected and refracted waves: (i) two reflected elastic waves in the first medium, one longitudinal (qP-wave) and the other transversal (qSV-wave), and a similar situation for the refracted waves in the second medium; (ii) two reflected potential electric waves and a similar situation for the refracted waves; (iii) two reflected thermal waves and a similar situation for the refracted waves. The amplitudes of the reflected and refracted waves are functions of the incident angle, of the thermal relaxation times and of the media elastic, electric, thermal constants. This study is relevant to signal processing, sound systems, wireless communications, surface acoustic wave devices and military defense equipment.     

A nth-order shear deformation theory for the free vibration analysis on the isotropic plates

In the present paper, a nth-order shear deformation theory is used to perform the free vibration analysis of the isotropic plates. The present nth-order shear deformation theory satisfies the zero transverse shear stress boundary conditions on the top and bottom surface of the plate. Reddy??s third order theory can be considered as a special case of present nth-order theory (n=3). The governing equations and boundary conditions are derived by the principle of virtual work. The governing differential equations of the isotropic plates are solved by the meshless radial point collocation method based on the thin plate spline radial basis function. The effectiveness of the present theory is demonstrated by applying it to free vibration problem of the square and circular isotropic plate.     

Approximate fundamental solutions and wave fronts for general anisotropic materials

The time-dependent differential equations of elastodynamics for homogeneous solids with a general structure of anisotropy are considered in the paper. A new method of computation of the fundamental solution for these equations is proposed. This method consists of the following. Applying the Fourier transformation with respect to space variables to these equations, we obtain a system of second order ordinary differential equations whose coefficients depend on Fourier parameters. Using the matrix transformations and properties of the coefficients, the Fourier image of the fundamental solution is computed. Finally, the fundamental solution is calculated by the inverse Fourier transformation to the obtained Fourier transform. The implementation and justification of the suggested method have been made by computational experiments in MATLAB. These experiments confirm the robustness of the suggested method. The visualization of the displacement components in general homogeneous anisotropic solids by modern computer tools allows us to see and evaluate the dependence between the structure of solids and the behavior of the displacement field. Our method allows users to observe the elastic wave propagation, arising from pulse point forces of the form e^mδ(x)δ(t), in monoclinic, triclinic and other anisotropic solids. The visualization of displacement components gives knowledge about the form of fronts of elastic wave propagation in Sodium Thiosulfate with monoclinic and Copper Sulphate Pentahydrate with triclinic structures of anisotropy.     

Analysis and computation for nonlinear elastic surface waves of permanent form

Linear elastic surface waves are nondispersive. All wavelengths travel at the Rayleigh wave speed c _R. This absence of frequency dispersion means that nonlinear waves of permanent form cannot be determined as a small perturbation from a sinusoidal wavetrain. By representing the general Rayleigh wave of the linear theory in terms of a pair of conjugate harmonic functions, waves which propagate without distortion are characterized as those having surface elevation profiles which satisfy a certain nonlinear functional equation. In the small-strain limit, this reduces to a quadratic functional equation. Methods for the analysis of this equation are presented for both periodic and nonperiodic waveforms. For periodic waveforms, the infinite system of quadratic equations for the Fourier coefficients of the profile is solved numerically in the case of a certain harmonic elastic material. Two distinct families of profiles having phase speed differing from the linearized Rayleigh wave speed are found. Additionally, two families of exceptional waveforms are found, describing profiles which travel at the Rayleigh wave speed.     

Phragmén-lindelöf type results for harmonic functions with nonlinear boundary conditions

This paper is concerned with investigating the asymptotic behavior of harmonic functions defined on a three-dimensional semi-infinite cylinder, where homogeneous nonlinear boundary conditions are imposed on the lateral surface of the cylinder. Such problems arise in the theory of steady-state heat conduction. The classical Phragmén-Lindelöf theorem states that harmonic functions which vanish on the lateral surface of the cylinder must either grow exponentially or decay exponentially with distance from the finite end of the cylinder. Here we show that the results are significantly different when the homogeneous Dirichlet boundary condition is replaced by the nonlinear heatloss or heat-gain type boundary condition. We show that polynomial growth (or decay) or exponential growth (or decay) may occur, depending on the form of the nonlinearity. Explicit estimates for the growth or decay rates are obtained.     

Free vibration analysis of thin rotating cylindrical shells using wave propagation approach

The wave propagation approach is extended to study the frequency characteristics of thin rotating cylindrical shells. Based on Sanders’ shell theory, the governing equations of motion, which take into account the effects of centrifugal and Coriolis forces as well as the initial hoop tension due to rotation, are derived. And, the displacement field is expressed in the form of wave propagation associated with an axial wavenumber k _m and circumferential wavenumber n. Using the wavenumber of an equivalent beam with similar boundary conditions as the cylindrical shell, the axial wavenumber k _m is determined approximately. Then, the relation between the natural frequency with the axial wavenumber and circumferential wavenumber is established, and the traveling wave frequencies corresponding to a certain rotating speed are calculated numerically. To validate the results, comparisons are carried out with some available results of previous studies, and good agreements are observed. Finally, the relative errors induced by the approximation using the axial wavenumber of an equivalent beam are evaluated with respect to different circumferential wavenumbers, length-to-radius ratios as well as thickness-to-radius ratios, and the conditions under which the analysis presented in this paper will be accurate are discussed.     

Boundary element modeling for defect characterization potential in a wave guide

Wave scattering analysis implemented by boundary element methods (BEM) and the normal mode expansion technique is used to study the sizing potential of two-dimensional shaped defects in a wave guide. Surface breaking half-elliptical shaped defects of three opening lengths (0.3, 6.35 and 12.7 mm) and through-wall depths of 10–90% on a 10 mm thick steel plate were considered. The reflection and transmission coefficients of both Lamb and shear horizontal (SH) waves over a frequency range 0.05–2 MHz were studied. A powerfully practical result was obtained whereby the numerical results for the S₀ mode Lamb wave and n₀ mode SH wave at low frequencies showed a monotonic increase in signal amplitude with an increase in the defect through-wall depth. At high frequency (usually above the cut-off frequency of the A₁ mode for Lamb waves and the n₁ mode for SH waves, respectively), the monotonic trend does not hold in general due to the energy redistribution to the higher order wave modes. Guided waves impinging onto an internal stringer-like an inclusion were also studied. Both the Lamb and SH waves were shown to be insensitive to the stringer internal inclusions at low frequency. Experiments with piezoelectric Lamb wave transducers and non-contact SH wave electro-magnetic acoustic transducers (EMAT) verified some of the theoretical results.     

Saint-Venant End Effects in Anti-plane Shear for Classes of Linear Piezoelectric Materials

This paper is concerned with further investigation of the effect of mechanical/electrical coupling on the decay of Saint-Venant end effects in linear piezoelectricity. Saint-Venant's principle and related results for elasticity theory have received considerable attention in the literature but relatively little is known about analogous issues in piezoelectricity. The current rapidly developing smart structures technology provides motivation for the investigation of such problems. The decay of Saint-Venant end effects is investigated in the context of anti-plane shear deformations for linear homogeneous piezoelectric solids. For a rather general class of anisotropic piezoelectric materials, the governing partial differential equations of equilibrium are a coupled system of second-order partial differential equations for the mechanical displacement u and electric potential ?. The traction boundary-value problem with prescribed surface charge can be formulated as an oblique derivative boundary-value problem for this elliptic system. Energy-decay estimates using differential inequality methods are used to study the axial decay of solutions on a semi-infinite strip subjected to non-zero boundary conditions only at the near end. This analysis is carried out for a rather general class of materials (the tetragonal ${\bar 4}$ crystal class). The boundary-value problem involves a full coupling of mechanical and electrical effects. There are four independent material constants appearing in the problem. An explicit estimated decay rate (a lower bound for the actual decay rate) is obtained in terms of two dimensionless piezoelectric parameters d ₀,r, the first of which provides a measure of the degree of piezoelectric coupling. The estimated decay rate is shown to be monotone decreasing with increasing values of the coupling parameter d ₀. In the limit as d ₀→0, we recover the exact decay rate for the purely mechanical case. Thus, for the tetragonal ${\bar 4}$ class of materials, piezoelectric end effects are predicted to penetrate further into the strip than their elastic counterparts, confirming recent results obtained in other contexts in linear piezoelectricity.     

Scattering of harmonic waves by a disk-shaped rigid inclusion in a three-dimensional elastic matrix

We consider a three-dimensional problem on the interaction of harmonic waves with a thin rigid movable inclusion in an infinite elastic body. The problem is reduced to solving a system of two-dimensional boundary integral equations of Helmholtz potential type for the stress jump functions on the opposite surfaces of the inclusion. We propose a boundary element method for solving the integral equations on the basis of the regularization of their weakly singular kernels. Using the asymptotic relations between the amplitude-frequency characteristics of the wave farzone field and the obtained boundary stress jump functions, we determine the amplitudes of the shear plane wave scattering by a circular disk-shaped inclusion for various directions of the wave incident on the inclusion and for a broad range of wave numbers.     

Surface waves in an exponentially graded, general anisotropic elastic material under the influence of gravity

In a recent paper Destrade [1] studied surface waves in an exponentially graded orthotropic elastic material. He showed that the quartic equation for the Stroh eigenvalue p is, after properly modified, a quadratic equation in p² with real coefficients. He also showed that the displacement and the stress decay at different rates with the depth x₂ of the half-space. Vinh and Seriani [2] considered the same problem and added the influence of gravity on surface waves. In this paper we generalize the problem to exponentially graded general anisotropic elastic materials. We prove that the coefficients of the sextic equation for p remain real and that the different decay rates for the displacement and the stress hold also for general anisotropic materials. A surface wave exists in the graded material under the influence of gravity if a surface wave can propagate in the homogeneous material without the influence of gravity in which the material parameters are taken at the surface of the graded half-space. As the wave number k → ∞, the surface wave speed approaches the surface wave speed for the homogeneous material. A new matrix differential equation for surface waves in an arbitrarily graded anisotropic elastic material under the influence of gravity is presented. Finally we discuss the existence of one-component surface waves in the exponentially graded anisotropic elastic material with or without the influence of gravity.     

Scattering and excitation of SH-surface waves by a protrusion at the mass-loaded boundary of an elastic half-space

A rigorous theory of the scattering and excitation of SH-surface waves by a protrusion at the mass-loaded boundary of an elastic half-space is presented. The boundary value problem (which is of the third kind) is solved by employing two suitably chosen Green functions. One of them is represented as a Fourier type of integral, the other is taken to be the Bessel function of the second kind and order zero. The procedure leads to a system of three, coupled, integral equations. This system is solved numerically. In case of an incident bulk wave, the amplitude of the launched surface wave is computed; in case of an incident surface wave, its transmission and reflection factor are computed. For both cases, an expression for the far-field radiation pattern of the scattered bulk wave is derived. A reciprocity relation is shown to exist between the amplitude of the launched surface wave and the far-field bulk wave radiation pattern. Numerical results are presented for a triangularly-prismatic protrusion; they are compared with the results pertaining to a corresponding indentation in the mass-loaded boundary, that have been obtained in a previous paper.     

Asymptotic crack-tip fields for dynamic crack growth in compressive power-law hardening materials

The effect of material compressibility on the stress and strain fields for a mode-I crack propagating steadily in a power-law hardening material is investigated under plane strain conditions. The plastic deformation of materials is characterized by the J₂ flow theory within the framework of isotropic hardening and infinitesimal displacement gradient. The asymptotic solutions developed by the present authors [Zhu, X.K., Hwang K.C., 2002. Dynamic crack-tip field for tensile cracks propagating in power-law hardening materials. International Journal of Fracture 115, 323–342] for incompressible hardening materials are extended in this work to the compressible hardening materials. The results show that all stresses, strains, and particle velocities in the asymptotic fields are fully continuous and bounded without elastic unloading near the dynamic crack tip. The stress field contains two free parameters σ_eq0 and s₃₃₀ that cannot be determined in the asymptotic analysis, and can be determined from the full-field solutions. For the given values of σ_eq0 and s₃₃₀, all field quantities around the crack tip are determined through numerical integration, and then the effects of the hardening exponent n, the Poisson ratio ν, and the crack growth speed M on the asymptotic fields are studied. The approximate behaviors of the proposed solutions are discussed in the limit of ν → 0.5 or n → ∞.     

Static analysis of anisotropic functionally graded magneto-electro-elastic beams subjected to arbitrary loading

This paper considers the analytical and semi-analytical solutions for anisotropic functionally graded magneto-electro-elastic beams subjected to an arbitrary load, which can be expanded in terms of sinusoidal series. For the generalized plane stress problem, the stress function, electric displacement function and magnetic induction function are assumed to consist of two parts, respectively. One is a product of a trigonometric function of the longitudinal coordinate (x) and an undetermined function of the thickness coordinate (z), and the other a linear polynomial of x with unknown coefficients depending on z. The governing equations satisfied by these z-dependent functions are derived. The analytical expressions of stresses, electric displacements, magnetic induction, axial force, bending moment, shear force, average electric displacement, average magnetic induction, displacements, electric potential and magnetic potential are then deduced, with integral constants determinable from the boundary conditions. The analytical solution is derived for beam with material coefficients varying exponentially along the thickness, while the semi-analytical solution is sought by making use of the sub-layer approximation for beam with an arbitrary variation of material parameters along the thickness. The present analysis is applicable to beams with various boundary conditions at the two ends. Two numerical examples are presented for validation of the theory and illustration of the effects of certain parameters.     

Flexural-gravity wave resistances due to a surface-moving line source on a fluid covered by a thin elastic plate

Analytical solutions for the flexural-gravity wave resistances due to a line source steadily moving on the surface of an infinitely deep fluid are investigated within the framework of the linear potential theory. The homogenous fluid, covered by a thin elastic plate, is assumed to be incompressible and inviscid, and the motion to be irrotational. The solution in integral form for the wave resistance is obtained by means of the Fourier transform and the explicitly analytical solutions are derived with the aid of the residue theorem. The dispersion relation shows that there is a minimal phase speed c_min, a threshold for the existence of the wave resistance. No wave is generated when the moving speed of the source V is less than c_min while the wave resistances firstly increase to their peak values and then decrease when V ? c_min. The effects of the flexural rigidity and the inertia of the plate are studied.     

Scattering of plane monochromatic waves from a heterogeneous inclusion of arbitrary shape in a poroelastic medium: An efficient numerical solution

Scattering of plane longitudinal monochromatic waves from a heterogeneous inclusion of arbitrary shape in an infinite poroelastic medium is considered. Wave propagation in the medium is described by Biot’s equations of poroelasticity. The scattering problem is formulated in terms of the volume integral equations for displacements of the solid skeleton and fluid pressure in the pore space in the region occupied by the inclusion. An efficient numerical method is applied to solve these equations. In the method, Gaussian approximating functions are used for discretization of the problem. For regular node grids, the matrix of the discretized problem has Toeplitz’s properties, and the Fast Fourier Transform technique can be used for the calculation of matrix–vector products. The latter accelerates substantially the process of iterative solution of the discretized problem. For material parameters of typical sedimentary rocks, the system of differential equations of poroelasticity contains a differential operator with a small parameter. As the result, the wave field in the inclusion region is split up into a slowly changing part, and boundary layer functions concentrated near the inclusion interface. The method of matched asymptotic expansions is used for the numerical solution in this case. For a spherical inclusion, the results of the numerical and matched asymptotic expansion methods are compared with a semi-analytical series solution. For a non-spherical heterogeneous inclusion, an example of the numerical solution is presented.     

On the T-matrix for water-wave scattering problems

A rigid cylinder of infinite length is floating in the free surface of deep water. The cylinder is held fixed and a given time-harmonic wave of small amplitude is incident upon it. The corresponding linear two-dimensional boundary-value problem for a velocity potential φ is treated using the null-field method, and an expression for the T-matrix is obtained. (The T-matrix connects the diffraction potential away from the cylinder to the given incident potential.) Fundamental properties of the T-matrix are derived from considerations of energy and reciprocity. For regular wavetrains incident from the right or from the left, there are well-known relations between the corresponding reflection and transmission coefficients; these relations are recovered by specialising the equations satisfied by the T-matrix. Two extensions to water of constant finite depth are described: one uses multipole potentials whilst the other uses Havelock wavemaker functions; this second approach also leads to a new method for treating the problem of waves in a semi-infinite channel with an end-wall of arbitrary shape.