李龍起，趙皓璆

(成都理工大學(xué) 環(huán)境與土木工程學(xué)院，四川 成都 610059)

0 引 言

斜樁基礎(chǔ)是樁基中部分或全部基樁傾斜的一類基礎(chǔ)形式，斜樁基礎(chǔ)的樁基在水平和豎直承載力方面均表現(xiàn)出巨大的優(yōu)越性。近年來，對豎直樁的承載性能研究已經(jīng)相對成熟，然而，由于承臺-樁-土之間相互作用產(chǎn)生的群樁效應(yīng)使得傾斜樁基與豎直樁基在豎直承載力方面表現(xiàn)出較大的差異性。目前在工程界，有關(guān)傾斜群樁的設(shè)計多借鑒豎直群樁。由于對傾斜群樁的工作特性認(rèn)識尚不足，也就無法對工程進(jìn)行相應(yīng)的優(yōu)化，出現(xiàn)設(shè)計過于保守或資源浪費的狀況。因此，有必要深入研究傾斜群樁在實際工程中的承載力學(xué)性能。

近年來，國內(nèi)外眾多學(xué)者對斜樁在豎直荷載作用下的受力變形特性進(jìn)行了一系列的試驗研究，并取得了一定的成果。李龍起等[1-2]通過室內(nèi)靜載模型試驗發(fā)現(xiàn)傾斜樁基中相鄰兩基樁之間的間距隨土體埋深的增大而增大，土體的遮攔作用減小，繞樁滑動作用變強，從而造成中樁與樁間土體承載能力不同程度的發(fā)揮；曹衛(wèi)平等[3-4]通過分析不同樁身傾角及長徑比，研究了豎直荷載下斜樁的承載性狀，并與豎直樁的承載性狀進(jìn)行了比較；鄧友生等[5]研究了超大群樁基礎(chǔ)的沉降特性，結(jié)果表明樁身的彈性壓縮變形和自重引起的彈性變形對計算基礎(chǔ)總沉降量有影響；劉金礪等[6-7]通過大量群樁試驗研究了群樁效應(yīng)與樁基沉降之間的關(guān)系；呂凡任等[8]對承受豎直荷載的斜樁和承受傾斜荷載的豎直樁分別進(jìn)行了研究，得出結(jié)論：傾角很小的斜樁和豎直樁在承受豎直荷載時兩者承載力沒有明顯的差別，而豎直樁承受小于10°的傾斜荷載時，水平位移的變化比豎直位移的變化明顯的多；任瑞虹[9]采用有限元研究了樁-土相互作用對斜樁承載性能的影響；鄭剛等[10]利用有限元軟件ABAQUS模擬豎直荷載下不同傾角、不同土質(zhì)斜樁的承載力；王云崗等[11]研究了斜樁不同的軸向剛度與側(cè)向剛度比值條件下，側(cè)向與軸向荷載的分擔(dān)情況；王曙光等[12]通過模型試驗研究了螺桿樁的豎直受壓承載機理，給出了計算螺桿灌注單樁的極限承載力公式；常遠(yuǎn)[13]采用有限元軟件研究了豎直荷載下單樁樁土荷載的分擔(dān)情況，提出了樁土荷載分擔(dān)與樁土面積比的關(guān)系式。

上述關(guān)于豎直荷載作用下斜樁承載性能的研究大部分都集中在承載力上，而對斜樁在不同影響因素下內(nèi)力及樁身變形性狀的研究較少。筆者采用有限元軟件ABAQUS對豎直分級加載下的斜樁進(jìn)行數(shù)值模擬，研究了不同傾角、不同樁間距條件下，斜樁的彎矩、樁身位移、樁側(cè)摩阻力的變化。研究成果可為類似的樁基工程在前期分析及設(shè)計施工時提供參考。

1 有限元模型及參數(shù)確定

ABAQUS是一款功能強大的工程模擬非線性元軟件，應(yīng)用領(lǐng)域十分廣泛。ABAQUS在模擬樁-土相互作用時，首先需預(yù)先平衡地應(yīng)力，使土體的初始位移沉降盡可能接近0；然后再施加荷載，研究樁身的承載性能及內(nèi)力的變化。

1.1 樁土材料本構(gòu)模型

本構(gòu)模型主要指材料的應(yīng)力-應(yīng)變關(guān)系。數(shù)值模擬中樁基采用彈性模型，土體采用摩爾-庫倫(Mohr-Coulomb)彈塑性模型。

1.2 有限元模型尺寸

