On formulas for the Rayleigh wave velocity in pre-stressed compressible solids

In this paper, formulas for the velocity of Rayleigh waves in compressible isotropic solids subject to uniform initial deformations are derived using the theory of cubic equation. They are explicit, have simple algebraic forms, and hold for a general strain energy function. Unlike the previous investigations where the derived formulas for Rayleigh wave velocity are approximate and valid for only small enough values of pre-strains, this paper establishes exact formulas for Rayleigh wave velocity being valid for any range of pre-strains. When the prestresses are absent, the obtained formulas recover the Rayleigh wave velocity formula for compressible elastic solids. Since obtained formulas are explicit, exact and hold for any range of pre-strains, they are good tools for evaluating nondestructively prestresses of structures.     

An approach for obtaining approximate formulas for the Rayleigh wave velocity

In this paper, we introduce an approach for finding analytical approximate formulas for the Rayleigh wave velocity for isotropic elastic solids and anisotropic elastic media as well. The approach is based on the least-square principle. To demonstrate its application, we applied it in order to obtain an explanation for Bergmann’s approximation, the earliest known approximation of the Rayleigh wave velocity for isotropic elastic solids, and used it to establish a new approximation. By employing this approach, the best approximate polynomials of the second order of the cubic power and the quartic power in the interval [0, 1] were found. By using the best approximate polynomial of the second order of the cubic power, we derived an approximate formula for the Rayleigh wave speed in isotropic elastic solids which is slightly better than the one given recently by Rahman and Michelitsch by employing Lanczos’s approximation. Also by using this second order polynomial, analytical approximate expressions for orthotropic, incompressible and compressible elastic solids were found. For incompressible case, it is shown that the approximation is comparable with Rahman and Michelitsch’s approximation, while for the compressible case, it is shown that our approximate formulas are more accurate than Mozhaev’s ones. Remarkably, by using the best approximate polynomials of the second order of the cubic power and the quartic power in the interval [0, 1], we derived an approximate formula of the Rayleigh wave velocity in incompressible monoclinic materials, where the explicit exact formulas of the Rayleigh wave velocity so far are not available.     

Rayleigh wave propagation in curved waveguides

A JWKB asymptotic expansion describing inplane elastic waves is used to approximate a Rayleigh-like wave guided within a curved elastic waveguide whose curvature is small and changes slowly over a wavelength. The two lowest eigenmodes in a curved guide, taken together, constitute the Rayleigh-like wave. It is shown that this wave lies in the shadows of four, closely spaced, virtual caustics, two caustics per constituent eigenmode. If the curvature becomes too large one or more of the caustics ceases to be virtual and enters the guide after which a Rayleigh-like wave cannot propagate. The overall disturbance is shown to have an amplitude that is modulated because the wavenumbers of the constituent eigenmodes differ by a small amount. Moreover, the disturbance is shown to propagate with a wavenumber that, to leading order, has a linear dependence on the curvature causing the phase to be modulated, as well. Passing from a thin guide to a very thick one suppresses the amplitude modulation, making the phase modulation evident. Propagation into an environment of increasing curvature, for both thin and thick, shallowly curved guides is studied so that the modulations may be observed.     

Implicit and explicit secular equations for Rayleigh waves in two-dimensional anisotropic media

This paper is concerned with the derivation of implicit and explicit secular equations for Rayleigh waves polarized in a plane of symmetry of an anisotropic linear elastic medium. It has been confirmed, in accord with Ting’s paper [2], that the Rayleigh waves propagate with no geometric dispersion. Numerical evaluations of both the implicit and explicit equations give the same values of Rayleigh wave velocities. In the case of orthotropic material (thin composites) it has been found that Rayleigh wave velocity depends significantly (as with bulk waves) on the directions of principal material axes. For the same material model the analytical solutions, based on implicit and explicit secular equations, were compared against the finite element and experimental data that had been published by Cerv et al. [4] in 2010. It emerged that the theory was in accordance with the experiment.     

