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.     

Stroh Formalism and Rayleigh Waves

The Stroh formalism is a powerful and elegant mathematical method developed for the analysis of the equations of anisotropic elasticity. The purpose of this exposition is to introduce the essence of this formalism and demonstrate its effectiveness in both static and dynamic elasticity. The equations of elasticity are complicated, because they constitute a system and, particularly for the anisotropic cases, inherit many parameters from the elasticity tensor. The Stroh formalism reveals simple structures hidden in the equations of anisotropic elasticity and provides a systematic approach to these equations. This exposition is divided into three chapters. Chapter 1 gives a succinct introduction to the Stroh formalism so that the reader could grasp the essentials as quickly as possible. In Chapter 2 several important topics in static elasticity, which include fundamental solutions, piezoelectricity, and inverse boundary value problems, are studied on the basis of the Stroh formalism. Chapter 3 is devoted to Rayleigh waves, for long a topic of utmost importance in nondestructive evaluation, seismology, and materials science. There we discuss existence, uniqueness, phase velocity, polarization, and perturbation of Rayleigh waves through the Stroh formalism.

---

The Table of Contents and Index are also provided as Electronic Supplementary Material for online readers at doi:      

Angular Dependence of Rayleigh-Wave Velocity in Prestressed Polycrystalline Media with Monoclinic Texture

This article is concerned with Rayleigh waves propagating along the free surface of a macroscopically homogeneous, prestressed half-space. In the meso-scale, the half-space in question is taken to be a textured polycrystalline aggregate of cubic crystallites, which has the normal to its free surface being a 2-fold axis of monoclinic sample symmetry. Under the theoretical framework of linear elasticity with initial stress, an angular dependence formula, which shows explicitly how the phase velocity of Rayleigh waves depends on the propagation direction, the prestress, and the crystallographic texture, is derived from a constitutive equation motivated by Hartig's law. This velocity formula includes terms which describe the effects of texture on acoustoelastic coefficients, and it is correct to within terms linear in the initial stress and in the anisotropic part of the incremental elasticity tensor. Since its derivation makes no presumption on the origin of the initial stress, this velocity formula is meant to be applicable when the prestress is induced by plastic deformations such as those incurred during the surface enhancement treatment of low plasticity burnishing. The angular dependence formula assumes a simpler form when the texture of the prestressed half-space is orthorhombic. This revised version was published online in July 2006 with corrections to the Cover Date.     

Wave propagation in thermally conducting linear fibre-reinforced composite materials

The propagation of plane waves in a fibre-reinforced, anisotropic, generalized thermoelastic media is discussed. The governing equations in x–y plane are solved to obtain a cubic equation in phase velocity. Three coupled waves, namely quasi-P, quasi-SV and quasi-thermal waves are shown to exist. The propagation of Rayleigh waves in stress free thermally insulated and transversely isotropic fibre-reinforced thermoelastic solid half-space is also investigated. The frequency equation is obtained for these waves. The velocities of the plane waves are shown graphically with the angle of propagation. The numerical results are also compared to those without thermal disturbances and anisotropy parameters.     

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

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

Rayleigh waves in an incompressible elastic half-space overlaid with a water layer under the effect of gravity

This paper is concerned with the propagation of Rayleigh waves in an incompressible isotropic elastic half-space overlaid with a layer of non-viscous incompressible water under the effect of gravity. The authors have derived the exact secular equation of the wave which did not appear in the literature. Based on it the existence of Rayleigh waves is considered. It is shown that a Rayleigh wave can be possible or not, and when a Rayleigh wave exists it is not necessary unique. From the exact secular equation the authors arrive immediately at the first-order approximate secular equation derived by Bromwich [Proc. Lond. Math. Soc. 30:98–120, 1898]. When the layer is assumed to be thin, a fourth-order approximate secular equation is derived and of which the first-order approximate secular equation obtained by Bromwich is a special case. Some approximate formulas for the velocity of Rayleigh waves are established. In particular, when the layer being thin and the effect of gravity being small, a second-order approximate formula for the velocity is created which recovers the first-order approximate formula obtained by Bromwich [Proc. Lond. Math. Soc. 30:98–120, 1898]. For the case of thin layer, a second-order approximate formula for the velocity is provided and an approximation, called global approximation, for it is derived by using the best approximate second-order polynomials of the third- and fourth-powers.     

