成层土中竖向受荷桩沉降性状简化计算方法

王立兴 吴文兵 梅国雄 闻敏杰 徐美娟

摘要：基于应力泡形虚土桩模型和考虑桩土界面初始临界摩阻力的改进双曲线模型，提出了一种可以考虑桩端应力扩散效应和桩土界面初始临界摩阻力影响的竖向受荷单桩沉降简化计算方法.首先，利用该计算方法分析讨论了模型参数对单桩沉降的影响程度，并给出了合理的取值范围；其次，将用该计算方法得到的结果与现场实测数据对比分析，验证了本文计算方法的合理性，并进一步分析研究了不同桩顶荷载作用下竖向受荷桩的桩身压缩与桩端土体沉降所占比重关系、桩身轴力分布规律以及桩侧摩阻力发挥特性等.本文的计算方法参数明确，应用简便，可以为桩基础设计提供更为准确的沉降值，具有广泛的工程应用前景.

关键词：竖向受荷桩；沉降性状；虚土桩模型；双曲线模型；桩身轴力；桩侧摩阻力

中图分类号：TU473.1文献标志码：A

Simplified Calculation Method for Settlement Behavior of Vertically Loaded Pile in Layered Soil

WANG Lixing1，WU Wenbing1，MEI Guoxiong1，2，WEN Minjie1，3，XU Meijuan2

（1. Faculty of Engineering，China University of Geosciences，Wuhan 430074，China；2. College of Civil Engineering and Architecture，Guangxi University，Nanning 530004，China；3. School of Civil Engineering and Architecture，Zhejiang Sci-Tech University，Hangzhou 310018，China）

Abstract：Based on the stress-bubble fictitious soil-pile model and the improved hyperbolic model considering the initial critical frictional resistance of the pile-soil interface，a simplified method for calculating the settlement of vertically loaded single pile is proposed，which can consider the effect of stress diffusion at the pile tip and the initial critical frictional resistance at the pile-soil interface. Firstly，the influence degree of the model parameters on the settlement of a single pile is analyzed and discussed by using the present calculation method，and the reasonable value ranges of model parameters are given. Then，the rationality of the present calculation method is verified by comparing the results obtained using the present calculation method with the field measured data. Furthermore，the proportion of pile compression and pile tip soil settlement，the distribution law of pile axial force and the characteristics of pile lateral friction resistance under different pile top loads are further analyzed and studied. The present calculation method has clear parameters and simple application，which can provide a more accurate settlement value for pile foundation design，thus having a wide engineering application prospect.

Key words：vertically loaded pile；settlement behavior；fictitious soil-pile model；hyperbolic model；pile axial force；pile lateral friction resistance

樁基础由于承载力高、稳定性好、适应性强等特点，已经被广泛应用于高层建筑、桥梁公路等构筑物的基础建设中.在桩基设计中，竖向受荷桩的承载机理和沉降特性分析是至关重要的.目前，基桩沉降计算方法有荷载传递法[1-2]、弹性理论法[3-4]、剪切位移法[5-6]和数值计算法[7-8]等，其中荷载传递法是Seed和Reese[1]基于桩身荷载传递机理提出的，可以灵活地选取不同的荷载传递函数模拟桩与桩侧土非线性相互作用特性及地基土的层状特性，得到了广泛应用.荷载传递法的准确性主要依赖于荷载传递函数的合理性，荷载传递函数即剪应力与桩土相对位移之间的关系，传统的荷载传递函数主要有理想弹塑性模型[9-10]、双折线硬化模型[2，11]、三折线模型[12-13]、双曲线模型[14-15]及侧阻软化模型[11，16-17].然而土与混凝土结构面的剪切试验[18-19]发现当剪应力小于某一个临界值时无相对位移产生，表明土与混凝土结构物之间应当存在初始临界摩阻力，赵明华等[20]和Zhang等[21]分别提出了考虑初始临界摩阻力的改进双曲线模型和侧阻软化模型用于桩基沉降计算，通过实例分析验证了考虑初始临界摩阻力的非线性荷载传递模型可以更准确地模拟桩土相互作用特性，具有更好的工程适用性.

