范慶來+張炳磊+賈小偉

摘 要：采用基于耦合歐拉拉格朗日算法的三維大變形有限元方法對(duì)上硬下軟黏土地基方形樁靴的貫入過程進(jìn)行數(shù)值模擬，系統(tǒng)地探討了上下土層強(qiáng)度之比、上層土厚度、土體重度及靜止側(cè)壓力系數(shù)等指標(biāo)對(duì)于樁靴貫入阻力曲線的影響，揭示了穿透現(xiàn)象發(fā)生時(shí)的土體變形機(jī)理。計(jì)算結(jié)果表明：上下土層強(qiáng)度之比與上層土厚度對(duì)于樁靴的貫入阻力曲線具有很大影響。靜止側(cè)壓力系數(shù)的減小在一定程度上導(dǎo)致樁靴承載力系數(shù)的降低，而上下土層重度的差異導(dǎo)致樁靴的承載力系數(shù)不滿足前人建議的方程。在上下土層重度差異不大以及上層土厚度確定的情況下，關(guān)于樁靴寬度、上層土重度與強(qiáng)度的無量綱化參數(shù)小于臨界值時(shí)將會(huì)發(fā)生穿透現(xiàn)象。

關(guān)鍵詞：方形樁靴；貫入阻力；穿透機(jī)理；黏土；大變形分析

中圖分類號(hào)：TU476 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1674-4764（2018）01-0062-06

Analysis of penetration and punch-through failure of square spudcan in two-layered clays

Fan Qinglai1， 2，Zhang Binglei2，Jia Xiaowei2

（1. College of Civil and Transportation Engineering， Hohai University， Nanjing 210098， P.R.China；2. Key Laboratory of Geotechnical Engineering， Ludong University， Yantai 264025， Shandong，P.R.China）

Abstract：Based on three-dimensional large deformation finite element method utilizing coupled Eulerian-Lagrangian algorithm， the penetration responses of square spudcan in strong over soft clays were numerically simulated. The effects on the penetration resistance curve were studied of the strength ratio of the bottom clay layer to the top clay layer， the relative thickness of the top layer， the unit weight of soil， and the coefficient of earth pressure at rest. The soil deformation mechanism of punch-through failure during penetration was investigated. The results show that the strength ratio of the bottom clay layer to the top clay layer and the relative thickness of the top layer significantly affects the penetration resistance of square spudcan footings. The peak bearing capacity factor of square spudcan decreases somewhat with the decrease of the coefficient of earth pressure at rest. The penetration resistance curve does not obey the available equation proposed by other authors due to the possible difference of unit weight of the top layer and the bottom layer. The square spudcan footings experience punch-through failure when the proposed dimensional parameter in terms of footing width， the unit weight and undrained strength of upper soil is less than a critical value in the condition that the thickness of upper clay layer is defined and the difference of unit weight of the layers is not appreciable.

Keywords：square spudcan； penetration resistance； punch-through failure； clay； large deformation analyses

