劉飛成,張建經(jīng),閆世杰

(西南交通大學(xué)土木工程學(xué)院,四川 成都 610031)

隨著中國東部沿海地區(qū)的快速發(fā)展,高速鐵路修建的需求也越來越大,而這一地區(qū)廣泛分布著承載能力差、壓縮性高的深厚軟土,易引起高速鐵路路堤的沉降超限,這也是各國的巖土工程師面臨的一大難題.樁基礎(chǔ)作為改良軟土地基的一種有效手段[1-2]已經(jīng)被廣泛運用在軟土地基的路堤修建中,并取得了良好的效果[3-4].

目前關(guān)于樁承式路堤的研究主要集中在承載能力分析[5-6]以及力的傳遞機制[7-8]上,對樁承式路堤的沉降研究相對較少.對一些軟土地基路堤的設(shè)計而言,路堤沉降控制相比于地基承載能力更為重要[9-10].Abusharar等[5]認(rèn)為目前尚沒有被廣泛認(rèn)可的針對樁承式路堤的沉降計算方法,需要進(jìn)行深入研究.路堤沉降與下部復(fù)合地基的沉降有著密切關(guān)系,準(zhǔn)確計算下部復(fù)合地基沉降對路堤沉降的計算具有重要的意義.

目前復(fù)合地基的沉降計算方法思想主要是把沉降分為加固區(qū)沉降和非加固區(qū)沉降兩部分來考慮.主要方法有規(guī)范方法、樁身壓縮量法、雙層應(yīng)力法和附加應(yīng)力法.其中：規(guī)范方法[11]將加固區(qū)看成一個整體并用復(fù)合模量表示其整體的壓縮能力,通過Boussinesq解計算應(yīng)力,采用分層總和法計算加固區(qū)和非加固區(qū)沉降；樁身壓縮量法[12]將復(fù)合地基的沉降看成是由樁身壓縮量和樁端刺入變形組成,該方法在確定樁的荷載分擔(dān)值時主要依賴于經(jīng)驗；雙層應(yīng)力法[13]認(rèn)為加固區(qū)土體和下臥層土體中的應(yīng)力可分別用Boussinesq解獲得,進(jìn)而分別求解加固區(qū)和下臥層的沉降,其樁間土頂面和下臥層頂面應(yīng)力水平的不確定性較大；附加應(yīng)力法[14-15]是基于Mindlin解推導(dǎo)土體中的附加應(yīng)力進(jìn)而求解土體的壓縮變形量,但仍是基于均質(zhì)土體進(jìn)行考慮的,這與實際情況明顯不同,誤差偏大.

基于此,本文提出了一種樁承式路堤沉降的計算方法,可分別計算樁頂沉降和樁間土沉降.將樁頂沉降分為樁端沉降、樁側(cè)摩阻力引起的沉降及樁身壓縮,采用球孔擴(kuò)張法計算樁端荷載,樁側(cè)摩阻力引起的沉降可通過Vesic解求得.將樁間土沉降分為樁間土表面荷載、樁側(cè)摩阻力和樁端力引起的沉降.樁端沉降和樁間土沉降均認(rèn)為是由地基土的附加應(yīng)力引起,路堤荷載作用下地基土任意點的附加應(yīng)力通過考慮多土層的非均質(zhì)地基土情況利用修正Mindlin解求得.將本文提出的計算方法所得結(jié)果與現(xiàn)場實例和數(shù)值分析結(jié)果進(jìn)行比較,以驗證該方法的合理性.

1　簡化計算方法的建立

為方便接下來的公式推導(dǎo),作出如下假設(shè):

(1) 路堤填土是均質(zhì)、各向同性的材料;(2) 每一層地基軟土是均質(zhì)、各向同性且僅發(fā)生豎向變形;(3) 樁身剛度足夠大不會產(chǎn)生破壞;(4) 不考慮土工格柵的影響,路堤填土的高度足夠形成完整土拱.

1.1　土拱效應(yīng)

1.1.1正方形布樁

