王立興 吳文兵 楊曉燕 張?jiān)迄i 王奎華

摘 要：基于虛土樁模型，研究了樁端應(yīng)力擴(kuò)散效應(yīng)對(duì)單樁沉降的影響。結(jié)合虛土樁模型和Boussinesq解，提出考慮樁端應(yīng)力泡形擴(kuò)散的應(yīng)力泡形虛土樁模型來考慮樁端土對(duì)樁的支承作用;利用荷載傳遞法，對(duì)樁周土采用雙折線模型，推導(dǎo)得到層狀地基中考慮樁端應(yīng)力泡性擴(kuò)散的單樁沉降解析解;討論了彈性極限位移、彈性抗剪切剛度系數(shù)、附加應(yīng)力值等參數(shù)對(duì)單樁沉降的影響;結(jié)合工程實(shí)測(cè)數(shù)據(jù)，將該方法與現(xiàn)有理論解進(jìn)行對(duì)比，驗(yàn)證了方法的合理性。應(yīng)力泡形虛土樁模型的最大優(yōu)勢(shì)在于能夠較嚴(yán)格地給出樁端應(yīng)力擴(kuò)散邊界，具有更廣泛的工程應(yīng)用價(jià)值。

關(guān)鍵詞：?jiǎn)螛冻两?應(yīng)力泡形虛土樁;層狀地基;應(yīng)力擴(kuò)散;荷載傳遞法

中圖分類號(hào)：TU473.1 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：2096-6717（2021）06-0065-09

Abstract： Based on the fictitious soil pile model， the effect of pile-end stress diffusion on the settlement of a single pile has been investigated theoretically. Firstly， a pile-end stress-bubble fictitious soil pile model considering the spread of the pile-end stress-bubble was proposed to simulate the supporting effect of the pile end soil on the pile body by combining the fictitious soil pile model and Boussinesq solution. Then， based on the load transfer method， an analytical solution of single pile settlement considering pile-end stress-bubble diffusion in layered foundation was derived by utilizing the double-broken-line model for the soil around the pile. Furthermore， the effects of parameters such as elastic limit displacement， elastic shear stiffness coefficient and additional stress value on the settlement of single pile were discussed based on the present solution. Finally， combining with the actual engineering data， the present method was validated by comparing with existing theoretical solution. The stress-bubble fictitious soil pile model is more meaningful for engineering applications on account of the obvious advantage that the stress diffusion boundary at pile end can be given more strictly.

Keywords：single pile settlement; stress-bubble fictitious soil pile; layered soils; stress diffusion; load transfer method

樁基沉降計(jì)算是樁基工程設(shè)計(jì)中非常重要的部分。學(xué)者們對(duì)樁基沉降進(jìn)行了廣泛而深入的研究，得到了很多有價(jià)值的單樁沉降計(jì)算方法，如荷載傳遞法[1-7]、剪切位移法[8-10]、彈性理論法[11-13]、規(guī)范法[14]和數(shù)值計(jì)算法[15-17]等，其中，荷載傳遞法因可通過不同的樁土相互作用模型來反映樁土共同作用、地基層狀特性而被學(xué)者廣泛運(yùn)用?？偟膩碚f，基于荷載傳遞法的樁土相互作用模型主要有理想彈塑性模型[1，4-6]、雙折線硬化模型[3，18]、三折線模型[19]及雙曲線模型[7]，通過選取不同的樁土相互作用模型，可以較準(zhǔn)確地得到不同土性中單樁的沉降，具有較好的工程適用性。

