海底橢圓管道在扭矩和水壓作用下的屈曲壓潰研究

余建星 ，安思宇 ，李志剛，段晶輝 ，孫震洲 ，蔣梅榮 ，韓夢(mèng)雪

(1. 天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072；2. 上海交通大學(xué)高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240；3. 中海油研究總院，北京 100028)

近年來，我國(guó)石油產(chǎn)量穩(wěn)步提升，而增量主要來自海洋石油開發(fā).油氣輸送的主渠道是深水海底管道，海底管道保證了整個(gè)系統(tǒng)的正常運(yùn)行，其穩(wěn)定性和安全性對(duì)于海洋油氣開發(fā)極為重要[1].海底管道在外部水壓的作用下有可能發(fā)生屈曲壓潰[2].壓潰后的管道極易發(fā)生屈曲傳播，這是由于海底管道壓潰壓力遠(yuǎn)大于其屈曲傳播壓力，能夠使屈曲沿管道軸向進(jìn)行傳播的外部水壓通常只有管道壓潰壓力的 18%～25%[3].因此，管道一旦發(fā)生屈曲壓潰，就會(huì)造成管道結(jié)構(gòu)的迅速破壞，導(dǎo)致重大的財(cái)產(chǎn)損失和海洋環(huán)境污染.

海底管道的壓潰壓力體現(xiàn)的是其極限承載力，而管道局部缺陷或損傷以及復(fù)雜載荷的作用都會(huì)對(duì)管道的極限承載力產(chǎn)生影響.在工程實(shí)際中，由于加工工藝水平限制和海洋環(huán)境載荷作用，管道截面通常不是完美的圓環(huán)，而是具有一定橢圓度的橢圓環(huán).因此，研究具有橢圓度缺陷的管道在扭矩和水壓作用下的屈曲破壞具有重要的意義.管道在海底服役前的加工和鋪管過程可能會(huì)受到不同程度扭矩的作用，管道發(fā)生扭轉(zhuǎn)變形[4].

一直以來，國(guó)內(nèi)外學(xué)者對(duì)于受復(fù)雜載荷管道的研究主要集中在受軸力和彎矩的管道屈曲特性[5-7]上.對(duì)于受扭矩和水壓作用的橢圓管道研究較少，Wang[8]、Arghavan等[9]對(duì)受扭矩作用的空心截面管進(jìn)行了理論研究，但沒有考慮外部水壓.Sharma等[10]對(duì)受到扭矩和外壓作用下的圓環(huán)截面管道進(jìn)行了理論計(jì)算，但沒有考慮管道的橢圓度缺陷.Zheleznov等[11]對(duì)受到扭矩和內(nèi)壓作用下的空心橢圓截面殼進(jìn)行了有限元分析，但殼的厚度很小，結(jié)論不能應(yīng)用于具有較大壁厚的海底管道.結(jié)合國(guó)內(nèi)外學(xué)者的研究來看，還沒有學(xué)者針對(duì)橢圓管道在扭矩和水壓作用下的屈曲破壞進(jìn)行深入的研究.

管道在服役前后的受力過程會(huì)經(jīng)歷兩個(gè)階段，第1個(gè)階段是受到較小扭矩作用的彈性變形階段，第 2個(gè)階段是受到較大外部水壓作用的塑性變形階段.本文基于 Timoshenko的彈性理論[12]，對(duì)橢圓管道在彈性階段受扭矩作用的位移方程進(jìn)行了求解.對(duì)于橢圓管道在塑性階段的屈曲壓潰，通過建立ABAQUS有限元仿真模型，利用弧長(zhǎng)法對(duì)受到扭矩和外壓作用下橢圓管道的屈曲破壞進(jìn)行數(shù)值模擬，研究管道壁厚、一致橢圓度以及外徑對(duì)于管道屈曲壓潰的影響.管道在服役過程中也有可能受到扭矩作用，故改變扭矩-水壓的加載路徑，探究管道在海底服役中受到扭矩和水壓作用下的屈曲壓潰.

1 彈性階段

1.1 海底橢圓管道彈性扭轉(zhuǎn)方程

在工程實(shí)際中，由于加工工藝水平限制、鋪管過程操作以及海底環(huán)境載荷作用，管道截面通常不是完美的圓環(huán)，而是具有一定橢圓度的等壁厚橢圓環(huán).對(duì)于較長(zhǎng)橢圓管道而言，這種橢圓度變形沿軸向幾乎是一致的，可稱為一致橢圓度.