Elastic Anisotropic Material with Purely Longitudinal and Transverse Waves

The simplest form of the matrix of elasticity moduli of an anisotropic material conducting purely longitudinal and transverse waves with an arbitrary direction of the wave normal is obtained. A generic solution of equations in displacements is represented in terms of three functions satisfying independent wave equations. In the case of planar deformation, this solution yields a complex representation coinciding with the Kolosov–Muskhelishvili formulas for an isotropic material. The formulas in the present work also determine an anisotropic material with Young's modulus identical for all directions, as in an isotropic medium.     

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.     

Initial stresses in elastic solids: Constitutive laws and acoustoelasticity

On the basis of the nonlinear theory of elasticity, the general constitutive equation for an isotropic hyperelastic solid in the presence of initial stress is derived. This derivation involves invariants that couple the deformation with the initial stress and in general, for a compressible material, it requires 10 invariants, reducing to 9 for an incompressible material. Expressions for the Cauchy and nominal stress tensors in a finitely deformed configuration are given along with the elasticity tensor and its specialization to the initially stressed undeformed configuration. The equations governing infinitesimal motions superimposed on a finite deformation are then used to study the combined effects of initial stress and finite deformation on the propagation of homogeneous plane waves in a homogeneously deformed and initially stressed solid of infinite extent. This general framework allows for various different specializations, which make contact with earlier works. In particular, connections with results derived within Biot's classical theory are highlighted. The general results are also specialized to the case of a small initial stress and a small pre-deformation, i.e. to the evaluation of the acoustoelastic effect. Here the formulas derived for the wave speeds cover the case of a second-order elastic solid without initial stress and subject to a uniaxial tension [Hughes and Kelly, Phys. Rev. 92 (1953) 1145] and are consistent with results for an undeformed solid subject to a residual stress [Man and Lu, J. Elasticity 17 (1987) 159]. These formulas provide a basis for acoustic evaluation of the second- and third-order elasticity constants and of the residual stresses. The results are further illustrated in respect of a prototype model of nonlinear elasticity with initial stress, allowing for both finite deformation and nonlinear dependence on the initial stress.     

Edge waves under initial stress

The frequency equation giving the phase velocity of a wave at the edge of a thick plate under initial stress is obtained. Some particular cases are discussed to derive (a) the velocity of edge waves in a thin plate and (b) the velocity of Rayleigh waves in a plate of infinite thickness under initial stress. The results are compared with those for zero initial stress.     

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.     

Stroh formulation for a generally constrained and pre-stressed elastic material

The Stroh formalism is essentially a spatial Hamiltonian formulation and has been recognized to be a powerful tool for solving elasticity problems involving generally anisotropic elastic materials for which conventional methods developed for isotropic materials become intractable. In this paper we develop the Stroh/Hamiltonian formulation for a generally constrained and prestressed elastic material. We derive the corresponding integral representation for the surface-impedance tensor and explain how it can be used, together with a matrix Riccati equation, to calculate the surface-wave speed. The proposed algorithm can deal with any form of constraint, pre-stress, and direction of wave propagation. As an illustration, previously known results are reproduced for surface waves in a pre-stressed incompressible elastic material and an unstressed inextensible fibre-reinforced composite, and an additional example is included analyzing the effects of pre-stress upon surface waves in an inextensible material.     

A Natural Generalization of Hypoelasticity and Eulerian Rate Type Formulation of Hyperelasticity

