砂土地基中混凝土實(shí)心方樁水平承載特性分析

楊華勛，王福臨，江 軍，張炳磊

(魯東大學(xué) 巖土工程重點(diǎn)實(shí)驗(yàn)室， 山東 煙臺(tái) 264025)

由于沿海地區(qū)地質(zhì)復(fù)雜，抗腐蝕要求很高，混凝土實(shí)心方樁已在海洋結(jié)構(gòu)工程中廣泛使用[1]。海洋環(huán)境中，樁基礎(chǔ)不僅承受上部結(jié)構(gòu)傳遞的豎向荷載，還承受海風(fēng)、波浪等作用下產(chǎn)生的橫向水平荷載。水平荷載的承擔(dān)一般是樁自身水平承載性能和樁周土體共同完成[2-6]。鄭剛等[7]基于數(shù)值模擬研究了豎向及水平荷載加載大小、次序?qū)螛冻休d性能的影響，在模擬中混凝土和鋼筋分別采用塑性損傷模型和彈性模型，模擬結(jié)果顯示柔性長(zhǎng)樁表現(xiàn)出與剛性、半剛性短樁不同的承載特性。蔡忠祥等[8]基于類(lèi)似計(jì)算方法探討了樁體配筋率、樁頭約束條件及豎向荷載大小對(duì)于小直徑混凝土灌注樁的影響。Conte等[9]、Ziccarelli等[10]通過(guò)現(xiàn)場(chǎng)原位試驗(yàn)與數(shù)值模擬結(jié)果的比較，表明在承受較大水平荷載時(shí)，考慮混凝土損傷和鋼筋屈服才能合理反映鋼筋混凝土單樁的受力性狀。范慶來(lái)等[11]采用三維非線性有限元方法對(duì)于水平荷載作用下圓形單樁力學(xué)響應(yīng)進(jìn)行了系統(tǒng)的數(shù)值分析，并根據(jù)計(jì)算結(jié)果提出了一種新的單樁水平承載力容許值的確定方法。以上水平承載樁的研究多是針對(duì)常規(guī)實(shí)心圓形樁、管樁，對(duì)于實(shí)心方樁發(fā)生混凝土保護(hù)層裂開(kāi)及鋼筋塑性屈服后的分析探討則不是很多。本文以混凝土實(shí)心方樁為研究對(duì)象，運(yùn)用大型通用有限元程序ABAQUS來(lái)構(gòu)建數(shù)值計(jì)算模型，研究樁周土強(qiáng)度參數(shù)、豎向荷載大小及樁體配筋率參數(shù)對(duì)于方樁橫向承載性能的影響。

1 有限元數(shù)值模型建立與驗(yàn)證

Ziccarelli等[10]針對(duì)混凝土實(shí)心方樁進(jìn)行了橫向承載力現(xiàn)場(chǎng)試驗(yàn)，得到了樁頂位移與水平力之間的關(guān)系等測(cè)試數(shù)據(jù)。本文為了驗(yàn)證有限元模型的合理性，對(duì)此試驗(yàn)進(jìn)行了分析與比對(duì)。在模型中，方樁的截面大小為1.2 m×2.8 m，埋置深度為11 m，樁頭高出地面1.5 m，樁身采用φ30鋼筋單層配筋，長(zhǎng)短邊為11×7均勻配置，保護(hù)層厚度取為4 cm。

試驗(yàn)場(chǎng)地地基土類(lèi)為砂土，分為兩層，其內(nèi)摩擦角分別為37°和42°，在計(jì)算中采用Mohr-Coulmb彈塑性模型模擬?；炷梁弯摻罘謩e使用彌散開(kāi)裂模型和理想彈塑性Mises模型[11]。地基土體和方樁樁體均使用縮減積分三維實(shí)體單元離散，鋼筋離散為桿單元，內(nèi)置于樁體網(wǎng)格中。樁土之間的接觸面采用的是硬接觸中的罰函數(shù)方法，界面摩擦系數(shù)取為f=tan(2φ/3)?？紤]到方樁的實(shí)際形狀，半無(wú)限地基土體截?cái)酁榱⒎襟w區(qū)域，平面上長(zhǎng)60 m，寬40 m，深度為20 m?？紤]到樁身臨近土體應(yīng)變較大，該處網(wǎng)格劃分較密，隨著離樁體的距離增大網(wǎng)格逐漸趨于稀疏。經(jīng)過(guò)試算，這樣的網(wǎng)格劃分情況滿足計(jì)算精度的要求。模型建立完成后，首先進(jìn)行地應(yīng)力平衡以消除土體自重產(chǎn)生的變形，然后逐級(jí)施加水平力進(jìn)行加載，進(jìn)而得到不同加載階段樁頭側(cè)移量。