筆者基于文獻(xiàn)[2]的室內(nèi)物理試驗?zāi)Ｐ?，在有限元中按等比尺建模，開展了樁基承載特性研究?；鶚都俺信_尺寸如下：基樁長1 180 mm，直徑36 mm，承臺長×寬×高=420 mm×280 mm×10 mm。試驗各土層剖面及基樁平面布置如圖1[2]。數(shù)值模擬中試驗用土的物理力學(xué)指標(biāo)及承臺、樁基的物理力學(xué)指標(biāo)見表1。

表1 試驗用土及承臺、樁基的物理力學(xué)指標(biāo)Table 1 Physical and mechanical indexes of test soil, pile caps and pile foundations

圖1 試驗土層剖面及基樁平面布置Fig. 1 Sectional view of test soil layer and layout of foundation piles

1.3 樁-土接觸面及參數(shù)

樁側(cè)與樁間土的接觸通過面面接觸來實現(xiàn)。由于樁體與土體的彈性模量相差較大，兩者的接觸面會產(chǎn)生相對位移，為保證模擬結(jié)果的精確性，以樁體為主表面，土體為從表面，接觸類型為有限滑移。樁側(cè)與樁間土之間的切向摩擦采用庫倫摩擦，法向采用硬接觸。通過反復(fù)調(diào)整樁土接觸面之間的摩擦參數(shù)，最終得到：當(dāng)摩擦系數(shù)μ=0.4時，有限元計算結(jié)果與文獻(xiàn)[2]的結(jié)果最為接近。

有限元模型采用三維模型，為防止土體邊界產(chǎn)生位移，對土體兩側(cè)水平方向及土體底部豎直和水平方向設(shè)置了邊界條件；數(shù)值分析中，為保證計算結(jié)果的精確性，樁側(cè)與樁間土接觸區(qū)域設(shè)置較密的網(wǎng)格，遠(yuǎn)離接觸面區(qū)域設(shè)置較稀疏的網(wǎng)格；傾斜樁基、地基土采用三維實體10節(jié)點2次四面體單元C3D10。有限元模型網(wǎng)格劃分如圖2。

圖2 有限元模型網(wǎng)格劃分Fig. 2 FEM meshing diagram

1.4 有限元模擬工況

為研究豎直荷載作用下，傾角α=0°、5°、8°、12°，樁間距l(xiāng)=2.5d、3.0d、4.0d、5.0d、6.0d(d為樁徑，單位：m)群樁的承載性能，設(shè)計了20種工況，見表2。

表2 有限元模擬工況Table 2 Finite element simulation conditions

2 數(shù)值模擬結(jié)果分析

2.1 數(shù)值模擬與室內(nèi)試驗結(jié)果對比

在傾角α=0°，樁間距l(xiāng)=3.0d時，各級豎直荷載下群樁承臺頂荷載-沉降(Q-S)曲線如圖3。

圖3 α=0°、l=3.0d時，承臺頂荷載-沉降曲線Fig. 3 Q-S curves on top of cap when α=0° and l=3.0d

由圖3可見，當(dāng)Q< 6 kN時，數(shù)值模擬與室內(nèi)試驗的Q-S曲線基本重合，呈線性關(guān)系；當(dāng)Q≥6 kN時，荷載與承臺頂沉降呈非線性關(guān)系，且在相同荷載作用下，數(shù)值模擬的沉降值均小于試驗得出來的沉降值。數(shù)值模擬的Q-S曲線整體上與試驗的Q-S曲線相吻合。

2.2 傾角、樁間距對承臺頂沉降的影響

不同傾角α、樁間距l(xiāng)下承臺頂荷載-沉降(Q-S)曲線如圖4。

由圖4可見：

1)當(dāng)荷載較小時，樁和土體均處于彈性狀態(tài)，Q-S曲線近似于一條斜直線，表現(xiàn)為線性關(guān)系；當(dāng)荷載增大到某一值時，靠近承臺底的土體由彈性轉(zhuǎn)變?yōu)樗苄?，且隨著荷載的繼續(xù)增大，塑性區(qū)不斷擴大，曲線表現(xiàn)出由線性向非線性變化。

2)無論是豎直群樁還是傾斜群樁，當(dāng)傾角α相同時，樁間距l(xiāng)越大，承臺頂沉降越小，如：α=5°時，l=5.0d所對應(yīng)的沉降值比l=2.5d、3.0d、4.0d所對應(yīng)的沉降值，分別減小了336.8%、238.0%、115.3%。分析原因是：隨著樁間距的增大，樁與樁之間應(yīng)力疊加效應(yīng)即群樁效應(yīng)減弱，使得樁間土的承載力提高，承臺頂?shù)某两迪鄳?yīng)減小。

