Development of a three-dimensional soil model for the dynamic analysis of end-bearing pile groups subjected to vertical loads

This paper presents a new analytical model for calculating the steady-state impedance of pile groups subjected to vertical dynamic loads. The derived solution allows considering effects of radially but also vertically propagating soil waves on the soil attenuation function, pile interaction factor, and pile group impedance. The proposed model provides accurate estimates of the soil stress field and of the response of the pile group in the low as well as in the high-frequency range, unlike earlier solutions based on the plane-strain model to describe the soil surrounding the piles, which ignores the vertical soil stress gradient. The latter assumption results in underestimating pile group impedance and overestimating radiation damping for frequencies lower than the cutoff frequencies of the system, which are explicitly captured with the proposed solution.     

竖向入射S波作用下单桩水平地震响应简化解析方法

基于连续介质模型并考虑桩-土运动相互作用,将单桩视为一维线弹性梁,研究了竖向入射S波作用下的单桩水平地震响应问题。将竖向入射S波模拟为基岩水平位移,基于平面应变模型建立的土体控制方程,推导出地震作用下土体水平动力阻抗函数表达式。将土体阻抗代入单桩控制方程并联立桩-土接触条件及桩顶和桩底的边界条件,得出了竖向入射S波作用下单桩的地震响应解析解。通过将所得解与已有文献理论解和有限元结果进行对比,验证了该方法的合理性。基于所得理论解进行参数分析得出：桩-土模量比的增加可以明显降低桩-土运动相互作用因子的最小值,而较大桩身长径比以及土体滞回阻尼对桩-土运动相互作用因子的影响较小;对于桩顶水平地震放大因子来说,桩-土模量比的增加仅在高共振频率处抑制其幅值,较大桩身长径比对其影响较小,而滞回阻尼比的增加会显著抑制共振频率处的幅值;桩身地震响应仅在较小桩径比时受桩-土模量比的影响明显,并随桩-土模量比的增加而降低。     

横观各向同性地基中管桩扭转振动响应解析解

考虑地基沉积过程中产生的竖向和水平向力学性质的差异,对横观各向同性地基中管桩扭转振动频域响应进行了理论研究。基于横观各向同性材料的本构关系以及桩-土耦合扭转振动,建立了桩土系统定解问题,通过Laplace变换和分离变量法求得了桩周土和桩芯土扭转振动位移形式解。通过桩-土接触面的连续条件,求得了管桩扭转频域响应解析解,并得到了桩顶复动刚度和速度导纳的表达式。将所得解退化到横观各向同性地基中实心桩解以及均匀地基中管桩解,并与已有文献进行了对比,验证了解的合理性。通过数值算例,分析了桩周土和桩芯土的横观各向同性力学参数对桩顶扭转复刚度及速度导纳的影响。     

桩在饱和土中水平振动的辐射阻尼简化算法

将适用于单相土－桩系统的Novak的平面应变假设和Gazetas的分象限假设推广到饱和土－桩系统中。通过简化后的微分方程推导了饱和土－桩系统辐射阻尼系数的解析解答。最后用算例说明了饱和土中两个压缩波的特性，比较了饱和土和单相土的单位桩长辐射阻尼系数，比较了简化算法和三维有限元方法求得的桩头幅射阻尼比，结果证明简化算法在一定条件下适用于计算层饱和土－桩系统的辐射阻尼。     

砂墙地基二维固结自由应变解

将含竖向排水体地基的三维固结变形问题等效为平面应变问题进行数值分析时,砂墙地基二维固结解析解答是合理确定其等效固结计算参数的重要依据。为辨析现有砂墙地基等应变固结近似解答的适用性,针对微单元土体严格的二维固结微分方程,考虑对地基固结有重要影响的井阻作用,以及涵盖完全透水和不完全透水的更一般边界面排水条件,推求得到了其自由应变解答。并与现有解答进行对比分析,同时,分析了泊松比效应以及水平和竖向排水对地基固结的影响。结果表明,现有砂墙地基的等应变固结解答虽然近似,但已有足够精确;砂墙地基以水平向固结为主,竖向固结几乎可以忽略不计;地基固结速率随着泊松比的增大而增大,在将竖向排水体等效为砂墙时,应考虑其作用影响。     

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.     

Boundary element analysis of axially loaded piles embedded in a multi-layered soil

In a field, piles are likely installed in a multi-layered soil. Analysis of axially loaded piles in a multi-layered soil is complicated and deserves more attention. A boundary element method is used in this study to analyze an axially loaded single pile in a multi-layered soil using the solution for vertical and horizontal axisymmetric ring loads in a multi-layered elastic medium. Good and reasonable agreement is obtained between the proposed and published solutions for a single pile in a homogenous soil, a finite soil, and a Gibson soil. The proposed solution is also used to evaluate an axially loaded single pile in a multi-layered (8 layers) soil.     

基于径向不变假定的现浇大直径管桩纵向振动响应频域解

