Consolidation analysis of saturated multi‐layered soils with anisotropic permeability caused by a point sink

This paper presents the analytical layer‐element method to analyze the consolidation of saturated multi‐layered soils caused by a point sink by considering the anisotropy of permeability. Starting from the governing equations of the problem, the solutions of displacements and stresses for a single soil layer are obtained in the Laplace–Hankel transformed domain. Then, the analytical layer‐element method is utilized to further derive the solutions for the saturated multi‐layered soils in the transformed domain by combining with the boundary conditions of the soil system and continuity conditions between adjacent layers. The actual solutions in the physical domain can be acquired by the inversion of Laplace–Hankel transform. Numerical results are carried out to show the accuracy and stability of the proposed method and evaluate the influence of sink depth and anisotropic permeability on excess pore pressure and surface settlement. Copyright © 2011 John Wiley & Sons, Ltd.     

Static analysis of vertically loaded pile group in multilayered transversely isotropic soils

In this paper, a coupling approach is presented to study the static responses of vertically loaded pile group embedded in multilayered transversely isotropic soils. The individual pile in pile group is modeled by the finite element method, while the analytical layer-element method is applied to represent the soil's behavior. Then, the interaction equation of piles and soils is obtained by considering the force equilibrium and displacement compatibility conditions and solved by a FORTRAN program. The results computed by the proposed approach compare favorably with those from some existing solutions and field test. Some typical parametric analysis cases are investigated to study the effect of soil anisotropy, pile stiffness ratio, and pile spacing on the behavior of vertically loaded pile group.     

考虑桩土侧移的被动桩中土拱效应数值分析

被动桩对侧向位移的土层起到遮拦作用的机制主要是土拱效应。采用土工有限元软件Plaxis Tunnel 3D 1.2,对堆载荷载作用下邻近桩基中的土拱效应产生机制和性状进行三维数值分析,指出目前被动桩中土拱效应二维有限元分析存在的问题。考虑桩土侧移与相对位移,再利用土工有限元软件Plaxis2D 8.2详细地研究了侧向土体位移大小、桩身水平位移大小、土体性质以及桩土接触面性质等影响因素对土拱效应性态和桩土荷载分担比的影响。     

A quasistatic analysis of a plate on consolidating layered soils by analytical layer‐element/finite element method coupling

In this paper, a coupling method between finite element and analytical layer‐elements is utilized to analyze the time‐dependent behavior of a plate of any shape and finite rigidity resting on layered saturated soils. Based on the integral transform techniques together with the aid of an order reduction method, an analytical layer‐element solution is derived from the governing equations for three‐dimensional Biot consolidation with respect to a Cartesian coordinate system and then extended to be the fundamental solution for the layered saturated soil under a point load. The Mindlin plate is modeled by eight‐noded isoparametric elements. The governing equations of the interaction between soil and plate in the Laplace‐Fourier transformed domain are deduced by referring to the coupling theory of FEM/BEM, and the final solution is obtained by applying numerical inversion. Numerical examples concerned with the time‐dependent response of a plate are performed to demonstrate the influence of soil and plate properties on the interaction process. Copyright © 2014 John Wiley & Sons, Ltd.     

A simplified analysis method for piled raft foundations in non‐homogeneous soils

A simplified method of numerical analysis based on elasticity theory has been developed for the analysis of axially and laterally loaded piled raft foundations embedded in non‐homogeneous soils and incorporated into a computer program “PRAB”. In this method, a hybrid model is employed in which the flexible raft is modelled as thin plates and the piles as elastic beams and the soil is treated as springs. The interactions between structural members, pile–soil–pile, pile–soil–raft and raft–soil–raft interactions, are approximated based on Mindlin's solutions for both vertical and lateral forces with consideration of non‐homogeneous soils. The validity of the proposed method is verified through comparisons with some published solutions for single piles, pile groups and capped pile groups in non‐homogeneous soils. Thereafter, the solutions from this approach for the analysis of axially and laterally loaded 4‐pile pile groups and 4‐pile piled rafts embedded in finite homogeneous and non‐homogeneous soil layers are compared with those from three‐dimensional finite element analysis. Good agreement between the present approach and the more rigorous finite element approach is demonstrated. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical studies of anchored sheet pile wall behavior constructed in cut and fill conditions