管道在服役前可能受到扭矩作用.考慮到在實(shí)際工程中，管道受到的扭矩通常較小，此時(shí)管道仍處于彈性變形階段或剛進(jìn)入塑性階段.因此可將橢圓管道簡(jiǎn)化為圖2所示等截面空心直桿，對(duì)其在受到扭矩作用下的自由彈性扭轉(zhuǎn)進(jìn)行分析.

圖2 空心橢圓截面桿件Fig.2 Hollow member bar with oval cross section

圖2 所示空心橢圓截面桿件，其外部橢圓長(zhǎng)軸為a，短軸為b，橢圓截面等壁厚，壁厚為t.

管道截面橢圓度的計(jì)算式為

式中：maxD 為橢圓最大直徑；minD 為橢圓最小直徑.

對(duì)于外部橢圓而言，其應(yīng)力函數(shù)1φ可取

式中m為待定常數(shù).

橢圓的相容方程為

式中C為常數(shù).由桿件端部邊界條件可得外部橢圓的扭矩[13]為

根據(jù) Timoshenko的彈性理論[12]，對(duì)于等壁厚的空心截面桿件，橢圓環(huán)中內(nèi)部橢圓可視作具有外部邊界實(shí)心橢圓軸的一條應(yīng)力線，線上任一點(diǎn)的剪應(yīng)力必然是沿著這條線的切向方向，這也就意味著移去內(nèi)部橢圓后，外部的應(yīng)力分布不受影響，應(yīng)力函數(shù)也可用于內(nèi)部橢圓.因此，對(duì)于一定的扭轉(zhuǎn)角，環(huán)形截面的扭矩等于實(shí)心部分的扭矩減去空心部分的扭矩.由此，可以得到內(nèi)部橢圓的扭矩為

由上可得空心橢圓桿件的扭矩為

故

將m代入式(3)，可得

因此，單位長(zhǎng)度桿件扭轉(zhuǎn)角為

位移分量為

式中u、v分別為結(jié)點(diǎn)水平位移和豎直位移.

式中w為結(jié)點(diǎn)軸向位移.對(duì)其分別積分，且不計(jì)剛體位移，可得

1.2 理論驗(yàn)證

通過大型商業(yè)有限元軟件ABAQUS對(duì)橢圓管道在扭矩作用下的準(zhǔn)靜態(tài)彈性變形進(jìn)行模擬來驗(yàn)證理論的正確性.建立圖 3所示外徑 325,mm、壁厚10,mm、管長(zhǎng) 6,m 且具有 0.5%一致橢圓度的橢圓管道.對(duì)其施加一定大小的扭矩后，將有限元軟件計(jì)算得到的位移結(jié)果與理論模型計(jì)算結(jié)果匯總于表1.

由表1可見，有限元模型的計(jì)算結(jié)果與理論計(jì)算的結(jié)果非常接近，可以認(rèn)為該理論模型較好反映了等壁厚橢圓管道在彈性階段受扭矩作用后的位移.

圖3 200,kN·m扭矩作用下橢圓管道變形Fig.3 Deformation of oval pipeline under torsion of 200,kN·m

表1 彈性階段理論模型驗(yàn)證Tab.1 Verification of theoretical model in elastic stage

2 塑性階段

上一節(jié)中對(duì)橢圓管道在服役前受扭矩作用的彈性位移方程進(jìn)行了求解.海底管道在深水服役中，會(huì)受到高達(dá)幾十 MPa的外部水壓作用，管道會(huì)產(chǎn)生塑性變形甚至發(fā)生屈曲，可使用 ABAQUS軟件對(duì)管道在塑性階段的變形和屈曲進(jìn)行數(shù)值模擬.

2.1 本構(gòu)關(guān)系

利用工程中得到廣泛應(yīng)用的Ramberg-Osgood本構(gòu)方程擬合得到材料的塑性應(yīng)力-應(yīng)變曲線.

Ramberg-Osgood 方程的原型[14]為

式中：ε為應(yīng)變；σ為應(yīng)力；K為強(qiáng)度系數(shù)；n為應(yīng)變硬化指數(shù)；E為材料的彈性模量.

引入關(guān)于材料屈服應(yīng)力yσ的參數(shù)α[15]為