上述研究采用荷載傳遞法計(jì)算樁基沉降時(shí)，均將樁端土對(duì)樁的支撐作用簡(jiǎn)化為彈簧模型，彈簧參數(shù)可結(jié)合樁端土硬化特性、非線性特性等根據(jù)實(shí)測(cè)的樁端反力位移曲線反演確定。然而，這種彈簧支撐模型存在一些不足：無法考慮樁端土體的層狀特性;無法建立樁端土層物理力學(xué)性質(zhì)與彈簧參數(shù)之間的相互關(guān)系;且無法反映持力層內(nèi)及以下土體性質(zhì)的變化對(duì)單樁沉降的影響[20]。為此，Wu等[20-21]提出了虛土樁模型，并用于樁基縱向和扭轉(zhuǎn)動(dòng)力特性分析。隨后，王奎華等[22]率先將虛土樁模型應(yīng)用于層狀地基中的單樁沉降計(jì)算分析，并提出采用錐型虛土樁模型來考慮樁端應(yīng)力擴(kuò)散效應(yīng)對(duì)單樁沉降的影響[23]，辛冬冬等[24]將虛土樁模型推廣到層狀地基中的群樁沉降分析。盡管虛土樁模型能很好地考慮樁端層狀特性對(duì)基樁靜動(dòng)力學(xué)特性的影響，且能在一定程度上考慮樁端應(yīng)力擴(kuò)散效應(yīng)，但虛土樁的合理長(zhǎng)度及應(yīng)力擴(kuò)散影響范圍到底多大仍不得而知。

為了更好地分析樁端土性質(zhì)對(duì)單樁沉降的影響，筆者基于Boussinesq解提出了一種應(yīng)力泡形虛土樁模型來考慮樁端應(yīng)力擴(kuò)散效應(yīng)，并基于荷載傳遞法對(duì)樁側(cè)土采用雙折線模型、對(duì)樁端土采用應(yīng)力泡形虛土樁模型，考慮樁端土中實(shí)際附加應(yīng)力的分布情況，推導(dǎo)得到了成層地基中的單樁樁頂沉降計(jì)算公式。

1 數(shù)學(xué)模型及假設(shè)條件

1.1 應(yīng)力泡形虛土樁模型的建立

根據(jù)Boussinesq解可知，樁端土體中附加應(yīng)力大致的分布規(guī)律為：在徑向上隨著徑向半徑的增加逐漸減小，在豎向上隨著深度的增加也逐漸減小，在整體空間中呈橢球型泡狀分布，即應(yīng)力泡。假設(shè)樁端對(duì)土體的作用力為均布荷載形式，根據(jù)Boussinesq解可以求出樁端土中應(yīng)力分布形式的解析式，然后建立應(yīng)力泡形虛土樁模型。建立如圖1所示的軸對(duì)稱坐標(biāo)系，基于Boussinesq解，通過積分可以求得圓形均布荷載作用下地基土體中任意一點(diǎn)的豎向應(yīng)力分布[25]。

式中：σ-為附加應(yīng)力分布線;σz為地基中任意點(diǎn)的豎向應(yīng)力;P為樁端平均應(yīng)力，即均布荷載;z為任意點(diǎn)的深度，即距樁端的豎向距離;R為任意點(diǎn)的徑向距離，即距均布荷載的水平距離;r0為樁體半徑，即圓形荷載的作用半徑;r為均布荷載微單元距離均布荷載中心（即樁軸心）的距離;θ為均布荷載微單元與任意點(diǎn)在水平面上投影的夾角。

將式（1）求解得到的相同數(shù)值的豎向應(yīng)力等值線繞z軸旋轉(zhuǎn)一周就形成一系列以z軸為中心的三維應(yīng)力等值面，將等值面以內(nèi)的土體看作是受荷主體，即為應(yīng)力泡形虛土樁的主體部分，如圖2所示。

1.2 樁土相互作用模型的建立

樁周土采用雙折線模型，如圖3所示，荷載傳遞函數(shù)為

式中：u為樁周土的位移;ub為樁周土的彈性極限位移;fs為單位長(zhǎng)度樁周側(cè)摩阻力;λ為理想彈塑性模型的彈性抗剪切剛度系數(shù)，kN/m2。

