非均匀介质散射问题的体积分方程数值解法

将非均匀介质视为某一均匀背景介质中的扰动，可建立用均匀背景介质格林函数作基本解的体积分方程．给出了配置法求解体积分方程的数值方法，首先解得扰动域内各点以速度扰动为权的波场函数，然后回代计算得到观测面上各接收点的散射波场．与边界元法和Ｂｏｒｎ近似法计算结果比较表明该方法具有很高的精度，可得到穿过非     

带形域内裂纹对SH波散射问题的边界元解

采用边界元法研究含裂纹的带形域各向同性弹性体，裂纹对SH波的散射问题，推导出带形域情况下不同边界条件的各种Green函数，导出了以裂纹张开位移为未知数的边界积分方程，计算出表面散射场和总位移．算例表明，利用所提供的格林函数和边界元格式解答带形域的散射问题比较方便．     

三维弹性问题Taylor展开多极边界元法的误差分析

Taylor展开多极边界元法有效的提高了边界元法的求解效率,使之可用于大规模问题的计算。然而,由于计算中对基本解进行了Taylor级数展开,与传统边界元方法相比计算精度有所下降。本文主要针对三维弹性问题Taylor展开多极边界元法的计算精度和误差进行研究。文中对两种方法的计算精度进行了比较;研究了核函数的Taylor展开性质;推导了三维弹性问题基本解的误差估计公式;给出了Taylor展开多极边界元法中远近场的划分原则。通过具体的算例,证明了该方法的正确性和误差估计公式的有效性,说明了影响Taylor展开多极边界元法求解精度的因素。     

带形域内裂编译地SH波散射问题的边界元解

采用边界元法研究含裂纹的带形域各向性弹性体，裂纹对ＳＨ波的散射问题，推导出带形域情况下不同边界条件的各种Ｇｒｅｅｎ函数，导出了以裂纹张开位移为未知数的边界积分方程，计算出表面散射场和总位移，算例表明，利用所提供的格林函数和边界元格式解答带形域的散射问题比较方便     

二维直角平面内固定圆形夹杂对稳态入射反平面剪切波的散射

利用复变函数法、多极坐标及傅立叶级数展开技术求解了二维直角平面内固定圆形夹杂对稳态入射反平面剪切（shearing horizontal, SH）波的散射问题。首先构造出介质内不存在夹杂时的入射波场和反射波场,然后建立介质内存在夹杂时由夹杂边界产生的能够自动满足直角边应力自由条件的散射波解,从而利用叠加原理写出介质内的总波场。利用夹杂边界处位移条件和傅立叶级数展开方法列出求解散射波中未知系数的无穷代数方程组,在满足计算精度的前提下通过有限项截断,得到相应有限代数方程组的解,最后通过算例具体讨论了二维直角平面水平边界点的位移幅度比和相位随量纲一波数、入射波入射角及夹杂位置的不同而变化的情况,结果表明了算法的有效实用性。     

一种用于计算板中缺陷兰姆波散射的混合边界元法

本文介绍了一种将传统边界元法和兰姆波的本征模式函数相结合的混合边界元法，用此方法解决兰姆波在板材中传播时遇到缺陷发生的散射问题，计算了兰姆波通过缺陷后的反射与透射系数。文中给出了若干兰姆波信号幅度与板中不同深度缺陷的相互关系。     

楔形空间中任意形状孔洞对平面SH波的散射：IBEM求解

采用间接边界元法(IBEM),求解了楔形空间中任意形状孔洞对平面SH波的散射问题。基于单层位势理论,首先在孔洞表面及其附近楔形空间表面上施加虚拟均布荷载,以构造散射波场。进而由自由表面零应力条件,建立方程求解得到虚拟荷载密度。总波场由楔形空间自由波场和散射波场叠加得到。研究表明：IBEM方法能够精确高效求解楔形空间中弹性波散射问题。楔形空间孔洞对波的散射特征依赖于波入射角、无量纲入射波频率、楔形夹角、孔洞位置及其形状;孔洞周围波的相干效应十分显著,空间表面位移幅值及孔边动应力集中因子比半空间情况放大一倍有余,分别达到8.5和10.0;该研究为楔形空间中更为复杂的P、SV波散射问题求解奠定了基础。     

