Plane waves in a hypoplastic half-space under periodic boundary conditions

The paper presents a numerical study of the propagation of plane waves in a half-space occupied by a granular material, with periodic boundary conditions for velocity or stresses prescribed at the boundary of the half-space. The constitutive behaviour of the material is described by a simplified hypoplastic equation which takes into account different values of the stiffness for different directions of deformation, and the coupling between shear and volumetric strains owing to dilatancy. These two features are responsible for a nonlinear character of longitudinal waves and for the generation of longitudinal motion by transverse disturbances. It is shown that longitudinal and transverse boundary disturbances produce qualitatively the same longitudinal waves at large distances from the boundary. As a longitudinal wave propagates, the amplitude of oscillations decreases and eventually vanishes, resulting in a single non-oscillating wave.Received: 10 September 2002, Accepted: 31 March 2003 Correspondence to: Y. A. Berezin     

具有表面效应的压电半空间中的表面波

纳米科技的快速发展使压电纳米结构在纳米机电系统中得到广泛应用,形成了诸如纳米压电电子学等新的研究方向.与传统的宏观压电材料相比,在纳米尺度下压电材料往往呈现出不同的力学特性,而造成这种差异的原因之一便是表面效应.本文基于Stroh公式、Barnett-Lothe积分矩阵及表面阻抗矩阵,研究计入表面效应的任意各向异性压电半空间中的表面波传播问题,导出了频散方程.针对横观各向同性压电材料,假设矢状平面平行于材料各向同性面,发现Rayleigh表面波和B-G波解耦,并得到各自的显式频散方程.结果表明,Rayleigh表面波的波速小于偏振方向垂直于表面的体波,而B-G波的波速小于偏振方向垂直于矢状平面的体波.以PZT-5H材料为例,用数值方法考察表面残余应力和电学边界条件对表面波频散特性的影响发现:表面残余应力只对第一阶Rayleigh波有明显的影响;电学开路情形的B-G波比电学闭路情形的B-G波传播快.本文工作可为纳米表面声波器件的设计或压电纳米结构的无损检测提供理论依据.     

Wave propagation at the boundary surface of elastic and initially stressed viscothermoelastic diffusion with voids media

The present investigation is concerned with the wave propagation at the boundary surface of elastic half-space and initially stressed viscothermoelastic diffusion with voids half-space. The longitudinal and transverse waves are incident obliquely at the plane interface between uniform elastic half-space and initially stressed viscothermoelastic diffusion with voids half-space. It is found that the amplitude ratios of various reflected and transmitted waves are functions of angle of incidence, frequency of incident wave and are influenced by the initial stress, diffusion, voids, elastic and viscoelastic properties of media. The expressions of amplitude ratios and energy ratios are obtained in closed form and computed numerically for a specific model. The variations of energy ratios with angle of incidence are shown graphically. The conservation of energy at the interface is verified.     

Free-field dynamic response of an inhomogeneous half-space

In this work, elastic wave propagation in the inhomogeneous half-space is solved by an analytical approach based on plane wave decomposition in conjunction with appropriate functional transformations for the displacement vector. Specifically, free-field motions are recovered at the surface of a half-space with either quadratic or exponential type of depth-dependent material parameters. The incident wave is a time harmonic, planar pressure wave and the resulting free-field motions are obtained in closed form, first for the full-space and then for the half-space by adding the reflected waves. Parametric studies show marked differences in the results when compared against the corresponding ones for a homogeneous background. Finally, sensitivities of the free-field waves on the basic characteristics of the underlying inhomogeneous material and of the incoming wave are investigated.     

Elastic waves in an electro-microelastic solid