圖1為得到的樁頂側(cè)移s與樁頂水平力Q關(guān)系曲線。經(jīng)過(guò)比較看到，當(dāng)水平荷載小于1.6 MN時(shí)，模型計(jì)算結(jié)果略小于試驗(yàn)結(jié)果，而當(dāng)水平荷載大于1.6 MN時(shí)，計(jì)算結(jié)果略大于試驗(yàn)結(jié)果。從總體來(lái)看，試驗(yàn)實(shí)測(cè)值與計(jì)算值比較吻合，從而表明本文所采用的的有限元模型是比較可靠的。

圖1橫向荷載施加時(shí)方樁頂部的側(cè)移曲線

圖2顯示的是不同水平荷載時(shí)，實(shí)心方樁整體的側(cè)向位移曲線。在水平荷載較小時(shí)，樁體基本上是一個(gè)剛體旋轉(zhuǎn)，旋轉(zhuǎn)中心的深度約為泥面以下10 m。當(dāng)樁頂施加的水平荷載較高時(shí)，樁體從樁頂?shù)?.5 m埋深處出現(xiàn)比較明顯的彎曲。這是因?yàn)樵谠搮^(qū)域混凝土逐漸產(chǎn)生累積損傷，鋼筋承受轉(zhuǎn)移而來(lái)的荷載而進(jìn)入屈服階段，因此誘發(fā)較大的彎曲變形。

圖2不同橫向載荷對(duì)實(shí)心方樁整體的側(cè)向位移影響

圖3顯示的是水平荷載為不同特征值時(shí)方樁損傷區(qū)域擴(kuò)展情況。

圖3方樁樁身隨著橫向荷載的增大

損傷區(qū)域(灰色區(qū)域)擴(kuò)展

在圖3(a)中，在橫向作用力小于2.23 MN時(shí)，方樁處于彈性狀態(tài)，樁體混凝土未產(chǎn)生損傷。當(dāng)水平荷載增加到3.01 MN時(shí)，如圖3(b)所示，在樁身的受拉側(cè)混凝土出現(xiàn)損傷區(qū)，而后隨著荷載的增大，樁體受拉側(cè)的損傷區(qū)逐步擴(kuò)展，如圖3(c)和3(d)所示的灰色區(qū)域。將數(shù)值分析獲得的損傷區(qū)域與Ziccarelli等的原位試驗(yàn)[10]進(jìn)行對(duì)比，除了有小范圍上移外，其他基本相符。因此在模型計(jì)算中合理考慮混凝土損傷與鋼筋屈服將更符合實(shí)際工況，尤其是水平荷載較大階段[7]。

2 水平承載特性影響因素分析

為了研究實(shí)心方樁水平承載特性，進(jìn)行如下變動(dòng)參數(shù)研究。實(shí)心方樁樁長(zhǎng)取為24.9 m，其中樁頭高出地面約0.9 m，其截面邊長(zhǎng)為0.5 m，樁身配筋情況為單層配筋，均勻排列，保護(hù)層厚度規(guī)定為4 cm。地基土持力層分兩層，其中上層土厚為7 m，下層砂土厚度為24 m。其它參數(shù)選取與前述情況相同，有限元模型網(wǎng)格劃分情況如圖4所示。

圖4有限元模型

2.1 樁周土體內(nèi)摩擦角的影響