另外，传统的荷载传递法进行基桩沉降计算时简化地用折线硬化模型[11，21]或双曲线模型[17，20]来表征桩端荷载与位移的关系，这样无法考虑桩端土体的层状特性和桩端应力扩散效应.为此，王奎华等[9-10]、Wu等[22]提出虚土桩模型来模拟桩端土对桩的支撑作用，并将虚土桩模型应用于层状地基中的单桩沉降计算分析和基桩动力特性研究[23-25].进一步，王立兴等[26]将应力泡形虚土桩模型应用到单桩沉降分析来考虑桩端应力扩散效应的影响，但是文中采用双折线函数模拟桩侧阻力具有一定局限性，因此，有必要进一步对采用侧阻非线性荷载传递模型结合端阻应力泡形虚土桩模型的成层土中竖向受荷桩的沉降特性展开研究.

综上，本文采用应力泡形虚土桩模型和考虑桩土界面初始临界摩阻力的改进双曲线模型，分别模拟桩端土体对桩体的支撑作用和桩侧摩阻力与桩土相对位移间的关系，并在上述两种模型的基础上提出一种考虑桩土界面初始临界摩阻力影响和桩端应力扩散效应的成层土中竖向受荷桩的沉降特性简化迭代计算方法.通过工程实例对比验证了本文提出计算方法的合理性，并进一步分析研究了不同桩顶荷载作用下，竖向受荷桩的桩顶荷载沉降曲线、桩身压缩与桩端沉降所占比重变化情况、桩身轴力分布规律以及桩侧摩阻力发挥特性等单桩沉降性状.

1计算模型及基本假设

1.1桩侧荷载传递模型

为考虑桩土界面初始临界摩阻力和桩土非线性相互作用，本文采用赵明华等[20]提出的改进双曲线荷载传递模型来模拟桩侧摩阻力与桩土相对位移之间的关系，如图1所示，该模型适用于正常固结黏土、中密砂与松砂等土体类型，其函数表达式为：

式中：τ（z）是深度z处的桩侧摩阻力；S（z）是深度z处的桩土相对位移；参数a、b为模型参数，其物理意义分别为曲线原点处的切线斜率的倒数和桩侧极限摩阻力τ的倒数；c（z）为深度z处桩土界面初始临界静摩阻力.参数a、b和c（z）的值可通过现场载荷试验结果的反分析确定或通过室内土-结构接触面剪切试验获得.

1.2桩端泡形虚土桩模型

辛公锋[27]通过对4根现场试桩实测数据的反演，发现桩在极限承载力下桩底土体大部分仍处于弹性状态，因此，可以采用Boussinesq解计算樁端附加应力范围，并利用弹性力学知识进行桩端土体沉降计算.假设桩端对土体的作用力为均布荷载形式，在桩端底面建立如图2所示的极坐标系，将Boussinesq解进行二重积分可以求得圆形均布荷载作用下桩端土体中任意一点的竖向应力分布为[25-26]：

对式（2）进行数值求解可以得到一组竖向附加应力等值线，将等值线绕z轴旋转一周就形成三维附加应力等值面，假设应力等值面外的土体内附加应力很小，将等值面内的土体看作受荷主体进行计算，如图3所示.

1.3桩土相互作用模型及相关假设

桩土相互作用体系简图如图3所示，根据地基土层的成层性情况将实体桩和虚土桩一共分为n段，其中实体桩分成m段，虚土桩分成n-m段.从实体桩顶部向下依次编号为1，2，…，i，…，m，m+1，…，n；各桩段厚度分别为l，l，…，l，…，l.实体桩总长为L，虚土桩总长为L，各段桩的埋深分别为z，z，…，z，…，z.桩顶作用的竖向荷载为P，产生的桩顶位移为S，荷载向下传递，假设第i段桩的桩顶轴力为P，相应的桩顶位移为S，桩侧摩阻力为τ.

为简化计算，本文桩土体系相关假设有：

1）实体桩桩身截面均匀，虚土桩以应力泡的形式存在，其边界由式（2）求得；

2）虚土桩的力学行为与实体桩相符合，参数取土体实际参数；

3）虚土桩的变形仅考虑弹性状态，虚土桩底部为刚性支撑不发生变形；

4）实体桩与虚土桩的各桩段界面之间为完全连续接触，满足连续性条件；

5）地基中同一土层的物理力学性质沿深度方向不变.

2单桩沉降计算流程

单桩的沉降由桩端土体的沉降和实心桩体的压缩两部分组成.本文采用应力泡形虚土桩模型来模拟桩底土体对实体桩的支撑作用，采用考虑初始临界摩阻力的双曲线荷载传递模型来模拟桩与桩周土体的相互作用，实体桩和虚土桩的压缩沉降量之和，即为所求实体桩的桩顶总沉降量.

