海底管道在剛性扣入式止屈器中的屈曲傳播理論研究

馬維林，余建星，張?chǎng)?/p>

?

海底管道在剛性扣入式止屈器中的屈曲傳播理論研究

馬維林1，2，余建星3，張?chǎng)?

(1.中國(guó)海洋大學(xué) 工程學(xué)院，山東 青島 266100；2.大連理工大學(xué) 海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024；3.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072)

為了從理論上研究海底管道在扣入式止屈器中的屈曲傳播問(wèn)題，建立了二維管道環(huán)模型，并運(yùn)用非線性離散彈簧模型來(lái)模擬管道內(nèi)壁以及管道外壁與止屈器內(nèi)壁之間的接觸作用。理論模型中考慮了結(jié)構(gòu)的非線性幾何大變形和彈塑性本構(gòu)關(guān)系，利用虛功原理建立了平衡方程，計(jì)算研究了海底管道在剛性扣入式止屈器中屈曲傳播時(shí)的外壓-截面積變化規(guī)律，然后基于Maxwell能量平衡法確定了扣入式止屈器的穿越壓力下限。通過(guò)將理論計(jì)算結(jié)果與全比例管道試驗(yàn)和有限元數(shù)值模擬進(jìn)行對(duì)比，證明了該理論方法的正確性，本理論模型為進(jìn)一步研究三維管道模型的屈曲傳播和止屈機(jī)理奠定了基礎(chǔ)。

海底管道；屈曲傳播；接觸作用；扣入式止屈器；穿越壓力

海底管道在海洋油氣開發(fā)中的結(jié)構(gòu)安全問(wèn)題不容忽視。隨著水深增加，靜水壓力也不斷升高，當(dāng)管道存在初始缺陷時(shí)，容易發(fā)生局部屈曲壓潰進(jìn)而發(fā)生屈曲傳播，從而造成大范圍的管道結(jié)構(gòu)破壞，為減小事故損失，通常在管道上安裝止屈器(如扣入式止屈器等)以阻止大范圍的屈曲傳播。

國(guó)內(nèi)外學(xué)者對(duì)海底管道屈曲傳播和止屈問(wèn)題開展了大量研究工作。在管道屈曲傳播研究方面，Dyau和Kyriakides[1]采用三維圓柱殼模型，對(duì)管道屈曲傳播的發(fā)生和發(fā)展過(guò)程進(jìn)行了模擬，不同徑厚比管道的屈曲傳播壓力計(jì)算值與試驗(yàn)值吻合較好；Kyriakides和Netto[2]運(yùn)用ABAQUS有限元軟件對(duì)管道的準(zhǔn)靜態(tài)和動(dòng)態(tài)屈曲傳播進(jìn)行了數(shù)值模擬研究；Xue[3]運(yùn)用ABAQUS有限元軟件分析了非均勻和受腐蝕海底管道在外壓作用下的局部屈曲和屈曲傳播特性。在止屈器研究方面，大多數(shù)學(xué)者采用有限元數(shù)值模擬以及試驗(yàn)手段來(lái)對(duì)帶止屈器管道的屈曲行為進(jìn)行研究。Kyriakides[4]運(yùn)用縮尺比試驗(yàn)方法研究了扣入式止屈器的止屈效率上下限問(wèn)題；Lee和Kyriakides[5,7]運(yùn)用ABAQUS有限元軟件對(duì)扣入式止屈器的止屈效率問(wèn)題進(jìn)行了模擬分析；Kyriakides和Lee[6-7]運(yùn)用ABAQUS有限元軟件對(duì)圍腔內(nèi)受限管道的屈曲傳播進(jìn)行了準(zhǔn)靜態(tài)模擬，并與試驗(yàn)結(jié)果進(jìn)行了對(duì)比?？梢钥闯?，國(guó)內(nèi)外學(xué)者對(duì)管道屈曲傳播和止屈問(wèn)題的理論研究工作并不多見(jiàn)，且國(guó)內(nèi)研究與國(guó)外研究還存在較大差距。