假定各物理量沿径向不变化,建立了低应变瞬态集中荷载作用下现浇大直径管桩振动响应的计算模型和波动方程。采用Laplace变换法,求得了波动方程的频域解析解,采用Fourier逆变换求得了时域响应。将文中解计算结果与三维频域解析解进行了对比分析,文中解的入射峰-反射峰时间差与三维频域解差别很小,说明对于PCC桩这种大直径薄壁管桩,采用径向不变假定对计算结果几乎没有影响。因此,在PCC桩低应变检测波形分析时,采用基于径向不变假定的二维解是完全合理的。将文中二维频域解的计算结果与二维时域解的结果进行了对比分析,结果表明,2种解在入射波、第一个桩底反射波、第二个桩底反射波峰值大小和到达时间非常吻合,这说明考虑环向位移与否对计算结果没有太大的影响,研究PCC桩低应变动力响应的问题时忽略水平方向的位移是完全可以的。     

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.     

地下综合管廊边界条件对地震动力响应影响数值分析

为研究综合管廊动力边界条件对地震动力响应的影响,以厦门地区的代表性土层为例,建立动力有限元数值模型,土体本构采用小应变硬化模型,分别设定固定边界、黏性边界和自由场3种人工边界条件,进行Rayleigh波和地震底部剪切波作用下的场地响应研究;并根据变形特征及拟绝对加速度反应谱（PSA）评价3种边界的有效性,提出综合管廊地震动力分析的优化动力边界组合方法。研究表明：在地震波（底部水平加速度时程）及Rayleigh波的作用下,由于考虑了黏性边界对外行波的吸收,但未考虑地震动的输入问题及边界外半无限介质的弹性恢复性能,边界会对模型内部土体的水平位移产生限制作用,使得场地内水平位移响应偏小,而采用自由场边界则基本不存在这种限制作用,表现出强烈的振荡;采用激励侧固定边界、远离激励侧黏性边界、其余侧自由场边界的优化组合动力边界,在Rayleigh波和底部加速度时程共同作用下,二者引起的动力响应交叉干扰较少,可按线性叠加处理;同时,黏性边界对地震波引起的动力响应有一定范围的吸收,自由场边界对Rayleigh波引起的动力响应也有一定范围的变形限制影响。研究成果可供地下综合管廊结构地震响应精细化数值模拟及抗震设计参考。     

被动桩水平位移荷载施加位置的探讨

被动桩是指一种由于土体水平位移而发生挠曲变形的桩。目前不论是室内试验还是数值计算,常在模型边界对土体施加水平位移荷载（位移边界条件）形成土体的位移场,用来研究被动桩的变形特征,但不同研究者采用的水平位移施加位置（边界位置）与被动桩的距离存在较大差异。结合一个工程案例,采用有限单元法对水平位移加载位置进行了单变量参数分析。研究表明,随着水平位移加载位置与被动桩距离的增加,桩身变形显著减小。这表明了施加位移边界后地基中存在明显的应力扩散。当该距离小于5.5倍桩径时,在模型边界施加同等大小的水平位移所需应力显著增加。同时讨论了水平荷载加载位置与被动桩距离的合理范围,认为应在5～8倍桩径范围之内。     

Study on low‐strain integrity testing of pipe‐pile using the elastodynamic finite integration technique

This study focuses on the three‐dimensional (3‐D) characteristics of wave propagation in pipe‐pile using elastodynamic finite integration technique. First, a real 3‐D pile‐soil model in cylindrical coordinate system is presented. Then, the governing equations are established. With the boundary and initial conditions, the numerical solution is obtained. The accuracy and feasibility of the self‐written code are further verified via comparing with the measured data. Velocity histories at different angles of pile top and pile tip are illustrated, and the snapshots reflecting the 3‐D characteristics of wave propagation are given. It shows that the interferences of Rayleigh waves can confuse the result interpretation for pile integrity testing. The increase of hammer contact time can effectively mitigate the interferences, and the interferences of Rayleigh waves are weakest at an angle of 90° from where hammer hits. Besides, surrounding soil can partly mitigate the wave interferences. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

边坡土体侧移模式对抗滑桩影响分析

在开挖、降雨或地震等外部因素作用下,边坡土体很容易进入局部或瞬态大变形乃至失稳滑动,使抗滑桩产生附加位移及弯矩。基于两阶段分析方法,采用Winkler模型模拟抗滑单桩与土之间的相互作用,建立单桩水平位移控制方程组,根据内力与位移的连续条件得到考虑不同土体侧移模式下求解桩身响应的矩阵解析表达式,并采用现场监测数据及Poulos弹性理论进行验证,证明该方法是合理可行的,并通过参数分析土体侧移对抗滑桩水平承载性状的影响程度。分析结果表明,土体侧移模式包括最大侧移值、分布形状及重心、侧移势等,对抗滑桩的挠度和弯矩均有显著影响,在工程设计中应予以充分重视。     

非均质土中桩端扩散虚土桩法的桩基纵向振动研究

