Nonlinear response of laterally loaded piles and pile groups

In spite of extensive studies on laterally loaded piles carried out over years, none of them offers an expedite approach as to gaining the nonlinear response and its associated depth of mobilization of limiting force along each pile in a group. To serve such a need, elastic–plastic solutions for free‐head, laterally loaded piles were developed recently by the author. They allow the response to be readily computed from elastic state right up to failure, by assigning a series of slip depths, and a limiting force profile. In this paper, equivalent solutions for fixed‐head (FixH) single piles were developed. They are subsequently extended to cater for response of pile groups by incorporating p‐multipliers. The newly established solutions were substantiated by existing numerical solutions for piles and pile groups. They offer satisfactory prediction of the nonlinear response of all the 6 single piles and 24 pile groups investigated so far after properly considering the impact of semi‐FixH restraints. They also offer the extent to ultimate state of pile groups via the evaluated slip depths. The study allows ad hoc guidelines to be established for determining input parameters for the solutions. The solutions are tailored for routine prediction of the nonlinear interaction of laterally loaded FixH piles and capped pile groups. Copyright © 2008 John Wiley & Sons, Ltd.     

Elastic models for nonlinear response of rigid passive piles

Recent study indicates that the response of rigid passive piles is dominated by elastic pile–soil interaction and may be estimated using theory for lateral piles. The difference lies in that passive piles normally are associated with a large scatter of the ratio of maximum bending moment over maximum shear force and induce a limiting pressure that is ~1/3 that on laterally loaded piles. This disparity prompts this study. This paper proposes pressure‐based pile–soil models and develops their associated solutions to capture response of rigid piles subjected to soil movement. The impact of soil movement was encapsulated into a power‐law distributed loading over a sliding depth, and load transfer model was adopted to mimic the pile–soil interaction. The solutions are presented in explicit expressions and can be readily obtained. They are capable of capturing responses of model piles in a sliding soil owing to the impact of sliding depth and relative strength between sliding and stable layer on limiting force prior to ultimate state. In comparison with available solutions for ultimate state, this study reveals the 1/3 limiting pressure (of the active piles) on passive piles was induced by elastic interaction. The current models employing distributed pressure for moving soil are more pertinent to passive piles (rather than plastic soil flow). An example calculation against instrumented model piles is provided, which demonstrates the accuracy of the current solutions for design slope stabilising piles. Copyright © 2014 John Wiley & Sons, Ltd.     

Response of rigid piles during passive dragging

This paper develops a three‐layer model and elastic solutions to capture nonlinear response of rigid, passive piles in sliding soil. Elastic solutions are obtained for an equivalent force per unit length p_s of the soil movement. They are repeated for a series of linearly increasing p_s (with depth) to yield the nonlinear response. The parameters underpinning the model are determined against pertinent numerical solutions and model tests on passive free‐head and capped piles. The solutions are presented in non‐dimensional charts and elaborated through three examples. The study reveals the following:

• On‐pile pressure in rotationally restrained, sliding layer reduces by a factor α, which resembles the p‐multiplier for a laterally loaded, capped pile, but for its increase with vertical loading (embankment surcharge), and stiffness of underlying stiff layer: α = 0.25 and 0.6 for a shallow, translating and rotating piles, respectively; α = 0.33–0.5 and 0.8–1.3 for a slide overlying a stiff layer concerning a uniform and a linearly increasing pressure, respectively; and α = 0.5–0.72 for moving clay under embankment loading.
• Ultimate state is well defined using the ratio of passive earth pressure coefficient over that of active earth pressure. The subgrade modulus for a large soil movement may be scaled from model tests.
• The normalised rotational stiffness is equal to 0.1–0.15 for the capped piles, which increases the pile displacement with depth.

The three‐layer model solutions well predict nonlinear response of capped piles subjected to passive loading, which may be used for pertinent design. Copyright © 2016 John Wiley & Sons, Ltd.     

Visco‐elastic load transfer models for axially loaded piles

