半无限空间中圆形孔洞周围SH波的散射

建立了求解在含有圆形孔洞的弹性半空间中SH波散射与圆形孔洞附近动应力集中问题的解析方法。利用SH波散射的对称性和多极坐标的方法，构造了一个可以预先满足半空间自由表面上应力自由边界条件的圆形孔洞对SH波散射的波函数。利用这一波函数，则可将该问题转化成对一个圆形孔洞散射的求解问题。该问题的解答最终又可归结为对一组无究代数方程组的求解问题，并可利用截断有限项的方法对其进行计算，最后给出了有关圆形也洞附近动应力集中问题的算例和数值结果，并讨论了波数与圆孔至自由边界距离变化对动应力集中的影响。     

半无限空间界面附近SH波对圆形衬砌的散射

建立了求解半无限空间中SH波对浅埋圆形衬砌结构的散射与动应力集中问题的解析方法。利用SH波散射的对称性和多极坐标的方法，在复平面上构造出了一个可以预先满足半空间自由表面上应力自由的边界条件的浅埋圆形衬砌对稳态SH波散射的波函数，并构造出衬砌内的散射波函数。然后根据衬砌周围的边界条件，将该问题转化为对一组无穷代数方程组的求解。最后给出了具体算例，并讨论了其数值结果。     

SH波冲击下浅埋任意形孔洞的动力分析

求解了稳态SH波垂直弹性半空间水平表面入射时,浅埋任意形孔洞的动力响应。采用复变函数和多极坐标方法构造了一个能够自动满足水平表面上应力自由边界条件的散射波函数。应用这一波函数,将半空间中的问题转化为求解一个全空间中任意形孔洞的散射问题,最终将问题归结为对一组无穷代数方程组的求解。作为对抗爆问题的研究,给出了浅埋椭圆孔和方孔附近的动应力集中系数的数值结果,并对算例进行了讨论。     

SH波在浅埋可移动圆柱形刚性夹杂处的散射与地震动

采用复函数、多极坐标法研究了含有可移动圆柱形刚性夹杂的弹性半空间对SH波的散射问题。构造一个能自动满足含可移动刚性圆柱的弹性半空间自由表面上应力为零的边界条件的散射波,应用可移动圆柱形刚性夹杂的运动条件来确定该散射波。最终则可将求解问题归结为求解一个无穷代数方程组,采用截断有限项的方法对其进行求解。给出了地表位移幅值的数值结果,并讨论了各种参数对它的影响。     

半无限空间中稳态P波在衬砌周围的散射

本文采用多级坐标和复变函数方法针对半无限空间中P波在一衬砌周围的散射问题提出了一种近似求解分析方法。具体做法是利用一个半径很大的圆来逼近半空间的直边界，将待解问题转化为全空间中P波在一圆孔和一衬砌周围的共同散射问题。我们预先写出问题波函数的一般形式解，利用边界条件并借助复数傅立叶级数展开把问题化为求解波函数中未知系数的无穷线性代数方程组，进而讨论了衬砌内外边界处动应力集中系数针对不同条件组合的分布和变化情况。算例结果表明：本文方法对于研究与P波（SV波）有关的散射问题是可行的。     

各向异性半空间中浅埋孔洞对地表 反平面运动的影响

本文采用复变函数方法研究了稳态水平剪切波(SH波)在各向异性弹性半空间中任意形状孔洞上的散射及其对地面运动的影响.在变换空间中构造出自动满足各向异性半空间水平表面上应力为零的散射波函数,并根据孔洞表面应力为零的边界条件,采用最小二乘法求解散射波函数的系数.用介质的各向异性性质来模拟地质条件,给出了SH波作用下含圆孔、椭圆孔和方孔的各向异性半空间表面位移幅值的数值结果,并分析了介质特性、孔洞的形状、埋深、入射波波数及入射角度等因素对地面运动的影响规律.数值结果表明:介质的各向异性对含有孔洞的半空间表面的地表位移具有显著的影响;沿一定角度的入射波在某一频段内所引起的地表位移幅值比各向同性介质的可能要大,且随着孔洞埋深的增加,地表位移的幅值逐渐减小.      