安裝扣入式止屈器管道屈曲分析的主要難點(diǎn)在于其同時(shí)涉及了三種非線性問(wèn)題，即：管道大變形的幾何非線性、管壁應(yīng)變較大區(qū)域的材料塑性行為非線性以及管壁與扣入式止屈器之間的接觸非線性。在加載過(guò)程中，三種非線性共同影響管道的變形行為，從而增加了問(wèn)題分析和求解的復(fù)雜性。本文對(duì)安裝扣入式止屈器管道的屈曲傳播問(wèn)題進(jìn)行理論研究，理論結(jié)果與ABAQUS有限元軟件的數(shù)值模擬結(jié)果進(jìn)行了對(duì)比；之后基于Maxwell能量平衡法，求出了扣入式止屈器穿越壓力的理論下限，并與全比例管道止屈試驗(yàn)結(jié)果進(jìn)行了對(duì)比分析。

1 管道屈曲傳播力學(xué)模型

1.1 管道模型

為了采用二維環(huán)模型來(lái)作為管道的力學(xué)模型，在三維管道的基礎(chǔ)上作出以下基本假定來(lái)進(jìn)行模型簡(jiǎn)化[8-9]：

1)管壁環(huán)向應(yīng)變沿管壁厚度方向滿足直法線假定條件；

2)管壁的徑向應(yīng)變忽略不計(jì)；

3)忽略扣入式止屈器與管道之間以及管道自身內(nèi)壁之間的相對(duì)切向運(yùn)動(dòng)摩擦；

4)管道變形過(guò)程中，材料保持各向同性；

5)僅考慮管道初始橢圓度幾何缺陷；

6)加載過(guò)程為準(zhǔn)靜態(tài)過(guò)程。

具有初始橢圓度幾何缺陷的二維管道環(huán)模型如圖1所示。

圖1 管道二維環(huán)模型坐標(biāo)系Fig.1 Coordinates of the 2D ring model

以管道中心為原點(diǎn)建立整體坐標(biāo)系，et、exp t分別為水平軸和豎直軸，用以描述管道的實(shí)際空間形狀；在管壁中性層上任意點(diǎn)處建立局部坐標(biāo)系，v、w分別表示管壁單元的環(huán)向和徑向位移，用以分析管道應(yīng)變。R是未發(fā)生變形管道的中性層半徑，中性層上任一點(diǎn)的位置可以用R和θ來(lái)表示；P表示管壁的內(nèi)外壓差。兩個(gè)坐標(biāo)系的關(guān)系表示為

(1)

(2)

1.2 接觸模型

假定扣入式止屈器為完全剛性，其截面形狀不發(fā)生變化。運(yùn)用非線性離散彈簧模型來(lái)模擬管道自身內(nèi)壁之間以及管道外壁與扣入式止屈器內(nèi)壁之間的接觸邊界，在管道內(nèi)外壁發(fā)生接觸的位置建立一系列離散的非線性彈簧模型，如圖2所示。彈簧接觸力與接觸距離之間呈非線性關(guān)系。

圖2 管壁上的接觸彈簧模型Fig.2 Spring models on pipeline′s inner and outer walls

對(duì)于管道內(nèi)壁的彈簧模型，設(shè)定與離散彈簧對(duì)應(yīng)的距離量為l，用以表征接觸面之間的接觸情況。計(jì)算加載過(guò)程中各內(nèi)彈簧所對(duì)應(yīng)的距離量l，中性層與相鄰界面距離為g0(等于壁厚t的一半)。當(dāng)l首次達(dá)到g0時(shí)，認(rèn)為該點(diǎn)開始發(fā)生接觸，而其余l(xiāng)>g0處的未接觸彈簧則暫時(shí)不予考慮。對(duì)于發(fā)生接觸的內(nèi)彈簧模型：當(dāng)lg0時(shí)，彈簧力為負(fù)，阻止已接觸管壁發(fā)生分離[10]。內(nèi)彈簧的彈簧力F與接觸距離l之間的關(guān)系式為

(3)

式中:F0為內(nèi)彈簧剛度系數(shù)，β為內(nèi)彈簧反作用力的增長(zhǎng)率。

(4)

1.3 基本方程

求解此二維環(huán)模型的屈曲問(wèn)題需要建立三類基本方程，即幾何方程、本構(gòu)方程和平衡方程，從而分別建立位移與應(yīng)變、應(yīng)力與應(yīng)變、位移與外力之間的關(guān)系。

