岳晨雨，馬震岳，王 偉，霍學(xué)平，王明輝

(大連理工大學(xué) 建設(shè)工程學(xué)部 水利水電研究所，遼寧 大連 116024)

接觸單元在樁—土接觸模擬中的應(yīng)用研究

岳晨雨，馬震岳*，王 偉，霍學(xué)平，王明輝

(大連理工大學(xué) 建設(shè)工程學(xué)部 水利水電研究所，遼寧 大連 116024)

針對(duì)坐落于軟巖地基上的水電站廠房地基摩擦樁加固措施，采用ANSYS中的面—面接觸單元模擬樁—土之間的接觸特性，并結(jié)合一中等直徑灌注樁，分析樁土模型的荷載—樁頂沉降曲線、接觸黏聚力—樁側(cè)/端摩阻力和樁徑—樁側(cè)/端摩阻力變化關(guān)系，總結(jié)荷載和參數(shù)作用與影響規(guī)律，為相似工程的樁土接觸數(shù)值模擬處理提供參考。

樁—土接觸；樁基礎(chǔ)；沉降曲線；接觸摩阻力

0 引 言

近年來，水電站建設(shè)的主戰(zhàn)場(chǎng)已逐漸轉(zhuǎn)移到水能豐富但地質(zhì)條件非常復(fù)雜的青藏高原地區(qū)，此處覆蓋層大都是軟土性質(zhì)的砂礫石土[1]?，F(xiàn)在面臨的最關(guān)鍵問題并不是水電站廠房結(jié)構(gòu)設(shè)計(jì)等技術(shù)問題，而是如何解決地基沉降尤其不均勻沉降等地基處理問題。我國(guó)工程建設(shè)中廣泛采用樁基礎(chǔ)，以提高軟土地基的承載力，但是人們并不能合理預(yù)測(cè)地基沉降量以及荷載傳遞規(guī)律，最終導(dǎo)致水電站主體結(jié)構(gòu)沉降量較大及其與從屬結(jié)構(gòu)之間的差異沉降量突出，從而對(duì)水電站廠房穩(wěn)定運(yùn)行產(chǎn)生不利影響。以往研究?jī)H給出樁頂荷載—沉降變形曲線、樁頂荷載—軸力變化曲線[2-3]，不能全面地說明樁土接觸之間的問題(接觸狀態(tài)，接觸摩擦，接觸壓力)，本文結(jié)合ANSYS軟件中面—面接觸單元，深入研究了樁頂荷載與接觸面摩阻力、樁端壓力的關(guān)系，為以后類似工程的技術(shù)設(shè)計(jì)提供參考。

1 模型選擇

1.1 土體與樁本構(gòu)模型選擇

ANSYS模型中的土體采用了Drucker-Prager(DP)材料，其材料特性包括黏聚力c、內(nèi)摩擦角φ和膨脹角φ。DP材料采用的是Drucker-Prager(德魯克- 普拉格) 屈服準(zhǔn)則。此屈服準(zhǔn)則是對(duì)Mohr-Coulumb(莫爾- 庫(kù)倫) 屈服準(zhǔn)則的近似，以此來修正VonMises(米賽斯) 屈服準(zhǔn)則。流動(dòng)準(zhǔn)則既可以使用相關(guān)流動(dòng)準(zhǔn)則，也可以使用不相關(guān)流動(dòng)準(zhǔn)則，材料是非強(qiáng)化材料，所以其屈服面并不隨材料的逐漸屈服而改變，因此沒有采用強(qiáng)化準(zhǔn)則[4]。其屈服強(qiáng)度隨著側(cè)限壓力的增加其塑性行為被假定為理想彈塑性。樁的本構(gòu)模型采用線彈性模型，主要是因?yàn)闃兜膹椥阅Ａ肯鄬?duì)土的壓縮模量較大，當(dāng)加載到極限值時(shí)，樁仍處于彈性階段。

1.2 樁—土體系模型

樁土接觸屬于邊界非線性問題，有限元建模時(shí)在樁與土之間引入面—面接觸單元模擬樁與土之間的剪力傳遞和相對(duì)滑移[5]。本文選擇剛體—柔體的接觸單元來模擬樁土交界面的非線性。在接觸問題中，一般設(shè)定較硬的面為目標(biāo)面，較軟的面為接觸面，所以目標(biāo)面在樁側(cè)面，接觸面在土側(cè)面，分別用Targe170、Conta174來模擬，該單元的特點(diǎn)是能夠進(jìn)行大變形計(jì)算，如較大尺寸的張開和滑移。

