A parallel finite–infinite element model for two-dimensional soil–structure interaction problems

A parallel soil–structure interaction (SSI) model is presented for applications on distributed computer systems. Substructring method is applied to the SSI system and a coupled finite–infinite element based parallel computer program is developed. In the SSI system, infinite elements are used to represent the soil which extends to infinity. In this case, a large finite element mesh is required to define the near field for reliable predictions. The resulting large-scale problems are solved on distributed computer systems in this study. The domain is represented by separated substructures and an interface. The number of substructures are determined by the available processors in the parallel platform. To avoid the formation of large interface equations, smaller interface equations are distributed to processors while substructure contributions are performed. This saves a lot of memory storage and computational effort. Direct solution techniques are used for the solution of interface and substructure equation systems. The program is investigated through some example problems. The example problems exposed the need for solving large-scale problems in order to reach better results. The results of the example problems demonstrated the benefits of the parallel SSI algorithm.     

Real-time dynamic hybrid testing for soil–structure interaction analysis

Simulating dynamic soil–structure interaction (SSI) problems is a challenge when using a shaking table because of the semi-infinity of soil foundations. This paper develops real-time dynamic hybrid testing (RTDHT) for SSI problems in order to consider the radiation damping effect of the semi-infinite soil foundation using a shaking table. Based on the substructure concept, the superstructure is physically tested and the semi-infinite foundation is numerically simulated. Thus, the response of the entire system considering the dynamic SSI is obtained by coupling the numerical calculation of the soil and the physical test of the superstructure. A two-story shear frame on a rigid foundation was first tested to verify the developed RTDHT system, in which the top story was modeled as the physical substructure and the bottom story was the numerical substructure. The RTDHT for a two-story structure mounted on soil foundation was then carried out on a shaking table while the foundation was numerically simulated using a lumped parameter model. The dynamic responses, including acceleration and shear force, were obtained under soft and hard soil conditions. The results show that the soil–structure interaction should be reasonably taken into account in the shaking table testing for structures.     

Real-time dynamic hybrid testing for soil–structure interaction analysis

Simulating dynamic soil–structure interaction (SSI) problems is a challenge when using a shaking table because of the semi-infinity of soil foundations. This paper develops real-time dynamic hybrid testing (RTDHT) for SSI problems in order to consider the radiation damping effect of the semi-infinite soil foundation using a shaking table. Based on the substructure concept, the superstructure is physically tested and the semi-infinite foundation is numerically simulated. Thus, the response of the entire system considering the dynamic SSI is obtained by coupling the numerical calculation of the soil and the physical test of the superstructure. A two-story shear frame on a rigid foundation was first tested to verify the developed RTDHT system, in which the top story was modeled as the physical substructure and the bottom story was the numerical substructure. The RTDHT for a two-story structure mounted on soil foundation was then carried out on a shaking table while the foundation was numerically simulated using a lumped parameter model. The dynamic responses, including acceleration and shear force, were obtained under soft and hard soil conditions. The results show that the soil–structure interaction should be reasonably taken into account in the shaking table testing for structures.     

Nonlinear 3D finite element analysis of soil–pile–structure interaction system subjected to horizontal earthquake excitation

The dynamic response of a seismic soil–pile–structure interaction (SSPSI) system is investigated in this paper by conducting nonlinear 3D finite element numerical simulations. Nonlinear behaviors such as non-reflecting boundary condition and soil–pile–structure interaction modeled by the penalty method have been taken into account. An equivalent linear model developed from the ground response analysis and the modified Drucker–Prager model are separately used for soil ground. A comparison of the two models shows that the equivalent linear soil model results in an underestimated acceleration response of the structure under this ground shaking and the soil behavior should be considered as a fully-nonlinear constitutive model in the design process of the SSPSI system. It was also observed that the dynamic response of the system is greatly affected by the nonlinearity of soil–pile interface and is not sensitive to the dilation angle of the soil. Furthermore, the effect of the presence of pile foundations on SSPSI response is also analyzed and discussed.     

Observed dynamic soil–structure interaction in scale testing of offshore wind turbine foundations

