Amplification of elastic waves by a three dimensional valley. Part 1: Steady state response

Steady state scattering of incident P, SV, SH and Rayleigh waves by general non-axisymmetric three dimensional dipping layers is investigated by using an indirect boundary integral equation method. Material of the half-space and the layer is assumed to be linear, weakly anelastic, homogeneous and isotropic. Systematic comparisons between three dimensional and two dimensional models demonstrate that the validity of a two dimensional approximation for a given basin shape may be affected strongly by changes in azimuthal angle of incidence, type of incident wave and frequency. The discrepancies of two dimensional modelling appear to be much more pronounced for the case of an incident SH wave. Another important feature of the results is the existence of strong coupling between P/SV and SH modes, which has no correspondence in two dimensional models. Such off-azimuthal mode conversions are particularly strong for an incident SH wave.     

Scattering of obliquely incident seismic waves by a cylindrical valley in a layered half-space

Three-dimensional scattering of seismic waves by a cylindrical alluvial valley embedded in a layered half-space is investigated by using the combination of the boundary integral representation and the finite element method. The surface displacements due to incident plane harmonic body waves (P, SV and SH) propagating at an arbitrary angle to the axis of the cylindrical valley are evaluated numerically for two semi-elliptical alluvial valleys. The presence of the layer is found to have a strong effect on the amplification of the surface displacements in some cases. The three-dimensional motion seems to be quite critical and may cause large amplification. The surface ground motion becomes significant when compared with corresponding free-field motion as the wavelengths become comparable to the characteristic length of the valley. The maximum amplification always occurs atop the valley. Numerical results show that the amplitude and the amplification pattern of the surface displacement strongly depend upon the frequency, the angle and the type of the incident waves.     

Seismic stability of jointed rock slopes under obliquely incident earthquake waves

Seismic stability of slopes has been traditionally analyzed with vertically propagated earthquake waves. However, for rock slopes, the earthquake waves might approach the outcrop still with a evidently oblique direction. To investigate the impact of obliquely incident earthquake excitations, the input method for SV and P waves with arbitrary incident angles is conducted, respectively, by adopting the equivalent nodal force method together with a viscous-spring boundary. Then, the input method is introduced within the framework of ABAQUS software and verified by a numerical example. Both SV and P waves input are considered herein for a 2D jointed rock slope. For the jointed rock mass, the jointed material model in ABAQUS software is employed to simulate its behavior as a continuum. Results of the study show that the earthquake incident angles have significance on the seismic stability of jointed rock slopes. The larger the incident angle, the greater the risk of slope instability. Furthermore, the stability of the jointed rock slopes also is affected by wave types of earthquakes heavily. P waves induce weaker responses and SV waves are shown to be more critical.     

Influences of type of wave and angle of incidence on seismic bending moments in pile foundations

When analysing the seismic response of pile groups, a vertically‐incident wavefield is usually employed even though it does not necessarily correspond to the worst case scenario. This work aims to study the influences of both the type of seismic body wave and its angle of incidence on the dynamic response of pile foundations. To this end, the formulation of SV, SH and P obliquely‐incident waves is presented and implemented in a frequency‐domain boundary element‐finite element code for the dynamic analysis of pile foundations and piled structures. Results are presented in terms of bending moments at cap level of single piles and 3 × 3 pile groups, both in frequency and in time domains. It is found that, in general, the vertical incidence is not the most unfavourable situation. In particular, obliquely‐incident SV waves with angles of incidence smaller than the critical one, a situation in which the mechanism of propagation of the waves in the soil changes and surface waves appear, yield bending moments much larger than those obtained for vertically‐incident wavefields. It is also shown that the influence of pile‐to‐pile interaction on the kinematic bending moments becomes significant for non‐vertical incidence, especially for P and SV waves. Copyright © 2013 John Wiley & Sons, Ltd.     

Effects of a covering layer in a circular-arc canyon on incident plane SV waves