1.3.1 幾何方程

利用Sanders[11]圓柱殼非線性大變形理論和基爾霍夫-勒夫假定，環(huán)模型環(huán)向方向的總應(yīng)變?yōu)?/p>

(5)

式中:Eθ為環(huán)向薄膜應(yīng)變，Kθ為環(huán)向彎曲應(yīng)變，z為沿中性層法線方向坐標(biāo)?？紤]管道橫截面大變形中的高階項(xiàng)，環(huán)向的薄膜應(yīng)變與彎曲應(yīng)變分別為

(6)

(7)

1.3.2 本構(gòu)方程

考慮管道屈曲過(guò)程中的大變形，管道會(huì)進(jìn)入彈塑性變形階段，采用J2塑性流動(dòng)理論來(lái)建立二維管道的環(huán)向和軸向應(yīng)力-應(yīng)變關(guān)系[12]，即

(8)

式中:

1.3.3 能量平衡方程

(9)

準(zhǔn)靜態(tài)加載過(guò)程中的外力功增量為[9]

(10)

(11)

管道環(huán)模型的應(yīng)變能增量為

(12)

內(nèi)彈簧應(yīng)變能增量的變分為所有已發(fā)生接觸內(nèi)彈簧應(yīng)變能之和的變分：

(13)

外彈簧應(yīng)變能增量的變分為所有已發(fā)生接觸外彈簧應(yīng)變能之和的變分：

(14)

2 理論結(jié)果與有限元數(shù)值模擬

2.1 理論結(jié)果

以管材API 5L X65為例，彈性模量為E=2.06×105MPa，屈服極限為σ0=448 MPa，分別計(jì)算了具有初始橢圓度的325 mm×10 mm、406 mm×10 mm兩種徑厚比管道的屈曲變形過(guò)程。

管道自身接觸發(fā)生之前以及管道自身接觸發(fā)生之后的管道中性層變形情況如圖3所示。

圖3 管道中性層變形情況Fig.3 Deformation of neutral layer

從圖3中可以看出，在管道自身接觸發(fā)生之前，管道截面按預(yù)想的變形方式在扣入式止屈器內(nèi)產(chǎn)生“狗骨式”變形；當(dāng)管道截面變形至自身上下內(nèi)壁發(fā)生接觸后，管道水平軸上的離散彈簧模型開始啟用。

2.2 有限元模擬對(duì)比

采用與上述理論計(jì)算相同的模型參數(shù)，運(yùn)用ABAQUS軟件進(jìn)行有限元數(shù)值模擬，如圖4所示。

圖4 管道和止屈器的有限元模型Fig.4 Finite element model of the pipeline and arrestor

管道外部運(yùn)用二維殼模型建立四分之一圓環(huán)來(lái)模擬扣入式止屈器，止屈器材料定義為彈性模量很大的剛性體；管道和止屈器網(wǎng)格劃分均采用非協(xié)調(diào)應(yīng)力單元(CPS4I)；模型采用流體加載，并在管道外部采用二維線模型建立與管道同心四分之一圓作為流體加載的外殼；模型運(yùn)用二維剛性面在管道水平軸上建立長(zhǎng)度小于管道內(nèi)徑的剛性面，用以模擬管道自身內(nèi)壁的接觸平面。

理論計(jì)算和有限元模擬的管道外壓變化曲線對(duì)比如圖5所示。

對(duì)于325 mm×10 mm、406 mm×10 mm兩種徑厚比的管道，理論計(jì)算和有限元模擬所得的外壓變化曲線在管道內(nèi)壁發(fā)生接觸后的上升段有較大差別，這可能是由于有限元單元的定義、劃分以及有限元模擬和理論計(jì)算接觸邊界的模型不同所致，但外壓的整體變化趨勢(shì)比較吻合，因此可以說(shuō)明該理論計(jì)算方法是正確的。

圖5 理論計(jì)算和有限元模擬的管道外壓變化曲線Fig.5 External pressure curves of the theoretical calculation and finite element simulation

3 扣入式止屈器的穿越壓力

3.1 穿越壓力的理論計(jì)算

