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.     

各向异性粘弹性多孔截止中平面波的传播

This study discusses wave propagation in perhaps the most general model of a poroelastic medium. The medium is considered as a viscoelastic, anisotropic and porous solid frame such that its pores of anisotropic permeability are filled with a viscous fluid. The anisotropy considered is of general type, and the attenuating waves in the medium are treated as the inhomogeneous waves. The complex slowness vector is resolved to define the phase velocity, homogeneous attenuation, inhomogeneous attenuation, and angle of attenuation for each of the four attenuating waves in the medium. A non-dimensional parameter measures the deviation of an inhomogeneous wave from its homogeneous version. An numerical model of a North-Sea sandstone is used to analyze the effects of the propagation direction, inhomogeneity parameter, frequency regime, anisotropy symmetry, anelasticity of the frame, and viscosity of the pore-fluid on the propagation characteristics of waves in such a medium.     

SH-type waves dispersion in an isotropic medium sandwiched between an initially stressed orthotropic and heterogeneous semi-infinite media

The present paper studies the propagation of shear waves (SH-type waves) in an homogeneous isotropic medium sandwiched between two semi infinite media. The upper half-space is considered as orthotropic medium under initial stress and lower half-space considered as heterogeneous medium. We have obtained the dispersion equation of phase velocity for SH-type waves. The propagation of SH-type waves are influenced by inhomogeneity parameters and initial stress parameter. The velocity of SH-type wave has been computed for different cases. We have also obtained the dispersion equation of phase velocity in homogeneous media in the absence of initial stress. The velocities of SH-type waves are calculated numerically as a function of kH (non-dimensional wave number) and presented in a number of graphs. To study the effect of inhomogeneity parameters and initial stress parameter we have plotted the velocity of SH-type wave in several figure. We have observed that the velocity of wave increases with the increase inhomogeneity parameters. We found that in both homogeneous and inhomogeneous media the velocity of SH-type wave increases with the increase of initial stress parameter. The results may be useful for the study of seismic waves propagation during any earthquake and artificial explosions.     

On the theory of plane inhomogeneous waves in anisotropic elastic media

In the theory of plane inhomogeneous elastic waves, the complex wave vector constituted by two real vectors in a given plane may be described with the aid of two complex scalar parameters. Either of those parameters may be taken as a free one in the characteristic condition assigned to the wave equation. This alternative underlies the two fundamental approaches in the theory, namely, one associated with the Stroh eigenvalue problem and the other with the generalized Christoffel eigenvalue problem. The two approaches are identical insofar as a partial nondegenerate wave solution (partial mode) is concerned, but they differ in the fundamental solution (wave packet) assembling, and their dissimilarity is also revealed in the presence of degeneracies, which may involve either of the two governing parameters or both of them. Therefore, use of both approaches is essential for studying the degeneracy phenomenon in the theory of inhomogeneous waves. The criteria for different types of degeneracy, related to a double eigenvalue of the Stroh matrix or the Christoffel matrix and at the same time to a repeated root of the characteristic condition, are formulated by appeal to the matrix algebra and to the theory of polynomial equations. On this basis, dimensions of the manifolds, associated with degeneracy of different types in the space of variables, are established for elastic media of unrestricted anisotropy. The relation to the boundary-value problems is discussed.     

Superposition of finite-amplitude waves propagating in a principal plane of a deformed Mooney–Rivlin material

Here we consider finite-amplitude wave motions in Mooney–Rivlin elastic materials which are first subjected to a static homogeneous deformation (prestrain). We assume that the time-dependent displacement superimposed on the prestrain is along a principal axis of the prestrain and depends on two spatial variables in the principal plane orthogonal to this axis. Thus all waves considered here are linearly polarized along this axis. After retrieving known results for a single homogeneous plane wave propagating in a principal plane, a superposition of an arbitrary number of sinusoidal homogeneous plane waves is shown to be a solution of the equations of motion. Also, inhomogeneous plane wave solutions with complex wave vector in a principal plane and complex frequency are obtained. Moreover, appropriate superpositions of such inhomogeneous waves are also shown to be solutions. In each case, expressions are obtained for the energy density and energy flux associated with the wave motion.     

位势流动和非均匀介质中声波方程的36种形式