Monopile foundations have been commonly used to support offshore wind turbine generators (WTGs), but this type of foundation encounters economic and technical limitations for larger WTGs in water depths exceeding 30 m. Offshore wind farm projects are increasingly turning to alternative multipod foundations (for example tetrapod, jacket and tripods) supported on shallow foundations to reduce the environmental effects of piling noise. However the characteristics of these foundations under dynamic loading or long term cyclic wind turbine loading are not fully understood. This paper summarises the results from a series of small scaled tests (1:100, 1:150 and 1:200) of a complete National Renewable Energy Laboratory (NREL) wind turbine model on three types of foundations: monopiles, symmetric tetrapod and asymmetric tripod. The test bed used consists of either kaolin clay or sand and up to 1.4 million loading cycles were applied. The results showed that the multipod foundations (symmetric or asymmetric) exhibit two closely spaced natural frequencies corresponding to the rocking modes of vibration in two principle axes. Furthermore, the corresponding two spectral peaks change with repeated cycles of loading and they converge for symmetric tetrapods but not for asymmetric tripods. From the fatigue design point of view, the two spectral peaks for multipod foundations broaden the range of frequencies that can be excited by the broadband nature of the environmental loading (wind and wave) thereby impacting the extent of motions. Thus the system lifespan (number of cycles to failure) may effectively increase for symmetric foundations as the two peaks will tend to converge. However, for asymmetric foundations the system life may continue to be affected adversely as the two peaks will not converge. In this sense, designers should prefer symmetric foundations to asymmetric foundations.     

Dynamic response of structures subjected to pounding and structure–soil–structure interaction

During strong earthquakes, adjacent structures with non-sufficient clear distances collide with each other. In addition to such a pounding, cross interaction of adjacent structures through soil can exchange the vibration energy between buildings and make the problem even more complex. In this paper, effects of both of the mentioned phenomena on the inelastic response of selected steel structures are studied. Number of stories varied between 3 and 12 and different clear distances up to the seismic codes prescribed value are considered. The pounding element is modeled within Opensees. A coupled model of springs and dashpots is utilized for through-the-soil interaction of the adjacent structures, for two types of soft soils. The pounding force, relative displacements of stories, story shears, and plastic hinge rotations are compared for different conditions as the maximum responses averaged between seven consistent earthquakes. As a result, simultaneous effects of pounding and structure–soil–structure interaction are discussed.     

Bounds for natural frequencies,Dunkerley's formula and application to soil–structure interaction

A new expression for a lower bound for the fundamental natural frequency of a structure is presented and its relationship with Dunkerley's formula is examined. A related improved estimate for the fundamental frequency is also obtained and applied to obtain an approximate expression for the fundamental frequency of a soil–structure system including a 1-DOF structure, a flat rigid foundation, and an elastic soil. The resulting approximation is then compared with the now standard result first obtained by Jennings and Bielak [6]     

The near-field method for dynamic analysis of structures on soft soils including inelastic soil–structure interaction

The problem of soil–structure interaction analysis with the direct method is studied. The direct method consists of explicitly modeling the surrounding soil to bedrock and the structure resting on the soil. For the soil medium, usually the traditional equivalent linear method with a reduced shear modulus and an increased damping ratio for the soil is used. However, this method does not work in the vicinity of foundation where the soil behavior is highly nonlinear because of presence of large strains. This research proposes a modified equivalent linear method with a further reduction of the soil shear modulus in the near-field of foundation that results in validity of using the equivalent linear method throughout. For regular short, intermediate and tall structures resting on such soft soils, a series of dynamic time-history analysis is implemented using earthquake records scaled to a sample design spectrum and the nonlinear structural responses are compared for different assumptions of soil behavior including the elasto-plastic Mohr–Coulomb, the traditional equivalent linear, and the proposed modified equivalent linear method. This analysis validates the proposed method.     

Conversion factors for design spectral accelerations including soil–structure interaction

In this paper maximum response of a single degree of freedom system resting on a flexible base is determined under consistent earthquakes and the results are presented as acceleration spectra including soil–structure interaction (SSI). Flexibility of base is modeled using frequency-dependent springs and dampers. The spring–damper coefficients are calculated for the desired natural mode of vibration of a multi-degree-of-freedom system. Consistency of earthquakes is maintained considering their magnitude, distance, local soil type, and return period. The latter parameter is accounted for by the use of earthquake categories identified by their similar spectral values. Ratio of spectral acceleration modification factors with SSI from this study to those calculated using the ASCE 7-10 procedure are determined for each case. Examination of the resulting curves shows that the mentioned code is conservative/non-conservative in estimation of spectral responses with SSI in certain cases for the lower/higher modes of vibration. The code’s procedure is modified using the developed curves for a conversion factor.     

Dynamic soil–structure interaction of monopile supported wind turbines in cohesive soil

Offshore wind turbines supported on monopile foundations are dynamically sensitive because the overall natural frequencies of these structures are close to the different forcing frequencies imposed upon them. The structures are designed for an intended life of 25 to 30 years, but little is known about their long term behaviour. To study their long term behaviour, a series of laboratory tests were conducted in which a scaled model wind turbine supported on a monopile in kaolin clay was subjected to between 32,000 and 172,000 cycles of horizontal loading and the changes in natural frequency and damping of the model were monitored. The experimental results are presented using a non-dimensional framework based on an interpretation of the governing mechanics. The change in natural frequency was found to be strongly dependent on the shear strain level in the soil next to the pile. Practical guidance for choosing the diameter of monopile is suggested based on element test results using the concept of volumetric threshold shear strain.     

