海底管道后挖溝過程的應(yīng)力解析分析方法

， ， ，經(jīng)天

(1.中國石油大學(xué)(北京)機械與儲運工程學(xué)院 北京 102249; 2.中石化天津液化天然氣有限責(zé)任公司 天津 300457)

·試驗研究·

海底管道后挖溝過程的應(yīng)力解析分析方法

夏夢瑩1，張宏1，吳鍇1，吳經(jīng)天2

(1.中國石油大學(xué)(北京)機械與儲運工程學(xué)院 北京 102249; 2.中石化天津液化天然氣有限責(zé)任公司 天津 300457)

海底管道在后挖溝過程中會產(chǎn)生較大的變形，此時管道內(nèi)產(chǎn)生較大的應(yīng)力，對管道安全的運行形成威脅。因此準(zhǔn)確給出挖溝過程中管道應(yīng)力計算方法具有重要意義。基于彈性梁理論，給出了挖溝過程中管道撓曲線微分方程，通過邊界條件與協(xié)調(diào)方程，給出了考慮幾何非線性的挖溝過程中管道應(yīng)力與變形解析計算方法。通過與有限元數(shù)值模型對比得到本解析方法結(jié)果與數(shù)值方法結(jié)果基本一致，可以用于工程中海底管道的應(yīng)力計算。

海底管道；挖溝敷設(shè)；應(yīng)力；解析方法；數(shù)值驗證

0 引 言

海底管道在海洋油氣開發(fā)過程中擔(dān)任著重要的作用。管道的挖溝埋設(shè)，是指用水下挖溝機械在海床上挖出一條指定深度、指定形狀的溝槽，并將海底管道埋藏在溝槽之中，以達到保護管道的目的[1]。這種敷設(shè)在我國近海得到了廣泛的應(yīng)用，挖溝敷設(shè)后的管道能夠更好的避免拋錨等第三方事故對管道的影響[2]。

由于海底管道挖溝過程中管道內(nèi)會產(chǎn)生較大的變形與應(yīng)力，從而保證管道在挖溝過程中的強度安全有著重要的實際意義。然而針對這方面的研究較少[3,4]。帥健[5]提出了一種管道挖溝狀態(tài)下應(yīng)力解析分析方法，但是其沒有考慮軸向變形對管道應(yīng)力的影響。劉秉德[6]基于AutoPIPE軟件開展了X60管道的挖溝過程的數(shù)值模擬分析，給出了最大彎曲應(yīng)力與變形的變化趨勢，然而其結(jié)果的適用性有限。總而言之，現(xiàn)有的研究對工程實際的應(yīng)用范圍有限。因此，本文給出一種通用的海底管道在挖溝過程中應(yīng)力分析的解析方法。方法基于彈性梁理論確定管道的彎曲變形，基于伸長量協(xié)調(diào)條件考慮軸向應(yīng)力的影響。基于本方法能夠準(zhǔn)確的得到海底管道挖溝過程中的應(yīng)力與變形結(jié)果。

1 管道力學(xué)模型

挖溝條件下管道變形形狀如圖1所示，溝深為Δ，管子在海水中的單位長度重量為q，AB段為上部懸空段，設(shè)其長度為l1；BC段為下部懸空段，設(shè)其長度為l2。假設(shè)基礎(chǔ)是剛性的，從而可以得到躺在剛性基礎(chǔ)上的部分管道存在彎矩為零的邊界條件。

最終可以得到彎曲段的受力如圖所示，在A、B、C點分別作用有三個集中力RA、RB和RC，這三個集中力和圖中的懸空長度l1和l2均為未知量。

圖1 挖溝條件下管道變形示意圖

1.1 彎曲變形分析

以A點為坐標(biāo)原點，建立坐標(biāo)系，對于AB段，彎曲微分方程為

(1)

代入段的邊界條件y1″|x=0=0,y1′|x=0=0,y1|x=0=0,y1|x=l1=0

對于BC段，彎曲微分方程為

EIy2″=RAx+RB(x-l1)-1/2qx2,l1≤x≤l1+l2

(2)

在B點處轉(zhuǎn)角連續(xù)，可以得得到：C3=0。

BC段邊界條件為：y2|x=l1=0,y2|x=l1+l2=-Δ,y2′|x=l1+l2=0,y2″|x=l1+l2=0

從而可以得到

RB=1/2l2q(l1+l2)(1/2l1+l2)

(3)

從而可以得到：l1=1.96(ΔEI/q)1/4

各段彎矩

(4)

從而可以得到最大彎矩作用在B處，其大小為

(5)

從而可以得到，最大彎矩隨管道的單位長度的水下重量、抗彎剛度和溝深的增大而增大。

同時根據(jù)對B點求矩的平衡可以得到：

1.2 管道軸向變形

1.2.1 幾何伸長量

管道的幾何伸長量主要來自于管道變形后曲線長度的增加，可以通過弧長積分法得到變形后的曲線長度，從而得到管道的幾何伸長量

(6)

1.2.2 物理伸長量

管道的物理伸長量由軸向變形積分得到，軸向變形主要由四段管道的變形產(chǎn)生，A點左側(cè)管道，AB段管道，BC段管道，C點右側(cè)段管道。其中AB段與BC段管道內(nèi)軸向力在B點處于B點處集中力形成一個平衡，受力關(guān)系如圖2所示。

圖2 B點受力平衡示意圖

