Scattering of torsional waves by a spherical cavity in a long circular elastic cylinder

Summary. The axisymmetric wave equation is solved for the problem of torsional elastic waves scattered by a spherical cavity located symmetrically in an infinitely long circular cylinder. Using Fourier transforms, the problem is reduced to the solution of an infinite system of simultaneous equations, which is suitable for the numerical solution. The numerical results on the transmission and reflection coefficients are shown for various values of sphere radius and frequency. Equation of energy conservation is utilized to check the numerical procedure.     

Propagation of waves from a spherical cavity with and without a shell embedment

Summary A spherical cavity in an infinite, elastic medium with and without a shell embedment is subjected to axisymmetric, non-torsional surface loads in the radial and meridional directions. The so-called Residual Variable Method (RVM) is used to obtain exact, closed-form solutions of the wave propagation problems. Some representative numerical results are presented graphically for the stresses created in two realistic loading situations.     

Boundary integral equation analysis for radiation and scattering of elastic waves in three dimensions

Abstract

In this paper, the boundary integral equation (BIE) method is employed to investigate the radiation and scattering of time‐harmonic elastic waves by obstacles of arbitrary shape embedded in an infinite medium. Based on the vector BIE, entirely free of Cauchy principal value integrals, an efficient numerical scheme using quadratic isoparametric boundary elements is proposed. Furthermore, the difficulty of non‐uniquess of a solution inherent with BIE formulations for exterior elastodynamic problems is studied numerically and analytically. The counterparts of the combined Helmholtz integral formulation method for elastodynamics together with the least‐square or Lagrange‐multiplier technique are derived and applied to overcome this difficulty successfully. In addition, the elastic‐wave fields radiated or scattered by either a spherical cavity or a rigid sphere in an infinite medium are calculated and the results are compared with the analytical solutions to demonstrate the accuracy and versatility of the proposed numerical scheme.     

基于快速多极子基本解方法(FMM-MFS)的弹性波二维散射模拟研究

针对弹性波二维散射问题,发展一种新的快速多极子基本解方法（FMM-MFS）。方法基于单层位势理论,通过在虚边界上设置膨胀波线源和剪切波线源以构造散射波场,从而避免了奇异性的处理和边界单元离散;结合快速多极子展开技术（FMM）,大幅度降低了计算量和存储量,突破了传统方法难以处理大规模散射问题的瓶颈。以全空间孔洞对P、SV波的二维散射为例,给出了具体求解步骤,并在个人计算机上实现了上百万自由度问题的快速精确计算。在方法效率和精度检验基础上,分别以单孔洞和随机孔洞群对平面波（P、SV波）的散射为例进行计算模拟,揭示了孔洞（群）周围弹性波散射的若干重要规律。     

Scattering of SH-waves by an interface cavity

Summary. The scattering of the SH-wave and dynamic stress concentrations near an arbitrary cavity situated at the planar interface separating two different elastic media are investigated. The total wave field can be obtained by superposition of the free field and the scattered field. The free field is composed of the incident, reflected and refracted waves. The scattered wave fields in adjacent media are expressed respectively, and the method of wave functions expansion is applied to obtain the solutions for these fields. The scattered wave functions can be expanded into Hankel-Fourier series with unknown coefficients. In solving for the unknown coefficients according to the boundary conditions for the total wave field at the interface and at the cavity wall, the non-orthogonality makes the system of equations for the unknown coefficients infinite and coupling each other. Another key point is to extend each scattered wave field from its own half-plane domain into the full plane domain by a certain way keeping the total wave field unchanged for the non-orthogonal Fourier integrals around the cavity. Finally, the scattering of the SH wave by an interface ellipse with different ratios between long and short axis is considered, and the distributions of dynamic stress concentration factors at the cavity wall are presented.     

Diffraction of a stress wave by a cylindrical cavity in an infinite elastic strip

The scattering of an incident plane shock wave by a cylindrical circular cavity in an infinite elastic strip is considered. In the Laplace transformed domain, boundary conditions at the plane surfaces and those at the circular hole are satisfied with the help of the Fourier transformation and the Schmidt method. A numerical Laplace inversion technique is taken to obtain the stresses in the physical space.     