Viscoelastic or creep behaviour can have a significant influence on the load transfer (t–z) response at the pile–soil interface, and thus on the pile load settlement relationship. Many experimental and theoretical models for pile load transfer behaviour have been presented. However, none of these has led to a closed‐form expression which captures both non‐linearity and viscoelastic behaviour of the soil. In this paper, non‐linear viscoelastic shaft and base load transfer (t–z) models are presented, based on integration of a generalized viscoelastic stress–strain model for the soil. The resulting shaft model is verified through published field and laboratory test data. With these models, the previous closed‐form solutions evolved for a pile in a non‐homogeneous media have been readily extended to account for visco‐elastic response. For 1‐step loading case, the closed‐form predictions have been verified extensively with previous more rigorous numerical analysis, and with the new GASPILE program analysis. Parametric studies on two kinds of commonly encountered loading: step loading, ramp (linear increase followed by sustained) loading have been performed. Two examples of the prediction of the effects of creep on the load settlement relationship by the solutions and the program GASPILE, have been presented. Copyright © 2000 John Wiley & Sons, Ltd.     

A Winkler model approach for vertically and laterally loaded piles in nonhomogeneous soil

An investigation is made to present analytical solutions provided by a Winkler model approach for the analysis of single piles and pile groups subjected to vertical and lateral loads in nonhomogeneous soils. The load transfer parameter of a single pile in nonhomogeneous soils is derived from the displacement influence factor obtained from Mindlin's solution for an elastic continuum analysis, without using the conventional form of the load transfer parameter adopting the maximum radius of the influence of the pile proposed by Randolph and Wroth. The modulus of the subgrade reaction along the pile in nonhomogeneous soils is expressed by using the displacement influence factor related to Mindlin's equation for an elastic continuum analysis to combine the elastic continuum approach with the subgrade reaction approach. The relationship between settlement and vertical load for a single pile in nonhomogeneous soils is obtained by using the recurrence equation for each layer. Using the modulus of the subgrade reaction represented by the displacement influence factor related to Mindlin's solution for the lateral load, the relationship between horizontal displacement, rotation, moment, and shear force for a single pile subjected to lateral loads in nonhomogeneous soils is available in the form of the recurrence equation. The comparison of the results calculated by the present method for single piles and pile groups in nonhomogeneous soils has shown good agreement with those obtained from the more rigorous finite element and boundary element methods. It is found that the present procedure gives a good prediction on the behavior of piles in nonhomogeneous soils. Copyright © 2011 John Wiley & Sons, Ltd.     

Embedded beam element with interaction surface for lateral loading of piles

This paper presents a numerical formulation of a three dimensional embedded beam element for the modeling of piles, which incorporates an explicit interaction surface between soil and pile. The formulation is herein implemented for lateral loading of piles but is able to represent soil–pile interaction phenomena in a general manner for different types of loading conditions or ground movements. The model assumes perfect adherence between beam and soil along the interaction surface. The paper presents a comparison of the results obtained by means of the present formulation and by means of a previously formulated embedded pile element without interaction surface, as well as reference semi‐analytical solutions and a fully 3D finite element (FE) model. It is seen that the proposed embedded element provides a better convergence behavior than a previously formulated embedded element and is able to reproduce key features of a full 3D FE model. Copyright © 2015 John Wiley & Sons, Ltd.     

Nonlinear analysis of laterally loaded rigid piles in cohesive soil

This article presents a method for the nonlinear analysis of laterally loaded rigid piles in cohesive soil. The method considers the force and the moment equilibrium to derive the system equations for a rigid pile under a lateral eccentric load. The system equations are then solved using an iteration scheme to obtain the response of the pile. The method considers the nonlinear variation of the ultimate lateral soil resistance with depth and uses a new closed‐form expression proposed in this article to determine the lateral bearing factor. The method also considers the horizontal shear resistance at the pile base, and a bilinear relationship between the shear resistance and the displacement is used. For simplicity, the modulus of horizontal subgrade reaction is assumed to be constant with depth, which is applicable to piles in overconsolidated clay. The nonlinearity of the modulus of horizontal subgrade reaction with pile displacement at ground surface is also considered. The validity of the developed method is demonstrated by comparing its results with those of 3D finite element analysis. The applications of the developed method to analyze five field test piles also show good agreement between the predictions and the experimental results. The developed method offers an alternative approach for simple and effective analysis of laterally loaded rigid piles in cohesive soil. Copyright © 2011 John Wiley & Sons, Ltd.     