直角平面内弹性圆夹杂对入射平面SH波的散射

利用复变函数法、多极坐标移动技术及傅立叶级数展开求解二维直角平面内圆形弹性夹杂对稳态入射平面SH波的散射问题。首先写出直角平面内不含夹杂时的入射波场和反射波场;其次建立直角平面内含夹杂时夹杂外的散射波解和夹杂内的驻波解,并利用叠加原理写出问题的总波场,借助夹杂边界处应力和位移的连续条件建立求解散射波解和驻波解中未知系数的无穷代数方程组并求解,通过算例具体讨论了直角平面水平边界点的位移幅度比和夹杂边界处径向应力集中系数随不同无量纲波数、入射角及圆孔位置的变化情况,结果表明了算法的有效实用性。     

杆件结构强迫振动分析时域基本解

基本解是边界元法、格林函数法、基本解法等数值方法必备的理论基础之一.大多问题的基本解已经被推导出来,但杆件结构强迫振动分析的时域基本解未见相关报道.现有方法在求解强迫振动问题时通常采用一定技巧将问题等效为静力问题,并利用静力基本解进行求解,这就增加了列式过程的复杂程度和方法的计算工作量.论文通过振型分解、杜哈梅积分等步骤推导出杆件结构强迫振动分析的时域基本解,该基本解可以直接用于结构的动力响应分析,提高运算效率.论文工作可以为边界元法等数值方法在进行杆系结构强迫振动分析时提供新的思路.     

直角平面区域内固定圆形刚性夹杂问题的Green函数解

利用复变函数法、多极坐标移动技术研究了直角平面区域内含有固定圆形夹杂时的反平面问题Green函数解.首先构造出不含夹杂的完整直角平面区域内满足边界应力条件的入射位移场;其次,建立直角平面区域内固定圆形夹杂对该入射场产生的满足直角边界应力自由条件的散射波解,并由叠加原理得到介质内的总波场.最后利用夹杂边界处的位移条件确定出散射波解中的未知系数,最终得到问题的Green函数解,还通过算例讨论了夹杂边界处的径向应力和环向应力随不同波数、角度和不同夹杂位置及不同点源位置的变化情况.算例结果表明了该文方法的有效实用性.     

Analysis of scattering waves in an elastic layered medium by means of the complete eigenfunction expansion form of the Green’s function

A scattering problem due to an object and a plane incident wave in an elastic layered half space is presented in this paper. The complete eigenfunction expansion form of the Green’s function developed by the author and the boundary integral equation method are introduced into the analysis. First, the complete eigenfunction expansion form of the Green’s function is investigated for its application to the scattering problem. A comprehensive explanation is also given for the fact that the complex Rayleigh wave modes exhibit standing waves. Next, a method for the analysis of scattering waves by means of the Green’s function is presented. The advantage of the present method is that the formulation itself is independent of the number of layers and the scattering waves can be decomposed into the modes for the spectra defined for the layered medium. Several numerical calculations are performed to examine the efficiency of the present method as well as the properties of the scattering waves. According to the numerical results, the complete eigenfunction expansion form of the Green’s function provides accurate values for application to a boundary element analysis. The spectral structure and radiation patterns of the scattering wave are presented and investigated. The differences in directionality can be found from the radiation patterns of the scattering waves decomposed into the modes for the spectra.     

Scattering of elastic waves by an arbitrary small imperfection in the surface of a half-space

