A Note on Perturbation Formulae for the Surface-Wave Speed Due To Perturbations in Material Properties

In a recent paper by Tanuma and Man, a two-term asymptotic formula was derived for the speed of surface waves propagating in an anisotropic elastic half-space whose elastic moduli differ only slightly from those for a (base) isotropic elastic material. This formula disagrees with that derived by Delsanto and Clark in an earlier paper using a different method. In this short note, we use a simple procedure to derive another two-term asymptotic formula for the surface-wave speed. Our formula takes the same compact form even if the base material is generally anisotropic. We show that when an error in the work of Delsanto and Clark is corrected, the three different methods do give equivalent results.      

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.     

Perturbation approach to elastic constitutive relations of polycrystals

Herein we obtain a formula for the effective elastic stiffness tensor C^eff of an orthorhombic aggregate of cubic crystallites by the perturbation method. The effective elastic stiffness tensor of the polycrystal gives the relationship between volume average stress and volume average strain. Under Voigt's model, Reuss’ model and Man's theory, the elastic constitutive relation accounts for the effect of the orientation distribution function (ODF) up to terms linear in the texture coefficients. However, the formula derived in this paper delineates the effect of crystallographic texture on elastic response and shows quadratic texture dependence. The formula is very simple. We also consider the influence of grain shape to elastic constitutive relations of polycrystals. Some examples are given to compare computational results of the formula with those given by Voigt's model, Reuss's model, the finite element method, and the self-consistent method. In Section 3, we also present an expression of the perturbation displacement field, in which Green's function for an orthorhombic aggregate of cubic crystallites is included.     

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.     

立方晶粒任意集合下的格林函数和有效弹性刚度张量

为了推导多晶体材料的有效弹性刚度张量,给出立方晶粒任意集合的格林函数封闭但近似的表达式,该格林函数表达式包含三个单晶弹性常数和多晶体材料五个织构系数,它考虑取向分布函数的影响直至织构系数的线性项,它适用于弱织构多晶体材料或具有弱各向异性晶粒的多晶体材料(如金属铝),它与Nishioka格林函数近似式的比较通过三个算例给出;Synge的格林函数积分式则直接通过数值计算完成,它可作为问题的精确解供参考.该文还简单介绍了多晶体材料有效弹性刚度张量的推导过程,并把所得结果和有限元计算结果进行比较.     

GREEN＇S FUNCTION AND EFFECTIVE ELASTIC STIFFNESS TENSOR FOR ARBITRARY AGGREGATES OF CUBIC CRYSTALS

A closed but approximate formula of Green‘s function for an arbitrary aggregate of cubic crystallites is given to derive the effective elastic stiffness tensor of the polycrystal. This formula, which includes three elastic constants of single cubic crystal and five texture coefficients,accounts for the effects of the orientation distribution function (ODF) up to terms linear in the texture coefficients. Thus it is expected that our formula would be applicable to arbitrary aggregates with weak texture or to materials such as aluminum whose single crystal has weak anisotropy.Three examples are presented to compare predictions from our formula with those from Nishioka and Lothe‘s formula and Synge‘s contour integral through numerical integration. As an application of Green‘s function, we briefly describe the procedure of deriving the effective elastic stiffness tensor for an orthorhombic aggregate of cubic crystallites. The comparison of the computational results given by the finite element method and our effective elastic stiffness tensor is made by an example.     

On the plastic wave speeds in rate-independent elastic–plastic materials with anisotropic elasticity

This paper presents a theoretical study of the speeds of plastic waves in rate-independent elastic–plastic materials with anisotropic elasticity. It is shown that for a given propagation direction the plastic wave speeds are equal to or lower than the corresponding elastic speeds, and a simple expression is provided for the bound on the difference between the elastic and the plastic wave speeds. The bound is given as a function of the plastic modulus and the magnitude of a vector defined by the current stress state and the propagation direction. For elastic–plastic materials with cubic symmetry and with tetragonal symmetry, the upper and lower bounds on the plastic wave speeds are obtained without numerically solving an eigenvalue problem. Numerical examples of materials with cubic symmetry (copper) and with tetragonal symmetry (tin) are presented as a validation of the proposed bounds. The lower bound proposed here on the minimum plastic wave speed may also be used as an efficient alternative to the bifurcation analysis at early stages of plastic deformation for the determination of the loss of ellipticity.     

大延伸非均匀介质中地震波全弹性散射理论Ⅰ--弹性波单次散射理论

所描述的工作聚焦于大延伸非均匀介质中非均匀弹性地震波散射问题的研究.应用Born近似及等效源原理,推出了来自连续横向无界非均匀层的弹性散射波的通解.这一工作是解决大延伸非均匀介质的弹性地震波多次散射问题的基础.在上述通解的基础上,建立了适用于大延伸非均匀介质的全弹性散射理论.该理论可包容小尺度非均匀体、大延伸非均匀介质全弹性波单次弱散射理论及标量波单次弱散射理论,因此可视其为一个更为广义和统一的弱散射理论.     