根據(jù)圍腔內(nèi)受限管道的屈曲傳播分析研究[6]，其外壓-時(shí)間的準(zhǔn)靜態(tài)變化曲線如圖6所示。

圖6 部分圍腔內(nèi)受限管道的外壓-時(shí)間準(zhǔn)靜態(tài)變化曲線Fig.6 External pressure-time history of a quasi-static test on a partially confined pipeline

圖6中PI表示管道的屈曲壓潰壓力，PP表示管道的屈曲傳播壓力，PIC可以認(rèn)為是圍腔內(nèi)受限管道的屈曲壓潰壓力，PPC表示圍腔內(nèi)受限管道的屈曲傳播壓力，也可以看成管道屈曲“穿越”長(zhǎng)圍腔的最小壓力。如果將長(zhǎng)圍腔縮短至扣入式止屈器的長(zhǎng)度，則PPC便是該扣入式止屈器的穿越壓力下限，而PIC則是穿越壓力上限。

圖7為管道外壓隨其截面變形的曲線(P-δυ)，認(rèn)為管道在外壓PPC作用下進(jìn)行準(zhǔn)靜態(tài)屈曲傳播，單位長(zhǎng)度管道發(fā)生屈曲傳播的外力功為

式中:A是管道未發(fā)生變形時(shí)的初始狀態(tài)，C是管道發(fā)生屈曲傳播后的最終狀態(tài)。

假定材料行為與路徑無(wú)關(guān)，此時(shí)管道應(yīng)變能變化便只與管道截面的初始狀態(tài)A和最終狀態(tài)C有關(guān)，而應(yīng)變能變化量應(yīng)等于外力功，即

(15)

對(duì)于彈性變形，材料行為與路徑無(wú)關(guān)，式(15)嚴(yán)格滿足；但對(duì)于塑性變形，材料行為與路徑相關(guān)，式(15)并不嚴(yán)格滿足，盡管如此，這里仍用該方法來(lái)進(jìn)行近似求解。

基于本文帶扣入式止屈器二維管道理論模型所求出的外壓-截面積變化曲線，可以運(yùn)用Maxwell作圖法[6-7]來(lái)確定扣入式止屈器的穿越壓力下限PPC。為滿足式(15)，可以將PPC的水平線畫在P-δυ圖中，保證面積A1等于面積A2即可。

圖7 P-δυ變化曲線和Maxwell作圖法Fig.7 Curve of P-δυ and Maxwell′s construction

3.2 扣入式止屈器的止屈試驗(yàn)

為了驗(yàn)證理論分析和有限元模擬結(jié)果的可靠性，筆者在天津大學(xué)深海壓力艙試驗(yàn)室[13-14]進(jìn)行了安裝扣入式止屈器管道的屈曲試驗(yàn)。

扣入式止屈器有多種形式，考慮到制作及試驗(yàn)安裝的便捷性，試驗(yàn)中采用螺栓夾緊型止屈器，如圖8所示。

圖8 螺栓夾緊型止屈器剖面Fig.8 Cross section of clamped arrestors

1)制作螺栓夾緊型止屈器，并將其安裝在管件中心位置，之后在管件兩端安裝固定密封法蘭，并在管件上每隔一定間距沿環(huán)向粘貼深海應(yīng)變片，然后將管件放入壓力艙內(nèi)部并固定。

2)壓力艙密封后，將艙體內(nèi)部完全注滿水，然后通過(guò)加壓系統(tǒng)進(jìn)行加壓，采集記錄水壓數(shù)據(jù)(如圖9所示)及變形數(shù)據(jù)。

3)當(dāng)艙體內(nèi)部水壓達(dá)到一定值時(shí)，由于管件左端存在初始橢圓度缺陷，因此管件左端首先發(fā)生局部屈曲，此為試驗(yàn)第1)階段：局部屈曲壓潰階段。

4)由于壓力艙內(nèi)部密封，管件發(fā)生局部屈曲后體積驟然減小，艙內(nèi)水容積增大，導(dǎo)致水壓瞬間大幅下降，然后繼續(xù)加壓，當(dāng)水壓上升到某一值時(shí)，便會(huì)發(fā)生屈曲傳播，此時(shí)艙內(nèi)水壓基本維持在一個(gè)穩(wěn)定值，此為試驗(yàn)第2)階段：止屈器前屈曲傳播階段。