3)加載前期，豎直樁基承臺頂沉降量大于傾斜樁基承臺頂沉降量，但兩者的沉降量相差不大；隨著豎直荷載的繼續(xù)增大，傾斜樁基承臺頂沉降量比豎直樁基承臺頂沉降量大，荷載越大，兩者的沉降量值相差越明顯。因此，考慮到工程中樁基承臺以承受豎直荷載為主，建議采用豎直群樁。

2.3 不同樁身傾角下各基樁樁側(cè)摩阻力沿樁身分布特征分析

豎直荷載作用下，樁土相互作用時會在樁側(cè)產(chǎn)生側(cè)摩阻力。沿樁某段的樁側(cè)摩阻力qsi可按式(1)計算：

(1)

式中：qsi為樁身第i-1截面與第i截面間的樁側(cè)摩阻力平均值，kPa；u為樁身周長，mm；li為樁身第i-1截面與第i截面間的樁長，m；Qi-1、Qi分別為樁身第i-1截面和第i截面處的軸力，N。

中樁、角樁、邊樁的樁側(cè)摩阻力qs隨樁體深度h變化的關(guān)系曲線如圖5。

圖5 豎直樁基和傾斜樁基樁側(cè)摩阻力與樁體深度關(guān)系曲線Fig. 5 Relationship curve between pile side friction and pile depth of vertical pile foundation and inclined pile foundation

由圖5可見：

1)各基樁樁頂下約0.2 m深度范圍內(nèi)的樁側(cè)摩阻力表現(xiàn)為負(fù)值，再往下，樁側(cè)摩阻力表現(xiàn)為正值，且隨豎直荷載的增加呈非線性增大。分析原因是：由于第一層土為淤泥質(zhì)黏土，土質(zhì)較差，壓縮模量小，因此，較小豎直荷載作用下樁周土體的沉降量大于樁體的沉降量，土體對樁體產(chǎn)生一個向下的摩擦力，即樁側(cè)負(fù)摩阻力；隨著豎直荷載的增大，樁體沉降量逐漸大于樁周土體沉降量，土體對樁體產(chǎn)生一個向上的摩擦力，即樁側(cè)正摩阻力；樁土之間的相對位移隨豎直荷載的增大而增大，從而樁側(cè)負(fù)摩阻力表現(xiàn)為非線性增大的趨勢。

2)相同荷載作用下，傾斜樁基的各基樁樁側(cè)摩阻力qs的工作特性不同于豎直樁基的，傾斜樁基表現(xiàn)為：qs傾斜中樁>qs傾斜角樁>qs傾斜邊樁；豎直樁基表現(xiàn)為：qs豎直角樁>qs豎直邊樁>qs豎直中樁。分析原因是：在傾斜樁基中，隨著樁入土深度的不斷增加，角樁、邊樁與中樁之間的水平距離不斷增大，樁體受土體的約束作用則不斷減小，群樁效應(yīng)減弱，土的承載力提高，此時中樁承擔(dān)大部分荷載，其樁側(cè)摩阻力最先發(fā)揮作用，其次是角樁，最后是邊樁；而在豎直樁基中，由于角樁受到土體的約束作用最小，荷載大部分由角樁最先開始承擔(dān)，其樁側(cè)摩阻力最先發(fā)揮作用，其次是邊樁，最后是中樁。

2.4 樁身豎直位移云圖分析

當(dāng)承臺中心承受的豎直荷載Q分別為4、8 kN時，在不同樁間距l(xiāng)下，豎直樁基及傾斜樁基的樁身豎直位移U剖面云圖如圖6。

圖6 豎直樁基及傾斜樁基樁身豎直位移云圖(單位:m)Fig. 6 Nephogram of vertical displacement of vertical pile foundation and inclined pile foundation

由圖6可見：

1)在豎直逐級加載的過程中，各樁基豎直方向上的變形均隨著承臺頂荷載的增大而增大，樁身上部率先產(chǎn)生壓縮變形，變形均從樁頂向樁端傳遞，且群樁基礎(chǔ)的最大豎直位移均出現(xiàn)在樁頂邊緣，最小位移均出現(xiàn)在樁端邊緣。

2)在豎直樁基中，中樁處于群樁的包圍之中，受鄰近樁的影響最大，由于應(yīng)力疊加，最易和土一起發(fā)生整體沉降，同時在外側(cè)基樁的遮攔作用下，角樁在豎直方向上的變形最小，邊樁次之，中樁最大；隨著樁間距增大，樁間土對中樁的攜帶作用減弱，使得各基樁豎直變形差異逐漸縮小。