The present study is concerned with the wave propagation in an electro-microelastic solid. The reflection phenomenon of plane elastic waves from a stress free plane boundary of an electro-microelastic solid half-space is studied. The condition and the range of frequency for the existence of elastic waves in an infinite electro-microelastic body are investigated. The constitutive relations and the field equations for an electro-microelastic solid are stemmed from the Eringen’s theory of microstretch elasticity with electromagnetic interactions. Amplitude ratios and energy ratios of various reflected waves are presented when an elastic wave is made incident obliquely at the stress free plane boundary of an electro-microelastic solid half-space. It has been verified that there is no dissipation of energy at the boundary surface during reflection. Numerical computations are performed for a specific model to calculate the phase speeds, amplitude ratios and energy ratios, and the results obtained are depicted graphically. The effect of elastic parameter corresponding to micro-stretch is noticed on reflection coefficients, in particular. Results of Parfitt and Eringen [Parfitt, V.R., Eringen, A.C., 1969. Reflection of plane waves from a flat boundary of a micropolar elastic half-space. J. Acoust. Soc. Am. 45, 1258–1272] have also been reduced as a special case from the present formulation.     

Scattering of pulsed Rayleigh surface waves by a cylindrical cavity

A pulsed Rayleigh surface wave of prescribed shape is incident on a cylindrical cavity which is parallel to both the plane free surface and the plane wave front. Multiple reflections at the cylindrical and plane free surface are considered and the resulting displacements and stress components are calculated in the surrounding of the cavity by approximately summing infinite double sums. Use is made of the stationary loading case simulated by a periodic train of wave pulses and its time Fourier series representation and of expansions of all incident and reflected waves in terms of cylindrical wave functions. For reflection, the free surface of the half-space is approximated by a fictitious convex (or concave) cylindrical surface of “large” radius. The wave pattern due to a single pulse loading is constructed from the stationary solution by enforcing homogeneous initial conditions in the half-space ahead of the single loading pulse and by prescribing a wide spacing in the periodically set-forth train of pulses. The numerical results for stresses and dynamic stress magnification factors are especially useful for the interpretation of recent measurements in dynamic photoelasticity.     

Inhomogeneous waves at the boundary of a generalized thermoelastic anisotropic medium

The Christoffel equation is derived for the propagation of plane harmonic waves in a generalized thermoelastic anisotropic (GTA) medium. Solving this equation for velocities implies the propagation of four attenuating waves in the medium. The same Christoffel equation is solved into a polynomial equation of degree eight. The roots of this equation define the vertical slownesses of the eight attenuating waves existing at a boundary of the medium. Incidence of inhomogeneous waves is considered at the boundary of the medium. A finite non-dimensional parameter defines the inhomogeneity of incident wave and is used to calculate its (complex) slowness vector. The reflected attenuating waves are identified with the values of vertical slowness. Procedure is explained to calculate the slowness vectors of the waves reflected from the boundary of the medium. The slowness vectors are used, further, to calculate the phase velocities, phase directions, directions and amounts of attenuations of the reflected waves. Numerical examples are considered to analyze the variations of these propagation characteristics with the inhomogeneity and propagation direction of incident wave. Incidence of each of the four types of waves is considered. Numerical example is also considered to study the propagation and attenuation of inhomogeneous waves in the unbounded medium.     

Rapid Sliding Indentation with Friction of a Pre-Stressed Thermoelastic Material

A rigid indentor travels with a constant speed over the surface of an isotropic thermoelastic half-space. Friction exists between the indentor and half-space, and the latter is initially in equilibrium at a uniform temperature under a uniform normal pre-stress. This pre-stress, below but near yield, is assumed to produce deformations that dominate the additional deformations due to indentation. Thus, the process is treated as one of small deformations superposed upon (relatively) large. The governing equations for the superposed deformation are those of nonisotropic coupled thermoelasticity. A steady-state two-dimensional study uses robust asymptotic analytical solutions to reduce the associated mixed boundary value problem to a classical singular integral equation which can be solved analytically. The solutions show that the pre-stress-induced de facto nonisotropy alters the values of the rotational and dilatational wave and Rayleigh speeds in the half-space and, in the case of a compressive pre-stress, generates a second, lower, critical speed. In addition, pre-stress generates noncritical sliding speeds at which the friction-dependent integral equation eigenvalue changes sign. For purposes of illustration, expressions for the half-space surface temperature change and its average over the contact zone, the equations necessary to determine contact zone size and location, the resultant moment on the indentor, and the maximum compressive stress on the contact zone are presented for a parabolic indentor. This revised version was published online in August 2006 with corrections to the Cover Date.     