至今，普遍應(yīng)用的面—面接觸模型中接觸摩擦力公式[6]為：

CS=CP×MμCP=CT×FKN

FKN=Mμ×PRES/SLTO

(1)

式中CS為接觸面間的摩擦力；CP為接觸壓力；CT為接觸面間的接觸滲透值；Mμ為接觸面間的摩擦系數(shù)，本文取0.3；FKN為接觸單元法向剛度；PRES為接觸單元接觸壓力；SLTO為接觸單元最大滑移因子。

傳統(tǒng)的解釋是：摩擦系數(shù)Mμ和法向剛度FKN控制了接觸的受力行為，它們的定義正確與否是決定接觸結(jié)果是否合理的關(guān)鍵。

為了使結(jié)構(gòu)能夠安全運(yùn)行，在建設(shè)水電站之前，首先需要確定廠址所在地基的最大承載力，按建筑樁基技術(shù)規(guī)范JPG94-2008提供的單樁豎向極限承載力標(biāo)準(zhǔn)值[7-8]，可按下式估算：

(2)

式中qs i κ為樁側(cè)第i層樁側(cè)第i層土的極限側(cè)阻力標(biāo)準(zhǔn)值；qp κ為極限端阻力標(biāo)準(zhǔn)值。

當(dāng)模型利用公式(1)計(jì)算出的接觸面間的摩阻力大于式(2)的值時(shí)，認(rèn)為結(jié)構(gòu)破壞。另外一種判斷方法是，樁基受荷載沉降，當(dāng)沉降超過規(guī)范規(guī)定的極限值，認(rèn)為結(jié)構(gòu)也達(dá)到極限狀態(tài)。

2 單樁與土體ANSYS計(jì)算模型

2.1 土層和樁的材料特性

本文引用王先軍等[6]的實(shí)驗(yàn)?zāi)Ｐ停肁NSYS重新模擬，創(chuàng)建幾何模型，劃分網(wǎng)格，設(shè)置接觸參數(shù)進(jìn)行計(jì)算，并與其實(shí)驗(yàn)結(jié)果和模擬結(jié)果比較，結(jié)果證明本文模擬結(jié)果可靠，同時(shí)繪出荷載—沉降曲線。模型各參數(shù)見表1、表2。

表1 樁參數(shù)表

表2 土的力學(xué)參數(shù)表

Potyondy和Acer等的研究表明，樁土界面之間的摩擦角φ是影響摩擦樁承載性能的關(guān)鍵因素，對(duì)于黏土取0.6～0.7(樁周土體的有效內(nèi)摩擦角)是比較合適的?，F(xiàn)場(chǎng)澆筑的樁如灌注樁，接觸面間的摩擦性能比較好，接觸面的c和φ值可以取與樁相鄰?fù)翆拥腸和φ的0.8倍左右，可以根據(jù)地質(zhì)條件作相應(yīng)的調(diào)整，本文取0.8[9]。

修正之后，接觸面的力學(xué)參數(shù)見表3。

表3 修正之后接觸面的力學(xué)參數(shù)表

2.2 樁—土幾何模型

由于模擬重點(diǎn)是樁—土接觸模擬，所以樁和土體均采用Solid45 單元。靜載作用下，樁和土體受力是空間軸對(duì)稱的，可利用其對(duì)稱性取1/4模型。根據(jù)以往經(jīng)驗(yàn)，沿水平向距樁徑中心5 m以外以及樁端3 m 以下的土體的變形很小，所以土體水平方向取5 m，豎向取11 m。土體圓柱體的底面及側(cè)面約束全部自由度，兩1/4 界面采用面對(duì)稱約束，頂面完全自由。

樁—土體系有限元網(wǎng)格模型見圖1，接觸面模型見圖2 。

圖1 樁土有限元模型縱剖面和橫剖面Fig.1 Longitudinal and transverse section of the FEM model of pile and soil

圖2 樁土接觸面Fig.2 Interface between pile and soil

圖2中c為接觸面之接的黏聚力，Pa；fmax為樁土接觸面之間的最大靜摩擦應(yīng)力，Pa。

3 計(jì)算結(jié)果分析

利用ANSYS軟件，在樁頂分級(jí)加靜載300、500、700、900、1 100、1 310 kN，計(jì)算得樁頂沉降模擬值和實(shí)驗(yàn)值見表4。