A three‐dimensional displacement approach for analysis of laterally loaded piles in nonhomogeneous soil

An analytical approach using the three‐dimensional displacement of a soil is investigated to provide analytical solutions of the horizontal response of a circular pile subjected to lateral loads in nonhomogeneous soil. The rocking stiffness coefficient of the pile shaft in homogeneous soil is derived from the analytical solution taking into account the three‐dimensional displacement represented in terms of scalar potentials in the elastic three‐dimensional analysis. The lateral stiffness coefficient of the pile shaft in nonhomogeneous soil is derived from the rocking stiffness coefficient taking into account the rocking rotation of a rigid pile shaft. The relationship between horizontal displacement, rotation, moment, and shear force of a pile subjected to horizontal loads in nonhomogeneous soil is obtainable in the form of the recurrence equation. The formulation of the lateral displacement and rotation of the pile base subjected to lateral loads in nonhomogeneous soils is presented by taking into account Mindlin's equation and the equivalent thickness for soil layers in the equivalent elastic method. There is little difference between lateral, rocking, and couple stiffness coefficients each obtained from both the two‐dimensional and three‐dimensional methods except for the case of Poisson's ratio near 0.5. The comparison of results calculated by the current method for a pile subjected to lateral loads in homogeneous and nonhomogeneous soils has shown good agreement with those obtained from analytical and numerical methods. Copyright © 2017 John Wiley & Sons, Ltd.     

Vertically loaded piles in non-homogeneous media

Analytical methods for the axial responses of piles can be classified under three broad categories of (1) simple but approximate analytical solutions, (2) one-dimensional numerical algorithms, (3) full axisymmetric analyses using boundary or finite element approaches. The first two categories rely on the so-called load transfer approach, with interaction between pile and soil determined by independent springs distributed along the pile shaft and at the pile base. The non-linear spring stiffness is related to the elastic–plastic properties of the actual soil partly by empirically based correlations and partly by theoretical arguments based on simplified models of the pile–soil system. This paper presents new closed-form solutions for the axial response of piles in elastic–plastic, non-homogeneous, media. The solutions fall in the first of the three categories above, and have been verified through extensive parametric studies using more rigorous one-dimensional and continuum analyses. The effect of non-homogeneity and partial slip on the load and displacement profiles along the pile shaft is explored, and comparisons are presented with experimental data. © 1997 John Wiley & Sons, Ltd.     

水平受荷长桩弹塑性计算解析解

当考虑桩侧土体非线性本构关系时对水平受荷桩的计算一般需采用数值方法,解析结果相对较少。基于Winkler地基模型和桩侧土体简化的弹塑性本构关系,对均质地基中水平荷载作用下桩头嵌固的长桩进行了解析推导,得到了桩身最大挠度及最大弯矩与荷载关系的统一解析表达式,并采用相同的方法求得高桩情形下桩头挠度的计算式。计算表明,联合荷载作用下桩身泥面处的挠度和转角不等于单个荷载作用时的线性叠加,采用常规的线性叠加法计算将偏于不安全。所求解析式借助计算器即可进行最大挠度和最大弯矩的计算,大大方便了工程的计算应用。     

Axial kinematic response of end‐bearing piles to P waves

Kinematic pile–soil interaction under vertically impinging seismic P waves is revisited through a novel continuum elastodynamic solution of the Tajimi type. The proposed model simulates the steady‐state kinematic response of a cylindrical end‐bearing pile embedded in a homogeneous viscoelastic soil stratum over a rigid base, subjected to vertically propagating harmonic compressional waves. Closed‐form solutions are obtained for the following: (i) the displacement field in the soil and along the pile; (ii) the kinematic Winkler moduli (i.e., distributed springs and dashpots) along the pile; (iii) equivalent, depth‐independent, Winkler moduli to match the motion at the pile head. The solution for displacements is expressed in terms of dimensionless transfer functions relating the motion of the pile head to the free‐field surface motion and the rock motion. It is shown that (i) a pile foundation may significantly alter (possibly amplify) the vertical seismic excitation transmitted to the base of a structure and (ii) Winkler moduli pertaining to kinematic loading differ from those for inertial loading. Simple approximate expressions for kinematic Winkler moduli are derived for use in applications. Copyright © 2013 John Wiley & Sons, Ltd.     

