Some analytical results on lateral dynamic stiffness for footings supported on hysteretic soil medium

The nonlinearity of the soil affects soil–structure interaction to a considerable extent. For a reliable and safe analysis of soil interaction effects on the dynamic response of structures, a more realistic and relatively straightforward method incorporating the nonlinear hysteretic nature of the underlying soil–foundation system needs to be developed. The present paper models the soil–foundation system as a single degree of freedom spring–dashpot system with nonlinear hysteresis in form of elasto-perfectly plastic behavior. Analytical results for the lateral dynamic stiffness on footing have been presented. An example study has been carried out in case of circular footings. It is shown how the analytical results can be used to get a preliminary idea of the lateral dynamic stiffness of footings on a soil medium prior to a detailed computational geo-mechanics analysis provided the static nonlinear load–deformation characteristic of the soil medium is known and can be modeled by a hysteretic elasto-plastic behavior. The corresponding results are presented in a graphical form. The results have been computed showing parametric variations with the change in the amplitude and dimensionless frequency of the non-dimensional excitation force. Analytical results are also presented for the asymptotic cases at low and very high values of dimensionless frequency parameter.     

Dynamic soil–structure interaction effects on the seismic response of asymmetric buildings

An approach is formulated for the linear analysis of three-dimensional dynamic soil–structure interaction of asymmetric buildings in the time domain, in order to evaluate the seismic response behaviour of torsionally coupled buildings. The asymmetric building is idealized as a single-storey three-dimensional system resting on different soil conditions. The soil beneath the superstructure is modeled as linear elastic solid elements. The contact surface between foundation mat and solid elements of soil is discretised by linear plane interface elements with zero thickness. An interface element is further developed to function between the rigid foundation and soil. As an example, the response of soil–structure interaction of torsionally coupled system under two simultaneous lateral components of El Centro 1940 earthquake records has been evaluated and the effects of base flexibility on the response behaviour of the system are verified.     

Seismic soil–structure interaction of massive flexible strip-foundations embedded in layered soils by hybrid BEM–FEM

A study on the seismic response of massive flexible strip-foundations embedded in layered soils and subjected to seismic excitation is presented. Emphasis is placed on the investigation of the system response with the aid of a boundary element–finite element formulation proper for the treatment of such soil–structure interaction problems. In the formulation, the boundary element method (BEM) is employed to overcome the difficulties that arise from modeling the infinite soil domain, and the finite element method (FEM) is applied to model the embedded massive flexible strip-foundation. The numerical solution for the soil–foundation system is obtained by coupling the FEM with the BEM through compatibility and equilibrium conditions at the soil–foundation and soil layer interfaces. A parametric study is conducted to investigate the effects of foundation stiffness and embedment on the seismic response.     

Seismic analysis of lock–soil–fluid systems by hybrid BEM–FEM

A study of soil–structure–fluid interaction (SSFI) of a lock system subjected to harmonic seismic excitation is presented. The water contained lock is embedded in layered soils supported by a half-space bedrock. The ground excitation is prescribed at the soil–bedrock interface. The response is numerically obtained through a hybrid boundary element (BEM) finite element method (FEM) formulation. The semi-infinite soil and the fluid are modeled by the BEM and the lock is modeled by the FEM. The equilibrium equation for the lock system is obtained by enforcing compatibility and equilibrium conditions at the fluid–structure, soil–structure and soil–layer interfaces under conditions of plane strain. To the authors’ knowledge this is the first study of a lock system that considers the effects of dynamic soil–fluid–structure interaction through a BEM–FEM methodology. A numerical example and parametric studies are presented to examine the effects of the presence of water, lock stiffness, and lock embedment on the response.     

Soil–structure interaction in the time domain using halfspace Green's functions

A time-domain formulation is proposed for the transient response analysis of general, three-dimensional structures resting on a homogeneous, elastic halfspace subjected to either external loads or seismic motions. The formulation consists of two parts: (a) the time domain formulation of the soil behaviour and (b) the coupling of the corresponding soil algorithms to the Finite Element Code ANSYS. As far as the structure is concerned, this coupling opens the way for the analysis of non-linear soil–structure interaction. The approach is based on halfspace Green's functions for displacements elicited by Heaviside time-dependent surface point loads. Hence, the spatial discretisation can be confined to the contact area between the foundation and the soil, i.e. no auxiliary grid beyond the foundation as for conventional boundary element formulations is required. The method is applied to analyse the dynamic response of a railway track due to a moving wheel set by demonstrating the influence of ‘through-the-soil coupling’.     

Investigation of seismic behavior of fluid–rectangular tank–soil/foundation systems in frequency domain