本文模擬值與實(shí)驗(yàn)數(shù)據(jù)、王先軍模型模擬值吻合較好。有一定差異，主要是接觸參數(shù)的取值設(shè)定和ANSYS計(jì)算的誤差等因素影響，但誤差滿足工程要求。因此本文應(yīng)用面—面接觸單元模擬樁土接觸的方法比較合理。王先軍主要利用ANSYS軟件模擬其所作的試驗(yàn)?zāi)Ｐ?，僅研究了樁軸力的變化規(guī)律。本文主要在王先軍的實(shí)驗(yàn)?zāi)Ｐ偷幕A(chǔ)上研究樁土接觸側(cè)/端黏聚力、樁徑和樁長(zhǎng)對(duì)接觸面最大摩擦應(yīng)力、接觸面總摩擦力、接觸端壓力以及樁土沉降的影響。

表4 樁頂在各級(jí)荷載作用下沉降模擬值和實(shí)驗(yàn)值對(duì)比

Table 4 Simulation and experiment value contrast of pile top’s settlement at all levels load

荷載值/kN實(shí)驗(yàn)值/mm王先軍模型模擬值/mm本文模擬值/mm與實(shí)驗(yàn)值相對(duì)誤差/%000003001.01.11.110 5002.12.01.99.57002.93.43.210.39005.05.35.48.011008.98.69.12.2131010.811.211.23.7

3.1 樁—土接觸側(cè)黏聚力的影響

將上述模型簡(jiǎn)化，幾何模型不變,樁周土均設(shè)置成1號(hào)土層的材料參數(shù)，樁的材料參數(shù)不變,并根據(jù)建筑樁基技術(shù)規(guī)范JPG94-2008，設(shè)置樁側(cè)樁土接觸面間最大靜摩擦力55 000 Pa,樁端樁土接觸面間最大靜摩擦力850 000 Pa。樁頂荷載，主要由樁側(cè)接觸面摩擦力和樁端壓力承擔(dān)，即由式(1)求得。當(dāng)樁側(cè)接觸面摩擦力達(dá)到最大靜摩擦力后，繼續(xù)加載，多余的力就沿樁軸向下傳遞，最終由樁端土層承擔(dān)[1-10]，這是導(dǎo)致樁沉降的主要原因。樁頂施加豎向荷載1 000 kPa,取樁側(cè)樁土接觸黏聚力分別為15.0、20.0、25.6、30.0、40.0、50.0、60.0、65.0 kPa，樁端樁土接觸黏聚力為25.6 kPa。計(jì)算得樁側(cè)摩阻力、樁端壓力隨樁側(cè)樁土接觸黏聚力的變化見表5。

由表5可見，樁側(cè)樁—土接觸黏聚力對(duì)樁土接觸摩擦力和樁端接觸壓力均有一定影響，隨著樁側(cè)黏聚力增大，樁側(cè)摩阻應(yīng)力逐漸增大，最后增大到55 kPa(樁土接觸最大靜摩擦力)不再變化；樁側(cè)總摩擦力也隨著黏聚力增大而增大，而樁端壓力逐漸減小，但變化趨勢(shì)相對(duì)樁側(cè)摩擦力較小。說明樁側(cè)黏聚力對(duì)樁側(cè)摩阻力影響較大。

3.2 樁—土接觸端黏聚力的影響

分別取樁端樁土接觸黏聚力為20.0、30.0、40.0、50.0、60.0 kPa，樁側(cè)樁土接觸黏聚力為25.6 kPa,其他參數(shù)與上述模型一樣，計(jì)算得樁側(cè)摩阻力、樁端壓力隨樁端樁土接觸黏聚力的變化見表6。

表5 樁側(cè)不同樁土接觸黏聚力對(duì)樁側(cè)阻力和樁端壓力的影響

Table 5 Effect of the different contact cohesions between pile and soil in the side of the pile on friction against the side and pressure against the end

樁土接觸黏聚力/kPa樁側(cè)最大摩阻力/Pa樁側(cè)總摩擦力/N樁端總壓力/N15.019014551342633920.024380557821797325.633848593271836230.038814600771822540.049837615781804550.055000623881782960.055000623891782765.0550006238917819

表6 樁端不同樁土接觸黏聚力對(duì)樁側(cè)阻力和樁端壓力的影響

Table 6 Effect of the different contact cohesion between pile and soil in the end of the pile on friction against the side and pressure against the end

樁土接觸端黏聚力/kPa樁側(cè)最大摩阻應(yīng)力/Pa樁側(cè)總摩擦力/N樁端總壓力/N20.033843592631833230.033851593781834840.033857594711837150.030237558991865960.0302415600218704

