成層土中豎向受荷樁沉降性狀簡化計(jì)算方法

王立興 吳文兵 梅國雄 聞敏杰 徐美娟

摘要：基于應(yīng)力泡形虛土樁模型和考慮樁土界面初始臨界摩阻力的改進(jìn)雙曲線模型，提出了一種可以考慮樁端應(yīng)力擴(kuò)散效應(yīng)和樁土界面初始臨界摩阻力影響的豎向受荷單樁沉降簡化計(jì)算方法.首先，利用該計(jì)算方法分析討論了模型參數(shù)對單樁沉降的影響程度，并給出了合理的取值范圍；其次，將用該計(jì)算方法得到的結(jié)果與現(xiàn)場實(shí)測數(shù)據(jù)對比分析，驗(yàn)證了本文計(jì)算方法的合理性，并進(jìn)一步分析研究了不同樁頂荷載作用下豎向受荷樁的樁身壓縮與樁端土體沉降所占比重關(guān)系、樁身軸力分布規(guī)律以及樁側(cè)摩阻力發(fā)揮特性等.本文的計(jì)算方法參數(shù)明確，應(yīng)用簡便，可以為樁基礎(chǔ)設(shè)計(jì)提供更為準(zhǔn)確的沉降值，具有廣泛的工程應(yīng)用前景.

關(guān)鍵詞：豎向受荷樁；沉降性狀；虛土樁模型；雙曲線模型；樁身軸力；樁側(cè)摩阻力

中圖分類號：TU473.1文獻(xiàn)標(biāo)志碼：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

樁基礎(chǔ)由于承載力高、穩(wěn)定性好、適應(yīng)性強(qiáng)等特點(diǎn)，已經(jīng)被廣泛應(yīng)用于高層建筑、橋梁公路等構(gòu)筑物的基礎(chǔ)建設(shè)中.在樁基設(shè)計(jì)中，豎向受荷樁的承載機(jī)理和沉降特性分析是至關(guān)重要的.目前，基樁沉降計(jì)算方法有荷載傳遞法[1-2]、彈性理論法[3-4]、剪切位移法[5-6]和數(shù)值計(jì)算法[7-8]等，其中荷載傳遞法是Seed和Reese[1]基于樁身荷載傳遞機(jī)理提出的，可以靈活地選取不同的荷載傳遞函數(shù)模擬樁與樁側(cè)土非線性相互作用特性及地基土的層狀特性，得到了廣泛應(yīng)用.荷載傳遞法的準(zhǔn)確性主要依賴于荷載傳遞函數(shù)的合理性，荷載傳遞函數(shù)即剪應(yīng)力與樁土相對位移之間的關(guān)系，傳統(tǒng)的荷載傳遞函數(shù)主要有理想彈塑性模型[9-10]、雙折線硬化模型[2，11]、三折線模型[12-13]、雙曲線模型[14-15]及側(cè)阻軟化模型[11，16-17].然而土與混凝土結(jié)構(gòu)面的剪切試驗(yàn)[18-19]發(fā)現(xiàn)當(dāng)剪應(yīng)力小于某一個臨界值時無相對位移產(chǎn)生，表明土與混凝土結(jié)構(gòu)物之間應(yīng)當(dāng)存在初始臨界摩阻力，趙明華等[20]和Zhang等[21]分別提出了考慮初始臨界摩阻力的改進(jìn)雙曲線模型和側(cè)阻軟化模型用于樁基沉降計(jì)算，通過實(shí)例分析驗(yàn)證了考慮初始臨界摩阻力的非線性荷載傳遞模型可以更準(zhǔn)確地模擬樁土相互作用特性，具有更好的工程適用性.

