多层弹性半空间中的地震波（一）

为了了解地震震源和地球介质的性质,很有必要对地震波的辐射、传播和衰减问题作仔细的分析。作为一种近似,可以暂且忽略地球的曲率,把传播地震波的地球介质视为多层半空间。为简便起见,地震波的衰减问题另作考虑。这样,便需要研究多层、均匀、各向同性和完全弹性半空间中地震震源辐射的地震波传播问题。 用哈斯克尔（Haskell）矩阵法解多层介质中弹性波的传播问题是很方便的。如果     

An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P,S and Rayleigh waves

A boundary integral formulation is presented and applied to model the ground motion on alluvial valleys under incident P, S and Rayleigh waves. It is based on integral representations for the diffracted and the refracted elastic waves using single-layer boundary sources. This approach is called indirect BEM in the literature as the sources' strengths should be obtained as an intermediate step. Boundary conditions lead to a system of integral equations for boundary sources. A discretization scheme based on the numerical and analytical integration of exact Green's functions for displacements and tractions is used. Various examples are given for two-dimensional problems of diffraction of elastic waves by soft elastic inclusion models of alluvial deposits in an elastic half-space. Results are displayed in both frequency and time domains. These results show the significant influence of locally generated surface waves in seismic response and suggest approximations of practical interest. For shallow alluvial valleys the response and its resonant frequencies are controlled by a coupling mechanism that involves both the simple one-dimensional shear beam model and the propagation of surface waves.     

Axisymmetric transient waves in transversely isotropic half-space

A transversely isotropic material in the sense of Green is considered. Using a series of potential functions proposed in [Eskandari-Ghadi M. A complete solution of the wave equations for transversely isotropic media. J Elasticity 2005; 81:1–19], the solutions of the transient wave equations within a half-space under surface load are obtained in the Laplace–Hankel domain for axisymmetric problems. The solutions are investigated in detail in the special case of a surface point force pulse varying with time as Heaviside function. Using Cagniard–De Hoop method, the inverse Laplace transform and inverse Hankel transform of the solutions are then obtained in the form of integrals with finite limits. For validity of the analytical results, the final formulations for surface waves are degenerated for an isotropic material and compared with the existing formulation obtained by Pekeris [The seismic surface pulse. Proc Natl Acad Sci USA 1955;41:469–80], to show that they are exactly the same. The numerical evaluations of the integrals for some transversely isotropic materials as well as an isotropic one are obtained. The present approach is then numerically verified by comparing a particular case of displacements for the surface of an isotropic half-space subjected to a point load of Heaviside function with the solutions obtained by Pekeris [The seismic surface pulse. Proc Natl Acad Sci USA 1955;41:469–80]. In addition, the wave equations for the mentioned medium are obtained on the vertical line directly under the applied surface load. The final formulations are degenerated for an isotropic material and compared with the existing formulation given in Graff [Wave motion in elastic solids. New York: Dover Publications Inc; 1975 [New Ed edition, November 1991]], to show that they are also exactly the same. Then equations are presented in graphical forms using an appropriate numerical evaluation.     

Scattering of vertically-incident P-waves by an embedded pile

An exact theoretical formulation is presented for the analysis of a thin-walled pile embedded in an elastic half-space under vertically-incident P-wave excitation. In the framework of three-dimensional elastodynamics and a shell theory, the axisymmetrical wave-scattering problem is shown to be reducible to a set of Fredholm boundary integral equations. With the incorporation of the singular characteristics of the wave-induced contact load distributions into the solution scheme, a computational boundary element method is developed for a rigorous treatment of the seismic soil-structure interaction problem. Typical results for the dynamic contact load distributions, displacements, complex-valued foundation input motion functions, and resonant pile foundation response are included for direct engineering applications.     

Rocking vibrations of a rigid circular foundation on poroelastic half-space to elastic waves

A study is carefully conducted for the rocking response of a rigid circular foundation resting on a poroelastic half-space when subjected to seismic waves under the framework of Biot’s theory. The free-field waves, rigid-body scattering field waves and radiation scattering field waves are introduced to consider the complex behavior of the soil owing to the scattering phenomena caused by the existence of the foundation. The contact surface between the soil and the foundation is supposed to be perfectly bonded and fully permeable. Combining with the divided wave fields, two sets of dual integral equations elaborating the mixed boundary-value conditions are established, and then reduced to Fredholm integral equations. Therefore, with a semi-analytical method, the expressions of the rocking displacements are obtained. The numerical results of the rocking vibration of the foundation for incident P, SV and Rayleigh waves are presented. The influences of certain parameters, such as the permeability of the soil, the incident angle, Poisson’s ratio and the mass of the foundation, on the rocking vibration of the foundation are explored and studied. Different reactions are found when the foundation is excited by different waves.     