Main purpose of this study is to evaluate the dynamic behavior of fluid–rectangular tank–soil/foundation system with a simple and fast seismic analysis procedure. In this procedure, interaction effects are presented by Housner's two mass approximations for fluid and the cone model for soil/foundation system. This approach can determine; displacement at the height of the impulsive mass, the sloshing displacement and base forces for the soil/foundation system conditions including embedment and incompressible soil cases. Models and equations for proposed method were briefly explained for different tank–soil/foundation system combinations. By means of changing soil/foundation conditions, some comparisons are made on base forces and sloshing responses for the cases of embedment and no embedment. The results showed that the displacements and base shear forces generally decreased, with decreasing soil stiffness. However, embedment, wall flexibility, and soil–structure interaction (SSI) did not considerably affect the sloshing displacement.     

Dynamic soil–structure interaction analysis via coupled finite-element–boundary-element method

In this paper, a study on the transient response of an elastic structure embedded in a homogeneous, isotropic and linearly elastic half-plane is presented. Transient dynamic and seismic forces are considered in the analysis. The numerical method employed is the coupled Finite-Element–Boundary-Element technique (FE–BE). The finite element method (FEM) is used for discretization of the near field and the boundary element method (BEM) is employed to model the semi-infinite far field. These two methods are coupled through equilibrium and compatibility conditions at the soil–structure interface. Effects of non-zero initial conditions due to the pre-dynamic loads and/or self-weight of the structure are included in the transient boundary element formulation. Hence, it is possible to analyse practical cases (such as dam–foundation systems) involving initial conditions due to the pre-seismic loads such as water pressure and self-weight of the dam. As an application of the proposed formulation, a gravity dam has been analysed and the results for different foundation stiffness are presented. The results of the analysis indicate the importance of including the foundation stiffness and thus the dam–foundation interaction.     

Numerical study of characteristics of underground blast induced surface ground motion and their effect on above-ground structures. Part II. Effects on structural responses

Studies of structural responses and damage to high-frequency blast motion are very limited. Current practice uses some empirical allowable ground vibration limits in assessing structural performance. These empirical limits overlook the physical parameters that govern structural response and damage, such as the ground motion characteristics and inherent structural properties. This paper studies the response of RC frame structures to numerically simulated underground blast-induced ground motions. The structural response and damage characteristics of frame structures to ground motions of different frequencies are investigated first. The effects of blast ground motion spatial variations and soil–structure interaction on structural responses are also studied. A suitable discrete model that gives accurate response prediction is determined. A damage index defined based on the accumulated plastic hinge rotation is used to predict structural damage level. Numerical results indicated that both the low structural vibration modes (global modes) and the first elemental vibration mode (local) might govern the dynamic structural responses depending on the ground motion frequency and structural response parameters under consideration. Both ground motion spatial variations and soil–structure interaction effects are prominent. Neglecting them might yield inaccurate structural response prediction. The overall structural response and damage are highly ground motion frequency dependent. Numerical results of structural damage are also compared with some test results obtained in a previous study and with code specifications. Discussions on the adequacy of the code allowable ground vibration limits on RC frame structures are also made.     

Real-time dynamic hybrid testing for soil–structure interaction analysis

Simulating dynamic soil–structure interaction (SSI) problems is a challenge when using a shaking table because of the semi-infinity of soil foundations. This paper develops real-time dynamic hybrid testing (RTDHT) for SSI problems in order to consider the radiation damping effect of the semi-infinite soil foundation using a shaking table. Based on the substructure concept, the superstructure is physically tested and the semi-infinite foundation is numerically simulated. Thus, the response of the entire system considering the dynamic SSI is obtained by coupling the numerical calculation of the soil and the physical test of the superstructure. A two-story shear frame on a rigid foundation was first tested to verify the developed RTDHT system, in which the top story was modeled as the physical substructure and the bottom story was the numerical substructure. The RTDHT for a two-story structure mounted on soil foundation was then carried out on a shaking table while the foundation was numerically simulated using a lumped parameter model. The dynamic responses, including acceleration and shear force, were obtained under soft and hard soil conditions. The results show that the soil–structure interaction should be reasonably taken into account in the shaking table testing for structures.     

Ambient vibration tests of a seven-story reinforced concrete building in Van Nuys, California, damaged by the 1994 Northridge earthquake