另外，傳統(tǒng)的荷載傳遞法進(jìn)行基樁沉降計(jì)算時簡化地用折線硬化模型[11，21]或雙曲線模型[17，20]來表征樁端荷載與位移的關(guān)系，這樣無法考慮樁端土體的層狀特性和樁端應(yīng)力擴(kuò)散效應(yīng).為此，王奎華等[9-10]、Wu等[22]提出虛土樁模型來模擬樁端土對樁的支撐作用，并將虛土樁模型應(yīng)用于層狀地基中的單樁沉降計(jì)算分析和基樁動力特性研究[23-25].進(jìn)一步，王立興等[26]將應(yīng)力泡形虛土樁模型應(yīng)用到單樁沉降分析來考慮樁端應(yīng)力擴(kuò)散效應(yīng)的影響，但是文中采用雙折線函數(shù)模擬樁側(cè)阻力具有一定局限性，因此，有必要進(jìn)一步對采用側(cè)阻非線性荷載傳遞模型結(jié)合端阻應(yīng)力泡形虛土樁模型的成層土中豎向受荷樁的沉降特性展開研究.

綜上，本文采用應(yīng)力泡形虛土樁模型和考慮樁土界面初始臨界摩阻力的改進(jìn)雙曲線模型，分別模擬樁端土體對樁體的支撐作用和樁側(cè)摩阻力與樁土相對位移間的關(guān)系，并在上述兩種模型的基礎(chǔ)上提出一種考慮樁土界面初始臨界摩阻力影響和樁端應(yīng)力擴(kuò)散效應(yīng)的成層土中豎向受荷樁的沉降特性簡化迭代計(jì)算方法.通過工程實(shí)例對比驗(yàn)證了本文提出計(jì)算方法的合理性，并進(jìn)一步分析研究了不同樁頂荷載作用下，豎向受荷樁的樁頂荷載沉降曲線、樁身壓縮與樁端沉降所占比重變化情況、樁身軸力分布規(guī)律以及樁側(cè)摩阻力發(fā)揮特性等單樁沉降性狀.

1計(jì)算模型及基本假設(shè)

1.1樁側(cè)荷載傳遞模型

為考慮樁土界面初始臨界摩阻力和樁土非線性相互作用，本文采用趙明華等[20]提出的改進(jìn)雙曲線荷載傳遞模型來模擬樁側(cè)摩阻力與樁土相對位移之間的關(guān)系，如圖1所示，該模型適用于正常固結(jié)黏土、中密砂與松砂等土體類型，其函數(shù)表達(dá)式為：

式中：τ（z）是深度z處的樁側(cè)摩阻力；S（z）是深度z處的樁土相對位移；參數(shù)a、b為模型參數(shù)，其物理意義分別為曲線原點(diǎn)處的切線斜率的倒數(shù)和樁側(cè)極限摩阻力τ的倒數(shù)；c（z）為深度z處樁土界面初始臨界靜摩阻力.參數(shù)a、b和c（z）的值可通過現(xiàn)場載荷試驗(yàn)結(jié)果的反分析確定或通過室內(nèi)土-結(jié)構(gòu)接觸面剪切試驗(yàn)獲得.

1.2樁端泡形虛土樁模型

辛公鋒[27]通過對4根現(xiàn)場試樁實(shí)測數(shù)據(jù)的反演，發(fā)現(xiàn)樁在極限承載力下樁底土體大部分仍處于彈性狀態(tài)，因此，可以采用Boussinesq解計(jì)算樁端附加應(yīng)力范圍，并利用彈性力學(xué)知識進(jìn)行樁端土體沉降計(jì)算.假設(shè)樁端對土體的作用力為均布荷載形式，在樁端底面建立如圖2所示的極坐標(biāo)系，將Boussinesq解進(jìn)行二重積分可以求得圓形均布荷載作用下樁端土體中任意一點(diǎn)的豎向應(yīng)力分布為[25-26]：

對式（2）進(jìn)行數(shù)值求解可以得到一組豎向附加應(yīng)力等值線，將等值線繞z軸旋轉(zhuǎn)一周就形成三維附加應(yīng)力等值面，假設(shè)應(yīng)力等值面外的土體內(nèi)附加應(yīng)力很小，將等值面內(nèi)的土體看作受荷主體進(jìn)行計(jì)算，如圖3所示.

1.3樁土相互作用模型及相關(guān)假設(shè)