Dynamic structure–soil–structure interaction between nearby piled buildings under seismic excitation by BEM–FEM model

The dynamic through–soil interaction between nearby pile supported structures in a viscoelastic half-space, under incident S and Rayleigh waves, is numerically studied. To this end, a three-dimensional viscoelastic BEM–FEM formulation for the dynamic analysis of piles and pile groups in the frequency domain is used, where soil is modelled by BEM and piles are simulated by one-dimensional finite elements as Bernoulli beams. This formulation has been enhanced to include the presence of linear superstructures founded on pile groups, so that structure–soil–structure interaction (SSSI) can be investigated making use of a direct methodology with an affordable number of degrees of freedom. The influence of SSSI on lateral spectral deformation, vertical and rotational response, and shear forces at pile heads, for several configurations of shear one-storey buildings, is addressed. Maximum response spectra are also presented. SSSI effects on groups of structures with similar dynamic characteristics have been found to be important. The system response can be either amplified or attenuated according to the distance between adjacent buildings, which has been related to dynamic properties of the overall system.     

Dynamic analysis of piles and pile groups embedded in homogeneous soils

A boundary element formulation for the dynamic analysis of axially and laterally loaded single piles and pile groups is presented. The piles are represented by compressible beam-column elements and the soil as a hysteretic elastic half-space. The governing equations of motion for the pile domain have been solved exactly for distributed periodic loading intensities. These solutions are then coupled with a numerical solution for the motion of the soil domain by satisfying equilibrium and compatibility at the pile-soil interface. The results obtained from the analysis compare favourably with those from alternative analyses, e.g. finite element, but at greatly reduced computational costs.     

Dynamic response of rectangular foundations to obliquely incident seismic waves

A study is made of the harmonic response of a rigid massless rectangular foundation bonded to an elastic half-space and subjected to the action of both external forces and obliquely incident plane seismic waves. The associated mixed boundary value problem is discretized and solved numerically. The results obtained indicate that the angle of incidence of the seismic wave has a marked effect on the nature and magnitude of the foundation response.     

多层弹性半空间中的地震波（二）

一、引言 地震面波的频散性质、地震波辐射的方向性等特性已经广泛地用于地壳和上地幔结构以及震源机制的研究中,并且取得了许多有用的成果.研究地震波如何从震源辐射出来、如何在实际介质中传播和衰减的这一问题,对于利用地震波确定地壳和上地幔结构以及震源的参数,是很有必要的.关于这一问题的研究,已经作过许多工作。已往的工作中,为了分析方便,往往采用简单的地壳-上地幔模型或简单的震源模型,或两者都相     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

Kinematic response of single piles to vertically incident P‐waves

Dynamic response of single piles to seismic waves is fundamentally different from the free‐field motion because of the interaction between the pile and the surrounding soil. Considering soil–pile interaction, this paper presents a new displacement model for the steady‐state kinematic response of single piles to vertically incident P‐waves on the basis of a continuum model. The governing equations and boundary conditions of the two undetermined functions in the model are obtained to be coupled by using Hamilton's principle. Then, the two unknown functions are decoupled and solved by an iterative algorithm numerically. A parametric study is performed to investigate the effects of the properties of the soil–pile system on the kinematic response of single piles. It is shown that the effects of the pile–soil modulus ratio, the slenderness ratio of the pile, and the frequency of the incident excitations are very significant. By contrast, the influence of soil damping on the kinematics of the system is slight and can be neglected. Copyright © 2013 John Wiley & Sons, Ltd.     

Dynamic response of 3-D rigid surface foundations by time domain boundary element method

The dynamic response of three-dimensional rigid surface foundations of arbitrary shape is numerically obtained. The foundations are placed on a linear elastic, isotropic and homogeneous half-space representing the soil medium and are subjected to either external dynamic forces or seismic waves of various kinds and directions, with a general transient time variation. The problem is formulated in the time domain by the boundary element method and the response is obtained by a time step-by-step integration. Two examples dealing with three-dimensional rectangular foundations are presented in detail, together with comparisons with other methods, in order to document the accuracy of the method. The main advantages of the proposed method are that, unlike frequency domain techniques, it provides directly the transient response and forms the basis for extension to the case of non-linear behaviour.     

A coupled fluid layer–rigid disk–poroelastic half-space vibration problem

