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

， ， ，經(jīng)天

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

·試驗(yàn)研究·

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

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

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

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

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

0 引 言

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

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

1 管道力學(xué)模型

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

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

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

1.1 彎曲變形分析

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

(1)

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

對(duì)于BC段，彎曲微分方程為

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

(2)

在B點(diǎn)處轉(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)

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

同時(shí)根據(jù)對(duì)B點(diǎn)求矩的平衡可以得到：

1.2 管道軸向變形

1.2.1 幾何伸長(zhǎng)量

管道的幾何伸長(zhǎng)量主要來(lái)自于管道變形后曲線長(zhǎng)度的增加，可以通過(guò)弧長(zhǎng)積分法得到變形后的曲線長(zhǎng)度，從而得到管道的幾何伸長(zhǎng)量

(6)

1.2.2 物理伸長(zhǎng)量

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

圖2 B點(diǎn)受力平衡示意圖

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

管道物理伸長(zhǎng)量計(jì)算方法如式(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)力的總大?。害襛xial=σservice+σbend+σ0。

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

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

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

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

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

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

3 結(jié) 論

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

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

[2] 孟 然. 海底管道后挖溝技術(shù)實(shí)驗(yàn)研究及數(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ìn)[J]. 江蘇船舶, 2010, 27(1):18-20.

[5] 帥 健. 海底管線開(kāi)溝埋設(shè)時(shí)的應(yīng)力分析[J]. 石油機(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

夏夢(mèng)瑩，女，1990年生，在讀博士研究生，主要從事油氣裝備失效分析與安全評(píng)價(jià)。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

馬小芳)