Seismic Reflection from an Interface Between an Elastic Solid and a Fractured Porous Medium with Partial Saturation

The theory of Tuncay and Corapcioglu (Transp Porous Media 23:237–258, 1996a) has been employed to investigate the possibility of plane wave propagation in a fractured porous medium containing two immiscible fluids. Solid phase of the porous medium is assumed to be linearly elastic, isotropic and the fractures are assumed to be distributed isotropically throughout the medium. It has been shown that there can exist four compressional waves and one rotational wave. The phase speeds of these waves are found to be affected by the presence of fractures, in general. Of the four compressional waves, one arises due to the presence of fractures in the medium and the remaining three are those encountered by Tuncay and Corapcioglu (J Appl Mech 64:313–319, 1997). Reflection and transmission phenomena at a plane interface between a uniform elastic half-space and a fractured porous half-space containing two immiscible fluids, are analyzed due to incidence of plane longitudinal/transverse wave from uniform elastic half-space. Variation of modulus of amplitude and energy ratios with the angle of incidence are computed numerically by taking the elastic half-space as granite and the fractured porous half-space as sandstone material containing non-viscous wetting and non-wetting fluid phases. The results obtained in case of porous half-space with fractures, are compared graphically with those in case of porous half-space without fractures. It is found that the presence of fractures in the porous half-space do affect the reflection/transmission of waves, which is responsible for raising the reflection and lowering the transmission coefficients.     

Reflection and refraction of attenuated waves at boundary of elastic solid and porous solid saturated with two immiscible viscous fluids

The propagation of elastic waves is studied in a porous solid saturated with two immiscible viscous fluids.The propagation of three longitudinal waves is represented through three scalar potential functions.The lone transverse wave is presented by a vector potential function.The displacements of particles in different phases of the aggregate are defined in terms of these potential functions.It is shown that there exist three longitudinal waves and one transverse wave.The phenomena of reflection and refraction due to longitudinal and transverse waves at a plane interface between an elastic solid half-space and a porous solid half-space saturated with two immiscible viscous fluids are investigated.For the presence of viscosity in pore-fluids,the waves refracted to the porous medium attenuate in the direction normal to the interface.The ratios of the amplitudes of the reflected and refracted waves to that of the incident wave are calculated as a nonsingular system of linear algebraic equations.These amplitude ratios are used to further calculate the shares of different scattered waves in the energy of the incident wave.The modulus of the amplitude and the energy ratios with the angle of incidence are computed for a particular numerical model.The conservation of the energy across the interface is verified.The effects of variations in non-wet saturation of pores and frequencies on the energy partition are depicted graphically and discussed.     

Finite-amplitude Love waves in a pre-stressed compressible elastic half-space with a double surface layer

The dispersive behavior of finite-amplitude time-harmonic Love waves propagating in a pre-stressed compressible elastic half-space overlaid with two compressible elastic surface layers of finite thickness is investigated. The half-space and layers are made of different pre-stressed compressible neo-Hookean materials. The dispersion relation which relates wave speed and wavenumber is obtained in explicit form. Results for the energy density and energy flux of the waves are also presented. The special case where the interfaces between the layers and the half-space are principal planes of the left Cauchy–Green deformation tensor is also investigated. Numerical results are presented showing the variation of the Love wave speed with the pre-stress and the propagation angle.     

Torsional surface waves in heterogeneous anisotropic half-space under initial stress

The present paper is concerned with the propagation of torsional surface waves in a heterogeneous anisotropic half-space under the initial compressive stress. The heterogeneity in the half-space is caused by the linear variation in rigidity, initial compressive stress and density. The solution part of the problem involves the use of Whittaker function. The dispersion equation has been obtained in a closed form, which shows the variation of phase velocity with corresponding wave number. Effects of anisotropy and initial stress have been shown by the means of graphs for different anisotropic materials. It has found that the phase velocity of torsional waves decreases with increment in initial stress and inhomogeneity. Obtained phase velocity of torsional surface wave is found to be less than the shear wave velocity, which agrees with the standard result.     