自升式海上鉆井平臺(tái)一般由上部平臺(tái)、升降裝置和樁腿等部分組成。這種平臺(tái)具有經(jīng)濟(jì)性和機(jī)動(dòng)性等優(yōu)點(diǎn)，已普遍應(yīng)用于近海油氣資源的勘探開采。自升式鉆井平臺(tái)的樁腿底端通常都設(shè)有倒圓錐形、方形及其它形式的樁靴。而海洋地基地質(zhì)條件比較復(fù)雜，很多地區(qū)的海床是由一層較硬的黏土層覆蓋在深厚軟黏土上。在這種上硬下軟地基條件下，樁靴很容易發(fā)生穿透破壞而導(dǎo)致上部平臺(tái)傾覆。Gourvenec等[1]、Merifield等[2]分別采用常規(guī)的小變形有限元方法和極限分析方法，提出了上硬下軟雙層黏土地基承載力系數(shù)，但是，這些研究主要針對(duì)建造在土層表面的淺基礎(chǔ)。張愛霞[3]根據(jù)假設(shè)的破壞模式，基于Meyerhof-Hanna沖剪理論，建立了一種雙層地基中方形樁靴承載力計(jì)算公式。Craig等[4]采用離心試驗(yàn)分析了樁靴在均勻黏土、砂土及上砂下黏雙層地基中的連續(xù)貫入過程，同時(shí)，建議了承載力估計(jì)公式。Lee等[5]隨后進(jìn)行了比較系統(tǒng)的離心試驗(yàn)，并給出了不同于現(xiàn)行SNAME規(guī)范的貫入阻力設(shè)計(jì)曲線。在數(shù)值模擬方面，由于樁靴基礎(chǔ)貫入導(dǎo)致土體產(chǎn)生大變形，傳統(tǒng)的基于總體拉格朗日（Total Lagrangian）列式或修正拉格朗日（Updated Lagrangian）列式的小變形有限元方法會(huì)產(chǎn)生網(wǎng)格畸變而導(dǎo)致計(jì)算中斷，無法跟蹤基礎(chǔ)連續(xù)貫入的整個(gè)過程。部分學(xué)者假設(shè)樁靴基礎(chǔ)預(yù)埋于地基中不同深度，仍然基于小變形有限元方法進(jìn)行研究，這雖然從一定程度上能夠反映出基礎(chǔ)貫入過程中承載力連續(xù)變化趨勢(shì)，但是，無法準(zhǔn)確預(yù)測(cè)穿透破壞時(shí)的入土深度。為了解決這個(gè)問題，Yu等[6]、Liu等[7]基于RITSS（remeshing and interpolation technique with small strain）大變形有限元方法，研究了方形樁靴貫入上硬下軟雙層黏土的過程，結(jié)果表明，在某些情況下，基礎(chǔ)連續(xù)貫入過程中的承載力隨著貫入深度變化，曲線同小變形計(jì)算得到的情況有較大差別，采用大變形有限元分析方法模擬基礎(chǔ)的連續(xù)貫入過程更加符合實(shí)際。最近，Tho等[8]、Qiu等[9]、任憲剛等[10]、鄭靜等[11]利用耦合歐拉拉格朗日（Coupled Eulerian-Lagrangian，簡稱CEL）有限元法對(duì)于自升式平臺(tái)樁靴貫入過程進(jìn)行了數(shù)值分析，在研究過程中采用拉格朗日單元模擬基礎(chǔ)，用歐拉網(wǎng)格代表土體所在區(qū)域，在土體發(fā)生大變形過程中，歐拉網(wǎng)格固定，因此，不產(chǎn)生單元扭曲畸變問題，模擬結(jié)果與相關(guān)試驗(yàn)比較吻合。endprint

采用CEL大變形有限元方法對(duì)于上硬下軟雙層黏土地基中方形樁靴貫入過程進(jìn)行了數(shù)值模擬，如圖1所示，比較系統(tǒng)地探討了上下土層強(qiáng)度之比sut/sub、上層土厚度與樁靴寬度之比H/B、上下土層重度γt、γb及靜止側(cè)壓力系數(shù)K0等指標(biāo)對(duì)于樁靴貫入阻力曲線的影響。如圖1所示，樁靴貫入過程中可能將基礎(chǔ)底面以下一部分上層黏土壓入下臥層，這部分上層黏土的厚度稱之為土塞厚度hp，D為貫入深度。

1 有限元模型

由于所研究問題的對(duì)稱性，因此，選取樁靴與地基耦合系統(tǒng)的1/4建立三維CEL有限元模型。對(duì)于地基土體所在區(qū)域采用8節(jié)點(diǎn)縮減積分Eulerian單元進(jìn)行離散，樁靴采用離散剛體Lagrangian單元模擬?；A(chǔ)寬度取B=5 m，為了避免邊界效應(yīng)的影響，地基土所在區(qū)域?qū)挾燃昂穸确謩e取為5B和6B。在有限元模型的底面邊界約束豎向和水平向自由度，在側(cè)面邊界上約束水平向自由度，基礎(chǔ)近場(chǎng)土體局部區(qū)域加密網(wǎng)格。Qiu等[9]通過計(jì)算和分析認(rèn)為樁靴表面光滑或者粗糙對(duì)于得到的貫入阻力曲線影響不大，因此，本文假定基礎(chǔ)表面光滑，基礎(chǔ)與地基土體之間的界面通過基于罰函數(shù)方法的廣義接觸算法來模擬[11]，如式（1）所示。

式中：Fp為界面上對(duì)應(yīng)點(diǎn)之間的接觸力；kp為罰剛度，其值與Lagrangian和Eulerian介質(zhì)材料特性有關(guān)，dp為穿入深度。這種算法允許Eulerian介質(zhì)穿入Lagrangian區(qū)域少量深度，因此，不如動(dòng)態(tài)接觸方法嚴(yán)格，但是，可以保證在困難情況下算法的收斂。值得注意的是，這種廣義接觸算法在Lagrangian和Eulerian單元之間不施加接觸約束，Lagrangian單元可以在不賦予介質(zhì)的Eulerian網(wǎng)格內(nèi)自由運(yùn)動(dòng)。

