Numerical simulations of shake-table experiment for dynamic soil-pile-structure interaction in liquefiable soils

A shake-table experiment on pile foundations in liquefi able soils composed of liquefi able sand and overlying soft clay is studied. A three-dimensional(3D) effective stress fi nite element(FE) analysis is employed to simulate the experiment. A recently developed multi-surface elasto-plastic constitutive model and a fully coupled dynamic inelastic FE formulation(u-p) are used to model the liquefaction behavior of the sand. The soil domains are discretized using a solid-fl uid fully coupled(u-p) 20-8 noded brick element. The pile is simulated using beam-column elements. Upon careful calibration, very good agreement is obtained between the computed and the measured dynamic behavior of the ground and the pile. A parametric analysis is also conducted on the model to investigate the effect of pile-pinning, pile diameter, pile stiffness, ground inclination angle, superstructure mass and pile head restraints on the ground improvement. It is found that the pile foundation has a noticeable pinning effect that reduces the lateral soil displacement. It is observed that a larger pile diameter and fi xed pile head restraints contribute to decreasing the lateral pile deformation; however, a higher ground inclination angle tends to increase the lateral pile head displacements and pile stiffness, and superstructure mass seems to effectively infl uence the lateral pile displacements.     

液化侧扩流场地桥梁群桩效应分析

基于u-p有限元公式模拟饱和砂土中水和土颗粒完全耦合效应,建立液化侧向流场地群桩动力反应分析的三维数值模型。模型中,砂土采用多屈服面弹塑性本构模型模拟、黏土采用多屈服面运动塑性模型模拟,群桩在计算过程中保持线弹性状态;采用20节点的六面体单元和考虑孔压效应的20-8节点分别划分黏土层和饱和砂层;选用剪切梁边界处理计算域的人工边界,模拟地震过程中土层的剪切效应;应用瑞利阻尼考虑体系的阻尼效应。随后对比分析2×2群桩中各单桩的地震反应规律,结果表明,各单桩的弯矩、位移时程规律基本一致,峰值弯矩及峰值位移出现时刻滞后于输入加速度峰值时刻,上坡向桩的弯矩和位移峰值大于下坡向的桩的反应值。接着通过改变桩间距研究群桩效应,随着桩间距增加,群桩中各单桩的弯矩最大值均出现在土层分界处,且各单桩的弯矩、桩顶位移逐渐增大。最后给出液化侧向流场地群桩效应的基本原因,得出该类场地群桩抗震设计的基本认识。     

Shake table test of soil-pile groups-bridge structure interaction in liquefiable ground

<正>This paper describes a shake table test study on the seismic response of low-cap pile groups and a bridge structure in liquefiable ground.The soil profile,contained in a large-scale laminar shear box,consisted of a horizontally saturated sand layer overlaid with a silty clay layer,with the simulated low-cap pile groups embedded.The container was excited in three E1 Centra earthquake events of different levels.Test results indicate that excessive pore pressure(EPP) during slight shaking only slightly accumulated,and the accumulation mainly occurred during strong shaking.The EPP was gradually enhanced as the amplitude and duration of the input acceleration increased.The acceleration response of the sand was remarkably influenced by soil liquefaction.As soil liquefaction occurred,the peak sand displacement gradually lagged behind the input acceleration;meanwhile,the sand displacement exhibited an increasing effect on the bending moment of the pile,and acceleration responses of the pile and the sand layer gradually changed from decreasing to increasing in the vertical direction from the bottom to the top.A jump variation of the bending moment on the pile was observed near the soil interface in all three input earthquake events.It is thought that the shake table tests could provide the groundwork for further seismic performance studies of low-cap pile groups used in bridges located on liquefiable groun.     

Soil–pile-bridge structure interaction in liquefying ground using shake table testing

The evaluation of seismic pile response is particularly useful for geotechnical engineers involved in the design of foundations in liquefying site. Shake table testing was performed to study the dynamic interactive behavior of soil–pile foundations in liquefying ground under different shaking frequency and amplitude. The soil profile consisted of a clayey layer over liquefiable sand over clay. The model was tested with a series of El Centro earthquake motions with peak accelerations ranging from 0.15g to 0.50g, and time step from 0.006 to 0.02 s. Representative data, including time histories of accelerations and excess pore pressure ratios that characterize the important aspects of soil–pile interaction in liquefying ground are presented. The shaking frequency has no significant effect on the magnitudes of excess pore pressure ratio, ground and pile accelerations and pile bending moments. Excess pore pressure ratio, ground acceleration and pile acceleration, and pile bending moment largely depend on the shaking amplitude.     