樁土相互作用體系簡圖如圖3所示，根據(jù)地基土層的成層性情況將實(shí)體樁和虛土樁一共分為n段，其中實(shí)體樁分成m段，虛土樁分成n-m段.從實(shí)體樁頂部向下依次編號為1，2，…，i，…，m，m+1，…，n；各樁段厚度分別為l，l，…，l，…，l.實(shí)體樁總長為L，虛土樁總長為L，各段樁的埋深分別為z，z，…，z，…，z.樁頂作用的豎向荷載為P，產(chǎn)生的樁頂位移為S，荷載向下傳遞，假設(shè)第i段樁的樁頂軸力為P，相應(yīng)的樁頂位移為S，樁側(cè)摩阻力為τ.

為簡化計(jì)算，本文樁土體系相關(guān)假設(shè)有：

1）實(shí)體樁樁身截面均勻，虛土樁以應(yīng)力泡的形式存在，其邊界由式（2）求得；

2）虛土樁的力學(xué)行為與實(shí)體樁相符合，參數(shù)取土體實(shí)際參數(shù)；

3）虛土樁的變形僅考慮彈性狀態(tài)，虛土樁底部為剛性支撐不發(fā)生變形；

4）實(shí)體樁與虛土樁的各樁段界面之間為完全連續(xù)接觸，滿足連續(xù)性條件；

5）地基中同一土層的物理力學(xué)性質(zhì)沿深度方向不變.

2單樁沉降計(jì)算流程

單樁的沉降由樁端土體的沉降和實(shí)心樁體的壓縮兩部分組成.本文采用應(yīng)力泡形虛土樁模型來模擬樁底土體對實(shí)體樁的支撐作用，采用考慮初始臨界摩阻力的雙曲線荷載傳遞模型來模擬樁與樁周土體的相互作用，實(shí)體樁和虛土樁的壓縮沉降量之和，即為所求實(shí)體樁的樁頂總沉降量.

首先推導(dǎo)虛土樁的沉降計(jì)算公式，根據(jù)靜力平衡條件，第i（m<i<n）段樁位移應(yīng)滿足式（3）：

式中：E是第i段虛土樁的彈性模量；A是第i段樁的截面積，A=πR：各個虛土樁段的半徑R由式（2）確定；f是第i段樁周單位摩阻力，滿足式（4）：

f=λS（4）

式中：λ為第i段樁彈性剪切剛度系數(shù)，且λ=τ/S.τ和S分別為第i段樁的極限側(cè)摩阻力和相應(yīng)的極限位移，由原位靜載試驗(yàn)獲得.

由各樁段分界面的軸力和位移的連續(xù)性條件，可以得到式（5）：

式中：P為第i+1段樁樁頂軸力；S為第i+1段樁樁頂位移；z為第i段樁底面到樁頂?shù)木嚯x，即各段樁的埋深.

第i段樁樁頂軸力和位移可由式（6）求解得到：

進(jìn)一步，由式（7）可得第i段樁樁頂剛度為：

由虛土樁端不發(fā)生變形，可以得到：

式（7）（8）組合可以得到第n段樁樁頂剛度：

根據(jù)上述推導(dǎo)，完整的基樁沉降預(yù)測流程可分為以下10個步驟：

1）將樁-土體系劃分為n段，如圖3所示；

2）假設(shè)第n段虛土樁的壓縮位移為S，則第n段虛土樁的軸力P可通過式（9）計(jì)算；

3）利用式（6）可遞歸得到每段虛土樁截面頂部的軸力和位移，實(shí)體樁端的軸力和位移等于虛土樁樁頂?shù)妮S力和位移；

4）假設(shè)第m段實(shí)體樁的壓縮位移為S，并將第m段實(shí)體樁樁頂位移值（S+S/2）代入式（1），就可得到第m段實(shí)體樁的樁周摩阻力τ；

5）利用式（10）計(jì)算第m段實(shí)心樁的軸力增量P：

P=lUτ（10）

式中：l是第m段實(shí)體樁的長度，U是第m段實(shí)體樁的截面周長，U=2πr，r是實(shí)體樁的半徑；

6）利用式（11）計(jì)算第m段樁的樁身壓縮S′：

