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.     

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.     

Analysis of lateral loading of pile groups using embedded beam elements with interaction surface

The numerical simulation of soil-pile interaction problems, by means of full 3D finite element models, involves a large number of degrees of freedom (DOF) and difficulties during the mesh generation process. In order to reduce the unknowns and simplify and properly analyze such class of geotechnical problems, the so-called embedded beam elements (EBE) have recently been developed. In a preceding contribution of the authors, an improved EBE formulation, which brings into play the soil-pile interaction surface, was proposed with the aim to localize material plasticity in the soil surrounding the pile. This embedded beam model couples two different finite elements, each described by distinct kinematics (ie, solid and beam). The coupling is incorporated in the formulation by means of kinematical constrains established over the solid and beam displacement fields on the interaction surface. One of the main advantages of the embedded elements is that the addition of beams structural members immersed within the 3D soil model does not represent a constraint for the solid mesh, which can be adopted independently from the beam mesh. In this paper, the lateral loading of pile groups is studied by means of the proposed EBE approach with elasto-plastic interfaces. In order to represent a rigid cap, a master node and a special set of kinematical restrictions are incorporated into the formulation. The paper presents results obtained by means of the present formulation compared against other well-established analysis methods and test results published in the literature, for both elastic and elasto-plastic cases.     

A non‐linear coupled finite element–boundary element model for the prediction of vibrations due to vibratory and impact pile driving

This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Interaction of a cylindrical thin‐walled pile with an inhomogeneous isotropic elastic medium

This paper presents a rigorous analysis for the static interaction of a cylindrical thin‐walled pile with an inhomogeneous isotropic elastic half‐space under vertical, horizontal, and torsional forces individually applied at the top of pile. The inhomogeneity is specified with the exponential variation of shear modulus along depth of the embedding medium, and the Poisson's ratio is assumed to be constant. By means of a set of Green's functions for pile and soil medium and satisfying the compatibility conditions between the 2 interacting media, the formulation is reduced to coupled Fredholm integral equations. Using the adaptive‐gradient elements, capable of capturing the singular stress transfer at both ends of the pile, a numerical procedure is developed and utilized for evaluating the relevant integral equations and studying the inhomogeneity effect on the soil‐pile interaction responses. The analysis results have been validated for different soil‐pile modulus ratios under axial load and for a Poisson's ratio of 0.3 under lateral load. The procedure does not consider the nonlinear behavior of the soil medium or plastic yielding in the pile section, and the impact of the unreliable results for the case of high Poisson's ratio is not examined.     

A revised solution for the horizontal vibration of an end‐bearing pile in viscoelastic soil

This note presents a new method to derive closed‐form expressions describing the horizontal response of an end‐bearing pile in viscoelastic soil subjected to harmonic loads at its head. The soil surrounding the pile is assumed as a linearly viscoelastic layer. The propagation of waves in the soil and pile is treated mathematically by three‐dimensional and one‐dimensional theories, respectively. Unlike previous studies of the problem, the formulation presented allows the governing equations of the soil to be solved directly, eliminating the need to introduce potential functions. Accordingly, the dynamic response of the pile is obtained by means of the initial parameter method, invoking the requirement for continuity at the pile–soil interface. It is demonstrated that the derived compact solution matches exactly an existing solution that utilises potential functions to formulate the problem. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of piles subjected to lateral soil movements using a three‐dimensional displacement approach

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 soil movements in nonhomogeneous soil. The lateral stiffness coefficient of the pile shaft in nonhomogeneous soil is derived from the rocking stiffness coefficient that is obtained from the analytical solution, taking into account the three‐dimensional displacement represented in terms of scalar potentials in the elastic three‐dimensional analysis. The relationship between horizontal displacement, rotation, moment, and shear force of a pile subjected to lateral soil movements in nonhomogeneous soil is obtainable in the form of the recurrence equation. For the relationship between the lateral pressure and the horizontal displacement, it is assumed that the behavior is linear elastic up to lateral soil yield, and the lateral pressure is constant under the lateral soil yield. The interaction factors between piles subjected to both lateral load and moment are calculated, taking into account the lateral soil movement. 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 the Mindlin equation and the equivalent thickness for soil layers in the equivalent elastic method. For lateral movement, lateral pressure, bending moment, and interaction factors, there are small differences between results obtained from the 1‐D and the 3‐D displacement methods except a very flexible pile. Copyright © 2015 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.     

Numerical analysis of the response of battered piles to inclined pullout loads