Results of two detailed ambient vibration surveys of a 7-story reinforced concrete building in Van Nuys, California, are presented. Both surveys were conducted after the building was severely damaged by the 17 January 1994, Northridge earthquake (M_L=5.3, epicenter 1.5 km west from the building site) and its early aftershocks. The first survey was conducted on 4 and 5 February 1994, and the second one on 19 and 20 April 1994, about one month after the 20 March aftershock (M_L=5.3, epicenter 1.2 km north–west from the building site). The apparent frequencies and two- and three-dimensional mode shapes for longitudinal, transverse and vertical vibrations were calculated. The attempts to detect the highly localized damage by simple spectral analyses of the ambient noise data were not successful. It is suggested that very high spatial resolution of recording points is required to identify localized column and beam damage, due to the complex building behavior, with many interacting structural components. The loss of the axial capacity of the damaged columns could be seen in the vertical response of the columns, but similar moderate or weak damage typically would not be noticed in ambient vibration surveys. Previous analysis of the recorded response of this building to 12 earthquakes suggests that, during large response of the foundation and piles, the soil is pushed sideways and gaps form between the foundation and the soil. These gaps appear to be closing during “dynamic compaction” when the building site is shaken by many small aftershocks. The apparent frequencies of the soil–foundation–structure system appear to be influenced significantly by variations in the effective soil–foundation stiffness. These variations can be monitored by a sequence of specialized ambient vibration tests.     

Structure, soil–structure response and effects of damage based on observations of horizontal-to-vertical spectral ratios of microtremors

The microtremor horizontal-to-vertical-spectral-ratio (HVSR) technique is widely used in the urban environment to assess the fundamental frequency response of the ground. Extensive literature exists about case histories using HVSR for microzonation in several cities, but no systematic studies have been devoted to check the presence of soil–structure interaction effects, and even less attention to study building behaviour after earthquake damage. To evaluate the above-mentioned effects, a series of experiments are reported in this article.We first made a series of microtremor measurements on buildings and civil structures to evaluate the reliability of fundamental frequency determinations. Then, we considered several case studies to evaluate the effect of soil–structure interaction in estimates of site response in the presence of tall buildings. Finally, an experiment on the frequency change due to damage was performed. It was possible to confirm that HVSR is able to detect building fundamental modes and once known the building frequency, it is also possible to detect the presence of soil–structure interaction. Thus, once the presence of the building natural frequency is identified, it is possible to infer the site response from free field measurements. We also found that the HVSR technique is equally useful for detecting structural damage by determining the frequency shift of the buildings.     

Wave propagation and soil–structure interaction on a cliff crest during the 1999 Athens Earthquake

During the 1999 Athens Earthquake the town of Adàmes, located on the eastern cliff of the Kifissos river canyon, experienced unexpectedly heavy damage. Despite the significant amplification potential of the slope geometry, topography effects cannot alone explain the uneven damage distribution within a 300 m zone behind the crest, characterized by a rather uniform structural quality. This paper illustrates the important role of soil stratigraphy, material heterogeneity, and soil–structure interaction on the characteristics of ground surface motion. For this purpose, we first perform elastic two-dimensional wave propagation analyses utilizing available geotechnical and seismological data, and validate our results by comparison with aftershock recordings. We then conduct non-linear time-domain simulations that include spatial variability of soil properties and soil–structure interaction effects, to reveal their additive contribution in the topographic motion aggravation.     

Effects of rainfall on soil–structure system frequency: Examples based on poroelasticity and a comparison with full-scale measurements

In this paper, a simple two-dimensional soil–structure interaction model, based on Biot's theory of wave propagation in fluid saturated porous media, is used to explain the observed increase of the apparent frequencies of Millikan library in Pasadena, California, during heavy rainfall and recovery within days after the rain. These variations have been measured for small amplitude response (to microtremors and wind excitation), for which Biot's linear theory is valid. The postulated hypothesis is that the observed increases in frequency are due to the water saturation of the soil. The theoretical model used to explore this hypothesis consists of a shear wall supported by a circular foundation embedded in a poroelastic half-space. This rigid foundation model may be appropriate only for the NS response of Millikan library. This paper presents results for the foundation stiffness, and for the system response for model parameters similar to those for Millikan library (located on alluvium with shear wave velocity of about 300 m/s). The foundation impedance matrix, foundation input motion and system response are compared for dry and fully saturated half-space, with permeable and impermeable foundation. The results show that for embedded foundations, the effects of saturation on the horizontal foundation stiffness are as significant as for the vertical stiffness, contrary to what has been known for surface foundations investigated by other authors. Further, the results suggest a 1–2% increase in system frequency of the first two modes of vibration, depending on the drainage condition along the foundation–soil interface. Such increases agree qualitatively with the observations.     

Numerical simulation of liquefaction effects on seismic SSI

The present paper deals with the influence of soil non-linearity, introduced by soil liquefaction, on the soil–foundation–structure interaction phenomena. The objective is to reveal the beneficial or unfavourable effects of the non-linear SSI on both structural drift and settlement of a given structure. Factors such as the signal modification due to liquefaction, and ratios of fundamental frequencies of soil, structure and signal may play an important role on the damage of the structure. The importance of each of these factors is evaluated through a significant parametric study. A 2D coupled finite element modelling is carried out using an elastoplastic multi-mechanism model to represent the soil behaviour. This paper presents the research work we did in the framework of the European Community project NEMISREF (New methods of mitigation of seismic risk on existing foundations, GRDI-40457), to study possible retrofitting measures using GEFDYN computational tools.     