Nonlinear analysis of laterally loaded rigid piles in cohesionless soil

In this paper, a method is developed for nonlinear analysis of laterally loaded rigid piles in cohesionless soil. The method assumes that both the ultimate soil resistance and the modulus of horizontal subgrade reaction increase linearly with depth. By considering the force and moment equilibrium, the system equations are derived for a rigid pile under a lateral eccentric load. An iteration scheme containing three main steps is then proposed to solve the system equations to obtain the response of the pile. To determine the ultimate soil resistance and the modulus of horizontal subgrade reaction required in the analysis, related expressions are selected by reviewing and assessing the existing methods. The degradation of the modulus of horizontal subgrade reaction with pile displacement at ground surface is also considered. The developed method is validated by comparing its results with those of centrifugal tests and three-dimensional finite element analysis. Applications of the developed method to laboratory model and field test piles also show good agreement between the predictions and the experimental results.     

分布荷载推力桩计算的p-y曲线法研究

为提高现有水平梯形分布荷载推力桩设计计算水平,提出了地基系数按非线性土抗力.水平位移(p-y)曲线抗力模式表达的水平梯形分布荷载推力桩位移和内力计算的有限差分数值分析方法,并详细推导了桩身位移的差分格式.基于这些公式,编制了全桩位移、内力及侧土抗力的计算和图形处理程序,可适用于滑坡抗滑桩和深基坑悬臂支护桩的设计计算.结合新型桩型现浇混凝土薄壁管桩(PCC桩)作为滑坡抗滑桩进行对比分析,算例表明:该方法方便可靠.当有限差分段划分得足够小时,可使数值解接近于真实解.     

Selection of p–y curves for the design of single laterally loaded piles

Two-dimensional finite element analysis has been used to find load–transfer relationships for translation of an infinitely long pile through undrained soil for a variety of soil-constitutive models. It has been shown that these load–transfer curves can be used as p–y curves in the analysis of single piles undergoing lateral pile head loading in undrained soils with non-linear stress–strain laws. Lateral pile response deduced from 2-D analysis input to the subgrade reaction method has been compared to the behaviour of a single pile analysed using three-dimensional finite element analysis. Good agreement between the two methods for non-linear soils suggests that the 2-D analysis may form a useful design method for calculation of p–y curves. Copyright © 1999 John Wiley & Sons, Ltd.     

Performance of seismically loaded shearwalls on nonlinear shallow foundations

If utilized, the energy dissipative capability of seismically loaded shallow foundations due to inelastic behavior can result in more economic design, provided the consequences, such as excessive deformations are accounted for. In this article, a Beam‐on‐Nonlinear‐Winkler‐Foundation (BNWF) model is used to assess the performance of shearwall‐foundation systems with different attributes, when subjected to ground motions of varied hazard levels. The numerical study indicates that the force and drift demands to the shearwall reduce significantly, when nonlinear foundation behavior is realized, while permanent settlement is well below the permissible limit. These results support the concept of shallow foundation capacity mobilization in future design. Copyright © 2010 John Wiley & Sons, Ltd.     

Soil resistances for laterally loaded rigid piles in multi-layered elastic soil

ABSTRACT

Short stubby piles like monopiles and large diameter drilled shafts undergo rigid body translation and rotation when subjected to a lateral force and/or a moment at the head. A method of analysis for these piles embedded in multi-layered elastic soil is developed using the variational principles of mechanics. Using this analysis, the soil resistance against pile movement can be rigorously related to the soil elastic constants, and the pile head displacement and rotation can be quickly calculated. The equilibrium equations for pile and soil displacements are obtained using the principle of virtual work and solved using an iterative algorithm. Pile responses obtained from the analysis match well with those obtained from three-dimensional finite element analyses in which the same inputs of loads, geometry, and material properties are given. Based on the new analysis, fitted equations for soil resistance parameters are developed, which can be used to directly calculate the pile head displacement and rotation without the use of the iterative algorithm. Numerical examples are provided that demonstrate how the method can be used to analyse practical problems.     

A modulus‐multiplier approach for non‐linear analysis of laterally loaded pile groups

