曹衛(wèi)平,李 賀,陶 鵬

(1.西安建筑科技大學土木工程學院,陜西 西安 710055; 2.陜西省巖土與地下空間工程重點實驗室,陜西 西安 710055)

輸電塔、碼頭和海洋工程中常用斜樁抵抗風、海浪、土壤壓力等較大的橫向荷載[1]。與直樁相比,斜樁具有更高的側向水平承載力[2],但由于斜樁樁身的軸線與建筑物所在平面并不垂直,導致斜樁在受到外部荷載作用時的承載變形特性與直樁有較大的不同。因此,對斜樁開展在水平荷載H和扭矩T耦合作用(H-T耦合作用)下的荷載變形特性和承載性能的研究對斜樁基礎的設計具有重要意義。

樁基承載變形問題是國內外有關學者廣泛關注的問題[3-5]。例如:Zhang等[6]通過離心模型試驗研究了循環(huán)荷載作用下樁的變形特性,發(fā)現隨著循環(huán)次數的增加,樁周土體的剛度逐漸減小;Ramadan等[7]通過對不同樁身傾角下受荷斜樁的樁身彎矩及位移分析,認為上拔荷載與水平荷載對樁身內力有著較明顯的影響;Meyerhof等[8]通過模型試驗研究了傾斜荷載作用下斜樁的變形性狀,給出了樁頂水平位移的影響曲線;呂凡任[9]通過模型試驗和現場試驗開展了對傾斜荷載作用下斜樁基礎的承載性狀研究;曹衛(wèi)平等[10-11]基于砂土中的斜樁模型試驗和數值模擬軟件,分別研究了豎向荷載和水平荷載作用下斜樁承載性狀,分析了樁身傾角和長徑比對斜樁承載力的影響;趙爽等[12]基于數值模擬軟件,研究了砂土中各方向水平荷載與傾斜面斜交作用下斜樁單樁的變形及承載特性。也有部分學者研究了H-T耦合作用下樁的承載特性。例如:鄒新軍等[13]進行了H-T耦合作用下直樁的模型試驗,并探討了施加荷載的順序對單樁的影響,得出無論是先施加水平荷載還是扭矩,都會使得單樁的扭矩或水平承載力降低,但先施加水平荷載對樁的承載力影響比先施加扭矩影響更大;郭沛翰等[14]在鄒新軍等[13]兩種荷載作用的基礎上添加了豎向荷載作用,完成了直樁室內砂土地基的模型試驗,得到了3種荷載共同作用下對單樁的影響曲線和承載力計算公式。

目前對斜樁的承載變形研究大多只分析單個荷載作用下的變形情況,對涉及兩個及多個荷載耦合作用下斜樁變形的研究較少,而H-T耦合作用對樁的影響研究大多只是針對直樁,對斜樁開展的相關研究較少。本文采用ABAQUS數值軟件對斜樁在H-T耦合作用下的承載特性進行數值模擬,分析H-T耦合作用下樁身傾角對斜樁承載特性的影響,研究斜樁樁-土接觸壓力及摩阻力沿樁身的分布特征。

1 有限元模型

1.1 數值模型

本文所建立的數值模型如圖1所示,為消除土體邊界效應,樁周土計算范圍取30D(D為樁徑),樁端底部土層厚度取10D[15-16]。樁頂承受H-T耦合作用,水平荷載和扭矩為等比例同時施加,其中水平荷載垂直作用于斜樁樁身傾斜面,即水平荷載作用于斜樁所在平面以外,扭矩方向為z軸正向。

圖1 有限元模型與計算網格

采用ABAQUS數值軟件對H-T耦合作用下的斜樁承載性狀進行模擬分析,選取線彈性模型模擬樁體,Drucker-Prager模型模擬土體變化。樁-土界面相互作用通過ABAQUS數值軟件中樁-土界面接觸算法實現,選擇樁側與樁端表面為主控面,土側與樁端土表面為從屬面。樁-土接觸界面法向作用采用硬接觸的方法模擬,切向作用采用罰函數模擬,樁-土接觸界面間的摩擦系數μ采用μ=tan(0.75φ)計算,其中φ為樁周土的內摩擦角[15]。樁身網格密度均勻布置,土體網格密度根據與樁身的距離由近至遠沿徑向由密到疏布置。樁身及土體單元均采用三維實體單元C3D8R[17]。

1.2 數值模型驗證