首先推导虚土桩的沉降计算公式，根据静力平衡条件，第i（m<i<n）段桩位移应满足式（3）：

式中：E是第i段虚土桩的弹性模量；A是第i段桩的截面积，A=πR：各个虚土桩段的半径R由式（2）确定；f是第i段桩周单位摩阻力，满足式（4）：

f=λS（4）

式中：λ为第i段桩弹性剪切刚度系数，且λ=τ/S.τ和S分别为第i段桩的极限侧摩阻力和相应的极限位移，由原位静载试验获得.

由各桩段分界面的轴力和位移的连续性条件，可以得到式（5）：

式中：P为第i+1段桩桩顶轴力；S为第i+1段桩桩顶位移；z为第i段桩底面到桩顶的距离，即各段桩的埋深.

第i段桩桩顶轴力和位移可由式（6）求解得到：

进一步，由式（7）可得第i段桩桩顶刚度为：

由虚土桩端不发生变形，可以得到：

式（7）（8）组合可以得到第n段桩桩顶刚度：

根据上述推导，完整的基桩沉降预测流程可分为以下10个步骤：

1）将桩-土体系划分为n段，如图3所示；

2）假设第n段虚土桩的压缩位移为S，则第n段虚土桩的轴力P可通过式（9）计算；

3）利用式（6）可递归得到每段虚土桩截面顶部的轴力和位移，实体桩端的轴力和位移等于虚土桩桩顶的轴力和位移；

4）假设第m段实体桩的压缩位移为S，并将第m段实体桩桩顶位移值（S+S/2）代入式（1），就可得到第m段实体桩的桩周摩阻力τ；

5）利用式（10）计算第m段实心桩的轴力增量P：

P=lUτ（10）

式中：l是第m段实体桩的长度，U是第m段实体桩的截面周长，U=2πr，r是实体桩的半径；

6）利用式（11）计算第m段桩的桩身压缩S′：

7）将第m段实体桩更新后的压缩位移S′和步骤4）中的压缩位移S值比较，如果S′与S之差不在允许公差（如1×10）内，用S′替换S，然后重复步骤4）到步骤6），直到两者之差位于允许公差范围内；

8）将计算得到的位移增量S和轴力增量P代入式（12），得到第m段实体桩顶部的轴力P和位移S：

9）對于第i段实体桩，重复步骤4）到步骤8），得

各实体桩桩顶轴力P和位移S，直至得到第1段桩桩顶轴力（即桩顶竖向荷载）P和位移S；

10）假定一组递增的第n段虚土桩压缩位移为S，然后重复步骤2）到步骤9），便可得到一系列荷载位移值，最终得到桩顶荷载沉降曲线.

3计算模型参数讨论

桩端土体的沉降对桩顶沉降有较大的贡献，为了突出分析端应力扩散效应对桩顶沉降的影响，下面分析中取单层均质地基，着重对泡形虚土桩模型参数进行讨论，对于成层地基情况，只需要根据实际土层情况进行赋值即可.假定桩土体系参数取值为：实体桩长L=20 m，桩身半径r=0.5 m，桩身弹性模量E=3×10MPa.地基土层为均质硬黏土层，土层的极限侧摩阻力τ=120 kPa，相应的极限位移s=8 mm.土层的双曲线模型计算参数a=0.03 mm/kPa，b=0.007/kPa，c=10 kPa.模型参数是参考本文工程案例的实际参数选取的合理近似值.

3.1竖向附加应力σ

竖向附加应力σ是采用应力泡形虚土桩模型分析单桩沉降的一个重要参数，选取不同桩端土体的竖向附加应力σ代表选取了不同长度（深度）的泡形虚土桩模型，进一步将竖向附加应力σ和任意深度值z代入式（2）进行数值求解可以得到不同深度处的径向宽度.不同边界的应力泡形虚土桩模型如图4 所示，基于本文计算方法竖向附加应力σ对桩顶沉降的影响如图5所示.

从图5可以看出当桩顶荷载较小（P<5 000 kN）时，桩端竖向附加应力σ对桩顶沉降影响很小，这是由于当桩顶荷载较小时传递到桩端的基底压力和桩端土体中的竖向附加应力相应较小，使得不同竖向附加应力对应的泡形虚土桩产生的压缩变形差值就会很小；当桩顶荷载继续增大时，其对桩顶沉降的影响程度会明显增加，相同荷载下随着桩端竖向附加应力的减小桩顶沉降值随之减小，这是由于竖向附加应力的减小意味着泡形虚土桩的边界增大，即承受基底压力的受荷土体体积增大，那么整体的桩端土体沉降值会减小，桩顶沉降也随之减小；另外随着竖向附加应力的增大，桩顶沉降增大的幅度不断减小并且最后收敛于σ=0.01 P附近，因此，在利用应力泡形虚土桩模型进行单桩沉降计算时建议桩端竖向附加应力取值范围在0.005P～0.05 P较为合理.