A study on photothermal waves in an unbounded semiconductor medium with cylindrical cavity

The theory of coupled plasma, thermal, and elastic waves was used to investigate the wave propagation on semiconductor material with cylindrical cavity during photo-thermoelastic process. An unbounded material, elastic semiconductor containing a cylindrical cavity with isotropic and homogeneous thermal and elastic properties has been considered. The inner surface of cavity is constrained, and the carrier density is photogenerated by an exponentially decaying pulse boundary heat flux. The eigenvalue approach, together with Laplace transform techniques, was used to obtain the analytical solutions. Numerical computations have been done for a silicon-like semiconductor material, and the results are presented graphically to estimate the effect of the coupling between the plasma, thermal, and elastic waves. The graphical results indicate that the thermal activation coupling parameter is an important phenomenon and has a great effect on the distribution of field quantities.     

开孔薄板弹性波散射与动应力集中

本文采用边界无法对开孔无限大薄板弹性波的散射与动应力集中问题进行理论分析和数值计算。基于动力学功的互等定理建立了薄板弯曲波动问题的边界积分方程,应用Ｍａｔｈｅｍａｔｉｃａ软件首次推导了各影响系数的计算公式．最后,给出了圆孔附近的动应力集中系数的数值结果。     

Three‐dimensional BEM for scattering of elastic waves in general anisotropic media

Based on the full‐space Green's functions, a three‐dimensional time‐harmonic boundary element method is presented for the scattering of elastic waves in a triclinic full space. The boundary integral equations for incident, scattered and total wave fields are given. An efficient numerical method is proposed to calculate the free terms for any geometry. The discretization of the boundary integral equation is achieved by using a linear triangular element. Applications are discussed for scattering of elastic waves by a spherical cavity in a 3D triclinic medium. The method has been tested by comparing the numerical results with the existing analytical solutions for an isotropic problem. The results show that, in addition to the frequency of the incident waves, the scattered waves strongly depend on the anisotropy of the media. Copyright © 2003 John Wiley & Sons, Ltd.     

带形介质内SH型导波对圆柱孔洞的动力分析

弹性板壳的反平面运动中由缺陷引起的应力集中问题可以按照带形介质中圆柱孔洞对SH型导波的散射问题来分析。首先,构造带形介质中相容导波的形式,即其满足上、下边界应力自由条件。之后,由波函数展开法给出圆柱孔洞散射波的级数表示,根据累次镜像叠加的方法构造由上、下边界反射所形成的相容的散射导波。最后,给定入射导波,由圆柱孔洞边界应力自由的条件来定解波函数级数的系数。数值算例求解了特定导波对圆柱孔洞的散射,给出了圆柱孔洞边沿的动应力分布,讨论了导波阶数、频率以及孔洞位置的影响。     

A method for linear elasto-static analysis of multi-layered axisymmetrical bodies using Hankel's transform

 A method to determine the distribution of stresses and displacements in an infinite, linear, elastic, multi-layered medium subjected to static axisymmetric loading is presented in this work. By using axisymmetric governing equations, Hankel's transform and matrix analysis, the methodology gives a clearly arranged way to calculate the stresses and displacements in the medium. A numerical method for Hankel's transform is employed to perform the calculation. Two representative examples are studied. The results can be utilized as a fundamental solution for boundary element methods for the linear, elasto-static, axisymmetric multi-layered problem with a little modification.     

A sweeping window method for detection of flaws using an explicit dynamic XFEM and absorbing boundary layers