According to the classical hypoelasticity theory, the hypoelasticity tensor, i.e. the fourth order Eulerian constitutive tensor, characterizing the linear relationship between the stretching and an objective stress rate, is dependent on the current stress and must be isotropic. Although the classical hypoelasticity in this sense includes as a particular case the isotropic elasticity, it fails to incorporate any given type of anisotropic elasticity. This implies that one can formulate the isotropic elasticity as an integrable-exactly classical hypoelastic relation, whereas one can in no way do the same for any given type of anisotropic elasticity. A generalization of classical theory is available, which assumes that the material time derivative of the rotated stress is dependent on the rotated Cauchy stress, the rotated stretching and a Lagrangean spin, linear and of the first degree in the latter two. As compared with the original idea of classical hypoelasticity, perhaps the just-mentioned generalization might be somewhat drastic. In this article, we show that, merely replacing the isotropy property of the aforementioned stress-dependent hypoelasticity tensor with the invariance property of the latter under an R-rotating material symmetry group R⋆ G ₀, one may establish a natural generalization of classical theory, which includes all of elasticity. Here R is the rotation tensor in the polar decomposition of the deformation gradient and G ₀ any given initial material symmetry group. In particular, the classical case is recovered whenever the material symmetry is assumed to be isotropic. With the new generalization it is demonstrated that any two non-integrable hypoelastic relations based on any two objective stress rates predict quite different path-dependent responses in nature and hence can in no sense be equivalent. Thus, the non-integrable hypoelastic relations based on any given objective stress rate constitute an independent constitutive class in its own right which is disjoint with and hence distinguishes itself from any class based on another objective stress rate. Only for elasticity, equivalent hypoelastic formulations based on different stress rates may be established. Moreover, universal integrability conditions are derived for all kinds of objective corotational stress rates and for all types of material symmetry. Explicit, simple, integrable-exactly hypoelastic relations based on the newly discovered logarithmic stress rate are presented to characterize hyperelasticity with any given type of material symmetry. It is shown that, to achieve the latter goal, the logarithmic stress rate is the only choice among all infinitely many objective corotational stress rates. This revised version was published online in August 2006 with corrections to the Cover Date.     

Rayleigh and Love surface waves in isotropic media with negative Poisson’s ratio

The behavior of Rayleigh surface waves and the first mode of the Love waves in isotropic media with positive and negative Poisson’s ratio is compared. It is shown that the Rayleigh wave velocity increases with decreasing Poisson’s ratio, and it increases especially rapidly for negative Poisson’s ratios less than ?0.75. It is demonstrated that, for positive Poisson’s ratios, the vertical component of the Rayleigh wave displacements decays with depth after some initial increase, while for negative Poisson’s ratios, there is a monotone decrease. The Rayleigh waves are characterized by elliptic trajectories of the particle motion with the change of the rotation direction at critical depths and by the linear vertical polarization at these depths. It is found that the elliptic orbits are less elongated and the critical depths are greater for negative Poisson’s ratios. It is shown that the stress distribution in the Rayleighwaves varies nonmonotonically with the dimensionless depth as (positive or negative) Poisson’s ratio varies. The stresses increase strongly only as Poisson’s ratio tends to?1. It is shown that, in the case of an incompressible thin covering layer, the velocity of the first mode of the Love waves strongly increases for negative Poisson’s ratios of the half-space material. If the thickness of the incompressible layer is large, then the wave very weakly penetrates into the halfspace for any value of its Poisson’s ratio. For negative Poisson’s ratios, the Love wave in a layer and a half-space is mainly localized in the covering layer for any values of its thickness and weakly penetrates into the half-space. For the first mode of the Love waves, it was discovered that there is a strong increase in the maximum of one of the shear stresses on the interface between the covering layer and the half-space as Poisson’s ratios of both materials decrease. For the other shear stress, there is a stress jump on the interface and a more complicated dependence of the stress on Poisson’s ratio on both sides of the interface.     

Effects of curvilinear anisotropy on radially symmetric stresses in anisotropic linearly elastic solids

