吳鳴，鄧露

（汕頭大學(xué)土木工程系，廣東汕頭515063）

軟黏土中超長(zhǎng)樁豎向荷載作用下的三維數(shù)值模擬

吳鳴，鄧露

（汕頭大學(xué)土木工程系，廣東汕頭515063）

我國(guó)沿海地區(qū)，軟弱土層較厚，且分布較廣，超長(zhǎng)樁的使用已經(jīng)成為必然，然而有關(guān)超長(zhǎng)樁在豎向荷載作用下的荷載傳遞機(jī)理及設(shè)計(jì)方法的研究遠(yuǎn)遠(yuǎn)落后于工程實(shí)際應(yīng)用.本文基于ABAQUS大型非線性有限元軟件，采用適合于軟黏土的劍橋模型，且考慮土的正常固結(jié)效應(yīng)，建立了豎向荷載作用下，樁土共同作用的三維數(shù)值模型，考慮了樁土之間產(chǎn)生的滑移和開裂，在此基礎(chǔ)上分析了樁長(zhǎng)，樁徑以及樁身壓縮量對(duì)超長(zhǎng)樁承載性能的影響.

超長(zhǎng)樁；劍橋模型；承載能力；有限元模擬

0 引言

我國(guó)沿海地區(qū)地質(zhì)條件極其復(fù)雜，如廣東、上海等沿海深厚軟土達(dá)150～400 m，80 m以上超長(zhǎng)樁已成為軟土地區(qū)超高層建筑和橋梁等建筑物樁基常用的一種樁型.隨著樁長(zhǎng)與樁徑的增加，其承載性能及其影響因素更加復(fù)雜，現(xiàn)行規(guī)范中仍沿用普通樁的設(shè)計(jì)理論.大量的試驗(yàn)分析表明，基樁的側(cè)阻和端阻發(fā)揮不同步，如果按照統(tǒng)一參數(shù)計(jì)算缺乏確切的安全度概念，這一問(wèn)題在超長(zhǎng)樁的設(shè)計(jì)中尤為突出；而超長(zhǎng)樁樁側(cè)阻在樁身中部常存在強(qiáng)化效應(yīng)，又使其設(shè)計(jì)存在一定的承載潛力，這種一方面基樁已具有一定承載潛力，另一方面設(shè)計(jì)人員仍試圖以增大樁長(zhǎng)提高基樁豎向承載力的矛盾，均源自于工程設(shè)計(jì)人員缺乏對(duì)超長(zhǎng)樁承載性狀的深入了解，因而根據(jù)現(xiàn)行規(guī)范對(duì)超長(zhǎng)樁的設(shè)計(jì)是不完善的，也是不合理的.實(shí)際工程中經(jīng)常出現(xiàn)因樁基沉降過(guò)大等引起的工程事故，或樁基設(shè)計(jì)中越來(lái)越保守的趨勢(shì)和嚴(yán)重浪費(fèi)的現(xiàn)象，因而開展超長(zhǎng)樁的承載性能與設(shè)計(jì)理論研究具有重要的工程意義和較大的經(jīng)濟(jì)價(jià)值，也是工程界的迫切要求.

單樁理論分析方法主要有：彈性理論法、荷載傳遞法、剪切位移法及數(shù)值分析法等，其中在眾多數(shù)值方法中，有限單元法發(fā)展得較為成熟，其可同時(shí)考慮眾多影響因素，如土的非線性、固結(jié)時(shí)間效應(yīng)以及動(dòng)力效應(yīng)等，許多學(xué)者為獲得各工況下的樁基承載性能，借助于有限元方法來(lái)模擬現(xiàn)場(chǎng)試驗(yàn)，得到了很多有意義的結(jié)論[1-4].然而在豎向荷載下，超長(zhǎng)樁樁土之間往往出現(xiàn)滑移，甚至脫開，樁周土往往處于塑性變形狀態(tài)，傳統(tǒng)的數(shù)值計(jì)算難以滿足如此要求，本文選用適合于樁土共同作用分析的ABAQUS大型非線性有限元軟件，采用適合于軟黏土的劍橋模型，且考慮土的正常固結(jié)效應(yīng)，建立了豎向荷載作用下，樁土共同作用的三維數(shù)值模型，考慮了樁土之間產(chǎn)生的滑移和開裂，在此基礎(chǔ)上分析了樁長(zhǎng)，樁徑以及樁身壓縮量對(duì)超長(zhǎng)樁承載性能的影響，得出了一些有益的結(jié)論，對(duì)于指導(dǎo)工程設(shè)計(jì)具有一定的意義.

