徐 峰，譚紅瑩，趙佳寧，王江宏，陳 曦

(海洋石油工程股份有限公司 天津 300451)

1 研究背景

深水系泊系統(tǒng)中，嵌入水底土中的錨鏈由于土壤提供的側(cè)向支撐而呈現(xiàn)典型的反懸鏈線狀，如圖1所示。該段錨鏈與土體間復(fù)雜的相互作用[1]，改變了傳遞至系泊基礎(chǔ)上的荷載大小與角度，并最終影響系泊基礎(chǔ)的承載力與失效模式。

圖1 錨固基礎(chǔ)受力示意圖Fig.1 Schematic diagram of force acting on anchored foundation

對于嵌入土中的錨鏈，國內(nèi)外學(xué)者已做過大量研究。Reese[2]最早研究了錨鏈與土體間的相互作用，提出了假定錨鏈形狀為多段圓弧構(gòu)成的簡單計算方法。Vivatrat等[3]建立了錨鏈微元段的微分方程，獲得嵌入土體錨鏈的反懸鏈線構(gòu)型。Degenkamp等[4]采用試驗方法研究了錨鏈系泊點與出泥點荷載之間的關(guān)系。Neubecker[5]和 Randolph[6]采取逐段微分的方法推導(dǎo)了錨鏈系泊點與出泥點荷載關(guān)系的計算公式，提出了嵌入土中錨鏈形態(tài)的計算公式。國振等[7]建立了二維錨泊線準(zhǔn)靜力分析模型，并對考慮錨泊線切入土體段的必要性進行了分析。Wang等[8]在二維和三維空間對錨鏈進行了準(zhǔn)靜力分析，并對錨鏈在土體中滑動可能對錨泊系統(tǒng)造成的影響進行了研究，指出該滑動對吸力錨的承載力有影響。

然而，對于嵌入土中反懸鏈段的計算，國內(nèi)尚未有成熟的模擬程序應(yīng)用于工程項目，通常還是由國外的工程公司完成該方面的分析。本文基于 Neubecker和Randolph提出的利用微分方式求取錨鏈力學(xué)響應(yīng)的方法[9]，建立了嵌入土中錨鏈段的數(shù)值計算模型，采用 MATLAB編制程序開展計算分析，并以工程項目為例，將數(shù)值分析結(jié)果與國外工程公司的計算結(jié)果進行對比，驗證該數(shù)值計算模型的可靠性。

2 嵌入土中錨鏈的力平衡方程

嵌入土中段錨鏈的微元段受力如圖2所示。

圖2 土體中錨鏈微元體的力平衡方程Fig.1 Force equilibrium equation of anchor chain microelement in soil

圖中，錨鏈微元段ds上荷載的減小值dT可通過下式來計算：

式中：T為錨鏈荷載；θ為錨鏈微元段方位角(θ=0°表示水平微元段)；ds為微元段長度；F為沿錨鏈切向的土壤阻力；wc為土中錨鏈單位長度的重量。

從一個錨鏈微元段到下一個錨鏈微元段的角度變化使用下式計算：

式中：Q為沿錨鏈法向的土壤阻力。

土壤阻力F和Q，分別采用下列公式計算：

式中：db為錨鏈單鏈環(huán)的名義直徑；En、Et分別為有效寬度和周長參數(shù)；f為沿錨鏈切向的單位摩阻力；q為沿錨鏈法向的單位端阻力。在粘性土中，(α為粘結(jié)系數(shù)，su為不排水抗剪強度，Nc為承載力系數(shù))。

3 MATLAB自編軟件

本文自編軟件基于2個基本假定：①錨鏈在基礎(chǔ)上的系泊點為固定點，即基礎(chǔ)承受錨鏈所傳遞的張力荷載時不會產(chǎn)生位移；②錨鏈為連續(xù)柔性，不受彎矩和扭轉(zhuǎn)，只受張緊荷載[10]。

3.1 軟件輸入?yún)?shù)

系泊點在泥面下的深度值 z(m)；泥面處的錨鏈荷載 T(N)；錨鏈與泥面的夾角θ(°)；錨鏈的軸向剛度 EA(N)；錨鏈單鏈環(huán)的名義直徑(m)；土體中錨鏈的單位長度重量(N/m)；錨鏈剖分單元長度(m)；土壤有關(guān)參數(shù)。

3.2 軟件輸出參數(shù)

