徐保照, 李 颯, 蘇洪路

( 1.天津港東疆建設(shè)開(kāi)發(fā)有限公司, 天津 300463; 2.天津大學(xué)建筑工程學(xué)院, 天津 300072; 3.天津二建建筑工程有限公司, 天津 300143 )

近幾十年來(lái)，海洋油氣資源的開(kāi)發(fā)成為現(xiàn)在能源開(kāi)發(fā)的熱點(diǎn)，海洋油氣的開(kāi)發(fā)促成了越來(lái)越多海洋平臺(tái)的建造逐步向深海發(fā)展。中國(guó)海域面積廣闊，土質(zhì)的分布也比較特殊，從一般的粘土、砂土到與粘土和砂土性質(zhì)不同的迭層土[1]，再到易破碎的珊瑚礁、鈣質(zhì)砂以及硬的鈣質(zhì)巖(礁灰?guī)r)，每種土都有各自的特殊性，鈣質(zhì)土中的鈣質(zhì)砂與一般的粘土及砂土最大的不同是，鈣質(zhì)砂在受力的情況下，表現(xiàn)出顆粒易碎、強(qiáng)度、壓縮及循環(huán)荷載下強(qiáng)度衰減等特性,在該土質(zhì)條件下需要慎重選擇錨固基礎(chǔ)。

重力錨是一種應(yīng)用比較早的錨固形式，主要靠自重下錨底與土之間的摩擦力來(lái)提供反力[2]，其結(jié)構(gòu)簡(jiǎn)單，由混凝土塊或碎金屬或其他高密度材料制成。對(duì)于含有礁灰?guī)r的鈣質(zhì)土，重力錨可以較好地適用于這種土質(zhì)條件。在國(guó)外有不少學(xué)者對(duì)重力錨進(jìn)行了相關(guān)的研究，Taylor R[3]設(shè)計(jì)了帶剪力鍵的重力錨，給出了重力錨在粘土中的承載力計(jì)算公式。Harris R E等[4]為波浪能轉(zhuǎn)換器(WECs)設(shè)計(jì)錨泊系統(tǒng)，討論了該體系所遇到疲勞荷載和磨損以及從工程造價(jià)和土質(zhì)適用上探討了傳統(tǒng)的錨泊系統(tǒng)對(duì)該體系的適用性，得出的結(jié)論是重力錨是比較優(yōu)越的。Taylor R等[5]介紹了熱帶島嶼或者火山環(huán)礁附近海域上的海洋熱能轉(zhuǎn)換器(OTEC)可能的錨固設(shè)計(jì)方案，根據(jù)環(huán)境條件，設(shè)計(jì)了拖曳嵌入式錨和帶有裙翼的重力錨，并探討其外形尺寸的選擇。Michael G S[6]為東南佛羅里達(dá)海洋再生能源體系設(shè)計(jì)了重力錨泊系統(tǒng)，主要設(shè)計(jì)了單錨的外形尺寸并進(jìn)行了強(qiáng)度驗(yàn)算。Jim B. Petersen[7]發(fā)明了用錨鏈連接在一起的兩個(gè)或者多個(gè)錨單元，用來(lái)增加錨抓力；同時(shí)，錨單元在安裝過(guò)程中可以轉(zhuǎn)動(dòng)，以確保錨以正確的方向安裝在海底。Ben U[8]在華盛頓州的SR520常青點(diǎn)浮橋和西雅圖登陸項(xiàng)目設(shè)計(jì)中，因地制宜地設(shè)計(jì)了軸鉆錨(Drilled shaft anchor)、重力錨、抓錨三種不同的錨來(lái)適應(yīng)不同的工況條件。在國(guó)內(nèi)，李懷亮等[9]對(duì)重力錨水平承載力進(jìn)行了相關(guān)的研究。對(duì)于鈣質(zhì)土的研究，近年來(lái)不少學(xué)者做了相關(guān)的研究[10-13]。