由表6可見，隨著樁端樁—土接觸黏聚力增大，上述特征值略有浮動(dòng)，但影響不明顯。比較表5、表6可見，樁—土接觸側(cè)黏聚力對(duì)上述特征值的影響明顯比樁端黏聚力影響大。若要有效提高樁基承載力，應(yīng)通過工程措施適當(dāng)提高樁側(cè)接觸黏聚力。

3.3 樁徑的影響

分別取樁半徑300、400、500 mm(規(guī)范規(guī)定可選取的半徑)，在樁頂加豎向荷載200 kPa,樁側(cè)/端樁土接觸黏聚力為25.6 kPa，幾何模型和其他材料參數(shù)與上述模型一致,計(jì)算得樁側(cè)摩阻力、樁端壓力隨樁徑的變化見表7。

樁徑大小的選取由上部荷載的大小，地質(zhì)資料，樁身材料特性等綜合來決定。此模型覆蓋層較薄，設(shè)置為短樁。根據(jù)工程經(jīng)驗(yàn)及所加荷載，選取3種樁徑，比較其對(duì)樁側(cè)最大摩阻應(yīng)力、樁側(cè)總摩擦力、樁端總壓力的影響。由表7可見，樁徑增大，樁側(cè)最大摩阻應(yīng)力和樁側(cè)總摩擦力減小、樁端總壓力增大，整體最大沉降略為減小，但總的來說，對(duì)短樁的上述特征值影響不大。

表7 不同樁徑對(duì)樁側(cè)阻力和樁端壓力的影響

Table 7 Effect of the different pile diameters on friction against the side and pressure against the end

樁半徑/mm樁側(cè)最大摩阻應(yīng)力/Pa樁側(cè)總摩擦力/N樁端總壓力/N整體最大沉降/mm50026591152727233.1040028131182426503.6630038691226517393.87

3.4 樁長(zhǎng)的影響

分別取樁長(zhǎng)2.5、5.7、8.0 m，樁徑定為0.3 m,在樁頂加豎向荷載200 kPa,樁側(cè)/端樁土接觸黏聚力為25.6 kPa，幾何模型和其他材料參數(shù)與上述模型一致，計(jì)算得樁側(cè)摩阻力、樁端壓力隨樁長(zhǎng)的變化見表8。

表8 不同樁長(zhǎng)對(duì)樁側(cè)阻力和樁端壓力的影響

Table 8 Effect of the different pile lengths on friction against the side and pressure against the end

樁長(zhǎng)/m樁側(cè)最大摩阻應(yīng)力/Pa樁側(cè)總摩擦力/N樁端總壓力/N整體最大沉降/mm2.583648986487010.825.747751142225205.718.038691226517393.87

由表8可見，隨著樁長(zhǎng)增大，樁側(cè)最大摩阻應(yīng)力減小、樁側(cè)總摩擦力增大、樁端總壓力減小。由計(jì)算結(jié)果可知，接觸面間的接觸摩擦應(yīng)力隨樁長(zhǎng)增加而減小，樁側(cè)總摩擦力反而增大，這是因接觸面積增加，導(dǎo)致分配到各個(gè)節(jié)點(diǎn)的摩擦應(yīng)力減小。樁端壓力是引起樁下土層變化的主要原因，隨樁長(zhǎng)增加，樁端壓力減小，最終的沉降變形也會(huì)減小,樁與樁周土之間的相對(duì)滑移減小，發(fā)生滑移破壞所需的外荷載增大，最終提高樁土地基的承載能力。

比較表7和表8可見，樁長(zhǎng)相對(duì)于樁徑對(duì)上述特征值影響較大。

結(jié)合理論和數(shù)值分析，說明利用ANSYS軟件模擬樁土之間的接觸是合理的，能夠恰當(dāng)?shù)啬M樁土之間的接觸摩擦、樁端壓力與樁側(cè)/端樁土接觸黏聚力和樁徑、樁長(zhǎng)關(guān)系。通過對(duì)加固后地基進(jìn)行模擬，計(jì)算樁基承載力和總摩阻力，若地基承載力達(dá)不到標(biāo)準(zhǔn)，可通過改變樁長(zhǎng)等措施，提高樁基承載力，使地基上方結(jié)構(gòu)安全運(yùn)行。

4 結(jié) 論