傳遞到系泊點的荷載大小(N)；錨鏈張力的水平分量(N)；錨鏈張力的豎向分量(N)；錨鏈形態(tài)的水平坐標(biāo)(m)；錨鏈形態(tài)的豎向坐標(biāo)(m)；系泊點處荷載與水平面的夾角(°)；未張緊錨鏈的初始長度(m)；受載伸張后的錨鏈長度(m)。

4 算例驗證

為保證數(shù)值計算模型結(jié)果的準(zhǔn)確度和可靠性，選取了以下幾個具有代表性的實際工程項目數(shù)據(jù)，采用MATLAB自編軟件進行計算，并與國外專業(yè)巖土工程公司的計算結(jié)果進行對比驗證。

4.1 算例一

某單點系泊項目，系泊點位于泥面下 14.0m 深度處，該位置以上土壤均為粘性土。泥面處的系泊荷載、與泥面的夾角、錨鏈以及土壤的具體參數(shù)如表1—表3所示。

表1 系泊荷載Tab.1 Mooring load

表2 錨鏈參數(shù)Tab.2 Anchor chain parameters

表3 土壤參數(shù)Tab.3 Soil parameters

MATLAB自編軟件的分析結(jié)果如下：系泊點處張力：2847.7kN；系泊點處荷載與水平面夾角：22.08°。

4.2 算例二

系泊點位于泥面下 4.0m的位置，其上土壤均為砂土。泥面處的系泊荷載、與泥面的夾角、錨鏈以及土壤的具體參數(shù)如表4—表6所示。

計算中考慮了1m的整體沖刷，則系泊點位置深度變?yōu)?3m，MATLAB自編軟件的分析結(jié)果如下：系泊點處張力：10746kN；系泊點處荷載與水平面夾角：29.5°。

表4 系泊荷載Tab.4 Mooring load

表5 錨鏈參數(shù)Tab.5 Anchor chain parameters

表6 土壤參數(shù)Tab.6 Soil parameters

4.3 算例三

系泊點位于泥面下 5.5m的位置，系泊點位置以上土壤為砂土與粘土互層。泥面處的系泊荷載等參數(shù)如表7—表9所示。

表7 系泊荷載Tab.7 Mooring load

表8 錨鏈參數(shù)Tab.8 Anchor chain parameters

表9 土壤參數(shù)Tab.9 Soil parameters

計算中考慮了1m的整體沖刷，則系泊點位置深度變?yōu)?4.5m，MATLAB自編軟件的分析結(jié)果如下：系泊點處張力 10977kN；系泊點處荷載與水平面夾角 28.3°。

4.4 結(jié)果對比

本文MATLAB自編軟件計算結(jié)果與國外工程公司計算結(jié)果對比如表10所示。

由上述對比結(jié)果可以看出，算例一和算例二系泊點處荷載與水平面夾角的數(shù)值計算結(jié)果與國外工程公司計算結(jié)果誤差在 4%以內(nèi)，計算結(jié)果都略小于國外工程公司結(jié)果；算例二(全砂土)計算的系泊點處荷載值誤差較大，約為 9%，但由于本文計算的系泊點處荷載大于國外工程公司結(jié)果，水平和豎向荷載分量都大于國外工程公司的計算結(jié)果，對設(shè)計而言是保守的；算例三的土壤情況一定要引起足夠重視，表層土為砂土，其下為比較軟的粘土，其承載能力受控于下層的軟粘土，這種現(xiàn)象國外工程公司稱之為“穿透現(xiàn)象”，數(shù)值計算和系泊點荷載與國外工程公司結(jié)果接近，但夾角偏差較大，還需進一步研究和完善。

表10 計算結(jié)果對比Tab.10 Comparisons of calculation results

5 結(jié) 論

①本文建立了一種利用微分方式求取嵌入土中錨鏈段錨鏈力學(xué)響應(yīng)的數(shù)值計算模型，并通過MATLAB編程實現(xiàn)可操作化。

②通過與國外工程公司計算結(jié)果對比，驗證了MATLAB自編軟件計算結(jié)果的準(zhǔn)確度和可靠性；對于均質(zhì)土，與國外工程公司計算結(jié)果誤差較小，可用于指導(dǎo)工程設(shè)計。

③對于存在軟弱下臥層的地層，由于可能會產(chǎn)生穿透現(xiàn)象，自編軟件計算的系泊點處荷載與水平面夾角誤差較大(約為 35%)。對于此地層情況的嵌入土體中的反懸鏈線分析還需進一步研究。