樁承式路堤土拱效應三維分析

摘 要：通過模型試驗和有限元分析，研究樁承式路堤中填土的破壞模式。結果表明，對于填土高度與樁凈間距之比H/s-a≤1的低路堤，填土中的破壞面為通過樁邊緣的豎直面；而對于H/s-a≥1.5的高路堤，破壞面為位于樁頂以上的泡狀曲面，可假設為對數螺旋曲面。基于模型試驗和有限元分析得到的破壞模式，建立了土拱效應的三維簡化分析方法，并與有限元和文獻收集得到的試驗數據進行了比較。結果表明，由于采用了合適的破壞面假設，該簡化方法的計算結果可靠合理。

關鍵詞：樁承式路堤；土拱效應；有限元；破壞模式

中圖分類號：TU 470.3 文獻標志碼：A 文章編號：16744764（2012）05007708

樁承式路堤是指在地基中設置剛性樁或半剛性樁來支承路堤填土荷載，其具有施工速度快、沉降小、穩定性高等特點，近年來在工程中得到了廣泛應用[12]。

在路堤荷載作用下，地基土的沉降大于樁的沉降，使得路堤填土中出現剪應力，荷載向樁頂轉移，從而減小了作用在地基土表面的壓力，這種現象稱為“土拱效應”。樁承式路堤的設計思路在于經濟合理的布置樁體，盡可能發揮路堤填土中的土拱效應，減少地基土承擔的荷載。另一方面，如果在樁頂布置土工格柵等加筋材料，加筋材料的選擇取決于其所承擔的荷載。因而合理評價土拱效應對樁承式路堤的設計至關重要[3]。

Terzaghi [4]根據活板門試驗研究了平面土拱效應，其假設填土中的破壞面為通過活動門邊緣的豎直面，作用在地基土表面的荷載等于滑動體的重力減去邊界上的摩阻力。Russell等[5]、Sloan等[6]將Terzaghi平面土拱效應理論推廣到三維情況。Naughton[7]認為填土中的破壞面為通過樁邊緣的對數螺旋曲線，該對數螺旋曲線的中心點為樁間土體的形心，作用在地基表面上的荷載為破壞面下方土體的重力，結果只適用于二維。Hewlett等[8]（以下簡稱HR方法）根據二維和三維的室內模型試驗，假設填土中存在支撐于樁頂的圓形土拱，以拱頂或拱腳土體單元的極限平衡條件為設計控制狀態。Zeaske等[9]根據模型槽試驗和數值分析，假設三維土拱為多個殼單元組成的系統，即所謂的“多拱模型”?！?D（〗 費 康，等：樁承式路堤土拱效應三維分析〖=〗

由于采用了不同的路堤填土破壞模式，各方法之間的計算結果差異較大，與實測數據也有所區別[1011]。筆者首先根據模型槽試驗和有限元數值分析對樁承式路堤填土中的三維破壞模式進行分析， 確定破壞面形狀；然后進行理論推導，得到三維土拱效應的簡化計算方法；最后與有限元計算結果和實測數據進行對比。1 破壞模式的模型槽試驗研究

1.1 試驗概況

1.1.1 試驗材料 試驗在一長100 cm、寬100 cm、高150 cm的鋼制模型槽中進行。為方便觀測，模型槽一側為鋼化玻璃。試驗用樁采用C20混凝土，截面為正方形，邊長a=10 cm，樁長L=40 cm，樁中心間距s=50 cm，為方便觀察樁頂土體變形，試驗為“半模試驗”，共使用了2根整樁和4根半樁，并將半樁與玻璃緊密貼合，具體布置方式見圖 1。

為了使得樁和地基土之間產生足夠大的差異變形，充分發揮土拱效應，試驗中地基土采用泡沫模擬， 厚度t=40 cm。根據室內側限壓縮試驗結果，在15%的應變范圍內，泡沫應力應變近似為線性關系，壓縮模量為28.2 kPa。

路堤填土采用標準石英砂模擬，砂土的比重ds為2.65，不均勻系數Cu為3.22，曲率系數Cc為117。填筑時的控制干密度ρd=1.5 g/cm3，相應密實度下由直剪試驗確定的摩擦角為27.3°。為減少模型槽側壁摩擦的影響，在砂土與玻璃壁之間涂抹凡士林，另外3邊鋪設塑料薄膜。由于樁的模量遠大于地基土，樁的變形很小，模型槽側壁摩擦對樁的性狀影響較小，因此試驗中未做額外的處理。