A simple and explicit algebraic expression for the Rayleigh wave velocity

In this note a factorization technique based on the theory of the Riemann problem is used to derive a compact algebraic formula for the velocity of Rayleigh waves. Unlike previous results based on rationalization and Cardan’s solution of a cubic, the present formulation leads to a formula for the velocity which is a continuous function of the Poission’s ratio and yet is simple enough to be of practical interest. The new formula also enables us to express the complex roots associated with the Rayleigh wave equation as simple functions of the Rayleigh wave velocity.     

On the distance to blow-up of acceleration waves in random media

We study the effects of material spatial randomness on the distance to form shocks from acceleration waves, , in random media. We introduce this randomness by taking the material coefficients and – that represent the dissipation and elastic nonlinearity, respectively, in the governing Bernoulli equation – as a stochastic vector process. The focus of our investigation is the resulting stochastic, rather than deterministic as in classical continuum mechanics studies, competition of dissipation and elastic nonlinearity. Quantitative results for are obtained by the method of moments in special simple cases, and otherwise by the method of maximum entropy. We find that the effect of even very weak random perturbation in and may be very significant on . In particular, the full negative cross-correlation between and $ results in the strongest scatter of , and hence, in the largest probability of shock formation in a given distance x. Received November 6, 2001 / Published online September 4, 2002 Dedicated to Professor Ingo Müller on the occasion of his 65th birthday Communicated by Kolumban Hutter, Darmstadt     

BOUNDARY INTEGRAL FORMULA OF ELASTIC PROBLEMS IN CIRCLE PLANE

IntroductionConcerningtheelasticplaneprobleminaunitcircle ,ZhengShenzhouandZhengXueliangdevelopedaboundaryintegralformulaofthestressfunction[1]:Φ(r,θ) =-( 1 -r2 ) 24π ∫2π0ν( φ)1 -2rcos(θ-φ) r2 dφ　　 12π∫2π011 -2rcos(θ-ω) r2 dω∫2π0μ( φ)1 -cos(ω-φ) dφ　　 1 -r22π∫2π0μ( φ)1 -2rcos(θ -φ) r2 dφ　　 ( 0 ≤r <1 ) ,( 1 )whereμ(θ) =Φ(r,θ) ｜r=1,ν(θ) = Φ n r=1= Φ r r=1.Intheformula ( 1 )theseconditemisastrongsingularintegral,itshouldbeunderstoodasanintegra…     

MULTIPLE SCATTERING AND DYNAMIC STRESS ANALYSIS OF ELASTIC WAVES IN A FIBER—REINFORCED COMPOSITE WITH INTERFACES

Based on the theory of elastic dynamics, multiple scattering of elastic waves and dynamic stress concentrations in fiber-reinforced composite are studied. The analytical expressions of elastic waves in different regions are presented. The mode coefficients of elastic waves are determined in accordance with the continuous conditions of displacement and stress on the boundary of the multi-interfaces. By using the addition theorem of Hankel functions, the formula of scattered wave fields in different local coordinates are transformed into those in one local coordinate to determine the unknown coefficients and dynamic stress concentration factors (DSCFs). The influences of the distance between two inclusions, material properties and structural size on the DSCFs near the interfaces are analyzed. As examples, the numerical results of DSCFs near the interfaces for two kinds of fiber-reinforced composites are presented and discussed. The project supported by the National Natural Science Foundation of China (19972018)     

弹塑性激波衰减规律的一种简便解法

本文提出了一个求解弹塑性激波衰减规律的级数解法,给出了各未知量级数展式的通项表达式。同时给出了塑性激波后方弹性卸载区内的应力和质点速度分布规律。该方法简单方便,级数收敛得很快。     

Nonlinear traveling waves in a compressible Mooney-Rivlin rod I. Long finite-amplitude waves

In literature, nonlinear traveling waves in elastic circular rods have only been studied based on single partial differential equation (pde) models, and here we consider such a problem by using a more accurate coupled-pde model. We derive the Hamiltonian from the model equations for the long finite-amplitude wave approximation, analyze how the number of singular points of the system changes with the parameters, and study the features of these singular points qualitatively. Various physically acceptable nonlinear traveling waves are also discussed, and corresponding examples are given. In particular, we find that certain waves, which cannot be counted by the single-equation model, can arise. The project supported by the Research Grants Council of the HKSAR, China (City U 1107/99P) and the National Natural Science Foundation of China (10372054 and 10171061)     

Equivalence Theorems on the Propagation of Small-Amplitude Waves in Prestressed Linearly Elastic Materials with Internal Constraints