對(duì)于不排水條件下的黏土地基，采用基于Tresca屈服準(zhǔn)則的理想彈塑性本構(gòu)模型，泊松比取為v=0.49，變形模量E=500su，su為黏土不排水抗剪強(qiáng)度。沒有特別聲明情況下，在本文分析中，上層土強(qiáng)度sut=100 kPa保持不變，通過改變下層土強(qiáng)度sub來實(shí)現(xiàn)不同的土層強(qiáng)度比sub/sut。

樁靴的貫入通過在剛體單元上施加速率邊界條件來實(shí)現(xiàn)，將樁靴基礎(chǔ)勻速壓入土層中，進(jìn)而得到基礎(chǔ)反力與貫入深度之間的關(guān)系。在工程中，樁靴的貫入速率一般在2 m/h的量級(jí)[8]，且基本上是一個(gè)擬靜力力學(xué)行為。而在ABAQUS中，CEL有限元采用了基于中心差分準(zhǔn)則的顯式動(dòng)力學(xué)解法，其穩(wěn)定時(shí)間步長主要取決于模型中單元尺寸大小與土的彈性性質(zhì)[12]。所以，若在有限元分析中采用與工程實(shí)際相同的貫入速率，將導(dǎo)致計(jì)算時(shí)間過長而不可接受。同時(shí)，貫入速率若取得過高的話，數(shù)值模擬結(jié)果將因?yàn)橐氲娘@著的慣性力而嚴(yán)重失真。為了不引起顯著的偽慣性效應(yīng)，同時(shí)縮短計(jì)算時(shí)間，經(jīng)過試算將貫入速率取c=0.25 m/s，這介于Tho等[8]、Qiu等[8]建議的取值范圍之間。

2 上層土厚度的影響

為了探討上層硬土層厚度對(duì)于樁靴貫入阻力的影響，暫時(shí)假定土層重度γt=γb=0，在下層土與上層土不排水強(qiáng)度之比sub/sut=0.25的情況下，分別取上層土厚度與基礎(chǔ)寬度之比為H/B=0、0.2、0.5、1和∞，計(jì)算得到的貫入阻力曲線如圖2所示，其中承載力系數(shù)Nc=F/Asut，F(xiàn)為基礎(chǔ)反力，A為方形樁靴基底面積。其中，H/B=0為不排水強(qiáng)度為sub的均質(zhì)土，H/B=∞實(shí)際上代表了不排水強(qiáng)度為sut的均質(zhì)土。

通過圖2可以看到，在H/B=0和∞，即均質(zhì)土情況下，貫入阻力隨著貫入深度的增加而增加，最終趨于深基礎(chǔ)的承載力極限值，不會(huì)發(fā)生穿透現(xiàn)象。而在H/B=0.2、0.5和1，即上硬下軟層狀土情況下，承載力系數(shù)隨著深度增加達(dá)到了峰值，而后開始出現(xiàn)下降，表現(xiàn)出一種脆性響應(yīng)，這表明方形樁靴在實(shí)際貫入時(shí)的預(yù)壓荷載一旦超過曲線上的峰值，就將產(chǎn)生較大沉降，這種穿透現(xiàn)象可能造成上部結(jié)構(gòu)的傾覆破壞。

根據(jù)Wang等[13]的研究，在貫入深度足夠大時(shí)，光滑矩形深基礎(chǔ)承載力系數(shù)將趨近于一個(gè)極限值Ncult。該值可采用Meyerhof[14]建議的解析解近似表達(dá)。

式中：B、L分別為基礎(chǔ)的寬度與長度，對(duì)于方形基礎(chǔ)，Ncult=9.53。對(duì)于H/B=∞情況，有限元計(jì)算得到的極限值Ncult=10.6，比Meyerhof的解析解高10%左右，這也說明了大變形有限元模型的可靠性。

本文著重考察穿透現(xiàn)象，同時(shí)，為了節(jié)省計(jì)算費(fèi)用，沒有對(duì)樁靴進(jìn)入下層土層之后的貫入特性進(jìn)行詳細(xì)分析。但是，通過圖2中的H/B=0、0.2、0.5和1情況仍然可以看到，這幾種情況下的貫入曲線在進(jìn)入下層土深度較大時(shí)，有趨向于極限值2.65的趨勢(shì)，這與 解析解也是比較接近的。

3 土層之間強(qiáng)度比的影響