Experimental and numerical investigation on the dynamic response of pile group in liquefying ground

The response of pile foundation in liquefiable sand reinforced by densification techniques remains a very complex problem during strong earthquakes. A shake-table experiment was carried out to investigate the behavior of a reinforced concrete low-cap pile group embedded in this type of ground. In this study, a three-dimensional (3D) finite element (FE) analysis of the experiment was conducted. The computed response of the soil-pile system was in reasonable agreement with the experimental results, highlighting some key characteristics. Then, a parametric study was performed to explore the influence of pile spacing, pile stiffness (EI), superstructure mass, sand permeability, and shaking characteristics of input motion on the behavior of the pile. The investigation demonstrated a stiffening behavior appearing in the liquefied mediumdense sand, and the pile group effect seemed negligible. Furthermore, the kinematic effect was closely connected with both EI and sand permeability. Nevertheless, the inertial effect was strongly influenced by the superstructure mass. Meanwhile, high frequency and large amplitude of the input motion could produced greater the pile’s moments. It is estimated that this case study could further enhance the current understanding of the behavior of low-cap pile foundations in liquefied dense sand.     

The effects of liquefaction-induced lateral spreading on pile foundations

Observations of pile foundation performance during previous earthquakes have shown that pile failure has been caused by lateral ground movements resulting from soil liquefaction. The recognition that lateral ground movements may play a critical role in pile performance during an earthquake has important implications for design and risk assessment, and requires that analytical models be devised to evaluate these potential problems.In this paper, parametric studies were conducted to estimate the maximum bending moments induced in piles subjected to lateral ground displacement. The results are summarized in charts using dimensionless parameters.The analyses reveal that the existence of a nonliquefiable layer at the ground surface can affect significantly the maximum bending moment of the pile. When a relatively thick nonliquefiable layer exists above a liquefiable layer, neither the material nonlinearity of the soil nor loss of soil stiffness within the liquefiable layer significantly affect the maximum bending moment. When the thickness of the liquefiable soils is greater than about three times that of an overlying intact layer, soil stiffness in the liquefiable layer must be chosen carefully when evaluating the maximum bending moment.     

Shake-table investigation of scoured RC pile-group-supported bridges in liquefiable and nonliquefiable soils

This paper presents results of one-g shake-table tests on scoured pile-group-supported bridge models in saturated (liquefiable) and dry (nonliquefiable) sands. The primary objective is to reveal the influence of liquefaction on seismic demands and failure mechanism of scoured bridges. To this end, two identical models, each consisting of a 2 × 2 reinforced concrete pile-group with a center-to-center spacing of 3 times pile diameter, a cap and a single pier with a lumped iron block, were constructed and embedded into saturated and dry sands, respectively, with the same scour depth of 4 times pile diameter. Typical test results, including excess pore pressure, acceleration and displacement demands are interpreted first, followed by the focus on curvature demands and associated seismic failure mechanism identification. Finally, inertial and kinematic effects on pile curvature demands are estimated using cross-correlation analyses. Results show that near-pile liquefied soils exhibit more remarkable dilation tendency as compared to far field. For bridges under the given scour depth, soil liquefaction tends to significantly affect the failure modes via transferring damage positions from pier bottom to pile head and meanwhile from underground pile to pile head. In addition, pile group effects appear to be significant in nonliquefiable soils while to be relatively inessential in liquefied soils. Moreover, the inertial effect is more prominent in nonliquefiable soils, while the kinematic effect itself generally appears to be more significant in liquefiable soils. The test results can be used to validate numerical models for future studies.     

Experimental investigation on seismic behavior of single piles in sandy soil

This paper describes a quasi-static test program featuring lateral cyclic loading on single piles in sandy soil. The tests were conducted on 18 aluminum model piles with different cross sections and lateral load eccentricity ratios, e/d, (e is the lateral load eccentricity and d is the diameter of pile) of 0, 4 and 8, embedded in sand with a relative density of 30% and 70%. The experimental results include lateral load-displacement hysteresis loops, skeleton curves and energy dissipation curves. Lateral capacity, ductility and energy dissipation capacity of single piles under seismic load were evaluated in detail. The lateral capacities and the energy dissipation capacity of piles in dense sand were much higher than in loose sand. When embedded in loose sand, the maximum lateral load and the maximum lateral displacement of piles increased as e/d increased. On the contrary, when embedded in dense sand, the maximum lateral load of piles decreased as e/d increased. Piles with a higher load eccentricity ratio experienced higher energy dissipation capacity than piles with e/d of 0 in both dense and loose sand. At a given level of displacement, piles with circular cross sections provided the best energy dissipation capacity in both loose and dense sand.     