By extending the procedure of linearization for constrained elastic materials in the papers by Marlow and Chadwick et al., we set up a linearized theory of constrained materials with initial stress (not necessarily based on a nonlinear theory). The conditions of propagation are characterized for small-displacement waves that may be either of discontinuity type of any given order or, in the homogeneous case, plane progressive. We see that, just as in the unconstrained case, the laws of propagation of discontinuity waves are the same as those of progressive waves. Waves are classified as mixed, kinematic, or ghost. Then we prove that the analogues of Truesdell"s two equivalence theorems on wave propagation in finite elasticity hold for each type of wave. This revised version was published online in August 2006 with corrections to the Cover Date.     

Rayleigh waves in an isotropic elastic half-space coated by a thin isotropic elastic layer with smooth contact

In the present paper, we are interested in the propagation of Rayleigh waves in an isotropic elastic half-space coated with a thin isotropic elastic layer. The contact between the layer and the half space is assumed to be smooth. The main purpose of the paper is to establish an approximate secular equation of the wave. By using the effective boundary condition method, an approximate, yet highly accurate secular equation of fourth-order in terms of the dimensionless thickness of the layer is derived. From the secular equation obtained, an approximate formula of third-order for the velocity of Rayleigh waves is established. The approximate secular equation and the formula for the velocity obtained in this paper are potentially useful in many practical applications.     

Dynamics of a prestressed stiff layer on an elastic half space: filtering and band gap characteristics of periodic structural models derived from long-wave asymptotics

Anti-plane and plane-strain, time-harmonic, small-amplitude vibrations of an elastic layer on an elastic half space are considered, superimposed upon a state of finite, uniform stress and strain. A (compressible) elastic material is considered, orthotropic with orthotropy axes aligned parallel and orthogonal both to the layer and the prestress principal directions. A non-uniform mass density is assumed in the layer. A formal long-wave asymptotic solution is derived under the assumptions of high contrast between the stiffnesses of the layer and the half space and between certain prestress components and the current elastic shear modulus.It is shown that (i) the layer asymptotically behaves as a beam subject to transversal and axial vibrations; (ii) the response of the half space can be found in a closed-form, under the assumption of plane wave motion (which becomes consistent when the density of the layer is uniform), otherwise it is represented by a hypersingular integral equation; (iii) if the nonlocality introduced by the hypersingular integral equation is restricted to an influence area of finite extent, the integral can be analytically approximated, so that a Winkler-type spring model representing the half space is rigorously derived. For uniform density of the layer, the constants defining the spring model are given as functions of the prestress and anisotropy parameters of the half space; and, finally, (iv) the asymptotic solution provides new analytical expressions for incremental displacement of the layer, which, compared to the exact numerical solution (also included), are shown to perform quite well, even for values of parameters much beyond the limits imposed by the asymptotic analysis.The asymptotic analysis allows us to explore, for the first time, dynamic properties of a periodic layer bonded to an elastic half space and subject to a uniform prestress state. We find that the system exhibits band gaps (ranges of forbidden frequencies) and that the prestress can be used as a parameter tuning the filtering properties of the structure, an effect which may have important consequences in the design of resonant devices.     

Evolution Equations for Plane Cubically Nonlinear Elastic Waves

Consideration is given to the nonlinear theory of elastic waves with cubic nonlinearity. This nonlinearity is separated out, and the interaction of four harmonic waves is studied. The method of slowly varying amplitudes is used. The shortened and evolution equations, the first integrals of these equations (Manley–Rowe relations), and energy balance law for a set of four interacting waves (quadruplet) are derived. The interaction of waves is described using the wavefront reversal scheme     

From imaging to material identification: A generalized concept of topological sensitivity

To establish a compact analytical framework for the preliminary stress-wave identification of material defects, the focus of this study is an extension of the concept of topological derivative, rooted in elastostatics and the idea of cavity nucleation, to 3D elastodynamics involving germination of solid obstacles. The main result of the proposed generalization is an expression for topological sensitivity, explicit in terms of the elastodynamic Green's function, obtained by an asymptotic expansion of a misfit-type cost functional with respect to the nucleation of a dissimilar elastic inclusion in a defect-free “reference” solid. The featured formula, consisting of an inertial-contrast monopole term and an elasticity-contrast dipole term, is shown to be applicable to a variety of reference solids (semi-infinite and infinite domains with constant or functionally graded elastic properties) for which the Green's functions are available. To deal with situations when the latter is not the case (e.g. finite reference bodies or those with pre-existing defects), an adjoint field approach is employed to derive an alternative expression for topological sensitivity that involves the contraction of two (numerically computed) elastodynamic states. A set of numerical results is included to demonstrate the potential of generalized topological derivative as an efficient tool for exposing not only the geometry, but also material characteristics of subsurface material defects through a local, point-wise identification of “optimal” inclusion properties that minimize the topological sensitivity at sampling location. Beyond the realm of non-invasive characterization of engineered materials, the proposed developments may be relevant to medical diagnosis and in particular to breast cancer detection where focused ultrasound waves show a promise of superseding manual palpation.