地震作用下预应力混凝土管桩运动响应三维数值分析

基于黏弹性人工边界,建立上部结构-桩-土的共同作用三维有限元模型,分析地震作用下预应力混凝土管桩的运动响应特性。分别针对预应力混凝土管桩的桩径、双层软硬土剪切波速比值、上覆土层厚度、上部结构荷载等影响因素进行数值计算。参数分析表明:在地震作用下,桩径的增大会导致桩身整体弯矩相应增加,特别是桩身土层分界面处增大明显;软硬土层剪切波速比及上覆土层厚度的增加,引起土层分界面处桩身峰值弯矩增加;固定桩头条件下,桩头与桩身软硬土层分界面处均会产生较大的运动弯矩;上部结构的惯性荷载对固定桩头的内力有着较大影响,对桩身深处段弯矩影响较小。本文研究结论可为预应力混凝土管桩抗震设计提供有益的理论参考。     

Pile‐head kinematic bending in layered soil

Kinematic effects at the head of a flexible vertical pile embedded in a two‐layer soil deposit are investigated by means of rigorous three‐dimensional elastodynamic finite‐element analyses. Both pile and soil are idealized as linearly viscoelastic materials, modelled by solid elements, without the restrictions associated with the use of strength‐of‐materials approximations. The system is analyzed by a time‐Fourier approach in conjunction with a modal expansion in space. Constant viscous damping is considered for each natural mode, and an FFT algorithm is employed to switch from frequency to time domain and vice versa in natural or generalized coordinates. The scope of the paper is to: (a) elucidate the role of a number of key phenomena controlling the amplitude of kinematic bending moments at the pile head; (b) propose a simplified semi‐analytical formula for evaluating such moments; and (c) provide some remarks about the role of kinematic bending in the seismic design of pile foundations. The results of the study provide a new interpretation of the interplay between interface kinematic moments and corresponding head moments, as a function of layer thickness, pile‐to‐soil stiffness ratio, and stiffness contrast between the soil layers. In addition, the role of diameter in designing against kinematic action, with or without the presence of an inertial counterpart, is discussed. Copyright © 2012 John Wiley & Sons, Ltd.     

Kinematic bending moments in pile foundations

In this paper the kinematic seismic interaction of single piles embedded in soil deposits is evaluated by focusing the attention on the bending moments induced by the transient motion. The analysis is performed by modeling the pile like an Euler–Bernoulli beam embedded in a layered Winkler-type medium. The excitation motion is obtained by means of a one-D propagation analysis. A comprehensive parametric analysis is carried out by varying the main parameters governing the dynamic response of piles like the soil properties, the bedrock location, the diameter and embedment in the bedrock of piles. On the basis of the parametric analysis, a new design formula for predicting the kinematic bending moments for both the cross-sections at the deposit–bedrock interface and at the pile head is proposed.     

A Study of Piles during Earthquakes: Issues of Design and Analysis