可以得到應(yīng)力-應(yīng)變關(guān)系表達(dá)式為

對(duì)于API 5,L X65鋼，材料的屈服應(yīng)力yσ的值可取為399.6,MPa，n值可取為10.7[16].

2.2 弧長(zhǎng)法基本原理

在 ABAQUS軟件中，采用弧長(zhǎng)法進(jìn)行橢圓管道受扭矩和水壓作用下屈曲壓潰的數(shù)值模擬.弧長(zhǎng)法是一種雙重增量計(jì)算法，需要同時(shí)在計(jì)算中考慮荷載增量和位移增量.

非線性方程組通?？杀硎緸?/p>

式中：P(x)為函數(shù)矩陣；F為外載荷矩陣.式(17)可寫為

式中：K為切線剛度矩陣；x為與K相關(guān)的位移矩陣.考慮到非線性方程組在迭代過程中存在殘差，引入一個(gè)荷載因子λ和迭代結(jié)束后存在的殘差R[17]，則方程組可改寫為

該方程組共有 j個(gè)方程， 1j+ 個(gè)未知量λ)，需附加控制方程

其中，Δl在每個(gè)增量步內(nèi)為固定值，也就意味著，每個(gè)增量步中迭代路徑沿半徑為Δl的圓弧進(jìn)行.

2.3 有限元建模

2.3.1 管道模型

考慮到受扭矩和水壓作用的管道并不能保持軸對(duì)稱的結(jié)構(gòu)，選擇建立管道的全圓模型，可以更好地模擬真實(shí)情況下管道的屈曲失效過程.

建模中還需設(shè)置與管道模型端面相耦合的參考點(diǎn)，以實(shí)現(xiàn)通過參考點(diǎn)對(duì)管道整體模型施加扭矩的目的.采用需考慮權(quán)重影響的連續(xù)分布耦合方式，在參考點(diǎn)至從節(jié)點(diǎn)的自由度傳遞過程中，權(quán)重系數(shù)公式為

式中：ri表示從節(jié)點(diǎn) i到參考點(diǎn)的距離；r0表示該距離的初始值.

2.3.2 有限元結(jié)果校驗(yàn)

天津大學(xué)建造的深海全尺寸壓力試驗(yàn)艙(如圖 4所示)可對(duì)具有各種一致橢圓度的橢圓管道進(jìn)行全尺寸試驗(yàn)[18]，得到橢圓管道的壓潰壓力.

圖4 深水全尺寸壓力艙Fig.4 Deep-water full-scale hyperbaric chamber

天津大學(xué)深水結(jié)構(gòu)試驗(yàn)室完成了一批Q235B材質(zhì)的全尺寸鋼管屈曲試驗(yàn)，將試驗(yàn)得到的管道壓潰壓力結(jié)果、有限元模型計(jì)算結(jié)果與以DNV F101規(guī)范[19]計(jì)算得到的屈曲壓力匯總于表 2.325,mm×10,mm管道屈曲壓潰及屈曲傳播的全尺寸試驗(yàn)結(jié)果如圖 5所示.

表2 有限元弧長(zhǎng)法計(jì)算結(jié)果驗(yàn)證Tab.2 Results verification of arc-length method

圖5 325,mm×10,mm管道屈曲壓潰及其傳播全尺寸試驗(yàn)Fig.5 Full-scale test of buckling collapse and propagation of pipeline with the parameter of 325,mm×10,mm

根據(jù)表2所示的結(jié)果對(duì)比，可看出弧長(zhǎng)法計(jì)算的結(jié)果與試驗(yàn)值和規(guī)范解的誤差不大.說明通過有限元模擬方法對(duì)橢圓管道的極限承載力進(jìn)行計(jì)算是可行的，可以用此模型對(duì)受到扭矩作用后橢圓管道的極限承載力進(jìn)行研究.

2.4 計(jì)算結(jié)果

2.4.1 敏感性因素分析

在驗(yàn)證有限元模擬計(jì)算管道壓潰壓力準(zhǔn)確性的基礎(chǔ)上，對(duì)可能影響管道受扭矩和水壓作用下壓潰壓力大小的因素，如管道壁厚、一致橢圓度、外徑等，進(jìn)行研究.