It has been known for some time that certain radial anisotropies in some linear elasticity problems can give rise to stress singularities which are absent in the corresponding isotropic problems. Recently related issues were examined by other authors in the context of plane strain axisymmetric deformations of a hollow circular cylindrically anisotropic linearly elastic cylinder under uniform external pressure, an anisotropic analog of the classic isotropic Lamé problem. In the isotropic case, as the external radius increases, the stresses rapidly approach those for a traction-free cavity in an infinite medium under remotely applied uniform compression. However, it has been shown that this does not occur when the cylinder is even slightly anisotropic. In this paper, we provide further elaboration on these issues. For the externally pressurized hollow cylinder (or disk), it is shown that for radially orthotropic materials, the maximum hoop stress occurs always on the inner boundary (as in the isotropic case) but that the stress concentration factor is infinite. For circumferentially orthotropic materials, if the tube is sufficiently thin, the maximum hoop stress always occurs on the inner boundary whereas for sufficiently thick tubes, the maximum hoop stress occurs at the outer boundary. For the case of an internally pressurized tube, the anisotropic problem does not give rise to such radical differences in stress behavior from the isotropic problem. Such differences do, however, arise in the problem of an anisotropic disk, in plane stress, rotating at a constant angular velocity about its center, as well as in the three-dimensional problem governing radially symmetric deformations of anisotropic externally pressurized hollow spheres. The anisotropies of concern here do arise in technological applications such as the processing of fiber composites as well as the casting of metals.     

Strong-contrast expansion correlation approximations for the effective elastic moduli of multiphase composites

Following the previous approach of Pham and Torquato (J Appl Phys 94:6591–6602, 2003) and Torquato (J Mech Phys Solids 45:1421–1448, 1997; Random heterogeneous media, Springer, Berlin, 2002), we derive the strong-contrast expansions for the effective elastic moduli K _e,G _e of d-dimensional multiphase composites. The series consists of a principal reference part and a fluctuation part (perturbation about a homogeneous reference or comparison material), which contains multi-point correlation functions that characterize the microstructure of the composite. We propose a three-point correlation approximation for the fluctuation part with an objective choice of the reference phase moduli, such that the fluctuation terms vanish. That results in the approximations for the effective elastic moduli of isotropic composites, which coincide with the well-known self-consistent and Maxwell approximations for two-phase composites having respective microstructures. Applications to some two-phase materials are given.     

川藏公路地质环境与整治改建方案的思考

川藏公路由于地质环境复杂、建设标准低、后遗病害多,抗灾能力差,泥石流、滑坡、山崩、雪害、水毁等自然灾害频繁发生,公路阻车断道严重。国家投入巨资进行整治改建,并取得了明显的效果,但由于自然环境特殊、影响因素复杂,许多特大型、大型工程地质病害问题还没有可行、可靠的解决方案。本文通过分析川藏公路沿线的地质环境和灾害特点,总结历年整治改建和经验的教训,提出川藏公路建设的途径、可能达到的目标和应采用的原则。     

无黏流体与固体界面的漏瑞利波

漏瑞利波存在于半无限无黏性流体和半无限固体媒质的界面处. 首先推导流固无限各向同性介质界面处漏瑞利波的特征方程和位移及应力的解析计算公式. 然后结合典型结构通过数值计算研究了漏瑞利波特性以及位移和应力在流体和固体中的分布规律. 数值计算结果表明漏瑞利波的相速度和衰减随流固密度比的增大而增大, 在流固界面上法向位移连续而切向位移不连续. 流固密度比对固体媒质中沿垂直于漏瑞利波的传播方向的位移、正应力和剪应力有比较大的影响,而对沿漏瑞利波的传播方向的正应力几乎没影响. 为利用漏瑞利波的无损检测与评价提供了理论基础.      

Rayleigh waves in anisotropic porous media and the polarization vector method

In this paper, the propagation of Rayleigh waves in orthotropic non-viscous fluid-saturated porous half-spaces with sealed surface-pores and with impervious surface is investigated. The main aim of the investigation is to derive explicit secular equations and based on them to examine the effect of the material parameters and the boundary conditions on the propagation of Rayleigh waves. By employing the method of polarization vector the explicit secular equations have been derived. These equations recover the ones corresponding to Rayleigh waves propagating in purely elastic half-spaces. It is shown from numerical examples that the Rayleigh wave velocity depends strongly on the porosity, the elastic constants, the anisotropy, the boundary conditions and it differs considerably from the one corresponding to purely elastic half-spaces. Remarkably, in the fluid saturated porous half-spaces, Rayleigh waves may travel with a larger velocity than that of the shear wave, a fact that is impossible for the purely elastic half-spaces.