Non-linear dynamic earth dam–foundation interaction using a BE–FE method

A general, rigorous, coupled Boundary Element–Finite Element (BE–FE) formulation is presented for non-linear seismic soil–structure interaction in two dimensions. The BE–FE method is applied to investigate the inelastic response of earth dams to transient SV waves. The dam body, consisting of heterogeneous materials modelled with a simple non-linear hysteretic model, is discretized with finite elements, whereas the elastic half-space is discretized with boundary elements. The study focuses on the combined effects of the material non-linearity and foundation flexibility. The results show the significant effect of the foundation flexibility in reducing the response through radiation of energy. For excitations with peak ground accelerations from 0·2gto 0·6g, the crest acceleration amplification ranges from 2·5 to 1·4 and seems to be comparable with field observations and results from other studies. Deamplification increasing with strain is reported at the lower part of the dam. The method is computationally powerful and can be used for efficient non-linear analysis of complex soil–structure systems. The efficiency of the BE–FE method allows further improvements with incorporation of a more advanced constitutive model and consideration of the generation and dissipation of pore-water pressures during the earthquake. © 1998 John Wiley & Sons, Ltd.     

Probabilistic evaluation of soil–foundation–structure interaction effects on seismic structural response

Complex seismic behaviour of soil–foundation–structure (SFS) systems together with uncertainties in system parameters and variability in earthquake ground motions result in a significant debate over the effects of soil–foundation–structure interaction (SFSI) on structural response. The aim of this study is to evaluate the influence of foundation flexibility on the structural seismic response by considering the variability in the system and uncertainties in the ground motion characteristics through comprehensive numerical simulations. An established rheological soil‐shallow foundation–structure model with equivalent linear soil behaviour and nonlinear behaviour of the superstructure has been used. A large number of models incorporating wide range of soil, foundation and structural parameters were generated using a robust Monte‐Carlo simulation. In total, 4.08 million time‐history analyses were performed over the adopted models using an ensemble of 40 earthquake ground motions as seismic input. The results of the analyses are used to rigorously quantify the effects of foundation flexibility on the structural distortion and total displacement of the superstructure through comparisons between the responses of SFS models and corresponding fixed‐base (FB) models. The effects of predominant period of the FB system, linear vs nonlinear modelling of the superstructure, type of nonlinear model used and key system parameters are quantified in terms of different probability levels for SFSI effects to cause an increase in the structural response and the level of amplification of the response in such cases. The results clearly illustrate the risk of underestimating the structural response associated with simplified approaches in which SFSI and nonlinear effects are ignored. Copyright © 2010 John Wiley & Sons, Ltd.     

Bridge–embankment interaction under transverse ground excitation

Seismic performance and dynamic response of bridge–embankments during strong or moderate ground excitations are investigated through finite element (FE) modelling and detailed dynamic analysis. Previous research studies have established that bridge–embankments exhibit increasingly flexible performance under high‐shear deformation levels and that soil displacements at bridge abutment supports may be significant particularly in the transverse direction. The 2D equation of motion is solved for the embankment, in order to evaluate the dynamic characteristics and to describe explicitly the seismic performance and dynamic response under transverse excitations accounting for soil nonlinearities, soil–structure interaction and imposed boundary conditions (BCs). Using the proposed model, equivalent elastic analysis was performed so as to evaluate the dynamic response of approach embankments while accounting for soil–structure interaction. The analytical procedures were applied in the case of a well‐documented bridge with monolithic supports (Painter Street Overcrossing, PSO) which had been instrumented and embankment participation was identified from its response records after the 1971 San Fernando earthquake. The dynamic characteristics and dynamic response of the PSO embankments were evaluated for alternative BCs accounting for soil–structure interaction. Explicit expressions for the evaluation of the critical embankment length L_c are provided in order to quantify soil contribution to the overall bridge system under strong intensity ground excitations. The dynamic response of the entire bridge system (deck–abutments–embankments) was also evaluated through simplified models that considered soil–structure interaction. Results obtained from this analysis are correlated with those of detailed 3D FE models and field data with good agreement. Copyright © 2007 John Wiley & Sons, Ltd.     

Nonlinear earthquake analysis of high arch dam–water–foundation rock systems

