傅文橋， 郭小剛， 陳雨琴， 李平均

(湘潭大學(xué) 土木工程與力學(xué)學(xué)院,湖南 湘潭 411105)

鉆孔灌注樁樁底沉渣對(duì)單樁承載性狀的影響分析*

傅文橋， 郭小剛*， 陳雨琴， 李平均

(湘潭大學(xué) 土木工程與力學(xué)學(xué)院,湖南 湘潭 411105)

伴隨鉆孔灌注樁在各個(gè)工程領(lǐng)域的廣泛應(yīng)用，樁底易產(chǎn)生一定厚度弱性沉渣的影響問題得到了很多學(xué)者的重視，并采用不同方法對(duì)樁底沉渣的影響性進(jìn)行了研究分析.憑借有限元軟件在巖土問題上的可操作性和全面性，本文采用數(shù)值模擬方法來分析樁底沉渣對(duì)單樁承載性狀的影響.為深入和具體了解沉渣的影響性，從沉渣厚度、沉渣物理力學(xué)性質(zhì)和灌注樁長(zhǎng)徑比三個(gè)方面進(jìn)行細(xì)致的研究分析.計(jì)算結(jié)果表明:沉渣厚度、沉渣彈性模量、沉渣內(nèi)摩擦角和樁長(zhǎng)徑比對(duì)單樁荷載-沉降曲線存在一定程度影響，其中沉渣厚度的影響存在一個(gè)臨界值，而沉渣粘聚力和泊松比的影響性并不明顯.這些結(jié)論可用于指導(dǎo)或輔助灌注樁設(shè)計(jì)、施工、樁底沉渣檢測(cè)和沉渣處理等.

沉渣；灌注樁；數(shù)值模擬；承載性狀

隨著經(jīng)濟(jì)建設(shè)的高速發(fā)展，鉆孔灌注樁由于具有較高的承載能力和控制沉降的可靠性且施工簡(jiǎn)便，得以廣泛應(yīng)用于建筑工程、公路橋梁工程和港口工程等.但在實(shí)際工程施工中，限于施工隊(duì)伍的素質(zhì)、施工工藝和技術(shù)水平等因素，灌注樁樁底往往極易形成一定厚度的性質(zhì)較樁底土弱的沉渣.從很多統(tǒng)計(jì)的工程資料中可以發(fā)現(xiàn)，樁底沉渣的存在會(huì)降低樁的承載力，對(duì)樁的承載力發(fā)揮有較大影響[1].這使得越來越多的學(xué)者投入到樁底沉渣對(duì)樁承載性狀的影響研究中去，研究表明樁底沉渣對(duì)樁承載力的影響一般可達(dá)10%～30%，甚至可能更高.

一直以來，很多學(xué)者做了相關(guān)內(nèi)容的研究：Wong[2]采用室內(nèi)離心機(jī)試驗(yàn)和FB-Pier程序的數(shù)值模擬對(duì)單樁和群樁中存在樁底沉渣等情況進(jìn)行了細(xì)致分析；Iskander[3]等通過靜載試驗(yàn)得到的數(shù)據(jù)結(jié)果來分析樁底夾土等對(duì)樁承載力的影響，研究認(rèn)為樁底沉渣的存在會(huì)降低樁的承載力，有較大影響；周紅波[4]運(yùn)用PLAXIS軟件進(jìn)行樁底沉渣和樁側(cè)泥皮的數(shù)值模擬，研究了樁底沉渣對(duì)樁的承載性狀的影響；袁振[5]采用荷載傳遞法來研究樁端條件對(duì)灌注樁荷載傳遞特性的影響.本文為獲得一個(gè)系統(tǒng)性、全面性、指導(dǎo)性的結(jié)論，在大型有限元軟件ANSYS平臺(tái)上進(jìn)行數(shù)值模擬，對(duì)樁底沉渣對(duì)單樁承載性狀的影響性進(jìn)行分析.

1 幾何模型