The seismic response of pile foundations is a very complex process involving inertial interaction between structure and pile foundation, kinematic interaction between piles and soils, seismically induced pore-water pressures (PWP) and the non-linear response of soils to strong earthquake motions. In contrast, very simple pseudo-static methods are used in engineering practice to determine response parameters for design. These methods neglect several of the factors cited above that can strongly affect pile response. Also soil–pile interaction is modelled using either linear or non-linear springs in a Winkler computational model for pile response. The reliability of this constitutive model has been questioned. In the case of pile groups, the Winkler model for analysis of a single pile is adjusted in various ways by empirical factors to yield a computational model for group response. Can the results of such a simplified analysis be adequate for design in all situations?The lecture will present a critical evaluation of general engineering practice for estimating the response of pile foundations in liquefiable and non-liquefiable soils during earthquakes. The evaluation is part of a major research study on the seismic design of pile foundations sponsored by a Japanese construction company with interests in performance based design and the seismic response of piles in reclaimed land. The evaluation of practice is based on results from field tests, centrifuge tests on model piles and comprehensive non-linear dynamic analyses of pile foundations consisting of both single piles and pile groups. Studies of particular aspects of pile–soil interaction were made. Piles in layered liquefiable soils were analysed in detail as case histories show that these conditions increase the seismic demand on pile foundations. These studies demonstrate the importance of kinematic interaction, usually neglected in simple pseudo-static methods. Recent developments in designing piles to resist lateral spreading of the ground after liquefaction are presented. A comprehensive study of the evaluation of pile cap stiffness coefficients was undertaken and a reliable method of selecting the single value stiffnesses demanded by mainstream commercial structural software was developed. Some other important findings from the study are: the relative effects of inertial and kinematic interactions between foundation and soil on acceleration and displacement spectra of the super-structure; a method for estimating whether inertial interaction is likely to be important or not in a given situation and so when a structure may be treated as a fixed based structure for estimating inertial loads; the occurrence of large kinematic moments when a liquefied layer or naturally occurring soft layer is sandwiched between two hard layers; and the role of rotational stiffness in controlling pile head displacements, especially in liquefiable soils. The lecture concludes with some recommendations for practice that recognize that design, especially preliminary design, will always be based on simplified procedures.     

SOIL–PILE–BRIDGE SEISMIC INTERACTION: KINEMATIC AND INERTIAL EFFECTS. PART I: SOFT SOIL

A substructuring method has been implemented for the seismic analysis of bridge piers founded on vertical piles and pile groups in multi-layered soil. The method reproduces semi-analytically both the kinematic and inertial soil–structure interaction, in a simple realistic way. Vertical S-wave propagation and the pile-to-pile interplay are treated with sufficient rigor, within the realm of equivalent-linear soil behaviour, while a variety of support conditions of the bridge deck on the pier can be studied with the method. Analyses are performed in both frequency and time domains, with the excitation specified at the surface of the outcropping (‘elastic’) rock. A parameter study explores the role of soil–structure interaction by elucidating, for typical bridge piers founded on soft soil, the key phenomena and parameters associated with the interplay between seismic excitation, soil profile, pile–foundation, and superstructure. Results illustrate the potential errors from ignoring: (i) the radiation damping generated from the oscillating piles, and (ii) the rotational component of motion at the head of the single pile or the pile-group cap. Results are obtained for accelerations of bridge deck and foundation points, as well as for bending moments along the piles. © 1997 by John Wiley & Sons, Ltd.     

Static equivalent method for the kinematic interaction analysis of single piles

This paper presents a static equivalent approach to estimate the maximum kinematic interaction effects on piles subjected to lateral seismic excitation. Closed-form expressions are reported for the evaluation of the maximum free-field soil movements and for the computation of maximum pile shear force and bending moments. Firstly, modal analysis, combined with a suitable damped response spectrum, is used to evaluate the maximum free-field response. Secondly, the pile is schematised as a Winkler's beam subjected to equivalent static forces defined according to soil vibration modal shapes and amplitude. The method may be applied by using response spectra suggested by National Standards or those obtained with accelerograms. The procedure proposed may be conveniently implemented in simple spreadsheets or in commercial finite element programs and easily used by practicing engineers. Method accuracy is demonstrated by comparing the results with those obtained with a more rigorous model. Good results may be achieved by considering only the first soil vibration mode making the procedure straightforward for practical design purposes.     

Influences of type of wave and angle of incidence on seismic bending moments in pile foundations

When analysing the seismic response of pile groups, a vertically‐incident wavefield is usually employed even though it does not necessarily correspond to the worst case scenario. This work aims to study the influences of both the type of seismic body wave and its angle of incidence on the dynamic response of pile foundations. To this end, the formulation of SV, SH and P obliquely‐incident waves is presented and implemented in a frequency‐domain boundary element‐finite element code for the dynamic analysis of pile foundations and piled structures. Results are presented in terms of bending moments at cap level of single piles and 3 × 3 pile groups, both in frequency and in time domains. It is found that, in general, the vertical incidence is not the most unfavourable situation. In particular, obliquely‐incident SV waves with angles of incidence smaller than the critical one, a situation in which the mechanism of propagation of the waves in the soil changes and surface waves appear, yield bending moments much larger than those obtained for vertically‐incident wavefields. It is also shown that the influence of pile‐to‐pile interaction on the kinematic bending moments becomes significant for non‐vertical incidence, especially for P and SV waves. Copyright © 2013 John Wiley & Sons, Ltd.     