1.1.2 試驗方法及儀器布置 試驗中先將模擬軟土地基的泡沫就位，再將樁插入到預留的樁孔中，隨后進行砂土的填筑。填筑過程中以5 cm為1層，分層填筑。每層填土之間鋪設薄層彩砂，并在玻璃外側繪制直線，通過比較變形前后的照片，可以直觀的觀察路段填土在不同填土高度下的變形模式。

為了評價填土中的土拱效應，試驗中在2根樁的樁頂各布置了一個土壓力盒（T1和T2），從樁間土的中心點往上以100 mm為間距共布置了10個土壓力盒（T3到T12），儀器具體布置見圖 1。

1.2 試驗結果分析

1.2.1 樁頂周圍填土的變形模式 填土高度為20 cm（Hs-a=0.5）和80 cm（Hs-a=2）時的變形模式及確定的破壞面分別見圖 2（a）和（b），圖中實線代表變形前，虛線代表變形后。結果表明，當填土較低時，不同高程處的彩砂在樁邊緣都出現了明顯的轉折，樁間填土的位移大于樁頂土體?？紤]到破壞面兩側的沉降變形變化劇烈，可認為破壞面為通過樁邊緣的豎直或略傾斜的平面， 且貫穿了整個填土層，這與Terzaghi方法中的假設一致。當填土較高時，樁頂以上局部范圍之內的土體沉降小于樁間填土的沉降，填土中伴隨有剪切作用；超過一定高度之后，填土中的沉降近似均勻分布，即存在著等沉面。同樣可根據填土的變形模式確定破壞面，在樁邊緣破壞面與豎直線成一角度向上傾斜發展，當達到一定高度后（本例中約為0.8 a）破壞面向樁中心線彎曲，整個破壞面形狀類似于泡狀，局限在樁頂1.8 a以上范圍內，破壞面的水平范圍最大為1.6 a。

1.2.2 填土中的豎向應力分布 圖3給出了實測的豎向應力分布。圖中的深度是由填土表面向下起算的距離。結果表明，當填土高度較低時，豎向應力沿深度單一增加。當填土高度較高時，豎向應力分布曲線從上到下大致可以分為沿深度增加、減小、再增加3個階段。

豎向應力沿深度的分布實質上是破壞模式的體現，圖 4將填土高度1.0 m時的豎向應力分布理論值與實測值進行了比較。若認為破壞模式符合Terzaghi方法，在不考慮等沉面的情況下，計算得到的豎向應力應沿深度逐漸減小，且減小的速率越來越小，這與填土高度較小時的豎向應力分布特點相近。當填土高度較高時，考慮到樁頂某一高度以上的填土中沉降基本相等，該高度以上的填土中沒有剪應力，Terzaghi方法中將其作為超載考慮。這種情況下計算得到的豎向應力分布曲線大致可分為2段，在等沉面高度以上為自重應力分布，以下則隨深度逐漸減小。由此可見，無論是否考慮等沉面，Terzaghi方法的計算結果與試驗實測曲線的分布形態都有所區別。

根據HR方法，填土中應存在球形土穹，對正方形布置的方樁土穹的內半徑為s-a2，外半徑為s2。在土穹頂部外側以上的填土中豎向應力為自重應力線性分布；在土穹中，豎向應力減??；在土穹以下，豎向應力又開始增加。按照HR的土拱效應理論，出現第2段曲線中的豎向應力的減小是因為填土中出現了土穹，豎向應力減小段的始點對應于土穹的外頂點，終點對應于土穹的內頂點。從試驗結果來看，反算的土穹內頂點的位置比HR方法明顯的要低，HR方法也存在一定的誤差。

值得注意的是，如果在Terzaghi方法中考慮圖 2（b）所示的曲線滑動面，即考慮滑動體的寬度隨高程變化，得到的豎向應力分布形態能較好的模擬實測規律（圖 5），這表明將破壞面看做局限在樁頂的泡狀曲面是合理的。

另外注意到，填土高度1.0 m時實測地基頂面豎向應力為7.4 kPa，約為填土自重應力的0.49。此時樁頂應力實測值為178.8 kPa（T1）和145.9 kPa（T2），若將樁頂應力取為T1和T2的平均值162.3 kPa，則樁和土共同承擔的荷載為13.6 kN，與填土總荷載理論值15.0 kN相差約9%，驗證了實測數據的可靠性。