Reflection and transmission of elastic waves at an elastic/porous solid saturated by two immiscible fluids

Wave propagation in a porous elastic medium saturated by two immiscible fluids is investigated. It is shown that there exist three dilatational waves and one transverse wave propagating with different velocities. It is found that the velocities of all the three longitudinal waves are influenced by the capillary pressure, while the velocity of transverse wave does not at all. The problem of reflection and refraction phenomena due to longitudinal and transverse wave incident obliquely at a plane interface between uniform elastic solid half-space and porous elastic half-space saturated by two immiscible fluids has been analyzed. The amplitude ratios of various reflected and refracted waves are found to be continuous functions of the angle of incidence. Expression of energy ratios of various reflected and refracted waves are derived in closed form. The amplitude ratios and energy ratios have been computed numerically for a particular model and the results obtained are depicted graphically. It is verified that during transmission there is no dissipation of energy at the interface. Some particular cases have also been reduced from the present formulation.     

Longitudinal waves at a micropolar fluid/solid interface

The possibility of plane wave propagation in a micropolar fluid of infinite extent has been explored. The reflection and transmission of longitudinal elastic wave at a plane interface between a homogeneous micropolar fluid half-space and a micropolar solid half-space has also been investigated. It is found that there can exist four plane waves propagating with distinct phase speeds in an infinite micropolar fluid. All the four waves are found to be dispersive and attenuated. The reflection and transmission coefficients are found to be the functions of the angle of incidence, the elastic properties of the half-spaces and the frequency of the incident wave. The expressions of energy ratios have also been obtained in explicit form. Frequency equation for the Stoneley wave at micropolar solid/fluid interface has also been derived in the form of sixth-order determinantal expression, which is found in full agreement with the corresponding result of inviscid liquid/elastic solid interface. Numerical computations have been performed for a specific model. The dispersion curves and attenuation of the existed waves in micropolar fluid have been computed and depicted graphically. The variations of various amplitudes and energy ratios are also shown against the angle of incidence. Results of some earlier workers have been deduced from the present formulation.     

Wave reflection and refraction in triclinic crystalline media

Summary In this paper, the reflection and refraction of a plane wave at an interface between two half-spaces composed of triclinic crystalline material is considered. It is shown that due to incidence of plane wave three types of waves, namely quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated governed by the propagation condition involving the acoustic tensor. A simple procedure has been presented for the calculation of all the three phase velocities of the quasi waves. It has been established that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation. Relations are established between directions of motion and propagation, respectively. The expressions for reflection and refraction coefficients of qP, qSV and qSH waves are obtained. Numerical results of reflection and refraction coefficients are presented for different types of anisotropic media and for different types of incident waves. Graphical representations have been made for incident qP waves, and for incident qSV and qSH waves numerical data are presented in tables.The work was completed while the author was visiting the University of Kaiserslautern, Department of Geomathematics as Visiting Professor. The Author is grateful to Professor Dr. W. Freeden for providing DAAD fellowship and all the facilities for conducting research, as well as to Dr. V.Michel for various discussions about the research work and also for all kinds of help during his stay at Kaiserslautern, Germany. This award is very gratefully acknowledged.     

Plane strain dynamics of elastic solids with intrinsic boundary elasticity, with application to surface wave propagation