在分析中，樁周土彈性模量取為E=100 MPa保持不變，上層土內(nèi)摩擦角φ進(jìn)行參數(shù)變動(dòng)，涵蓋了從松砂到密砂的情況，分別取32.5°、35°、37.5°和40°進(jìn)行計(jì)算。圖5為不同情況下水平荷載Q與樁頂側(cè)移s之間的關(guān)系曲線。

圖5土體內(nèi)摩擦角對(duì)樁頂位移的影響

從圖5中可以看出，內(nèi)摩擦角對(duì)樁土相互作用影響比較顯著，內(nèi)摩擦角增大，樁基水平承載力明顯提高。這也意味著，在相同荷載水平下，上層土體內(nèi)摩擦角φ增大，樁頂水平位移大大降低。例如當(dāng)水平荷載是300 kN時(shí)，φ=37.5°情況下的樁頭側(cè)移量比φ=32.5°工況降低12.2 mm。

圖6提取了不同內(nèi)摩擦角φ情況下樁體受拉側(cè)最大彎矩截面處鋼筋應(yīng)力隨著水平荷載變化的Q-σs曲線。根據(jù)范慶來(lái)等的建議方法[11]，根據(jù)圖6能夠判斷出方樁水平臨界荷載和極限荷載，如表1所示。由于實(shí)心方樁目前主要用在房建工程中，根據(jù)相應(yīng)規(guī)范其水平位移限值取為10 mm。結(jié)合圖5，可以將樁頭位移達(dá)到10 mm時(shí)的水平荷載作為其容許值Q10mm，一起列在表1中以便于比較。通過(guò)比較可以看到，在所考慮內(nèi)摩擦角φ范圍內(nèi)，臨界荷載的大小與Q10mm相差不大，并且Q10mm值比臨界荷載稍大。根據(jù)侯勝男等[12]的建議，在這種情況下，需要取較小值，即臨界荷載作為水平承載力容許值，取極限荷載作為承載力極限值。通過(guò)表1可見(jiàn)，內(nèi)摩擦角從φ=32.5°增加到φ=40°時(shí)，方樁的承載力容許值增加了38.2%，承載力極限值增加了37.1%。

圖6 內(nèi)摩擦角對(duì)鋼筋應(yīng)力Q-σs的影響

2.2 豎向荷載水平的影響

實(shí)際工程中，方樁基礎(chǔ)承受豎向荷載和水平荷載的共同作用，豎向荷載對(duì)其水平承載特性有較大影響。保持彈性模量E=100 MPa，內(nèi)摩擦角φ=35°不變，在樁頂施加豎向荷載分量分別是0、0.2、0.4、0.6和0.8倍的豎向承載力極限值(對(duì)該模型單純施加豎向荷載，確定出方樁豎向承載力極限值Vult=2 050 MN)，再逐級(jí)加載水平荷載，得到實(shí)心方樁樁頂側(cè)移與樁頭水平荷載之間的Q-s關(guān)系曲線如圖7所示。

圖7豎向作用力對(duì)樁頂位移的影響

由圖7可以看到，當(dāng)荷載小于150 kN時(shí)，Q-s關(guān)系曲線處在線彈性區(qū)間，豎向荷載水平對(duì)于方樁樁頂側(cè)移基本不產(chǎn)生影響。隨著水平荷載的增大，豎向荷載水平的變化對(duì)于Q-s曲線產(chǎn)生了比較明顯的影響。值得注意的是，增大豎向荷載水平，相同水平荷載引起的樁頂?shù)膫?cè)移量趨于降低，這個(gè)結(jié)論與已有研究[13-16]是相統(tǒng)一的。這主要是因?yàn)樨Q向荷載的變大使得樁周臨近土體中側(cè)向應(yīng)力變大，進(jìn)而導(dǎo)致樁土界面上法向應(yīng)力變大，因此樁側(cè)受到的摩阻力也相應(yīng)變大，最終促使樁的水平承載特性有所提高。

