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.     

A lumped parameter model for time‐domain inertial soil‐structure interaction analysis of structures on pile foundations

The paper presents a lumped parameter model for the approximation of the frequency‐dependent dynamic stiffness of pile group foundations. The model can be implemented in commercial software to perform linear or nonlinear dynamic analyses of structures founded on piles taking into account the frequency‐dependent coupled roto‐translational, vertical, and torsional behaviour of the soil‐foundation system. Closed‐form formulas for estimating parameters of the model are proposed with reference to pile groups embedded in homogeneous soil deposits. These are calibrated with a nonlinear least square procedure, based on data provided by an extensive non‐dimensional parametric analysis performed with a model previously developed by the authors. Pile groups with square layout and different number of piles embedded in soft and stiff soils are considered. Formulas are overall well capable to reproduce parameters of the proposed lumped system that can be straightforwardly incorporated into inertial structural analyses to account for the dynamic behaviour of the soil‐foundation system. Some applications on typical bridge piers are finally presented to show examples of practical use of the proposed model. Results demonstrate the capability of the proposed lumped system as well as the formulas efficiency in approximating impedances of pile groups and the relevant effect on the response of the superstructure.     

A time‐domain seismic SSI analysis method for inelastic bridge structures through the use of a frequency‐dependent lumped parameter model

It is highlighted in the past that the soil–structure interaction phenomenon can produce a significant alteration on the response of a bridge structure. A variety of approaches has been developed in the past, which is capable of tackling the soil–structure interaction problem from different perspectives. The popular approach of a discretized truncated finite element model of the soil domain is not always a numerically viable solution, especially for computationally demanding simulations such as the probabilistic fragility analysis of a bridge structure or the real time hybrid simulation. This paper aims to develop a complete modeling procedure that is capable of coping with the soil–structure interaction problem of inelastic bridge structures through the use of a frequency dependent lumped parameter assembly. The proposed procedure encounters accuracy and global stability issues observed on past methods while maintaining the broad applicability of the method by any commercial FEM software. A case study of an overpass bridge structure under earthquake excitations is illustrated in order to verify the proposed method. Copyright © 2015 John Wiley & Sons, Ltd.     

Effects of seismic isolation on the seismic response of a California high‐speed rail prototype bridge with soil‐structure and track‐structure interactions

With the launch of the high‐speed train project in California, the seismic risk is a crucial concern to the stakeholders. To investigate the seismic behavior of future California High‐Speed Rail (CHSR) bridge structures, a 3D nonlinear finite‐element model of a CHSR prototype bridge is developed. Soil‐structure and track‐structure interactions are accounted for in this comprehensive numerical model used to simulate the seismic response of the bridge and track system. This paper focuses on examining potential benefits and possible drawbacks of the a priori promising application of seismic isolation in CHSR bridges. Nonlinear time history analyses are performed for this prototype bridge subjected to two bidirectional horizontal historical earthquake ground motions each scaled to two different seismic hazard levels. The effect of seismic isolation on the seismic performance of the bridge is investigated through a detailed comparison of the seismic response of the bridge with and without seismic isolation. It is found that seismic isolation significantly reduces the deck acceleration and the force demand in the bridge substructure (i.e., piers and foundations), especially for high‐intensity earthquakes. However, seismic isolation increases the deck displacement (relative to the pile cap) and the stresses in the rails. These findings imply that seismic isolation can be promisingly applied to CHSR bridges with due consideration of balancing its beneficial and detrimental effects through using appropriate isolators design. The optimum seismic isolator properties can be sought by solving a performance‐based optimum seismic design problem using the nonlinear finite‐element model presented herein. Copyright © 2016 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.     

Response of nonlinear soil‐MDOF structure systems subjected to distinct frequency‐content components of near‐fault ground motions