Probabilistic buckling analysis of axially loaded piles in liquefiable soils

Recent research has demonstrated that axial load alone can cause a slender pile to fail by forming a plastic hinge, if soil surrounding the pile liquefies in an earthquake. This failure mechanism is due to buckling instability. Lateral loads from lateral spreading of the surrounding soil and/or inertia and imperfection inherent in a pile result in increased deflection, which can promote more rapid collapse. These effects are secondary to the basic requirement so that axially loaded piles passing through liquefiable soils should be checked against Euler's buckling in addition to the bending mechanism of failure, i.e. incorporation of P−Δ effect. While fewer large diameter piles are currently being used in modern construction practice (which performed better in liquefiable areas rather than multiple small diameter piles), there are many pile-supported structures where buckling considerations were not taken into account and therefore may need retrofitting. This paper develops a probabilistic tool which can be used for assessing the likelihood of a buckling failure of existing piled foundations due to a scenario earthquake. This tool can equally serve as a valuable decision-support tool for implementing earthquake risk mitigation measures. A case study is presented to show the applicability of the method.     

Study of pore pressure variation during liquefaction using two constitutive models for sand

Numerical analyses of liquefiable sand are presented in this paper. Liquefaction phenomenon is an undrained response of saturated sandy soils when they are subjected to static or dynamic loads. A fully coupled dynamic computer code is developed to predict the liquefaction potential of a saturated sandy layer. Coupled dynamic field equations of extended Biot's theory with u–P formulation are used to determine the responses of pore fluid and soil skeleton. Generalized Newmark method is employed for integration in time. The soil behavior is modelled by two constitutive models; a critical state two-surface plasticity model, and a densification model. A class ‘B’ analysis of a centrifuge experiment is performed to simulate the dynamic response of level ground sites. The results of the numerical analyses demonstrate the capability of the critical sate two-surface plasticity model in producing pore pressures that are consistent with observations of the behavior of liquefiable sand in the centrifuge test.     

Liquefaction-induced deformation of earthen embankments on non-homogeneous soil deposits under sequential ground motions

Damage of embankments during earthquakes is widely attributed to the liquefaction of foundation soil. Previous studies have investigated the dynamic response of embankments by mainly considering uniform sand foundation and a single earthquake event. However, the foundation of an embankment consists of many sublayers of soil from liquefiable sand to relatively impermeable layer, and during earthquakes a mainshock may trigger numerous aftershocks within a short time which may have the potential to cause additional damage to soil structures. Accordingly, the investigation of liquefaction-induced deformation of earthen embankments on various liquefiable foundation conditions under mainshock–aftershock sequential ground motions is carried out by a series of dynamic centrifuge tests in this study. The liquefiable foundation includes uniform sand profile, continuous layered soil profile, and non-homogeneous soil profiles. Effects of various foundation conditions on embankment deformations are compared and analyzed. From the test results, it is found that the embankment resting on non-homogeneous soil deposits suffer more damage compared to the uniform sand foundation of same relative density. The test results also suggest that the sequential ground motions have a significant effect on the accumulated deformation of embankment.     

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.     

双向地震作用下液化水平和倾斜场地-桩基-桥梁结构地震反应的差异研究

根据已经完成的液化侧向扩展场地-群桩基础-上部结构体系大型振动台试验,在有限元软件OpenSees中建立了可液化倾斜场地振动台试验的有限元模型。通过与试验结果对比,验证了数值模型的可靠性。基于此,建立了典型水平和倾斜液化场地-桩基-桥梁结构体系的数值模型,讨论了双向地震作用下水平和倾斜场地体系地震响应的差异,结果表明:相比水平场地,倾斜场地超孔隙水压力在峰值阶段波动幅度更大,土体的侧向位移增加明显,尤其是在饱和砂土中部位置;倾斜场地中桩基础的破坏程度更大,可液化层中部桩基曲率最大可增大约13倍,桩身水平位移显著增加;而水平场地桥墩曲率比倾斜场地桥墩曲率大,建议在液化场地桩基设计中应考虑场地倾斜带来的影响。      

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.     