Hewelett等[16]根據(jù)模型試驗觀測的結(jié)果認(rèn)為,在正方形布樁情況下,樁間將形成一個半球形土拱和4個平面土拱,如圖1所示.其中半球形土拱以相鄰對角樁的樁間距為直徑,相鄰?fù)凉按嬖诓糠种丿B.土拱的極限狀態(tài)僅可能出現(xiàn)在兩個位置:土拱拱頂和拱腳.取樁的荷載分?jǐn)偙容^小的情況作為實際工作狀態(tài).

(1) 拱頂達(dá)到極限狀態(tài)

根據(jù)拱頂單元的豎向力平衡,土拱單元的法向應(yīng)力為[16]

(1)

(2)

式中:r為路堤填土的重度；H為路堤填土高度;Kp為朗肯被動土壓力系數(shù);s為樁間距；γ為路堤填土的重度;φs為路堤填土的內(nèi)摩擦角,對于黏性土而言,其等效內(nèi)摩擦角φeq,s=arctan(tanφs+c/γH)，其中，c為樁端土的黏聚力.

因此,樁間土中間處所受到的荷載為

(3)

式中:a為樁帽寬度.

式(3)求得的實際上是樁間土壓力的最大值,而樁間土應(yīng)力并非均勻分布,尤其對軟土地基而言.考慮地基表面荷載分布不均勻性[17],引入均勻系數(shù)μ，則樁間土表面等效均布荷載為

σbs=μσs.

(4)

樁的荷載承擔(dān)比為

(5)

圖1　土拱模型Fig.1　Arching model

(2) 拱腳達(dá)到極限狀態(tài)

當(dāng)拱腳單元處于極限狀態(tài)時,土拱內(nèi)邊界的法向應(yīng)力認(rèn)為等于Kpσbs.根據(jù)拱腳土體單元的法向力平衡可得到土拱的法向應(yīng)力為[16]

(6)

單個樁帽受到的合力由附近4個平面應(yīng)變土拱產(chǎn)生,每個平面土拱均占據(jù)1/4樁帽面積.樁帽受到的合力可通過對樁帽上的不同半徑處應(yīng)力進(jìn)行積分得到,即

(7)

根據(jù)豎向荷載平衡,可得到

Pt+σbs(s2-a2)=γHs2.

(8)

由(8)可得到σbs,進(jìn)而可求得樁的荷載分?jǐn)偙?

1.1.2三角形布樁

圖2　三角形布樁的土拱模型Fig.2　Arch supported by piles arranged in triangle pattern

(1) 拱頂達(dá)到極限狀態(tài)

(9)

因此樁間土中間處所受到的荷載為

(10)

由式(4)、(5)和(10)可求得樁的荷載分?jǐn)偙?

(2) 拱腳達(dá)到極限狀態(tài)

單個樁帽受到的合力由附近6個平面應(yīng)變土拱產(chǎn)生,每個平面土拱均占據(jù)1/6樁帽面積.因此樁帽受到的合力也可通過對樁帽上的不同半徑處應(yīng)力進(jìn)行積分得到,即

(11)

同樣根據(jù)式(8)可求得樁間土表面等效均布荷載σbs,進(jìn)而可求得樁的荷載分?jǐn)偙?

1.2　樁頂沉降計算

單樁樁頂沉降Spt可認(rèn)為由3部分組成:樁端力引起的樁端沉降Spb,樁身壓縮量Spc和樁側(cè)阻力引起的沉降Sps,即

Spt=Spd+Spc+Sps.

(12)

1.2.1樁端Vesic球孔擴(kuò)張理論

Yasufuku和Hydea[18]通過試驗觀察到樁端的向下刺入類似于球孔擴(kuò)張過程.因此,引入Vesic球孔擴(kuò)張理論[19]進(jìn)行樁端阻力計算,則樁端應(yīng)力可表示為

(13)

式中:σ0為樁端平均應(yīng)力;φ為樁端土的內(nèi)摩擦角.

群樁基礎(chǔ)的樁間距一般都小于6倍樁徑,考慮到相鄰樁之間的橫向相互擠壓,對樁端下刺具有抑制作用,樁端阻力可近似計算為[20]

(14)

1.2.2樁端沉降計算

對于均質(zhì)彈性半空間體內(nèi)深度h處作用集中力P時,深度z處任意點的豎向應(yīng)力可表示為[21]