1 軟黏土的本構(gòu)模型

軟土具有靈敏度高，低強(qiáng)度，低透水性等特性，從而使地基易產(chǎn)生不均勻沉降.劍橋塑性模型與多孔介質(zhì)彈性模型共同使用，可以很好的描述黏土的特性，可以考慮土的固結(jié)，以及考慮到孔隙比隨土層的變化，以及黏土的透水性.

多孔介質(zhì)彈性模型是一種非線性的各向同性彈性模型，該模型體積應(yīng)力應(yīng)變關(guān)系：

e0是初始孔隙比，P0是初始平均應(yīng)力，pelt是彈性狀態(tài)的拉應(yīng)力極限值，Jel是彈性體積應(yīng)變，κ是對(duì)數(shù)體積模量.

剪切應(yīng)力應(yīng)變關(guān)系為：

G為彈性剪切模量，eel為彈性偏應(yīng)變.

劍橋塑性模型采用聯(lián)合流動(dòng)法則，即塑性勢(shì)與屈服函數(shù)為同一函數(shù)：

2 樁土接觸面的摩擦模型

樁身材料是剛性材料，而土層除了未風(fēng)化的巖石外基本上屬于柔性材料，由于在豎向荷載作用下，兩種材料變形不一致，在其接觸面上產(chǎn)生剪應(yīng)力.研究接觸面的變形，主要從兩個(gè)方面來(lái)研究：一是接觸面的相互作用力學(xué)模型，它主要包含接觸面的法向作用和切向作用；二是接觸面單元，它是有限元分析時(shí)用來(lái)模擬接觸面變形的特殊單元.接觸面的力學(xué)模型要適應(yīng)相應(yīng)的接觸面單元，而接觸面單元需要表達(dá)接觸面的變形，兩者是相互聯(lián)系的.

ABAQUS對(duì)于接觸面的模擬，采用主-從接觸算法，提供了點(diǎn)-面接觸對(duì)和面-面接觸對(duì).由于面-面離散方法需要分析整個(gè)從面的接觸行為，其計(jì)算代價(jià)比點(diǎn)-面離散方法要高，因此，樁-土接觸面設(shè)置點(diǎn)-面接觸對(duì)，將剛性材料的樁身界面設(shè)置為主控面（master surface），土界面設(shè)置為從屬面（slave surface），采用小滑動(dòng)（small sliding）跟蹤方法.ABAQUS提供了硬接觸與軟接觸兩大類的接觸面法向模型，樁與土的接觸面主要考慮其剪應(yīng)力，不傳遞法向應(yīng)力，因此模擬時(shí)采用硬接觸模型，但是ABAQUS中的接觸面是沒(méi)有厚度只有長(zhǎng)度的，摩擦模型采用的是庫(kù)倫模型，是理想彈塑性的.1968年，Goodman等人提出了無(wú)厚度接觸單元[5]，這與ABAQUS提供的單元類似，只是摩擦模型不同，通常將剪應(yīng)力與剪切變形之間的關(guān)系模擬為非線性彈性；在此基礎(chǔ)上，1984年，Desai等人提出了薄層四邊形單元[6]；但是有厚度接觸單元在數(shù)值模擬時(shí)相應(yīng)參數(shù)的設(shè)計(jì)理論并不成熟，一些學(xué)者在此基礎(chǔ)上做了進(jìn)一步研究[7-8].

因此，本文有關(guān)于樁土接觸面的模擬采用Goodman無(wú)厚度接觸單元.