This paper is devoted to investigate the effects of near‐fault ground motions on the seismic responses of nonlinear MDOF structures considering soil‐structure interaction (SSI). Attempts are made to take into account the effects of different frequency‐content components of near‐fault records including pulse‐type (PT) and high‐frequency (HF) components via adopting an ensemble of 54 near‐fault ground motions. A deep sensitivity analysis is implemented based on the main parameters of the soil‐structure system. The soil is simulated based on the Cone model concept, and the superstructure is idealized as a nonlinear shear building. The results elucidate that SSI has approximately increasing and mitigating effects on structural responses to the PT and HF components, respectively. Also, a threshold period exists above which the HF component governs the structural responses. As the fundamental period of the structure becomes shorter and structural target ductility reduces, the contribution of the HF component to the structural responses increases, elaborately. Soil flexibility makes the threshold period increase, and the effect of the PT component becomes more significant than the HF one. In the case of soil‐structure system, slenderizing the structure also increases this threshold period and causes the PT component to be dominant. Copyright © 2013 John Wiley & Sons, Ltd.     

Collapse capacity of soil‐structure systems under pulse‐like earthquakes

In this paper, a comprehensive study is carried out to examine the possibility of dynamic instability produced in soil‐structure systems using an ensemble of 50 pulse‐like records. A number of structural models with various vibration periods varying from 0.1 to 2 s are used in this study. The superstructure is simulated as a non‐linear SDOF oscillator with a two‐segment backbone curve having negative post‐yield stiffness. The soil is idealized based on the cone model concept widely used for practical purposes. The results of this investigation demonstrate that as the pulse period increases, the collapse relative lateral strength ratio decreases and probability of dynamic instability enhances. Moreover, soil flexibility makes the system dynamically more unstable, and as the non‐dimensional frequency increases, the collapse relative lateral strength ratio highly reduces. Additionally, the aspect ratio has insignificant effects on the collapse relative lateral strength ratio. Furthermore, comparison of the collapse relative lateral strength ratios resulting from pulse‐like motions with those obtained from studies under non‐pulse‐like motions (Miranda and Akkar; FEMA 440) for fixed‐base conditions shows that high‐velocity pulses exacerbate the dynamic instability problem and decrease the collapse relative lateral strength ratio. Copyright © 2014 John Wiley & Sons, Ltd.     

Probabilistic analysis of soil‐structure interaction effects on the seismic performance of structures

This paper revisits the phenomenon of dynamic soil‐structure interaction (SSI) with a probabilistic approach. For this purpose, a twofold objective is pursued. First, the effect of SSI on inelastic response of the structure is studied considering the prevailing uncertainties. Second, the consequence of practicing SSI provisions of the current seismic design codes on the structural performance is investigated in a probabilistic framework. The soil‐structure system is modeled by the sub‐structure method. The uncertainty in the properties of the soil and the structure is described by random variables that are input to this model. Monte Carlo sampling analysis is employed to compute the probability distribution of the ductility demand of the structure, which is selected as the metrics for the structural performance. In each sample, a randomly generated soil‐structure system is subjected to a randomly selected and scaled ground motion. To comprehensively model the uncertainty in the ground motion, a suite of 3269 records is employed. An extensive parametric study is conducted to cover a wide range of soil‐structure systems. The results reveal the probability that SSI increases the ductility demand of structures designed based on the conventional fixed‐based assumption but built on flexible soil in reality. The results also show it is highly probable that practicing SSI provisions of modern seismic codes increase the ductility demand of the structure. Copyright © 2016 John Wiley & Sons, Ltd.     

Earthquake soil‐structure interaction of nuclear power plants,differences in response to 3‐D, 3 × 1‐D,and 1‐D excitations

In soil‐structure interaction modeling of systems subjected to earthquake motions, it is classically assumed that the incoming wave field, produced by an earthquake, is unidimensional and vertically propagating. This work explores the validity of this assumption by performing earthquake soil‐structure interaction modeling, including explicit modeling of sources, seismic wave propagation, site, and structure. The domain reduction method is used to couple seismic (near‐field) simulations with local soil‐structure interaction response. The response of a generic nuclear power plant model computed using full earthquake soil‐structure interaction simulations is compared with the current state‐of‐the‐art method of deconvolving in depth the (simulated) free‐field motions, recorded at the site of interest, and assuming that the earthquake wave field is spatially unidimensional. Results show that the 1‐D wave‐field assumption does not hold in general. It is shown that the way in which full 3‐D analysis results differ from those which assume a 1‐D wave field is dependent on fault‐to‐site geometry and motion frequency content. It is argued that this is especially important for certain classes of soil‐structure systems of which nuclear power plants subjected to near‐field earthquakes are an example.     