扭矩的加載路徑為：先對(duì)管道施加扭矩(T)，之后施加均布外部水壓(P)，以模擬管道服役前及受到扭矩作用后在海底服役中受到水壓壓潰的工程實(shí)際.該路徑可稱為 T-P路徑.需要注意的是，由于需要保證管道在 T-P加載路徑下的扭轉(zhuǎn)變形處于彈性階段或剛進(jìn)入塑性階段，以避免管道產(chǎn)生扭轉(zhuǎn)屈曲破壞，影響服役中的管道在外部水壓作用下的壓潰壓力計(jì)算結(jié)果，所以在有限元模擬對(duì)管道施加的扭矩不宜過大.同時(shí)也應(yīng)該考慮在極端工況中可能出現(xiàn)的較大扭矩對(duì)管道屈曲壓潰的影響.因此扭矩的最大值可取為管道抗扭極限承載力的 80%.對(duì)于具有不同壁厚、不同一致橢圓度的管道，其抗扭極限承載力不同，應(yīng)對(duì)施加的扭矩大小進(jìn)行相應(yīng)的調(diào)整.

2.4.2 T-P路徑下管道壁厚、一致橢圓度對(duì)管道壓潰壓力的影響

計(jì)算采用外徑為325,mm的工程常見管道，管道材料選擇 API 5,L X65型鋼.設(shè)置管道壁厚分別為10,mm、15,mm、20,mm.對(duì)于每種不同壁厚的管道，分別設(shè)置1%、2%和5%,3種一致橢圓度.在保持管道外徑、施加扭矩大小和扭矩加載路徑相同的情況下，分別計(jì)算具有不同壁厚和不同一致橢圓度管道在T- P路徑下的壓潰壓力.計(jì)算結(jié)果如圖6所示.

由圖 6可見，對(duì)于具有相同壁厚的橢圓管道，其一致橢圓度越大，管道在 T-P路徑下的壓潰壓力越小.同時(shí)，對(duì)于具有相同壁厚的橢圓管道，施加的扭矩越大，管道的壓潰壓力通常越低，并且這種變化不是線性的，扭矩越大，壓潰壓力的下降速度越快.對(duì)于壁厚較大、一致橢圓度較小的橢圓管道，壓潰壓力的變化曲線會(huì)出現(xiàn)極值點(diǎn)，在極值點(diǎn)后管道的壓潰壓力呈上升趨勢(shì).這說明扭矩對(duì)管道的影響并非是始終消極的，在管道壁厚較大、一致橢圓度較小時(shí)，管道扭轉(zhuǎn)變形有可能增加管道的極限承載力.

另一方面，對(duì)于具有相同初始一致橢圓度的橢圓管道，其壁厚越大，在相同扭矩作用后壓潰壓力也越大，這說明壁厚越大的橢圓管道的極限承載力越強(qiáng).壁厚較小時(shí)壓潰壓力隨壁厚增大而增長(zhǎng)的速率較快，而壁厚較大時(shí)壓潰壓力的增長(zhǎng)速率較慢.與此同時(shí)，對(duì)于一致橢圓度較小的管道，其壓潰壓力隨扭矩增大而下降的速率越慢.對(duì)于一致橢圓度較小、壁厚較大的管道，仍然可以得到與上文相同的結(jié)果.

圖6 不同壁厚、不同一致橢圓度管道的壓潰壓力Fig.6 Collapse pressures of pipeline with different thicknesses and uniform ovalities

圖7 橢圓管道變形形式Fig.7 Forms of deformation of oval pipeline

在 T-P路徑下，管道的變形會(huì)有 3種形式.由圖7(a)可以看到，在扭矩較小時(shí)，具有一致橢圓度缺陷的管道的初始?jí)簼⑽恢貌皇鞘冀K位于模型的端面，其位置是不確定的.這是因?yàn)榕ぞ剡€沒有使管道端部產(chǎn)生足夠程度的變形，這種變形對(duì)于管道缺陷程度并沒有明顯影響，這就導(dǎo)致了管道的初始?jí)簼⑽恢檬遣淮_定的，可能會(huì)發(fā)生在管道沿軸向方向的任何一個(gè)位置.而當(dāng)扭矩逐漸增大，管道會(huì)在 2個(gè)位置同時(shí)發(fā)生壓潰，如圖 7(b)所示.這說明在扭矩的作用下，管道的端面發(fā)生了一定程度的變形，這種缺陷與管道本身存在的一致橢圓度缺陷相疊加后，與管道一致橢圓度缺陷相接近.如圖 7(c)所示，當(dāng)扭矩繼續(xù)增大，管道端面的變形也隨之增大，管道端面的極限承載力有了明顯的下降，已經(jīng)低于一致橢圓度段管道的壓潰壓力，因此管道端面也就成了管道壓潰的初始位置.