若不考慮兩個(gè)方向上剪應(yīng)力的相互影響，接觸面的本構(gòu)關(guān)系為：

克拉夫和鄧肯認(rèn)為剪應(yīng)力和相對(duì)剪切位移之間符合雙曲線關(guān)系，則將ks1、ks2表示為：

式中K1、K2、Rf、n為非線性指標(biāo)；δ是接觸面的界面摩擦角；rw為水的容重；pa為大氣壓力.

3 有限元模型及計(jì)算參數(shù)

ABAQUS是比較成熟的大型非線性有限元分析軟件，擁有能夠真實(shí)反映土體性狀的本構(gòu)模型，還具有強(qiáng)大的接觸面處理功能，能夠很好的模擬土與樁之間的相互作用，并且提供了二次開發(fā)接口，用戶可以自定義材料特性，單元特性，以及接觸面模型等等.而用ABAQUS模擬樁土相互作用，最主要的就是能夠很好的模擬樁土的材料性能以及樁土接觸面的特點(diǎn).

本計(jì)算模型樁長(zhǎng)為40 m，樁徑為1 m（目前對(duì)于超長(zhǎng)樁的定義l/d≥40），樁周土取樁徑的10倍，取模型的1/4進(jìn)行計(jì)算，如圖1，樁采用線彈性模型，彈性模量E＝28 Gpa，泊松比ν＝0.2，土采用劍橋模型，相關(guān)計(jì)算參數(shù)如表1，樁土接觸面采用goodman單元，接觸面的界面摩擦角δ為10度.

表1 土的模型參數(shù)[9]

考慮土體的正常固結(jié)效應(yīng)，以及孔隙比沿深度線性分布等影響，得到豎向荷載作用下樁土共同作用的變形后的應(yīng)力模型圖（如圖2）、荷載-位移曲線圖（如圖3）、不同時(shí)刻樁側(cè)摩阻力圖（如圖4）以及樁端土體的孔壓場(chǎng)（圖5）.

從圖3可以看出，荷載與位移層曲線變化，這與理想的Q-S曲線比較接近，超長(zhǎng)單樁的承載力達(dá)到了極限值，在極限承載力之前，荷載與位移呈線性變化，而且位移隨著荷載的增加而緩慢增加，隨后，位移隨著荷載的增加速率加快，這就說(shuō)明此時(shí)，樁以較快的速度刺入樁端土中.

從圖4可以看出，計(jì)算的摩阻力最大值為125.4 kPa，不同時(shí)刻的樁側(cè)摩阻力都是在接近樁端的時(shí)候達(dá)到最大值，但是在未接近樁端時(shí)，摩阻力比較小，在樁深20 m之前呈緩慢增長(zhǎng)速度，在之后呈緩慢遞減趨勢(shì)，這主要是與樁身壓縮變形有關(guān)，樁身壓縮量是隨著荷載的增加而增大.這就說(shuō)明，超長(zhǎng)樁屬于摩擦型樁，樁側(cè)摩阻力是從樁頂開始，隨著荷載的增大，沿著樁身長(zhǎng)度逐漸變大，在接近樁底時(shí)，達(dá)到極限摩阻力.

由圖5及計(jì)算結(jié)果表明，樁端下的土體由于受到豎向壓力，產(chǎn)生了正的孔壓，最大值約為410 kPa；而在靠近樁端上部的土體中，由于下部土體受壓有向上滑動(dòng)的趨勢(shì)，因而產(chǎn)生負(fù)的孔壓，最大值約為-150 kPa.樁土接觸面的孔壓也得到了比較滿意的結(jié)果，在樁土之間沒(méi)有發(fā)生流體的流動(dòng)，與模擬不排水條件相吻合.

圖1 樁土有限元模型

圖2 變形后的應(yīng)力模型圖

圖3 荷載-沉降曲線

圖4 不同時(shí)刻的樁側(cè)摩阻力

4 影響因素分析

1）樁長(zhǎng)對(duì)單樁承載性能的影響

