Impedances of embedded rigid strip foundations

This paper is concerned with the dynamic response of rigid strip foundations of arbitrary geometry embedded in a homogeneous elastic half-space. The embedded rigid foundation is modelled by an equivalent domain in a uniform half-space which is subjected to an appropriate body force field. The components of the impedance matrix are determined through the solution of a linear simultaneous equation system which is established by invoking rigid body displacements of discrete locations within the equivalent domain and appropriate equilibrium consideration. It is found that high numerical efficiency and flexibility can be achieved using the body force model when compared to boundary integral formulations through the selection of appropriate displacement influence functions and a ‘parent domain’ in the analysis. Numerical results are presented to illustrate the influence of the embedment ratio, frequency of excitation, foundation geometry and Poisson's ratio on the vertical, horizontal, rocking and coupled impedances of a single embedded foundation. The effect on the impedance due to the presence of an adjacent embedment is investigated for various distances between foundations and embedment ratios.     

2-D transient soil–surface foundation interaction and wave propagation by time domain BEM

This paper shows an effective implementation of the half-plane Green function for surface strip impulses (Lamb's problem), which was previously developed in a closed form by the authors, into the time-domain boundary element method for the analysis of related initial boundary value problems. The time-stepping algorithm utilizing Heaviside step function makes the solution process free from the Rayleigh wave front singularity. Illustrative analyses performed include that: First, the response due to an impulsive uniform strip loading is dealt with in order to check the accuracy of the present solution and to interpret the associated wave motion in the medium. Second, a rigid massless strip surface foundation is analysed when subjected to various impulsive loadings in vertical, horizontal and rotational directions to observe which wave is most concerned with the respective foundation motion. The field response is also of interest with respect to distance attenuation. Third, the dynamic cross-interaction between active and passive foundations through soil is investigated when multiple strip foundations are placed separately on a half-space with a certain distance.     

Vibration of surface foundations of arbitrary shapes

Presented is a systematic procedure for generating impedance (or compliance) matrices for foundations with arbitrary shapes, resting on an elastic half-space medium. A technique to decompose prescribed harmonic tractions on the half-space medium is employed to solve analytically the differential wave equations in cylindrical coordinates. However, the interaction stresses due to the vibration of a foundation with arbitrary shape are described in rectangular coordinates, and assumed to be piecewise constant in the region of the arbitrary shape. A coordinate transformation matrix is introduced for the piecewise constant tractions in order to use the solution of the differential wave equations in cylindrical coordinates. Finite element modelling is assumed in rectangular coordinates for the foundation itself. The impedance matrix is then obtained for the finite element model, using a variational principle and the reciprocal theorem. A simple example of a rigid square plate resting on a half-space medium and subjected to vertical excitation is used to demonstrate the efficiency and effectiveness of the procedure. Some numerical aspects are investigated and some possible extensions of the procedure are also discussed.     

Dynamic interaction between rigid foundations in a layered half-space

A computationally efficient boundary integral equation technique to calculate the dynamic response of a group of rigid surface foundations bonded to a layered viscoelastic half-space and subjected to external forces and seismic waves is presented. The technique relies on an iterative scheme which minimizes in-core memory requirements and takes advantage of any geometrical symmetry of the foundations. Extensive results for the case of two rigid square foundations placed at different separations and bonded to a viscoelastic half-space are presented. It was found that the choice of discretization of the foundations has a marked effect on the calculated impedance functions for extremely small separations. Illustrative results for a case of several closely-spaced foundations bonded to a layered half-space are also presented.     

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 response of rigid foundations of arbitrary shape

An approximate numerical procedure for calculation of the harmonic force-displacement relationships for a rigid foundation of arbitrary shape placed on an elastic half-space is presented. This procedure is used to evaluate the vertical, rocking and horizontal compliance functions for rigid rectangular foundations and the vertical compliance for a rigid square foundation with an internal hole. Several comparisons between the results obtained by the proposed approach and other methods are also presented.     

Discrete models for vertical vibrations of surface and embedded foundations

A five-parameter discrete model that approximates the dynamic force4isplacement relationship for rigid foundations undergoing vertical vibrations on a uniform elastic half-space is presented. The model involves a combination of two springs, two viscous dampers and a mass. Values of the parameters for circular, square and rectangular foundations placed on the surface or embedded in an elastic half-space are listed. The parameters are obtained by minimizing the discrepancy between the force4isplacement relation for the model and that obtained by solution of the mixed boundary-value problem of the rigid foundation on an elastic half-space. The definition of an appropriate input motion to represent wave excitation is also discussed. The input motion to the discrete model differs from the input motion that should be used in a continuum model.     

Dynamic stiffness matrices for two circular foundations

A systematic procedure is presented for generating dynamic stiffness matrices for two independent circular foundations on an elastic half-space medium. With the technique reported in References 1–3, the analytic solution of three-dimensional (3D) wave equations satisfying the prescribed traction due to the vibration of one circular foundation can be found. Since there are two analytic solutions for two prescribed tractions due to the vibrations of two circular foundations, the principle of superposition must be used to obtain the total solution. The interaction stresses (prescribed tractions) are assumed to be piecewise linear in the r-directions of both cylindrical co-ordinates for the two circular foundations. Then, the variational principle and the reciprocal theorem are employed to generate the dynamic stiffness matrices for the two foundations. In the process of employing the variational principle, a co-ordinate transformation matrix between two cylindrical co-ordinate systems is introduced. Some numerical results of dynamic stiffness matrices for the interaction of two identical rigid circular foundations are presented in order to show the effectiveness and efficiency of the present method, and some elaborations for its future extensions are also discussed.     

埋置基础扭转振动的实用化计算

根据弹性半空间理论及基础振动试验和数值计算的最新成果,采用方程对等法推演出埋置块体基础扭转振动的实用化计算公式。采用该实用化计算公式和弹性半空间理论分别对复杂形状的埋置块体基础在谐和扰力矩作用下的扭转振动动力响应进行计算,并将计算结果进行了比较,两者结果吻合,振幅误差为9.13%。该实用化计算公式具有概念明确,计算简单的优点,对于任意形状、任意埋置状况的动力设备块体基础都实用,对于结构工程的地震响应计算具有参考意义。     

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.     

Application of a second-order paraxial boundary condition to problems of dynamics of circular foundations on a porous layered half-space

A half-space finite element and a consistent transmitting boundary in a cylindrical coordinate system are developed for analysis of rigid circular (or cylindrical) foundations in a water-saturated porous layered half-space. By means of second-order paraxial approximations of the exact dynamic stiffness for a half-space in plane-strain and antiplane-shear conditions, the corresponding approximation for general three-dimensional wave motion in a Cartesian coordinate system is obtained and transformed in terms of cylindrical coordinates. Using the paraxial approximations, the half-space finite element and consistent transmitting boundary are formulated in a cylindrical coordinate system. The development is verified by comparison of dynamic compliances of rigid circular foundations with available published results. Examination of the advantage of the paraxial condition vis-á-vis the fixed condition shows that the former achieves substantial gain in computational effort. The developed half-space finite element and transmitting boundary can be employed for accurate and effective analysis of foundation dynamics and soil–structure interaction in a porous layered half-space.     

Dynamic vertical compliance of strip foundations in layered soils

In this paper, results of a detailed investigation on the dynamic response of rigid strip foundations in viscoelastic soils under vertical excitation are presented. An advanced boundary element algorithm developed by incorporating isoparametric quadratic elements and a sophisticated self-adapting numerical integration scheme has been used for this investigation. Foundations supported on three types of soil profiles, namely, homogeneous half-space, stratum-over-half-space and stratum-over-bedrock are considered. The influence of material properties like Poisson's ratio and material damping as well as the influence of geometrical properties such as depth of embedment and layer thickness are studied. The effect of the type of contact at the soil-foundation interface is also investigated.     

Dynamic interaction between flexible strip foundations

This paper is concerned with the dynamic response of a system of flexible strip foundations resting in smooth contact with a homogeneous isotropic viscoelastic half space. An arbitrary number of foundations with different flexibilities and geometries subjected to time-harmonic distributed loadings are considered in the formulation. The response of each strip foundation is governed by the classical plate theory and its transverse deflection is represented by an admissible function containing a set of generalized co-ordinates. A coupled variational-Green's function scheme is employed to establish the equations of motion of the strip foundation system. The numerical stability and convergence of the solution scheme are established. The influence of the foundation flexibility, distance between adjacent foundations and frequency of motion on the response of the foundation system is investigated in the numerical study.     

Effects of the dynamic cross-interaction in the seismic analysis of multiple embedded foundations

A boundary element formulation of the substructure deletion method is presented for the seismic analysis of the dynamic cross-interaction between multiple embedded foundations. This approach is particularly suitable for three-dimensional foundations of any arbitrary geometrical shape and spatial location, since it requires only the discretization of the foundations’ surfaces. The surrounding soil is represented by a homogeneous viscoelastic half-space while the foundations are assumed to be rigid and subjected to incoming SH-, P-, and SV-waves arbitrarily inclined in both the horizontal and vertical planes. The proposed methodology is tested for the case of two identical embedded square foundations for different values of the foundations’ embedment and distance. The effects of the cross-interaction are outlined in the components of the impedance matrix and of the foundation input motion. © 1997 John Wiley & Sons, Ltd.     

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.     

Dynamics of rigid strip foundations embedded in orthotropic elastic soils

This study is concerned with the dynamic response of an arbitrary shaped rigid strip foundation embedded in an orthotropic elastic soil. The foundation is subjected to time-harmonic vertical, horizontal and moment loadings. The boundary-value problem related to an embedded foundation is analysed by using the indirect boundary integral equation method. The kernel functions of the integral equations are displacement and traction Green's functions of an anisotropic elastic half plane. Exact analytical solutions are used for the Green's functions. The boundary integral equation is solved by using numerical techniques. Selected numerical results are presented for the impedances of rectangular and semi-circular rigid strip foundations embedded in four types of anisotropic soils. A discussion on the influence of soil anisotropy and frequency of excitation on the impedances is presented. The versatility of the analysis is demonstrated by considering the through soil interaction between two semi-circular strip foundations.     

Dynamic interaction between 3-D rigid surface foundations and comparison with the ATC-3 provisions

The dynamic behaviour of a system of three-dimensional, massless, rigid, surface foundations of arbitrary shape perfectly bonded to the elastic half-space is numerically studied with the frequency domain boundary element method. This method employs the dynamic Green's function for the surface of the half-space and this results in a discretization of only the soil-foundation interfaces. In addition, use of isoparametric quadratic quadrilateral boundary elements increases the accuracy of the method, which is confirmed by comparison with other known numerical solutions. Externally applied loads, harmonically varying with time, are considered. The through the soil coupling effect between the foundations as a function of distance and frequency is assessed through extensive parametric studies involving two and four rigid foundations being isolated or interconnected. It is found that the assertion of ATC-3 regulations that omission of coupling effects leads to conservative results is not always correct for all frequencies.     

Dynamic stiffness of rigid rectangular foundations on the half-space

The dynamic soil–structure interaction of a rigid rectangular foundation with the subsoil represents a mixed-boundary value problem. This problem is formulated in terms of a system of coupled Fredholm integral equations of the first kind. The subsoil is modelled by a homogeneous, linear-elastic and isotropic half-space which is perfectly bonded to the rigid, rectangular foundation. An approximate solution for the resultant loads between the foundation and the half-space due to a unit forced displacement or rotation is obtained using the Bubnov–Galerkin method. Using this method the displacement boundary value conditions are exactly satisfied and the contact stress distributions between the foundation and the half-space are approximated by series expansions of Chebyshev polynomials. This method provides a simple means of studying the soil-structure interaction of rectangular foundations with different inertia properties.     

Cone models for nearly incompressible soil

Cones can be used to model soil in a unified strength-of-materials approach. For the vertical and rocking motions involving predominantly compressional-extensional deformation, the corresponding dilatational wave velocity tends to infinity for Poisson's ratio approaching 1/2. Based on the rigorous solution for the dynamic stiffness of a rigid disk for all frequencies, whereby the partition of the power among P-, S- and Rayleigh waves is also discussed, two special features are necessary for the vertical and rocking motions for nearly incompressible soil with Poisson's ratio between 1/3 and 1/2: (1) The appropriate wave velocity is selected as twice the shear wave velocity and not as the dilatational wave velocity; (2) A trapped mass which increases linearly with Poisson's ratio is introduced. The trapped mass can be assigned to the base mat, allowing the cone model to be constructed in the same way for all Poisson's ratios. The realization of cone models for surface foundations on a homogeneous half-space and on a layer on a flexible half-space and for embedded and pile foundations is addressed.