由于幾何形狀、材料性質(zhì)、位移邊界條件和區(qū)域應(yīng)力都滿足軸對(duì)稱條件，故有限元分析可采用軸對(duì)稱模型.有限元模型建立時(shí)，關(guān)于土體邊界范圍的取值，有些文獻(xiàn)[7,8]取有限元模型的計(jì)算區(qū)域深度為2L(L為樁長(zhǎng))，寬度為1L；Desai[9]在進(jìn)行砂土中單樁承載性狀分析時(shí)，取計(jì)算區(qū)域深度為1.5L，寬度為30R(R為樁半徑)；Randolph[10]建議取計(jì)算區(qū)域深度為2.5L，寬度為50R；文獻(xiàn)[11]認(rèn)為樁端以下深度要延伸20R～30R，模型寬度要超過20R～30R，或者超過樁長(zhǎng)L；胡鋮波[6]等用有限元軟件對(duì)樁土模型進(jìn)行分析，提出一種樁土模型所需的合理邊界，樁側(cè)土體邊界取24R，樁端土體深度取30R.參照以上國(guó)內(nèi)外文獻(xiàn)并經(jīng)ANSYS試算結(jié)果分析，本文取有限元模型寬度為24R，深度為30R和1.5L的較大值，這樣不僅能滿足計(jì)算結(jié)果的精度要求和提高計(jì)算效率，又能合理考慮到土的承載能力.其二維軸對(duì)稱模型和擴(kuò)展三維模型分別如圖1和圖2所示.

2 土體本構(gòu)模型和樁土接觸面處理

建模時(shí)樁體采用線彈性實(shí)體單元，而土體和沉渣采用ANSYS中自帶的D-P模型.M-C強(qiáng)度準(zhǔn)則能較好地描述巖土材料的破壞行為，但是它在三維應(yīng)力空間中的屈服面存在角點(diǎn)的奇異性，即在偏平面上的破壞曲線是一個(gè)不規(guī)則的六角形，導(dǎo)致數(shù)值計(jì)算時(shí)收斂緩慢.因此產(chǎn)生了許多光滑的曲線來逼近M-C強(qiáng)度曲線的模型，而其中，Drucker-Prager準(zhǔn)則[12，13]就是在偏平面上用圓來逼近.

Mohr-Coulomb強(qiáng)度準(zhǔn)則可表示為：

τ=σntanφ+c，

(1)

其中，τ為剪應(yīng)力；σn為受剪面上的法向應(yīng)力；c為凝聚力；φ為內(nèi)摩擦角.

將式(1)用應(yīng)力不變量可表示為：

(2)

其中，I1為主應(yīng)力的第1不變量；J2偏應(yīng)力第2不變量.

將式(2)整理后可得：

(3)

當(dāng)取式(3)中θσ=-π/6，即可得到D-P強(qiáng)度準(zhǔn)則的表達(dá)式：

(4)

樁土接觸選擇ANSYS中面-面接觸，與土體和沉渣接觸的樁體表面設(shè)置為TARGE169，與樁體接觸的土體面和與樁體接觸的沉渣面設(shè)置為CONTA172.

3 位移邊界條件和初始地應(yīng)力模擬

參考其他文獻(xiàn)[14，15]的位移邊界條件及考慮本文模擬的實(shí)際情況，采用的有限元模型中土體底部邊界和土體四周邊界施加全約束.為方便有限元模擬的進(jìn)行，采用土體自重應(yīng)力來近似考慮實(shí)際土體的初始應(yīng)力.而ANSYS軟件中設(shè)置的應(yīng)力文件存入和讀取為初始地應(yīng)力的模擬提供了便利之處.

4 利用樁底沉渣的單樁數(shù)值模擬

4.1 工程實(shí)例模擬分析

