Interactive behavior of soil–pile-superstructure system in transient state to liquefaction by means of large shake table tests

Shaking table tests were conducted by means of a large-scale laminar box with 4 m in length, 2 m in width and 2 m in height in order to investigate behavior of a soil-pile-superstructure system in liquefiable ground. A model two-storey structure, supported by a pile group, was set in a saturated sand deposit, and subjected to a sinusoidal base motion with increasing amplitude. Discussions are focused on the transient behavior until soil liquefaction occurs. Main interests are characteristics of springs used in a sway-rocking model and a multi-freedom lumped mass (MFLM) model that are frequently used in soil–pile interaction analysis. The spring constant in the sway-rocking model is represented by restoring force characteristics at the pile head, and that in the MFLM system is represented by an interaction spring connecting the pile to the free field. The transient state prior to soil liquefaction is shown to be important in the design of a pile because dynamic earth pressure shows peak response in this state. The reduction of the stiffness due to excess porewater generation and strain dependent nonlinear behavior is evaluated.     

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.     

Seismic response of a bridge–soil–foundation system under the combined effect of vertical and horizontal ground motions

Different levels of model sophistication have recently emerged to support seismic risk assessment of bridges, but mostly at the expense of neglecting the influence of vertical ground motions (VGMs). In this paper, the influence of VGMs on bridge seismic response is presented and the results are compared with the case of horizontal‐only excitations. An advanced finite element model that accounts for VGMs is first developed. Then, to investigate the effect of soil–structure interaction (SSI) including liquefaction potential, the same bridge with soil‐foundation and fixed boundary conditions is also analyzed. Results show that the inclusion of the VGMs has a significant influence on the seismic response, especially for the axial force in columns, normal force of bearings, and the vertical deck bending moments. However, VGMs do not have as much influence on the seismic demand of the pile cap displacements or pile maximum axial forces. Also, the significant fluctuation of the column axial force can reduce its shear and flexural capacity, and a heightened reversal of flexural effects may induce damage in the deck. In addition, relative to the fixed base case, SSI effects tend to reduce response quantities for certain ground motions while increasing demands for others. This phenomenon is explained as a function of the frequency content of the ground motions, the shift in natural vertical periods, and the VGM spectral accelerations at higher modes. Moreover, the mechanisms of liquefaction are isolated relative to SSI effects in nonliquefiable soils, revealing the influence of liquefaction on bridge response under VGMs. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of dynamic soil–structure interaction in shaking table testing

This paper provides an insight into the numerical simulation of soil–structure interaction (SSI) phenomena studied in a shaking table facility. The shaking table test is purposely designed to confirm the ability of the numerical substructure technique to simulate the SSI phenomenon. A model foundation–structure system with strong SSI potential is embedded in a dry bed of sand deposited within a purpose designed shaking-table soil container. The experimental system is subjected to a strong ground motion. The numerical simulation of the complete soil–foundation–structure system is conducted in the linear viscoelastic domain using the substructure approach. The matching of the experimental and numerical responses in both frequency and in time domain is satisfying. Many important aspects of SSI that are apparent in the experiment are captured by the numerical simulation. Furthermore, the numerical modelling is shown to be adequate for practical engineering design purposes.     

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.     

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     

Prediction of free field vibrations due to pile driving using a dynamic soil–structure interaction formulation

This paper presents a numerical model for the prediction of free field vibrations due to vibratory and impact pile driving. As the focus is on the response in the far field where deformations are relatively small, a linear elastic constitutive behaviour is assumed for the soil. The free field vibrations are calculated by means of a coupled FE–BE model based on a subdomain formulation. First, the case of vibratory pile driving is considered, where the contributions of different types of waves are investigated for several penetration depths. In the near field, the soil response is dominated by a vertically polarized shear wave, whereas in the far field, body waves are importantly attenuated and Rayleigh waves dominate the ground vibration. Second, the case of impact pile driving is considered. A linear wave equation model is used to estimate the impact force during the driving process. Apart from the response of a homogeneous halfspace, it is also investigated how the soil stratification influences the ground vibration for the case of a soft layer on a stiffer halfspace. When the penetration depth is smaller than the layer thickness, the layered medium has no significant influence on ground vibrations. However, when the penetration depth is larger than the layer thickness, the influence of the layered medium becomes more significant. The computed ground vibrations are finally compared with field measurements reported in the literature.     

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.     

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.     

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.     

Three-dimensional non-linear soil–building interaction analysis in the lakebed zone of Mexico city during the hypothetical Guerrero earthquake

The 1985 Michoacan earthquake (M=8·1) caused very severe damage to mid-rise buildings in the lakebed zone of Mexico City, which is approximately 400 km from the epicentre in the Pacific Ocean. In the present study, we perform a three-dimensional (3-D) non-linear soil–building interaction analysis for several types of low- to high-rise buildings during the hypothetical Guerrero earthquake, and try to understand the real cause of heavy damage to mid-rise buildings in the lakebed zone during the 1985 Michoacan earthquake. We make a reasonable estimation of the input earthquake motions and the local site effects. The non-linear soil-building interaction analysis explains the damage pattern observed during the 1985 earthquake, although other analyses do not. We realize that all the factors from the earthquake source to the building superstructure must be taken into account adequately. © 1998 John Wiley & Sons, Ltd.     

Inelastic dynamic analysis of RC bridges accounting for spatial variability of ground motion,site effects and soil–structure interaction phenomena. Part 2: Parametric study