1 模型的可靠性驗(yàn)證

1.1 理論計(jì)算

為了對(duì)串錨進(jìn)行相關(guān)分析，首先需要對(duì)建立的計(jì)算模型進(jìn)行可靠性驗(yàn)證。重力錨的結(jié)構(gòu)示意圖如下：

圖1 重力錨結(jié)構(gòu)示意圖Fig.1 The structure diagram of gravity anchor

此圖為錨的斷面圖，錨主體為長(zhǎng)方體型，即蹲塊式重力錨，重力錨有安裝剪力鍵和無(wú)剪力鍵的情況，各個(gè)物理參數(shù)在下面進(jìn)行詳細(xì)說(shuō)明。

模型的驗(yàn)證采用如下方法進(jìn)行：分別用理論及室內(nèi)模型試驗(yàn)和有限元法計(jì)算重力錨的水平極限承載力，然后對(duì)比三種方法的計(jì)算結(jié)果，如果有限元計(jì)算的結(jié)果和理論計(jì)算以及室內(nèi)試驗(yàn)結(jié)果的誤差在容許的范圍內(nèi)，則認(rèn)為模型是可靠的，反之則不可靠。目前相關(guān)規(guī)范[2]給出，在非粘性土中，錨在水平荷載作用下，錨底土體發(fā)生排水剪切破壞，水平極限承載力按下式計(jì)算：

Qul=μ[(Wb+γbADf-Fve)cosβ-Fhsinβ]+Rp。

(1)

式中：Qul為錨水平極限承載力；μ為錨底與土體的摩擦系數(shù)或者是在有剪力鍵的前提下土與土之間的摩擦系數(shù)，計(jì)算公式為μ=tan(φ-5°)(在沒(méi)有剪力鍵的鋼質(zhì)或混凝土的錨底需要折減5°，有剪力鍵則不用折減)；Wb為錨在水下的重力；Fve為設(shè)計(jì)荷載與錨鏈荷載在豎直方向上的分力(豎直向上為正)；Fh為設(shè)計(jì)荷載與錨鏈荷載在水平方向上的分力(沿錨底坡向方向?yàn)樨?fù))；γbADf為剪力鍵間土體的浮重；γb為土的浮重度； A為錨底面積；Df為剪力鍵入泥深度(海底泥面以下為正)；β為海底泥面坡角；φ是錨底土體內(nèi)摩擦角；Rp為錨底前沿的被動(dòng)土壓力，Rp可由下式計(jì)算：

(2)

式中Kp為被動(dòng)土壓力系數(shù)，計(jì)算公式為：

Kp=tan2(φ/2+45°)。

Talory對(duì)式(1)進(jìn)行了簡(jiǎn)化，當(dāng)β=0時(shí)，

Qul=μ[(Wb-Fve)]+Rp。

(3)

對(duì)于平底錨(無(wú)剪力鍵)，在水平荷載作用下，公式(3)可簡(jiǎn)化為

Qul=μWb。

(4)

此次驗(yàn)證采用單個(gè)重力錨，并且沒(méi)有錨底的剪力鍵。對(duì)于重力錨下面土體的計(jì)算相關(guān)參數(shù)，可以根據(jù)室內(nèi)試驗(yàn)結(jié)果。采用室內(nèi)試驗(yàn)方法對(duì)鈣質(zhì)砂的成分進(jìn)行測(cè)試。樣本的鈣質(zhì)砂中CaCO3含量都大于95%，鈣質(zhì)砂的重度為19 kN/m3。試驗(yàn)用鈣質(zhì)砂的顆粒級(jí)配見(jiàn)圖2，采用此級(jí)配的鈣質(zhì)砂在不同的相對(duì)密實(shí)度下進(jìn)行直剪試驗(yàn)，得到鈣質(zhì)砂內(nèi)摩擦角隨相對(duì)密實(shí)度變化的關(guān)系見(jiàn)圖3。根據(jù)上述試驗(yàn)結(jié)果，確定重力錨下土體的計(jì)算參數(shù)。