The construction of sheet pile walls may involve either excavation of soils in front or backfilling of soils behind the wall. These construction procedures generate different loading conditions in the soil and therefore different wall behavior should also be expected. The conventional methods, which are based on limit equilibrium approach, commonly used in the design of anchored sheet pile walls do not consider the method of construction. However, continuum mechanics numerical methods, such as finite element method, make it possible to incorporate the construction method during the analyses and design of sheet pile walls. The effect of wall construction type for varying soil conditions and wall heights were investigated using finite element modeling and analysis. The influence of construction method on soil behavior, wall deformations, wall bending moments, and anchor forces were investigated. The study results indicate that walls constructed by backfill method yield significantly higher bending moments and wall deformations. This paper presents the results of the numerical parametric study performed and comparative analyses of the anchored sheet pile walls constructed by different construction methods.     

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.     

Axial and lateral response of pile groups embedded in nonhomogeneous soils

A numerical method of analysis based on elasticity theory is presented for the analysis of axially and laterally loaded pile groups embedded in nonhomogeneous soils. The problem is decomposed into two systems, namely the group piles acted upon by external applied loads and pile–soil interaction forces, and a layered soil continuum acted upon by a system of pile–soil interaction forces at the imaginary positions of the piles. The group piles are discretized into discrete elements while the nonhomogeneous soil behaviour is determined from an economically viable finite element procedure. The load–deformation relationship of the pile group system is then determined by considering the equilibrium of the pile–soil interaction forces, and the compatibility of the pile and soil displacements. The influence of soil nonlinearity can be studied by limiting the soil forces at the pile–soil interface, and redistributing the ‘excess forces’ by an ‘initial stress’ process popular in elasto-plastic finite element analysis. The solutions from this approach are compared with some available published solutions for single piles and pile groups in homogeneous and nonhomogeneous soils. A limited number of field tests on pile groups are studied, and show that, in general, the computed response compares favourably with the field measurements.     

A study of the effect of driving on pre‐bored piles

A finite element model for pile‐driving analysis is developed and used to investigate the behaviour of pre‐bored piles, which are then driven the last 1.25 or 2.25 m to their final design depth. The study was conducted for the case of saturated clays. The model traces the penetration of the pile into the soil and accommodates for large deformations. The non‐linear behaviour of the clay in this study is predicted using the bounding‐surface‐plasticity model, as applied to isotropic cohesive soils. The details of the 3‐D numerical modelling and computational schemes are presented. A significant difference was observed in the pile displacement during driving, and in the computed soil resistance at the pile tip, particularly at the earliest driving stages. No difference in soil resistance at the soil pile interface along the pile shaft was detected between the pre‐bored piles whether driven 1.25 or 2.25 m. Copyright © 2002 John Wiley & Sons, Ltd.     

灌注桩桩土相互作用试验及有限元模拟研究

以灌注桩桩土相互作用的原位试验为基础, 结合西安土层结构、性质; 采用有限单元法对黄土地区的砼灌注桩桩土相互作用进行仿真模拟, 研究了桩土相互作用、荷载的传递规律、桩土相对位移与桩侧摩阻力的关系。通过桩载试验资料对比, 得出三维有限单元法的模拟结果与实测值相近, 该研究进一步揭示了桩土相互作用的实质, 为最优桩长的选取和单桩承载力的确定奠定了基础。     

Effect of compressible parameters on vertical vibration of an elastic pile in multilayered poroelastic media

This paper mainly investigates the influences of compressible parameters on the vertical vibration of a pile embedded in layered poroelastic soil media. The pile is treated as a 1D elastic bar by the finite element method, and fundamental solutions for the layered poroelastic soils due to a vertical dynamic load are obtained by the analytical layer element method. Based on the compatibility conditions, the pile-soil dynamic interaction problem is solved. The numerical scheme has been compiled into a Fortran program for numerical calculation. Influences of the pile-soil stiffness ratio, compressible parameters, vibration frequency and the soil stratification are discussed.     