2.4.3 管道外徑對(duì)管道壓潰壓力的影響

考慮到管道外徑會(huì)對(duì)其壓潰壓力產(chǎn)生影響，故改變管道外徑.敏感性因素計(jì)算應(yīng)保持相同的徑厚比，計(jì)算中采用徑厚比均為26的325,mm×12.5,mm以及104,mm×4,mm 管道，對(duì)于每種管道，分別設(shè)置 0.5%、1.0%、1.5%、2.0% 4種一致橢圓度，計(jì)算并對(duì)比2種管道在 T-P路徑下的壓潰壓力，計(jì)算結(jié)果如圖8所示.

圖8 4種一致橢圓度、2種外徑的管道的壓潰壓力對(duì)比Fig.8 Comparison of collapse pressure of pipeline with two different diameters and four kinds of ovality

由圖8可知，在具有相同一致橢圓度時(shí)，104,mm×4,mm管道與325,mm×12.5,mm管道相比，在不受扭矩作用或受到較小扭矩作用時(shí)的壓潰壓力更大.同時(shí)，對(duì)于104,mm×4,mm管道，扭矩達(dá)到10,000,N·m時(shí)就可以使管道的壓潰壓力發(fā)生明顯下降.而對(duì)于325,mm×12.5,mm 管道，扭矩達(dá)到 100,000,N·m時(shí)，管道壓潰壓力仍沒有明顯下降，這說明外徑越大的管道對(duì)于扭矩就具有越強(qiáng)的承載力.對(duì)于實(shí)際工程而言，如果管道在服役前后受到扭矩作用且扭矩較小，那么外徑相對(duì)較小的管道極限承載力更強(qiáng)，能適應(yīng)更大的外部水壓.另一方面，如果管道在服役前后受到扭矩且扭矩較大，那么應(yīng)該選用外徑較大的管道，由于其抵抗扭矩的能力更強(qiáng)，能提高管道在受到扭矩作用后的極限承載力，具有更好的安全性.

3 改變扭矩-水壓加載路徑

管道在海底服役中還可能會(huì)受到海洋環(huán)境載荷造成的扭矩.對(duì)于這種路徑的加載過程，管道先受到恒定外部水壓，之后受到扭矩的作用.可在有限元軟件中采用 P - T加載路徑，即先對(duì)管道施加扭矩，之后施加均布外部水壓，以模擬管道在海底服役過程中受到扭矩作用導(dǎo)致壓潰的工程實(shí)際.

P- T路徑中恒定外部水壓的數(shù)值，根據(jù)橢圓管道在 T - P路徑中壓潰壓力的數(shù)值而進(jìn)行選取.即對(duì)于相同的管道壓潰壓力，通過比較2種加載路徑所需施加扭矩的大小，研究 P - T加載路徑下扭矩對(duì)管道極限承載力的影響.

以具有 2%一致橢圓度的 325,mm×10,mm管道為例，在相同壓潰壓力下2種加載路徑下橢圓管道的扭矩大小差異如表3所示.

表3 325,mm×10,mm管道不同加載路徑扭矩差異Tab.3 Differences of torsion in load paths of pipeline with the parameter of 325,mm×10,mm

對(duì)于外徑為 325,mm、壁厚為 10,mm 的管道，分別設(shè)置 0.5%、1.0%、1.5%、2.0%、5.0% 5種一致橢圓度.將具有不同一致橢圓度的 325,mm×10,mm管道在2種加載路徑下的扭矩差異曲線如圖9所示.