5)屈曲傳播遇到止屈器后便被抑制，當(dāng)艙內(nèi)水壓繼續(xù)上升至一定水壓值后，屈曲傳播便會(huì)穿越止屈器繼續(xù)傳播，此時(shí)所對(duì)應(yīng)的水壓值即為該止屈器的穿越壓力，此為試驗(yàn)第3)階段：屈曲傳播穿越止屈器階段。該階段中管道外壁與扣入式止屈器內(nèi)壁間的接觸如圖10所示。

6)在屈曲傳播穿越止屈器后，水壓會(huì)有小幅下降而后上升至穩(wěn)定值，此為試驗(yàn)第4)階段：止屈器后屈曲傳播階段。

圖9 扣入式止屈器試驗(yàn)中的水壓變化曲線Fig.9 Variation of water pressure in the slip-on arrestor′s experiment

圖10 管道外壁與扣入式止屈器內(nèi)壁間的接觸Fig.10 Contact between pipeline′s outer wall and slip-on arrestor′s inner wall

3.3 理論結(jié)果與試驗(yàn)結(jié)果對(duì)比分析

試驗(yàn)測(cè)出的扣入式止屈器穿越壓力值和3.1節(jié)中的穿越壓力理論值對(duì)比如表1所示?？梢钥闯?，對(duì)于同一徑厚比的管道來(lái)說(shuō)，扣入式止屈器的厚度比長(zhǎng)度對(duì)穿越壓力的影響更為顯著，止屈器越厚，穿越壓力越大；另外從文獻(xiàn)[6]可知，當(dāng)扣入式止屈器較長(zhǎng)時(shí)，管道會(huì)在PPC外壓下在止屈器內(nèi)部發(fā)生屈曲傳播而穿越止屈器。

由于本文理論采用的是帶扣入式止屈器的二維管道模型，因此，止屈器長(zhǎng)度可以看作無(wú)限長(zhǎng)，所求出的理論值只是穿越壓力的下限值PPC。

表1 扣入式止屈器穿越壓力試驗(yàn)值和理論值

4 結(jié)論

1)通過(guò)理論計(jì)算和有限元數(shù)值模擬對(duì)比，說(shuō)明本文的二維管道環(huán)模型和接觸彈簧模型可用于研究帶扣入式止屈器海底管道在外部靜水壓力作用下的截面變形及外壓變化情況。

2)由于假定扣入式止屈器為完全剛性，因此，所求出的穿越壓力理論值更適用于厚度較大的止屈器；接下來(lái)應(yīng)進(jìn)一步研究止屈器厚度對(duì)止屈效果的影響機(jī)制。

3)運(yùn)用Maxwell能量平衡法求得的穿越壓力下限值，由于未考慮材料塑性行為的路徑相關(guān)性，所求出的穿越壓力理論值略小于試驗(yàn)值。

[1]DYAU J Y, KYRIAKIDES S. On the propagation pressure of long cylindrical shells under external pressure[J]. International journal of mechanical sciences, 1993, 35(8): 675-713.

[2]KYRIAKIDES S, NETTO T A. On the dynamics of propagating buckles in pipelines[J]. International journal of solids and structures, 2000, 37(46/47): 6843-6867.

[3]XUE Jianghong. A non-linear finite-element analysis of buckle propagation in subsea corroded pipelines[J]. Finite elements in analysis and design, 2006, 42(14/15): 1211-1219.

[4]KYRIAKIDES S. Efficiency limits for slip-on type buckle arrestors for offshore pipelines[J]. Journal of engineering mechanics, 2002, 128(1): 102-111.

[5]LEE L H, KYRIAKIDES S. On the arresting efficiency of slip-on buckle arrestors for offshore pipelines[J]. International journal of mechanical sciences, 2004, 46(7): 1035-1055.

[6]KYRIAKIDES S, LEE L H. Buckle propagation in confined steel tubes[J]. International journal of mechanical sciences, 2005, 47(4/5): 603-620.

[7]LEE L H. On the design of slip-on buckle arrestors for offshore pipelines[D]. Austin: The University of Texas at Austin, 2007.