声波方程是对大多数声学问题进行数学描述的出发点. 那些得到 广泛应用的经典波动方程及对流波动方程都存在苛刻的适用条件, 即仅适用于描述处于静态或匀速运动状态的定常 均匀介质中的线性无耗散声波. 然而, 很多实际场合并不满足这些严格的适用条件. 本文对经典声波方程和对流声波 方程进行推广, 导出了编号为W1$\sim$W36的36种不同形式的声波方程, 涵盖了处于静止、势流或旋涡流状态下的非均匀 和/或非定常介质中的声波传播问题. 所考虑的声波传播情形包括: (1) 线性波, 即具有小梯度(小振幅)性质; (2)非线性波, 即具有陡峭梯度性质, 包括``波纹'(小振幅大梯度)或者大振幅波. 本文仅考虑非耗散声波, 即排除了由剪切、体积黏度及热传导所引起的耗散. 对具有匀熵或等熵(熵沿流线守恒)性质的均匀介质和非均匀介质中的声传播进行了研究但非等熵(即耗散)情况除外; 另外, 对非定常介质中的 声波问题也进行了分析. 所涉及的介质可以处于静止、匀速运动状态, 或者是非匀速的和/或非定常的平均流动, 包括: (1)低Mach数的势平均流(即不可压缩的平均态), 或高速势平均流(即非均匀可压缩的平均流); ② 变截面管 道中的准一维传播, 包括无平均流的号管和具有低或高Mach数平均流的喷管; 或③平面的、空间的、或轴对称的单 向剪切平均流. 本文没有探讨其他类型的旋涡平均流(将与耗散及其他情形一起留待下一步研究), 例如, 可能与剪切效应相结合的轴对称旋转平均流. 通过对流体力学的一般方程进行消元处理或根据声学变分原理, 导出了36种波动方程, 对一些波动方程还采用这两种方法进行相互校验. 尽管声波方程的36种形式没有涵盖非线性、非均匀与非定常及非匀速运动介质 这3个效应的所有可能的组合情形, 但它们的确包括了孤立状态下的各种效应, 并包括了多种多重效应组合的 情形. 虽然经典波动方程和对流波动方程仅适用于处于静止(或匀速运动)的均匀定常介质中的线性无耗散声波, 但它们在 相关文献中已被广泛采用; 本文给出的36种声波方程提供了它们多种有用的推广形式. 在许多实际应用中, 经典波动方 程和对流波动方程仅是粗略的近似, 声波方程的更一般形式可提供更令人满意的理论模型. 本文每节末尾给出了这些应用 的众多范例. 在这篇评论文章中引用了240篇参考文献.     

径向非均匀介质中圆形夹杂的动应力分析

基于复变函数理论,研究了径向非均匀弹性介质中均匀圆夹杂对弹性波的散射问题. 介质的非均匀性体现在介质密度沿着径向按幂函数形式变化且剪切模量是常数. 利用坐标变换法将变系数的非均匀波动方程转为标准亥姆霍兹(Helmholtz) 方程. 在复坐标系下求得非均匀基体和均匀夹杂同时存在的位移和应力表达式. 通过具体算例分析了圆夹杂周边的动应力集中系数(DSCF). 结果表明:基体与夹杂的波数比和剪切模量比,基体的参考波数和非均匀参数对动应力集中有较大的影响.      

Free-field wave solutions in a half-plane exhibiting a special-type of continuous inhomogeneity

This work presents closed-form solutions for free-field motions in a continuously inhomogeneous half-plane that include contributions of incident waves as well as of waves reflected from the traction-free horizontal surface. Both pressure and vertically polarized shear waves are considered. Furthermore, two special types of material inhomogeneity are studied, namely (a) a shear modulus that varies quadratically with respect to the depth coordinate and (b) one that varies exponentially with the said coordinate. In all cases, Poisson’s ratio is fixed at one-quarter, while both shear modulus and material density profiles vary proportionally. Next, a series of numerical results serve to validate the aforementioned models, and to show the differences in the wave motion patterns developing in media that are inhomogeneous as compared to a reference homogeneous background. These results clearly show the influence of inhomogeneity, as summarized by a single material parameter, on the free-field motions that develop in the half-plane. It is believed that this type of information is useful within the context of wave propagation studies in non-homogeneous continua, which in turn find applications in fields as diverse as laminated composites, geophysical prospecting, oil exploration and earthquake engineering.     

