Three-dimensional dynamic soil–structure interaction analysis in the time domain

A new numerical procedure is proposed for the analysis of three-dimensional dynamic soil–structure interaction in the time domain. In this study, the soil is modelled as a linear elastic solid, however, the methods developed can be adapted to include the effects of soil non-linearities and hysteretic damping in the soil. A substructure method, in which the unbounded soil is modelled by the scaled boundary finite-element method, is used and the structure is modelled by 8–21 variable-number-node three-dimensional isoparametric or subparametric hexahedral curvilinear elements. Approximations in both time and space, which lead to efficient schemes for calculation of the acceleration unit-impulse response matrix, are proposed for the scaled boundary finite-element method resulting in significant reduction in computational effort with little loss of accuracy. The approximations also lead to a very efficient scheme for evaluation of convolution integrals in the calculation of soil–structure interaction forces. The approximations proposed in this paper are also applicable to the boundary element method. These approximations result in an improvement over current methods. A three-dimensional Dynamic Soil–Structure Interaction Analysis program (DSSIA-3D) is developed, and seismic excitations (S-waves, P-waves, and surface waves) and externally applied transient loadings can be considered in analysis. The computer program developed can be used in the analysis of three-dimensional dynamic soil–structure interaction as well as in the analysis of wave scattering and diffraction by three-dimensional surface irregularities. The scattering and diffraction of seismic waves (P-, S-, and Rayleigh waves) by various three-dimensional surface irregularities are studied in detail, and the numerical results obtained are in good agreement with those given by other authors. Numerical studies show that the new procedure is suitable and very efficient for problems which involve low frequencies of interest for earthquake engineering. Copyright © 1999 John Wiley & Sons Ltd     

A discontinuous Galerkin method for seismic soil–structure interaction analysis in the time domain

A technique for modeling transient wave propagation in unbounded media is extended and applied to seismic soil–structure interaction analysis in the time domain. The technique, based on the discontinuous Galerkin method, requires lower computational cost and less storage than the boundary element method, and the time‐stepping scheme resulting from Newmark's method in conjunction with the technique is unconditionally stable, allowing for efficient and robust time‐domain computations. To extend the technique to cases characterized by seismic excitation, the free‐field motion is used to compute effective forces, which are introduced on the boundary of the computational domain containing the structure and the soil in the vicinity of the structure. A numerical example on a dam–foundation system subjected to seismic excitation demonstrates the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

A hybrid time frequency domain formulation for non-linear soil–structure interaction

Methods that combine frequency and time domain techniques offer an attractive alternative for solving Soil–Structure-interaction problems where the structure exhibits non-linear behaviour. In the hybrid-frequency-time-domain procedure a reference linear system is solved in the frequency domain and the difference between the actual restoring forces and those in the linear model are treated as pseudo-forces. In the solution scheme explored in this paper, designated as the hybrid-time-frequency-domain (HTFD) procedure, the equations of motion are solved in the time domain with due consideration for non-linearities and with the unbounded medium represented by frequency-independent springs and dampers. The frequency dependency of the impedance coefficients is introduced by means of pseudo-forces evaluated in the frequency domain at the end of each iteration. A criterion of stability for the HTFD approach is derived analytically and its validity is sustained numerically. As is often the case, the criterion takes the form of a limit of unity on the spectral radius of an appropriately defined matrix. Inspection of the terms in this matrix shows that convergence can be guaranteed by suitable selection of the reference impedance. The CPU times required to obtain converged solutions with the HTFD are found, in a number of numerical simulations, to be up to one order of magnitude less than those required by the alternative hybrid-frequency-time-domain approach. © 1998 John Wiley & Sons, Ltd.     

Discrete model for dynamic through-the-soil coupling of 3-D foundations and structures