為驗證數值模型材料本構關系選取和計算參數設置的合理性,分別選取呂凡任等[18]和王書行[19]進行的現場載荷試驗和模型試驗進行驗證分析。

呂凡任等[18]進行了軟土地基上斜樁的現場試驗研究,試驗用斜樁樁長為15m,樁徑為0.25m,樁身傾角為10°。在進行數值模擬時,參考文獻[13-14],取鋼樁密度為7.9t/m3;各土層的泊松比均取為0.3,土體剪脹角取6.5°?，F場試驗與數值模擬的對比結果如圖2所示。

圖2 現場試驗與數值模擬結果對比

王書行[19]開展了飽和粉土地基中扭矩作用下直樁室內模型試驗。試驗用直樁材料為鋼管,樁長為6m(入土深度為3.41m),樁徑為114mm,壁厚為4.5mm;土層密度為1.7t/m3,有效內摩擦角為30°,有效黏聚力為1kPa,取粉土地基的壓縮模量為9.9MPa,泊松比為0.3,模型土體彈性模量取壓縮模量的4倍。模型試驗與數值模擬的對比結果如圖3所示。

圖3 模型試驗與數值模擬結果對比

圖2和圖3分別為荷載(Q)-沉降(s)曲線和樁頂扭矩(T)-扭轉角(ψ)曲線,結果表明,數值模擬結果與現場試驗和模型試驗的結果吻合較好,說明本文建立的數值模型是合理的,能較好地反應樁基與地基土之間的相互作用。

1.3 數值模擬方案

本文基于呂凡任等[18]的研究背景,對H-T耦合作用下斜樁的承載變形特性進行有限元數值模擬分析,模擬樁樁徑為0.5m,樁長為15m,密度為2.4t/m3,泊松比為0.2,彈性模量為31.5GPa。土層參數如圖4所示,各層土的泊松比均取為0.3,土體剪脹角可取為6.5°[11]。

圖4 數值模擬土體示意圖

斜樁數值模擬共分為4組,分別模擬樁身傾角、樁-土剛度比、荷載偏心距、樁身長徑比等4種影響因素對斜樁承載變形特性的響應,模擬方案如表1所示,除第一組包含4根樁外,其他3組均只含3根樁。

表1 數值模擬方案

2 直樁與不同樁身傾角斜樁的承載變形特性

2.1 樁頂扭轉角和位移

圖5為直樁和不同樁身傾角斜樁樁頂扭矩-扭轉角曲線,當曲線出現急劇上升變化時即可認為樁頂扭矩達到極限扭矩而導致樁身破壞。由圖5可見,樁身傾角越大,斜樁樁頂極限扭矩和極限扭轉角越大。定義α、β分別為斜樁對直樁極限扭矩和極限扭轉角的提高幅度,表2為不同樁身傾角下斜樁的極限扭矩和極限扭轉角相對直樁的提高幅度。圖5和表2的結果表明,與直樁相比,斜樁能提高樁身受扭承載力,且樁身傾角越大,提高幅度越大。

表2 樁身傾角對斜樁樁頂扭轉承載變形特性的影響

圖5 樁頂扭矩-扭轉角曲線

從表2和圖5中還可以看出,通過增大樁身傾角可以增加樁身極限扭矩,但同時要注意樁身傾角越大的斜樁其樁頂扭轉角相比其他斜樁也越大,如20°斜樁的極限扭轉角相對直樁的提高幅度為63%。在實際工程中可以通過增大一定的樁身傾角來增加樁的受扭承載力,但是要考慮扭轉角的極限值,同時更要考慮更大變形的不利影響。

直樁和不同樁身傾角斜樁樁頂水平荷載(H)-位移(s)曲線如圖6所示。結合圖5并分析曲線拐點可知,在樁身達到最大扭矩前,樁頂水平位移與樁頂水平荷載呈線性正相關關系;在樁頂水平荷載相同的情況下,樁頂水平位移與樁身傾角呈正相關關系。

圖6 樁頂水平荷載-位移曲線

斜樁樁頂在受到H-T耦合作用時,參考趙春風等[20]研究成果和JGJ106—2014《建筑基樁檢測技術規(guī)范》,取泥面處單樁水平位移達到0.1D時的水平荷載為單樁水平極限承載力。由表2、圖6以及鄒新軍等[21]的試驗結果可知,斜樁達到極限扭矩時,相應水平荷載作用下的樁頂水平位移均遠小于0.1D,即此時尚未達到斜樁的水平極限承載力,因此在H-T耦合作用下,斜樁樁身為受扭破壞,即相對水平極限承載力而言樁身的受扭承載力較低。