用波函数展开法求解界面圆孔的SH波散射问题

研究了ＳＨ波绕界面圆孔散射的波函数展开理论。由入射,反射和透射波组成的自由波场与孔的散射场叠加成总波场。将定义于两个半平面的散射场按一定方式延拓到全平面。通过Ｈａｎｋｅｌ－Ｆｏｕｒｉｅｒ型积分来满足界面与孔缘处的边界条件,得到了确定特定系数的封闭代数方程组。给出了具体算例,计算了孔缘处的动应力系数分布。     

Scattering of SH waves by a scalene triangular hill with a shallow cavity in half-space

In this study, a theoretical approach is used to investigate the scattering problem of circular holes under a scalene triangle on the surface. The wave displacement function is obtained by solving the Helmholtz equation that meets the zero-stress boundary conditions by adopting the method of separation of variables. Based on the complex function, multi-polar coordinate method, and region-matching technique, algebraic equations are established at auxiliary boundaries and free boundaries conditions in a complex domain. The auxiliary circle is used to solve the singularity of the reflex angle at the triangle corner. Then, according to sample statistics, the least squares method is used instead of the Fourier expansion method to solve the undetermined coefficient of the algebraic equations by discrete boundary. Numerical results show that the continuity of the auxiliary boundaries and the accuracy of the zero-stress boundaries are adequate, and the displacement of the free surface and the stress of the circular hole are related to the shape of the triangle, the position of the circular hole, the direction of the incident wave, and the frequency content of the excitation. Finally, time-domain responses are calculated by FFT based on the frequency domain theory, and the results reveal the wave propagation mechanism in a complicated structure.

     

SH波入射时半圆形凸起地形附近浅埋圆形衬砌结构的动应力分析

本文给出了地下圆形衬砌结构与地面上的半圆形凸起地形对垂直于地面入射的SH波散射问题的解答。方法是将求解区域分割成两部分。其一为包含半圆形凸起地形在内的圆形区域Ⅰ，其二为带有一个半圆形凹陷和一个圆形衬砌结构的弹性半空间Ⅱ，半圆形凹陷部分为其公共边界，在区域Ⅰ和Ⅱ中分别构造其位移解，然后再通过移动坐标，使其满足“公共边界”上的条件和地下衬砌的边界条件，建立起求解该问题的无穷代数方程组。最后，本文给出了算例，并讨论了数值结果，给出了圆形衬砌结构周边上的动应力集中系数变化规律。     

SH波入射时多个半圆形沉积谷地附近浅埋圆形孔洞的动力分析

本文采用辅助函数的思想,利用复变函数和多级坐标的方法给出了SH波入射条件下多个半圆形沉积谷地附近浅埋圆形孔洞动力分析问题的解答。将整个求解区域分割成两部分来处理,区域I为多个半圆形沉积谷地,区域II为浅埋圆形孔洞附近带半圆形凹陷的半无限弹性空间。在区域I和II中分别构造位移解,并在二个区域的“公共边界”上实施位移应力的连续条件,建立求解该问题的无穷代数方程组。最后,本文给出了算例和数值结果,并对其进行了讨论。     

SH-wave diffraction by a semi-circular hill revisited: A null-field boundary integral equation method using degenerate kernels

Following the success of seismic analysis of a canyon [1], the problem of SH-wave diffraction by a semi-circular hill is revisited using the null-field boundary integral equation method (BIEM). To fully utilize the analytical property in the null-field boundary integral equation approach in conjunction with degenerate kernels for solving the semi-circular hill scattering problem, the problem is decomposed into two regions to produce circular boundaries using the technique of taking free body. One is the half-plane problem containing a semi-circular boundary. This semi-infinite problem is imbedded in an infinite plane with an artificial full circular boundary such that degenerate kernel can be fully applied. The other is an interior problem bounded by a circular boundary. The degenerate kernel in the polar coordinates for two subdomains is utilized for the closed-form fundamental solution. The semi-analytical formulation along with matching boundary conditions yields six constraint equations. Instead of finding admissible wave expansion bases, our null-field BIEM approach in conjunction with degenerate kernels have five features over the conventional BIEM/BEM: (1) free from calculating principal values, (2) exponential convergence, (3) elimination of boundary-layer effect, (4) meshless and (5) well-posed system. All the numerical results are comparing well with the available results in the literature. It is interesting to find that a focusing phenomenon is also observed in this study.     

