Impedance matrices for circular foundation embedded in layered medium

A numerical scheme is developed in the paper for calculating torsional, vertical, horizontal, coupling and rocking impedances in frequency domain for axial-symmetric foundations embedded in layered media. In the scheme, the whole soil domain is divided into interior and exterior domains. For the exterior domain, the analytic solutions with unknown coefficients are obtained by solving three-dimensional (3D) wave equations in cylindrical coordinates satisfying homogeneous boundary conditions. For the interior domain, the analytical solutions are also obtained by solving the same 3D wave equations satisfying the homogeneous boundary conditions and the prescribed boundary conditions. The prescribed conditions are the interaction tractions at the interfaces between embedded foundation and surrounding soil. The interaction tractions are assumed to be piecewise linear. The piecewise linear tractions at the bottom surface of foundation will be decomposed into a series of Bessel functions which can be easily fitted into the general solutions of wave equations in cylindrical coordinates. After all the analytic solutions with unknown coefficients for both interior and exterior domains are found, the variational principle is employed using the continuity conditions (both displacements and stresses) at the interfaces between interior and exterior domains, interior domain and foundation, and exterior domain and foundation to find impedance functions.     

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

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

3-D dynamic foundation-soil-foundation interaction on layered soil

In this work the interaction between adjacent rigid, surface foundations resting on a viscoelastic layered soil medium is studied. A 3-D frequency domain BEM formulation in conjunction with infinite space fundamental solutions and the so called `successive stiffness method', initially developed for elastostatics and adapted here for the solution of elastodynamic problems, are used for the simulation of a layered soil medium. As a result, a discretization of the soil-foundation interface and the surrounding free surface as well as the soil layers' interfaces is necessary. However, it is shown in this work that reasonably accurate results can be obtained by using a substantially reduced discretization scheme involving only a small portion of the free surface surrounding the foundation and the corresponding interfaces of the soil layers. The presented numerical results demonstrate the importance of the dynamic foundation-soil-foundation interaction phenomenon which becomes even more pronounced where the supporting soil medium is made up of relatively shallow layers close to its free surface.     

Transmitting boundary for water-saturated transversely isotropic strata based on the u–U formulation

A new transmitting boundary in a cylindrical coordinate system has been developed for modeling the elastic waves radiating out to an infinite boundary in water-saturated transversely isotropic soil strata over a rigid bedrock. The saturated soil strata are assumed to consist of a porous material and modeled as a transversely isotropic two-phase medium, based on the u−U formulation. The newly developed transmitting boundary is combined with the finite-elements model of the near-field region, using the same u−U formulation, and applied to the study of the dynamics of a rigid circular foundation in porous isotropic or transversely isotropic layered strata, either fully or partly saturated with water. The verification and application examples give valuable insights into new and interesting aspects of the dynamic behavior of rigid circular foundations in fully or partly saturated two-phase ground in terms of permeability, transverse anisotropy, and ground-water table level.     

Dynamic response of a multi-layered poroelastic medium

An exact stiffness matrix method is presented to evaluate the dynamic response of a multi-layered poroelastic medium due to time-harmonic loads and fluid sources applied in the interior of the layered medium. The system under consideration consists of N layers of different properties and thickness overlying a homogeneous half-plane or a rigid base. Fourier integral transform is used with respect to the x-co-ordinate and the formulation is presented in the frequency domain. Fourier transforms of average displacements of the solid matrix and pore pressure at layer interfaces are considered as the basic unknowns. Exact stiffness (impedance) matrices describing the relationship between generalized displacement and force vectors of a layer of finite thickness and a half-plane are derived explicitly in the Fourier-frequency space by using rigorous analytical solutions for Biot's elastodynamic theory for porous media. The global stiffness matrix and the force vector of a layered system is assembled by considering the continuity of tractions and fluid flow at layer interfaces. The numerical solution of the global equation system for discrete values of Fourier transform parameter together with the application of numerical quadrature to evaluate inverse Fourier transform integrals yield the solutions for poroelastic fields. Numerical results for displacements and stresses of a few layered systems and vertical impedance of a rigid strip bonded to layered poroelastic media are presented. The advantages of the present method when compared to existing approximate stiffness methods and other methods based on the determination of layer arbitrary coefficients are discussed.     

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

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

三维层状黏弹性半空间中球面SH、P和SV波源自由场

采用刚度矩阵方法结合Hankel积分变换,求解了层状黏弹性半空间中球面SH、P和SV波的自由波场.首先,在柱坐标系下建立层状黏弹性半空间的反轴对称(柱面SH波)和轴对称(柱面P-SV波)情况精确动力刚度矩阵.进而由Hankel变换将空间域内的球面波展开为波数域内柱面波的叠加,然后将球面波源所在层的上下端面固定,求得固定层内的动力响应和固定端面反力,将固端反力反向施加到层状黏弹性半空间上,采用直接刚度法求得固端反力的动力响应,叠加固定层内和固端反力动力响应,求得波数域内球面波源动力响应.最后由Hankel积分逆变换求得频率-空间域内球面波源自由场,时域结果由傅里叶逆变换求得.文中验证了方法的正确性,并以均匀半空间和基岩上单一土层中球面SH、P和SV波为例分别在频域和时域内进行了数值计算分析.研究表明基岩上单一土层中球面波自由场与均匀半空间情况有着本质差异;基岩上单一土层中球面波位移频谱峰值频率与场地固有频率相对应,基岩面的存在使得基岩上单一土层地表点的位移时程非常复杂,振动持续时间明显增长;阻尼的增大显著降低了动力响应的峰值,同时也显著减少了波在土层的往复次数.     

Vertical vibration of three-dimensional rigid foundations on layered media

A numerical method of analysis is presented for the determination of the steady-state vertical vibration of rigid foundations with arbitrary three-dimensional geometries resting on the surface of a layered soil medium. The method utilizes the flexibility concept applied to steady-state periodic problems and it is solved in the frequency domain. The accuracy of the method is verified by comparison with several published solutions for massless, smooth rigid rectangular foundations on a homogeneous, isotropic elastic half-space. Parametric solutions are presented to study the dynamic behaviour of massless, smooth rigid rectangular foundations on a homogeneous, elastic stratum.     

Dynamic stress concentration near a fluid-filled permeable borehole induced by general modal vibrations of an internal cylindrical radiator

Stress distribution in the vicinity of a permeable cylindrical cavity surface (borehole wall) arising due to modal vibrations of an internal cylindrical radiator of infinite extent is studied. Biot phenomenological model is used to represent the behavior of sound in the fluid-saturated elastic porous medium and closed-form solution in the form of an infinite series is developed. A numerical example for the infinite cylindrical surface excited in vibrational modes of zeroth and first order while immersed in a water-filled cavity embedded within a water-saturated Ridgefield Sandstone environment is presented and several limiting cases are examined. Effects of axial and radial vibration frequencies, porosity, frame stiffness, and interface permeability condition on stress distribution at the borehole surface are presented and discussed.     

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

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

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

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

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

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

饱和土中深埋圆柱形衬砌洞室对瞬态平面波的散射

基于Biot饱和多孔介质动力学理论,运用Laplace变换和波函数展开法,根据饱和土体与衬砌结构交界面的连续条件和衬砌结构内边界上的应力自由条件,得到饱和土中深埋圆柱形衬砌洞室对瞬态平面P波和SV波散射问题的解答,该解答可以退化成为饱和土中深埋圆柱形空穴或弹性夹塞物的情形,并很容易转换成为对稳态波散射的解.通过与已有的相关问题的解析解答进行对比,验证了该解答的正确性.同时利用Laplace逆变换的数值方法,给出了饱和土和衬砌中应力和位移场在时域内的数值解,通过算例,分析了衬砌厚度、刚度对衬砌内边界处应力集中因子的影响.     

Surface motion of a cylindrical hill of circular-arc cross-section for incident plane SH waves

A closed-form solution of two-dimensional scattering of plane SH waves by a cylindrical hill of circular-arc cross-section in a half-space is presented using the wave functions expansion method. The solution is reduced to solving a set of infinite linear algebraic equations, using the auxiliary functions and the exterior region form of the Graf's addition theorem. Numerical solutions are obtained by truncation of the infinite equations and their accuracies are demonstrated by convergence of the numerical results and by the extent to which the numerical results fit the exact boundary conditions with increasing the truncation order. The numerical results for some typical cases are then presented for checking accuracies of various numerical methods. The effects of the height-to-width ratio of the hill on surface ground motion are finally illustrated.     

Scattering of SV waves by a canyon in a fluid-saturated, poroelastic layered half-space, modeled using the indirect boundary element method

The scattering of SV waves by a canyon in a fluid-saturated, poroelastic layered half-space is modeled using the indirect boundary element method in the frequency domain. The free-field responses are calculated to determine the displacements and stresses at the surface of the canyon, and fictitious distributed loads are then applied at the surface of the canyon in the free field to calculate the Green's functions for displacements and stresses. The amplitudes of the fictitious distributed loads are determined from the boundary conditions, and the displacements arising from the waves in the free field and from the fictitious distributed loads are summed to obtain the solution. The effects of fluid saturation, boundary conditions, porosity, and soil layers on the surface displacement amplitudes and phase shifts are discussed, and some useful conclusions are obtained. It is shown that the surface displacement amplitudes due to saturation and boundary conditions, different porosities, or the presence of a soil layer can be very dissimilar, and large phase shifts can be observed. The resulting wavelengths for an undrained saturated poroelastic medium are slightly longer than those for a drained saturated poroelastic medium; and are longer for a drained saturated poroelastic medium than those for a dry poroelastic medium. As porosity increases, the wavelengths become longer; and a layered half-space produces longer wavelengths than a homogeneous half-space.     

A new analytical approach to reconstruct the acceleration time history at the bedrock base from the free surface signal records

Acceleration time histories of earthquake events are typically measured in seismic stations that are placed close to the soil top surface. These acceleration records are often used as input data for seismic analysis. It may be used for base excitation in seismic analysis of above ground structures with shallow foundations.. However it may not be used for seismic analysis of underground structures, or even for above ground buildings with deep foundations and several underground stories. The required base excitation data of the latter should have been measured below the top surface, at a level that may be determined according to the specific analyzed building geometry or at the bedrock below. If the acceleration time history at the bedrock would have been known, the seismic wave propagation through the soil medium, from the bedrock towards the top surface, could have been carried out and the base excitation of the buried structure could be determined. Since there is no data on the acceleration time history at the bedrock, and the only given data is the acceleration records at the top surface, the goal of this paper is to provide an exact reverse analysis procedure to determine the unknown acceleration time history at the bedrock that would exactly produce the measured acceleration time history at the top surface. Once this goal is achieved, seismic analysis of buried structures may be carried out with the determined acceleration record at the bedrock as input. This paper presents an analytical exact solution of the inverse problem for determination of the acceleration, velocity and displacement time histories at the bedrock base of a layered geological medium that are compatible with the given acceleration record at the soil top surface. This new proposed method is based on analytical solutions of the initial-boundary value problems of the linear wave equation in the case of a layered medium. The relationship between waves in one layer and waves in another adjacent layer is derived considering the continuity of stresses and displacements at the common interface between the layers. The efficiency and accuracy of the proposed method is demonstrated through several examples involving the nonstationary response of the free surface. The case of the San Fernando Earthquake is studied. Excellent agreement is achieved between the recorded free surface time history and the reconstructed signal. This excellent agreement is obtained due to the exact analytical method used in deriving the inverse problem solution. This exact analytical method allows one to obtain an acceleration (velocity/displacement) distribution along all the layers at any time.     

Dynamic response of axisymmetric embedded foundations

The boundary element method is used to obtain dynamic stiffness functions of rigid cylindrical foundations embedded in a uniform or layered viscoelastic half-space. Dynamic stiffness functions of hemispherical foundations embedded in a uniform half-space are also computed. The direct integral equation formulation is used in combination with the complete space point load fundamental solution that is integrated numerically along the azimuthal coordinate. The approach is easy to implement because of the simplicity of the fundamental solution. The numerical results obtained by this method for cylindrical and hemispherical foundations are very close to corresponding published results obtained by different procedures. A parametric study shows the important effects of the Poisson's ratio on the dynamic stiffness functions of cylindrical foundations embedded in a uniform viscoelastic half-space. The effect of the bedrock compliance on the stiffness functions is also shown in the case of cylindrical foundations embedded in a soil layer that rests on a bedrock.     

Dynamic plate-soil interaction - finite and infinite, flexible and rigid plates on homogeneous, layered or Winkler soil

An integral method to calculate the solution of a homogeneous or layered soil due to a harmonic point load is described. An infinite plate at the surface of the soil can be introduced in this integration in wavenumber domain, too. Finite structures on the soil are calculated by a combined finite element and boundary element method, which makes use of the point load solution of the soil. The compliance functions for a vertical point load and some vibration modes are calculated for realistic parameters of the plate and the soil and for a wide range of frequencies. The influence of the stiffness of the soil and the foundation is investigated, showing that the soil mainly affects the low-frequent response whereas the structural properties are more important at higher frequencies. A rigid approximation of flexible plates is only found at low frequencies, if the elastic length is used as the radius of a rigid disk. At higher frequencies, a characteristic behaviour of the flexible plate of approximately is observed, what is in clear contrast to the compliance of rigid foundations. A plate on a visco-elastic support (Winkler soil) shows similar displacements as a plate on a homogeneous half-space, but the maximal stresses between the plate and the soil are considerably smaller which is found to be more realistic for a plate on a layered soil. For practical applications, a normalized diagram and some explicit formulas of the exact and the approximate solutions of an infinite plate on a homogeneous half-space are given, which is a useful model to represent the soil-structure interaction of flexible foundations.     

饱和层状场地中任意凹陷地形对入射SV波的散射

利用饱和土层的精确动力刚度矩阵和动力格林(Green)函数,采用间接边界元法,在频域内求解了层状饱和场地中任意凹陷地形对入射SV波的散射问题.通过自由场反应分析,求得凹陷地形表面各点的位移和各单元的应力响应;同样计算了虚拟分布荷载的格林影响函数,求得相应的位移和应力响应;根据边界条件确定虚拟分布荷载,将自由场位移响应和...     

The influence of the degree of saturation on dynamic response of a cylindrical lined cavity in a nearly saturated medium

In most previous studies on the dynamic response of a long cylindrical cavity subjected to internal transient dynamic loads, the porous medium was usually assumed to be completely saturated by ground water. In practice, however, the full saturation condition does not always exist. In this paper the surrounding soil and the lining of the cavity are respectively treated as a nearly saturated porous medium and an elastic material, and the governing equations for the dynamic problem are derived. A set of exact solutions are obtained in the Laplace transform domain for three types of transient loads, i.e. suddenly applied constant load, gradually applied step load and triangular pulse load. By utilizing a reliable numerical method of inverse Laplace transforms, the time-domain solutions are then presented. The influence of the degree of saturation of the surrounding soil on the dynamic response of the lined cavity is examined for numerical examples.