(15)

式中:R1、R2為計算半徑,R1=(D2+(z-h)2)1/2,R2=(D2+(z+h)2)1/2;ν為土體的泊松比;D為所求點與集中力作用點的水平距離.

將樁端荷載看成圓形均布荷載,采用圓柱坐標(biāo)系,力的作用點坐標(biāo)為(r0,θ0,L),任意一點的坐標(biāo)為(r1,θ1,z),兩點的水平距離為

(16)

在圓形均布荷載選取一個面積元荷載dP=σptrdφdr,其在地基內(nèi)任意點引起的附加應(yīng)力為

(17)

樁端荷載在地基內(nèi)任意點的附加應(yīng)力可通過在圓形面積內(nèi)積分求得,即

(18)

對于群樁基礎(chǔ)而言,顯然與半無限彈性空間體這一前提條件不符.因此本文認(rèn)為在單根樁的處理范圍內(nèi)滿足Mindlin解的前提條件.對于正方形布置的群樁,僅考慮單樁處理范圍內(nèi)的樁和樁間土荷載的影響.

根據(jù)樁間土表面荷載是均勻分布的假定,樁間土表面荷載在地基內(nèi)任意點產(chǎn)生的附加應(yīng)力為

(19)

樁端下部土層內(nèi)任意點的附加豎向應(yīng)力為

σz=σz1+σz2.

(20)

對于成層地基而言,可將每個土層看成半無限空間體,對每個土層分別采用Mindlin公式求解,令第i-1層土層下邊界的附加應(yīng)力為σ(i-1)b,第i層土層上邊界的附加應(yīng)力為σit,根據(jù)邊界附加應(yīng)力連續(xù)條件,可得到

σ(i-1)b=σiz.

(21)

根據(jù)式(19)、(20)和(21)即可求得成層土地基任意點的附加應(yīng)力，隨后利用分層總和法可得到樁端荷載和樁間土表面荷載在樁端持力層產(chǎn)生的附加應(yīng)力引起的持力層壓縮,即為樁端的沉降量Spb.

1.2.3樁身壓縮量

一般來說,樁身壓縮量的計算是先將樁分成若干段,假定整根樁的樁身橫截面積不變,樁身彈性模量不變,分別計算每一段樁身的壓縮量,然后將各段的樁身壓縮量相加即可得到整根樁的樁身壓縮量,但這需要對樁身的軸力分布有明確認(rèn)識,通常情況下較難實現(xiàn).張乾青[22]通過現(xiàn)場試樁靜載試驗發(fā)現(xiàn),分段計算得到的全樁長樁身壓縮量略大于直接利用樁頂荷載和樁端力計算得到全樁長樁身壓縮量,極限荷載下分段計算的樁身壓縮量接近于樁身實際壓縮量,兩者的比值為0.86～0.93.本文為簡化計算,直接利用樁頂荷載和樁端力計算全樁長樁身壓縮量，即

(22)

式中:Ap為樁身橫截面積;Ep為樁身彈性模量;d為樁的直徑;χ為修正系數(shù),本文建議χ=1.10.

1.2.4樁側(cè)阻力引起的沉降

將樁周土體假設(shè)為彈性各向同性的彈簧單元,其性質(zhì)與土體的變形模量和泊松比有關(guān),則樁側(cè)阻力引起的沉降Sps為[23]:

(23)

(24)

式中:L為樁長;Es為地基土的壓縮模量.

若地基土含有多層土層,則地基土的復(fù)合壓縮模量可通過式(25)求解[24].

(25)

式中:h1(Es1)、h2(Es2)、…、hn(Esn)為成層地基的第1層、第2層、…、第n層的土層厚度(彈性模量).

1.3　樁間土沉降計算

樁間土沉降可認(rèn)為由樁間土表面荷載、樁側(cè)摩阻力和樁端力產(chǎn)生的附加應(yīng)力引起,求解三者在地基中引起的附加應(yīng)力是計算樁間土沉降的前提.樁端力和樁間土表面荷載產(chǎn)生的附加應(yīng)力可分別采用式(14)、(15)進(jìn)行求解.