Introduction The effect of local topographies on seismic wave propagation is one of the most attractive topics in seismology. It may be resolved by either an analytical method or a numerical method. Here, the analytical method is the wave function expansion method; numerical methods include the finite difference, finite element, boundary integration equation, and discrete wave number method, etc. The advantage of these numerical methods is that they can be applied to local inhomogeneity and ir…     

圆弧形凹陷地形表面覆盖层对入射平面SV波的影响

利用Fourier Bessel级数展开法，给出了表面具有覆盖层的圆弧形凹陷地形对入射平面SV波散射问题的一个解析解，并利用该解分析了不同深宽比凹陷地形表面覆盖层刚度和厚度对入射SV波的影响．结果表明，凹陷地形表面覆盖层的存在，即使厚度很薄，也会显著增加对入射SV波的放大作用，该放大作用可达到单一凹陷地形的2.5倍以上．覆盖层刚度和厚度的变化对入射平面SV波也具有很大影响．      

Structural response for stochastic kinematic interaction

The influence of stochastic kinematic interaction (SKI) on structural response is investigated in this paper. The SKI is evaluated through a computational model based on the boundary element method (BEM) formulated in the frequency domain. The singular integrals required in the computation of BEM are evaluated in a closed form. It is assumed that the foundation input motion (FIM) is the result of the superposition of many plane, stationary, correlated stochastic SH‐, P‐ and SV‐waves travelling within a homogeneous viscoelastic soil at different angles. The results obtained indicate that the effect of SKI on the foundation response is qualitatively similar to that of wave passage. Both effects involve a reduction of translational components of the response at intermediate and high frequencies and creation of a rotational response component at intermediate frequencies, which decreases at high frequencies. While, it is found that the SKI decreases the maximum response of structures built on embedded rigid strip foundations excited by SH‐ and P‐waves, it increases the maximum response for SV‐waves, except when the natural frequency of the structure is less than 0.5 Hz and for short structures excited by shallowly incident SV‐waves. Copyright © 2001 John Wiley & Sons, Ltd.     

Influence of scattering of SH-waves in dynamic interaction of shear wall with soil layers

The dynamic soil-structure interaction of a shear wall embedded in elastic isotropic and homogeneous soil layers underlain by bedrock, subjected to SH waves, is modeled in the present article. The soil layers consist of irregular interfaces and it has been shown that the scattering due to the roughness of the layers has significant effect on the displacement of both the foundation and the shear wall. To demonstrate the phenomena indirect boundary element method(IBEM) has been used on the basis of its validation in previous problems of similar type. The system response is compared with the analytical solution of the same type of model for vertically propagating incident SH waves. It is observed that for the low frequency of wave, displacement is abruptly high, and as a result the combination of shear wall and foundation perceives resonance. The thickness of the soil layer, mass of the shear wall, stiffness of the bedrock and the soil layers all affects the system frequency and displacement.     

Ground motion amplification due to elastic inclusions in a half-space

Scattering of plane harmonic SH, P, SV and Rayleigh waves by several inclusions of arbitrary shape, completely embedded into an elastic half-space, is considered. Perfect bonding between the half-space and the inclusions is assumed. The problem is investigated for linear, isotropic and homogeneous elastic materials. The displacement field is evaluated throughout the elastic medium so that the continuity conditions between the half-space and the inclusions are satisfied in mean-square sense. Numerical results of the surface displacement field are evaluated for single and two elliptic inclusions. The results show the following: (a) presence of a subsurface inhomogeneity may lead to large amplifications of the surface ground motion; (2) different surface displacement patterns emerge for different incident waves; (3) the presence of an additional inclusion may change significantly the surface displacement response of a single inclusion; (4) the surface motion extremes strongly depend upon (i) angle of incidence; (ii) frequency of incident field; (iii) embedment depth of the inclusions; (iv) separation distance between the inclusions; (v) material properties of the half-space and the inclusions; and (vi) location of observation point on the surface of the half-space.     

弹性层状半空间中凸起地形对入射平面SH波的放大作用