An efficient discrete model for predicting the dynamic through-the-soil interaction between adjacent rigid, surface foundations supported by a homogeneous, isotropic and linear elastic half-space is presented. The model utilizes frequency-independent springs and dashpots, and the foundation mass, for the consideration of soil–foundation interaction. The through-the-soil coupling of the foundations is attained by frequency-independent stiffness and damping functions, developed in this work, that interconnect the degrees of freedom of the entire system of foundations. The dynamic analysis of the resulting coupled system is performed in the time domain and includes the time lagging effects of coupled dynamic input due to wave propagation using an appropriate modification of the Wilson-θ method. The basic foundation interaction model is also extended to the evaluation of coupled building-foundation systems. © 1998 John Wiley & Sons, Ltd.     

System identification for evaluating soil–structure interaction effects in buildings from strong motion recordings

Parametric system identification is used to evaluate seismic soil–structure interaction effects in buildings. The input–output strong motion data pairs needed for evaluations of flexible- and fixed-base fundamental mode parameters are derived. Recordings of lateral free-field, foundation, and roof motions, as well as foundation rocking, are found to be necessary for direct evaluations of modal parameters for both cases of base fixity. For the common situation of missing free-field or base rocking motions, procedures are developed for estimating the modal parameters that cannot be directly evaluated. The accuracy of these estimation procedures for fundamental mode vibration period and damping is verified for eleven sites with complete instrumentation of the structure, foundation, and free-field. © 1998 John Wiley & Sons, Ltd.     

Seismic response and damage analysis of buildings supported on flexible soils

The investigation reported in this paper studies the effects of soil–structure interaction (SSI) on the seismic response and damage of building–foundation systems. A simple structural model is used for conducting a parametric study using a typical record obtained in the soft soil area of Mexico City during the 1985 earthquake. Peak response parameters chosen for this study were the roof displacement relative to the base and the hysteretic energy dissipated by the simple structural model. A damage parameter is also evaluated for investigating the SSI effects on the seismic damage of buildings. The results indicate that in most cases of inelastic response, SSI effects can be evaluated considering the rigid‐base case and the SSI period. Copyright © 2000 John Wiley & Sons, Ltd.     

Local and global shape functions in a boundary element formulation for the calculation of traffic induced vibrations

This paper compares the use of local and global shape functions in a boundary element method that is used in a prediction model for traffic induced vibrations. The boundary element formulation describes the interaction problem between a linear elastic layered half-space and a longitudinally invariant structure representing a road or a railway track. The boundary element formulation in the frequency–wavenumber domain is obtained by means of a weighted residual method. Constant element shape functions, as well as Legendre and Chebyshev shape functions are considered. Their effect on both accuracy and computational effort is investigated. The presence of a singularity in the Chebyshev based shape functions allows to obtain a better approximation for the soil tractions. The theory is applied to road traffic induced vibrations where the response is calculated in a large number of output points.     

A practical method to transform frequency dependent impedance to time domain

In order to perform time history earthquake response analyses with consideration to both the dynamic soil–structure interaction and the non‐linear behaviour of the structure, it is important to transform the soil impedance in the frequency domain to the impulse response in the time domain. In this paper, a new transform method with high practicality is proposed. First, the formulation of the proposed transform method is described. Next, the validity of the method is examined using an example problem whose impulse response is analytically obtained. Then, the impedance of the rigid foundation on 2‐layered soil is transformed to the time domain, and the characteristics of the impulse response are investigated. Finally, time history earthquake response analyses of a structure on the soil using the obtained impulse response are carried out. The validity and the efficiency of the proposed method are confirmed through these investigations. Copyright © 2005 John Wiley & Sons, Ltd.     

Recursive evaluation of interaction forces and property matrices from unit-impulse response functions of unbounded medium based on balancing approximation

Starting from the unit-impulse response matrix of the unbounded medium, a discrete-time formulation permitting the recursive evaluation of the interaction forces and a continuous-time formulation yielding property matrices corresponding to a model with a finite number of degrees of freedom are discussed. This is achieved using the balancing approximation method which is easily automated, guarantees stability and leads to highly accurate results. © 1998 John Wiley & Sons Ltd.     