T , the first of two articles, is concerned with the scattering of elastic waves by arbitrary surface-breaking or near surface defects in an isotropic half-plane. We present an analytical solution, by the method of matched asymptotic expansions, when the parameter , which is the ratio of a typical length scale of the imperfection to the incident radiation's wavelength, is small. The problem is formulated for a general class of small defects, including cracks, surface bumps and inclusions, and for arbitrary incident waves. As a straightforward example of the asymptotic scheme we specialize the defect to a two-dimensional circular void or protrusion, which breaks the free surface, and assume Rayleigh wave excitation ; this inner problem is exactly solvable by conformal mapping methods. The displacement field is found uniformly to leading order in , and the Rayleigh waves which are scattered by the crack are explicitly determined. In the second article we use the method given here to tackle the important problem of an inclined edge-crack. In that work we show that the scattered field can be found to any asymptotic order in a straightforward manner, and in particular the Rayleigh wave coefficients are given to O(²).     

Application of the reciprocity theorem to scattering of surface waves by a cavity

Scattering of surface waves by a cylindrical cavity at the surface of a homogenous, isotropic, linearly elastic half-space is analyzed in this paper. In the usual manner, the scattered field is shown to be equivalent to the radiation from a distribution of tractions, obtained from the incident wave on the surface of the cavity. For the approximation used in this paper, these tractions are shifted to tractions applied to the projection of the cavity on the surface of the half-space. The radiation of surface waves from a normal and a tangential line load, recently determined by the use of the reciprocity theorem, is employed to obtain the field scattered by the cavity from the superposition of displacements due to the distributed surface tractions. The vertical displacement at some distance from the cavity is compared with the solution of the scattering problem obtained by the boundary element method (BEM) for various depths and widths of the cavity. Comparisons between the analytical and BEM results are graphically displayed. The limitations of the approximate approach are discussed based on the comparisons with the BEM results.     

Scattering of elastic waves by a small inclined surface-breaking crack

I we examine the scattering of Rayleigh waves by an inclined two-dimensional plane surface-breaking crack in an isotropic elastic half-plane. We follow the method already introduced by the authors (A and W , 1992a, J. Mech. Phys. Solids 40, 1683) to obtain an analytical solution when the parameter , which is the ratio of a typical length scale of the imperfection to the incident radiation's wavelength, is small. The procedure employed is the method of matched asymptotic expansions, which involves determining the scattered wave field both away from and near the crack. The outer solution is specialized from the general expansion in the first part of this work (A and W , 1992a, J. Mech. Phys. Solids 40, 1683), and the inner problem is exactly solved by the Wiener-Hopf technique. The displacement field and scattered Rayleigh waves are found uniformly to third order in , and concluding remarks are made about the general method as well as the results presented here.     

多个椭圆柱波浪力的一种解析解

波浪在大尺寸结构表面产生不可忽略的散射波, 该散射波在多柱体体系中继续传播, 并在同体系中的其他柱体上产生高次散射波. 本文基于椭圆坐标系和绕射波理论首先推导了波浪作用下椭圆单柱体产生的散射波压力公式, 随后考虑该散射波在多柱体系中的传播, 将其视为第二次入射波, 推导出柱体上第二次散射波压力公式, 同理可以推导出高次散射波压力公式, 最后得到椭圆多柱体波浪力解析解, 并用数值解验证了本文解析方法的正确性. 本文以双柱体和四柱体体系为例, 分析了不同参数(波数、净距、波浪入射角度等)下, 高次散射波对柱体上波浪作用的影响. 结果表明: 波数较大的情况下, 高次散射波引起柱体上的波浪力不能忽略; 结构间距较大的情况下, 虽然高次波的作用有减小的趋势但仍然明显; 高次散射波来自多个柱体对入射波的散射, 柱体数目的增加后, 高次波的影响会增加, 结构所受的高次波作用因参数变化而起的波动会变剧烈; 高次波对上游柱体波浪力的贡献较对下游柱体的贡献大.      

Propagation of surface waves in a half-space with vertical,inclined or curved interfaces