圖2 鈣質(zhì)砂的顆粒級(jí)配Fig.2 The particle gradation of calcareous sand

圖3 鈣質(zhì)砂內(nèi)摩擦角隨相對(duì)密實(shí)度的變化關(guān)系Fig.3 The relationship of internal friction angle and relative density of calcareous sand

1.2 有限元模型建立

為了分析重力錨的水平承載機(jī)理，采用有限元方法對(duì)承受水平拉力的重力錨進(jìn)行分析。計(jì)算采用ABAQUS 有限元軟件，土體單元類(lèi)型采用C3D8R 實(shí)體單元，土體的計(jì)算范圍根據(jù)錨體的外形尺寸確定，計(jì)算深度與計(jì)算寬度均取約為5倍錨的相應(yīng)高度和寬度以消除邊界的影響，土體尺寸選定為20 m × 20 m × 10 m(長(zhǎng)×寬×高)。土體采用M-C(摩爾-庫(kù)倫) 本構(gòu)模型，土重度取19 kN/m3(浮重度9 kN/m3)，內(nèi)摩擦角取35°，相當(dāng)于密實(shí)鈣質(zhì)砂的內(nèi)摩擦角，彈性模量30 MPa，泊松比0.25。為了提高模擬精度，土體單元采用漸進(jìn)加密的單元?jiǎng)澐?。重力錨干重84.78 t，錨的尺寸為3.0 m×3.0 m×1.2 m，對(duì)于重力錨尺寸的要求相關(guān)規(guī)范[2]有詳細(xì)說(shuō)明，其材質(zhì)為鋼材，其本構(gòu)模型采用性彈性體，其密度為7 850 kg/m3(考慮重力錨在海底受到浮力，計(jì)算時(shí)采用水下密度6 850 kg/m3)，彈性模量取2.1×105MPa，泊松比0.3。

重力錨與土體之間設(shè)置接觸面，接觸面采用surface-to-surface 類(lèi)型，在接觸屬性定義切向和法向接觸面特性，對(duì)于切向特性，平底錨采用土體內(nèi)摩擦角折減5°取正切；對(duì)于法向特性，采用硬接觸。土體底部約束三個(gè)方向上的位移，四周立面約束水平方向上的位移。在計(jì)算分析過(guò)程中，均考慮重力錨應(yīng)用于懸鏈線(xiàn)式系泊系統(tǒng)，即重力錨只承受水平荷載，F(xiàn)ve = 0，海底地面平坦，即β = 0，水平力施加在錨前面形心位置。

根據(jù)提取的錨上水平反力和位移對(duì)應(yīng)關(guān)系，可繪制有限元法所計(jì)算的重力錨水平極限承載力-位移圖見(jiàn)圖4。

圖4 單錨位移-水平極限承載力關(guān)系Fig.4 The relationship of displacement-horizontal bearing capacity of single gravity anchor

從上圖可以看出，重力錨在水平荷載作用下，錨底土體發(fā)生滑動(dòng)破壞，圖中水平位移-承載力關(guān)系曲線(xiàn)有明顯的拐點(diǎn)，從圖中可得到重力錨的水平極限承載力為390.201 kN。

1.3 模型試驗(yàn)