[8]余建星, 李智博, 杜尊峰, 等. 深海管道非線性屈曲理論計(jì)算方法[J]. 海洋工程, 2013, 31(1): 54-60.

YU Jianxing, LI Zhibo, DU Zunfeng, et al. Theoretical calculation method of the nolinear buckling of deepsea pipes[J]. The ocean engineering, 2013, 31(1): 54-60.

[9]楊源. 海底管道屈曲破壞與止屈器研究[D]. 天津: 天津大學(xué), 2013. YANG Yuan. The research of submarine pipeline buckling failure and buckle arrestor[D]. Tianjin: Tianjin University, 2013.

[10]DYAU J Y. Localization and propagation of instabilities in long cylindrical shells under external pressure[D]. Austin: The University of Texas at Austin, 1993.

[11]SANDERS J L JR. Nonlinear theories for thin shells[J]. Quarterly of applied mathematics, 1963, 21(1): 21-36.

[12]黃克智, 黃永剛. 固體本構(gòu)關(guān)系[M]. 北京: 清華大學(xué)出版社, 1999. HUANG Kezhi, HUANG Yonggang. Constitutive relationship of solid[M]. Beijing: Tsinghua University Press, 1999.

[13]卞雪航, 吳海欣, 余楊, 等. 深水海底管道屈曲傳播試驗(yàn)設(shè)計(jì)[J]. 海洋技術(shù), 2011, 30(4): 93-95, 131.

BIAN Xuehang, WU Haixin, YU Yang, et al. Scheme design for deep water pipeline buckling propagation experiment[J]. Ocean technology, 2011, 30(4): 93-95, 131.

[14]余建星, 卞雪航, 余楊, 等. 深水海底管道全尺寸壓潰試驗(yàn)及數(shù)值模擬[J]. 天津大學(xué)學(xué)報(bào), 2012, 45(2): 154-159. YU Jianxing, BIAN Xuehang, YU Yang, et al. Full-scale collapse test and numerical simulation of deepwater pipeline[J]. Journal of Tianjin University, 2012, 45(2): 154-159.

Theoretical study of buckling propagation of subsea pipelines with rigid slip-on arrestors

MA Weilin1, 2, YU Jianxing3, ZHANG Xin1

(1. College of Engineering, Ocean University of China, Qingdao, Shandong 266100, China; 2. State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian, Liaoning 116024, China; 3. State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China)

In order to theoretically study buckling propagation of subsea pipelines with slip-on arrestors, a two-dimensional ring model was used to represent the pipeline and a nonlinear discrete spring model was adopted to simulate the contact between a pipeline′s inner walls and between a pipeline′s outer wall and an arrestor′s inner wall. The theoretical model considered a large geometric deformation and the elastic-plastic constitutive relation of material, a balance equation was established on the basis of the principle of virtual work. The variation of external pressure with respect to the pipeline′s cross-sectional area was analyzed, and the lower bound of the crossing pressure of the slip-on arrestor was calculated on the basis of Maxwell′s energy balance method. By comparing the theoretical results with the full-scale pipeline test and finite element numerical simulation, the theoretical method proved to be correct and reliable. The theoretical model presented in this paper lays the foundation for further research on buckling propagation and the arresting mechanism of three-dimensional models of pipelines.

subsea pipeline; buckling propagation; contact effect; slip-on arrestor; crossover pressure

2015-08-23.

日期：2016-08-29.

國(guó)家自然科學(xué)基金項(xiàng)目(51609222)；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)(201513039)；大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(LP1505)；天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(HESS-1602).

馬維林(1986- )，男，講師，博士后.

馬維林，E-mail：maweilinfreedom@163.com.

10.11990/jheu.201508043

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20160829.1421.028.html

P756.2

A

1006-7043(2016)10-1353-06

馬維林，余建星，張?chǎng)? 海底管道在剛性扣入式止屈器中的屈曲傳播理論研究[J]. 哈爾濱工程大學(xué)學(xué)報(bào)， 2016, 37(10): 1353-1358.

MA Weilin, YU Jianxing, ZHANG Xin. Theoretical study of buckling propagation of subsea pipelines with rigid slip-on arrestors[J]. Journal of Harbin Engineering University, 2016, 37(10): 1353-1358.