2.1 有限元分析模型

2.1.1 幾何模型及邊界條件 計算中取一根樁及其控制范圍進行分析。考慮到對稱性，所采用分析模型的平面示意圖如圖 6所示，圖中的OAB為1/4樁，ACDEB為樁間地基土的表面。

不同摩擦角下的計算結果也表現出相同的規律，以s=2.5 m，H=6.0 m為例（圖 8），不同摩擦角下的破壞面都顯示出相同的特性，即破壞面局限在樁頂上部填土內，且可近似看成對數螺旋面；摩擦角越大，破壞面范圍也越大。

3 三維土拱效應簡化分析方法

根據模型試驗及有限元分析結果，樁承式路堤中填土的破壞模式可分為低路堤和高路堤破壞模式2種，當路堤填土較小時（Hs-d≤1），填土中的破壞面為通過樁邊緣的豎直面，可仍采用Terzaghi土拱效應理論分析分析。當路堤填土較高時（Hs-d≥1.5），填土中的破壞面局限在樁頂局部范圍內，可近似看成對數螺旋曲面，本章對高路堤破壞模式進行推導，進而得到相應的土拱效應分析方法。

3.1 破壞面形狀

建立如圖 9所示的軸對稱坐標系，原點O為樁中心點，r軸水平向右，z軸豎直向上。OA為樁的半徑d2。假設破壞面ABCD是以O′為中心的對數螺旋曲面，BE為通過O′點的豎直線，對數螺旋曲線上任一點D的切線方向與水平方向的夾角為δ，D點距O′點的距離R為R1exp（ω-α1）tan φ，R1是O′C的長度，α1、α2、α3和ω意義如圖所示。O′C的長度為R2=R1exp（α2+α3）tan φ。

圖 10（a）給出了應力折減系數Sr與填土高度和樁凈間距之比n=Hs-d的變化關系，圖中的m=ds-d，反映了不同樁間距的影響。由圖可見，雖然簡化分析方法計算得到的應力折減系數略高于有限元計算結果，但整體變化規律是一致的。隨著Hs-d增加Sr減小，但降低的速率逐漸放緩，填土高度較高后（Hs-d>4），Sr基本保持不變。而對不同m值，m值越小，Sr越大，地基土需要承擔更多的荷載，因而需要合理的設置樁間距及樁帽尺寸。同時也注意到，當Hs-d<1.5時，m對Sr的影響不明顯，這是因為此時路堤填土中的破壞模式主要是低填土路堤模式，土拱效應的程度只取決于破壞面上發揮的摩擦力，Sr主要與Hs-d有關。

圖 10（b）給出了s=2.5 m，H=6 m時應力折減系數Sr與摩擦角φ之間的關系。結果表明，有限元方法和該簡化方法得到的規律是一致的。隨著摩擦角的增加，Sr減小的很明顯，當φ=20°時，簡化分析方法計算得到Sr=0.84；當φ=40°，Sr快速降低為0.25，更多的荷載通過土拱效應傳遞給樁。

取值，這2點仍然存在不確定因素。整體上來看，對Hs-d≥1.5的路堤，由于采用了更合適的破壞模式，筆者提出的簡化分析方法計算結果最好。如對試驗1，實測的應力折減系數為0.49，該方法計算值為055，在所有方法的計算結果中最為接近。若考慮樁頂荷載，試驗1中的樁頂應力實測值（T1和T2的平均值）為162.3 kPa，對應的極限荷載為1623 kN，與由該方法式（15）計算得到的1.792 kN相差約94%，吻合度也較好。這是因為應力折減系數Sr反映的是樁間土體分擔荷載的情況，Sr相同意味著土體承擔的荷載相同，也就意味著樁承擔的荷載也相同。

另外，也需要指出，該方法求解的是極限情況下的土拱效應。因而對試驗7～9，該方法的計算結果明顯偏小。這是因為試驗7～9是現場測試數據，地基土有一定的支撐作用，樁和樁間土的差異沉降未能使得土拱效應完全發生。而試驗1、試驗2～4中地基土采用泡沫模擬，試驗5～6地基土采用軟弱粘土，試驗10～11中采用了類似活動門試驗的做法，地基土表面可任意沉降。在這些試驗中樁和樁間土之間都產生了足夠的差異沉降，填土達到了極限情況。5 結 論