A shake table investigation of dynamic behavior of pile supported bridges in liquefiable soil deposits

Bridges are a part of vital infrastructure,which should operate even after a disaster to keep emergency services running.There have been numerous bridge failures during major past earthquakes due to liquefaction.Among other categories of failures,mid span collapse(without the failure of abutments)of pile supported bridges founded in liquefiable deposits are still observed even in most recent earthquakes.This mechanism of collapse is attributed to the effects related to the differential elongation of natural period of the individual piers during liquefaction.A shake table investigation has been carried out in this study to verify mechanisms behind midspan collapse of pile supported bridges in liquefiable deposits.In this investigation,a typical pile supported bridge is scaled down,and its foundations pass through the liquefiable loose sandy soil and rest in a dense gravel layer.White noise motions of increasing acceleration magnitude have been applied to initiate progressive liquefaction and to characterize the dynamic features of the bridge.It has been found that as the liquefaction of the soil sets in,the natural frequency of individual bridge support is reduced,with the highest reduction occurring near the central spans.As a result,there is differential lateral displacement and bending moment demand on the piles.It has also been observed that for the central pile,the maximum bending moment in the pile will occur at a higher elevation,as compared to that of the interface of soils of varied stiffness,unlike the abutment piles.The practical implications of this research are also highlighted.     

Numerical analyses of pile performance in laterally spreading frozen ground crust overlying liquefiable soils

Lateral spread of frozen ground crust over liquefied soil has caused extensive bridge foundation damage in the past winter earthquakes. A shake table experiment was conducted to investigate the performance of a model pile in this scenario and revealed unique pile failure mechanisms. The modelling results provided valuable data for validating numerical models. This paper presents analyses and results of this experiment using two numerical modeling approaches: solid-fluid coupled finite element (FE) modeling and the beam-on-nonlinear-Winkler-foundation (BNWF) method. A FE model was constructed based on the experiment configuration and subjected to earthquake loading. Soil and pile response results were presented and compared with experimental results to validate this model. The BNWF method was used to predict the pile response and failure mechanism. A p-y curve was presented for modelling the frozen ground crust with the free-field displacement from the experiment as loading. Pile responses were presented and compared with those of the experiment and FE model. It was concluded that the coupled FE model was effective in predicting formation of three plastic hinges at ground surface, ground crust-liquefiable soil interface and within the medium dense sand layer, while the BNWF method was only able to predict the latter two.     

A simplified method for analyzing the fundamental frequency of monopile supported offshore wind turbine system design

Preliminary design of offshore wind turbines requires high precision simplified methods for the analysis of the system fundamental frequency. Based on the Rayleigh method and Lagrange's Equation, this study establishes a simple formula for the analysis of system fundamental frequency in the preliminary design of an offshore wind turbine with a monopile foundation. This method takes into consideration the variation of cross-section geometry of the wind turbine tower along its length, with the inertia moment and distributed mass both changing with diameter. Also the rotational flexibility of the monopile foundation is mainly considered. The rigid pile and elastic middle long pile are calculated separately. The method is validated against both FEM analysis cases and field measurements, showing good agreement. The method is then used in a parametric study, showing that the tower length Lt, tower base diameter d0, tower wall thickness δt, pile diameter db and pile length Lb are the major factors influencing the fundamental frequency of the offshore wind turbine system. In the design of offshore wind turbine systems, these five parameters should be adjusted comprehensively. The seabed soil condition also needs to be carefully considered for soft clay and loose sand.     

Factors influencing the behavior of buried pipelines subjected to earthquake faulting

Seismic ground faulting is a severe hazard for continuous buried pipelines. Over the years, researchers have attempted to understand pipe behavior, most frequently via numerical modeling and simulation. However, there has been little, if any, physical modeling and tests to verify the numerical modeling approaches and assumptions. This paper presents results of five pairs of centrifuge tests designed to investigate the influence of various factors on the behavior of buried high-density polyethylene (HDPE) pipelines subjected to strike-slip faulting. Parameters considered are the soil moisture content, fault offset rate, relative burial depth (H/D), and pipe diameter. The centrifuge test results show that pipe behavior, specifically pipe strain, is nominally not affected by the soil moisture content and fault offset rate when the pipe is subjected to strike-slip faulting. On the other hand, the burial depth ratio (H/D) and pipe diameter influence peak pipe strain, and in some cases, the ground soil failure pattern.     

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.