Scattering of plane sh waves by a semi-elliptical canyon

In this paper we analyse the two-dimensional scattering and diffraction of plane SH waves by a semi-elliptical canyon. The exact series solution of the problem, for general angle of incidence of the plane SH waves, has been used to examine the dependence of surface amplifications inside and near the canyon. The nature of ground motion has been found to depend on two key parameters: (a) The angle of incidence. (b) The ratio of the canyon width to the wave length of incident SH waves. For short incident waves surface displacement amplitudes change rapidly from one point to another, while for the long waves and shallow canyons displacement amplitudes display only minor departure from the uniform half-space amplification of 2. For shallow canyons and long incident waves, the angle of incidence introduces only minor changes into the overall behaviour of surface amplitudes. For deep canyons and nearly grazing incidences, a prominent shadow zone is realized behind the canyon.     

Application of the weighted residual method to diffraction by 2-D canyons of arbitrary shape: I. Incident SH waves

A solution for the two-dimensional scattering and diffraction of plane SH waves by canyons of arbitrary shape in an elastic half space is presented. The wave field for arbitrary geometry in this paper is computed numerically by the method of weighted residues (moment method). The wave displacement field computed by the present residual method for the case of a semi-circular canyon was shown to agree analytically and numerically with that computed by the exact closed form series solution. The same observations about ground amplifications, their dependence on frequencies and orientations of the incident waves, can be stated here for canyons of arbitrary shape as previously made for circular canyons.     

Scattering of plane harmonic SH,SV, P and rayleigh waves by non-axisymmetric three-dimensional canyons: A wave function expansion approach

Scattering of elastic waves by three-dimensional canyons embedded within an elastic half-space is investigated by using a wave function expansion technique. The geometry of the canyon is assumed to be non-axisymmetric. The canyon is subjected to incident plane Rayleigh waves and oblique incident SH, SV and P waves. The unknown scattered wavefield is expressed in terms of spherical wave functions which satisfy the equations of motion and radiation conditions at infinity, but they do not satisfy stress-free boundary conditions at the half-space surface. The boundary conditions are imposed locally in the least-squares sense at several points on the surface of the canyon and the half-space. Through a comparative study the validity and limitations of two-dimensional approximations (antiplane strain and plane strain models) have been examined. It is shown that scattering of waves by three-dimensional canyons may cause substantial change in the surface displacement patterns in comparison to the two-dimensional models. These results emphasize the need for three-dimensional modelling of realistic problems of interest in strong ground motion seismology and earthquake engineering.     

Three-dimensional analysis of spatially varying ground motions around a uniform canyon in a homogeneous half-space

Department of Civil Engineering, University of California, Berkeley, CA 94720, U.S.A. A direct boundary element method to determine the three-dimensional seismic response of an infinitely-long canyon of arbitrary but uniform cross-section cut in a homogeneous viscoelastic half-space is presented. The seismic excitation is represented by P, SV, SH or Rayleigh waves at arbitrary angles with respect to the axis of the canyon. The accuracy of the procedure and implementing computer program is demonstrated by comparison with previous solutions for the limiting case of two-dimensional response, recently obtained three-dimensional response results for infinitely-long canyons, and three-dimensional boundary method solutions presented in this paper for finite canyons.     

The use of 3D modeling for the prediction of the seismic demands on the gravity dams

Seismic behavior of gravity dams has long been evaluated using a representative two‐dimensional (2D) system. Formulated for the gravity dams built in wide canyons, the assumption is nevertheless utilized extensively for almost all concrete dams due to the established procedures as well as the expected computational costs of a three‐dimensional model. However, a significant number of roller‐compacted concrete dams, characterized as such systems, do not conform to the basic assumptions of these methods by violating the conditions on canyon dimensions and joint‐spacing/details. Based on the premise that the 2D modeling assumption is overstretched for practical purposes in a variety of settings, the purpose of this study is to critically evaluate the use of 2D modeling for the prediction of the seismic demands on these systems. Using a rigorous soil–structure interaction approach, the difference between the two and three‐dimensional response for gravity dams was investigated first in the frequency domain for a range of canyon widths and foundation to dam moduli ratios. Then, the time domain differences between the crest displacements and the maximum principal stress were obtained using 70 different ground motions in order to show the possible bias introduced into the analysis results due to the modeling approach. The results of the study show that even for relatively wide canyons, the 2D analysis can lead to misleading predictions. Copyright © 2017 John Wiley & Sons, Ltd.     

Diffraction of plane sh waves in a half-space