A nonlinear finite element model for earthquake response analysis of arch dam–water–foundation rock systems is proposed in this paper. The model includes dynamic dam–water and dam–foundation rock interactions, the opening of contraction joints, the radiation damping of semi‐unbounded foundation rock, the compressibility of impounded water, and the upstream energy propagating along the semi‐unbounded reservoir. Meanwhile, a new equivalent force scheme is suggested to achieve free‐field input in the model. The effects of the earthquake input mechanism, joint opening, water compressibility, and radiation damping on the earthquake response of the Ertan arch dam (240 m high) in China are investigated using the proposed model. The results show that these factors significantly affect the earthquake response of the Ertan arch dam. Such factors should therefore be considered in the earthquake response analysis and earthquake safety evaluation of high arch dams. Copyright © 2011 John Wiley & Sons, Ltd.     

Two-dimensional effective stress analysis of liquefaction of irregular ground including soil–structure interaction

An effective stress method is presented for analysis of seismic response and liquefaction of irregular ground including soil–structure interaction, based on an implicit–explicit finite element method. A pore water pressure is computed with iteration from the total stress considering an undrained condition. The simulated pore water pressure is in reasonably good agreement with the experimental data. The proposed method of analysis is compared with other well-known methods for a one-dimensional model, which is in good agreement. The present effective stress method is also applied to liquefaction problems involving a two-dimensional soil–structure model. The structure is modelled by not only a rigid model but also as a multi-degree-of-freedom system with bi-linear springs. The numerical results are considered to be significant from the viewpoint of earthquake engineering.     

Linear analysis of concrete arch dams including dam–water–foundation rock interaction considering spatially varying ground motions

The available substructure method and computer program for earthquake response analysis of arch dams, including the effects of dam–water–foundation rock interaction and recognizing the semi‐unbounded size of the foundation rock and fluid domains, are extended to consider spatial variations in ground motions around the canyon. The response of Mauvoisin Dam in Switzerland to spatially varying ground motion recorded during a small earthquake is analyzed to illustrate the results from this analysis procedure. Copyright © 2009 John Wiley & Sons, Ltd.     

Influence of frequency‐dependent soil–structure interaction on the fragility of R/C bridges

Bridge performance under earthquake loading can be significantly influenced by the interaction between the structure and the supporting soil. Even though the frequency dependence of the interaction mentioned in this study has long been documented, the simplifying assumption that the dynamic stiffness is dominated by the mean or predominant excitation frequency is still commonly made, primarily as a result of the associated numerical difficulties when the analysis has to be performed in the time domain. This study makes use of the advanced lumped parameter models recently developed 1 in order to quantify the impact of the assumption on the predicted fragility of bridges mentioned in this study. This is achieved by comparing the predicted vulnerability for the case of a reference, well studied, actual bridge using both conventional, frequency‐independent, Kelvin–Voigt models and the aforementioned lumped parameter formulation. Analysis results demonstrate that the more refined consideration of frequency dependence of soil–structure interaction at the piers and the abutments of a bridge not only leads to different probabilities of failure for given intensity measures but also leads to different hierarchy and distribution of damage within the structure for the same set of earthquake ground motions even if the overall probability of exceeding a given damage state is the same. The paper concludes with the comparative assessment of the effect for different soil conditions, foundation configurations, and ground motion characteristics mentioned in this study along with the relevant analysis and design recommendations. Copyright © 2016 John Wiley & Sons, Ltd.     

Seismic soil–structure interaction in multi‐span bridges: Application to a railway bridge

The effects of soil‐structure interaction on the seismic response of multi‐span bridges are investigated by means of a modelling strategy based on the domain decomposition technique. First, the analysis methodology is presented: kinematic interaction analysis is performed in the frequency domain by means of a procedure accounting for radiation damping, soil–pile and pile‐to‐pile interaction; the seismic response of the superstructure is evaluated in the time domain by means of user‐friendly finite element programs introducing suitable lumped parameter models take into account the frequency‐dependent impedances of the soil–foundation system. Second, a real multi‐span railway bridge longitudinally restrained at one abutment is analyzed. The input motion is represented by two sets of real accelerograms: one consistent with the Italian seismic code and the other constituted by five records characterized by different frequency contents. The seismic response of the compliant‐base model is compared with that obtained from a fixed‐base model. Pile stress resultants due to kinematic and inertial interactions are also evaluated. The application demonstrates the importance of performing a comprehensive analysis of the soil–foundation–structure system in the design process, in order to capture the effects of soil‐structure interaction in each structural element that may be beneficial or detrimental. Copyright © 2011 John Wiley & Sons, Ltd.     