A modulus‐multiplier approach, which applies a reduction factor to the modulus of single pile p–y curves to account for the group effect, is presented for analysing the response of each individual pile in a laterally loaded pile group with any geometric arrangement based on non‐linear pile–soil–pile interaction. The pile–soil–pile interaction is conducted using a 3D non‐linear finite element approach. The interaction effect between piles under various loading directions is investigated in this paper. Group effects can be neglected at a pile spacing of 9 times the pile diameter for piles along the direction of the lateral load and at a pile spacing of 6 times the pile diameter for piles normal to the direction of loading. The modulus multipliers for a pair of piles are developed as a function of pile spacing for departure angle of 0, 90, and 180sup>/sup> with respect to the loading direction. The procedure proposed for computing the response of any individual pile within a pile group is verified using two well‐documented full‐scale pile load tests. Copyright © 2006 John Wiley & Sons, Ltd.     

Transverse seismic kinematics of single piles by a modified Vlasov model

Within the framework of soil–pile interaction, a novel displacement scheme for the transverse kinematic response of single piles to vertically propagating S waves is proposed on the basis of the modified Vlasov foundation model. The displacement model contains a displacement function along the pile axis and an attenuation function along the radial direction. The governing equations and boundary conditions of the two undetermined functions are obtained in a coupled form by using Hamilton's principle. An iterative algorithm is adopted to decouple and solve the two unknown functions. In light of the governing equation of the pile kinematics, a mechanical model is proposed to evaluate the present method on a physical basis considering material damping. The coefficient of the equivalent Winkler spring is derived explicitly as function of the displacement decay parameter γ and soil Poisson's ratio. A parametric study is performed to investigate the effects of the soil–pile system properties on the kinematic response of single piles. The results show that the dimensionless pile length controls the transverse kinematics of piles. In terms of the theory of beams on elastic foundation, the classification limits of the dimensionless pile length may be π ∕ 4 and π, respectively. Copyright © 2014 John Wiley & Sons, Ltd.     

Expanded superposition method for impedance functions of inclined‐pile groups

This paper presents a superposition method expanded for computing impedance functions (IFs) of inclined‐pile groups. Closed‐form solutions for obtaining horizontal, vertical, and rocking IFs, estimated by using pile‐to‐pile interaction factors, are proposed. IFs of solitary inclined piles, crossed IFs, and explicit incorporation of compatibility conditions for pile‐head movements are also appropriately taken into consideration. All of these factors should be known in advance and will be computed and shown for the most relevant cases. The accuracy of the proposed closed‐form solutions is verified for 2 × 2 and 3 × 3 square inclined‐pile groups embedded in an isotropic viscoelastic homogeneous half‐space soil medium, with hysteretic damping. The pile‐to‐pile interaction factors are computed by means of a three‐dimensional time‐harmonic boundary elements–finite elements coupling formulation. The results indicate that the IFs obtained from the proposed method are in good agreement with those obtained from the coupling formulation. Furthermore, crossed vertical‐rocking IFs of solitary piles need to be appropriately considered for obtaining rocking IFs when the number of piles is small. Copyright © 2015 John Wiley & Sons, Ltd.     

Beam–solid contact formulation for finite element analysis of pile–soil interaction with arbitrary discretization

This paper presents an embedded beam formulation for discretization independent finite element (FE) analyses of interactions between pile foundations or rock anchors and the surrounding soil in geotechnical and tunneling engineering. Piles are represented by means of finite beam elements embedded within FEs for the soil represented by 3D solid elements. The proposed formulation allows consideration of piles and pile groups with arbitrary orientation independently from the FE discretization of the surrounding soil. The interface behavior between piles and the surrounding soil is represented numerically by means of a contact formulation considering skin friction as well as pile tip resistance. The pile–soil interaction along the pile skin is considered by means of a 3D frictional point‐to‐point contact formulation using the integration points of the beam elements and reference points arbitrarily located within the solid elements as control points. The ability of the proposed embedded pile model to represent groups of piles objected to combined axial and shear loading and their interactions with the surrounding soil is demonstrated by selected benchmark examples. The pile model is applied to the numerical simulation of shield driven tunnel construction in the vicinity of an existing building resting upon pile foundation to demonstrate the performance of the proposed model in complex simulation environments. Copyright © 2014 John Wiley & Sons, Ltd.