基于前面有限元模型及樁土的力學(xué)性能相同條件下，不同樁長(zhǎng)情況下的荷載-位移曲線分別如圖6（樁徑1.5 m）、圖7（樁徑2.0 m）.

由圖6-7可以看出，隨著樁長(zhǎng)的增大，單樁極限承載力越大.在施加同樣的位移（displacement=100 mm）條件下，樁長(zhǎng)越長(zhǎng)，樁頂承受的荷載越大，不過(guò)，這種增長(zhǎng)的趨勢(shì)隨著樁長(zhǎng)的增加而減少.

圖5 樁端土體的孔壓場(chǎng)

圖6 不同樁長(zhǎng)的荷載-位移曲線對(duì)比圖（樁徑1.5m）

圖7 不同樁長(zhǎng)的荷載-位移曲線對(duì)比圖（樁徑2 m）

2）樁徑對(duì)單樁承載性能的影響

基于前面有限元模型及樁土的力學(xué)性能相同條件下，不同樁徑情況下的荷載-位移曲線分別如圖8（樁長(zhǎng)60 m）、圖9（樁長(zhǎng)90 m）以及圖10（樁長(zhǎng)120 m）.

由圖8-10可以看出隨著樁徑的增大，單樁極限承載力增大.在施加同樣的位移（displacement=100 mm）條件下，樁徑越大，樁頂承受的荷載越大，不過(guò)這種趨勢(shì)隨著樁長(zhǎng)的增大而減少.

隨著樁長(zhǎng)及樁徑的增大，單樁承載力也增大，但是這種增長(zhǎng)的趨勢(shì)隨著樁長(zhǎng)及樁徑的增加而減少，所以，合理的選擇超長(zhǎng)樁的長(zhǎng)徑比很重要，在實(shí)際工程中，避免造成浪費(fèi).

圖8 不同樁徑的荷載-位移曲線對(duì)比圖（樁長(zhǎng)60m）

圖9 不同樁徑的荷載-位移曲線對(duì)比圖（樁長(zhǎng)90m）

3）樁身壓縮量的討論

超長(zhǎng)樁在豎向荷載作用下，其樁頂沉降主要是由樁本身的壓縮量以及樁側(cè)阻力引起的樁端土沉降兩部分構(gòu)成，因此對(duì)于樁身壓縮量的研究是有必要的.

基于前面有限元模型及樁土的力學(xué)性能相同條件下，不同樁徑情況下的荷載-樁身壓縮量曲線分別如圖11（樁徑1.5 m）、圖12（樁徑2.0 m）.

由圖11、圖12可以看出，樁身壓縮量基本上與樁頂荷載成正比，且隨著樁長(zhǎng)的增加而增大.樁長(zhǎng)越長(zhǎng)，樁身壓縮量的增大速率也呈現(xiàn)增漲趨勢(shì).

從上圖分析得到，超長(zhǎng)樁的樁身壓縮量占總沉降量的比例隨著樁長(zhǎng)的增加而增加，隨著樁徑的增大而減少，樁長(zhǎng)為60 m時(shí)，其樁身壓縮量在10 mm以內(nèi)，但是，樁長(zhǎng)120 m時(shí)，卻達(dá)到了55 mm左右.可見(jiàn)，樁長(zhǎng)越長(zhǎng)，對(duì)于樁身壓縮量的研究就更重要了，在實(shí)際工程中，超長(zhǎng)樁的設(shè)計(jì)就需要考慮樁身壓縮量.

樁本身的壓縮量是構(gòu)成超長(zhǎng)樁在豎向荷載作用下樁頂沉降的重要部分，所占比例隨著樁長(zhǎng)的增加而增大，尤其是樁長(zhǎng)達(dá)到120 m（樁徑1.5 m）時(shí)，比例達(dá)到了50%以上.

圖10 不同樁徑的荷載-位移曲線對(duì)比圖（樁長(zhǎng)120m）

圖11 不同樁長(zhǎng)的荷載-樁身壓縮量曲線對(duì)比圖（D=1.5m）

5 結(jié)論