3.2虚土桩分段精度

进一步讨论虚土桩的分段精度对桩顶沉降计算的影响.本文中虚土桩分段精度用虚土桩分段厚度l来进行量化分析，基于本文计算方法虚土桩分段厚度l对桩顶沉降的影响如图6所示.

从图6可以看出，当桩顶荷载较小（P<6 000 kN）时，虚土桩分段精度l对桩顶沉降影响很小；随着桩顶荷载的增大虚土桩分段精度l对桩顶沉降影响增加，并且随着分段精度l的增大，其对桩顶沉降影响幅度不断减小并且最后收敛于l=1m附近，因此，在利用应力泡形虚土桩模型进行单桩沉降计算时建议虚土桩分段精度%，取值在0.5～1.5 m较为合理.

3.3虚土桩弹性模量

利用土体弹性模量E的变化分析桩端土体性质对桩顶沉降的影响，基于本文计算方法桩端土体弹性模量E对桩顶沉降的影响如图7所示.

从图7可以看出，当桩顶荷载较小（P<6 000 kN）时，桩端土弹性模量E对桩顶沉降影响很小；随着桩顶荷载的增大桩端土弹性模量E对桩顶沉降影响程度增加，具体表现为桩端土弹性模量E越大，桩顶沉降量越小，这是因为桩端土弹性模量E越大表示桩端土体性质越好，抗压缩变形能力越强，因此整体桩身沉降量会减小.

4工程实例分析

为了验证本文所述基桩沉降计算方法的准确性，利用本文计算方法对文献［28］中的1号试桩进行沉降计算，将计算结果与文献［28］中的实测曲线进行比较以验证本文计算方法的合理性，并进一步利用本文计算方法对竖向受荷桩进行沉降性状分析.

试验场地地势平坦，地基土层主要由冲洪积相的亚黏土、黏土、粗砂组成.1号试桩的实测桩径为1.09 m，桩身弹性模量E为2.43×10 MPa，桩长为25 m，桩端持力层为粗砂层，根据文献［28］的实测数据，粗砂层的极限侧摩阻力为183 kPa，相应的极限位移为11.6 mm，粗砂层的弹性模量E取35 MPa.试桩的桩侧各土层双曲线模型计算参数如表1所示，模型计算参数是赵明华等[20]根据文献［28］的实测数据拟合得到的.

根据前文计算模型参数讨论分析，在利用本文计算方法计算1号试桩桩基沉降时竖向附加应力取0.01，虚土桩分段厚度l取0.5 m，可以满足计算精度要求.图8即是利用本文计算方法所得结果与文献［28］中的实测数据的比较图示，计算结果包括桩顶沉降曲线、桩端沉降曲线以及桩身压缩曲线.

从图8可以看出，利用本文计算方法所得到的桩顶沉降曲线和桩端沉降曲线与文献［28］中的实测曲线有较好的一致性，说明本文竖向受荷桩的沉降计算方法是合理可行的，可以满足工程应用的计算精度要求.另外从计算得到的桩端沉降曲线和桩身压缩曲线可以清晰地看出在不同桩顶荷载作用下桩端沉降和桩身压缩对整个桩顶沉降的贡献比重，当桩顶荷载较小时（如P<5 000 kN），桩顶沉降主要是由桩身压缩造成，桩端土体的沉降量较小；然而随着桩顶荷载的不断增加，传递到桩端的基底压力随之增加，造成桩端土体的压缩量逐渐增大，直到桩端土体沉降量成为桩顶沉降的主要部分.

图9和图10是利用本文计算方法所得到的1号试桩的桩身轴力分布曲线和桩侧摩阻力分布曲线. 从图9可以看出桩身轴力沿桩身逐渐减小，随着桩顶荷载增加传递到桩端的基底压力随之增加.从图10可以看出桩侧摩阻力整体上大致随深度不断增加，其中由于第2层粗砂层不具有黏聚力使得桩侧摩阻力略微减小，相应地，图9中桩身轴力在第2层粗砂层中减小幅度变缓.

5结论