1）三維試驗和有限元分析結果表明，樁承式路堤的破壞模式取決于路堤填土的高度，填土高度小于樁凈間距時，破壞面從樁邊緣向上豎直延伸至填土表面；隨著填土高度的增加，破壞面逐漸演變為局限在樁頂的泡狀曲面。

2）根據試驗和有限元分析結果所假設的高路堤破壞模式，建立了土拱效應三維分析方法，并與有限元數值分析和試驗結果進行了對比，結果吻合較好，驗證了破壞模式和分析方法的可靠性。

3）所提出的三維簡化分析方法適用于求解土拱效應完全發揮時的應力折減系數，未達到極限狀態時樁、土之間的荷載分擔需進一步研究。

參考文獻：

[1]費康， 劉漢龍. 樁承式加筋路堤設計理論研究進展[J]. 水利水電科技進展， 2008， 28（4）： 8994.

FEI Kang， LIU Hanlong. Advances in research on design theory of geosynthetic reinforced pile supported embankment[J]. Advances in Science and Technology of Water Resources， 2008， 28（4）： 8994.

[2]費康， 劉漢龍， 高玉峰. 路堤下現澆薄壁管樁復合地基工作特性分析[J]. 巖土力學， 2004， 25（9）： 13901396.

FEI Kang， LIU Hanlong， GAO Yufeng. Performance of castinsitu concrete thinwall pipe pile composite foundation under road embankment[J]. Rock and Soil Mechanics， 2004， 25（9）： 13901396.

[3]Potts V J， Zdravkovic L. Finiteelement study of arching behavior in reinforced soils[J]. Ground Improvement， 2010， 163（4）： 217229.

[4]Terzaghi K. Theoretical soil mechanics[M]. New York： Wiley， 1943.

[5]Russell D， Pierpoint N. An assessment of design methods for piled embankments[J]. Ground Engineering， 1997，30（11）：3944.

[6]Sloan J A， Filz G M， Collin J G. A generalized formulation of the Adapted Terzaghi Method of arching in columnsupported embankments[C]//Proc. of GeoFrontiers， GSP 211 Advanced in Geotechnical Engineering， Dallas， Texas， USA， 2011：798805.

[7]Naughton P J. The significance of critical height in the design of piled embankments[C]//Proc. GeoDenver， GSP 172 Soil Improvement， Denver Colorado， USA， 2007：1323.

[8]Hewlett J， Randolph M F. Analysis of piled embankments[J]. Ground Engineering， 1988， 21（3）： 1218.

[9]Zaeske D， Kempfert H G. Berechnung und wirkungsweise von unbewehrten und bewehrten mineralischen Tragschichten über punktund linienfrmigen Traggliedern[J]. Bauingenieur， 2002 （77）： 8086.

Zaeske D， Kempfert H G.Calculation and mode of unreinforced and reinforced base layers on point and line shaped support members[J]. Civil Engineering， 2002 （77）： 8086.

[10]費康， 劉漢龍. 樁承式加筋路堤的現場試驗及數值分析[J]. 巖土力學， 2009， 30（4）： 10041012.

FEI kang， LIU Hanlong. Field test study and numerical analysis of a geogridreinforced and pilesupported embankment[J]. Rock and Soil Mechanics， 2009， 30（4）： 10041012.

[11]Ellis E， Aslam R. Arching in piled embankments： comparison of centrifuge tests and predictive methods[J]. Ground Engineering， 2009， 42（6）：3438.

[12]費康，王軍軍，陳毅. 樁承式路堤土拱效應的試驗和數值研究[J]. 巖土力學， 2011， 32（7）： 19751983.

FEI kang， WANG Junjun，CHEN Yi. Experimental and numerical studies of soil arching in piled embankment[J]. Rock and Soil Mechanics， 2011， 32（7）： 19751983.

[13]陳云敏，賈寧，陳仁朋. 樁承式路堤土拱效應分析[J]. 中國公路學報，2004，17（4）：16.

CHEN Yunmin， JIA Ning，CHEN Renpeng. Soil arch analysis of pilesupported embankments[J]. China Journal of Highway and Transport，2004，17（4）：16.

[14]Demerdash M A. An experimental study of piled embankments incorporating geosynthetic basal reinforcement[D]. Newcastle upon Tyne：Newcastle University，1996.

（編輯 胡英奎）doi：10.3969/j.issn.16744764.2012.05.013