樁端土采用應(yīng)力泡形虛土樁模型，把樁端至基巖之間的應(yīng)力泡形土體看成“土樁”，即所謂虛土樁，其參數(shù)分別取各實(shí)際土層參數(shù)，而變形按類似于樁的平面假定。根據(jù)樁端土層成層性情況及樁端土體性質(zhì)和受力環(huán)境是否變化而把實(shí)體樁和虛土樁一共分為n（n>m）段，其中實(shí)體樁分成m段，虛土樁分成N（即N=n-m）段，從實(shí)體樁頂部向下依次編號(hào)為1、2、…、i、…、m、m+1、…、n，各樁段厚度分別為l1、l2、…、li、…、ln，實(shí)體樁總長(zhǎng)為L(zhǎng)1，虛土樁總長(zhǎng)為L(zhǎng)2，各段樁的埋深分別為z1、z2、…、zi、…、zn。計(jì)算時(shí)，可以將虛土樁各樁段再細(xì)分為很薄的圓柱體進(jìn)行近似計(jì)算，計(jì)算精度由單元數(shù)N控制，計(jì)算中可以直接將樁土體系按照m+N段進(jìn)行計(jì)算。具體的樁體分段與樁土相互作用體系簡(jiǎn)圖如圖2所示。由式（1）可求得虛土樁各深度樁段對(duì)應(yīng)的半徑Ri，即各段實(shí)體樁和虛土樁的半徑分別為

1.3 樁土體系相關(guān)假設(shè)

1）樁土體系呈線彈性狀態(tài)，實(shí)體樁樁身截面均勻，虛土樁以應(yīng)力泡的形式存在。

2）樁體與土體之間沒有相對(duì)滑移，實(shí)體樁與虛土樁的各樁段交界面之間為完全連續(xù)接觸，滿足連續(xù)性條件。

3）樁單元中任意點(diǎn)的位移只與該點(diǎn)的樁周側(cè)摩阻力有關(guān)，虛土樁端部為剛性支撐。

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

5）虛土樁的力學(xué)行為與樁體相符合，參數(shù)為土體參數(shù)。

2 樁頂沉降計(jì)算方程建立及求解

根據(jù)荷載在樁身的傳遞機(jī)理，樁體的截面軸力及位移隨深度的增加，逐漸減小，且隨著樁頂荷載的增加，截面軸力及位移逐漸向樁底發(fā)展。因此，隨著荷載的增加，樁周土由淺入深逐漸進(jìn)入塑性階段。根據(jù)樁土體系所處狀態(tài)的不同，分3個(gè)階段分別給出樁頂荷載和沉降的計(jì)算公式。需要說明的是，由于樁基受力后通常僅對(duì)樁周徑向一定范圍內(nèi)的土體產(chǎn)生影響，這里所指的樁周土體即樁基應(yīng)力影響范圍內(nèi)的土體。

首先推導(dǎo)荷載傳遞函數(shù)基本微分方程，根據(jù)樁上任意一個(gè)單元體的靜力平衡條件，可以得到

式中：E為該樁段彈性模量;A為該樁段截面面積;P為該樁段頂部荷載。

對(duì)式（5）兩端求導(dǎo)，并將式（4）代入，可得荷載傳遞法的基本微分方程

2.1 樁周土全部處于彈性階段

當(dāng)樁頂荷載較小時(shí)，樁周土全部處于彈性階段，根據(jù)式（4）～式（6），此時(shí)第i個(gè)樁段截面位移u（x）應(yīng)滿足方程

式中：Pi為第i段樁體頂部荷載;Si為第i段樁體頂部位移;zi為第i段樁體截面到樁頂?shù)木嚯x;Ei為第i段樁體彈性模量;Ai為第i段樁體橫截面積;λi為第i段樁體彈性抗剪切剛度系數(shù)。

各虛土樁段面積Ai值計(jì)算采用式（8）。

根據(jù)式（12）和式（14）可以遞推得到樁體頂部剛度K1，進(jìn)而根據(jù)樁頂荷載P1，就可以求得樁頂沉降值S1。