如圖 9所示，對(duì)于具有不同一致橢圓度的325,mm×10,mm管道而言，其壓潰壓力越小，2種加載路徑得到的扭矩差異普遍越大，且管道壓潰壓力相同時(shí)， T-P加載路徑中施加的扭矩普遍大于 P-T加載路徑中施加的扭矩.需要注意的是，當(dāng)扭矩很小時(shí)，2種加載路徑得到的扭矩值出現(xiàn)較大差異，甚至出現(xiàn)了接近-100%的誤差，這說明此時(shí) P-T路徑扭矩遠(yuǎn)大于 T-P路徑扭矩.而當(dāng)扭矩較大時(shí)，差異則在30%以內(nèi)且呈平穩(wěn)增長(zhǎng)趨勢(shì).由此可知， P-T加載路徑下的扭矩變化幅值較小，極小值相對(duì)于 T-P加載路徑極小值要大，而極大值相對(duì)于 T-P加載路徑極大值要小.這說明 P-T加載路徑下的橢圓管道的初始抗扭能力較強(qiáng)，但對(duì)于扭矩更加敏感，一旦扭矩超過了一定值，就可以使橢圓管道的壓潰壓力迅速下降.在實(shí)際工程中，管道在服役前和服役中受到的扭矩通常較小.通過上述分析可知，扭矩較小時(shí)，服役前的管道更容易發(fā)生屈曲壓潰，因此應(yīng)該避免管道在生產(chǎn)、運(yùn)輸和鋪設(shè)過程中受到扭矩作用，影響其極限承載力.而對(duì)于服役中的管道，在受到外部水壓作用后，其極限承載力對(duì)于扭矩更加敏感，當(dāng)扭矩較大時(shí)更容易發(fā)生屈曲壓潰.因此在鋪設(shè)管道時(shí)，應(yīng)注意避開海洋環(huán)境載荷多變且劇烈的地點(diǎn)，同時(shí)要注重海洋環(huán)境監(jiān)測(cè)，提前做好管道的加固等防范性工作.

圖9 10,mm壁厚管道扭矩路徑差異曲線Fig.9 Difference curve of torsion in two load paths of pipeline with the thickness of 10,mm

與此同時(shí)，對(duì)于具有不同一致橢圓度的管道而言，橢圓度越大，2種加載路徑得到的扭矩差異越小.這說明管道的一致橢圓度越大，其對(duì)扭矩的敏感性就越弱，反之管道的一致橢圓度越小，其對(duì)扭矩的敏感性就越強(qiáng).從實(shí)際工程角度考慮，管道的初始橢圓度通常較小，一般范圍為 0.1%～2.0%.而對(duì)于一致橢圓度小的管道，其極限承載力比一致橢圓度大的管道更強(qiáng)，但對(duì)扭矩更加敏感，因此同樣需要避免管道在服役前和服役中受到扭矩作用影響其極限承載力.

4 結(jié) 論

通過對(duì)海底橢圓管道的理論計(jì)算，得到了其在彈性階段的位移方程.利用有限元弧長(zhǎng)法計(jì)算得到了受扭矩和水壓作用的橢圓管道在壁厚、一致橢圓度、外徑等敏感性因素影響下的壓潰壓力.并改變扭矩-水壓加載路徑，研究橢圓管道在海底服役過程中受到扭矩作用的屈曲壓潰.

(1)對(duì)于等壁厚橢圓管道，在服役前受到扭矩作用的彈性位移可通過理論計(jì)算得到.管道在該過程中水平方向和豎直方向的變形位移較大，而沿管道軸向方向的翹曲位移較小.

(2)利用有限元中的弧長(zhǎng)法計(jì)算得到的橢圓管道極限承載力是準(zhǔn)確的.利用該方法對(duì)不同尺度、不同加載路徑下管道壓潰壓力的大小進(jìn)行計(jì)算，結(jié)果發(fā)現(xiàn)管道在不同尺度和不同加載情況下的極限承載力有較大差別，管道的屈曲破壞形式也具有不同的特點(diǎn).

(3)對(duì)于海底管道的設(shè)計(jì)，為了提高其安全性和穩(wěn)定性，應(yīng)考慮管道在加工過程中和在服役前后可能會(huì)形成的一致橢圓度以及可能會(huì)承受的扭矩，并根據(jù)實(shí)際生產(chǎn)需求選擇恰當(dāng)壁厚及外徑的海底管道.

［1］Yu Jianxing，Wu Mengning，Sun Zhenzhou，et al.Research progress of buckling propagation experiment of deep-water pipelines[J].Transactions of Tianjin University，2016，22(4)：285-293.

［2］余建星，馬維林，陳飛宇，等. 海底管道止屈器形式與設(shè)計(jì)方法研究概述[J]. 海洋工程，2013，31(4)：100-105.Yu Jianxing，Ma Weilin，Chen Feiyu，et al. Research on types and design method of buckling arrestors for subsea pipelines[J].The Ocean Engineering，2013，31(4)：100-105(in Chinese).