Scattering of antiplane shear waves (SH) in two dimensions by surface and near-surface defects in a homogeneous, isotropic elastic semi-infinite medium has been studied. Attention has been focused here in the range of medium to long wavelengths. A combined finite element and analytical technique has been used to study the problems of scattering by semi-circular and triangular canyons. The results for the former case are compared with the known exact solution and those for the latter case are compared with some available approximate solutions. Finally a problem of multiple scattering by a triangular canyon and a nearby circular tunnel is studied. Numerical results are presented showing the effects of multiple scattering and different angles of incidence. These results are of interest in earthquake engineering.     

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

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

汶川地震震损水库土坝动力反应与几何坝形的经验关系

在5·12汶川地震后的四川省水库土坝震害调查成果的基础上,选取汶川地震中受损的有完整资料的96座水库土坝为研究对象,应用等效线性模型对土坝进行了二维动力反应分析。选择3条有代表性的汶川地震实测记录,以三水准峰值加速度输入,得到土坝的动力反应（放大系数、最大动剪应力）与土坝几何形状（宽高比、上游坡比、坝高）间的经验关系。结果表明,土坝放大系数和最大动剪应力随着宽高比和上游坡比的增大而减小,随着坝高的增大而增大;输入波的频谱和峰值强度均对土坝动力反应与其几何形状的经验关系有重要影响。     

Ground motion at canyons of arbitrary shape under incident sh waves

A method for calculating the two-dimensional scattering of incident SH waves by canyons of arbitrary shape is presented. The problem is formulated in terms of a Fredholm integral equation of the first kind with the integration path outside the boundary. Point-source discretization and a least-squares scheme are used. Numerical results are compared with the known analytic solution for a semi-cylindrical canyon. Spatial variations of surface amplitudes are computed for triangular and half-cycle sinusoidal canyons as well.     

含峭壁V形峡谷对地震SH波散射的解析解

地表地形常引起地震动的局部放大,这是由于地震波传播至局部地形时产生了散射现象.本文利用波函数展开方法和区域匹配技术,提出了含峭壁V形峡谷对平面SH波散射问题的解析解,并进行了退化验证.通过频域内的参数分析,揭示了峭壁深度、入射波频率和角度等因素对峡谷场地地面运动的影响规律,发现上部峭壁会增强峡谷对地震动的地形放大效应....     

RESPONSE OF EARTH DAMS TO P AND SV WAVES USING COUPLED FE-BE FORMULATION

A study of the effects of dam–foundation interaction on the response of earth dams to obliquely incident P and SV waves is presented. Emphasis is placed on the effects of the foundation flexibility and the spatial variability of the ground motion. The study is based on a rigorous hybrid numerical formulation that combines the efficiency and versatility of the Finite Element Method (FEM) and the ability of Boundary Element Method (BEM) to account for the radiation conditions at the far field. The developed hybrid method is very powerful and can be used efficiently to obtain accurate solutions of problems of complex geometry, material heterogeneity and, for time-domain analyses, material nonlinearity. The 2-D frequency-domain formulation is used here to investigate the response of infinitely long earth dams to obliquely incident P and SV waves. By accounting rigorously for the energy radiated back into the half-space, the study demonstrates the dramatic effect of the flexibility of the foundation rock in reducing the overall response of the dam. The effects of the spatial variability of the ground motion for P and SV waves travelling across the width of the dam are also important, but somewhat less pronounced than those of the foundation flexibility.     

圆弧形凹陷饱和土场地对平面P波散射问题的解析解

以Biot 饱和多孔介质动力学理论为基础，利用Fourier Bessel级数展开法, 得到饱和多孔介质半空间中圆弧形凹陷地形对平面P波的散射问题的解析解. 数值计算给出地表位移幅值，分析了入射波波长、入射角、圆弧高宽比对地表位移幅值的影响，并与现存的在单相介质情况下得到的结论进行对比.     

Effect of reservoir bottom on earthquake response of concrete dams

The absorption of hydrodynamic pressure waves at the reservoir bottom has dominant effects on the structural response of the dam when subjected to ground motion. In the present study, a model is proposed for the absorption effects of the reservoir bottom in the earthquake analysis of dams. The model utilizes the wave reflection coefficient approach and is based on the solution of the wave equation in a sediment layer of viscoelastic material with a constant thickness overlying an elastic, semi-infinite foundation. Numerical studies were conducted to evaluate the effect of the sediment layer thickness and material properties as well as the effect of reflection of waves from the underlying rock. It is shown that the current approach of assuming the wave reflection coefficient at the reservoir bottom based on the characteristics of the sediment material and excluding the effect of the reflected waves from the underlying rock, may significantly underestimate the seismic response of the dam.     

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.     

圆弧形凹陷在SH波入射下瞬态反应的解析解