重力錨構(gòu)件尺寸為3.0 m×3.1 m×1.2 m。模型幾何比尺為1∶15，所以重力錨模型試件尺寸取：0.20 m×0.20 m×0.10 m。試驗(yàn)裝置由流速控制系統(tǒng)、量測(cè)系統(tǒng)、計(jì)算機(jī)數(shù)據(jù)采集系統(tǒng)等幾部分組成。試驗(yàn)過(guò)程基本原理為：設(shè)定水流流速，通過(guò)調(diào)節(jié)變頻電動(dòng)機(jī)轉(zhuǎn)速?gòu)亩淖兯昧髁?，同時(shí)在管道內(nèi)布置畢托管檢測(cè)實(shí)際流速，直到達(dá)到水槽設(shè)定水流流速。重力錨后端布置拉線(xiàn)式位移傳感器，監(jiān)測(cè)重力錨運(yùn)動(dòng)；在重力錨前端的纜繩上布置拉力傳感器。待水流流速平穩(wěn)后，利用數(shù)據(jù)采集器將傳感器上的信號(hào)傳送到電腦上。本次試驗(yàn)主要通過(guò)水壓加荷。通過(guò)定滑輪在水槽末端懸掛一塑料水箱，利用水泵向水箱內(nèi)勻速加水，實(shí)現(xiàn)對(duì)重力錨試件的勻速加載。

試驗(yàn)過(guò)程按照相關(guān)的規(guī)范[14-17]操作進(jìn)行，錨靜止在水槽鈣質(zhì)砂上，等土中超靜水孔壓完全消失后，開(kāi)始均勻施加水平荷載，至錨移動(dòng)后觸發(fā)行程開(kāi)關(guān)，加載停止，錨在水平荷載下移動(dòng)一定距離后將停下(見(jiàn)圖5)，此時(shí)可認(rèn)為錨下面土體已破壞，由于水流速度很小，對(duì)錨的阻力很小，可以忽略，故錨移動(dòng)以前的荷載為錨水平承載力極限值。

圖5 錨在水平荷載作用下滑動(dòng)破壞Fig.5 Sliding failure of gravity with horizontal load

根據(jù)計(jì)算機(jī)采集的數(shù)據(jù)，將錨的水平荷載和位移關(guān)系繪制如下圖。

圖6 錨位移-荷載關(guān)系Fig.6 The relationship between the displacement and horizontal bearing capacity of gravity anchor

從圖中可以看出，水平荷載與位移的關(guān)系曲線(xiàn)有明顯的轉(zhuǎn)折點(diǎn)，水平極限承載力為120.21 N。

模型試驗(yàn)中錨的總干重25.004 kg，按照相似三定理的關(guān)系，原型錨的質(zhì)量為84.38 t，按公式(4)可計(jì)算得模型錨水平極限承載力為125.87 N，有限元計(jì)算的結(jié)果也按相似三定理?yè)Q算，其結(jié)果為115.62 N，相對(duì)誤差分析見(jiàn)下表1。

表1 不同方法所得錨水平極限承載力Table 1 horizontal ultimate bearing capacity of gravity with different method

Note：①Computing method；②Horizontal ultimate bearing capacity；③Relative deviation (FEM compared with theoretical and model test respectively)；④Theoretical calculation；⑤Finite element method；⑥Model test.

從表中可以看出，不論是與理論計(jì)算還是模型試驗(yàn)結(jié)果相比，有限元法所計(jì)算結(jié)果的相對(duì)誤差都在10%以?xún)?nèi)，因此，可以確信有限元模擬計(jì)算的結(jié)果是可靠的。

2 串錨承載力機(jī)理特性

串錨是將多個(gè)重力錨單元用錨鏈串聯(lián)起來(lái)，本次主要研究由兩個(gè)重力錨單元串聯(lián)起來(lái)的串錨承載力特性，兩個(gè)錨單元的尺寸仍采用上面重力錨的尺寸，錨鏈的長(zhǎng)度分別為10、20、40、60和80 m，計(jì)算的土質(zhì)條件仍不變，對(duì)應(yīng)土體尺寸分別為40 m×20 m×10 m、 60 m×20 m×10 m、80 m×20 m×10 m、100 m×20 m×10 m、120 m×20 m×10 m。計(jì)算模型見(jiàn)下圖：