本文采用应力泡形虚土桩模型和考虑桩土界面初始临界摩阻力的改进双曲线模型，提出了一种可以考虑桩端应力扩散效应和桩土界面初始临界摩阻力影响的竖向受荷单桩沉降简化计算方法.通过将利用本文方法得到的计算结果和实测曲线进行对比，表明本文单桩沉降计算方法是合理可行的.

采用本文提出的单桩沉降计算方法可以严格地给出桩端应力扩散边界，并且进一步直观地分析桩端土体性质变化对桩顶沉降的影响.通过对计算模型参数的分析讨论，在利用应力泡形虚土桩模型进行单桩沉降计算时，建议桩端附加应力在0.005P～0.05P之间取值，虚分段精度在0.5～1.5 m之间取值较为合理.

利用本文提出的单桩沉降计算方法可以方便地得到桩顶沉降曲线、桩端沉降曲线、桩身压缩曲线以及桩身轴力和桩侧摩阻力分布曲线，对分析研究竖向受荷桩的沉降性状具有重要意义，同时算法参数明确，计算简单，具有广泛的工程应用前景.

参考文献

[1]SEED H B，REESE L C. The action of soft clay along friction piles [J]. Transactions of the American Society of Civil Engineers，1957，122（1）：731-754.

[2]陳龙珠，梁国钱，朱金颖，等.桩轴向荷载-沉降曲线的一种解析算法[J].岩土工程学报，1994，16（6）：30-38.

CHEN L Z，LIANG G Q，ZHU J Y，et al. Analytical calculation of axial loadig-settlement curve of piles [J]. Chinese Journal of Geotechnical Engineering，1994，16（6）：30-38. （In Chinese）

[3]POULOS H G. Analysis of the settlement of pile groups [J]. Geotechnique，1968，18（4）：449-471.

[4]BUTTERFIELD R，BANERJEE P K. The elastic analysis of compressible piles and pile groups [J]. Geotechnique，1971，21（1）：43-60.

[5]RANDOLPH M F，WROTH C P. An analysis of the vertical deformation of pile groups [J]. Geotechnique，1979，29（4）：423-439.

[6]高盟，高广运，杨成斌，等.层状地基群桩沉降计算的剪切位移解析算法[J].岩土力学，2010，31（4）：1072-1077.

GAO M，GAO G Y，YANG C B，et al. Analytical solution for settlement of group piles in layered ground based on shear displacement method[J]. Rock and Soil Mechanics，2010，31（4）：1072-1077. （In Chinese）

[7]COMODROMOS E M，PAPADOPOULOU M C，RENTZEPERIS I K Pile foundation analysis and design using experimental data and 3D numerical analysis[J]. Computers and Geotechnics，2009，36（5）：819-836.

[8]周红波.桩侧泥皮和桩底沉渣对钻孔桩承载力影响的数值模拟[J].岩土力学，2007，28（5）：956-960.

ZHOU H B Numerical simulation of effects of mud cake and bottom sediment on bearing capacity of bored pile[J]. Rock and Soil Mechanics，2007，28（5）：956-960. （In Chinese）

[9]王奎华，吕述晖，吴文兵，等.层状地基中基于虚土桩模型的单桩沉降计算方法[J].工程力学，2013，30（7）：75-83.

WANG K H，LU S H，WU W B，et al. A new calculation method for the settlement of single pile based on virtual soil-pile model in layered soils[J]. Engineering Mechanics，2013，30（7）：75-83. （In Chinese）

[10]王奎华，罗永健，吴文兵，等.层状地基中考虑桩端应力扩散的单桩沉降计算[J].浙江大学学报（工学版），2013，47（3）：472-479

WANG K H，LUO Y J，WU W B，et al. Calculation method for settlement of single pile considering stress dispersion of pile end soil[J]. Journal of Zhejiang University （Engineering Science），2013，47（3）：472-479.（In Chinese）

[11]张乾青，李术才，李利平，等.考虑侧阻软化和端阻硬化的群桩沉降简化算法[J].岩石力学与工程学报，2013，32（3）：615-624.

ZHANG Q Q，LI S C，LI L P，et al. Simplified method for settlement prediction of pile groups considering skin friction softing and end resistance hardening[J]. Chinese Journal of Rock Mechanics and Engineering，2013，32（3）：615-624.（In Chinese）

[12]刘杰，张可能，肖宏彬.考虑桩侧土软化时单桩荷载-沉降关系的解析算法[J].中国公路学报，2003，16（2）：61-64.