This paper proposes a coupled fluid layer–foundation–poroelastic half-space vibration model to study how still water affects foundations operating underwater. As an example, we consider the problem of the vertical vibration of a rigid disk on a poroelastic half-space covered by a fluid layer having a finite depth. The solution of the disk vibration problem is obtained using the boundary conditions at the free surface of the fluid layer and the boundary conditions at the fluid layer–poroelastic medium interface. The solution is expressed in terms of dual integral equations that are converted into Fredholm integral equations of the second kind and solved numerically. Selected numerical results for the vertical dynamic impedance coefficient are examined based on different water depths, poroelastic materials, disk permeabilities and frequencies of excitation. Based on the numerical results, it is proposed that the hydrodynamic pressure caused by the foundation vibration is the intrinsic reason that the existence of a fluid layer has such a great effect on the dynamic characteristics of the foundation. In many cases, the hydrodynamic pressure caused by the foundation vibration cannot be ignored when designing dynamic underwater foundations. These results are helpful in understanding the dynamic response of foundations under still water without water waves, such as foundations in pools, lakes and reservoirs.     

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

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

Seismic response of foundations embedded in a layered half-space

An integral equation technique to determine the response of foundations embedded in a layered viscoelastic half-space when subjected to various types of seismic waves is presented. The technique is validated by comparison with previous results for rigid hemispherical and cylindrical foundations embedded in a uniform half-space. Illustrative results for rigid cylindrical foundations embedded in layered media are also presented.     

Dynamic stiffness and kinematic response of single piles in inhomogeneous soil

The effect of soil inhomogeneity on dynamic stiffness and kinematic response of single flexural elastic piles to vertically-propagating seismic SH waves is explored. A generalized parabolic function is employed to describe the variable shear wave propagation velocity in the inhomogeneous stratum. A layered soil with piece-wise homogeneous properties is introduced to approximate the continuous inhomogeneity in the realm of a Beam-on-Dynamic-Winkler-Foundation model. The problem is treated numerically by means of a layer transfer-matrix (Haskell–Thompson) formulation, and validated using available theoretical solutions and finite-element analyses. The role of salient model parameters such as pile-head fixity conditions, pile-to-soil stiffness ratio, surface-to-base shear wave velocity ratio and rate of inhomogeneity is elucidated. A new normalization scheme for inertial and kinematic response of such systems is presented based on an average Winkler wavenumber. With reference to long piles in moderately inhomogeneous soils, results indicate that: (a) kinematic pile response is essentially governed by a single dimensionless frequency parameter accounting for pile-to-soil stiffness ratio, pile slenderness and soil inhomogeneity and (b) definition of a characteristic pile wavelength allows an approximate estimation of pile elastodynamic response for preliminary design or analysis. Issues related to domain discretization and Winkler moduli are discussed.     

Soil–pile interaction during the passage of rayleigh waves: An analytical solution

A simple analytical solution is presented to calculate the single-pile response when excited by the passage of Rayleigh seismic waves. Closed-form expressions for the horizontal and vertical displacement distributions are presented for piles with finite or infinite length. The analytical results for both free-head and fixed-head piles are obtained through a dynamic Winkler model, with realistic frequency-dependent ‘springs’ and ‘dashpots’. The results of the presented method are in excellent agreement with results of a rigorous solution. It is shown that in vertical motion, the differences between pile and soil displacements are far more significant than in horizontal motion, and therefore, further work is needed to investigate the importance of pile-soil-pile interaction (group effects), because of the vertical component of Rayleigh seismic waves.     

Response of pile foundations to rayleigh waves and obliquely incident body waves

The paper presents results of a study on the harmonic response of piles and pile groups embedded in a halfspace to various forms of seismic waves. These include the Rayleigh wave as well as obliquely incident P, SV and SH waves. The pertinent mixed boundary value problems of pile-soil-pile interaction are solved by a numerical model of the boundary integral nature. All modes of foundation vibrations, i.e. translational, rocking and torsional, are included in the model. The results presented are used to highlight the salient features of the seismic response of piles. In addition, the influence of certain pile-soil parameters, such as pile rigidity and pile spacing, on the seismic behaviour of pile foundations is investigated.     

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

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

Seismic response of a cylindrical shell embedded in a layered viscoelastic half-space. I: Formulation

A method to obtain the three-dimensional harmonic response of a infinitely long cylindrical shell of circular cross-section embedded in a layered viscoelastic half-space and subjected to harmonic plane waves impinging at an oblique angle with respect to the axis of the shell is presented. The procedure combines an indirect integral representation for the field in the exterior half-space with a model of the pipeline or tunnel based on Donnell shell theory. The integral representation for the soil is based on the use of moving Green's functions for the layered viscoelastic half-space. The accuracy of the formulation is tested by comparison of results obtained by using different discretizations. Extensive comparisons with previous two- and three-dimensional results for the case of a shell embedded in a uniform half-space and some new numerical results for a shell embedded in a multilayered half-space are presented in a companion paper.