This paper presents a three‐dimensional finite element analysis of the response of battered piles to the combined lateral and vertical pullout loads. Analyses are carried out using an elastoplastic constitutive law based on the non‐associated Mohr–Coulomb criterion. The influence of the contact condition at the pile–soil interface is also investigated. Analyses show that the load's inclination with regard to the pile's axis affects both the lateral and axial response of the battered piles. Analyses also show that the pullout capacity of battered piles is affected by the pile's inclination regarding the vertical axis as well as the load's inclination regarding the pile's axis. The investigation of the influence of the contact condition at the soil–pile interface shows that the possibility of sliding at the soil–pile interface affects the response of battered piles subjected to loads with low inclination regarding the pile's axis. Copyright © 2008 John Wiley & Sons, Ltd.     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

Improved numerical algorithms for frictional contact in pile penetration analysis

This paper presents a numerical formulation for frictional contact problems associated with pile penetration. The frictional contact at the soil–pile interface is formulated using the theory of hardening/softening plasticity, so that advanced models for the interface can be dealt with. A smooth discretisation of the pile surface is proposed using Bézier polynomials. An automatic load stepping scheme is proposed, which features an error control algorithm and automatic subincrementation of the load increments. The numerical algorithms are then used to analyse the installation process of pushed-in axial piles. It is shown that the smooth discretisation of the pile surface is effective in reducing the oscillation in the predicted pile resistances and the automatic load stepping scheme outperforms the classical Newton–Raphson scheme for this type of problem.     

Analysis of pile-soil interaction under lateral loading using infinite and finite elements

In order to gain a better understanding of pile-soil interaction under lateral loading, this paper presents a numerical analysis which combines the infinite and finite element method. Interest is focused on the group effect on ultimate lateral soil resistance. Firstly, a single isolated pile is analysed and reasonably good agreement is found between existing analytical solutions and results obtained by the present method. A limited parametric study is also presented and some parameters influencing the ultimate lateral soil resistance are identified. The analysis of pile groups is then considered and it is shown that the group effect tends to reduce pile capacity when the spacings between piles are within the practical ranges. The extent of the reduction depends on the arrangement of piles within the group.     

Cyclic macro‐element for soil–structure interaction: material and geometrical non‐linearities

This paper presents a non‐linear soil–structure interaction (SSI) macro‐element for shallow foundation on cohesive soil. The element describes the behaviour in the near field of the foundation under cyclic loading, reproducing the material non‐linearities of the soil under the foundation (yielding) as well as the geometrical non‐linearities (uplift) at the soil–structure interface. The overall behaviour in the soil and at the interface is reduced to its action on the foundation. The macro‐element consists of a non‐linear joint element, expressed in generalised variables, i.e. in forces applied to the foundation and in the corresponding displacements. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. Mechanisms of yielding and uplift are modelled through a global, coupled plasticity–uplift model. The cyclic model is dedicated to modelling the dynamic response of structures subjected to seismic action. Thus, it is especially suited to combined loading developed during this kind of motion. Comparisons of cyclic results obtained from the macro‐element and from a FE modelization are shown in order to demonstrate the relevance of the proposed model and its predictive ability. Copyright © 2001 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.     

Three‐dimensional finite element analyses of passive pile behaviour

Piles may be subjected to lateral soil pressures as a result of lateral soil movements from nearby construction‐related activities such as embankment construction or excavation operations. Three‐dimensional finite element analyses have been carried out to investigate the response of a single pile when subjected to lateral soil movements. The pile and the soil were modelled using 20‐node quadrilateral brick elements with reduced integration. For compatibility between the soil–pile interface elements, 27‐node quadrilateral brick elements with reduced integration were used to model the soil around the pile adjacent to the soil–pile interface. A Mohr–Coulomb elastic–plastic constitutive model with large‐strain mode was assumed for the soil. The analyses indicate that the behaviour of the pile was significantly influenced by the pile flexibility, the magnitude of soil movement, the pile head boundary conditions, the shape of the soil movement profile and the thickness of the moving soil mass. Reasonable agreement is found between some existing published solutions and those developed herein. Copyright © 2005 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.     

碎石土土体性状对桩基水平承载能力的影响

桩基作为一种广泛使用的基础形式,其在如地震荷载、风荷载下的水平承载力确定十分重要。桩基水平承载能力主要由桩身及土体性质控制,而土体性质中密实度和含石量是影响碎石土场地上桩基水平承载能力的主要因素。为明确碎石土场地土性对桩基水平承载能力的影响,在室内开展5组不同土体密实度和含石量下的水平场地单桩水平静载试验,以此模型试验数据计算所得的临界荷载、极限荷载、水平扩散角、计算宽度和水平抗力系数的比例系数m值等为指标评价土体密实度和含石量对桩基水平承载能力的影响。试验结果表明增大密实度和含石量对各指标均有提升效果,但相比之下,密实度对各项指标的提升效果均明显优于含石量。综合考量下,碎石土密实度和含石量的增大均对桩基水平承载力有提升作用,但土体密实度的提升作用更明显。故对水平承载力要求较高的碎石土上桩基工程的选址中,应优先选择密实度较大的场地。     

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.