7）將第m段實(shí)體樁更新后的壓縮位移S′和步驟4）中的壓縮位移S值比較，如果S′與S之差不在允許公差（如1×10）內(nèi)，用S′替換S，然后重復(fù)步驟4）到步驟6），直到兩者之差位于允許公差范圍內(nèi)；

8）將計(jì)算得到的位移增量S和軸力增量P代入式（12），得到第m段實(shí)體樁頂部的軸力P和位移S：

9）對于第i段實(shí)體樁，重復(fù)步驟4）到步驟8），得

各實(shí)體樁樁頂軸力P和位移S，直至得到第1段樁樁頂軸力（即樁頂豎向荷載）P和位移S；

10）假定一組遞增的第n段虛土樁壓縮位移為S，然后重復(fù)步驟2）到步驟9），便可得到一系列荷載位移值，最終得到樁頂荷載沉降曲線.

3計(jì)算模型參數(shù)討論

樁端土體的沉降對樁頂沉降有較大的貢獻(xiàn)，為了突出分析端應(yīng)力擴(kuò)散效應(yīng)對樁頂沉降的影響，下面分析中取單層均質(zhì)地基，著重對泡形虛土樁模型參數(shù)進(jìn)行討論，對于成層地基情況，只需要根據(jù)實(shí)際土層情況進(jìn)行賦值即可.假定樁土體系參數(shù)取值為：實(shí)體樁長L=20 m，樁身半徑r=0.5 m，樁身彈性模量E=3×10MPa.地基土層為均質(zhì)硬黏土層，土層的極限側(cè)摩阻力τ=120 kPa，相應(yīng)的極限位移s=8 mm.土層的雙曲線模型計(jì)算參數(shù)a=0.03 mm/kPa，b=0.007/kPa，c=10 kPa.模型參數(shù)是參考本文工程案例的實(shí)際參數(shù)選取的合理近似值.

3.1豎向附加應(yīng)力σ

豎向附加應(yīng)力σ是采用應(yīng)力泡形虛土樁模型分析單樁沉降的一個重要參數(shù)，選取不同樁端土體的豎向附加應(yīng)力σ代表選取了不同長度（深度）的泡形虛土樁模型，進(jìn)一步將豎向附加應(yīng)力σ和任意深度值z代入式（2）進(jìn)行數(shù)值求解可以得到不同深度處的徑向?qū)挾?不同邊界的應(yīng)力泡形虛土樁模型如圖4 所示，基于本文計(jì)算方法豎向附加應(yīng)力σ對樁頂沉降的影響如圖5所示.

從圖5可以看出當(dāng)樁頂荷載較?。≒<5 000 kN）時，樁端豎向附加應(yīng)力σ對樁頂沉降影響很小，這是由于當(dāng)樁頂荷載較小時傳遞到樁端的基底壓力和樁端土體中的豎向附加應(yīng)力相應(yīng)較小，使得不同豎向附加應(yīng)力對應(yīng)的泡形虛土樁產(chǎn)生的壓縮變形差值就會很小；當(dāng)樁頂荷載繼續(xù)增大時，其對樁頂沉降的影響程度會明顯增加，相同荷載下隨著樁端豎向附加應(yīng)力的減小樁頂沉降值隨之減小，這是由于豎向附加應(yīng)力的減小意味著泡形虛土樁的邊界增大，即承受基底壓力的受荷土體體積增大，那么整體的樁端土體沉降值會減小，樁頂沉降也隨之減小；另外隨著豎向附加應(yīng)力的增大，樁頂沉降增大的幅度不斷減小并且最后收斂于σ=0.01 P附近，因此，在利用應(yīng)力泡形虛土樁模型進(jìn)行單樁沉降計(jì)算時建議樁端豎向附加應(yīng)力取值范圍在0.005P～0.05 P較為合理.

3.2虛土樁分段精度

進(jìn)一步討論虛土樁的分段精度對樁頂沉降計(jì)算的影響.本文中虛土樁分段精度用虛土樁分段厚度l來進(jìn)行量化分析，基于本文計(jì)算方法虛土樁分段厚度l對樁頂沉降的影響如圖6所示.