Simplified model for torsional foundation vibrations

A simplified model is presented to simulate unbounded soil for torsional foundation vibration problems. Based on the criterion of equivalent displacement response, a group of equivalent models are developed for a foundation-soil system. An optimal equivalent model is then determined to represent the best simplified model. The parameters of the optimal equivalent model may be obtained by a much easier and more efficient method than lengthy optimization techniques used by most existing lumped-parameter models. The dynamic torsional responses of the foundation-soil system using the optimal equivalent models are very consistent with those obtained by the half-space theory and by the existing models. With fewer parameters, the optimal equivalent model is also found to be as accurate as most existing models. This proposed method may be effectively applied to practical torsionally vibrating problems involving soil–structure interaction.     

Seismic behavior of a post-tensioned integral bridge including soil–structure interaction (SSI)

Current practice usually pays little attention to the effect of soil–structure interaction (SSI) on seismic analysis and design of bridges. The objective of this research study is to assess the significance of SSI on the modal with geometric stiffness and seismic response of a bridge with integral abutments that has been constructed using a new bridge system technology. Emphasis is placed on integral abutment behavior, since abutments together with piers are the most critical elements in securing the integrity of bridge superstructures during earthquakes. Comparison is made between analytical results and field measurements in order to establish the accuracy of the superstructure–abutment model. Sensitivity studies are conducted to investigate the effects of foundation stiffness on the overall dynamic and seismic response of the new bridge system.     

Dynamic response of pile groups with different configurations

A general methodology is outlined for a complete seismic soil—pile-foundation—structure interaction analysis. A Beam-on-Dynamic-Winkler-Foundation (BDWF) simplified model and a Green's-function-based rigorous method are utilized in determining the dynamic response of single piles and pile groups. The simplified model is validated through comparisons with the rigorous method. A comprehensive parameter study is then performed on the effect of pile group configuration on the dynamic impedances of pile foundations. Insight is gained into the nature of dynamic pile—soil—pile interaction. The results presented herein may be used in practice as a guide in obtaining the dynamic stiffness and damping of foundations with a large number of piles.     

Response of L-shaped rigid foundations embedded in a uniform half-space to traveling seismic waves

A simplified indirect boundary element method is applied to compute the impedance functions for L-shaped rigid foundations embedded in a homogeneous viscoelastic half-space. In this method, the waves generated by the 3D vibrating foundation are constructed from radiating sources located on the actual boundary of the foundation. The impedance functions together with the free-field displacements and tractions generated along the soil–foundation interface are used to calculate the foundation input motion for incident P, S and Rayleigh waves. This is accomplished by application of Iguchi's averaging method which, in turn, is verified by comparison with results obtained rigorously using the relation between the solutions of the basic radiation (impedance functions) and scattering (input motions) problems. Numerical results are presented for both surface-supported and embedded foundations. It is shown how the seismic response of L-shaped foundations with symmetrical wings differs from that of enveloping square foundations. The effects of inclination and azimuth of the earthquake excitation are examined as well. These results should be of use in analyses of soil–structure interaction to account for the traveling wave effects usually overlooked in practice.     

Plain strain soil–structure interaction model for a building supported by a circular foundation embedded in a poroelastic half-space

A simple theoretical model for soil–structure interaction in water saturated poroelastic soils is presented, developed to explore if the apparent building–foundation–soil system frequency changes due to water saturation. The model consists of a shear wall supported by a rigid circular foundation embedded in a homogenous, isotropic poroelastic half-space, fully saturated by a compressible and inviscid fluid, and excited by in-plane wave motion. The motion in the soil is governed by Biot's theory of wave propagation in fluid saturated porous media. Helmholtz decomposition and wave function expansion of the two P-wave and the S-wave potentials is used to represent the motion in the soil. The boundary conditions along the contact surface between the soil and the foundation are perfect bond (i.e. welded contact) for the skeleton, and either drained or undrained hydraulic condition for the fluid (i.e. pervious or impervious foundation). For the purpose of this exploratory analysis, the zero stress condition at the free surface is relaxed in the derivation of the foundation stiffness matrix, which enables a closed form solution. The implications of this assumption are discussed, based on published comparisons for the elastic case. Also, a closed form representation is derived for the foundation driving forces for incident plane (fast) P-wave or SV wave. Numerical results and comparison with the full-scale measurements are presented in the companion paper, published in this issue.