A model for the 3D kinematic interaction analysis of pile groups in layered soils

The paper presents a numerical model for the analysis of the soil–structure kinematic interaction of single piles and pile groups embedded in layered soil deposits during seismic actions. A finite element model is considered for the pile group and the soil is assumed to be a Winkler‐type medium. The pile–soil–pile interaction and the radiation problem are accounted for by means of elastodynamic Green's functions. Condensation of the problem permits a consistent and straightforward derivation of both the impedance functions and the foundation input motion, which are necessary to perform the inertial soil–structure interaction analyses. The model proposed allows calculating the internal forces induced by soil–pile and pile‐to‐pile interactions. Comparisons with data available in literature are made to study the convergence and validate the model. An application to a realistic pile foundation is given to demonstrate the potential of the model to catch the dynamic behaviour of the soil–foundation system and the stress resultants in each pile. Copyright © 2009 John Wiley & Sons, Ltd.     

Kinematic interaction factors of deep foundations with inclined piles

The beneficial or detrimental role of battered piles on the dynamic response of piled foundations has not been yet fully elucidated. In order to shed more light on this aspect, kinematic interaction factors of deep foundations with inclined piles, are provided for single‐battered piles, as well as for 2 × 2 and 3 × 3 groups of piles subjected to vertically incident plane shear S waves. Piles are modelled as linear‐elastic Bernoulli beams, whereas soil is assumed to be a linear, isotropic, homogeneous viscoelastic half‐space. Different pile group configurations, pile‐soil stiffness ratios, and rake angles are considered. The relevance and main trends observed in the influence of the rake angle on the kinematic interaction factors of the analysed foundations are inferred from the presented results. An important dependence of the kinematic interaction factors on the rake angle is observed together with the existence of an inclination angle at which cap rotation and excitation become out of phase in the low‐to‐mid frequency range. The existence of a small batter angle that provides minimum cap rotation is also shown. Copyright © 2014 John Wiley & Sons, Ltd.     

Kinematic response of single piles to vertically incident P‐waves

Dynamic response of single piles to seismic waves is fundamentally different from the free‐field motion because of the interaction between the pile and the surrounding soil. Considering soil–pile interaction, this paper presents a new displacement model for the steady‐state kinematic response of single piles to vertically incident P‐waves on the basis of a continuum model. The governing equations and boundary conditions of the two undetermined functions in the model are obtained to be coupled by using Hamilton's principle. Then, the two unknown functions are decoupled and solved by an iterative algorithm numerically. A parametric study is performed to investigate the effects of the properties of the soil–pile system on the kinematic response of single piles. It is shown that the effects of the pile–soil modulus ratio, the slenderness ratio of the pile, and the frequency of the incident excitations are very significant. By contrast, the influence of soil damping on the kinematics of the system is slight and can be neglected. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical model for the dynamic analysis of inclined pile groups

The paper presents a numerical model for the dynamic analysis of pile groups with inclined piles in horizontally layered soil deposits. Piles are modelled with Euler–Bernoulli beams, while the soil is supposed to be constituted by independent infinite viscoelastic horizontal layers. The pile–soil–pile interaction as well as the hysteretic and geometric damping is taken into account by means of two‐dimensional elastodynamic Green's functions. Piles cap is considered by introducing a rigid constraint; the condensation of the problem permits a consistent derivation of both the dynamic impedance matrix of the soil–foundation system and the foundation input motion. These quantities are those used to perform inertial soil–structure interaction analyses in the framework of the substructure approach. Furthermore, the model allows evaluating the kinematic stress resultants in piles resulting from waves propagating in the soil deposit, taking into account the pile–soil–pile interactions. The model validation is carried out by performing accuracy analyses and comparing results in terms of dynamic impedance functions, kinematic response parameters and pile stress resultants, with those furnished by 3D refined finite element models. To this purpose, classical elastodynamic solutions are adopted to define the soil–pile interaction problem. The model results in low computational demands without significant loss of precision, compared with more rigorous approaches or refined finite element models. Copyright © 2015 John Wiley & Sons, Ltd.     