某辦公樓為多高層框架結(jié)構(gòu)，樁基礎(chǔ)采用樁型為鉆孔灌注樁，設(shè)計(jì)樁徑為600 mm，樁身混凝土等級(jí)C20，單樁極限承載力為900 kN，樁長(zhǎng)27 m.樁基施工結(jié)束后按規(guī)范要求進(jìn)行單樁靜載荷試驗(yàn)，發(fā)現(xiàn)13號(hào)樁單樁極限承載力達(dá)不到設(shè)計(jì)要求，屬不合格樁.根據(jù)施工現(xiàn)場(chǎng)和地質(zhì)資料推斷樁下部有較厚的沉渣存在，取芯檢驗(yàn)結(jié)果發(fā)現(xiàn)樁底存在大約300 mm的軟質(zhì)沉渣，導(dǎo)致單樁承載力降低.其樁及地質(zhì)環(huán)境具體參數(shù)見表1.

表1 樁及土體的設(shè)置參數(shù)

對(duì)此工程實(shí)例進(jìn)行數(shù)值模擬，得到樁底有沉渣和無沉渣時(shí)的荷載-沉降曲線，并與實(shí)測(cè)的荷載-沉降數(shù)據(jù)對(duì)比，如圖3所示.

由圖3可見數(shù)值模擬很好地反映了工程實(shí)際情況，并發(fā)現(xiàn)樁底存在的沉渣對(duì)樁的承載力有較大的降低影響和對(duì)沉降有較大的增大影響，樁底沉渣的存在使單樁極限承載力降低了20%左右.

4.2 沉渣厚度對(duì)單樁承載性狀的影響

采用4.1節(jié)的樁土模型參數(shù)，但沉渣的彈性模量改用7.5 MPa，沉渣厚度分別為0 mm、50 mm、100 mm、200 mm、300 mm、400 mm、600 mm和800 mm，得到荷載-沉降曲線如圖4所示.從圖4可以看出沉渣厚度變化對(duì)單樁荷載-沉降曲線有著較大的影響.在相同樁頂荷載下，樁頂沉降隨著沉渣厚度的變大而增大.

圖5 不同沉渣厚度下的極限承載力損失率
Fig.5 Loss rate of ultimate capacity with different thickness

圖5為不同沉渣厚度下極限承載力的損失率,即各沉渣厚度下相對(duì)于沉渣厚度為0 mm時(shí)的極限承載力損失率.隨著樁底沉渣厚度的增加，單樁極限承載力逐漸衰減，但沉渣達(dá)到一定厚度后，極限承載力損失率趨于穩(wěn)定.當(dāng)沉渣厚度從100 mm增大到200 mm時(shí)，極限承載力損失率突增到8.2%，損失率接近100 mm沉渣時(shí)的2倍，符合《建筑樁基技術(shù)規(guī)范》(JGJ 94—2008)中對(duì)摩擦樁沉渣厚度不應(yīng)大于100 mm的要求.沉渣厚度繼續(xù)增大到300 mm，極限承載衰減速率變緩，其極限承載力損失率為9.5%，在10%以內(nèi).當(dāng)沉渣厚度從300 mm增加到400 mm，極限承載力減少的趨勢(shì)加大，其極限承載力損失率為11.3%.但沉渣厚度超過400 mm后，極限承載力損失率增大趨勢(shì)明顯放緩，趨于收斂，其極限承載力損失率在11.3%～12%之間，這表明沉渣影響存在臨界厚度，在300～400 mm之間.

4.3 沉渣物理力學(xué)性質(zhì)的影響

同樣采用4.1節(jié)的樁土模型，主要考慮樁底沉渣的彈性模量Ec、泊松比νc、粘聚力cc和內(nèi)摩擦角φc因素的各單因素影響性分析，避免各參數(shù)交叉分析時(shí)產(chǎn)生的不確定性.樁底沉渣厚度取為300mm.樁底沉渣彈性模量的取值為1.5MPa(0.06)、5MPa(0.2)、10MPa(0.4)、15MPa(0.6)、20MPa(0.8)、25MPa(1.0)，沉渣泊松比取值為0.2(0.8)、0.25(1.0)、0.3(1.2)，沉渣粘聚力取值為0.3(0.1)、1(0.33)、2(0.67)、3(1.0)，沉渣內(nèi)摩擦角取值為5°(0.15)、10°(0.31)、15°(0.46)、30°(1.0)，其中括號(hào)里的數(shù)據(jù)表示樁底沉渣參數(shù)值與相應(yīng)樁端下臥土層參數(shù)值的比值，即Ec/Ed、νc/νd、cc/cd、tanφc/tanφd，其中Ed、νd、cd和φd表示樁端下臥土層的彈性模量、泊松比、粘聚力和內(nèi)摩擦角.