A nonlinear macroelement model for the seismic analysis of masonry buildings

The macroelement technique for modelling the nonlinear response of masonry panels is particularly efficient and suitable for the analysis of the seismic behaviour of complex walls and buildings. The paper presents a macroelement model specifically developed for simulating the cyclic in‐plane response of masonry walls, with possible applications in nonlinear static and dynamic analysis of masonry structures. The model, starting from a previously developed macroelement model, has been refined in the representation of flexural–rocking and shear damage modes, and it is capable of fairly simulating the experimental response of cyclic tests performed on masonry piers. By means of two internal degrees of freedom, the two‐node macroelement permits to represent the coupling of axial and flexural response as well as the interaction of shear and flexural damage. Copyright © 2013 John Wiley & Sons, Ltd.     

A macro‐element model for non‐linear soil–shallow foundation–structure interaction under seismic loads: theoretical development and experimental validation on large scale tests

In this paper, different formulations of a macro‐element model for non‐linear dynamic soil‐structure interaction analyses of structures lying on shallow foundations are first reviewed, and secondly, a novel formulation is introduced, which combines some of the characteristics of previous approaches with several additional features. This macro‐element allows one to model soil‐footing geometric (uplift) and material (soil plasticity) non‐linearities that are coupled through a stiffness degradation model. Footing uplift is introduced by a simple non‐linear elastic model based on the concept of effective foundation width, whereas soil plasticity is treated by means of a bounding surface approach in which a vertical load mapping rule is implemented. This mapping is particularly suited for the seismic loading case for which the proposed model has been conceived. The new macro‐element is subsequently validated using cyclic and dynamic large‐scale laboratory tests of shallow foundations on dense sand, namely: the TRISEE cyclic tests, the Public Works Research Institute and CAMUS IV shaking table tests. Based on this comprehensive validation process against a set of independent experimental results, a unique set of macro‐element parameters for shallow foundations on dense sand is proposed, which can be used to perform predictive analyses by means of the present model. Copyright © 2011 John Wiley & Sons, Ltd.     

Vertical earthquake response of megawatt‐sized wind turbine with soil‐structure interaction effects

The dynamic response of a wind turbine on monopile is studied under horizontal and vertical earthquake excitations. The analyses are carried out using the finite element program SAP2000. The finite element model of the structure is verified against the results of shake table tests, and the earthquake response of the soil model is verified against analytical solutions of the steady‐state response of homogeneous strata. The focus of the analyses in this paper is the vertical earthquake response of wind turbines including the soil‐structure interaction effects. The analyses are carried out for both a non‐homogeneous stratum and a deep soil using the three‐step method. In addition, a procedure is implemented which allows one to perform coupled soil‐structure interaction analyses by properly tuning the damping in the tower structure. The analyses show amplification of the ground surface acceleration to the top of the tower by a factor of two. These accelerations are capable of causing damage in the turbine and the tower structure, or malfunctioning of the turbine after the earthquake; therefore, vertical earthquake excitation is considered a potential critical loading in design of wind turbines even in low‐to‐moderate seismic areas. Copyright © 2015 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.     

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

On the performance of lumped parameter models with gyro‐mass elements for the impedance function of a pile‐group supporting a single‐degree‐of‐freedom system

Lumped parameter models with a so called “gyro‐mass” element (GLPMs) have been proposed recently in response to a strong demand for efficiently and accurately representing frequency‐dependent impedance functions of soil–foundation systems. Although GLPMs are considered to be powerful tools for practical applications in earthquake engineering, some problems remain. For instance, although GLPMs show fairly close agreement with the target impedance functions, the accuracy of the transfer functions and the time‐histories of dynamic responses in structural systems comprising GLPMs have never been verified. Furthermore, no assessment has been performed on how much difference appears in the accuracy of dynamic responses obtained from GLPMs and those from conventional Kelvin–Voigt models comprising a spring and a dashpot arranged in parallel with various frequency‐independent constants. Therefore, in this paper, these problems are examined using an example of 2×4 pile groups embedded in a layered soil medium, supporting a single‐degree‐of‐freedom system subjected to ground motions. The results suggest that GLPMs are a new option for highly accurate computations in evaluating the dynamic response of structural systems comprising typical pile groups, rather than conventional Kelvin–Voigt models. Copyright © 2011 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.     