针对现有地震波散射问题的解析解均在频域内给出,其适用的频带较窄,不能获取到可靠的瞬态反应解析解的问题,本文利用宽频带频域解析解进一步获取了不同空间点地震反应的宽频带频响函数,并基于傅里叶变换得出SH波入射下圆弧形凹陷地形瞬态动力响应的解析解.以上述为基础,研究了输入脉冲与实际地震动时凹陷的瞬态反应的基本特征,并重点分析...     

SH Wave Number Green’s Function for a Layered,Elastic Half-Space. Part I: Theory and Dynamic Canyon Response by the Discrete Wave Number Boundary Element Method

We present a closed-form frequency-wave number (ω – k) Green’s function for a layered, elastic half-space under SH wave propagation. It is shown that for every (ω – k) pair, the fundamental solution exhibits two distinctive features: (1) the original layered system can be reduced to a system composed by the uppermost superficial layer over an equivalent half-space; (2) the fundamental solution can be partitioned into three different fundamental solutions, each one carrying out a different physical interpretation, i.e., an equivalent half-space, source image impact, and dispersive wave effect, respectively. Such an interpretation allows the proper use of analytical and numerical integration schemes, and ensures the correct assessment of Cauchy principal value integrals. Our method is based upon a stiffness-matrix scheme, and as a first approach we assume that observation points and the impulsive SH line-source are spatially located within the uppermost superficial layer. We use a discrete wave number boundary element strategy to test the benefits of our fundamental solution. We benchmark our results against reported solutions for an infinitely long circular canyon subjected to oblique incident SH waves within a homogeneous half-space. Our results show an almost exact agreement with previous studies. We further shed light on the impact of horizontal strata by examining the dynamic response of the circular canyon to oblique incident SH waves under different layered half-space configurations and incident angles. Our results show that modifications in the layering structure manifest by larger peak ground responses, and stronger spatial variability due to interactions of the canyon geometry with trapped Love waves in combination with impedance contrast effects.     

Effects of reservoir bottom absorption on earthquake response of concrete gravity dams

A procedure is presented to analyse the response of concrete gravity dams due to horizontal and vertical earthquake ground motion components considering dam-water interaction and partial absorption of hydrodynamic pressure waves at the reservoir bottom into the foundation medium. The effects of reservoir bottom absorption on the hydrodynamic force on a rigid dam are examined first. The harmonic response of an idealized dam cross-section is presented for a wide range of parameters characterizing the properties of the dam, the impounded water and the foundation medium. Based on these frequency response functions the effects of dam-water interaction and of reservoir bottom absorption in the response of dams due to horizontal and vertical components of ground motion are investigated.     

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

Soil–structure interaction for a SDOF oscillator supported by a flexible foundation embedded in a half-space: Closed-form solution for incident plane SH-waves

A closed-form analytical solution is presented for the dynamic response of a SDOF oscillator, supported by a flexible foundation embedded in an elastic half-space, and excited by plane SH waves. The solution is obtained by the wave function expansion method. The solution is verified for the special case of a rigid foundation by comparison with published results. The model is used to investigate the effect of the foundation flexibility on the system response. The results show that the effect is significant for both foundation response and structural relative response. For a system with more flexible foundation, the radiation damping is smaller, the foundation response is larger, especially for obliquely incident waves, while the structural relative response is smaller, and the system frequency shifts towards lower frequencies. This simple model may be helpful to obtain insight into the effects of soil–structure interaction for a slim structure on an extended flexible foundation.     

Elasto-plastic earthquake shear-response of one-dimensional earth dam models

A simplified analysis procedure for the non-linear hysteretic earthquake-response of earth dams is presented. The dam is modelled as a one-dimensional hysteretic shear-wedge subjected to base excitation. The hysteretic stress-strain behaviour of the dam materials is modelled by using elasto-plastic constitutive equations based on multi-surface kinematic plasticity theory. The method is based on a Galerkin formulation of the equations of motion in which the solution is expanded using eigenmodes of the linearized problem defined over the spatial domain occupied by the dam. The technique is applied to analyse the non-linear dynamic response of an earth dam subject to two very different input ground motions. The following investigations are presented: (i) comparison between the results obtained using two soil models depicting different nonlinear properties, (ii) comparison between the results of the one-mode and the multi-mode solution expansions, (iii) comparison with the results obtained through an elaborate finite element representation of the dam, and finally, (iv) comparison with the results obtained through the Makdisi-Seed¹¹ iterative procedure for earth dam analysis. The comparisons show that the proposed technique can be used to determine adequately the transient earthquake response of long earth dams. Furthermore, the efficiency and low computational cost make the technique very attractive; it can easily and systematically be extended to two- and three-dimensional calculations of earth dam response.