2.2 樁身位移和內力

選取H=3.2kN、T=4kN·m,得到H-T耦合作用下不同樁身傾角時樁身側向水平位移沿樁身相對深度(樁身沿z軸入土深度與總樁長L的比值)的分布曲線如圖7所示(位移以沿y軸正向為正)。從圖7可以看出,隨著樁身相對深度的增大,直樁和不同樁身傾角斜樁樁身側向水平位移發(fā)展區(qū)均分布在樁頂部以下樁身相對深度0.3以內。另外,最大側向位移均出現在樁頂,隨后隨著深度逐漸減小。直樁和不同樁身傾角斜樁的反彎區(qū)基本都處于3D～9D之間,斜樁反彎區(qū)沿樁身相對深度范圍分布與直樁相近,范圍基本一致,且反彎區(qū)最大位移與樁頂水平位移相比很小,可忽略不計。這表明在H-T耦合作用下,樁身側向水平位移主要在樁周上部區(qū)域有變化。

圖7 樁身側向水平位移沿樁身相對深度的分布曲線

圖8為H=3.2kN、T=4kN·m時,在H-T耦合作用下直樁和不同樁身傾角斜樁樁身扭矩沿樁身相對深度的分布曲線。由圖8可知,直樁和不同樁身傾角斜樁的扭矩沿樁身分布,且隨著樁身相對深度的增大,樁身扭矩從樁頂的最大值一直逐漸減小,到樁端減到最小,接近于0。這是由于樁側摩阻力發(fā)揮抵抗作用使得樁身扭矩逐漸減小,在相同的深度下,樁身傾角越大,樁所受的扭矩越小。因此工程上可以通過提高一定的樁身傾角來增強樁身受扭能力。

圖8 樁身扭矩沿樁身相對深度的分布曲線

圖9為H=3.2kN、T=4kN·m時,在H-T耦合作用下直樁和不同樁身傾角斜樁樁身彎矩(M)沿樁身相對深度的分布曲線。由圖9可知,斜樁的樁身彎矩均大于直樁,且樁身傾角越大,樁身彎矩越大。直樁和不同樁身傾角斜樁的樁身彎矩都是隨著樁身相對深度的增大先迅速增大到最大值,20°斜樁樁身彎矩最大,最大值出現在1.8D左右處,隨后樁身彎矩沿樁身相對深度急劇減小為0。與圖8結果對比可知,在樁身相對深度0.15范圍內,樁身彎矩對樁的影響較大,而扭矩則沿全樁均有分布。

圖9 樁身彎矩沿樁身相對深度的分布曲線

2.3 樁-土接觸壓力

由于樁頂水平荷載作用在斜樁樁身傾斜平面外,同時樁頂受到扭矩的影響,故導致樁-土接觸壓力不是軸對稱的,因此選擇沿斜樁樁身深度方向上的A、B、C、D 4條路徑(圖10)對樁-土接觸壓力進行有限元模擬。

圖10 樁-土接觸壓力路徑

圖11為15°斜樁在樁頂扭矩為4kN·m時樁-土接觸壓力(p)沿樁身相對深度的分布曲線,從圖11可知,不同路徑側的樁-土接觸壓力有著較為明顯的不同:在樁頂向下樁身相對深度0.3范圍內,由于H-T耦合作用,沿B、C路徑側土體有脫離樁體的趨勢,故B、C路徑側樁-土接觸壓力小于靜止土壓力,而沿A、D路徑側土體有受到擠密的趨勢,故A、D路徑側的靜止土壓力小于樁-土接觸壓力。這表明在H-T耦合作用下,B、C路徑側土體為卸荷,A、D路徑側土體受擠壓。從卸荷側來說,路徑C側樁-土接觸壓力小于B側,這說明相比較路徑B側來說C側土體脫離樁體的趨勢更大。從擠壓側來說,路徑A側的樁-土接觸壓力大于D側,這說明相比路徑D側來說A側土體受到的擠壓效應更大。

圖11 15°斜樁樁-土接觸壓力沿樁身相對深度的分布曲線

2.4 樁側環(huán)向摩阻力