樁側(cè)摩阻力的分布對計算樁側(cè)摩阻力在地基內(nèi)任意點引起的附加應(yīng)力至關(guān)重要.當(dāng)路堤填土荷載較小時,樁土相對位移較小,樁側(cè)摩阻力處于彈性階段,此時樁側(cè)摩阻力和樁土相對位移呈線性關(guān)系;當(dāng)路堤填土荷載較大時,樁側(cè)摩阻力與樁土相對位移呈非線性關(guān)系.在本文分析中,為簡化計算,假定樁側(cè)摩阻力沿深度呈三角形分布.

樁側(cè)某處的摩阻力在地基內(nèi)任意點引起的附加應(yīng)力為

(26)

樁側(cè)總的摩阻力在地基內(nèi)任意點引起的附加應(yīng)力為

(27)

對于成層土而言,可采用式(20)、(21)和(27)求解.同理,根據(jù)樁間土表面荷載、樁側(cè)摩阻力和樁端力產(chǎn)生的附加應(yīng)力,可通過分層總和求得樁間土沉降.

2　計算方法的驗證

為驗證本文計算方法的合理性,選取兩個樁承式路堤算例進(jìn)行分析,將本文方法的計算結(jié)果與現(xiàn)場實測和數(shù)值分析結(jié)果進(jìn)行對比.

2.1　算例1

現(xiàn)場實測數(shù)據(jù)采用Briancon等[25]對樁承式路堤沉降的監(jiān)測數(shù)據(jù).測試地點位于巴黎東北部20 km 處的一處鐵路路基,軟土層主要由低塑性的粘土和砂質(zhì)粘土,具有中度壓縮性,軟土層下部為致密礫石層.僅采用樁基礎(chǔ)加固路段的路堤填高均為5 m,樁的直徑均為0.30 m,間距為2 m的正方形布樁,各樁樁端均嵌入致密礫石層0.3 m.土體物理力學(xué)參數(shù)參考文獻(xiàn)[25-26],列于表1、2.路堤高度足夠形成完整土拱,計算樁長取8.8 m.

由于參數(shù)未給出壓縮模量,因此根據(jù)上述參數(shù),進(jìn)行壓縮模量的求解.一般來說,對于土的壓縮曲線,采用半對數(shù)坐標(biāo)表示成e-lgp曲線,其分為兩段:曲線段和斜直線段.斜線方程可表示為

(28)

(29)

式中:Cc為e-lgp曲線的直線段壓縮指數(shù);pc為壓縮曲線的轉(zhuǎn)折點對應(yīng)的壓力;e為壓縮曲線的轉(zhuǎn)折點對應(yīng)的孔隙比;p為作用在土體上的壓力;e0為壓力為p時對應(yīng)的孔隙比；av為土的壓縮系數(shù).

令mv為土體的體積壓縮系數(shù),體積壓縮系數(shù)與壓縮模量存在倒數(shù)關(guān)系.根據(jù)體應(yīng)變與垂直壓力成正比,土體體積壓縮系數(shù)與壓縮指數(shù)關(guān)系如式(30).

(30)

將式(28)和(29)代入到式(30)中,可得

(31)

因此可求得下部致密礫石的壓縮模量為6.38 MPa,軟土的等效壓縮模量為1.66 MPa.根據(jù)土拱效應(yīng)計算得到兩種狀態(tài)下樁的荷載分?jǐn)偙确謩e為0.736(拱頂極限狀態(tài))和0.373(拱腳極限狀態(tài)).取拱腳極限狀態(tài)進(jìn)行分析,Pt=137.9 kN,σbs=60.2 kPa.根據(jù)式(14)計算得到σpb=753.7 kPa,樁端沉降為7.21 mm,樁側(cè)摩阻力引起的沉降為0.37 mm,樁身壓縮量為0.04 mm.將計算結(jié)果與現(xiàn)場實測數(shù)據(jù)以及文獻(xiàn)[26]中的數(shù)值分析數(shù)據(jù)列于表3.由表3可知,本文提出的計算方法計算得到的樁頂沉降與現(xiàn)場實測相當(dāng)接近.

表1　土體物理力學(xué)參數(shù)Tab.1　Physical and mechnical parameters of soil