根据桩端土应力扩散的规律,建立了桩端扩散虚土桩模型。基于该模型对非均质土中桩-土纵向耦合振动进行研究。利用复刚度传递多圈层平面应变模型,得到桩与虚土桩桩侧土的剪切复刚度。结合边界条件、初始条件和连续条件,对扩散虚土桩和实体桩动力方程从底层往顶层逐层进行求解,得到桩顶动力响应的频域解析解和时域半解析解。通过对桩端扩散虚土桩扩散角、扩散层厚度、桩侧土非均质性和桩长的影响进行计算分析,得到基于扩散虚土桩法桩-土纵向振动响应特性。研究结论可为桩基础动力设计和动态检测提供理论依据。     

Poroelastic model for pile–soil interaction in a half‐space porous medium due to seismic waves

In this paper, frequency domain dynamic response of a pile embedded in a half‐space porous medium and subjected to P, SV seismic waves is investigated. According to the fictitious pile methodology, the problem is decomposed into an extended poroelastic half‐space and a fictitious pile. The extended porous half‐space is described by Biot's theory, while the fictitious pile is treated as a bar and a beam and described by the conventional 1‐D structure vibration theory. Using the Hankel transformation method, the fundamental solutions for a half‐space porous medium subjected to a vertical or a horizontal circular patch load are established. Based on the obtained fundamental solutions and free wave fields, the second kind of Fredholm integral equations describing the vertical and the horizontal interaction between the pile and the poroelastic half‐space are established. Solution of the integral equations yields the dynamic response of the pile to plane P, SV waves. Numerical results show the parameters of the porous medium, the pile and incident waves have direct influences on the dynamic response of the pile–half‐space system. Significant differences between conventional single‐phase elastic model and the poroelastic model for the surrounding medium of the pile are found. Copyright © 2007 John Wiley & Sons, Ltd.     

Load transfer approach for laterally loaded piles

A two‐parameter model has been proposed previously for predicting the response of laterally loaded single piles in homogenous soil. A disadvantage of the model is that at high Poisson's ratio, unreliable results may be obtained. In this paper, a new load transfer approach is developed to simulate the response of laterally loaded single piles embedded in a homogeneous medium, by introducing a rational stress field. The approach can overcome the inherent disadvantage of the two‐parameter model, although developed in a similar way. Generalized solutions for a single pile and the surrounding soil under various pile‐head and base conditions were established and presented in compact forms. With the solutions, a load transfer factor, correlating the displacements of the pile and the soil, was estimated and expressed as a simple equation. Expressions were developed for the modulus of subgrade reaction for a Winkler model as a unique function of the load transfer factor. Simple expressions were developed for estimating critical pile length, maximum bending moment, and the depth at which the maximum moment occurs. All the newly established solutions and/or expressions, using the load transfer factor, offer satisfactory predictions in comparison with the available, more rigorous numerical approaches. The current solutions are applicable to various boundary conditions, and any pile–soil relative stiffness. Copyright © 2001 John Wiley & Sons, Ltd.     

Soil‐pile interaction in the pile vertical vibration considering true three‐dimensional wave effect of soil

The dynamic response of an end bearing pile embedded in a linear visco‐elastic soil layer with hysteretic type damping is theoretically investigated when the pile is subjected to a time‐harmonic vertical loading at the pile top. The soil is modeled as a three‐dimensional axisymmetric continuum in which both its radial and vertical displacements are taken into account. The pile is assumed to be vertical, elastic and of uniform circular cross section. By using two potential functions to decompose the displacements of the soil layer and utilizing the separation of variables technique, the dynamic equilibrium equation is uncoupled and solved. At the interface of soil‐pile system, the boundary conditions of displacement continuity and force equilibrium are invoked to derive a closed‐form solution of the vertical dynamic response of the pile in frequency domain. The corresponding inverted solutions in time domain for the velocity response of a pile subjected to a semi‐sine excitation force applied at the pile top are obtained by means of inverse Fourier transform and the convolution theorem. A comparison with two other simplified solutions has been performed to verify the more rigorous solutions presented in this paper. Using the developed solutions, a parametric study has also been conducted to investigate the influence of the major parameters of the soil‐pile system on the vertical vibration characteristics of the pile. Copyright © 2013 John Wiley & Sons, Ltd.     

基于传递矩阵法的层状土中管桩水平动力阻抗分析

在考虑土体分层特性的基础上,分别建立了管桩桩周土体和桩芯土体的水平振动控制方程。通过引入势函数并考虑桩周土和桩芯土径向位移和环向位移的边界条件及其奇偶性,求得了管桩-土动力相互作用的刚度系数和阻尼系数。将土体模拟为连续分布的弹簧-阻尼器,并考虑桩芯土和桩周土的作用,建立了层状土中管桩的水平振动方程。借助初参数法和传递矩阵法求解了管桩的水平振动,得到了管桩桩顶的水平动力阻抗。通过数值分析,得到了土层剪切模量、管桩壁厚、桩周土和桩芯土剪切模量比、土层厚度等对管桩桩顶动力阻抗的影响规律。土层剪切模量、管桩壁厚、桩周土和桩芯土剪切模量比对层状土中管桩水平振动的影响主要在低频处,土层厚度在较宽的频率范围内对管桩水平振动有影响;管桩壁越厚,桩周土的剪切模量越大时,管桩水平动力阻抗的绝对值越大。     

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.