The methodology for dealing with spatial variability of ground motion, site effects and soil–structure interaction phenomena in the context of inelastic dynamic analysis of bridge structures, and the associated analytical tools established and validated in a companion paper are used herein for a detailed parametric analysis, aiming to evaluate the importance of the above effects in seismic design. For a total of 20 bridge structures differing in terms of structural type (fundamental period, symmetry, regularity, abutment conditions, pier‐to‐deck connections), dimensions (span and overall length), and ground motion characteristics (earthquake frequency content and direction of excitation), the dynamic response corresponding to nine levels of increasing analysis complexity was calculated and compared with the ‘standard’ case of a fixed base, uniformly excited, elastic structure for which site effects were totally ignored. It is concluded that the dynamic response of RC bridges is indeed strongly affected by the coupling of the above phenomena that may adversely affect displacements and/or action effects under certain circumstances. Evidence is also presented that some bridge types are relatively more sensitive to the above phenomena, hence a more refined analysis approach should be considered in their case. Copyright @ 2003 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.     

Inelastic dynamic analysis of RC bridges accounting for spatial variability of ground motion,site effects and soil–structure interaction phenomena. Part 1: Methodology and analytical tools

During strong ground motion it is expected that extended structures (such as bridges) are subjected to excitation that varies along their longitudinal axis in terms of arrival time, amplitude and frequency content, a fact primarily attributed to the wave passage effect, the loss of coherency and the role of local site conditions. Furthermore, the foundation interacts with the soil and the superstructure, thus significantly affecting the dynamic response of the bridge. A general methodology is therefore set up and implemented into a computer code for deriving sets of appropriately modified time histories and spring–dashpot coefficients at each support of a bridge with account for spatial variability, local site conditions and soil–foundation–superstructure interaction, for the purposes of inelastic dynamic analysis of RC bridges. In order to validate the methodology and code developed, each stage of the proposed procedure is verified using recorded data, finite‐element analyses, alternative computer programs, previous research studies, and closed‐form solutions wherever available. The results establish an adequate degree of confidence in the use of the proposed methodology and code in further parametric analyses and seismic design. Copyright © 2003 John Wiley & Sons, Ltd.     

Fluid–structure interaction analysis of 3-D rectangular tanks by a variationally coupled BEM–FEM and comparison with test results

A variationally coupled BEM–FEM is developed which can be used to analyse dynamic response, including free-surface sloshing motion, of 3-D rectangular liquid storage tanks subjected to horizontal ground excitation. The tank structure is modelled by the finite element method and the fluid region by the indirect boundary element method. By minimizing a single Lagrange function defined for the entire system, the governing equation with symmetric coefficient matrices is obtained. To verify the newly developed method, the analysis results are compared with the shaking-table test data of a 3-D rectangular tank model and with the solutions by the direct BEM–FEM. Analytical studies are conducted on the dynamic behaviour of 3-D rectangular tanks using the method developed. In particular, the characteristics of the sloshing response, the effect of the rigidity of adjacent walls on the dynamic response of the tanks and the orthogonal effects are investigated. © 1998 John Wiley & Sons, Ltd.     

Response analysis of the Higashi–Kobe Bridge and surrounding soil in the 1995 Hyogoken–Nanbu Earthquake

This paper presents results of observation and analysis of the response of one of the longest cable-stayed bridges in the world to the Hyogoken–Nanbu (Kobe) Earthquake of 17 January 1995. It is determined that interaction of the foundations of the bridge towers with the supporting soil plays a decisive role in the overall structural behaviour. The key factor governing the changes of the soil properties at this site is pore-water pressure buildup, which results in liquefaction of the saturated surface soil layers under large dynamic loads. Models of the soil and structure are created and initially validated by accurately simulating the system response to a small earthquake. Soil parameters reflecting the pore-water pressure buildup in the strong earthquake are determined by an advanced non-linear effective stress analysis, combining the Ramberg–Osgood model of stress–strain dependence with a pore pressure model based on shear work concept. They are utilized to investigate and simulate the interaction of the foundation and the supporting soil using the program SASSI with the flexible volume substructuring approach. The results show a good agreement with the observations and have useful implications to the scientific and engineering practice. © 1998 John Wiley & Sons, Ltd.     

Simulation analyses of buildings damaged in the 1995 Kobe,Japan, earthquake considering soil–structure interaction

A series of studies was conducted on three buildings of steel reinforced concrete structures with RC shear walls damaged in the 1995 Hyogo-ken Nanbu earthquake. These buildings are located in an area where structural damage centred around. Two of these buildings suffered severe damage, while the third was not structurally damaged. Our studies deal with site inspections, including micro-tremor measurement of buildings, the evaluation of input motions, and the response analyses considering soil–structure interaction. The results of simulation analyses of the two severely damaged buildings correspond to their actual damage state. From the response analyses of the one slender building with no structural damage, it was concluded that uplifting is the main reason it did not suffer any structural damage. Through these studies, the importance of soil–structure interaction and effective input motion is fully understood. Copyright © 1999 John Wiley & Sons, Ltd.     

Use of statistical analysis to formulate conceptual models of geochemical behavior: water chemical data from the Botucatu aquifer in São Paulo state, Brazil

Water chemical data from the Botucatu Sandstone aquifer in the São Paulo State part of the Paraná Basin, Brazil, was evaluated using geochemical methods and two statistical analyses: cluster analysis and factor analysis. The results were used to develop a conceptual geochemical model, in which three geochemical regions were identified, and their chemical behavior was modeled. The characteristic chemicals, changing from the recharge area to the center of the basin, are: SiO₂—(HCO₃⁻ and Ca²⁺)—(Na⁺, CO₃²⁻, and SO₄²⁻). The distribution of the chemicals is interpreted as controlled by different water–rock interaction processes in the different regions. In the recharge area, dissolution of alkali–feldspar minerals in the sandstone is the main reaction observed; in the mid-section of the basin, calcite dissolution results in high calcium and bicarbonate concentration; in the center of the basin, leakage from underlying layers becomes the governing factor.