圖6 不同沉渣彈性模量下的荷載-沉降曲線
Fig.6 Load-settlement curve with different elasticity modulus

圖6為樁底沉渣不同彈性模量下的荷載-沉降曲線.從圖中可以看出沉渣彈性模量對(duì)樁的承載力和沉降都存在明顯的影響.隨著樁底沉渣彈性模量的減小，相同荷載作用下產(chǎn)生的沉降量增大.當(dāng)沉渣的彈性模量減小到一定值后，沉降的增加幅度越來越明顯，沉渣彈性模量的影響越為顯著.當(dāng)沉渣的彈性模量從與樁端下臥層土的彈性模量一致的25MPa減小至10MPa，即從1.0Ed降到0.4Ed，荷載-沉降曲線基本無變化，但當(dāng)沉渣的彈性模量繼續(xù)減小至小于10MPa后，即小于0.4Ed，荷載-沉降曲線開始變化明顯.

樁底沉渣的彈性模量不同時(shí)相應(yīng)的極限承載力見圖7.圖中曲線反映了沉渣的彈性模量對(duì)樁極限承載力的影響，可見沉渣彈性模量在10～25MPa時(shí)，即0.4Ed～1.0Ed時(shí)，樁的極限承載力變化不大，最大的衰減幅值僅為2.2%.而沉渣彈性模量降低到5MPa時(shí)，即為0.2Ed時(shí)，相當(dāng)于沉渣彈性模量為25MPa時(shí)的樁的極限承載力，降低幅度達(dá)7%.當(dāng)沉渣的彈性模量繼續(xù)降低到1.5MPa時(shí)，承載力降低率甚至達(dá)到15.6%.可見樁底沉渣的彈性模量減小到0.2Ed～0.4Ed后，沉渣的彈性模量對(duì)樁的極限承載力有較大影響，影響性類似于指數(shù)衰減.

圖8為不同樁底沉渣泊松比下的荷載-沉降曲線.可見在沉渣的泊松比為0.8νd～1.2νd時(shí)，沉渣泊松比的變化對(duì)荷載-沉降曲線的影響很小.顯然其極限承載力很接近，相同樁頂荷載下的沉降相差在1mm之內(nèi).

不同沉渣粘聚力下的荷載-沉降曲線如圖9所示.由圖可見樁底沉渣的粘聚力對(duì)荷載-沉降曲線的影響并不明顯，不同粘聚力取值下荷載-沉降曲線基本重合，說明沉渣粘聚力的變化對(duì)單樁的承載力和沉降基本沒有影響.

樁底沉渣不同內(nèi)摩擦角下的荷載-沉降曲線如圖10所示.樁底沉渣的內(nèi)摩擦角的變化對(duì)荷載-沉降曲線有較明顯的影響，類似于沉渣的彈性模量對(duì)荷載-沉降曲線的影響.隨著內(nèi)摩擦角的減小，單樁的極限承載力降低且沉降增加.在內(nèi)摩擦角為5°，即tanφc/tanφd為0.15時(shí)的荷載-沉降曲線已顯然和其他三條曲線分離.這是因?yàn)閺椥阅Ａ亢蛢?nèi)摩擦角會(huì)影響到土體的抗剪強(qiáng)度.

4.4 樁長(zhǎng)徑比的影響