(①2# anchor unit；②A(yíng)nchor chain；③1# anchor unit；④Horizontal force；⑤Soil.)

圖7 串錨有限元模型示意圖
Fig.7 Finite element model of multiple gravity anchors

在A(yíng)BAQUS中，錨鏈常用truss單元來(lái)模擬，truss單元為桁架單元，所有的單元的節(jié)點(diǎn)均為鉸接點(diǎn)，只有平動(dòng)自由度，這就意味著truss單元沒(méi)有轉(zhuǎn)動(dòng)方向的剛度，只能受拉或受壓。對(duì)于本次模擬的錨鏈，只受拉，故可以用truss單元來(lái)模擬，錨鏈的等效截面積為0.007 85 m2，本次研究的錨鏈材質(zhì)為鋼質(zhì)，錨鏈直接與前后兩個(gè)錨立面形心連接以簡(jiǎn)化計(jì)算模型，連接方式為錨鏈兩端的節(jié)點(diǎn)與前后錨單元的立面耦合(Coupling)。為了模擬可操作性以及橫向可對(duì)比性，錨鏈的初始狀態(tài)均為剛好繃直(相對(duì)水平荷載，錨鏈自重可以忽略)。水平荷載施加在1號(hào)錨前面形心位置，與單錨情況相同。分別提取1號(hào)錨單元受力點(diǎn)的水平荷載和位移，繪制水平荷載和位移關(guān)系圖如下圖所示：

圖8 不同錨鏈長(zhǎng)度的串錨水平承載力與水平位移關(guān)系

從上圖可看出，串錨的水平承載力與位移關(guān)系和單錨類(lèi)似，圖中所有曲線(xiàn)有明顯的拐點(diǎn)，錨鏈長(zhǎng)度為10 m的串錨水平承載力最小，從圖中可以得到不同長(zhǎng)度的錨鏈的串錨水平承載力，見(jiàn)下表2。

表2 不同錨鏈長(zhǎng)度的串錨水平承載力與水平位移Table 2 the relationship between horizontal bearing capacity and displacement with different length of chain

注：錨鏈長(zhǎng)度為10 m的串錨水平極限承載力小于700 kN，故水平反力為700 kN時(shí)錨底土體已破壞，和其他工況的串錨水平位移沒(méi)有比較意義。

Note:①Chain length of multiple anchors；②Horizontal ultimate bearing capacity；③Horizontal displacement under horizontal ultimate bearing capacity；④Horizontal displacement under horizontal reaction 700 kN.

很明顯，總體對(duì)比起來(lái)，隨錨鏈的長(zhǎng)度增加，串錨的水平承載力有增加的趨勢(shì)，當(dāng)錨鏈達(dá)到一定長(zhǎng)度，承載力基本不變，總體上來(lái)講，錨鏈越長(zhǎng)，串錨達(dá)到水平極限承載力時(shí)發(fā)生的位移也越大，達(dá)到同樣的水平承載力時(shí)，同樣錨鏈越長(zhǎng)，錨發(fā)生的位移越大。對(duì)比單錨和串錨水平承載力，見(jiàn)下表3。

表3 串錨承載力與單錨對(duì)比Table 3 the horizontal bearing capacity of single gravity anchor and multiple gravity anchor

注：這里的折減系數(shù)是串錨水平承載力與單錨水平承載力2倍的比值。

Note:Reduction coefficient=Horizontal ultimate bearing capacity of China anchor÷(Horizontal ultimate bearing capacity of single anchor ×2).

從上表3可以看出，所有串錨的水平承載力均小于單錨水平承載力的2倍，也就是說(shuō)串錨的承載力不是單個(gè)錨承載力的簡(jiǎn)單相加，造成這種情況的原因可以從水平承載力機(jī)理上解釋。

首先可以研究單錨承載力機(jī)理，下面兩個(gè)圖為單錨在水平荷載作用下發(fā)生滑動(dòng)破壞后土體的等效塑性應(yīng)變?cè)茍D(PEEQ)和塑性屈服圖(AC YIELD)。