Seismic effectiveness of tuned mass dampers considering soil–structure interaction

This paper presents how soil–structure interaction affects the seismic performance of Tuned Mass Dampers (TMD) when installed on flexibly based structures. Previous studies on this subject have led to inconsistent conclusions since the soil and structure models employed considerably differ from each other. A generic frequency-independent model is used in this paper to represent a general soil–structure system, whose parameters cover a wide spectrum of soil and structural characteristics. The model structure is subjected to a stationary random excitation and the root-mean-square responses of engineering interest are used to measure the TMD's performance. Extensive parametric studies have shown that strong soil–structure interaction significantly defeats the seismic effectiveness of TMD systems. As the soil shear wave velocity decreases, TMD systems become less effective in reducing the maximum response of structures. For a structure resting on soft soil, the TMD system can hardly reduce the structural seismic response due to the high damping characteristics of soil–structure systems. The model structure is further subjected to the NS component of the 1940 El Centro, California earthquake to confirm the TMD's performance in a more realistic environment. Copyright © 1999 John Wiley & Sons Ltd.     

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     

Input and system identification of the Hualien soil–structure interaction system using earthquake response data

This paper presents an input and system identification technique for a soil–structure interaction system using earthquake response data. Identification is carried out on the Hualien large‐scale seismic test structure, which was built in Taiwan for international joint research. The identified quantities are the input ground acceleration as well as the shear wave velocities of the near‐field soil regions and Young's moduli of the shell sections of the structure. The earthquake response analysis on the soil–structure interaction system is carried out using the finite element method incorporating the infinite element formulation for the unbounded layered soil medium and the substructured wave input technique. The criterion function for the parameter estimation is constructed using the frequency response amplitude ratios of the earthquake responses measured at several points of the structure, so that the information on the input motion may be excluded. The constrained steepest descent method is employed to obtain the revised parameters. The simulated earthquake responses using the identified parameters and input ground motion show excellent agreement with the measured responses. Copyright © 2003 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.     

Soil–pile–structure interaction: experimental outcomes from shaking table tests

An effective way to study the complex seismic soil‐structure interaction phenomena is to investigate the response of physical scaled models in 1‐g or n‐g laboratory devices. The outcomes of an extensive experimental campaign carried out on scaled models by means of the shaking table of the Bristol Laboratory for Advanced Dynamics Engineering, University of Bristol, UK, are discussed in the present paper. The experimental model comprises an oscillator connected to a single or a group of piles embedded in a bi‐layer deposit. Different pile head conditions, that is free head and fixed head, several dynamic properties of the structure, including different masses at the top of the single degree of freedom system, excited by various input motions, e.g. white noise, sinedwells and natural earthquake strong motions recorded in Italy, have been tested. In the present work, the modal dynamic response of the soil–pile–structure system is assessed in terms of period elongation and system damping ratio. Furthermore, the effects of oscillator mass and pile head conditions on soil–pile response have been highlighted, when the harmonic input motions are considered. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Kinematic soil–structure interaction for long-span cable-supported bridges

A general procedure is presented to study the dynamic soil–structure interaction effects on the response of long-span suspension and cable-stayed bridges subjected to spatially varying ground motion at the supporting foundations. The foundation system is represented by multiple embedded cassion foundations and the frequency-dependent impedance matrix for the multiple foundations system takes into account also the cross-interaction among adjacent foundations through the soil. To illustrate the potential implementation of the analysis, a numerical example is presented in which the dynamic response of the Vincent–Thomas suspension bridge (Los Angeles, CA) subjected to the 1987 Whittier earthquake is investigated. Although both kinematic and inertial effects are included in the general procedure, only the kinematic effects of the soil–structure interaction are considered in the analysis of the test case. The results show the importance of the kinematic soil–foundation interaction on the structural response. These effects are related to the type, i.e. SH-, SV-, P- or Rayleigh waves and to the inclination of the seismic wave excitation. Moreover, rocking components of the foundation motion are emphasized by the embedment of the foundation system and greatly alter the structural response.     

Investigation of the dynamic behaviour of a storage tank with different foundation types focusing on the soil‐foundation‐structure interactions using centrifuge model tests