The paper reviews the author's theoretical studies on reflection and transmission of surface waves at vertical and near-vertical boundaries. For estimation of reflection and transmission coefficients a new approximate method based on the Green function technique is proposed. The method is valid if a body waves field arised at the boundary is approximately orthogonal to the fields of reflected and transmitted normal modes. It may be used for normal as well as for oblique incidence. A Green function in each quarterspace is assumed to be the same as in a halfspace with the same vertical distribution of elastic parameters. For large angles of incidence when the reflected wave field is strong enough, the Green function may be corrected for the reflected wave. The results of calculation for some models by means of this method are discussed.     

Modeling nonlinear Rayleigh wave fields generated by angle beam wedge transducers—A theoretical study

Nonlinear Rayleigh wave fields generated by an angle beam wedge transducer are modeled in this study. The calculated area sound sources underneath the wedge are used to model the fundamental Rayleigh sound fields on the specimen surface, which are more accurate than the previously used line sources with uniform or Gaussian amplitude distributions. A general two-dimensional nonlinear Rayleigh wave equation without parabolic approximation is introduced and the solutions are obtained using the quasilinear theory. The second harmonic Rayleigh wave due to material nonlinearity is given in an integral expression with these fundamental Rayleigh waves radiated by the wedge transmitter acting as a forcing function. Multi-Gaussian beam (MGB) models are employed to simplify these integral solutions and to extract the diffraction and attenuation correction terms explicitly. The effect of nonlinearity of generating sources on the second harmonic Rayleigh wave fields is taken into consideration; simulation results show that it will affect the magnitude and diffraction correction of the second harmonic waves in the region close to the Rayleigh wave sound sources. This research provides a theoretical improvement to alleviate the experimental restriction on analyzing the effects of diffraction, attenuation and source nonlinearity when using angle beam wedge transducers as transmitters.     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

Full field analysis of an anti-plane interface crack subjected to dynamic body forces

In this study, the transient full field response of an interface crack between two different media subjected to dynamic body force at one material is investigated. For time t < 0, the bimaterial medium is stress free and at rest. At t = 0, a concentrated anti-plane dynamic point loading is applied at the medium as shown in Fig. 1. The total wave field is due to the effect of this point loading and the scattering of the incident waves by the interface crack. An alternative methodology that is different from the conventional superposition method is used to construct the reflected, refracted and diffracted wave fields. A useful fundamental solution is proposed in this study and the full field solution is determined by superposition of the fundamental solution in the Laplace transform domain. The proposed fundamental problem is the problem of applying an exponentially distributed traction (in the Laplace transform domain) on the interfacial crack faces. The Cagniard–de Hoop method of Laplace inversion is used to obtain the transient solution in time domain. Exact transient closed form solutions for stresses and stress intensity factors are obtained. Numerical results for the time history of stresses and stress intensity factors during the transient process are discussed in detail.     

Scattering of guided waves by circumferential cracks in composite cylinders

A somewhat generalized numerical procedure is used in this paper to study the problem of wave scattering by circumferential cracks in composite pipes. The study is motivated by the need to develop a model for the quantitative, ultrasonic non-destructive evaluation of cracks in pipes. For this purpose, a stiffness-based Rayleigh–Ritz type approach is employed first to obtain the approximate wave numbers and wave modes. Using the wave function expansions of the incident and scattered fields in the axial direction and decomposing the problem into separate symmetric and anti-symmetric problems, a three-dimensional wave scattering problem is reduced to two, independent two-dimensional problems over the circular cross-section. Both these problems can be reduced further to quasi-one-dimensions by discretizing the cross-section into finite elements and using a transfer matrix approach in the circumferential direction. This simplification greatly reduces the computational time. A comparison of the results for an isotropic pipe demonstrates the reliability and accuracy of the modified numerical procedure. Numerical results for the reflection and transmission coefficients of different incident wave modes are also presented for a 2-ply composite pipe with a crack. The crack may have an arbitrary circumferential length and radial depth. Simple extrapolations from one wave to another wave, separately incident on a crack, are demonstrated to be impossible due to different mode conversions by the crack.