An analytical study on the Rayleigh wave generation in a stratified structure

In this paper, an analytical method is used to investigate the Rayleigh wave generation in a stratified structure and the wave generation in a dry sandy layer constrained between the couple stress and inhomogeneous orthotropic half-spaces. This study is devoted to analyzing the impact of various effective parameters associated with the media on the phase velocities of the wave. The displacement components for each medium are derived by implementing the separable variable method. The frequency eq...     

Growing crack as superposition of surface Rayleigh waves

A crack, symmetrically propagating in elastic material, was considered as superposition of surface Rayleigh waves. The self-similar growth of face loaded crack was considered in detail. Exact expressions of deformation and stress fields in the crack’s surrounding were found and asymptotic behavior of stress near crack’s tips was also obtained. A condition that determines the crack’s velocity of self-similar propagation was found.     

一种求解瑞利波散射问题的修正边界元方法

边界元法的一大优势是用于求解半空间等无限域问题,但在求解弹性波的传播问题中,传统边界元法采用的全平面或全空间格林函数会在截断边界处产生虚假的反射回波,这会引起散射场求解的误差.为了避免这种误差,论文在传统边界元法基础上提出一种修正边界元法,该修正方法主要包括:以瑞利波形式的远端散射场代替原本因截断而舍去的部分,通过互易定理建立单位瑞利波和全平面格林函数的积分方程,求得修正系数,并代入修正边界元矩阵,计算出瑞利波的散射场.文中基于该方法计算了无缺陷平面的瑞利波(与解析解的误差为1.24×10-5),并运用该方法计算了不同缺陷的散射场.由文中对比结果表明,论文所提修正边界元法可准确求解瑞利波散射场,为基于表面波的缺陷反演问题研究提供了有效的正演途径.     

超声瑞利波测试材料表面载荷应力的实验研究

实验研究了利用超声瑞利波传播特性测试材料表面应力的方法.在测试过程中采用回振法测声速,以分辨应力作用引起声速的徽小变化;并设计夹具实现探头与试样之间的稳定耦合和弹性接触,精简了实验机构,提高了实验精度.研究结果表明,在应力作用下,平行和垂直于应力方向传播的瑞利波发生了相应的变化:应力较小时,声速随应力变化较快;随着应力...     

Accurate 3-D finite element simulation of elastic wave propagation with the combination of explicit and implicit time-integration methods

One of the big issues in finite element solutions of wave propagation problems is the presence of spurious high-frequency oscillations that may lead to divergent results at mesh refinement. The paper deals with the extension of the new two-stage time-integration technique developed in our previous papers to the solution of wave propagation problems with explicit time-integration methods.The explicit central difference method is used for accurate time-integration of the semi-discrete system of elastodynamics at the stage of basic computations and allows spurious high-frequency oscillations. To filter these oscillations, pre- or/and post-processing (the filtering stage) is applied using a few time increments of the implicit time-continuous Galerkin method with large numerical dissipation.A special calibration procedure is used for the selection of the minimum necessary amount of numerical dissipation (in terms of a time increment) at the filtering stage. In contrast to existing approaches that use a time-integration method with the same dissipation (or artificial viscosity) for all time increments, the new technique yields accurate and non-oscillatory results for wave propagation problems without interaction between user and computer code. The solutions of 3-D wave propagation and impact problems show the effectiveness of the new approach.     

Reflection and transmission of obliquely incident Rayleigh surface waves by a viscoelastic interphase

The reflection and transmission of obliquely incident Rayleight surface waves by an interphase between two quarter spaces of identical or different materials, have been investigated. The mechanical behavior of the interphase is represented by a thin viscoelastic layer. By using the full space Green's functions due to a spatially harmonic line load, the mathematical statement of the 3-dimension problem is reduced to a 2-dimension system of singular integral equations. The far-field behavior of the scattered waves leads to the definition of reflection and transmission coefficients,R andT. The system of the singular integral equations are solved forR andT with the boundary element method. The results are presented for selected values of the elastic constants of the joined quarter spaces, the parameters of the interphase and the incident angles of Rayleigh surface waves.     