圖8顯示的是樁頂在水平荷載350 kN時(shí)，不同豎向荷載產(chǎn)生的方樁樁身彎矩圖?？梢钥吹剑Q向荷載變大時(shí)，樁身的最大正彎矩值隨之變大，其位置在泥面以下大約2.5 m處，基本保持不變，負(fù)彎矩在大約泥面以下7.5 m處達(dá)到峰值，而且豎向荷載水平對(duì)負(fù)彎矩峰值影響不很明顯。

圖8豎向作用力對(duì)于方樁樁身彎矩的作用

通過(guò)圖7、圖8綜合分析，豎向荷載水平能在一定程度上改善實(shí)心方樁的水平承載特性，但考慮到豎向荷載分量的存在會(huì)導(dǎo)致樁體正彎矩峰值明顯變大，增加配筋量，因此在實(shí)際工程中，砂性土地基上方樁的豎向荷載水平最優(yōu)值為0.4Vult～0.6Vult，這與砂土地基上淺基礎(chǔ)承載性狀有類(lèi)似之處[17]。

2.3 樁體配筋率的影響

隨著樁頭水平荷載的增加，樁周土體將發(fā)生累積塑性應(yīng)變，受拉側(cè)混凝土逐漸出現(xiàn)損傷裂隙而退出工作，其承擔(dān)的應(yīng)力將轉(zhuǎn)移到該區(qū)域鋼筋來(lái)負(fù)擔(dān)。因此方樁基礎(chǔ)配筋率對(duì)于方樁基礎(chǔ)的水平承載特性至關(guān)重要。樁周土體彈性模量E=100 MPa，內(nèi)摩擦角φ=35°，不同的配筋率如表2所示。

表2 方樁的配筋率

圖9給出的是不同的配筋率情況下樁頂側(cè)移Q-s曲線。從圖9可以看到，增大樁體配筋率，可以有效降低相同水平荷載時(shí)的樁頂側(cè)移，改善實(shí)心方樁承載性能。在樁頂水平荷載較小時(shí)，樁身受拉側(cè)混凝土沒(méi)有出現(xiàn)裂縫，配筋率對(duì)Q-s曲線影響較小，但水平荷載較大時(shí)，增加配筋率能夠有效地抑制樁頭側(cè)移。例如在Q=350 kN時(shí)，配筋率由1.16%增加到3.03%時(shí)，砂土中方樁樁頭位移量減小了36.4%。

圖9不同配筋率情況下方樁樁頂位移曲線

圖10給出了不同配筋率時(shí)鋼筋的Q-σs曲線。從圖10中可見(jiàn)，增加樁身配筋率，會(huì)很大程度上提高實(shí)心方樁極限荷載值。對(duì)于配筋率ρ較小的情況，如配筋率ρ=1.16%時(shí)，臨界荷載為126 kN，Q10mm=135 kN，所以此時(shí)水平承載力容許值受混凝土損傷控制。然而對(duì)于配筋率較大的情況，例如配筋率ρ=3.03%時(shí)，臨界荷載是221 kN，Q10mm=196 kN，所以對(duì)于高配筋率的方樁來(lái)說(shuō)，其承載力受水平位移限值決定[18]。

圖10配筋率對(duì)于鋼筋應(yīng)力的影響

3 結(jié) 論

(1) 相同水平荷載作用下，增大內(nèi)摩擦角可有效增強(qiáng)實(shí)心方樁的水平承載特性，從φ=32.5°增加到φ=40°時(shí)，方樁的水平荷載容許值增長(zhǎng)了38.2%。

(2) 豎向荷載增大會(huì)降低相同水平荷載情況下實(shí)心方樁樁頭位移，但同時(shí)導(dǎo)致實(shí)心方樁樁身正彎矩峰值變大，因此砂性土中方樁豎向荷載水平建議取為豎向極限荷載的0.4倍～0.6倍。

(3) 方樁配筋率由1.16%變化到3.03%時(shí)，方樁樁頂處位移量減小了36.42%，樁頂?shù)膫?cè)移得到顯著降低。低配筋率方樁的水平承載力容許值可取為臨界荷載，該值受樁體開(kāi)裂控制，較高配筋率方樁則由樁頭水平位移限值來(lái)控制。