A lattice Boltzmann-finite element model for two-dimensional fluid–structure interaction problems involving shallow waters

In this paper, a numerical method for the modeling of shallow waters interacting with slender elastic structures is presented. The fluid domain is modeled through the lattice Boltzmann method, while the solid domain is idealized by corotational beam finite elements undergoing large displacements. Structure dynamics is predicted by using the time discontinuous Galerkin method and the fluid–structure interface conditions are handled by the Immersed Boundary method. An explicit coupling strategy to combine the adopted numerical methods is proposed and its effectiveness is tested by computing the error in terms of the energy that is artificially introduced at the fluid–solid interface.     

Simulation of soil–foundation–structure interaction of Hualien large-scale seismic test using dynamic centrifuge test

Understanding the soil–structure interaction (SSI) mechanism is crucial in the seismic design of nuclear power plant (NPP) containment systems. Although the numerical analysis method is generally used in seismic design, there is a need for experimental verification for the reliable estimation of SSI behavior. In this study a dynamic centrifuge test was performed to simulate the SSI behavior of a Hualien large-scale seismic test (LSST) during the Chi-Chi earthquake. To simulate the soil profile and dynamic soil properties of the Hualien site, a series of resonant column (RC) tests was performed to determine the model soil preparation conditions, such as the compaction density and the ratio of soil–gravel contents. The variations in the shear wave velocity (V_S) profiles of the sand, gravel, and backfill layers in the model were estimated using the RC test results. During the centrifuge test, the V_S profiles of the model were evaluated using in-flight bender element tests and compared with the in-situ V_S profile at Hualien. The containment building model was modeled using aluminum and the proper scaling laws. A series of dynamic centrifuge tests was performed with a 1/50 scale model using the base motion recorded during the Chi-Chi-earthquake. In the soil layer and foundation level, the centrifuge test results were similar to the LSST data in both the time and frequency domains, but there were differences in the structure owing to the complex structural response as well as the material damping difference between the concrete in the prototype and aluminum in the model. In addition, as the input base motion amplitude was increased to a maximum value of 0.4g (prototype scale), the responses of the soil and containment model were measured. This study shows the potential of utilizing dynamic centrifuge tests as an experimental modeling tool for site specific SSI analyses of soil–foundation–NPP containment system.     

Pile–soil–pile interaction in pile groups with batter piles under dynamic loads

A simple analytical solution is developed for computing the dynamic interaction tensor for floating pile groups with batter piles. For this purpose, the governing differential equations are derived for an unloaded batter floating pile closely spaced to another loaded pile with the same properties. The reaction of soil against pile deformation is modeled by the springs and dashpots along the length of the pile. The soil is assumed linear viscous elastic and the pile behavior is linear elastic. The closed form solutions of governing equations are obtained using appropriate boundary conditions. The results are verified and compared with rigorous and approximate analytical solutions for vertical piles. The proposed method can be readily applied by engineers in the design of pile groups with batter piles.     

The effect of soil–structure interaction on the seismic risk to buildings

This paper studies the effect of soil–structure interaction (SSI) on the seismic risk estimates of buildings. Risk, in this context, denotes the probability distribution of seismic monetary loss due to structural and nonstructural damage. The risk analysis here uncovers the probability that SSI is beneficial, detrimental, or uninfluential on seismic losses. The analyses are conducted for a wide range of buildings with different structural systems, numbers of stories, and foundation sizes on various soil types. A probabilistic approach is employed to account for prevailing sources of uncertainty, i.e., those in ground motion and in the properties of the soil–structure system. In this approach, probabilistic models are employed to predict the response, damage, and repair cost of buildings. To properly account for the ground motion uncertainty, a suite of nearly 7000 accelerograms recorded on soil is employed. It is concluded that structures on very soft soils are extremely likely to incur smaller losses due to SSI, which is in line with the common belief that SSI is a favorable effect for such systems. However, the results for buildings on moderately soft soils reveal a considerable probability, up to 0.4, that SSI has an adverse effect on the structure and increases the seismic losses.     

Effects of soil–structure interaction on the dynamic properties and seismic response of piled structures

This paper presents a simple and stable procedure for the estimation of periods and dampings of piled shear buildings taking soil–structure interaction into account. A substructuring methodology that includes the three-dimensional character of the foundations is used. The structure is analyzed as founded on an elastic homogeneous half-space and excited by vertically incident S waves. The strategies proposed in the literature to estimate the period and damping are revised, and a modified strategy is proposed including crossed impedances and all damping terms. Ready-to-use graphs are presented for the estimation of flexible-base period and damping in terms of their fixed-base values and the system configuration. Maximum shear forces together with base displacement and rocking peak response are also provided. It is shown that cross-coupled impedances and kinematic interaction factors need to be taken into account to obtain accurate results for piled buildings.     