表2　路堤填土、致密礫石和樁的參數(shù)Tab.2　Parametes of soil and pile

表3不同計算方法得到的樁和土沉降
Tab.3Calculated and observed settlement of pile and soilmm

不同方法本文理論計算現(xiàn)場實測結(jié)果數(shù)值分析結(jié)果樁頂7．689．9樁間土92．410567．0

2.2　算例2

Chen等[8]利用數(shù)值分析方法計算了兩種工況下樁承式路堤.工況1:樁長20 m,地基土為25 m厚的軟土層,軟土層下部為硬土層;工況2:樁長為27 m,軟土層為25 m,軟土層下部為硬土層.將計算結(jié)果與其提出的計算模型進(jìn)行對比,計算模型的幾何參數(shù)以及材料物理力學(xué)參數(shù)如表4所示.理論計算結(jié)果和數(shù)值分析結(jié)果如表5所示.

表4　模型的幾何參數(shù)和材料物理力學(xué)參數(shù)Tab.4　Material properties used in the numerical modeling and the proposed analytical solution

表5不同計算方法得到的沉降
Tab.5The settlement by theoretical methods and numerical anaysismm

位置工況本文理論數(shù)值分析Chen模型樁頂1131．0148．0141．0216．312．021．0樁間土1152．0168．0165．0237．432．043．0

由表5可知：對工況1,本文方法與Chen的數(shù)值分析和計算模型結(jié)果吻合良好,說明本文方法和Chen計算模型在土層為均質(zhì)條件下可較準(zhǔn)確的預(yù)測樁頂沉降；對工況2,本文計算方法結(jié)果與數(shù)值分析結(jié)果仍吻合較好,說明本文方法能較準(zhǔn)確地預(yù)估路堤荷載下懸浮樁和嵌固樁復(fù)合地基的沉降.

3　結(jié)束語

本文提出了一種軟土地基樁承式路堤沉降的簡化計算方法.考慮土拱效應(yīng)對路堤荷載分配的影響,求得樁頂荷載和樁間土表面荷載,在此基礎(chǔ)上將樁間土沉降分為樁間土表面荷載、樁側(cè)摩阻力和樁端力引起的沉降,樁頂沉降分為樁端沉降、樁側(cè)摩阻力引起的沉降及樁身壓縮.樁間土沉降和樁端沉降均認(rèn)為是由地基土中附加應(yīng)力引起的,且成層地基土任意點的附加應(yīng)力可通過改進(jìn)的Mindlin解求得.通過將計算與現(xiàn)場實測結(jié)果進(jìn)行對比發(fā)現(xiàn),計算與實測結(jié)果吻合良好,說明本文方法的可靠性.

參考文獻(xiàn):

[1]MAGAN J P.Methods to reduce the settlement of embankments on soft clay:a review[J].Vertical & Horizontal Deformations of Foundatio Embankments,2015， 1(40):77-91.

[2]LIU H L,NG C W W,FEI K.Performance of a geogrid-reinforced and pile-supported highway embankment over soft clay case study[J].Journal of Geotechnical and Geoenvironmental Engineering,2007,133(12):1483-1493.

[3]陳云敏,賈寧,陳仁朋.樁承式路堤土拱效應(yīng)分析[J].中國公路學(xué)報,2004,17(4):1-6.

CHEN Yunmin,JIA Ning,CHEN Renpeng.Soil arch analysis of pile-supported embankments[J].Chinese Journal of Highway and Transport,2004,17(4):1-6.

[4]LU W H,MIAO L C.A simplified 2-D evaluation method of the arcing effect for geosynthetic-reinforced and pile-supported embankments[J].Computers and Geotechnics,2015,65(65):97-103.

[5]ABUSHARAR S W,ZHENG J J,CHEN B G,et al.A simplified method for analysis of a piled embankment reinforced with geosynthetics[J].Geotextiles and Geomembranes,2009,27(1):39-52.

[6]EEKELEN S J M V,BEZUIJEN A,TOL A F V.Validation of analytical models for the design of basal reinforced piled embankments[J].Geotextiles and Geomembranes,2015,43(1):56-81.