In this paper, in a development of the static theory derived by Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853), we establish the equations of motion for a non-linearly elastic body in plane strain with an elastic surface coating on part or all of its boundary. The equations of (linearized) incremental motions superposed on a finite static deformation are then obtained and applied to the problem of (time-harmonic) surface wave propagation on a pre-stressed incompressible isotropic elastic half-space with a thin coating on its plane boundary. The secular equation for (dispersive) wave speeds is then obtained in respect of a general form of incompressible isotropic elastic strain-energy function for the bulk material and a general energy function for the coating material. Specialization of the form of strain-energy function enables the secular equation to be cast as a quartic equation and we therefore focus on this for illustrative purposes. An explicit form for the secular equation is thereby obtained. This involves a number of material parameters, including residual stress and moment in the properties of the coating. It is shown how this equation relates to previous work on waves in a half-space with an overlying thin layer set in the classical theory of isotropic elasticity and, in particular, the significant effect of omission of the rotatory inertia term, even at small wave numbers, is emphasized. Corresponding results for a membrane-type coating, for which the bending moment, inertia and residual moment terms are absent, are also obtained. Asymptotic formulas for the wave speed at large wave number (high frequency) are derived and it is shown how these results influence the character of the wave speed throughout the range of wave number values. A bifurcation criterion is obtained from the secular equation by setting the wave speed to zero, thereby generalizing the bifurcation results of Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853) to the situation in which residual stress and moment are present in the coating. Numerical results which show the dependence of the wave speed on the various material parameters and the finite deformation are then described graphically. In particular, features which differ from those arising in the classical theory are highlighted.     

Effect of rigid boundary on propagation of torsional surface waves in porous elastic layer

The paper presents the effect of a rigid boundary on the propagation of torsional surface waves in a porous elastic layer over a porous elastic half-space using the mechanics of the medium derived by Cowin and Nunziato (Cowin, S. C. and Nunziato, J. W. Linear elastic materials with voids. Journal of Elasticity, 13(2), 125–147 (1983)). The velocity equation is derived, and the results are discussed. It is observed that there may be two torsional surface wave fronts in the medium whereas three wave fronts of torsional surface waves in the absence of the rigid boundary plane given by Dey et al. (Dey, S., Gupta, S., Gupta, A. K., Kar, S. K., and De, P. K. Propagation of torsional surface waves in an elastic layer with void pores over an elastic half-space with void pores. Tamkang Journal of Science and Engineering, 6(4), 241–249 (2003)). The results also reveal that in the porous layer, the Love wave is also available along with the torsional surface waves. It is remarkable that the phase speed of the Love wave in a porous layer with a rigid surface is different from that in a porous layer with a free surface. The torsional waves are observed to be dispersive in nature, and the velocity decreases as the oscillation frequency increases.     

Visual observations of the flow around a half-submerged oscillating circular cylinder

Visual observations are made on the flow around a horizontal circular cylinder which is half-submerged in still water and forced to oscillate vertically. The ends of the cylinder have great influence on the wave pattern and flow field. Progressive plane waves are generated at small forcing amplitudes, but cross-waves are superimposed on the progressive plane waves at large forcing amplitudes. The wavelength of the cross-waves in the direction parallel to the cylinder axis increases with the forcing amplitude. The crests of the cross-waves are in parallel lines which are oblique to the cylinder axis. The angle at which the parallel lines meet the cylinder axis decreases as the forcing amplitude is increased. Three kinds of steady flows are induced in the water: surface flow, undersurface flow, and vertical jet.     

Inhomogeneous plane waves in compressible viscous fluids

The propagation of inhomogeneous plane waves in a compressible viscous fluid is considered. The frequency and the slowness vector are both allowed to be complex. There are seen to be two types of solutions: (a) two transverse waves, which involve no density or pressure fluctuations, (b) a longitudinal wave, which involves no fluctuations in vorticity. For each type, a propagation condition is obtained giving the (complex) squared length of the slowness vector as a function of frequency. Each depends also on the viscosities. It is seen how to recover the incompressible case as the limit in which the inviscid acoustic wave speed tends to infinity. Each wave is shown to be linearly stable for real frequencies. These waves are attenuated in space and time but nevertheless it is possible to define constant weighted mean values (over a cycle of the propagating part of the wave) of the energy density, energy flux and dissipation. The energy-dissipation equation and the propagation conditions are used to derive relationships between these constant weighted means, some of which are generalizations to compressible fluids of previously known results for incompressible fluids. Explicit expressions in terms of frequency are given for the weighted means.