复合地基承载特性的弹塑性分析

对桩及承台采用线弹性有限元模型,对承台下桩周土采用弹塑性有限元模型,对群桩以外的土体采用线弹性无限元模型,在桩土接触面上设置接触面单元,利用三维弹塑性有限元对桩%D土%D承台相互作用进行了分析。得出了如下结论 :承台下桩顶反力总体表现出角桩最大,边桩次之,中桩最小的分布规律,随着作用在承台上的荷载增大,桩顶反力趋于均匀分布,承台下桩侧摩阻力是由桩端向桩顶逐渐发展的,承台对桩上部侧摩擦阻力存在削弱作用。为了验证本文方法的可行性,对承台下有九桩的情况进行了静载试验,将试验结果与本文计算结果进行了比较。     

Analysis of piled raft with nodular pile subjected to vertical load using a Winkler model approach

An investigation is made to present analytical solutions provided by a Winkler model approach for analysis of piled rafts with nodular pile subjected to vertical loads in nonhomogeneous soils. The vertical stiffness coefficient along a piled raft with the nodular pile in nonhomogeneous soils is derived from the displacement given by the Mindlin solution for elastic continuum analysis. The vertical stiffness coefficients for the bases of the raft and the nodular part in the nodular pile in a soil are expressed by the Muki solution for the 3‐D elastic analysis. The relationship between settlement and vertical load on the pile base is presented considering the Mindlin solution and the equivalent thickness in the equivalent elastic method. The interaction factor between the shaft of the nodular pile and the soil is expressed taking into account the Mindlin solution and the equivalent elastic modulus. The relationship between settlement and vertical load for a piled raft with the nodular pile in nonhomogeneous soils is obtained by using the recurrence equation of influence factors of the pile for each layer. The percentage of each load carried by both nodular pile and raft subjected to vertical load is represented through the vertical influence factors proposed here. Comparison of the results calculated by the present method for piled rafts with nodular piles in nonhomogeneous soils has shown good agreement with those obtained from the finite element method and a field test. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of pile‐raft foundations under complex loads in layered soils

Considering there is hardly any concerted effort to analyze the pile‐raft foundations under complex loads (combined with vertical loads, horizontal loads and moments), an analysis method is proposed in this paper to estimate the responses of pile‐raft foundations which are subjected to vertical loads, horizontal loads and moments in layered soils based on solutions for stresses and displacements in layered elastic half space. Pile to pile, pile to soil surface, soil surface to pile and soil surface to soil surface interactions are key ingredients for calculating the responses of pile‐raft foundations accurately. Those interactions are fully taken into account to estimate the responses of pile‐raft foundations subject to vertical loads, horizontal loads and moments in layered soils. The constraints of the raft on vertical movements, horizontal movements and rotations of the piles as well as the constraints of the raft on vertical movements and horizontal movements of the soils are considered to reflect the coupled effect on the raft. The method is verified through comparisons with the published methods and FEM. Then, the method is adopted to investigate the influence of soil stratigraphy on pile responses. The study shows that it is necessary to consider the soil non‐homogeneity when estimating the responses of pile‐raft foundations in layered soils, especially when estimating the horizontal responses of pile‐raft foundations. The horizontal loads and the moments have a significant impact on vertical responses of piles in pile‐raft foundations, while vertical loads have little influence on horizontal responses of piles in pile‐raft foundations in the cases of small deformations. The proposed method can provide a simple and useful tool for engineering design. Copyright © 2013 John Wiley & Sons, Ltd.     

改进的薄层单元在桩-土动力相互作用中的应用

在地震荷载作用过程中,由于桩身材料性能与周围土体性质差异较大,土体在桩-土接触面上发生张开或滑移,该强非线性接触行为直接影响接触面附近土体与桩体应力状态,从而影响上部结构地震响应水平。真实桩-土接触面由于桩身表面混凝土约束关系,具有一定厚度且存在相关体变规律的接触带。通过有限元软件ABAQUS用户自定义单元UEL程序编写改进Desai薄层接触单元,在Desai薄层接触单元中加入Rayleigh阻尼项以模拟地震作用时桩-土强非线性接触行为能量耗散过程。接触面法向与两切向本构关系采用双曲线模型。规定了Desai薄层接触单元在桩-土接触面上的行为模式以模拟土体在接触面黏结、滑动、张开、再闭合等接触状态。建立精细三维桩-土-结构体系动力相互作用模型以研究改进Desai单元对上部结构峰值动力响应水平影响,为工程结构抗震设计时程分析提供参考依据。     