假定樁周土為單一的均質(zhì)土，其樁、樁周土和沉渣的各指標(biāo)參數(shù)如表2所示，沉渣的厚度固定為300 mm，樁的長(zhǎng)徑比取值為10(9/0.9)、20(18/0.9)、30(27/0.9)、45(27/0.6)、60(36/0.6)和75(45/0.6)，其中括號(hào)中的數(shù)據(jù)表示樁長(zhǎng)與樁徑之比.

表2 各材料參數(shù)

在樁的不同長(zhǎng)徑比情況下，豎向荷載作用下的荷載-沉降曲線如圖11所示.從圖中可見樁的長(zhǎng)徑比變化會(huì)影響沉渣對(duì)樁荷載-沉降曲線的作用.不同長(zhǎng)徑比時(shí)，樁底沉渣的存在會(huì)使單樁的沉降加大.但不難看出，隨著樁長(zhǎng)徑比的增大，有無樁底沉渣下的荷載-沉降曲線趨于緊密，兩者之間的變化幅度越來越小，尤其在樁長(zhǎng)徑比為75時(shí)的荷載-沉降曲線在陡降段前近于重合.說明隨著樁長(zhǎng)徑比變大，樁底沉渣對(duì)單樁荷載-沉降曲線的影響越來越小.

圖11 不同長(zhǎng)徑比下的荷載-沉降曲線
Fig.11 Load-settlement curve with different length diameter ratio

不同長(zhǎng)徑比下的單樁極限承載力損失率如圖12所示.由圖可見當(dāng)樁徑一致時(shí)，樁長(zhǎng)增加，長(zhǎng)徑比變大，單樁極限承載力損失率減小.當(dāng)長(zhǎng)徑比較小時(shí)，如長(zhǎng)徑比為10時(shí)，極限承載力損失率超過15%，表明此時(shí)沉渣對(duì)單樁的承載力有較大的影響.但當(dāng)長(zhǎng)徑比大于30時(shí)，極限承載力損失率在5%以內(nèi).長(zhǎng)徑比繼續(xù)增大到60以上后，極限承載力損失率可控制在2%以內(nèi)，表明此時(shí)沉渣對(duì)單樁承載力的影響可以忽略.

圖12 不同長(zhǎng)徑比下的極限承載力損失率
Fig.12 Loss rate of ultimate capacity with different length diameter ratio

5 結(jié) 語

本文對(duì)樁底沉渣對(duì)單樁承載性狀的影響展開了研究分析，并得到以下結(jié)論：

(1) 樁底沉渣厚度的影響存在一個(gè)臨界值，在300～400 mm之間.沉渣厚度大于100 mm后，單樁極限承載力損失率會(huì)有一個(gè)突增.工程中應(yīng)嚴(yán)格按照規(guī)范要求控制沉渣厚度.

(2) 沉渣彈性模量的減小導(dǎo)致單樁極限承載力類似于指數(shù)形式衰減，同樣造成沉降增大.隨沉渣內(nèi)摩擦角減小，單樁極限承載力減小，樁頂沉降變大.但沉渣的粘聚力和內(nèi)摩擦角對(duì)單樁承載性狀的影響并不明顯.因此在工程中進(jìn)行樁底沉渣檢測(cè)時(shí)，可主要考慮沉渣彈性模量和內(nèi)摩擦角的大小對(duì)單樁承載性狀的影響，而不必太注重沉渣泊松比和粘聚力參數(shù).

(3) 隨著樁長(zhǎng)徑比變大，樁底沉渣對(duì)單樁承載力和樁頂沉降的影響變小.樁長(zhǎng)徑比大于30后，單樁極限承載力損失率已很小，所以工程中灌注樁的長(zhǎng)徑比大于30時(shí)，可以適當(dāng)合理忽略沉渣對(duì)單樁承載性狀產(chǎn)生的微弱不良影響.

[1] 樓曉明, 陳強(qiáng)華, 俞有煒, 等. 鉆孔灌注樁承載力異常現(xiàn)象分析[J]. 巖土工程學(xué)報(bào), 2001, 23(5): 547-551.