Scattering of plane P,SV or Rayleigh waves by a shallow lined tunnel in an elastic half space

Based on the plane complex variable theory and the image technique, an analytical solution is presented for scattering of plane harmonic P, SV or Rayleigh waves by a shallow lined circular tunnel in an elastic half space. The major contribution of this study is the treatment of the orthogonality of the boundary conditions along the half surface and the cavity wall. In terms of the image technique, the scattered waves by the half surface are simulated as transmitting from the image source of the origin of the tunnel. Using two different conformal mapping functions, we obtained the complex-valued stresses and displacements of the elastic medium and the liner in the image domain, respectively. The boundary value problem results in a set of infinite algebraic equations. The accuracy of the present approach is verified by comparing the present solution results with the available published data. Parametric study indicates that the embedment depth, the shear modulus and the thickness of the liner have significant influences on the dynamic response of the liner and the medium.     

曲线坐标在弹性波散射研究中的应用──SH波对不等深度凹陷地形的散射

本文利用曲线坐标的分析方法研究ＳＨ波对凹陷地形散射的边值问题。首先给出了在曲线坐标系中波动方程解的逼近函数序列，并使此解满足由曲线坐标所描述的具有不同深度的凹陷地形上的边界条件，最后将待解的问题归结为对一组无穷代数方程组的求解．本文提供的方法有如下的特点：（１）可将通常求解凹陷地形对ＳＨ波散射的波动方程的混合边值问题（指在不同边界区域上给出不同的应力分量的边界条件）化为一应力边值问题；（２）本文所使用的曲线坐标，可描述具有不同深度与宽度比的拟半圆形凹陷地形散射问题，并给出一般解；（３）本文方法可推广用于求解Ｐ（－ＳＶ）波对凹陷地形的散射问题。作为算例，本文给出了ＳＨ波对具有不同深度与宽度比的凹陷地形散射的数值结果，并进行了分析比较。     

Focusing phenomenon and near-trapped modes of SH waves

In this study, the null-field boundary integral equation method (BIEM) and the image method are used to solve the SH wave scattering problem containing semi-circular canyons and circular tunnels. To fully utilize the analytical property of circular geometry, the polar coordinates are used to expand the closed-form fundamental solution to the degenerate kernel, and the Fourier series is also introduced to represent the boundary density. By collocating boundary points to match boundary condition on the boundary, a linear algebraic system is constructed. The unknown coefficients in the algebraic system can be easily determined. In this way, a semi-analytical approach is developed. Following the experience of near-trapped modes in water wave problems of the full plane, the focusing phenomenon and near-trapped modes for the SH wave problem of the half-plane are solved, since the two problems obey the same mathematical model. In this study, it is found that the SH wave problem containing two semi-circular canyons and a circular tunnel has the near-trapped mode and the focusing phenomenon for a special incident angle and wavenumber. In this situation, the amplification factor for the amplitude of displacement is over 300.     

Scattering of SH waves from a partially debonded rigid elliptic cylinder

Scattering of SH waves by an embedded rigid elliptic cylinder of finite length, which is partially debonded from elastic soil is studied. The debonding regions are modeled as multiple elliptic arc-shaped interface cracks with non-contacting faces. The scattered wave field is expressed as a Mathieu function expansion with unknown coefficients. The mixed boundary conditions of the problem lead to a set of singular integral equations of the first type in terms of the dislocation density functions of the cracks. A quadrature method is used to solve these integral equations numerically. The results for dynamic stress intensity factors, far-field pattern of the displacement and scattering cross sections are presented. We particularly discuss the effects of the ratio of the short radius to long radius of the cylinder.     

SH波入射时半圆形凸起与凹陷地形的地震动