Nonlinear seismic behaviour of wall-frame dual systems accounting for soil–structure interaction

The aim of this paper is to study the effects of soil–structure interaction on the seismic response of coupled wall-frame structures on pile foundations designed according to modern seismic provisions. The analysis methodology based on the substructure method is recalled focusing on the modelling of pile group foundations. The nonlinear inertial interaction analysis is performed in the time domain by using a finite element model of the superstructure. Suitable lumped parameter models are implemented to reproduce the frequency-dependent compliance of the soil-foundation systems. The effects of soil–structure interaction are evaluated by considering a realistic case study consisting of a 6-storey 4-bay wall-frame structure founded on piles. Different two-layered soil deposits are investigated by varying the layer thicknesses and properties. Artificial earthquakes are employed to simulate the earthquake input. Comparisons of the results obtained considering compliant base and fixed base models are presented by addressing the effects of soil–structure interaction on displacements, base shears, and ductility demand. The evolution of dissipative mechanisms and the relevant redistribution of shear between the wall and the frame are investigated by considering earthquakes with increasing intensity. Effects on the foundations are also shown by pointing out the importance of both kinematic and inertial interaction. Finally, the response of the structure to some real near-fault records is studied. Copyright © 2012 John Wiley & Sons, Ltd.     

A model for the seismic analysis of a periodic viaduct when considering the pile–soil–structure interaction

In this study, a new model is developed for the aseismic design of a periodic viaduct when the pile–soil–structure interaction is considered. To account for the influence of the pile–soil–structure interaction, a wavenumber domain boundary element method (WDBEM) model for the periodic pile row supporting the viaduct is developed using the sequence Fourier transform as well as the boundary element method for the elastic medium. By using the WDBEM model for the pile row, the transfer matrices for the beams and piers, the joint conditions at the beam–beam–pier (BBP) junction as well as the periodicity condition for the viaduct, the wavenumber domain response of the periodic viaduct to spatially harmonic waves is determined. Based on the wavenumber domain response of the viaduct, the space-domain response of the viaduct to an arbitrary seismic wave can be obtained by invoking the inverse sequence Fourier transform method. Numerical results show that when the periodic viaduct is exposed to the spatially harmonic wave, resonances may occur at the bounding frequencies of the passbands of the characteristic waves of the viaduct. Also, it is found that the coincidence between the traveling seismic wave and characteristic waves of the viaduct will generate additional resonant frequencies located in passbands of the characteristic waves.     

Seismic soil–pile–structure kinematic and inertial interaction—A new beam approach

The main purpose of this study is to investigate the accuracy of an advanced beam model for the soil–pile–structure kinematic and inertial interaction and demonstrate its efficiency and advantages compared to other commonly used beam or solid models. Within this context, a Beam on Nonlinear Winkler Foundation model is adopted based on the Boundary Element Method (BEM), accounting for the effects induced by geometrical nonlinearity, rotary inertia and shear deformation, employing the concept of shear deformation coefficients. The soil nonlinearity is taken into consideration by means of a hybrid spring configuration consisting of a nonlinear (p–y) spring connected in series to an elastic spring–damper model. The nonlinear spring captures the near-field plastification of the soil while the spring–damper system (Kelvin–Voigt element) represents the far-field viscoelastic character of the soil. An extensive case study is carried out on a pile-column–deck system of a bridge, found in two cohesive layers of sharply different stiffness and subjected to various earthquake excitations, providing insight to several phenomena. The results of the proposed model are compared with those obtained from a Beam-FE solution as well as from a rigorous fully three-dimensional (3-D) continuum FE scheme.     

Effects of the site dynamic characteristics on soil–structure interaction (I): Incident SH-Waves

In this paper, the role that the site dynamic characteristics play in soil–structure interaction is studied on a simple model in which the site is represented as a soil layer over bedrock (half-space), and using the indirect boundary-element method (IBEM). For the purpose of comparison with published analytical solutions, the structure is represented as a shear wall supported by a semi-circular rigid foundation, subjected to incident plane SH waves. The accuracy of the method is verified, numerical results are analyzed, and the model response is compared with earthquake observations at the Hollywood Storage Building. It is shown that the effects of dynamic soil–structure interaction may become more significant near the characteristic frequencies of the site, and that the resonance of the system shifts to lower frequencies. The thickness of the soil layer, the stiffness of the bedrock, and the mass and the stiffness of the superstructure all influence the values of the system frequencies and system amplitudes.