Seismic design of yielding structures on flexible foundations

This paper introduces a simple method to consider the effects of inertial soil–structure interaction (SSI) on the seismic demands of a yielding single‐degree‐of‐freedom structure. This involves idealizing the yielding soil–structure system as an effective substitute oscillator having a modified period, damping ratio, and ductility. A parametric study is conducted to obtain the ratio between the displacement ductility demand of a flexible‐base system and that of the corresponding fixed‐base system. It is shown that while additional foundation damping can reduce the overall response, the effects of SSI may also increase the ductility demand of some structures, mostly being ductile and having large structural aspect ratio, up to 15%. Finally, a design procedure is provided for incorporation of the SSI effects on structural response. Copyright © 2015 John Wiley & Sons, Ltd.     

The role of soil in the collapse of 18 piers of Hanshin Expressway in the Kobe earthquake

An investigation is presented of the collapse of a 630 m segment (Fukae section) of the elevated Hanshin Expressway during the 1995 Kobe earthquake. The earthquake has, from a geotechnical viewpoint, been associated with extensive liquefactions, lateral soil spreading, and damage to waterfront structures. Evidence is presented that soil–structure interaction (SSI) in non‐liquefied ground played a detrimental role in the seismic performance of this major structure. The bridge consisted of single circular concrete piers monolithically connected to a concrete deck, founded on groups of 17 piles in layers of loose to dense sands and moderate to stiff clays. There were 18 spans in total, all of which suffered a spectacular pier failure and transverse overturning. Several factors associated with poor structural design have already been identified. The scope of this work is to extend the previous studies by investigating the role of soil in the collapse. The following issues are examined: (1) seismological and geotechnical information pertaining to the site; (2) free‐field soil response; (3) response of foundation‐superstructure system; (4) evaluation of results against earlier studies that did not consider SSI. Results indicate that the role of soil in the collapse was multiple: First, it modified the bedrock motion so that the frequency content of the resulting surface motion became disadvantageous for the particular structure. Second, the compliance of soil and foundation altered the vibrational characteristics of the bridge and moved it to a region of stronger response. Third, the compliance of the foundation increased the participation of the fundamental mode of the structure, inducing stronger response. It is shown that the increase in inelastic seismic demand in the piers may have exceeded 100% in comparison with piers fixed at the base. These conclusions contradict a widespread view of an always‐beneficial role of seismic SSI. Copyright © 2005 John Wiley & Sons, Ltd.     

Seismic response analysis of highway overcrossings including soil–structure interaction

This paper presents a systematic procedure for the seismic response analysis of highway overcrossings. The study employs an elementary stick model and a more sophisticated finite element formulation to compute response quantities. All dynamic stiffnesses of approach embankments and pile groups are approximated with frequency‐independent springs and dashpots that have been established elsewhere. A real eigenvalue analysis confirms the one‐to‐one correspondence between modal characteristics obtained with the three‐dimensional finite element solutions and the result of the simpler stick‐model idealization. A complex eigenvalue analysis yields modal damping values in the first six modes of interest and shows that modal damping ratios assume values much higher than those used by Caltrans. The validity of the proposed method is examined by comparing the computed time response quantities with records from the Meloland Road and Painter Street overcrossings located in southern and northern California, respectively. The proposed procedure allows for inexpensive parametric analysis that examines the importance of considering soil–structure interaction at the end abutments and centre bent. Results and recommendations presented by past investigations are revisited and integrated in comprehensive tables that improve our understanding of the dynamic characteristics and behaviour of freeway overcrossings. The study concludes with a step‐by‐step methodology that allows for a simple, yet dependable dynamic analysis of freeway overcrossings, that involves a stick model and frequency‐independent springs and dashpots. Copyright © 2002 John Wiley & Sons, Ltd.     

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