從圖6可以看出，當(dāng)樁頂荷載較小（P<6 000 kN）時，虛土樁分段精度l對樁頂沉降影響很小；隨著樁頂荷載的增大虛土樁分段精度l對樁頂沉降影響增加，并且隨著分段精度l的增大，其對樁頂沉降影響幅度不斷減小并且最后收斂于l=1m附近，因此，在利用應(yīng)力泡形虛土樁模型進(jìn)行單樁沉降計(jì)算時建議虛土樁分段精度%，取值在0.5～1.5 m較為合理.

3.3虛土樁彈性模量

利用土體彈性模量E的變化分析樁端土體性質(zhì)對樁頂沉降的影響，基于本文計(jì)算方法樁端土體彈性模量E對樁頂沉降的影響如圖7所示.

從圖7可以看出，當(dāng)樁頂荷載較?。≒<6 000 kN）時，樁端土彈性模量E對樁頂沉降影響很??；隨著樁頂荷載的增大樁端土彈性模量E對樁頂沉降影響程度增加，具體表現(xiàn)為樁端土彈性模量E越大，樁頂沉降量越小，這是因?yàn)闃抖送翉椥阅Ａ縀越大表示樁端土體性質(zhì)越好，抗壓縮變形能力越強(qiáng)，因此整體樁身沉降量會減小.

4工程實(shí)例分析

為了驗(yàn)證本文所述基樁沉降計(jì)算方法的準(zhǔn)確性，利用本文計(jì)算方法對文獻(xiàn)［28］中的1號試樁進(jìn)行沉降計(jì)算，將計(jì)算結(jié)果與文獻(xiàn)［28］中的實(shí)測曲線進(jìn)行比較以驗(yàn)證本文計(jì)算方法的合理性，并進(jìn)一步利用本文計(jì)算方法對豎向受荷樁進(jìn)行沉降性狀分析.

試驗(yàn)場地地勢平坦，地基土層主要由沖洪積相的亞黏土、黏土、粗砂組成.1號試樁的實(shí)測樁徑為1.09 m，樁身彈性模量E為2.43×10 MPa，樁長為25 m，樁端持力層為粗砂層，根據(jù)文獻(xiàn)［28］的實(shí)測數(shù)據(jù)，粗砂層的極限側(cè)摩阻力為183 kPa，相應(yīng)的極限位移為11.6 mm，粗砂層的彈性模量E取35 MPa.試樁的樁側(cè)各土層雙曲線模型計(jì)算參數(shù)如表1所示，模型計(jì)算參數(shù)是趙明華等[20]根據(jù)文獻(xiàn)［28］的實(shí)測數(shù)據(jù)擬合得到的.

根據(jù)前文計(jì)算模型參數(shù)討論分析，在利用本文計(jì)算方法計(jì)算1號試樁樁基沉降時豎向附加應(yīng)力取0.01，虛土樁分段厚度l取0.5 m，可以滿足計(jì)算精度要求.圖8即是利用本文計(jì)算方法所得結(jié)果與文獻(xiàn)［28］中的實(shí)測數(shù)據(jù)的比較圖示，計(jì)算結(jié)果包括樁頂沉降曲線、樁端沉降曲線以及樁身壓縮曲線.

從圖8可以看出，利用本文計(jì)算方法所得到的樁頂沉降曲線和樁端沉降曲線與文獻(xiàn)［28］中的實(shí)測曲線有較好的一致性，說明本文豎向受荷樁的沉降計(jì)算方法是合理可行的，可以滿足工程應(yīng)用的計(jì)算精度要求.另外從計(jì)算得到的樁端沉降曲線和樁身壓縮曲線可以清晰地看出在不同樁頂荷載作用下樁端沉降和樁身壓縮對整個樁頂沉降的貢獻(xiàn)比重，當(dāng)樁頂荷載較小時（如P<5 000 kN），樁頂沉降主要是由樁身壓縮造成，樁端土體的沉降量較??；然而隨著樁頂荷載的不斷增加，傳遞到樁端的基底壓力隨之增加，造成樁端土體的壓縮量逐漸增大，直到樁端土體沉降量成為樁頂沉降的主要部分.