We present a sweeping window method in elastodynamics for detection of multiple flaws embedded in a large structure. The key idea is to measure the elastic wave propagation generated by a dynamic load within a smaller substructural detecting window domain, given a sufficient number of sensors. Hence, rather than solving the full structure, one solves a set of smaller dynamic problems quickly and efficiently. To this end, an explicit dynamic extended FEM with circular/elliptical void enrichments is implemented to model the propagation of elastic waves in the detecting window domain. To avoid wave reflections, we consider the window as an unbounded domain with the option of full‐infinite/semi‐infinite/quarter‐infinite domains and employ a simple multi‐dimensional absorbing boundary layer technique. A spatially varying Rayleigh damping is proposed to eliminate spurious wave reflections at the artificial model boundaries. In the process of flaw detection, two phases are proposed: (i) pre‐analysis—identification of rough damage regions through a data‐driven approach, and (ii) post‐analysis‐–identification of the true flaw parameters by a two‐stage optimization technique. The ‘pre‐analysis’ phase considers the information contained in the ‘pseudo’ healthy structure and the scattered wave signals, providing an admissible initial guess for the optimization process. Then a two‐stage optimization approach (the simplex method and a damped Gauss–Newton algorithm) is carried out in the ‘post‐analysis’ phase for convergence to the true flaw parameters. A weighted sum of the least squares, of the residuals between the measured and simulated waves, is used to construct the objective function for optimization. Several benchmark examples are numerically illustrated to test the performance of the proposed sweeping methodology for detection of multiple flaws in an unbounded elastic domain. Copyright © 2015 John Wiley & Sons, Ltd.     

Elastodynamic guided wave scattering in infinite plates

The scattering of elastic guided waves by defects in two‐dimensional infinite plates is analysed in the plane and antiplane cases, corresponding, respectively, to Lamb and SH modes. A hybrid boundary element–finite element technique is used, where the defect neighbourhood is discretized with quadratic boundary elements and the radiation condition in the plate is satisfied through a normal mode expansion. A semi‐analytical finite element technique is applied in the infinite plate to calculate its dispersion curves and normal modes. This hybrid technique, which showed excellent performance in the solution of Lamb wave reflection at the edge of semi‐infinite plates, is extended in this paper to a wider range of problems, such as Lamb mode scattering by delaminations and surface defects and SH mode interaction with step discontinuities. Copyright © 2003 John Wiley & Sons, Ltd.     

Theoretical and experimental investigations of lateral modes in 1-3piezocomposites

Materials with a periodic microstructure show resonances caused by the elastic wave Bragg diffraction. This paper presents a simple approach to describe these resonances (called lateral resonances) in 1-3 piezoelectric composite materials which have a 2-D periodicity. Our model is based on the analysis of the propagation of transverse waves in a 2-D periodic medium of infinite thickness and takes into account the periodic and interfacial boundary conditions. This model predicts the displacement field vectors and frequencies of lateral resonances from which the phase velocity of lateral waves is determined. The theoretical and experimental variations of this velocity versus the ceramic rod width to pitch ratio are compared. It is shown that the first lateral mode frequency is maximum when the ceramic volume fraction is around 0.65. Theoretical predictions of the mechanical displacement at the composite surface are compared with measurements obtained by an interferometric laser technique. A good agreement is observed, showing that lateral waves are mainly vertically polarized     

Reflection and refraction of plane waves at interface between two piezoelectric media

This paper analyses reflection and refraction of plane waves at a perfect interface between two anisotropic piezoelectric media. The equations of elastic waves, quasi-static electric field, and constitutive relationships for the piezoelectric media are derived. A solution based on the inhomogeneous wave theory is developed to address the inconsistency between the numbers of independent wave modes in the media and the numbers of interfacial boundary conditions to obtain accurate reflection and refraction coefficients in case of strong piezoelectric media, where all the elastic and electric continuity conditions across the interface are satisfied simultaneously. The study shows that there exist independent and zero energy wave modes satisfying the general Snell’s law and propagating along the interface for any incident wave angle. These waves can be treated as pseudo surface waves. It is further found that all the reflection/refraction waves including the pseudo surface waves obey the energy conservation law at the interface boundary. In addition, the analysis also reveals that the reflection and refraction elastic waves can turn into pseudo surface waves at some critical incident angles.     

反平面冲击荷载作用下圆柱形洞室的动力响应