[2] WONG Y W. Behaviors of large diameter bored pilegroups with defects[D].Hong Kong:Department of Civil Engineering, The Hong Kong University of Science and Technology, 2004.

[3] ISKANDER M, ROY D, KELLEY S, et al. Drilled shaft defects: detection, and effects on capacity in varved clay[J]. Journal of Geotechnical and Geoenviron-mental Engineering,2003, 129(12):1 128-1 137.

[4] 周紅波. 樁側(cè)泥皮和樁底沉渣對(duì)鉆孔樁承載力影響的數(shù)值模擬[J]. 巖土力學(xué),2007(5):956-960.

[5] 袁振,陳錦劍,王建華. 樁端條件對(duì)灌注樁荷載傳遞特性的影響[J]. 巖土力學(xué),2006(8):1 398-1 402.

[6] 胡鋮波,梅嶺,梅國(guó)雄，等. 樁土模型中土體邊界選取的有限元分析[J]. 建筑科學(xué),2009(9):18-20,29.

[7] 紀(jì)淑鵬. 彈塑性土體中豎向受荷單樁數(shù)值分析研究[J]. 工業(yè)建筑,2006,S1:755-757，760.

[8] 戚科駿,張?jiān)栖?王旭東，等. 考慮大變形的靜力壓樁數(shù)值模擬[J]. 巖土力學(xué),2005,S1:129-132.

[9] DESAI C S. Numerical design-analysis for piles in sands[J]. Journal of Geotechnical Engineering Division, 1974,100(6):613-633.

[10] RANDOLPH M F, POULOS H G. Estimating of the flexibility of offshore pile group[C]. Proc 2nd Int Conf Numerical Method in Offshore Piling Austin,1982:313-328.

[11] 錢家歡，殷宗澤. 土工原理與基礎(chǔ)[M].2版.北京：中國(guó)水利水電出版社，1998.

[12] HJIAJ M,FORTIN J,DE SAXCé G. A complete stress update algorithm for the non-associated Drucker-Prager model including treatment of the apex[J]. International Journal of Engineering Science,2003,4 110：1 109-1 143.

[13] 周鳳璽,李世榮. 廣義Drucker-Prager強(qiáng)度準(zhǔn)則[J]. 巖土力學(xué),2008(3):747-751.

[14] 劉自由,林杭,江學(xué)良. 樁頂豎向荷載作用下樁土響應(yīng)的數(shù)值分析[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011(2):508-513.

[15] 董志高,付婷婷,洪曉珍,等. 灌注樁樁側(cè)泥皮性質(zhì)對(duì)樁承載性狀影響研究[J]. 建筑科學(xué),2014(9):16-21.

責(zé)任編輯：羅 聯(lián)

Influence of Bottom Sediment on the Bearing Capacity of Bored Cast-in-Place Pile

FUWen-qiao,GUOXiao-gang*,CHENYu-qin,LIPing-jun

(College of Civil Engineering and Mechanics，Xiangtan University，Xiangtan 411105 China)

Many scholars pay attention to harmful effects of sediment which always exist at bottom of pile as bored cast-in-place pile is widely used in any engineering field. They also study effects of bottom sediment by using different methods. This paper rely on finite element software which is workable and comprehensive on dealing with geotechnical problems to research the influence of bottom sediment on the bearing capacity of bored cast-in-place pile. It also analysis the effects from three sides including the thickness of sediment, the physical properties of sediment and length diameter ratio of pile. Thickness, elasticity modulus，internal friction angle of sediment and length diameter ratio of pile have obvious effects, but cohesive force and poisson ratio of sediment have no obvious effects. Effect of thickness also has a critical value.

sediment；cast-in-place pile；numerical simulation；bearing capacity

2015-01-14

郭小剛(1960— )，男，湖南 湘潭人，教授.E-mail:teenhero@163.com

TU473.1

A

1000-5900(2015)02-0053-08