对Wolf理论进行拓展,使之可解决凸起地形对波的散射问题,进而利用间接边界元法,求解了弹性层状半空间中凸起地形对入射平面SH波的放大作用问题。本文模型的显著特点之一是考虑了层状半空间的动力特性以及层状半空间和凸起地形的阻尼;特点之二是计算精度高。文中以基岩上单一土层中半圆凸起地形对入射平面SH波的放大作用为例进行了数值计算分析。研究表明,基岩上单一土层中凸起地形对入射平面SH波放大作用和均匀半空间中凸起地形有着本质的差别;土层动力特性不仅影响凸起地形地表位移的幅值,还会影响地表位移的频谱;阻尼会显著降低凸起地形对高频波的放大作用。     

A note on near-field site amplification effects of ground motion from a radially inhomogeneous valley

To improve the understanding of the near-field soil and topographic amplification effects, an analytical solution by the authors for the scattering of plane SH waves by a radially inhomogeneous semi-cylindrical valley is extended to the case of a line source of cylindrical SH waves. Upon confirmation of its accuracy with past exact solutions for a homogeneous and an inhomogeneous semi-cylindrical valley under far-field plane SH waves, the extended solution is used to calculate the ground motion amplification factors for both the homogeneous and inhomogeneous valleys subjected to near-field waves. A comprehensive parametric study is conducted with respect to the location of the wave source, the dimensionless frequency of the incident waves, and the inhomogeneity degree of the covering soil layer. It is found that more amplifications and reductions of ground motions will occur within a certain range in and around the valley as the sources are located further. Consistent with the far-field case, it is confirmed that an increase of the degree of inhomogeneity of the covering soil layer generally amplifies the ground motions significantly.     

地震波斜入射下层状场地中地下综合管廊地震响应分析

利用黏弹性人工边界和等效地震荷载时域波动输入方法,结合土层和半空间的精确动力刚度矩阵,实现了地震波斜入射下层状场地地下综合管廊地震反应分析,建立了不同场地条件下地下综合管廊分析模型。计算结果表明:地震波倾斜入射情况下,综合管廊结构地震响应与垂直入射时具有显著差异,一般SV波以30°临界角附近入射时结构地震反应最为剧烈;地下综合管廊动应力集中主要分布在管廊角部、中柱上下端;成层土波速结构变化对地下综合管廊地震反应亦具有显著影响。总体上看:当穿越软夹层时管廊结构地震反应更为剧烈,且覆盖层越厚,管廊结构内力幅值越大。因此地下综合管廊结构抗震设计宜考虑地震波倾斜入射及场地土层性质的影响。     

Amplification of elastic waves by a three dimensional valley. Part 2: Transient response

Transient response of three dimensional dipping layers of different shapes subjected to incident P, SV, SH and Rayleigh waves is investigated. The time domain response is constructed from steady state solutions through the Fourier synthesis. An indirect boundary integral equation method is applied to calculate the required steady state solutions. The material of the half-space and the layer is assumed to be linear, weakly inelastic, homogeneous and isotropic. Numerical results show that the maximum amplification of motion is strongly dependent upon the type of incident wave, the shape of the basin and signal frequency. The change in the shape of the valley from hemispherical to semi-prolate causes a significant increase in the amplitude of surface waves near the edges; however, the maximum amplification of motion near the centre of the valley decreases. This phenomenon is especially apparent for the case of an incident P wave. In comparison to the corresponding two dimensional responses, the amplitude of motion near the centre of the valley is in general higher for three dimensional models.     

Response of a circular foundation on a uniform half-space to elastic waves

An alternative technique to obtain the dynamic response of a massless rigid circular foundation resting on a uniform elastic half-space when subjected to harmonic plane waves is presented. The technique relies on the use of an integral representation involving the free-field ground motion and the contact tractions obtained in the course of calculating the dynamic force–displacement relationship of the foundation for external forces. Tables listing the translational and rotational components of the response of the foundation for non-vertically incident SH, P, SV and Rayleigh waves are presented.     

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%.     