对无限弹性土体内圆柱形洞室在反平面冲击荷载的作用下的瞬态响应进行了探讨。用残余变量的方法结合拉普拉斯变换及其逆变换,在频域内的得到土体位移和应力一般表达式,并采用拉普拉斯逆变换的数值方法,给出了问题的数值解。在时域内分析了无限弹性土体内圆柱形孔洞在沿轴向方向的冲击荷载作用下的土体动力响应的变化规律,并将计算结果与静力情况下的作了比较：反平面荷载作用后,波向外发散传播,并沿半径方向衰减,衰减速度较静力情况下的衰减要慢,且距离波源越远二者差别越大;波到达后,该点土体的应力和位移均发生与外荷载相同的三角形变化：先增大,随后减小最后保持为零不变。     

A spectral alternating method for elastostatic problems with multiple spherical cavities

The problem of an infinite elastic solid containing an arbitrary number of non-overlapping spherical cavities of arbitrary sizes and locations and with arbitrary boundary tractions is considered. The numerical procedure to solve this boundary value problem is based on a spectral method in which the boundary displacements and tractions are represented by truncated series of surface spherical harmonics. By using an alternating method, the problem containing multiple spherical cavities is replaced by a sequence of problems for a single spherical cavity with the boundary conditions adjusted iteratively to account for the cavity interactions. A least squares approximation is used to determine the unknown coefficients of surface spherical harmonics representing these cavity interactions. Several examples are given to illustrate the approach. An erratum to this article can be found at     

Elastic waves in a fluid-loaded, semi-infinite axisymmetric rod

In a fluid-loaded, semi-infinite axisymmetric rod, a free shear stress boundary condition on the circular cross-sectional end introduces complicated, nondispersive waves in the solid. They are composed of a pulse wave, which has the same waveform as the transmitted one and travels at speed c1, and different kinds of pulse trains, each of which travels along the rod at the speed of either c1 or square root of 2c2, where c1 and c2 are the propagating speeds of the longitudinal and transversal bulk waves, respectively. Furthermore, one can conclude from the solutions to the boundary conditions that c1 and square root of 2c2 are the only phase speeds of nondispersive waves. Frequency equations associated with these waves are established, and the solutions are solved and discussed analytically and numerically. The acoustic field in the fluid is also fully discussed, and it is more complicated than a single outgoing Hankel function as described for an infinite rod. The acoustic energy coupling between the solid and the fluid and the end reflection and transmission are quantified as well. In the end, experimental examinations of the echo spectra, using an aluminum rod immersed in the water and air, fully confirm the numerical solutions to the frequency equations.     

Non-reflecting boundary conditions for plane waves in anisotropic elasticity and poroelasticity

The paper presents exact non-reflecting boundary conditions for transient plane waves in an anisotropic elastic solid for oblique incidence. The boundary conditions are expressed through the eigenvectors of the acoustic tensor and are written in impedance form as a relation between the velocity vector and the traction vector. The approach is extended to anisotropic fluid-saturated porous solids. Exact plane-wave non-reflecting boundary conditions are derived for transient non-dissipative waves in a medium with infinite or zero permeability, and for steady-state dissipative waves.     

Expansion formulae for wave structure interaction problems with applications in hydroelasticity

Alternate derivations of the expansion formulae for wave structure interaction problems are obtained in case of water of infinite depth and utilized to analyze the hydroelastic behavior of large floating structures. Considering the boundary value problem associated with Laplace equation having higher order boundary condition on the horizontal boundary and a Dirichlet type boundary condition on the vertical boundary in a quarter plane, Fourier sine transform is applied in the horizontal direction to convert the problem to a Sturm-Liouville type boundary value problem associated with non-homogeneous ordinary differential equation (ODE) in the transformed variable. Finally, inverting the transformed functions and applying the regularity criterion of the transformed function, the required expansion formula is derived. The expansion formula thus derived is extended to deal with similar boundary value problems having Neumann type boundary condition. The expansion formulae are applied to (i) analyze oblique scattering of flexural gravity waves by an articulated floating elastic plate and (ii) study the effect of compression on the oblique scattering of flexural gravity waves by a line discontinuity in a large floating ice sheet in water of infinite depth, which find applications in marine technology and arctic engineering, respectively. The present derivations of the expansion formulae are very simple and straightforward and can be easily used to study a large class of problems in the area of fluids and structures in mathematical physics and engineering.