[7]ZHANG C L,JIANG G L,LIU X F,et al.Arching in geogrid-reinforced pile-supported embankments over silty clay of medium compressibility:field data and analytical solution[J].Computers and Geotechnics,2016,77:11-25.

[8]CHEN R P,CHEN Y M,HAN J,et al.A theoretical solution for pile-supported embankments on soft soils under one-dimensional compression[J].Canadian Geotechnical Engineering,2008,45(6):611-623.

[9]HAN J,OZTOPRAK S,PARSON R L,et al.Numerical analysis of foundation columns to support widening of embankments[J].Computers and Geotechnics,2007,34(6):435-448.

[10]SHAHIRA H,PAK A.Estimating liquefaction-induced settlement of shallow foundation by numerical approach[J].Computers and Geotechnics,2010,37(3):267-279.

[11]中國建筑科學(xué)研究院.JGJ79—2012 中華人民共和國行業(yè)標(biāo)準(zhǔn)建筑地基處理技術(shù)規(guī)范[S].北京:中國建筑工業(yè)出版社，2012.

[12]龔曉南.復(fù)合地基理論及工程應(yīng)用[M].北京:中國建筑工業(yè)出版社,2002:150-151.

[13]池躍君,宋二祥,陳肇元.剛性樁復(fù)合地基沉降計算方法的探討及應(yīng)用[J].土木工程學(xué)報,2003,36(11):19-23.

CHI Yuejun,SONG Erxiang,CHEN Zhaoyuan.An simplified method for settlement of rigid-pile composite foundation[J].China Civil Engineering Journal,2003,36(11):19-23.

[14]閆寶杰,宋修廣,岳越群.粉噴樁復(fù)合地基附加應(yīng)力的Boussinesq-Mindlin聯(lián)合解法[J].巖土力學(xué),2007,28(6):1255-1258.

YAN Baojie,SONG Xiuguang,YUE Yuequn.Boussinesq-Mindlin united solution method for DJMP composite foundation additional stress[J].Rock and Soil Mechanics,2007,28(6):1255-1258.

[15]陳永輝,王新泉,劉漢龍,等.Y型樁樁側(cè)摩阻力產(chǎn)生附加應(yīng)力的分析計算[J].巖土力學(xué),2008，29(11):2905-2911.

CHEN Yonghui,WANG Xinquan,LIU Hanlong,et al.Analysis and calculation of additional stress due to skin friction of Y-shaped vibro-pile[J].Rock and Soil Mechanics,2008， 29(11):2905-2911.

[16]HELWETT W J,RANDOLPH M F.Analysis of piled embankment[J].Ground Engineering,1988,21(3):12-18.

[17]LOW B K,TANG S K,CHOA V.Arching in piled embankments[J].Journal of Geotechnical Engineering,1994,120(11):1917-1938.

[18]YASUFUKU N， HYDE A F.Pile end-bearing capacity in crushable sands[J].Geotechnique,1995,45(4):663-676.

[19]VESIC A S.Bearing capacity of deep foundations in sand[J].Highway Research Record,1963,10(39):112-153.

[20]劉俊飛.鐵路CFG復(fù)合地基沉降控制機理與計算方法研究[D].成都:西南交通大學(xué),2011.

[21]MINDLIN R D.Force at a point in the interior of a semiinfinite solid[J].Physics,1936,7(5):195-202.

[22]張乾青.軟土地基樁基受力性狀和沉降特性試驗與理論研究[D].杭州:浙江大學(xué),2012.

[23]VESIC A S.Design of pile foundation[M].Washington:National Cooperative Highway Research Program Synthesis of Practice,1977:31-32.

[24]ZHUANG Y,WANG K Y,LIU H L.A simplified model to analyze the reinforced pile embankments[J].Geotextiles and Geomembranes,2014,42(2):154-165.

[25]BRIANCON L,SIMON B.Performance of pile-supported embankment over soft soil full-scale experiment[J].Journal of Geotechnical and Geoenvironmental Engineering,2011,138(4):551-561.

[26]ROWE R K,LIU K W.Three-dimensional finite element modelling of a full-sacle geosynthetic-reinforced,piled-supported embankment[J].Canadian Geotechnical Engineering,2015,52(12):2041-2054.