有限弹性变形地壳中的地震波

地壳由半无限大的基岩上一层厚度为H^-的表土层组成，入射地震波为垂直的SH波，产生水平地面运动。当浅源大地震发生时，在极震区以外行波传播产生地面运动将使地壳介质有非线性的有限弹性变形。用小参数摄动法使非线性控制方程为线性化的小参数各阶控制方程，得出头两阶线性控制方程的解析解。     

A plane strain soil-structure interaction model

A plane strain model for dynamic soil-structure interaction problems under harmonic state is presented. The boundary element method is used to study the response of a homogeneous isotropic linear elastic soil. The far field displacement at the free surface is approximated by an outgoing Rayleigh wave. The finite element method is used to describe the response of the building, of the foundation and possibly of a finite part of the inhomogeneous non-linear soil. Two coupling procedures are described. The model is applied to a problem previously studied in the antiplane case. Incident P, SV and Rayleigh waves are considered. The results show an amplification and an attenuation of the structure motion with frequency when incident Rayleigh waves and P, SV body waves are respectively considered.     

含黏滞流体各向异性孔隙介质中弹性波的频散和衰减

基于Biot理论,考虑液相的黏弹性变形和固液相接触面上的相对扭转,提出了含黏滞流体VTI孔隙介质模型.从理论上推导出,在该模型中除存在快P波、慢P波、SV波、SH波以外,还将存在两种新横波－慢SV波和慢SH波.数值模拟分析了6种弹性波的相速度、衰减、液固相振幅比随孔隙度、频率的变化规律以及快P波、快SV波的衰减随流体性质、渗透率、入射角的变化规律.结果表明慢SV波和慢SH波主要在液相中传播,高频高孔隙度时,速度较高;大角度入射时,快P波衰减表现出明显的各向异性,而快SV波的衰减则基本不变;储层纵向和横向渗透率存在差异时,快SV波衰减大的方向渗透率高.     

弹性波逆散射微扰论

本文将普遍声逆散射微扰论应用于弹性波层析成像问题,在Born变换下推出了以旋转角为补偿参数的各阶微扰重建公式,实现了对非均匀各向同性散射体内3个参数（质量密度ρ和两个Lamé系数λ,μ）的同时重建. 对于层析成像问题,在弹性波的传播过程中P波与SV波有耦合,但它们不会和SH波发生耦合,于是可以得到3个形式相对简单的标量方程. 在Born变换下,在散射波中引入微扰参数,将散射体的3个参数分别按该微扰参数展开,然后利用二维自由空间的Green函数分别得到散射的P波、SV波和SH波的积分表示. 最后,经一维傅氏变换后,得到Born变换下散射体3个参数的各阶微扰重建公式.     

Scattering of elastic waves by a general anisotropic basin. Part 1: a 2D model

Scattering of incident plane harmonic pseudo P‐, SH‐, and SV‐waves by a two‐dimensional basin of arbitrary shape is investigated by using an indirect boundary integral equation approach. The basin and surrounding half‐space are assumed to be generally anisotropic, homogeneous, linearly elastic solids. No material symmetries are assumed. The unknown scattered waves are expressed as linear combinations of full‐space time‐harmonic two‐dimensional Green functions. Using the Radon transform, the Green functions are obtained in the form of finite integrals over a unit circle. An algorithm for the accurate and efficient numerical evaluation of the Green functions is discussed. A detailed convergence and parametric analysis of the problem is presented. Excellent agreement is obtained with isotropic results available in the literature. Steady‐state surface ground motion is presented for semi‐circular basins with generally anisotropic material properties. The results show that surface motion strongly depends upon the material properties of the basin as well as the angle of incidence and frequency of the incident wave. Significant mode conversion can be observed for general triclinic materials which are not present in isotropic models. Comparison with an isotropic basin response demonstrates that anisotropy is very important for assessing the nature of surface motion atop basins. Copyright © 2001 John Wiley & Sons, Ltd.