Anti-plane （SH） waves diffraction by cavity： analytical an underground semi-circular solution

Diffraction of a two-dimensional (2D) semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called “improved cosine halfrange expansion” algorithm exhibits an excellent performance in reducing displacement residual errors at two rim points of concern. The governing equations are developed in a manner that minimizes the residues of the boundary conditions. Detailed derivation and analysis procedures as well as truncation of infinite linear governing equations are presented. The semi-circular cavity model presented in this paper, due to its simple profile, is expected to be used in seismic wave propagation studies as a benchmark for examining the accuracies of various analytical or numerical methods for mixed-boundary wave propagation problems.     

Out-of-plane (SH) soil-structure interaction: Semi-circular rigid foundation revisited

The model studied in this paper presents an extension of previous work for a shear wall on a semi-circular rigid foundation in an isotropic homogeneous and elastic half-space. The objective is to develop a soil-structure interaction model that can later be applied to the case of a flexible foundation. As shown in the Introduction below, Luco considered the case of a rigid foundation subjected to vertical incident plane SH waves, and Trifunac extended the solution for the same rigid foundation subjected to SH waves but for arbitrary angles of the incidence. In this paper, a new approach and model are presented for the same semi-circular rigid foundation with a tapered-shape (instead of rectangular) superstructure. The analytical expression for the deformation of the semi-circular rigid foundation below this tapered shear wall with soil-structure interaction in an isotropic homogeneous and elastic half-space is thus derived. Results are then compared with those of Trifunac discussed in the section below. This problem formulation can and will later be extended in the case of a flexible foundation that is semi-circular or arbitrarily shaped.     

Soil-structure interaction on shallow rigid circular foundation: plane SH waves from far-field earthquakes

A closed-form wave function analytic solution of two-dimensional scattering and diffraction of incident plane SH-waves by a fl exible wall on a rigid shallow circular foundation embedded in an elastic half-space is presented. This research generalizes the previous solution by Trifunac in 1972, which tackled only the semi-circular foundation, to arbitrary shallow circular-arc foundation cases, and is thus comparatively more realistic. Ground surface displacement spectra at higher frequencies are also obtained. As an analytical series solution, the accuracy and error analysis of the numerical results are also discussed. It was observed from the results that the rise-to-span ratio of the foundation profi le, frequency of incident waves, and mass ratios of different media(foundation-structure-soil) are the three primary factors that may affect the surface ground motion amplitudes near the structure.     

层状半空间中洞室对平面SH波的放大作用

利用间接边界元法,求解了弹性层状半空间中洞室对入射平面SH波的放大作用问题,并以基岩上单一土层中洞室对入射平面SH波的放大作用为例进行了数值计算分析。本文模型的特点之一是考虑了层状场地的动力特性,因而更接近于实际工程;特点之二是计算精度非常高。研究表明,层状半空间中洞室对波的放大作用与均匀半空间中情况有着本质的差别;层状半空间中洞室附近地表动力响应由土层动力特性和洞室对波的散射二者共同决定。土层动力特性不仅影响洞室附近地表位移的幅值,还会影响地表位移的频谱。在土层的前几阶共振频率附近,随着基岩与土层剪切波速比的增大,土层的影响随之增大,而随着土层厚度的增加,土层的影响随之减小,并逐渐趋于均匀半空间情况。     

风化半圆形河谷对柱面SH波的散射解析解

风化河谷地震效应对建在此类场地上工程结构(如大坝、桥梁)的动力稳定性具有潜在的威胁,为揭示河谷风化层对地面运动的影响,利用波函数展开法推导风化半圆形河谷对线源柱面SH波散射问题的解析解。计算不同震源位置条件下风化半圆形河谷的地表位移反应和地面运动放大因子,分析地面运动的幅值和形态,发现随着震源距离的增大,在河谷附近地震动的放大和衰减的交替更频繁。     

Diffraction of plane SV waves by a cavity in poroelastic half-space

This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green’s functions of compressive and shear wave sources are derived based on Biot’s theory. The scattered waves are constructed using fi ctitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verifi cation of the accuracy is performed by: (1) checking the satisfaction extent of the boundary c...     

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

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

Scattering and diffraction of plane P-waves in a 2-D elastic half-space Ⅱ:shallow arbitrary shaped canyon