研究了弹性半空间中半圆形凸起与凹陷相连地形对SH波的散射问题。将整个求解区域分割为2部分,在其中分别构造满足边界条件的位移解,通过移动坐标使之满足“公共边界”以及半圆形凹陷表面上的边界条件,从而建立起求解该问题的无穷代数方程组。最后,给出了地表位移幅值的数值结果以及凸起地形顶点和凹陷地形最低点处位移幅值的反应谱并进行了讨论。     

Scattering and reflection of SH waves around a slope on an elastic wedged space

The scattering and reflection of SH waves by a slope on an elastic wedged space is investigated. A series solution is obtained by using the wave function expansion method. The slope on a wedged space is divided into two subregions by an artificial, auxiliary circular arc. The wave fields with unknown complex coefficients within each sub-region are derived. Applying Graf addition theorem, the scattered waves in the sub-regions are expressed in a global coordinate system. Fourier transform is adopted to derive a consistent form of standing waves in the inner region using the orthogonality of the cosine functions. The boundary-valued problem is solved by stress and displacement continuity along the artificial, auxiliary arc to obtain the unknown complex coefficients. Parametric studies are next performed to investigate how the topography from the slope on the wedged space will affect the scattering and diffraction, and hence the amplification and de-amplification of the SH waves. Numerical results show that the surface motions on the slope of the wedged space is influenced greatly by the topography. Amplification of the surface motions near the slope vertex is significant. The corresponding phases along the wedged space surfaces are consistent with the direction that the SH waves are propagating.     

SH波对圆弧形凸起地形的散射

本文采用“契合”的方法,给出了一个求解平面ＳＨ波对圆弧形凸起形散射的新方法。利用包括圆绵形凸起上边界线在内的一个圆域中预先构造的满足凸起边界应力为零。其余边界应力待定的级数解和其余下的具有圆弧形凹陷的半空间中的解答。通过在其结合面上完成“契合”的过程中分别确定出圆域和半空间听解答。给出了问题的最终结果。利用上述方法,问题的求解仍归结为对一个无穷代数方程组的求解。最后,本文给出了数值结果并对其进行了     

Scattered Surface Wave Energy in the Seismic Coda

One of the many important contributions that Aki has made to seismology pertains to the origin of coda waves (Aki, 1969; Aki and Chouet, 1975). In this paper, I revisit Aki's original idea of the role of scattered surface waves in the seismic coda. Based on the radiative transfer theory, I developed a new set of scattered wave energy equations by including scattered surface waves and body wave to surface wave scattering conversions. The work is an extended study of Zeng et al. (1991), Zeng (1993) and Sato (1994a) on multiple isotropic-scattering, and may shed new insight into the seismic coda wave interpretation. The scattering equations are solved numerically by first discretizing the model at regular grids and then solving the linear integral equations iteratively. The results show that scattered wave energy can be well approximated by body-wave to body wave scattering at earlier arrival times and short distances. At long distances from the source, scattered surface waves dominate scattered body waves at surface stations. Since surface waves are 2-D propagating waves, their scattered energies should in theory follow a common decay curve. The observed common decay trends on seismic coda of local earthquake recordings particular at long lapse times suggest that perhaps later seismic codas are dominated by scattered surface waves. When efficient body wave to surface wave conversion mechanisms are present in the shallow crustal layers, such as soft sediment layers, the scattered surface waves dominate the seismic coda at even early arrival times for shallow sources and at later arrival times for deeper events.     

Dynamic stress concentration and scattering of SH-wave by interface elliptic cylindrical cavity

In this paper, the dynamic stress concentration and scattering of SH-waves by bi-naterial structures that possess an interface elliptic cavity are investigated. First, by using the complex function method, the Green‘s function is constructed. This yields the solution of the displacement field for an elastic half space with a semi-elliptic canyon impacted by an anti-plane harmonic line source loading on the horizontal surface. Then, the problem is divided into an upper and lower half space along the horizontal interface, regarded as a harmony model. In order to satisfythe integral continuity condition,the unknown anti-plane forces are applied to the interface. The integral equations with unknown forces can be established through the continuity condition, and after transformation, the algebraic equations are solved numerically. Finally, the distribution of the dynamic stress concentration factor (DSCF) around the elliptic cavity is given and the effect of different parameters on DSCF is discussed.