Nonlinear response of surface soil and NTT building due to soil–structure interaction during the 1995 Hyogo-ken Nanbu (Kobe) earthquake

The 1995 Hyogo-ken Nanbu (Kobe) earthquake brought about enormous damage to structures in the Hanshin and Awaji areas. In this paper the importance of investigating the relationship between ground motion and structural damage is pointed out.

Strong seismic motion was observed at the NTT (Nippon Telegraph and Telephone) Building during this earthquake. The structural damage to this building was relatively slight. In order to evaluate the relationship between ground motion and structural damage, it is necessary to assess the effects of the soil–structure interaction. In this study, the seismic response of the building and of the surface soil were evaluated by means of a nonlinear soil–structure interaction analysis using FEM.

It was found that, the nonlinearity of surface soil near the building had a great effect on the soil–structure interaction, especially the rocking of the building.     

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.     

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.     

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.     

Resistivity mapping using the VLF-MT method around surface fault ruptures of the 1995 Hyogo-ken Nanbu earthquake, Japan

Abstract Distinctive fault ruptures, the Nojima Fault and Ogura Fault, appeared along the northwestern coast of Awaji Island at the time of the 1995 Hyogo-ken Nanbu earthquake (Kobe earthquake). In order to delineate the shallow resistivity structures around the faults just after they formed, Very Low Frequency Magnetotelluric (VLF-MT) surveys were made at five sites along the Nojima Fault and at one site along the Ogura Fault. Fourteen transects were made at the one site on the Ogura Fault, and another transect covers the area between the two faults. Changes in apparent resistivity or phase, or both, commonly occur when crossing the surface location of one of the faults, except for the northern transects at OGR-0 on the Ogura Fault. Apparent resistivity values of less than 100 Ωm were observed for Tertiary and Quaternary sediments and values larger than 200 Ωm for granitic rocks. The resistivity structures are related to the morphological characteristics of the fault ruptures. Remarkably conductive zones (less than 10 Ωm in apparent resistivity and 30–40 m in width) were found where the surface displacement is distinct and prominent along a single fault plane. If remarkably conductive zones were formed at the time of the 1995 Hyogo-ken Nanbu earthquake, the results provide a good constraint on the dimensions of a conductive zone near the surface that was made by one earthquake. Alternatively, if characteristic resistivity structures existed prior to the earthquake, the conductive zone was probably formed by some tens of earthquakes in relatively modern times. In this case, this phenomenon is inferred to be a concentration of fracturing in a narrow zone and is associated with the formation of clay minerals, which enhance rock conductivity.     

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.     

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.     

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.     

In situ stress measurements in NIED boreholes in and around the fault zone near the 1995 Hyogo-ken Nanbu earthquake, Japan

Abstract The 1995 Hyogo-ken Nanbu (Kobe) earthquake, M 7.2, occurred along the north-east–south-west trending Rokko–Awaji Fault system. Three boreholes of 1001 m, 1313 m and 1838 m deep were drilled in the vicinity of the epicenter of the earthquake. Each borehole is located at characteristic sites in relation to active faults and the aftershock distribution. In particular, the Nojima–Hirabayashi borehole [Hirabayashi National Research Institute for Earth Science and Disaster Prevention (NIED) drilling] in Awaji Island was drilled to a depth of 1838 m, approximately 320 m southeast from the surface rupture of the Nojima Fault, and it crosses fracture zones below a depth of 1140 m. In situ stress measurements by the hydraulic fracturing method were conducted in these boreholes within 1.5 years after the earthquake. Measurement results suggest the following: (i) Differential stress values are very small, approximately 10 MPa at a depth of 1000 m at each site; (ii) the orientation of maximum horizontal compression is almost the same in the boreholes, perpendicular to the surface trace of the faults, north-west–south-east; (iii) fault types estimated from the state of stress differ among these sites; and (iv) the differential stress value just beneath the fault fracture zone decreases abruptly to one-half of that above the fault zone in the Hirabayashi NIED drilling. These features support the idea that the shear stress along the Rokko–Awaji Fault system decreased to a low level just after the earthquake.     

Influence of the soil–structure interaction on the fundamental period of buildings

This paper includes an investigation of the influence of the soil–structure interaction (SSI) on the fundamental period of buildings. The behaviour of both the soil and the structure is assumed to be elastic. The soil‐foundation system is modelled using translational and rotational discrete springs. Analysis is first conducted for one‐storey buildings. It shows that the influence of the SSI on the fundamental frequency of building depends on the soil–structure relative rigidity K_ss. Analysis is then extended for multi‐storey buildings. It allows the generalization of the soil–structure relative rigidity K_s to such complex structures. Charts are proposed for taking into account the influence of the SSI in the calculation of the fundamental frequency of a wide range of buildings. Copyright © 2007 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.     

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