為了探討土層之間強(qiáng)度比對(duì)于樁靴貫入阻力的影響，在上層土厚度與基礎(chǔ)寬度之比H/B=0.2、0.5和1情況下，分別取sub/sut=0.1、0.25、0.5、0.7和1，計(jì)算得到的貫入阻力曲線如圖3所示。

通過圖3可知，對(duì)于無重土，在H / B=0.2、0.5和1 3種情況下，當(dāng)sub/sut<0.5時(shí)，貫入阻力曲線的Nc值均在達(dá)到峰值后，出現(xiàn)明顯的下降，即產(chǎn)生了比較明顯的穿透現(xiàn)象，而在sub/sut>0.5時(shí)，不會(huì)發(fā)生穿透問題。sub/sut=0.5時(shí)，在上層土厚度較大情況下，例如H/B=0.5和1，發(fā)生了穿透現(xiàn)象，但在上層厚度較小情況下，此時(shí)硬土層中尚未來得及產(chǎn)生峰值，樁靴已進(jìn)入下層土，因此沒有明顯的穿透問題，如圖3（a）所示。

根據(jù)Wang等[13]的研究，基礎(chǔ)進(jìn)入上硬下軟地基中，承載力主要受兩方面因素的影響：基底下土體的抗剪強(qiáng)度與貫入過程中被壓迫到基礎(chǔ)周圍的土體所提供的強(qiáng)度及超載作用。這兩個(gè)因素共同作用，決定了貫入阻力曲線的趨勢(shì)。對(duì)于無重土，不考慮基底以上土體的超載效應(yīng)。圖4給出了H/B=1、sub/sut=0.1時(shí)，不同貫入深度情況下土體變形模式，其中深色區(qū)域代表硬土，淺色區(qū)域代表軟土。圖4（a）對(duì)應(yīng)著貫入阻力曲線上的峰值，此時(shí)基底下土塞厚度hp基本保持hp/B=1，再加上推擠開的土體所提供的強(qiáng)度貢獻(xiàn)，因此承載力達(dá)到了最大值，但隨后基底下硬土層發(fā)生側(cè)向塑性流動(dòng)，其厚度逐漸變薄，在D/B=0.5時(shí)hp/B=0.9，漸變?yōu)榘雮€(gè)倒置的橄欖球形狀，在D/B=1.2時(shí)hp/B=0.75，所提供的承載力明顯降低，盡管此時(shí)基礎(chǔ)周圍被迫擠出的土體所提供的強(qiáng)度貢獻(xiàn)有所增加，但總體上導(dǎo)致承載力系數(shù)出現(xiàn)下降。endprint

4 土重的影響

隨著方形樁靴的貫入，基底以上土體重力所產(chǎn)生的超載效應(yīng)對(duì)于樁靴的承載力有著較大影響。為了研究土重的影響，在H/B=1、sub/sut=0.25條件下，分別取如下幾種情況進(jìn)行計(jì)算：1）γt=γb=10 kN/m3，γtB/sut=0.5；2）γt=γb=20 kN/m3，γtB/sut=1；3）γt=γb=10 kN/m3，γtB/sut=1；4）γt=γb=10 kN/m3，γtB/sut=2；5） γt=γb=10 kN/m3，γtB/sut=2.5，計(jì)算過程中土體初始地應(yīng)力靜止側(cè)壓力系數(shù)取K0=1，所得到的計(jì)算結(jié)果如圖5所示。圖5中同時(shí)給出了Yu等[6]建議的方形基礎(chǔ)貫入阻力預(yù)測(cè)公式，如式（3）所示。

通過比較可知，在上下層土重相同情況下，當(dāng)γtB/sut≤1時(shí)，計(jì)算結(jié)果與已有解析解在深度D/B<0.8情況下比較吻合，而當(dāng)γtB/sut>1時(shí)，解析解則低估了方形樁靴承載力。另外，從圖5中可以看到，土重對(duì)于穿透現(xiàn)象的發(fā)生具有明顯的抑制作用，這主要是因?yàn)楸粔浩鹊交A(chǔ)周圍的土體的超載作用對(duì)于承載力起到了補(bǔ)償作用。當(dāng)無量綱系數(shù)γtB/sut>1時(shí)，貫入阻力曲線沒有出現(xiàn)所謂的脆性特征。

在有限元計(jì)算中，工況（2）與工況（3）所得貫入阻力曲線幾乎完全一致，這也說明了在上下層土重相同情況下，將無量綱系數(shù)γtB/sut作為穿透現(xiàn)象判別依據(jù)的合理性。