Scattering and Diffraction of elastic in-plane P- and SV- waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong-motion seismologists for over forty years. The case of out-of-plane SH waves on the same elastic canyon that is semi-circular in shape on the half-space surface is the first such problem that was solved by analytic closed form solutions over forty years ago by Trifunac. The corresponding case of in-plane P- and SV-waves on the same circular canyon is a much more complicated problem because, the in-plane P- and SV- scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by the author in the work of Lee and Liu. This paper uses the technique of Lee and Liu of defining these stress-free scattered waves to solve the problem of the scattered and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape.     

Diffraction of obliquely incident waves by a cylindrical cavity embedded in a layered viscoelastic half-space

The three-dimensional harmonic response in the vicinity of an infinitely long, cylindrical cavity of circular cross-section buried in a layered, viscoelastic half-space is obtained when the half-space is subjected to homogeneous plane waves and surface waves impinging at an oblique angle with respect to the axis of the cavity. The solution is obtained by an indirect boundary integral method based on the use of moving Green's functions for the viscoelastic half-space. Numerical results describing the motion on the ground surface and the motion and stresses on the wall of the cavity are presented for obliquely incident P-, SV-, SH- and Rayleigh waves with different horizontal angles of incidence.     

饱和土半空间中地下圆形衬砌洞室对平面SV波的散射

在Biot饱和多孔介质动力学理论的基础上,利用Fourier-Bessel级数展开法,通过对舍有衬砌洞室的局部场地进行波场分析,得到饱和土半空间中圆柱形衬砌洞室对平面SV波的散射问题的解析解。经验证,本文得到的解可以退化为半空间单相介质的情况。通过与已有的单相弹性介质半空间中圆柱形衬砌洞室对平面SV波散射问题的解析解的对比,验证了此解的正确性。在解析解的基础上,数值计算给出洞口动应力集中放大系数,分析了入射频率和孔洞埋深对柱面上的应力集中因子的影响。     

Two-dimensional scattering and diffraction of P- and SV-waves around a semi-circular canyon in an elastic half-space: An analytic solution via a stress-free wave function

A well-defined boundary-valued problem of wave scattering and diffraction in elastic half-space should have closed-form analytic solutions. This two-dimensional (2-D) scattering around a semi-circular canyon in elastic half-space subjected to seismic plane and cylindrical waves has long been a challenging boundary-value problem. In all cases, the diffracted waves will consist of both longitudinal (P-) and shear (S-) rotational waves. Together at the half-space surface, these in-plane longitudinal P- and shear SV-waves are not orthogonal over the infinite half-space flat-plane boundary. Thus, to simultaneously satisfy both the zero normal and shear stresses at the flat-plane boundary, some approximation of the geometry and/or wave functions often has to be made, or in some cases, relaxed (disregarded). This paper re-examines this two-dimensional (2-D) boundary-value problem from an applied mathematics points of view and redefines the proper form of the orthogonal cylindrical-wave functions for both the longitudinal P- and shear SV-waves so that they can together simultaneously satisfy the zero-stress boundary conditions at the half-space surface. With the zero-stress boundary conditions satisfied at the half-space surface, the most difficult part of the problem will be solved, and the remaining boundary conditions at the finite-canyon surface are then comparatively less complicated to solve. This is now a closed-form analytic solution of the 2-D boundary-valued problem satisfying the half-space zero-stress boundary conditions exactly.     

弹性层状半空间中无限长洞室对斜入射 平面SH波的三维散射(Ⅱ) 数值结果与分析

以基岩上单一土层场地为例, 计算分析了在斜入射平面SH波作用下弹性层状半空间中无限长洞室附近的地表位移. 研究表明, 层状半空间中地下洞室对波的散射与均匀半空间情况存在显著差别. 层状场地由于考虑了场地自身的动力特性, 使得洞室附近地表位移幅值的空间变化更为复杂, 基岩与土层刚度比、 土层厚度对散射效应均有着重要影响. 随着基岩与土层刚度比的增大, 地表位移幅值整体上逐渐增大; 随着土层厚度的增大, 土层对地表位移幅值的影响逐渐减小. 在频域解答的基础上, 给出了层状半空间中洞室对斜入射SH波散射的时域解答, 并以Ricker波为例进行了数值计算.      

Analytic wave series solution of out-of-plane(SH) waves diffraction by an almost semi-circular shallow cylindrical hill