抗滑桩极限阻力的有限元分析

应用平面应变有限元法分析了不同强度、应力状态下土体的单桩极限抗滑阻力,并与现有的极限阻力经验公式的计算结果进行了比较,论证了结果的合理性.对不同桩径及不同土体力学参数情况下的桩极限阻力进行了分析,研究了桩径及土体参数等因素对桩极限阻力的影响.为进一步利用该方法研究单桩极限抗滑阻力的计算公式提供了基础.     

Analysis of rectangular piles subjected to lateral loads in nonhomogeneous soil using a Winkler model approach

An analytical approach using a Winkler model based on two lateral soil displacement components in a three‐dimensional soil is investigated to provide analytical solutions of horizontal response of a rectangular pile subjected to lateral loads in nonhomogeneous soil. The two lateral displacement components of a soil surrounding the rectangular pile are represented by the Fourier series of displacement potential functions in the elastic three‐dimensional analysis. The lateral stiffness coefficient of the rectangular pile shaft in nonhomogeneous soil is derived from the rocking stiffness coefficient taking into account rocking rotation of a rigid pile shaft. The relationship between horizontal displacement, rotation, moment, and shear force for the rectangular pile subjected to horizontal loads in nonhomogeneous soil is obtainable in the form of the recurrence equation. The formulation of lateral displacement and rotation for a rectangular pile subjected to lateral loads on the pile base in nonhomogeneous soil is proposed by taking into account Mindlin's equation and the equivalent thickness for soil layers in the equivalent elastic method. The difference of lateral behavior between square and circular piles subjected to lateral loads is insignificant. The effect of aspect ratio of the rectangular pile on the lateral behavior is great for the lower stiffness ratio between pile and soil and the larger length–equivalent diameter ratio. The effect of the value of Poisson's ratio of soil on lateral stiffness coefficient is relatively small except Poisson's ratio close to 0.5. The comparison of the results calculated by the current method for a rectangular pile subjected to lateral loads in nonhomogeneous soil has shown good agreement with those obtained from the analytical methods and the finite element method. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical modeling of dredging effect on berthing structure

Piles and diaphragm wall-supported berthing structure on marine soils are loaded laterally from horizontal soil movements generated by dredging. The literature on the adequacy of the finite element method modeling of berthing structure to analyze their behavior during dredging is limited. This paper describes a finite element approach for analyzing the lateral response of pile and diaphragm wall during dredging. Piles are represented by equivalent sheet-pile walls and a plane strain analysis using the finite element method is performed. Results from the finite element method are compared with full-scale field test data. Full-scale field test was conducted on a bearing structure to measure the lateral deflection on pile and diaphragm wall for their full length using inclinometer during dredging in sequence. The finite element method results are in good agreement with full-scale field results. Conclusions are drawn regarding application of the analytical method to study the effect of dredging on piles and diaphragm wall-supported berthing structures.     

砂土中单桩竖向抗压承载机制的离散元分析

将单桩视为平面问题,采用二维离散单元法分析了砂土中单桩竖向抗压承载机制。首先利用离散元模拟地基的形成、挖孔灌注桩的成桩过程以及地基土的双轴压缩试验;然后通过离散元模拟单桩竖向抗压静载试验,分析其承载机制。结果表明：随着桩的沉降增大,桩端阻力一直增大,桩侧摩阻力先增大后达到稳定值;单位桩长侧摩阻力从桩顶到桩端呈非线性增长趋势;当荷载达到桩的极限承载力时,桩端发生刺入破坏,导致桩侧与桩端附近土体的转动场梯度、应力场梯度增大,主应力发生大角度旋转。     

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.