對(duì)于上下層土體重度不同γt≠γb的情況，在H/B=1、sut=100 kPa與sub/sut=0.25條件下，主要考察了兩種情況：1）γt=10 kN/m3，γb=20 kN/m3；2）γt=20 kN/m3，γb=10 kN/m3；計(jì)算結(jié)果如圖6所示。為了比較起見，將γt=γb=10 kN/m3、γt=γb=20 kN/m3兩種情況的曲線也列在圖6中。通過比較可以看到，在貫入深度范圍內(nèi)，γt≠γb兩種情況下的曲線幾乎重合，并且介于γt=γb=10 kN/m3和γt=γb=20 kN/m3兩種情況之間。通過局部放大可見，γt=20 kN/m3，γb=10 kN/m3情況的承載力略大于γt=10 kN/m3，

γb=20 kN/m3情況，也產(chǎn)生了比較明顯的穿透特征，因此式（3）以及無量綱系數(shù)γtB/sut在上下層土重相差不大時(shí)可以作為穿透現(xiàn)象判別依據(jù)，但在上下層土重相差較大時(shí)不適用。

5 靜止側(cè)壓力系數(shù)的影響

對(duì)于黏土地基承載力進(jìn)行小變形有限元計(jì)算時(shí)，通常認(rèn)為土體初始應(yīng)力場(chǎng)靜止側(cè)壓力系數(shù)K0對(duì)計(jì)算結(jié)果沒有影響[15-16]。為了探討靜止側(cè)壓力系數(shù)對(duì)于方形樁靴貫入阻力的影響，在H/B=1、sub/sut=0.25、γt=γb=10 kN/m3及γtB/sut=0.5條件下，分別取K0=0.5和1.0兩種情況進(jìn)行計(jì)算，結(jié)果如圖7所示。通過圖7可以看到，兩種情況下的曲線在D/B>0.8后逐漸趨于重合，但是K0=0.5時(shí)的穿透荷載比K0=1情況略低6%，從而導(dǎo)致其貫入阻力曲線的脆性特征有所減弱。根據(jù)前述穿透現(xiàn)象發(fā)生時(shí)的土體變形機(jī)理解釋，樁靴貫入阻力主要來自于基底下土體的抗剪強(qiáng)度與貫入過程中被壓迫到基礎(chǔ)周圍的土體所提供的強(qiáng)度及超載作用。側(cè)壓力系數(shù)K0的減小降低了受壓迫向基礎(chǔ)周圍擠出的土體的側(cè)向約束作用，因此，導(dǎo)致承載力系數(shù)有所下降。

6 結(jié)論

1）上下土層強(qiáng)度比與上層土厚度對(duì)于樁靴的貫入阻力曲線具有很大影響。

2）靜止側(cè)壓力系數(shù)的減小在一定程度上導(dǎo)致樁靴承載力系數(shù)的降低，其貫入阻力曲線的脆性特征有所減弱。

3）上下土層重度的差異導(dǎo)致樁靴的貫入阻力曲線不滿足前人建議的方程。無量綱系數(shù)γtB/sut值在上下層土重度相差不大時(shí)，可作為穿透現(xiàn)象判別依據(jù)。

參考文獻(xiàn)：

[1] GOURVENEC S， RANDOLPH M F， KINGSNORTH O. Undrained bearing capacity of square and rectangular footings[J]. International Journal of Geomechanics， 2006， 6（3）： 147-157.

[2] MERIFIELD R S， NGUYEN V Q. Two- and three-dimensional bearing-capacity solutions for footings on two-layered clays[J]. Journal of Geomechanics and Geoengineering， 2006， 1（2）： 151-162.

[3] 張愛霞. 方形樁靴自升式平臺(tái)地基極限承載力及結(jié)構(gòu)抗穿刺研究[D]. 遼寧 大連： 大連理工大學(xué)， 2013.

ZHANG A X. On ultimate bearing capacity of jack-up with square spudcan and application to structure investigation against punch-through[D]. Dalian，Liaoning： Dalian University of Technology， 2013. （in Chinese）

[4] CRAIG W H， CHUA K. Deep penetration of spud-can foundations on sand and clay [J]. Geotechnique， 1990， 40（4）： 541-556.

[5] LEE K K， CASSIDY M J， RANDOLPH M F. Bearing capacity on sand overlying clay soils： experimental and finite-element investigation of potential punch-through failure [J]. Geotechnique， 2013， 63（15）： 1271-1284.endprint