Influence of dynamic soil–structure interaction on the nonlinear response and seismic reliability of multistorey systems

A set of reinforced concrete structures with gravitational loads and mechanical properties (strength and stiffness) representative of systems designed for earthquake resistance in accordance with current criteria and methods is selected to study the influence of dynamic soil–structure interaction on seismic response, ductility demands and reliability levels. The buildings are considered located at soft soil sites in the Valley of Mexico and subjected to ground motion time histories simulated in accordance with characteristic parameters of the maximum probable earthquake likely to occur during the system's expected life. For the near‐resonance condition the effects of soil–structure interaction on the ductility demands depend mainly on radiation damping. According to the geometry of the structures studied this damping is strongly correlated with the aspect ratio, obtained by dividing the building height by its width. In this way, for structures with aspect ratio greater than 1.4 the storey and global ductility demands increase with respect to those obtained with the same structures but on rigid base, while for structures with aspect ratio less than 1.4 the ductility demands decrease with respect to those for the structures on rigid base. For the cases when the fundamental period of the structure has values very different from the dominant ground period, soil–structure interaction leads in all cases to a reduction of the ductility demands, independently of the aspect ratio. The reliability index β is obtained as a function of the base shear ratio and of the seismic intensity acting on the nonlinear systems subjected to the simulated motions. The resulting reliability functions are very similar for systems on rigid or on flexible foundation, provided that in the latter case the base rotation and the lateral displacement are removed from the total response of the system. Copyright © 2006 John Wiley & Sons, Ltd.     

Amplitude dependence of equivalent modal parameters in monitored buildings during earthquake swarms

This paper investigates the dynamic response of three sample buildings belonging to the Seismic Observatory for Structures, the Italian network for the permanent seismic monitoring of strategic structures, managed by the Italian Department of Civil Protection. The case studies cover different building types that could loosely represent the Italian building stock, with a special emphasis on cultural heritage and masonry structures. Observed under a low‐intensity seismic swarm comprising about 30 aftershocks after a main event, the three buildings are analysed through an input–output, model‐driven linear dynamic identification procedure, depicting the relation between the shaking level at the site and the variation of the equivalent structural modal parameters, while keeping into account the effects of soil–structure interaction. Finite element models will be used to investigate one of the case studies and to compare the law of variation of the structural modal parameters with respect to simplified models proposed by technical standards. Copyright © 2017 John Wiley & Sons, Ltd.     

Free‐vibration analysis of a three‐dimensional soil–structure system

A procedure which involves a non‐linear eigenvalue problem and is based on the substructure method is proposed for the free‐vibration analysis of a soil–structure system. In this procedure, the structure is modelled by the standard finite element method, while the unbounded soil is modelled by the scaled boundary finite element method. The fundamental frequency, and the corresponding radiation damping ratio as well as the modal shape are obtained by using inverse iteration. The free vibration of a dam–foundation system, a hemispherical cavity and a hemispherical deposit are analysed in detail. The numerical results are compared with available results and are also verified by the Fourier transform of the impulsive response calculated in the time domain by the three‐dimensional soil–structure–wave interaction analysis procedure proposed in our previous paper. The fundamental frequency obtained by the present procedure is very close to that obtained by Touhei and Ohmachi, but the damping ratio and the imaginary part of modal shape are significantly different due to the different definition of damping ratio. This study shows that although the classical mode‐superposition method is not applicable to a soil–structure system due to the frequency dependence of the radiation damping, it is still of interest in earthquake engineering to evaluate the fundamental frequency and the corresponding radiation damping ratio of the soil–structure system. Copyright © 2001 John Wiley & Sons, Ltd.     

Effects of foundation embedment during building–soil interaction