A closed-form wave equation analytic solution of two-dimensional scattering and diffraction of outof-plane(SH) waves by an almost semi-circular shallow cylindrical hill on a flat, elastic and homogeneous half space is proposed by applying the discrete Fourier series expansions of sine and cosine functions. The semi-circular hill problem is discussed as a special case for the new formulated equation.Compared with the previous semi-circular cases solutions, the present method can give surface displacement amplitudes which agrees well with previous results. Although the proposed equation can only solve the problem of SH-waves diffracted by almost semi-circular shallow hills, the stress and displacement residual amplitudes are numerical insignificantly everywhere. Moreover, the influences of the depth-towidth ratio(a parameter defined in this paper to evaluate the shallowness of the topography of hills) on ground motions are presented and summarized. The limitations and errors of truncation from Graf's addition theorem and Fourier series equations in the present paper are also discussed.     

平面P波在饱和半空间中洞室周围的散射(I):解析解

利用波函数展开法给出了入射平面P波在饱和半空间中圆形洞室周围散射问题的一个解析解。半空间假定为无粘性流体饱和介质,满足Biot理论。采用一种基于实验数据的孔隙率和模量之间的线性关系来确定Biot模型中的介质参数。解答考虑了透水边界和非透水边界两种情况。对边界条件进行了数值检验,结果表明,随着级数截断项数的增大,边界残量衰减很快。解答为进一步研究入射波频率和角度、边界渗透条件、孔隙率、泊松比等参数对散射的影响奠定了基础。     

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

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

Estimation of Elastic Moduli in a Compressible Gibson Half-space by Inverting Rayleigh-wave Phase Velocity

A Gibson half-space model (a non-layered Earth model) has the shear modulus varying linearly with depth in an inhomogeneous elastic half-space. In a half-space of sedimentary granular soil under a geostatic state of initial stress, the density and the Poisson’s ratio do not vary considerably with depth. In such an Earth body, the dynamic shear modulus is the parameter that mainly affects the dispersion of propagating waves. We have estimated shear-wave velocities in the compressible Gibson half-space by inverting Rayleigh-wave phase velocities. An analytical dispersion law of Rayleigh-type waves in a compressible Gibson half-space is given in an algebraic form, which makes our inversion process extremely simple and fast. The convergence of the weighted damping solution is guaranteed through selection of the damping factor using the Levenberg-Marquardt method. Calculation efficiency is achieved by reconstructing a weighted damping solution using singular value decomposition techniques. The main advantage of this algorithm is that only three parameters define the compressible Gibson half-space model. Theoretically, to determine the model by the inversion, only three Rayleigh-wave phase velocities at different frequencies are required. This is useful in practice where Rayleigh-wave energy is only developed in a limited frequency range or at certain frequencies as data acquired at manmade structures such as dams and levees. Two real examples are presented and verified by borehole S-wave velocity measurements. The results of these real examples are also compared with the results of the layered-Earth model.     

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.     

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

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

Scattering of anti-plane (SH) waves by a semi-elliptical hill: II—deep hill

The ground motion and wave propagation resulting from the scattering of anti-plane (SH) waves by a deep semi-elliptical hill on an elastic half-space are investigated. The solution is the closed-form one based on the method of wave-function expansion, and the expansion of elliptical sine and cosine over the half-range [−π/2, π/2]. The analysis shows that the interference of ground motion depends on several factors such as frequency of the incident waves, incident angle, and the aspect ratio of elliptical hills. Similar to part I: Shallow Hill, numerical results when semi-elliptical tends to semi-circular are compared with the work done by Lee et al. [17].     

In-plane soil–structure interaction in layered,fluid-saturated,poroelastic half-space I: Structural response

Linear in-plane soil–structure interaction in two dimensions (2D) is studied in fluid-saturated, poroelastic, layered half-space using the Indirect Boundary Element Method (IBEM). The structure is a shear wall supported by a rigid embedded foundation. Exact stiffness matrices for the soil layer and half-space, and Green׳s functions of uniformly distributed loads and pore pressure on an inclined line are derived. Results of the system response in the frequency domain are presented for the special case of single soil layer over bedrock, semi-circular foundation and zero seepage force. The effects of water saturation, soil porosity, depth of soil layer, rigidity contrast between layer and bedrock are investigated in the frequency domain for incident plane P- and SV waves. The results suggest that water saturation may cause increase of the system frequency by more than 10%.