圖9 單錨錨底土體等效塑性應(yīng)變?cè)茍D和塑性屈服圖Fig.9 PEEQ and AC YIELD of the soil under single gravity anchor

上面的剖面是沿錨水平軸線(xiàn)縱向剖開(kāi)，荷載施加方向?yàn)樗较蛴?。為了清楚看到土體破壞后各項(xiàng)指標(biāo)，故將錨移除。從上圖可以明顯看出，在水平荷載作用下，錨底土體發(fā)生滑動(dòng)破壞，破壞面為弧形，靠近前端，破壞面較深。

圖10分別是錨鏈長(zhǎng)度10、20、40、60和80 m串錨錨底土體等效塑性應(yīng)變?cè)茍D。圖10均是右邊是1號(hào)錨錨底土體，左邊是2號(hào)錨錨底土體，由于土體過(guò)于狹長(zhǎng)，只截取了錨底一定范圍內(nèi)的土體。

從上圖(a)到(e)可看出，串錨的每個(gè)錨單元破壞模式基本上和單錨相同，都是錨底土體發(fā)生滑動(dòng)破壞，不同的是，前面的錨單元錨底土體塑性應(yīng)變值較后面的錨單元大，破壞面也較后面的大；隨著錨鏈長(zhǎng)度的增加，這種差距變的較小。其原因是：在水平荷載作用下，由于錨鏈的存在，每個(gè)錨單元受力是不一樣的，前面的錨單元受力大于后面的錨單元，也就是說(shuō)，隨著水平荷載作用的增加，前面的錨單元會(huì)先發(fā)生錨底土體滑動(dòng)破壞，導(dǎo)致水平滑動(dòng)失效，隨后將力傳給后面的錨單元，后面的最后也發(fā)生滑動(dòng)破壞。因此導(dǎo)致串錨的水平承載力不是單個(gè)錨水平承載力相加。對(duì)于不同錨鏈長(zhǎng)度的串錨，在一定范圍內(nèi)錨鏈越長(zhǎng)，破壞時(shí)前面錨單元發(fā)生的位移越大，這樣后面的錨單元能夠充分發(fā)揮承載力效應(yīng)。

圖10 不同錨鏈長(zhǎng)度串錨錨底土體等效塑性應(yīng)變?cè)茍DFig.10 PEEQ of the soil under multiple gravity anchor with different chain length

3 結(jié)論

通過(guò)有限元對(duì)串錨在水平荷載作用下破壞過(guò)程的模擬，可以得出下面幾點(diǎn)結(jié)論：

(1)串錨在水平荷載作用下，首先是前面施加荷載的錨單元先發(fā)生滑動(dòng)破壞，隨后是后面的錨單元發(fā)生滑動(dòng)破壞，即前面錨單元受力較后面的大，滑動(dòng)破壞模式和相應(yīng)的單錨基本相同。

(2)串錨的水平承載力與錨鏈長(zhǎng)度有關(guān)，在一定范圍內(nèi)，增加錨鏈的長(zhǎng)度可以提高串錨的水平承載力，超過(guò)某一長(zhǎng)度后，錨鏈的長(zhǎng)度對(duì)串錨的承載力影響不大。

(3)串錨的水平承載力不是單個(gè)錨的水平承載力相加，其水平承載力小于相應(yīng)的單個(gè)錨水平承載力之和，在進(jìn)行串錨水平承載力設(shè)計(jì)計(jì)算時(shí)應(yīng)給與適當(dāng)?shù)目紤]。

針對(duì)串錨的破壞模式以及錨鏈長(zhǎng)度對(duì)水平承載力的影響，在此文的參考下，需要用模型試驗(yàn)來(lái)進(jìn)一步確認(rèn)錨鏈最佳長(zhǎng)度，這正是下一步所要做的研究工作。