非均匀吸收介质中地震波的广义传播矩阵解法

本文在复频域内,通过应用混合变量粘弹性波方程和线性常微分方程组的指数矩阵解法,给出了一种计算非均匀吸收介质中地震波传播的广义传播矩阵解法。该方法适用于各种粘弹性模型,可模拟任意震源及所产生的各种体波、面波,数值结果表明具有很高的计算精度。     

Low-frequency dilatational wave propagation through unsaturated porous media containing two immiscible fluids

An analytical theory is presented for the low-frequency behavior of dilatational waves propagating through a homogeneous elastic porous medium containing two immiscible fluids. The theory is based on the Berryman–Thigpen–Chin (BTC) model, in which capillary pressure effects are neglected. We show that the BTC model equations in the frequency domain can be transformed, at sufficiently low frequencies, into a dissipative wave equation (telegraph equation) and a propagating wave equation in the time domain. These partial differential equations describe two independent modes of dilatational wave motion that are analogous to the Biot fast and slow compressional waves in a single-fluid system. The equations can be solved analytically under a variety of initial and boundary conditions. The stipulation of “low frequency” underlying the derivation of our equations in the time domain is shown to require that the excitation frequency of wave motions be much smaller than a critical frequency. This frequency is shown to be the inverse of an intrinsic time scale that depends on an effective kinematic shear viscosity of the interstitial fluids and the intrinsic permeability of the porous medium. Numerical calculations indicate that the critical frequency in both unconsolidated and consolidated materials containing water and a nonaqueous phase liquid ranges typically from kHz to MHz. Thus engineering problems involving the dynamic response of an unsaturated porous medium to low excitation frequencies (e.g., seismic wave stimulation) should be accurately modeled by our equations after suitable initial and boundary conditions are imposed.     

Wave propagation in transversely isotropic porous piezoelectric materials

Wave propagation in porous piezoelectric material (PPM), having crystal symmetry 6 mm, is studied analytically. Christoffel equation is derived for the propagation of plane harmonic waves in such a medium. The roots of this equation give four complex wave velocities which can propagate in such materials. The phase velocities of propagation and the attenuation quality factors of all these waves are described in terms of complex wave velocities. Phase velocities and attenuation of the waves in PPM depend on the phase direction. Numerical results are computed for the PPM BaTiO₃. The variation of phase velocity and attenuation quality factor with phase direction, porosity and the wave frequency is studied. The effects of anisotropy and piezoelectric coupling are also studied. The phase velocities of two quasi dilatational waves and one quasi shear waves get affected due to piezoelectric coupling while that of type 2 quasi shear wave remain unaffected. The phase velocities of all the four waves show non-dispersive behavior after certain critical high frequency. The phase velocity of all waves decreases with porosity while attenuation of respective waves increases with porosity of the medium. The characteristic curves, including slowness curves, velocity curves, and the attenuation curves, are also studied in this paper.     

上覆单相弹性层的饱和地基上弹性圆板轴对称竖向振动分析

基于弹性层和饱和土半空间轴对称弹性波动方程，运用Hankel积分变换方法，得到它们在变换空间内的解．进而由层间完全接触条件及圆板底面的混合边值条件，构造一组描述上覆单相弹性层的饱和地基上弹性圆板轴对称竖向振动的对偶积分方程．将该对偶积分方程化为易于数值计算的第二类Fredholm积分方程，求得地基表面动力柔度系数随无量纲频率的变化曲线及弹性圆板的相对位移幅值．     

On theory of generalized thermoelastic solids with voids and diffusion

Governing equations of thermoelastic diffusion material with voids are modified with the help of Lord and Shulman theory of generalized thermoelasticity. These governing equations are then solved in two-dimension to show the existence of four coupled longitudinal waves and a shear wave. The complex absolute values of the speeds of the coupled longitudinal waves are computed numerically against the frequency for Magnesium material. The reflection of these plane waves from a stress free thermally insulated boundary is also studied, where the dependence of the reflection coefficients on angle of incidence is shown graphically for the incidence of coupled longitudinal wave only. The speeds and reflection coefficients of plane waves are also computed numerically in the absence of voids and diffusion parameters, which are shown graphically to observe the effects of voids and diffusion.     