A fundamental dual‐zone continuum theory for dynamic soil–structure interaction

A continuum theory for an improved characterization of dynamic soil–structure interaction in the framework of three‐dimensional elastodynamics is presented. Effective in demonstrating the importance of integrating free‐field and near‐field effects under general soil and foundation conditions, a compact two‐zone delineation of the soil medium is proposed as a quintessential mechanics perspective for this class of problems. Sufficient to deliver a practical resolution of some perennial analytical and experimental conflicts, a fundamental formulation commensurate to a gradated unification of the homogenization approach and any sole free‐field inhomogeneous representation is developed and implemented computationally. Specialized to the problem of a rigid circular footing on sand, a nominal set of dynamic contact stress distributions and related impedance functions by the dual‐zone theory is included for theoretical and engineering evaluation. Through its comparison with benchmark analytical solutions and relevant physical measurements, the usage of the underlying conceptual platform as an advanced yet practical foundation for general dynamic soil–structure interaction is illustrated. Copyright © 2010 John Wiley & Sons, Ltd.     

Inelastic displacement ratios for soil–structure systems allowed to uplift

The simultaneous effects of soil–structure interaction, foundation uplift and inelastic behavior of the superstructure on total displacement response of soil–structure systems are investigated. The superstructure is modeled as an equivalent single‐degree‐of‐freedom system with bilinear behavior mounted on a rigid foundation resting on distributed tensionless Winkler springs and dampers. It is well known that the behavior of soil–structure systems can be well described using a limited number of nondimensional parameters. Here, by introducing two new parameters, the concept is extended to inelastic soil–structure systems in which the foundation is allowed to uplift. An extensive parametric study is conducted for a wide range of the key parameters through nonlinear time history analyses. It is shown that while uplifting soil–structure systems experience excessive displacements, in comparison with systems that are not allowed to uplift, ductility demand in the superstructure generally decreases owing to foundation uplift. A new inelastic displacement ratio (IDR) is proposed in conjunction with a nonlinear static analysis of uplifting soil–structure systems. Simplified expressions are also provided to estimate the proposed IDR. Copyright © 2014 John Wiley & Sons, Ltd.     

Optimum lateral load pattern for seismic design of elastic shear‐buildings incorporating soil–structure interaction effects

Recently, several new optimum loading patterns have been proposed by researchers for fixed‐base systems while their adequacy for soil–structure systems has not been evaluated yet. Through intensive dynamic analyses of multistory shear‐building models with soil–structure interaction subjected to a group of 21 artificial earthquakes adjusted to soft soil design spectrum, the adequacy of these optimum patterns is investigated. It is concluded that using these patterns the structures generally achieve near optimum performance in some range of periods. However, their efficiency reduces as soil flexibility increases especially when soil–structure interaction effects are significant. In the present paper, using the uniform distribution of damage over the height of structures, as the criterion, an optimization algorithm for seismic design of elastic soil–structure systems is developed. The effects of fundamental period, number of stories, earthquake excitation, soil flexibility, building aspect ratio, damping ratio and damping model on optimum distribution pattern are investigated. On the basis of 30,240 optimum load patterns derived from numerical simulations and nonlinear statistical regression analyses, a new lateral load pattern for elastic soil–structure systems is proposed. It is a function of the fundamental period of the structure, soil flexibility and structural slenderness ratio. It is shown that the seismic performance of such a structure is superior to those designed by code‐compliant or recently proposed patterns by researchers for fixed‐base structures. Using the proposed load pattern in this study, the designed structures experience up to 40% less structural weight as compared with the code‐compliant or optimum patterns developed based on fixed‐base structures. Copyright © 2012 John Wiley & Sons, Ltd.