3)在傾斜樁基中，樁型采用角樁與邊樁背離中心樁傾斜布置，沿樁身向下豎直方向上，中心樁與周邊樁樁距與樁徑比逐漸增大；隨著斜樁傾斜角度的增大，群樁中各基樁相互作用逐漸變小，應(yīng)力疊加逐漸被削弱；樁間距l(xiāng)=2.5d、3.0d時，由于樁間距小，角樁與邊樁排列緊密，中樁與外側(cè)基樁之間排列相對彌散，樁周土體對邊樁的攜帶作用大于整個群樁對中樁的攜帶作用，因此，豎直位移U角樁

無論是豎直樁基還是傾斜樁基，均有以下現(xiàn)象：隨著樁間距的增大，群樁中各基樁相互作用減弱，樁土之間相互約束減弱，各基樁變形性狀相似。說明樁型的布置和樁間距的變化是影響傾斜群樁基礎(chǔ)各基樁豎直變形的主導(dǎo)因素，不同樁間距基樁豎直變形會發(fā)生顯著的變化。

2.5 傾角、樁間距對樁身彎矩的影響

圖7為樁間距l(xiāng)=3.0d、5.0d時，傾斜樁基的角樁樁身彎矩M隨樁體深度h的變化曲線。

圖7 樁間距l(xiāng)=3.0d、5.0d時，傾斜樁基的角樁樁身彎矩M隨樁體深度h的分布Fig. 7 Distribution of bending moment M of corner pile of inclined pile foundation with pile depth h, when l=3.0d, 5.0d

由圖7可見：

1)樁頂下0.09 m深度范圍內(nèi)，樁基下側(cè)受拉，樁頂產(chǎn)生較大的正彎矩，隨著傾斜角度的增加，樁頂附近彎矩急劇增加，樁身彎矩沿深度逐漸減小。

2)在樁身0.09 m處，開始出現(xiàn)反彎矩。主要原因是：由于上部豎直荷載的作用，傾斜樁基樁身整體在沉降的過程中，不僅樁體產(chǎn)生豎直位移，而且樁身背離樁基中心向兩側(cè)土體移動并發(fā)生彎曲，樁身上側(cè)受拉，產(chǎn)生負(fù)彎矩。

3)當(dāng)樁身傾角α=12°時，在樁身0.67～0.90 m處再次出現(xiàn)反彎矩，彎矩為正。對于超長傾斜樁基，當(dāng)樁周土較硬，樁體入土深度較大時，樁身產(chǎn)生撓曲變形，且變形隨荷載的增大而增大。樁端與土體相當(dāng)于固接，限制了樁端的撓曲變形，樁端處產(chǎn)生負(fù)彎矩。

4)樁間距相同時，樁身傾角α對樁身彎矩突變位置并沒有影響；相同荷載作用下，樁間距l(xiāng)=5.0d時斜樁彎矩突變的位置，比l=3.0d時的位置略淺；樁頂彎矩隨著樁間距的增加而減小，樁身中部彎矩繞樁中軸線呈上大下小的螺旋分布，螺旋狀隨基樁傾角的增大而增大，同時樁身下部出現(xiàn)較明顯的反彎曲現(xiàn)象，樁底彎矩均接近于0。

3 結(jié) 論

筆者通過有限元軟件模擬分析，研究了樁身傾角α、樁間距l(xiāng)對樁基豎直承載特性的影響。得到主要結(jié)論如下：

1)在較小豎直荷載作用下，豎直樁基與傾斜樁基的Q-S曲線均為線性，兩種樁的沉降值相差不大；隨著豎直荷載的增大，兩者的Q-S曲線逐漸向下彎曲，呈現(xiàn)非線性狀態(tài)，該階段傾斜樁基沉降量比豎直樁基的大。當(dāng)豎直荷載、傾角相同時，承臺頂沉降值隨著樁間距的增大而減小。

2)在相同豎直荷載作用下，各基樁樁側(cè)摩阻力qs大小關(guān)系，對于傾斜樁基為：qs傾斜中樁>qs傾斜角樁>qs傾斜邊樁；對于豎直樁基為：qs豎直角樁>qs豎直邊樁>qs豎直中樁。

3)在相同豎直荷載作用下，豎直樁基中，各基樁豎直位移U大小關(guān)系為：U中樁>U邊樁>U角樁，位移差隨著樁間距的增大而減??；傾斜樁基中，各基樁豎直位移U大小關(guān)系，對于較小樁間距為：U角樁

4)在豎直荷載作用下，樁身彎矩曲線整體或局部呈現(xiàn)“S”型曲線。當(dāng)樁間距相同時，樁身彎矩突變位置與傾角無關(guān)。增大樁間距有利于減小樁頂彎矩。