Damage‐based seismic planar pounding analysis of adjacent symmetric buildings considering inelastic structure–soil–structure interaction

In cities and urban areas, building structures located at close proximities inevitably interact under dynamic loading by direct pounding and indirectly through the underlying soil. Majority of the previous adjacent building pounding studies that have taken the structure–soil–structure interaction (SSSI) problem into account have used simple lumped mass–spring–dashpot models under plane strain conditions. In this research, the problem of SSSI‐included pounding problem of two adjacent symmetric in plan buildings resting on a soft soil profile excited by uniaxial earthquake loadings is investigated. To this end, a series of SSSI models considering one‐directional nonlinear impact elements between adjacent co‐planar stories and using a method for direct finite element modeling of 3D inelastic underlying soil volume has been developed to accurately study the problem. An advanced inelastic structural behavior parameter, the seismic damage index, has been considered in this study as the key nonlinear structural response of adjacent buildings. Based on the results of SSSI and fixed base case analyses presented herein, two main problems are investigated, namely, the minimum building separation distance for pounding prevention and seismic pounding effects on structural damage in adjacent buildings. The final results show that at least three times, the International Building Code 2009 minimum distance for building separation recommended value is required as a clear distance for adjacent symmetric buildings to prevent the occurrence of seismic pounding. At the International Building Code‐recommended distance, adjacent buildings experienced severe seismic pounding and therefore significant variations in storey shear forces and damage indices. Copyright © 2016 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.     

Estimation of strong motion distribution in the 1995 Kobe earthquake based on building damage data

The 1995 Kobe earthquake caused unprecedented damage to buildings and civil infrastructures in the city of Kobe and its surrounding areas. In order to evaluate the structural damage in this area due to the earthquake, it is important to estimate the distribution of earthquake ground motion. However, since the number of strong ground motion records is not enough in the heavily damaged areas, it is necessary to estimate the distribution using other data sources. In this paper, the fragility curves for low‐rise residential buildings were constructed using the recorded motions and the building damage data from the intensive field survey by the AIJ and CPIJ group. The fragility curves obtained were then employed to estimate the strong motion distribution in the district level for Kobe and the surrounding areas during the earthquake. The results may be useful to investigate the various damages caused by the earthquake. Copyright © 2001 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.     

Bound of earthquake input energy to soil–structure interaction systems

A new critical excitation method is developed for soil–structure interaction systems. In contrast to previous studies considering amplitude nonstationarity only, no special constraint of input motions is needed on nonstationarity. The input energy to the soil–structure interaction system during an earthquake is introduced as a new measure of criticality. It is demonstrated that the input energy expression can be of a compact form via the frequency integration of the product between the input component (Fourier amplitude spectrum) and the structural model component (so-called energy transfer function). With the help of this compact form, it is shown that the formulation of earthquake input energy in the frequency domain is essential for solving the critical excitation problem and deriving a bound on the earthquake input energy for a class of ground motions. The extension of the concept to MDOF systems is also presented.     

Outline of the fault zone drilling project by NIED in the vicinity of the 1995 Hyogo-ken Nanbu earthquake, Japan

Abstract Three boreholes, 1001 m, 1313 m and 1838 m deep, were drilled by the National Research Institute for Earth Science and Disaster Prevention (NIED) in the vicinity of the epicenter of the 1995 Hyogo-ken Nanbu (Kobe) earthquake to investigate tectonic and material characteristics near and in active faults. Using these boreholes, an integrated study of the in situ stress, heat flow, and material properties of drill cores and crustal resistivity was conducted. In particular, the Nojima–Hirabayashi borehole was drilled to a depth of 1838 m and directly intersected the Nojima Fault, and three possible fault strands were detected at depths of 1140 m, 1313 m and 1800 m. Major results obtained from this study include the following: (i) shear stress around the fault zone is very small, and the orientation of the maximum horizontal compression is perpendicular to the surface trace of faults; (ii) from the results of a heat flow study, the lower cut-off depth of the aftershocks was estimated to be roughly 300°C; (iii) cores were classified into five types of fault rocks, and an asymmetric distribution pattern of these fault rocks in the fracture zones was identified; (iv) country rock is characterized by a very low permeability and high strength; and (v) resistivity structure can be explained by a model of a fault extending to greater depths but with low resistivity.     

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.