首先假設(shè)AB段軸力大小為F0，從而可以得到B點管道軸向分量為：FB=RBsin(α)=RBsin(atan(yB′))，因此BC段軸力FOB=F0+FB。

管道物理伸長量計算方法如式(7)：

(7)

根據(jù)兩者協(xié)調(diào)條件可以得到管道的軸向力大小F0，由于管道受拉所以舍棄負(fù)值，從而可以得到：

F0=

1.3 管道軸向應(yīng)力

管道內(nèi)的軸向應(yīng)力包括，軸力引起的應(yīng)力分量σ0=F0/A；溫度與內(nèi)壓引起的管道的軸向應(yīng)力σserive=-Eα(T2-T1)+vσh；彎矩引起的軸向應(yīng)力σbend=MR/I。三者求和可以得到管道的軸向應(yīng)力的總大小：σaxial=σservice+σbend+σ0。

2 模型準(zhǔn)確性驗證

為了驗證本文方法的準(zhǔn)確性，這里采用實際工程參數(shù)進行計算，并將計算結(jié)果與有限元模型結(jié)果進行對比。工況參數(shù)如下：管徑500 mm，壁厚6 mm，摩擦系數(shù)0.5，熱脹系數(shù)1.2e-5，單位長度管道重量720 N/m，管溝深度1 m。

有限元模型使用ABAQUS軟件建立，使用PIPE21單元模擬管道，使用離散剛體模擬海床，使用通用接觸定義兩者接觸關(guān)系，所有參數(shù)與采用上面給出的工況參數(shù)，計算分為兩個載荷步，第一個載荷步為工作荷載施加載荷步，第二個載荷步為管溝沖溝載荷步。在第一個載荷步施加內(nèi)壓，溫度等載荷，第二個載荷步中將右側(cè)土塊下移模擬沖溝產(chǎn)生的溝底深度。計算結(jié)果如圖3所示。

圖3 挖溝條件下管道響應(yīng)數(shù)值計算模型

圖4給出了本文解析結(jié)果與有限元結(jié)果的對比圖。結(jié)果表明，管道在B點處應(yīng)力最大，B點管頂軸向應(yīng)力為軸向拉應(yīng)力，管底為軸向壓應(yīng)力。本文方法與數(shù)值計算方法吻合很好，從工程角度，兩者誤差幾乎可以忽略。從而表明了本文給出方法的準(zhǔn)確性和合理性。

圖4 變形段位移與應(yīng)力結(jié)果

3 結(jié) 論

本文基于彈性梁理論給出了一種海底管道挖溝過程中應(yīng)力與變形求解的解析分析方法，方法基于彈性梁理論，得到了管道的撓曲線微分方程，基于邊界條件得到了撓曲線方程的解析解，并由此得到了管道內(nèi)彎曲應(yīng)力。在此基礎(chǔ)上，通過管道幾何伸長量和物理伸長量的協(xié)調(diào)條件，得到了管道內(nèi)的軸向力應(yīng)力。最終通過彎曲應(yīng)力與軸向應(yīng)力疊加得到管道內(nèi)的軸向總應(yīng)力。通過與有限元數(shù)值模型對比得到本解析方法結(jié)果與數(shù)值方法結(jié)果基本一致，可以給工程實際提供一定的參考。

[1] 方維章. 國外海底管線的鋪設(shè)和開溝埋管技術(shù)[J]. 海洋工程, 1983(3):82-88.

[2] 孟 然. 海底管道后挖溝技術(shù)實驗研究及數(shù)值模擬[D]. 天津大學(xué), 2012.

[3] PAULIN M，HUMBY D，COCKER J，et al． Trenching of Pipelines for Protection in Ice Environments [C]. SPE Arctic and Extreme Environments Conference and Exhibition, SPE, 2011: 1-7.

[4] 秦衛(wèi)華. 海底管道挖溝技術(shù)的改進[J]. 江蘇船舶, 2010, 27(1):18-20.

[5] 帥 健. 海底管線開溝埋設(shè)時的應(yīng)力分析[J]. 石油機械, 1997(8):39-40.

[6] 劉秉德, 倪玲英, 詹燕民,等. 海底管道挖溝的數(shù)值模擬[J]. 管道技術(shù)與設(shè)備, 2014(6):5-7.

AnalyticalStressAnalysisMethodforSubmarinePipelineduringTrenching

XIAMengying1,ZHANGHong1,WUKai1,WUJingtian2

(1.CollegeofMechanicalandTransportationEngineering,ChinaUniversityofPetroleum-Beijing,Beijing102249,China; 2.SINOPECTianjinLiquefiedNaturalGasCo.Ltd.,Tianjin, 300457,China)

During the after trenching process, large deformation will occur in the pipe, which can induce large stress leading to the potential failure. Thus an accurate stress estimation method for pipelines under this load condition was necessary. Based on elastic beam theory, the pipe deflection equations were proposed. The analytical stress and deformation results were derived by the boundary conditions and the coordination condition between the geometrical elongation and the physical one. Finally, this proposed method was validated accuracy for calculating pipe stress by comparing with a finite element model.

submarine pipeline； trenching lay； stress； analytical method； numerical validation

夏夢瑩，女，1990年生，在讀博士研究生，主要從事油氣裝備失效分析與安全評價。E-mail：xiamengying322@163.com

TE 973

A

2096-0077(2017)06-0041-03

10.19459/j.cnki.61-1500/te.2017.06.011

2017-04-20

馬小芳)