LIU J，ZHANG K N，XIAO H B. Analytical method of loadsettlement relation on single pile under pile-side softening[J]. China Journal of Highway and Transport，2003，16（2）：61-64. （In Chinese）

[13]LIU J，XIAO H B，TANG J，et al. Analysis of load-transfer of single pile in layered soil[J].Computers and Geotechnics，2004，31（2）：127-135.

[14]李鏡培，陈浩华，李林，等.楔形单桩与群桩非线性荷载-沉降曲线计算方法[J].哈尔滨工业大学学报，2017，49（12）：102-109.

LI J P，CHEN H H，LI L，et al.A calculation approach for nonlinear load-settlement curve of single tapered pile and tapered pile group [J]. Journal of Harbin Institute of Technology，2017，49 （12）：102-109.（In Chinese）

[15]WANG ZJ，XIEXY，WANG J C. A new nonlinear method for vertical settlement prediction of a single pile and pile groups in layered soils[J]. Computers and Geotechnics，2012，45：118-126.

[16]ZHANG Q Q，ZHANG Z M. A simplified nonlinear approach for single pile settlement analysis[J]. Canadian Geotechnical Journal，2012，49（11）：1256-1266.

[17]WANG L X，WU W B，ZHANG Y P，et al.Nonlinear analysis of single pile settlement based on stress bubble fictitious soil pile model[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2022，46（7）：1187-1204.

[18]殷宗澤，朱泓，许国华.土与结构材料接触面的变形及数学模拟[J].岩土工程学报，1994，16（3）：14-22.

YIN Z Z，ZHU H，XU G H. Numerical simulation of the deformation in the interface between soil and structural material[J]. Chinese Journal of Geotechnical Engineering，1994，16（3）：14 - 22. （In Chinese）

[19]YIN Z Z，ZHU H，XU G H. A study of deformation in the interface between soil and concrete[J]. Computers and Geotechnics，1995，17（1）：75-92.

[20]赵明华，刘苏，尹平保，等.考虑桩土界面初始临界摩阻力影响的基桩沉降计算方法[J].中南大学学报（自然科学版），2013，44（8）：3425-3431.

ZHAO M H，LIU S，YIN P B，et al.Settlement calculation method of piles considering effect of initial critical friction on pile-soil interface [J]. Journal of Central South University （Science and Technology），2013，44（8）：3425-3431.（In Chinese）

[21]ZHANG Q X，ZOU J F. Simplified approach for settlement analysis of vertically loaded pile[J]. Journal of Engineering Mechanics，2017，143（11）：04017124.

[22]WU W B，WANG K H，ZHANG Z Q，et al. A new approach for time effect analysis of settlement for single pile based on virtual soil-pile model[J].Journal of Central South University，2012，19 （9）：2656-2662.

[23]WU W B，WANG K H，MA S J，et al. Longitudinal dynamic response of pile in layered soil based on virtual soil pile model[J]. Journal of Central South University，2012，19（7）：1999-2007.

[24]WU W B，LIU H，EL NAGGAR M H，et al. Torsional dynamic response of a pile embedded in layered soil based on the fictitious soil pile model [J]. Computers and Geotechnics，2016，80：190-198.

[25]王宁，王奎华，房凯.考虑应力扩散时桩端土对桩体阻抗的影响[J].哈尔滨工业大学学报，2012，44（12）：89-94.

WANG N，WANG K H，FANG K. Influence of stress diffusion of soil surrounding pile tip on the pile impedance[J]. Journal of Harbin Institute of Technology，2012，44（12）：89-94 （In Chinese）

[26]王立兴，吴文兵，杨晓燕，等.层状地基中考虑桩端应力泡形扩散的单桩沉降计算方法[J].土木与环境工程学报，2021，43 （6）：65-73

WANG LX，WU WB，YANG X Y，et al. A calculation method for settlement of single pile considering stress bulb dispersion of pile end soil[J]. Journal of Civil and Environmental Engineering，2021，43（6）：65-73 （In Chinese）

[27]辛公锋.大直径超长桩侧阻软化试验与理论研究[D].杭州：浙江大学，2006.

XIN G F Test and theory study on shaft resistance softening of large diameter and super-long piles[D]. Hangzhou：Zhejiang University，2006 （In Chinese）

[28]洪毓康，陈强华.钻孔灌注桩的荷载传递性能[J].岩土工程学报，1985，7（5）：22-35

HONG Y K，CHEN Q H. Load transfer behaviour of bored piles [J]. Chinese Journal of Geotechnical Engineering，1985，7（5）：22-35 （In Chinese）