［3］Kyriakides S. Propagating instabilities in structures[J].Advances in Applied Mechanics，1993，30(114)：67-189.

［4］Longva V，Saevik S. On prediction of torque in flexible pipe reeling operations using a Lagrangian-Eulerian FE framework[J].Marine Structures，2016，46：229-254.

［5］Merodio J，Ogden R W. Extension，inflation and torsion of a residually stressed circular cylindrical tube[J].Continuum Mechanics and Thermodynamics，2016，28(1)：157-174.

［6］Bardi F C，Kyriakides S. Plastic buckling of circular tubes under axial compression(Part I)：Experiments[J].International Journal of Mechanical Sciences，2006，48(8)：830-841.

［7］祝恩淳，Calladine C R. 圓柱薄殼在局部軸向壓力作用下的屈曲[J]. 工程力學(xué)，2003，20(2)：168-170.Zhu Enchun，Calladine C R. Buckling of thin cylindrical shells under locally axial compression[J].Engineering Mechaincs，2003，20(2)：168-170(in Chinese).

［8］Wang C Y. Torsion of tubes of arbitrary shape[J].International Journal of Solids & Structures，1998，35(s7/8)：719-731.

［9］Arghavan S，Hematiyan M R. Torsion of functionally graded hollow tubes[J].European Journal of Mechanics，2009，28(3)：551-559.

［10］Sharma R，Aggarwal A K，Sharma S. Collapse pressure analysis in torsion of a functionally graded thick-walled circular cylinder under external pressure[J].Procedia Engineering，2014，86(4)：738-747.

［11］Zheleznov L P，Kabanov V V，Boiko D V. Nonlinear deformation and stability of supported oval cylindrical shells under torsion and bending with internal pressure[J].Russian Aeronautics(Iz VUZ)，2010，53(4)：466-469.

［12］Timoshenco S P，Goodier J N.Theory of Elasticity[M].2nd ed.New York：McGraw-Hill，1951.

［13］Muskhelishvili N I.Some Basic Problems of the Mathematical Theory of Elasticity[M]. Netherlands：Noordhoff，Groningen，2009.

［14］劉玉卿，張振永，楊 帆. 管道過量塑性變形值Ramberg-Osgood方程研究[J]. 中國(guó)石油和化工標(biāo)準(zhǔn)與質(zhì)量，2014(11)：239-240.Liu Yuqing，Zhang Zhenyong，Yang Fan. Study on Ramberg-Osgood equation of over plastic deformation[J].China Petroleun and Chemical Standard andQuality，2014(11)：239-240(in Chinese).

［15］余建星，李 驍，徐立新，等. 深水海底管道凹陷機(jī)理[J]. 天津大學(xué)學(xué)報(bào)：自然科學(xué)與工程技術(shù)版，2015，48(11)：1009-1013.Yu Jianxing，Li Xiao，Xu Lixin，et al. Denting mechanism of deepwater pipeline[J].Journal of Tianjin University：Science and Technology，2015，48(11)：1009-1013(in Chinese).

［16］陳飛宇，余建星，趙羿羽，等. 復(fù)雜載荷條件下有缺陷海底管道非線性屈曲分析[J]. 中南大學(xué)學(xué)報(bào)：自然科學(xué)版，2015，46(7)：2701-2706.Chen Feiyu，Yu Jianxing，Zhao Yiyu，et al. Nonlinear buckling of subsea pipes with imperfection under complex loads[J].Journal of Central South University：Science and Technology，2015，46(7)：2701-2706(in Chinese).

［17］Crisfield M A. An arc-length method including line searches and accelerations[J].International Journal for Numerical Methods in Engineering，1983，19(9)：1269-1289.

［18］楊 源. 海底管道屈曲失效機(jī)理及其防護(hù)研究[D]. 天津：天津大學(xué)建筑工程學(xué)院，2013.Yang Yuan. The Research of Submarine Pipeline Buckling Failure and Buckle Arrestor[D]. Tianjin：School of Civil Engineering，Tianjin University，2013(in Chinese).

［19］Det Norske Veritas(DNV). Offshore Standard DNV-OSF101 Submarine Pipeline Systems[S]. Hφvik：DNV，2007.