Kinematic internal forces in deep foundations with inclined piles

This paper presents a parametric study that looks into the influence of pile rake angle on the kinematic internal forces of deep foundations with inclined piles. Envelopes of maximum kinematic bending moments, shear forces and axial loads are presented along single inclined piles and 2 × 2 symmetrical square pile groups with inclined elements subjected to an earthquake generated by vertically incident shear waves. Inclination angles from 0° to 30° are considered, and three different pile–soil stiffness ratios are studied. These results are obtained through a frequency–domain analysis using a boundary element–finite element code in which the soil is modelled by the boundary element method as a homogeneous, viscoelastic, unbounded region, and the piles are modelled by finite elements as Euler–Bernoulli beams. The rotational kinematic response of the pile foundations is shown to be a key factor on the evolution of the kinematic internal forces along the foundations. Copyright © 2015 John Wiley & Sons, Ltd.     

Artificial neural network application to estimate kinematic soil pile interaction response parameters

Six artificial neural network (ANN) models are developed to predict various response parameters of kinematic soil pile interaction. These responses include (1) kinematic response factors for free and fixed head piles in homogenous soil layer to derive foundation input motion (2) normalized bending moment at fixed head of pile in homogenous soil layer (3) normalized kinematic pile moment at the interface of two soil layers of sharply different soil stiffnesses. These ANN models represent simple solutions that can be implemented in a simple calculator capable of matrix operation and bypass the site response analysis and the complex wave diffraction analysis. The data required for ANN training is generated using beam on dynamic Winkler formulation (BDWF). Fifty percent of the data is used to train the ANN models while remaining 50% is used to test the ANN models. The trained ANN models show good agreement with BDWF results.     

Experimental study of the effect of inclined pile on the seismic behavior of pile group

Seismic behavior of inclined piles has been considered detrimental for years. However, recent researches show that battered piles can have a beneficial effect. In this framework, a series of centrifuge tests on an inclined pile group is performed. The analysis is based on the comparative response of two 2×1 simplified pile groups: one with vertical piles and the other with one vertical and one inclined pile. The response of these pile groups to repeated earthquakes or sinusoidal inputs is analyzed through the response frequencies, the envelop curves of bending moment profiles, the axial loads measured in both piles and the kinematic response of the cap. Results highlight that the effect of inclined pile is highly influenced by the frequency content of the input. In addition, the inclined pile induces non-negligible residual bending moments, higher horizontal stiffness at the pile cap and larger rotation.     

Seismic soil–structure interaction in multi‐span bridges: Application to a railway bridge

The effects of soil‐structure interaction on the seismic response of multi‐span bridges are investigated by means of a modelling strategy based on the domain decomposition technique. First, the analysis methodology is presented: kinematic interaction analysis is performed in the frequency domain by means of a procedure accounting for radiation damping, soil–pile and pile‐to‐pile interaction; the seismic response of the superstructure is evaluated in the time domain by means of user‐friendly finite element programs introducing suitable lumped parameter models take into account the frequency‐dependent impedances of the soil–foundation system. Second, a real multi‐span railway bridge longitudinally restrained at one abutment is analyzed. The input motion is represented by two sets of real accelerograms: one consistent with the Italian seismic code and the other constituted by five records characterized by different frequency contents. The seismic response of the compliant‐base model is compared with that obtained from a fixed‐base model. Pile stress resultants due to kinematic and inertial interactions are also evaluated. The application demonstrates the importance of performing a comprehensive analysis of the soil–foundation–structure system in the design process, in order to capture the effects of soil‐structure interaction in each structural element that may be beneficial or detrimental. Copyright © 2011 John Wiley & Sons, Ltd.     