圖9和圖10是利用本文計(jì)算方法所得到的1號試樁的樁身軸力分布曲線和樁側(cè)摩阻力分布曲線. 從圖9可以看出樁身軸力沿樁身逐漸減小，隨著樁頂荷載增加傳遞到樁端的基底壓力隨之增加.從圖10可以看出樁側(cè)摩阻力整體上大致隨深度不斷增加，其中由于第2層粗砂層不具有黏聚力使得樁側(cè)摩阻力略微減小，相應(yīng)地，圖9中樁身軸力在第2層粗砂層中減小幅度變緩.

5結(jié)論

本文采用應(yīng)力泡形虛土樁模型和考慮樁土界面初始臨界摩阻力的改進(jìn)雙曲線模型，提出了一種可以考慮樁端應(yīng)力擴(kuò)散效應(yīng)和樁土界面初始臨界摩阻力影響的豎向受荷單樁沉降簡化計(jì)算方法.通過將利用本文方法得到的計(jì)算結(jié)果和實(shí)測曲線進(jìn)行對比，表明本文單樁沉降計(jì)算方法是合理可行的.

采用本文提出的單樁沉降計(jì)算方法可以嚴(yán)格地給出樁端應(yīng)力擴(kuò)散邊界，并且進(jìn)一步直觀地分析樁端土體性質(zhì)變化對樁頂沉降的影響.通過對計(jì)算模型參數(shù)的分析討論，在利用應(yīng)力泡形虛土樁模型進(jìn)行單樁沉降計(jì)算時，建議樁端附加應(yīng)力在0.005P～0.05P之間取值，虛分段精度在0.5～1.5 m之間取值較為合理.

利用本文提出的單樁沉降計(jì)算方法可以方便地得到樁頂沉降曲線、樁端沉降曲線、樁身壓縮曲線以及樁身軸力和樁側(cè)摩阻力分布曲線，對分析研究豎向受荷樁的沉降性狀具有重要意義，同時算法參數(shù)明確，計(jì)算簡單，具有廣泛的工程應(yīng)用前景.

參考文獻(xiàn)

[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].巖土工程學(xué)報(bào)，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]高盟，高廣運(yùn)，楊成斌，等.層狀地基群樁沉降計(jì)算的剪切位移解析算法[J].巖土力學(xué)，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]周紅波.樁側(cè)泥皮和樁底沉渣對鉆孔樁承載力影響的數(shù)值模擬[J].巖土力學(xué)，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ì)算方法[J].工程力學(xué)，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]王奎華，羅永健，吳文兵，等.層狀地基中考慮樁端應(yīng)力擴(kuò)散的單樁沉降計(jì)算[J].浙江大學(xué)學(xué)報(bào)（工學(xué)版），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]張乾青，李術(shù)才，李利平，等.考慮側(cè)阻軟化和端阻硬化的群樁沉降簡化算法[J].巖石力學(xué)與工程學(xué)報(bào)，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]劉杰，張可能，肖宏彬.考慮樁側(cè)土軟化時單樁荷載-沉降關(guān)系的解析算法[J].中國公路學(xué)報(bào)，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ì)算方法[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào)，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]殷宗澤，朱泓，許國華.土與結(jié)構(gòu)材料接觸面的變形及數(shù)學(xué)模擬[J].巖土工程學(xué)報(bào)，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ì)算方法[J].中南大學(xué)學(xué)報(bào)（自然科學(xué)版），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]王寧，王奎華，房凱.考慮應(yīng)力擴(kuò)散時樁端土對樁體阻抗的影響[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2012，44（12）：89-94.

WANG N，WANG K H，F(xiàn)ANG 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]王立興，吳文兵，楊曉燕，等.層狀地基中考慮樁端應(yīng)力泡形擴(kuò)散的單樁沉降計(jì)算方法[J].土木與環(huán)境工程學(xué)報(bào)，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]辛公鋒.大直徑超長樁側(cè)阻軟化試驗(yàn)與理論研究[D].杭州：浙江大學(xué)，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]洪毓康，陳強(qiáng)華.鉆孔灌注樁的荷載傳遞性能[J].巖土工程學(xué)報(bào)，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）