（1）超長(zhǎng)樁屬于摩擦型樁，樁側(cè)摩阻力是從樁頂開始，隨著荷載的增大，沿著樁身長(zhǎng)度逐漸變大，在接近樁底時(shí)，達(dá)到極限摩阻力.

（2）隨著樁長(zhǎng)及樁徑的增大，單樁承載力也增大，但是這種增長(zhǎng)的趨勢(shì)隨著樁長(zhǎng)的增加而減少，所以，合理的選擇超長(zhǎng)樁的長(zhǎng)徑比很重要，在實(shí)際工程中，避免造成浪費(fèi).

（3）樁本身的壓縮量是構(gòu)成超長(zhǎng)樁在豎向荷載作用下樁頂沉降的重要部分，所占比例隨著樁長(zhǎng)的增加而增大，尤其是樁長(zhǎng)達(dá)到120 m（樁徑1.5 m）時(shí)，比例達(dá)到了50%以上.可見(jiàn)，樁長(zhǎng)越長(zhǎng)，對(duì)于樁身壓縮量的研究就更重要了，在此類型土中，超長(zhǎng)樁的設(shè)計(jì)就需要考慮樁身壓縮量.

圖12 不同樁長(zhǎng)的荷載-樁身壓縮量曲線對(duì)比圖（D=2.0m）

參考文獻(xiàn)

[1]李穎.超長(zhǎng)單樁豎向承載力影響因素有限元分析[J].濟(jì)南大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，25（3）：287-290.

[2]張忠苗，張乾青.樁端土強(qiáng)度對(duì)樁側(cè)阻力影響的研究[J].巖土工程學(xué)報(bào)，2010，32（增2）：59-63.

[3]徐亞利.黃土地基中大直徑超長(zhǎng)樁受力特性的非線性有限元分析[J].江蘇建材，2010（3）：28-30.

[4]何繼周，張東星.超長(zhǎng)樁與土共同作用性狀非線性有限元分析[J].施工技術(shù)，2009，38（增）：29-32.

[5]GOODMAN R E，TAYLOR R L，BREKKE T L.A model for the mechanics of jointed rock[J].J of the Soil Mechanics and Foundations，1968，94（3）：660-677.

[6]DESAI C S，ZAMAN MM，LIGHTHER J G，et al.Thin-layer element for interfaces and joints[J].Int J for Numerical Methods in Geomechanics，1984，8（1）：19-43.

[7]杜成斌，任青文.用于接觸面模擬的三維非線性接觸單元[J].東南大學(xué)學(xué)報(bào)（自然科學(xué)版），2001，31（4）：92-96.

[8]簡(jiǎn)政，劉用暖，李九紅.樁-土三維等厚度接觸單元模型及其工作性狀[J].水利學(xué)報(bào)，2002（11）：113-118.

[9]費(fèi)康，張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京：中國(guó)水利水電出版社，2009.

Three Dimensional Numerical Simulation of Super-Long Pile in Soft Clay Under Vertical Load

WU Ming,DENG Lu
（Department of Civil Engineering,Shantou University,Shantou 515063,Guangdong,China）

In coastal areas,weak soil layer is thick and has a wide distribution.The research of the bearing capacity of super-long pile in soft clay under the action of vertical load is important. Based on the large-scale nonlinear software of ABAQUS,the three-dimensional model of the interaction of pile and soil under the action of vertical load is established.The Cam-Clay is more suitable for the soil constitutive model and the normal soil consolidation effect is considered.The influence of pile length,pile diameter and pile body compression on the bearing capacity of super-long pile is analyzed.

super-long pile;cam-clay;bearing capacity;finite element simulation

TU473

A

1001-4217（2016）04-0074-07

2015-12-14

吳鳴（1975—），男（漢族），浙江省東陽(yáng)市人，副教授，研究方向：樁基設(shè)計(jì)理論，土-結(jié)構(gòu)動(dòng)力相互作用。E-mail：mwu@stu.edu.cn

汕頭市科技計(jì)劃項(xiàng)目（2011-153；2012-167）