1st Ishihara Lecture: An overview of the behavior of pile foundations in liquefiable and non-liquefiable soils during earthquake excitation

The seismic response of a pile foundation is usually analyzed by approximate methods in practice. These methods typically neglect one or more of the important factors that affect seismic response such as inertial interaction, kinematic interaction, seismic pore water pressures, soil nonlinearity, cross stiffness coupling and dynamic pile to pile interaction. A nonlinear 3-D analysis is used to show how all these factors affect pile response, to demonstrate some of the consequences of using various approximate methods and to provide a comprehensive overview of how pile foundations behave during earthquakes in liquefiable and non-liquefiable soils.     

Soil–pile–structure interaction: experimental outcomes from shaking table tests

An effective way to study the complex seismic soil‐structure interaction phenomena is to investigate the response of physical scaled models in 1‐g or n‐g laboratory devices. The outcomes of an extensive experimental campaign carried out on scaled models by means of the shaking table of the Bristol Laboratory for Advanced Dynamics Engineering, University of Bristol, UK, are discussed in the present paper. The experimental model comprises an oscillator connected to a single or a group of piles embedded in a bi‐layer deposit. Different pile head conditions, that is free head and fixed head, several dynamic properties of the structure, including different masses at the top of the single degree of freedom system, excited by various input motions, e.g. white noise, sinedwells and natural earthquake strong motions recorded in Italy, have been tested. In the present work, the modal dynamic response of the soil–pile–structure system is assessed in terms of period elongation and system damping ratio. Furthermore, the effects of oscillator mass and pile head conditions on soil–pile response have been highlighted, when the harmonic input motions are considered. Copyright © 2015 John Wiley & Sons, Ltd.     

Three-dimensional nonlinear analysis for seismic soil–pile-structure interaction

A three-dimensional method of analysis is presented for the seismic response of structures constructed on pile foundations. An analysis is formulated in the time domain and the effects of material nonlinearity of soil on the seismic response are investigated. A subsystem model consisting of a structure subsystem and a pile-foundation subsystem is used. Seismic response of the system is found using a successive-coupling incremental solution scheme. Both subsystems are assumed to be coupled at each time step. Material nonlinearity is accounted for by incorporating an advanced plasticity-based soil model, HiSS, in the finite element formulation. Both single piles and pile groups are considered and the effects of kinematic and inertial interaction on seismic response are investigated while considering harmonic and transient excitations. It is seen that nonlinearity significantly affects seismic response of pile foundations as well as that of structures. Effects of nonlinearity on response are dependent on the frequency of excitation with nonlinearity causing an increase in response at low frequencies of excitation.     

Transient kinematic pile bending in two-layer soil

The dynamic response of piles to seismic loading is explored by means of an extensive parametric study based on a properly calibrated Beam-on-Dynamic-Winkler-Foundation (BDWF) model. The investigated problem consists of a single vertical cylindrical pile, modelled as an Euler–Bernoulli beam, embedded in a subsoil consisting of two homogeneous viscoelastic layers of sharply different stiffness resting on a rigid stratum. The system is subjected to vertically propagating seismic S waves, in the form of a transient motion imposed on rock outcrop. Several accelerograms recorded in Italy are employed as input motions in the numerical analyses. The paper highlights the severity of kinematic pile bending in the vicinity of the interface separating the two soil layers. In addition to factors already investigated such as layer stiffness contrast, relative soil–pile stiffness, interface depth and intensity of ground excitation, the paper focuses on additional important factors, notably soil material damping, stiffness of Winkler springs and frequency content of earthquake excitation. Existing predictive equations for assessing kinematic pile bending at soil layer interfaces are revisited and new regression analyses are performed. A synthesis of findings in terms of a set of simple equations is provided. The use of these equations is discussed through examples.     

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

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