式（12）和（14）中所需要的計(jì)算參數(shù)均可以由實(shí)測(cè)的土層參數(shù)直接或間接計(jì)算得到。因此，根據(jù)式（11）～式（15）就可以求得樁頂沉降值。

2.2 樁周土部分處于塑性階段

樁周土部分處于塑性階段時(shí)，樁土的相互作用體系簡(jiǎn)圖如圖4所示，塑性樁段的長(zhǎng)度為L(zhǎng)a，由于塑性樁段的長(zhǎng)度La與樁頂荷載和樁周土的性質(zhì)有關(guān)，在實(shí)際工程應(yīng)用中應(yīng)當(dāng)根據(jù)現(xiàn)場(chǎng)試驗(yàn)確定樁周土處于塑性階段的深度La。

根據(jù)假設(shè)，當(dāng)樁周土部分處于塑性階段時(shí)，可以依據(jù)樁周土所處的狀態(tài)把樁分為兩個(gè)部分，即塑性段和彈性段，先按照式（11）～式（15）計(jì)算彈性樁段的剛度，再按照式（18）計(jì)算塑性樁段的剛度，就可以得到單樁樁頂沉降值。

2.3 樁周土全部處于塑性階段

樁周土全部處于塑性階段時(shí)，如圖4所示，即La=L1+L2?？汕蟮玫趎段樁的樁頂荷載及位移為

當(dāng)樁周土全部處于塑性階段時(shí)，結(jié)合式（18）和式（19）遞推就可以得到單樁樁頂沉降值。

3 計(jì)算參數(shù)的選取及影響

3.1 彈性極限位移ub

對(duì)于彈性極限位移，很多學(xué)者都有研究，目前有經(jīng)驗(yàn)法[26-27]、經(jīng)驗(yàn)公式法[28]及應(yīng)變計(jì)測(cè)試方法[29]。為了在實(shí)際工程應(yīng)用中方便取值，建議采用應(yīng)變計(jì)測(cè)試方法進(jìn)行計(jì)算。

式中：s0為樁頂沉降量;ubi為第i段樁土體系的極限位移;li為第i段樁的長(zhǎng)度;εi為靜載試驗(yàn)中在樁頂作用承載力設(shè)計(jì)值時(shí)第i個(gè)測(cè)試面的應(yīng)變。

通過靜載試驗(yàn)，將加載到極限承載力時(shí)測(cè)試得到的樁頂沉降量和應(yīng)變值代入式（20），就可以求得樁土彈性極限位移。

3.2 彈性抗剪切剛度系數(shù)λ

與彈性極限位移類似，為了在實(shí)際工程應(yīng)用中方便取值，建議彈性抗剪切剛度系數(shù)采用應(yīng)變計(jì)測(cè)試方法進(jìn)行計(jì)算。

式中：τui為第i段樁的極限側(cè)摩阻力;τi與ui為靜載試驗(yàn)?zāi)骋患?jí)荷載作用下計(jì)算得到的第i段樁的側(cè)摩阻力與相對(duì)位移。

在雙折線模型中，λ為斜線段斜率，知道任一點(diǎn)的側(cè)摩阻力和相對(duì)位移就可以計(jì)算得到抗剪切剛度系數(shù)。

3.3 附加應(yīng)力值σz

由式（1）可以看出，附加應(yīng)力值σz是利用應(yīng)力泡形虛土樁模型分析單樁沉降的一個(gè)重要參數(shù)，選取不同的附加應(yīng)力值σz意味著選用了范圍大小不同的應(yīng)力泡形虛土樁模型，如圖5所示。樁端土體的沉降對(duì)樁頂沉降有較大的貢獻(xiàn)，因此，選取范圍不同的應(yīng)力泡形虛土樁模型對(duì)樁頂沉降影響很大。

3.3.1 算例1：樁端土為軟弱土