Perturbation Formula for Phase Velocity of Rayleigh Waves in Prestressed Anisotropic Media

Herein we consider Rayleigh waves propagating along the free surface of a macroscopically homogeneous, anisotropic, prestressed half-space. We adopt the formulation of linear elasticity with initial stress and assume that the deviation of the prestressed anisotropic medium from a comparative ‘unperturbed’, unstressed and isotropic state, as formally caused by the initial stress and by the anisotropic part of the incremental elasticity tensor, be small. No assumption, however, is made on the material anisotropy of the incremental elasticity tensor. With the help of the Stroh formalism, we derive a first-order perturbation formula for the shift of phase velocity of Rayleigh waves from its comparative isotropic value. Our perturbation formula does not agree totally with that which was derived some years ago by Delsanto and Clark, and we provide another argument as further support for our version of the formula. According to our first-order formula, the anisotropy-induced velocity shifts of Rayleigh waves, taken in totality of all propagation directions on the free surface, carry information only on 13 elastic constants of the anisotropic part of the incremental elasticity tensor. The remaining eight elastic constants are those which would become zero if were monoclinic with the two-fold symmetry axis normal to the free surface of the material half-space in question.     

On the Rayleigh Wave Speed in Orthotropic Elastic Solids

Recently, a formula for the Rayleigh wave speed in an isotropic elastic half-space has been given by Malischewsky and a detailed derivation given by the present authors. This study deals with the generalization of this formula to orthotropic elastic materials and Malischewsky’s formula is recovered as a special case. The formula is obtained using the theory of cubic equations and is expressed as a continuous function of three dimensionless material parameters.     

Theoretical study of the generation of screw dislocations for Rayleigh and Lamb waves in isotropic solids

Phase singularities are generic structures which occur in all wave fields, and they are characterised by an inability to assign a value to the phase. Screw dislocations are a particular kind of phase singularity where the phase possesses a helical structure, with the singularity at the centre of the helix. In this paper we show that it is possible to generate screw dislocations on the surface of elastic isotropic solids by means of the interference of three Rayleigh waves or three Lamb waves. The dispersive character of Lamb waves leads to more complicated behaviour, which may in turn result in greater potential for applications.     

Dispersion and localisation in structured Rayleigh beams

This paper brings a comparative analysis between dynamic models of couple-stress elastic materials and structured Rayleigh beams on a Winkler foundation. Although physical phenomena have different physical origins, the underlying equations appear to be similar, and hence mathematical models have a lot in common. In the present work, our main focus is on the analysis of dispersive waves, band-gaps and localised waveforms in structured Rayleigh beams. The Rayleigh beam theory includes the effects of rotational inertia which are neglected in the Euler–Bernoulli beam theory. This makes the approach applicable to higher frequency regimes. Special attention is given to waves in pre-stressed Rayleigh beams on elastic foundations.     

界面凹坑上纵波－瑞利波转换的实验研究

应用超声技术对界面半球形凹坑上的纵波－瑞利波的波型转换进行了实验研究。在一个特制的钢质模型上，频率为ＩＭＨｚ的纵波从６个不同方向依次入射到３个不同直径的表面凹坑上，然后在８个不同的方向上接收瑞利波。所接收的散射信号被数字化后进行了ＦＦＴ运算。根据所得实验数据得知，当纵波垂直入射凹坑时所转换的瑞利波的总能量最大；当凹坑直径与入射波长相等时波型转换率最大等结论。其中一些结果为反卷积运算和已知数值计算所验证。     

Dispersion and Attenuation of Surface Waves in a Fluid-Saturated Porous Medium