This paper proposes a dynamic centrifuge model test method for the accurate simulation of the behaviours of a liquid storage tank with different types of foundations during earthquakes. The method can be used to determine the actual stress conditions of a prototype storage‐tank structure. It was used in the present study to investigate the soil‐foundation‐structure interactions of a simplified storage tank under two different earthquake motions, which were simulated using a shaking table installed in a centrifuge basket. Three different types of foundations were considered, namely, a shallow foundation, a slab on the surface of the ground connected to piles and a slab with disconnected piles. The test results were organised to compare the ground surface and foundation motions, the slab of foundation and top of structure motions and the horizontal and vertical motions of the slab, respectively. These were used to establish the complex dynamic behaviours of tank models with different foundations. The effects of soil–foundation–structure interaction with three foundation conditions and two different earthquake motions are focused and some important factors, that should be considered for future designs are also discussed in this research. Copyright © 2017 John Wiley & Sons, Ltd.     

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.     

Probabilistic estimation of inelastic displacement demands in soil–structure interaction systems on cohesive soils

Inelastic displacement ratios (IDRs) of nonlinear soil–structure interaction (SSI) systems located at sites with cohesive soils are investigated in this study. To capture the effects of inelastic cyclic behavior of the supporting soil, the Beam on Nonlinear Winkler Foundation (BNWF) model is used. The superstructure is modeled using an inelastic single-degree-of-freedom (SDOF) system model. Nonlinear SSI systems representing various combinations of unconfined compressive strengths and shear wave velocities are considered in the analysis. A set of strong ground motions recorded at sites with soft to stiff soils is used for considering the record-to-record variability of IDRs. It is observed that IDRs for nonlinear SSI systems are sensitive to the strength and the stiffness properties of both the soil and the structure. For the case of SSI systems on the top of cohesive soils, the compressive strength of the soil has a significant impact on the IDRs, which cannot be captured by considering only the shear wave velocity of the soil. Based on the results of nonlinear time-history analysis, a new equation is proposed for estimating the mean and the dispersion of IDRs of SSI systems depending on the characteristic properties of the supporting soil, dimensions of the foundation, and properties of the superstructure. A probabilistic framework is presented for the performance-based seismic design of SSI systems located at sites with cohesive soils.     

Soil–pile–structure interaction under SH wave excitation

A continuum model for the interaction analysis of a fully coupled soil–pile–structure system under seismic excitation is presented in this paper. Only horizontal shaking induced by harmonic SH waves is considered so that the soil–pile–structure system is under anti‐plane deformation. The soil mass, pile and superstructure were all considered as elastic with hysteretic damping, while geometrically both pile and structures were simplified as a beam model. Buildings of various heights in Hong Kong designed to resist wind load were analysed using the present model. It was discovered that the acceleration of the piled‐structures at ground level can, in general, be larger than that of a free‐field shaking of the soil site, depending on the excitation frequency. For typical piled‐structures in Hong Kong, the amplification factor of shaking at the ground level does not show simple trends with the number of storeys of the superstructure, the thickness and the stiffness of soil, and the stiffness of the superstructure if number of storeys is fixed. The effect of pile stiffness on the amplification factor of shaking is, however, insignificant. Thus, simply increasing the pile size or the superstructure stiffness does not necessarily improve the seismic resistance of the soil–pile–structure system; on the contrary, it may lead to excessive amplification of shaking for the whole system. Copyright © 2003 John Wiley & Sons, Ltd.     

H-adaptivity applied to liquefiable soil in nonlinear analysis of soil–pile interaction

This paper is concerned with application of the h-adaptive finite element method to dynamic analysis of a pile in liquefiable soil considering large deformation. In finite element analysis of pile behavior in liquefiable soil during an earthquake, especially considering large deformation of liquefied ground, error due to discretization in the zone near the pile becomes very large. Our purpose was to refine the approximation of the finite element method. The updated Lagrangian formulation and a cyclic elasto-plastic model based on the kinematic hardening rule were adopted to deal with the nonlinearity of the soil. The mixed finite element and finite difference methods together with the u-p formulation and Biot's two-phase mixture theory were used. To improve the accuracy and increase the efficiency of finite element analysis, an h-adaptive scheme that included a posteriori error estimation and h-version mesh refinement was applied to the analysis. The calculated results of effective stress were smoothed locally by the extrapolation method and smoothed stress was used to calculate the L₂ norm of the effective stress error in the last step of the calculation of each time increment. The mesh was refined by a fission procedure based on the indication of the error estimate As a numerical example, a soil–pile interaction system loaded cyclically was analyzed by our method.