Squirt-Flow in Fluid-Saturated Porous Media: Propagation of Rayleigh Waves

Propagation of attenuated waves is studied in a squirt-flow model of porous solid permeated by two different pore regimes saturated with same viscous fluid. Presence of soft compliant microcracks embedded in the grains of stiff porous rock defines the double-porosity formation. Microcracks and pores respond differently to the compressional effect of a propagating wave, which induces the squirt-flow from microcracks to pores. Elastodynamics of constituent particles in porous aggregate is represented through a single-porosity formulation, which involves the frequency-dependent complex moduli. This formulation is deduced as a special case of double-porosity formation allowing the wave-induced flow of pore-fluid. This squirt-flow model of porous solid supports the attenuated propagation of two compressional waves and one shear wave. Superposition of these body waves, subject to stress-free surface, defines the propagation of Rayleigh wave. This wave is governed by a complex irrational dispersion equation, which is solved numerically after rationalising into an algebraic equation. For existence of Rayleigh wave, a complex solution of the dispersion equation should represent a leaky wave, which decays for propagation along any direction in the semi-infinite medium. A numerical example is solved to analyse the effects of squirt-flow on phase velocity, attenuation and polarisation of the Rayleigh waves, for different combinations of parameters. Numerical results suggest the existence of an additional (second) Rayleigh wave in the squirt-flow model of dissipative porous solids.     

Lamb wave scattering by a symmetric pair of surface-breaking cracks in a plate

The scattering problem of a Lamb wave incident on a symmetric pair of surface-breaking transverse cracks in a plate is considered. The Lamb wave is assumed to be obliquely incident on the crack plane. Since the cracks are part-through, the scattered field will contain reflected as well as transmitted waves. The energy of the incoming wave is partitioned into reflected and transmitted wave modes. Energy coefficients of the reflected and transmitted waves are calculated as a function of incident frequency and crack depth. The incidence angle of the incoming wave is also treated as a parameter. Both the reflected and transmitted wave fields are considered as linear superpositions of all real and complex wave modes in the plate. Decomposition of modes is achieved with the help of an orthogonality condition based on the principle of reciprocal work. Continuity of displacement and stress fields is imposed at the crack plane. Energy coefficients for reflection and transmission are obtained from the mode amplitudes. Energy coefficients are shown to be a strong function of incident frequency and crack depth. Experiments are conducted with a PZT transducer network interacting with a symmetric pair of machined cracks in an aluminum plate. Trends predicted by the analysis are reflected in the experimental results.     

Propagation velocities of elastic waves in saturated soils

Based on the wave equations established by the authors, the characteristics of propagation velocities of elastic waves in saturated soils are analyzed and verified by ultrasonic test in laboratory and seismic survey in the field. The results provide theoretical basis for the determination of physical and mechanical parameters of saturated soils using propagation velocities of elastic waves, especially P-wave Velocity.     

On a modification of the wave equation for a layered medium

A new form of the wave equation in inhomogeneous media is presented which does not contain derivatives of the medium parameters in its coefficients. Hence this equation can be used not only for the case of smooth but also for the case of abrupt changes of the parameters with the coordinates. The equation can be used for waves of different nature.

To illustrate the advantages of the new form of the wave equation four problems have been solved. They are: scattering of a plane sound wave by weak inhomogeneities; excitation of a lateral wave; the symmetry of the plane-wave transmission coefficient with respect to inversion of the path of the wave; and plane were reflection from a thin inhomogeneous layer.     

Anti-Plane Shear Waves Scattering from a Partially Debonded Magneto-Electro-Elastic Circular Cylindrical Inhomogeneity

This article presents an analysis of the scattering of anti-plane shear waves from a single cylindrical inhomogeneity partially bonded to an unbounded magneto-electro-elastic matrix. The magneto-electric permeable boundary conditions are adopted. The crack opening displacement is represented by Chebyshev polynomials and a system of equations is derived and solved for the unknown coefficients. Some examples are calculated and the results are illustrated. The results show that the COD increases when the piezomagnetic coefficient of the inhomogeneity bonded to the piezoelectric matrix becomes larger, and that the COD decreases when the piezomagnetic coefficient of the matrix with the piezoelectric inhomogeneity increases.