An investigation is conducted of propagation of surface waves in a porous medium consisting of a microscopically incompressible solid skeleton in which a microscopically incompressible liquid flows within the interconnected pores, and particularly the case where the solid skeleton deforms linear elastically. The frequency equations of Rayleigh- and Love-type waves are derived relating the dependence of wave numbers, being complex quantities, on frequency, as a result those waves are dispersive as well as inhomogeneous. Nevertheless, the amplitudes of both surface waves attenuate along the surface of the porous medium, whereas they decay exponentially receding from the surface of the medium.     

Propagation of Rayleigh waves on curved surfaces

The propagation and properties of Rayleigh waves on curved surfaces are investigated theoretically. The Rayleigh wave dispersion equation for propagation on a curved surface is derived as a parabolic equation, and its penetration depth is analyzed using the curved surface boundary. Reciprocity is introduced to model the diffracted Rayleigh wave beams. Simulations of Rayleigh waves on some canonical curved surfaces are carried out, and the results are used to quantify the influence of curvature. It is found that the velocity of the surface wave increases with greater concave surface curvature, and a Rayleigh wave no longer exists once the surface wave velocity exceeds the bulk shear wave velocity. Moreover, the predicted wave penetration depth indicates that the energy in the Rayleigh wave is transferred to other modes and cannot propagate on convex surfaces with large curvature. A strong directional dependence is observed for the propagation of Rayleigh waves in different directions on surfaces with complex curvatures. Thus, it is important to include dispersion effects when considering Rayleigh wave propagation on curved surfaces.     

A novel method for investigation of acoustic and elastic wave phenomena using numerical experiments

The emergence of new types of composite materials, the depletion of existing hydrocarbon deposits, and the increase in the speed of trains require the development of new research methods based on wave scattering. Therefore, it is necessary to determine the laws of wave scattering in inhomogeneous media. We propose a method that combines the advantages of a numerical simulation with an analytical study of the boundary value problem of elastic and acoustic wave equations. In this letter we present the results of the study using the proposed method: the formation of a response from a shear wave in an acoustic medium and the formation of shear waves when a vertically incident longitudinal wave is scattered by a vertical gas-filled fracture. We have obtained a number of analytical expressions characterising the scattering of these wave types.     

Macroscopic and microscopic examination of the relationship between crack velocity and path and Rayleigh surface wave speed in single crystal silicon

The speed of Rayleigh surface waves, denoted C_R, is the accepted upper limit for Mode I crack velocity in monolithic solids. In the current contribution, we discuss several critical issues associated with the velocity of Rayleigh surface waves and crack velocity in single crystal (SC) brittle solids, and the global and local influence of C_R on crack path selection in particular.Recent cleavage experiments in SC silicon showed that crack velocity at certain cleavage planes and crystallographic orientations cannot exceed a small fraction of C_R, and thereafter the crack deflects to other cleavage planes. Indeed, C_R defined by the continuum mechanics ignores atomistic phenomena occurring during rapid crack propagation, and therefore is limited in predicting the crack velocity. Examination of these anomalies shows that this limitation lies in microstructural lattice arrangement and in anisotropic phonon radiation during rapid crack propagation. Globally, C_R has no influence on the crack deflection phenomenon. However, the misfit in C_R between the original plane of propagation and the deflected plane generates local instabilities along the deflection zone.     

On the relation between particle velocity and shock wave speed for thermoelastic materials

Results of shock-wave experiments in solids often suggest a nearly-linear relation between the particle velocity behind the shock and the shock wave speed. The present note reconsiders the question of whether thermoelastic material models may be consistent with such observations. Emphasis is placed on the role played by the response of the material in severe compression, as distinguished from its response for small or moderate deformations. The details are illustrated for materials of Mie-Grüneisen type. Received 21 December 2001 / Accepted 22 April 2002 Published online 8 July 2002