為了考慮樁端應(yīng)力泡范圍不同的影響，取上硬下軟雙層土作為算例1，實(shí)體樁深度范圍內(nèi)為硬黏土，虛土樁范圍內(nèi)為軟弱土，土層參數(shù)來源于文獻(xiàn)[30]，如表1所示。

樁身彈性模量Ep=3×104 MPa，樁身半徑r0=0.5 m，樁頂荷載P取2 000、3 000、4 000 kN，附加應(yīng)力值σz取0.1P～0.000 01P，虛土樁部分按照0.5 m分層計(jì)算，得到樁頂?shù)某两抵等鐖D6所示。

如圖6所示，在樁側(cè)土為硬土，樁端土為軟土的雙層地基中，附加應(yīng)力范圍對(duì)樁頂沉降的影響不明顯，沉降變化不大。整體上，隨著附加應(yīng)力范圍的增大（即應(yīng)力泡形虛土樁范圍增大），樁頂沉降隨之減小，但減小幅度不大。

3.3.2 算例2：樁端土為硬黏土或砂性土

取上軟下硬雙層土作為算例2，實(shí)體樁深度范圍內(nèi)為軟弱土，虛土樁范圍內(nèi)為硬黏土或砂性土，各參數(shù)取值如表2所示。

樁身彈性模量Ep=3×104 MPa，樁身半徑r0=0.5 m，樁頂荷載P取2 000、3 000、4 000 kN，附加應(yīng)力值σz取0.1P～0.000 01P時(shí)，虛土樁部分按照0.5 m分層計(jì)算，得到樁頂?shù)某两抵等鐖D7所示。

如圖7所示，在樁側(cè)土為軟土，樁端土為硬土的雙層地基中，附加應(yīng)力范圍對(duì)樁頂沉降的影響明顯，尤其是附加應(yīng)力分布線σ-從0.1～0.001變化時(shí)，沉降變化較大，特別是當(dāng)附加應(yīng)力分布線σ-從0.003～0.004變化時(shí)會(huì)有一個(gè)沉降突變，這是由于應(yīng)力泡的影響深度大于樁端土的厚度，即應(yīng)力泡下端部分變成了基巖，而基巖不發(fā)生沉降，因此樁頂沉降突然減小。但整體上，隨著附加應(yīng)力范圍的增大（即應(yīng)力泡形虛土樁范圍增大），樁頂沉降隨之減小，并且當(dāng)附加應(yīng)力范圍的增大到一定值后，樁頂沉降變化幅值變小。

綜合上述分析發(fā)現(xiàn)，當(dāng)樁周土性好時(shí)，樁以側(cè)阻為主，應(yīng)力泡形虛土樁范圍的大小對(duì)樁頂沉降量影響不大;而當(dāng)樁周土性差時(shí)，樁以端承為主，應(yīng)力泡形虛土樁范圍的大小對(duì)樁頂沉降量影響明顯;且當(dāng)應(yīng)力泡的影響深度大于樁端土層厚度時(shí)，會(huì)有沉降值的突變。實(shí)際工程中，應(yīng)針對(duì)不同類型的樁周土及樁端土的性質(zhì)與樁端土的厚度，選擇合適的應(yīng)力泡范圍，當(dāng)樁端土厚度并非過大時(shí)，可以選擇應(yīng)力泡的影響深度等于樁端土厚度時(shí)的應(yīng)力泡作為虛土樁模型的邊界進(jìn)行計(jì)算。

4 工程實(shí)例分析

王奎華等[23]提出了考慮應(yīng)力擴(kuò)散角的樁頂沉降計(jì)算模型，并通過計(jì)算兩根試樁的沉降與實(shí)測(cè)值對(duì)比驗(yàn)證了其實(shí)用性。為驗(yàn)證該方法的合理性及實(shí)用性，對(duì)文獻(xiàn)[23]中的T3、T4樁進(jìn)行擬合，并將擬合結(jié)果與原文實(shí)測(cè)值及其計(jì)算值做對(duì)比分析，與文獻(xiàn)[23]保持一致，采用樁周土處于彈性階段公式計(jì)算。樁周土及樁端土的物理力學(xué)性質(zhì)參數(shù)來源于文獻(xiàn)[23]，如表3所示。