A numerical solution for evaluating the effects of foundation embedment on the effective period and damping and the response of soil–structure systems is presented. A simple system similar to that used in practice to account for inertial interaction effects is investigated, with the inclusion of kinematic interaction effects for the important special case of vertically incident shear waves. The effective period and damping are obtained by establishing an equivalence between the interacting system excited by the foundation input motion and a replacement oscillator excited by the free-field ground motion. In this way, the use of standard free-field response spectra applicable to the effective period and damping of the system is permitted. Also, an approximate solution for total soil–structure interaction is presented, which indicates that the system period is insensitive to kinematic interaction and the system damping may be expressed as that for inertial interaction but modified by a factor due to kinematic interaction. Results involving both kinematic and inertial effects are compared with those obtained for no soil–structure interaction and inertial interaction only. The more important parameters involved are identified and their influences are examined over practical ranges of interest. © 1998 John Wiley & Sons, Ltd.     

The influence of dynamic soil–structure interaction on traffic induced vibrations in buildings

This paper deals with the dynamic response of buildings due to traffic induced wave fields. The response of a two-storey single family dwelling due to the passage of a two-axle truck on a traffic plateau is computed with a model that fully accounts for the dynamic interaction between the soil and the structure. The results of three cases where the structure is founded on a slab foundation, a strip foundation and a box foundation are calculated and a comprehensive analysis of the dynamic structural response is performed. A methodology is also proposed to calculate the structural response, neglecting the effects of dynamic soil–structure interaction. A comparison with the results of calculations where dynamic soil–structure interaction is accounted for shows that a good approximation is obtained in the case of a rigid structure resting on a soft soil.     

Effects of pulse period of near‐field ground motions on the seismic demands of soil–MDOF structure systems using mathematical pulse models

In this paper, the effects of pulse period associated with near‐field ground motions on the seismic demands of soil–MDOF structure systems are investigated by using mathematical pulse models. Three non‐dimensional parameters are employed as the crucial parameters, which govern the responses of soil–structure systems: (1) non‐dimensional frequency as the structure‐to‐soil stiffness ratio; (2) aspect ratio of the superstructure; and (3) structural target ductility ratio. The soil beneath the superstructure is simulated on the basis of the Cone model concept. The superstructure is modeled as a nonlinear shear building. Interstory drift ratio is selected as the main engineering demand parameter for soil–structure systems. It is demonstrated that the contribution of higher modes to the response of soil–structure system depends on the pulse‐to‐interacting system period ratio instead of pulse‐to‐fixed‐base structure period ratio. Furthermore, results of the MDOF superstructures demonstrate that increasing structural target ductility ratio results in the first‐mode domination for both fixed‐base structure and soil–structure system. Additionally, increasing non‐dimensional frequency and aspect ratio of the superstructure respectively decrease and increase the structural responses. Moreover, comparison of the equivalent soil–SDOF structure system and the soil–MDOF structure system elucidates that higher‐mode effects are more significant, when soil–structure interaction is taken into account. In general, the effects of fling step and forward directivity pulses on activating higher modes of the superstructure are more sever in soil–structure systems, and in addition, the influences of forward directivity pulses are more considerable than fling step ones. Copyright © 2013 John Wiley & Sons, Ltd.     

Stability of the time‐domain analysis method including a frequency‐dependent soil–foundation system

A number of methods have been proposed that utilize the time‐domain transformations of frequency‐dependent dynamic impedance functions to perform a time‐history analysis. Though these methods have been available in literature for a number of years, the methods exhibit stability issues depending on how the model parameters are calibrated. In this study, a novel method is proposed with which the stability of a numerical integration scheme combined with time‐domain representation of a frequency‐dependent dynamic impedance function can be evaluated. The method is verified with three independent recursive parameter models. The proposed method is expected to be a useful tool in evaluating the potential stability issue of a time‐domain analysis before running a full‐fledged nonlinear time‐domain analysis of a soil–structure system in which the dynamic impedance of a soil–foundation system is represented with a recursive parameter model. Copyright © 2015 John Wiley & Sons, Ltd.