1)本文按建筑樁基技術(shù)規(guī)范JPG94-2008，設(shè)置樁土接觸面間的最大靜摩擦力，結(jié)合不同結(jié)構(gòu)沉降變形極限值，從而確定該模型的最大承載力。最大承載力主要由樁土接觸的各參數(shù)決定，對(duì)實(shí)際工程進(jìn)行模擬時(shí)，要想得到較為合理的結(jié)果，必須獲得準(zhǔn)確的樁土物理參數(shù)。本模型模擬結(jié)果，為不同地基參數(shù)的選定以及對(duì)外載與樁土變形、樁側(cè)/端摩阻力等變化規(guī)律的認(rèn)識(shí)提供支持，同時(shí)為提高地基承載力提供了一種預(yù)估方法。

2)根據(jù)ANSYS模擬可見樁側(cè)接觸黏聚力對(duì)樁土接觸摩擦力以及樁端壓力影響較大，適當(dāng)改善樁側(cè)接觸黏聚力可提高樁基承載力。

3)短樁半徑的改變對(duì)地基承載力的影響相對(duì)較小。增加樁的長(zhǎng)度可提高樁土地基的承載力，此結(jié)論同樣適用于水電站廠房下的群樁地基。

[1]中國(guó)水力發(fā)電工程學(xué)會(huì)，中國(guó)水電工程顧問集團(tuán)公司，中國(guó)水利水電建設(shè)集團(tuán)公司.中國(guó)水力發(fā)電科學(xué)技術(shù)發(fā)展報(bào)告[R].北京:中國(guó)水力發(fā)電工程學(xué)會(huì),2012.

[2]肖宏彬,鐘輝虹,張亦靜,等.單樁荷載—沉降關(guān)系數(shù)值分析模擬方法[J].巖土力學(xué),2002,23(5):592-596.

[3]杜 蓉,張建友,隋麗麗,等.樁基承載力的ANSYS有限元分析[J].水利與建筑工程學(xué)報(bào),2010，(4):213-215.

[4]張學(xué)言，閆澍旺.巖土彈塑性力學(xué)基礎(chǔ)[M].天津:天津大學(xué)出版社,2004:78-81.

[5]宰金珉,宰金璋.高層建筑基礎(chǔ)分析與設(shè)計(jì)[M].北京:中國(guó)建筑工業(yè)出版社,1994.

[6]王先軍,周文宇,蔣 鑫.ANSYS在模擬樁土接觸中的應(yīng)用[J].森林工程,2006，(3):49-51.

[7]JGJ94-2008,建筑樁基技術(shù)規(guī)范[S].

[8]邱 敏,鹿 寧,鮑呈蒼,等.樁基在多布水電站廠房深覆蓋層地基中的應(yīng)用[J].西北水電,2014，(1):44-48.

[9]陳育民,徐鼎平.FLAC/FLAC3D基礎(chǔ)與工程實(shí)例[M].北京:中國(guó)水利水電出版社,2008:140-141.

[10]李作勤.摩擦樁的荷載傳遞及承載力的一些問題[J].巖土力學(xué),1990，(4):1-12.

Application of the interface elements between pile and soil

YUE Chen-Yu,MA Zhen-Yue*,WANG Wei,HUO Xue-Ping,WANG Ming-Hui

(Institute of Hydraulic and Hydropower，F(xiàn)aculty of Infrastructure Engineering，Dalian University of Technology,Dalian 116024,Liaoning,China)

Surface to surface interface element of ANSYS is cited to simulate the contact between pile and soil,and combines with a medium diameter filling pile.The load-settlement curve in the top of the pile,the contactcohesion-contact friction change rule and the pile diameter-contact friction change rule are analyzed.Summarizing the rule,the important information about pile-soil processing and the settlement forecast are provided for the same zone.

contact between pile and soil;pile foundation;settlement curve;contact friction

10.13524/j.2095-008x.2015.03.037

2015-05-18;

2015-05-31

國(guó)家自然科學(xué)基金資助項(xiàng)目(51379030)

岳晨雨(1989-)，女，河南商丘人，碩士研究生，研究方向：軟巖地基水電站廠房結(jié)構(gòu)穩(wěn)定分析，E-mail:YCYdlut@163.com；*通訊作者：馬震岳(1962-)，男，河南南陽(yáng)人，教授，博士研究生導(dǎo)師，研究方向：水輪發(fā)電機(jī)組動(dòng)力學(xué)、水工水電建筑物結(jié)構(gòu)分析，E-mail:dmzy@dlut.edu.cn。

TU473.12

A

2095-008X(2015)03-0001-05