圖12為樁頂在H-T耦合作用下,直樁和不同樁身傾角斜樁的樁側環(huán)向摩阻力(qsi)沿樁身相對深度的分布曲線,可以看出,斜樁樁側環(huán)向摩阻力分布與直樁明顯不同:直樁樁側環(huán)向摩阻力一開始為0,而斜樁在樁頂以下一開始就有較大的摩阻力;直樁樁側環(huán)向摩阻力沿樁身相對深度逐漸增大,而斜樁樁側環(huán)向摩阻力沿樁身相對深度增大先減小,在樁身相對深度0.3處又反向沿樁身相對深度逐漸增大。同一深度位置處,樁側環(huán)向摩阻力與樁身傾角成正比,即20°斜樁的樁側環(huán)向摩阻力最大。這表明樁側環(huán)向摩阻力的發(fā)揮程度與樁身傾角有關,傾角越大,發(fā)揮程度越大,且樁身上部區(qū)段的發(fā)揮程度要大于樁身下部區(qū)段。此外,對于樁身中下部區(qū)段(z/L=0.3～0.9),樁側環(huán)向摩阻力沿樁身相對深度增大的一個原因是樁-土接觸壓力增大。

圖12 樁側環(huán)向摩阻力沿樁身相對深度的分布曲線

2.5 不同因素對斜樁樁頂扭矩-扭轉角曲線的影響

圖13～15分別為不同樁-土剛度比、荷載偏心距和樁身長徑比影響下15°斜樁的樁頂扭矩-扭轉角曲線,圖中曲線均為陡升型,若以曲線急劇變化點對應的扭矩值作為斜樁的極限扭矩,可得到對應于不同影響因素的斜樁的最大扭轉角。由此可知,在一定樁頂扭轉角下,樁-土剛度比和荷載偏心距越大,斜樁所承受的扭矩越小,表明斜樁受扭承載力與樁-土剛度比和荷載偏心距呈負相關關系,但其影響程度較小。而在同一樁頂扭轉角下,斜樁極限扭矩與樁身長徑比呈正相關關系且影響較大,短樁的極限扭矩增幅較小,長樁的極限扭矩增幅較大。對于短樁,樁-土界面發(fā)生破壞,樁側環(huán)向摩阻力為其提供了受扭承載力,故短樁的極限扭矩增幅較小。對于長樁,土體內部發(fā)生破壞,其受扭承載力主要由樁側土體的內部剪切作用和樁-土界面的扭轉剪切作用共同提供。綜上,樁-土剛度比和荷載偏心距對斜樁極限扭矩的影響不大,而樁身長徑比則對其影響較為顯著。

圖14 不同荷載偏心距時的樁頂扭矩-扭轉角曲線

如圖15所示,對同一傾角(15°)的斜樁來說,若以同一樁頂扭轉角作為極限值,則明顯可以看出樁身長徑比對受扭承載力的影響,即樁身長徑比越大的斜樁其極限扭矩也越大,樁頂扭轉角給定的前提下,扭矩與樁身長徑比呈正相關關系。從圖15還可以看出,雖然樁身長徑比較大的樁其極限扭矩也較大,但其相應的樁頂扭轉角也較大,即極限扭矩提高的同時也使得極限扭轉角增大,因此在實際工程中可以通過增大樁身長徑比來提高樁身受扭能力,但是要注意樁頂扭轉角的極限值不能過大。

圖15 不同樁身長徑比時的樁頂扭矩-扭轉角曲線

3 結 論

a.H-T耦合作用下,斜樁樁身表現為受扭破壞,斜樁相比直樁有更高的受扭承載力,且樁身傾角越大,承載力越高。

b.在樁頂H-T耦合作用下,斜樁樁身扭矩沿樁身分布,不存在臨界樁長的問題。斜樁樁身扭矩小于直樁,最大值出現在樁頂,沿樁身向下逐漸減小,在樁頂下同一截面位置處,斜樁樁身扭矩與樁身傾角呈負相關關系。

c.在樁頂H-T耦合作用下,相比直樁,斜樁更有利于樁側環(huán)向摩阻力的發(fā)揮,樁身傾角越大,樁側環(huán)向摩阻力的發(fā)揮程度越大,且樁身上部區(qū)段的發(fā)揮程度要大于樁身下部區(qū)段。

d.斜樁極限扭矩與樁-土剛度比和荷載偏心距呈負相關關系,但兩者對極限扭矩影響程度較小;而樁身長徑比與斜樁極限扭矩呈正相關關系且影響程度較大。