T3試樁：樁直徑為800 mm，樁長(zhǎng)61.8 m;T4試樁：樁直徑為1 000 mm，樁長(zhǎng)73.8m，樁身彈性模量Ep=3.3×107 kPa。采用本文解計(jì)算各樁頂沉降量并與文獻(xiàn)[23]計(jì)算結(jié)果及實(shí)測(cè)值對(duì)比，見圖8及圖9。利用該方法計(jì)算時(shí)，T3樁與T4樁的虛土樁部分取至第15層粉土底部，根據(jù)式（1）計(jì)算虛土樁影響深度剛好為第15層粉土底部時(shí)的T3試樁附加應(yīng)力值為0.000 36P，T4試樁附加應(yīng)力值為0.001 2P，并分別選取此時(shí)的應(yīng)力泡范圍作為虛土樁邊界進(jìn)行計(jì)算。

與兩根試樁的實(shí)測(cè)值與文獻(xiàn)[23]的擴(kuò)散角法計(jì)算值對(duì)比可以發(fā)現(xiàn)，該方法考慮應(yīng)力泡形擴(kuò)散的計(jì)算結(jié)果與文獻(xiàn)[23]考慮應(yīng)力擴(kuò)散角的計(jì)算結(jié)果較為接近，與實(shí)測(cè)結(jié)果具有較好的一致性。與文獻(xiàn)[23]的錐形虛土樁模型相比，該模型的最大優(yōu)勢(shì)在于能夠較嚴(yán)格地給出樁端應(yīng)力擴(kuò)散邊界，因而具有更廣泛的工程應(yīng)用價(jià)值。

5 結(jié)論

1）基于虛土樁模型并結(jié)合Boussinesq解，提出了考慮樁端應(yīng)力擴(kuò)散效應(yīng)的應(yīng)力泡形虛土樁模型，并建立了基于該模型單樁沉降計(jì)算的解析解。該解能夠計(jì)算樁側(cè)土體處于彈性及塑性階段時(shí)的樁頂荷載沉降曲線，參數(shù)明確，應(yīng)用簡(jiǎn)便，可以用于工程現(xiàn)場(chǎng)單樁沉降量的計(jì)算分析，為樁基礎(chǔ)設(shè)計(jì)提供較為準(zhǔn)確的沉降值。

2）當(dāng)樁周土性較好時(shí)，樁側(cè)摩阻力發(fā)揮主要作用，樁周土承擔(dān)主要荷載，應(yīng)力泡范圍的大小對(duì)樁頂沉降量影響不明顯，可取較小范圍的應(yīng)力泡形虛土樁;當(dāng)樁周土性較差時(shí)，樁端阻力發(fā)揮主要作用，樁端土承擔(dān)主要荷載，應(yīng)力泡范圍的大小對(duì)樁頂沉降量影響明顯，應(yīng)針對(duì)不同類型的樁周土及樁端土的性質(zhì)與樁端土的厚度，選擇合適的應(yīng)力泡范圍。

3）通過分析不同附加應(yīng)力取值對(duì)單樁樁頂沉降量的影響發(fā)現(xiàn)，當(dāng)應(yīng)力泡的影響深度大于樁端土厚度時(shí)，會(huì)發(fā)生沉降量的突變。因此，在工程應(yīng)用中，可以選取應(yīng)力泡的影響深度恰好是樁端土厚度時(shí)的應(yīng)力泡形虛土樁模型進(jìn)行計(jì)算。

4）該方法與錐形虛土樁模型對(duì)兩根試樁的沉降計(jì)算結(jié)果與實(shí)測(cè)結(jié)果對(duì)比表明，該方法的計(jì)算結(jié)果與實(shí)測(cè)結(jié)果有較好的一致性。同時(shí)，該模型的最大優(yōu)勢(shì)在于能夠較嚴(yán)格地給出樁端應(yīng)力擴(kuò)散邊界，因而具有更廣泛的工程應(yīng)用價(jià)值。

參考文獻(xiàn)：

[1]BOLTON SEED H， 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]KRAFT L M， RAY R P， KAKAAKI T. Theoretical t-z curves[J]. Journal of Geotechnical Engineering Division， ASCE， 1981， 107（11）： 1543-1561.

[3]陳龍珠， 梁國(guó)錢， 朱金穎， 等. 樁軸向荷載沉降曲線的一種解析算法[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）

[4]趙明華， 鄒丹， 鄒新軍. 群樁沉降計(jì)算的荷載傳遞法[J]. 工程力學(xué)， 2006， 23（7）： 119-123.

ZHAO M H， ZOU D， ZOU X J. Settlement calculation of pile groups by load transfer method[J]. Engineering Mechanics， 2006， 23（7）： 119-123. （in Chinese）

[5]徐禮閣， 王奎華， 張鵬， 等. 層狀地基中考慮樁側(cè)水泥土與土的靜鉆根植樁單樁沉降計(jì)算[J]. 中南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2016， 47（3）： 868-874.

XU L G， WANG K H， ZHANG P， et al. Calculation method for settlement of single static drill rooted pile considering cemented soil and soil around pile in layered soil[J]. Journal of Central South University （Science and Technology）， 2016， 47（3）： 868-874. （in Chinese）

[6]周佳錦， 龔曉南， 王奎華， 等. 層狀地基中靜鉆根植竹節(jié)樁單樁沉降計(jì)算[J]. 巖土力學(xué)， 2017， 38（1）： 109-116.

ZHOU J J， GONG X N， WANG K H， et al. A simplified approach to calculating settlement of a single pre-bored grouting planted nodular pile in layered soils[J]. Rock and Soil Mechanics， 2017， 38（1）： 109-116. （in Chinese）

[7]李鏡培， 陳浩華， 李林， 等. 楔形單樁與群樁非線性荷載沉降曲線計(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）

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

[9]高盟， 高廣運(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）

[10]潘軍， 楊小平， 劉庭金. 考慮樁身非線性壓縮的單樁沉降解析算法[J]. 地下空間與工程學(xué)報(bào)， 2017， 13（6）： 1579-1584.

PAN J， YANG X P， LIU T J. Analytical solution for single pile settlement considering non-linear compression of pile shaft[J]. Chinese Journal of Underground Space and Engineering， 2017， 13（6）： 1579-1584. （in Chinese）

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

[12]張瑜， 吳江斌， 蔡永昌. 考慮樁徑和樁側(cè)摩阻力分布的改進(jìn)Geddes樁基沉降分析方法[J]. 巖土力學(xué)， 2015， 36（Sup2）： 406-412.

ZHANG Y， WU J B， CAI Y C. A modified Geddes s method for estimating settlement of pile foundation considering influences of pile diameter and distribution pattern of shaft resistance[J]. Rock and Soil Mechanics， 2015， 36（Sup2）： 406-412. （in Chinese）

[13]趙明華， 馬思齊， 肖堯， 等. 基于厚壁圓筒模型的筋箍碎石樁沉降計(jì)算[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2019， 46（1）： 101-108.

ZHAO M H， MA S Q， XIAO Y， et al. Settlement calculation of geogrid-encased stone columns based on thick-walled cylinder model[J]. Journal of Hunan University （Natural Sciences）， 2019， 46（1）： 101-108. （in Chinese）

[14]建筑樁基技術(shù)規(guī)范： JGJ 94—2008[S]. 北京： 中國(guó)建筑工業(yè)出版社， 2008.

Technical code for building pile foundations： JGJ 94-2008[S]. Beijing： China Architecture & Building Press， 2008. （in Chinese）

[15]SHENG D C， EIGENBROD K D， WRIGGERS P. Finite element analysis of pile installation using large-slip frictional contact[J]. Computers and Geotechnics， 2005， 32（1）： 17-26.

[16]賈煜， 宋福貴， 王炳龍， 等. 基于改進(jìn)荷載傳遞法計(jì)算降水引起的基樁沉降[J]. 巖土力學(xué)， 2015， 36（1）： 68-74， 82.

JIA Y， SONG F G， WANG B L， et al. Modified load transfer method for calculation of foundation pile settlement due to dewatering[J]. Rock and Soil Mechanics， 2015， 36（1）： 68-74， 82. （in Chinese）

[17]林智勇， 戴自航. 由單樁載荷試驗(yàn)推算群樁沉降的相互作用系數(shù)法[J]. 巖土工程學(xué)報(bào)， 2016， 38（1）： 155-162.

LIN Z Y， DAI Z H. Interaction factor method for piles group settlement by static load tests of single pile[J]. Chinese Journal of Geotechnical Engineering， 2016， 38（1）： 155-162. （in Chinese）

[18]張乾青， 李術(shù)才， 李利平， 等. 考慮側(cè)阻軟化和端阻硬化的群樁沉降簡(jiǎn)化算法[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）

[19]劉杰， 張可能， 肖宏彬. 考慮樁側(cè)土軟化時(shí)單樁荷載沉降關(guān)系的解析算法[J]. 中國(guó)公路學(xué)報(bào)， 2003， 16（2）： 61-64.

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

[20]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.

[21]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.

[22]王奎華， 呂述暉， 吳文兵， 等. 層狀地基中基于虛土樁模型的單樁沉降計(jì)算方法[J]. 工程力學(xué)， 2013， 30（7）： 75-83， 99.

WANG K H， LV 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， 99. （in Chinese）

[23]王奎華， 羅永健， 吳文兵， 等. 層狀地基中考慮樁端應(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）

[24]辛冬冬， 張樂文， 宿傳璽. 基于虛土樁模型的層狀地基群樁沉降研究[J]. 巖土力學(xué)， 2017， 38（8）： 2368-2376， 2394.

XIN D D， ZHANG L W， SU C X. Settlement research of pile groups in layered soils based on virtual soil-pile model[J]. Rock and Soil Mechanics， 2017， 38（8）： 2368-2376， 2394. （in Chinese）

[25]王寧， 王奎華， 房凱. 考慮應(yīng)力擴(kuò)散時(shí)樁端土對(duì)樁體阻抗的影響[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào)， 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]POULOS H G. Modified calculation of pile-group settlement interaction[J]. Journal of Geotechnical Engineering， 1988， 114（6）： 697-706.

[27]洪毓康， 陳強(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）

[28]RANDOLPH M F， WORTH C P. Analysis of deformation of vertically loaded piles[J]. Journal of Soil Mechanics and Foundations Division， ASCE， 1978（12）： 56-61.

[29]董金榮， 林勝天， 戴一鳴. 大口徑鉆孔灌注樁荷載傳遞性狀[J]. 巖土工程學(xué)報(bào)， 1994， 16（6）： 123-131.

DONG J R， LIN S T， DAI Y M. The load tranfer behavior of large diameter cast-in-situ pile in crushed pebble stratum[J]. Chinese Journal of Geotechnical Engineering， 1994， 16（6）： 123-131. （in Chinese）

[30]羅永健. 考慮虛土樁擴(kuò)散角時(shí)基樁沉降計(jì)算理論及工程應(yīng)用[D]. 杭州： 浙江大學(xué)， 2013.

LUO Y J. Settlement calculation theory of single pile considering cone angle of fictitious pile and its application[D]. Hangzhou： Zhejiang University， 2013. （in Chinese）

（編輯 黃廷）