考慮裂紋影響的海洋管道極限彎曲載荷分析

，，

(1.北京石油化工學(xué)院 機(jī)械工程學(xué)院，北京 102617；2.中國(guó)石油大學(xué)(北京) 海洋油氣研究中心，北京 102249)①

管道在卷管鋪設(shè)過(guò)程中受到隨機(jī)載荷、冷彎、地面運(yùn)動(dòng)等作用，不可避免地發(fā)生了塑性變形。管道上卷時(shí)受到的彎曲載荷較大，塑性變形也較大。如果沒(méi)有準(zhǔn)確預(yù)估彎曲載荷的大小，海管在上卷過(guò)程中很可能會(huì)因局部的壓潰，甚至斷裂而造成整根海管失效。尤其是對(duì)于焊接中產(chǎn)生的裂紋，如圖1所示，在彎曲載荷作用下裂紋逐漸變大，且截面發(fā)生了橢圓化變形，使裂紋迅速擴(kuò)展，從而使得海管失效速度加快。

裂紋會(huì)對(duì)海管的極限彎曲載荷產(chǎn)生一定的影響，準(zhǔn)確預(yù)估海管的極限彎曲載荷對(duì)保證海管安全運(yùn)行有重要的意義。

圖1 含有裂紋的海管受到彎曲載荷作用發(fā)生變形

對(duì)含有裂紋的管道進(jìn)行極限彎曲載荷研究，最早的是Miller[1]和Zahoor[2]，他們分別提出了含有裂紋管道的極限彎曲載荷表達(dá)式，被后來(lái)的很多學(xué)者引用，并以此為基礎(chǔ)對(duì)裂紋問(wèn)題做進(jìn)一步的探討和研究。其中，Chattopadhyay對(duì)含裂紋的管道研究的較多，包括用試驗(yàn)和解析方法研究含裂紋管道的三點(diǎn)彎曲問(wèn)題、含縱向裂紋的管道受到面內(nèi)彎矩作用的塑性壓潰問(wèn)題、含周向裂紋的管道受到面內(nèi)彎矩和內(nèi)壓聯(lián)合作用的壓潰問(wèn)題等[3-8]。Yahiaooui等總結(jié)了最近幾年關(guān)于含有縱向和周向裂紋的碳鋼管道受到彎曲載荷作用的試驗(yàn)和有限元模擬的研究，并對(duì)這些研究成果做了比較，指出以前研究人員計(jì)算得到的數(shù)據(jù)較為保守，可為實(shí)際工程提供參考[9]。Lei在文獻(xiàn)[10]中提出了關(guān)于含有軸向裂紋的厚壁筒受到內(nèi)壓作用時(shí)的極限載荷求解方法，用這種方法求解的結(jié)果與有限元計(jì)算結(jié)果吻合良好。在國(guó)內(nèi)，對(duì)含有裂紋的管道的研究也較多。帥健等給出了基于J積分理論的含裂紋管道的失效評(píng)定曲線的基本方程，并用有限元方法計(jì)算裂紋管道的J積分[11]。韓克江等針對(duì)現(xiàn)有的斷裂力學(xué)方法應(yīng)用于管道大范圍屈服斷裂評(píng)估的局限性，分別從雙參數(shù)斷裂力學(xué)、基于約束校正的斷裂韌性測(cè)試、基于應(yīng)變的斷裂評(píng)估和基于應(yīng)變的失效評(píng)估4個(gè)方面詳細(xì)地介紹了管道大范圍屈服斷裂評(píng)估的研究現(xiàn)狀，指出目前基于約束校正的管道斷裂韌性測(cè)試的主要方法是用SENT試件方法和表觀斷裂韌性方法[12]。周晶等分析了海底管線全壽命安全運(yùn)行的關(guān)鍵問(wèn)題[13]。金偉良等介紹了結(jié)構(gòu)失效評(píng)估圖(FAD)技術(shù)，詳細(xì)闡述了BS7910：1999規(guī)范中3個(gè)不同等級(jí)的斷裂評(píng)估方法和疲勞評(píng)估方法[14]。對(duì)于卷管鋪設(shè)方法中所使用的海管，Netto等研究了其疲勞壽命評(píng)估方法[15]，并在文獻(xiàn)[16]中研究了含有環(huán)向焊縫且有預(yù)應(yīng)力的管道的疲勞問(wèn)題。Tkaczyk等研究了受循環(huán)彎曲載荷作用的管道的疲勞壽命評(píng)估[17]，并在文獻(xiàn)[18]中研究了帶有屈服平臺(tái)的管道材料的非線性斷裂評(píng)估方法。

綜合以上研究?jī)?nèi)容，關(guān)于含有裂紋管道的研究較多，但大多集中在關(guān)于裂紋擴(kuò)展、疲勞裂紋研究等問(wèn)題，而對(duì)于含裂紋且有初始彎曲曲率的管道受彎曲載荷作用的研究較少。本文針對(duì)該問(wèn)題進(jìn)行了研究，首先總結(jié)已有的無(wú)裂紋管道的極限彎曲載荷表達(dá)式，并簡(jiǎn)述其適用范圍，包括本文所推導(dǎo)的無(wú)裂紋管道的極限彎曲載荷表達(dá)式。然后，列出已有的比較有代表性的含裂紋的管道的極限彎曲載荷表達(dá)式。本文用有限元軟件ABAQUS模擬有初始裂紋且有初始彎曲曲率的海管受彎曲載荷作用的多種工況，并將計(jì)算結(jié)果與已有的表達(dá)式結(jié)果進(jìn)行比較，得出本文的有關(guān)含有初始裂紋的彎曲海管的極限彎曲載荷表達(dá)式，并對(duì)其影響因素做了進(jìn)一步探討。

1 極限彎曲載荷理論

有裂紋的海管受到彎曲載荷作用后，將會(huì)有拉伸應(yīng)力(e-f弧線)和壓縮應(yīng)力(c-d弧線)產(chǎn)生，如圖1所示。這種拉伸和壓縮應(yīng)力使海管的截面發(fā)生橢圓化變形，而裂紋在cd弧線段上面，若裂紋尺寸小于cd弧線尺寸，則裂紋也是在受壓的一側(cè)，且不會(huì)張開(kāi)，在海管彎曲成一定曲率之前不會(huì)影響海管的穩(wěn)定性。反之，裂紋將會(huì)使極限彎曲載荷發(fā)生變化，會(huì)對(duì)海管的穩(wěn)定性產(chǎn)生一定的影響，這與無(wú)裂紋海管的極限彎曲載荷有一定關(guān)系。

1.1 無(wú)裂紋海管的極限彎曲載荷理論

海管在受到彎曲載荷作用時(shí)，其極限彎曲載荷即為海管塑性失穩(wěn)或者局部壓潰的載荷，塑性失穩(wěn)載荷的特點(diǎn)為載荷單調(diào)增加時(shí)的最大載荷，在該處海管發(fā)生了塑性變形。對(duì)于極限彎曲載荷的計(jì)算方法，在很多文獻(xiàn)中都有涉及，文獻(xiàn)[19]中提出含有初始彎曲曲率的管道的極限彎曲載荷為：

(1)

h=tRb/R2

(2)

式中：σy為管道的屈服強(qiáng)度；h為管道彎曲性質(zhì)參數(shù)；Rb為管道的彎曲半徑；R為管道平均半徑；D為管道外徑；t為管道壁厚。

Calladine在文獻(xiàn)[20]中得出含有初始彎曲曲率的管道的極限彎曲載荷為：

(3)

以上2個(gè)表達(dá)式均基于小變形分析，且管道材料均假設(shè)為理想彈性模型。對(duì)于大變形分析，Goodall在文獻(xiàn)[21]中提出了另外一種表達(dá)形式，即：

(4)

其中：

(5)

式中：ν為泊松比；E為彈性模量。

Touboul基于實(shí)驗(yàn)提出了塑性壓潰的2種彎矩表達(dá)式[22]。

對(duì)于圖2，極限彎曲載荷表達(dá)式為：

M0=0.715h2/3(4R2tσy)

(6)

對(duì)于圖3，極限彎曲載荷表達(dá)式為：

M0=0.722h1/3(4R2tσy)

(7)

圖2 受正向彎曲載荷作用的海管

圖3 受反向彎曲載荷作用的海管

本文在文獻(xiàn)[23]基礎(chǔ)上提出了含有初始彎曲曲率海管的極限彎曲載荷表達(dá)式為：

(8)

式中：Δ為橢圓度；k為管道彎曲曲率；t為壁厚；r為管道彎曲半徑；k0為初始彎曲曲率。

(9)

然而，對(duì)于含有初始裂紋的海管，無(wú)法用上面的式子進(jìn)行計(jì)算。故需要對(duì)含有裂紋的海管進(jìn)行分析和研究。

1.2 有裂紋海管極限彎曲的載荷理論

含裂紋的截面的彎矩可由式(10)求得。

(10)

式中：M為彎矩；σ為局部應(yīng)力；Z為含裂紋的截面模量。

隨著加載彎矩的逐漸變大，局部應(yīng)力σ隨著管道材料的硬化也在逐漸發(fā)生變化，其值與加載方式有關(guān)。

對(duì)于含有裂紋和初始彎曲曲率的管道的極限彎曲載荷，Miller和Zahoor最先提出了各自的彎曲載荷表達(dá)式，并被后來(lái)的學(xué)者們作為研究的基礎(chǔ)和參考。

針對(duì)管道部分周向貫穿裂紋，Miller[1]提出的極限彎曲載荷表達(dá)式為：

(11)

M0=0.935(4R2tσf)h2/3

(12)

h=tRb/R2

(13)

式中：3θ為總的周向裂紋角。

針對(duì)部分周向貫穿裂紋，Zahoor[2]提出的極限彎曲載荷表達(dá)式為：

(14)

式(14)適用于：

(15)

式中：a為裂紋長(zhǎng)度的一半；Dm為管道的平均彎曲半徑。

Zahoor還提出了含有縱向貫穿裂紋的極限彎曲載荷表達(dá)式。

(16)

然而，這些公式并沒(méi)有區(qū)分彎曲載荷的施加方式(如圖2和圖3)。式(8)是基于其施加方式的不同，并考慮海管在變形中截面發(fā)生了橢圓化變形而進(jìn)行推導(dǎo)的，本文將基于此，對(duì)含有裂紋的海管的極限彎曲載荷做進(jìn)一步的研究。假設(shè)：

ML=M0X

(17)

其中，M0由式(8)得。X可通過(guò)后面的有限元模擬計(jì)算，并對(duì)計(jì)算結(jié)果進(jìn)行擬合得出。

2 有限元計(jì)算分析

本文對(duì)卷管鋪設(shè)過(guò)程中的海管采用大型通用非線性有限元軟件ABAQUS進(jìn)行數(shù)值模擬。海管外徑219 mm，壁厚12 mm，長(zhǎng)度為40 mm，初始彎曲曲率為-0.01 m-1，即，施加彎矩并使海管的彎曲曲率與初始的彎曲曲率方向相反，海管的部分周向貫穿裂紋角的一半為π/16，材料參數(shù)選取X65鋼。在Assembly模塊中通過(guò)劃分實(shí)體，得到所要的初始裂紋的形狀，在Interaction模塊里Special-Crack-Assign Seam定義初始裂紋。對(duì)裂紋附近區(qū)域劃分了較為密集的網(wǎng)格，如圖4所示。在模擬中，施加位移載荷來(lái)控制海管的變形，不考慮裂紋的擴(kuò)展。當(dāng)局部塑性應(yīng)變達(dá)到一定程度，即可判定海管失穩(wěn)，計(jì)算結(jié)果如圖5所示。

從圖5的海管變形圖中可以看出，裂紋處海管發(fā)生的變形最大，且有應(yīng)力集中現(xiàn)象。隨著彎矩載荷的變大，裂紋處最先發(fā)生局部失穩(wěn)，故裂紋對(duì)海管的極限彎矩載荷有較大影響。

圖4 海管有限元模型

圖5 管道塑性應(yīng)變累計(jì)分布

本文建立了裂紋半角為π/16、π/8、3π/16、π/4、5π/16、3π/8、7π/16、π/2的有限元模型，并分別進(jìn)行計(jì)算。表1列出了式(11)、式(14)的計(jì)算結(jié)果，以及本文所做的有限元計(jì)算結(jié)果，并將計(jì)算結(jié)果繪制在圖6中。

表1 含裂紋海管的極限彎矩載荷計(jì)算結(jié)果

注：1) 裂紋半角為π/16、π/8、3π/16等。

2)ML為含裂紋管道極限彎曲載荷。

3)M0為無(wú)裂紋管道極限彎曲載荷。

圖6 極限彎曲載荷與裂紋角度對(duì)應(yīng)關(guān)系

從圖6中可看出，3種方法得到的計(jì)算結(jié)果均表明，隨著裂紋角的增大，極限彎曲載荷逐漸減小。Zahoor結(jié)果與有限元計(jì)算結(jié)果略高于Miller結(jié)果。當(dāng)裂紋角超過(guò)π/2后，有限元計(jì)算結(jié)果將迅速減小，這主要是由于海管局部發(fā)生失穩(wěn)，極限彎曲載荷會(huì)迅速減小。有限元計(jì)算結(jié)果與實(shí)際的海管變形情況更加接近，其可靠度更高一些。對(duì)有限元計(jì)算結(jié)果進(jìn)行擬合，得出：

(18)

從式(18)可看出，有裂紋海管的極限彎曲載荷與裂紋角有關(guān)，且和無(wú)裂紋海管的極限彎曲載荷有關(guān)。無(wú)裂紋海管的極限彎曲載荷又受海管管徑、壁厚、初始彎曲曲率等參數(shù)的影響。故有裂紋海管的極限彎曲載荷也與管徑、壁厚、初始彎曲曲率有關(guān)。

3 影響彎曲載荷變化的主要因素

針對(duì)式(18)，分析有裂紋海管的極限彎曲載荷隨影響因素變化的規(guī)律時(shí)，首先求出無(wú)裂紋海管的極限彎曲載荷，這主要通過(guò)式(8)和式(9)求得。海管幾何參數(shù)的變化主要是管外徑和壁厚的變化，本計(jì)算中固定海管管徑為219 mm，只對(duì)壁厚的變化進(jìn)行研究。除了海管的幾何參數(shù)外，海管的初始彎曲曲率也對(duì)海管的極限彎曲載荷產(chǎn)生影響，故對(duì)海管的初始彎曲曲率進(jìn)行研究，負(fù)值表示海管加載時(shí)海管產(chǎn)生彎曲曲率的方向與海管最初的彎曲曲率方向相反，正值表示相同(例如圖2與圖3)。對(duì)t=12、14、16、18 mm，和k0=-0.01、-0.02、-0.03、-0.04、0.01、0.02、0.03、0.04 m-1分別進(jìn)行了計(jì)算，計(jì)算工況如表2所示。表3列出了無(wú)裂紋海管的極限彎曲載荷值，并將這些數(shù)據(jù)用Origin軟件進(jìn)行繪圖，如圖7所示。

表2 計(jì)算工況

注：k0為初始彎曲曲率；t為壁厚。

表3 無(wú)裂紋海管的極限彎曲載荷M0計(jì)算結(jié)果

圖7 無(wú)裂紋海管的極限彎曲載荷

從圖7中可以看出，無(wú)裂紋海管的極限彎曲載荷隨著壁厚的增大而增大，且增長(zhǎng)趨勢(shì)較快。隨初始彎曲曲率的增大而減小，但變化趨勢(shì)不是很明顯。由此可見(jiàn)，增加管壁厚有利于提高無(wú)裂紋海管的抗彎能力。

3.1 壁厚

對(duì)于具有相同壁厚的海管，其極限彎曲載荷由于受到裂紋的影響將會(huì)大幅減小。本文對(duì)初始彎曲曲率為-0.01 m-1的海管進(jìn)行了研究。壁厚t=12、14、16、18 mm，裂紋半角分別為θ=π/16、π/8、3π/16、π/4、5π/16、3π/8、7π/16、π/2，根據(jù)式(18)進(jìn)行計(jì)算，計(jì)算結(jié)果如表4所示。將表4中數(shù)據(jù)繪圖，如圖8所示。

表4 k0=-0.01 m-1時(shí)含裂紋海管的ML隨壁厚變化計(jì)算結(jié)果

圖8 k0=-0.01 m-1時(shí)的有裂紋海管的極限彎曲載荷與壁厚、裂紋半角的關(guān)系

從圖8中可以看出，極限彎曲載荷隨裂紋半角的增大而減小，隨壁厚的增大而增大。在裂紋為零處，極限彎曲載荷變化趨勢(shì)與無(wú)裂紋海管的極限彎曲載荷變化趨勢(shì)一致。壁厚較大的海管，極限彎曲載荷也較大，但其隨裂紋半角的增大，極限彎曲載荷變小的趨勢(shì)更快一些。對(duì)于裂紋半角較大的海管，其極限彎曲載荷隨壁厚的增大而增大，但增長(zhǎng)趨勢(shì)比較緩慢。

3.2 初始彎曲曲率

初始彎曲曲率不僅影響無(wú)裂紋海管的極限彎曲載荷，也對(duì)有裂紋海管的極限彎曲載荷產(chǎn)生一定的影響，在壁厚t=12 mm時(shí)，研究初始彎曲曲率和裂紋半角的大小對(duì)極限彎曲載荷的影響，初始彎曲曲率的取值為-0.01，-0.02，-0.03，-0.04，0.01，0.02，0.03，0.04，裂紋半角分別為θ=π/16、π/8、3π/16、π/4、5π/16、3π/8、7π/16、π/2，根據(jù)(18)式進(jìn)行計(jì)算，計(jì)算結(jié)果如表5所示，圖9為其相應(yīng)的三維圖。

表5 t=12 mm時(shí)含裂紋海管的ML隨初始彎曲曲率變化的計(jì)算結(jié)果

圖9 t=12mm時(shí)，海管極限彎曲載荷與裂紋半角、初始彎曲曲率的對(duì)應(yīng)關(guān)系

從圖9中可看出，對(duì)于有同一裂紋半角的海管，隨著初始彎曲曲率的增大，其極限彎曲載荷減小。無(wú)論是正向的初始彎曲曲率，還是負(fù)向的初始彎曲曲率，都沒(méi)有改變極限彎曲載荷隨裂紋半角的增大而減小的曲線變化規(guī)律。在裂紋半角接近π/2時(shí)，極限彎曲載荷隨初始彎曲曲率變化的范圍較小，這主要是由于裂紋角越大，海管的局部越容易產(chǎn)生應(yīng)力集中，其局部越容易失穩(wěn)，使得海管的極限彎曲載荷越小。

4 結(jié)論

本文總結(jié)了有關(guān)無(wú)裂紋海管的極限彎曲載荷的理論計(jì)算方法，其中包括本文所提出的理論計(jì)算公式。在此基礎(chǔ)上，分析了Miller和Zahoor有關(guān)含裂紋海管的極限彎曲載荷理論計(jì)算方法。用有限元方法進(jìn)行了計(jì)算，得到了多組裂紋半角對(duì)應(yīng)的極限彎曲載荷值，并用Origin軟件進(jìn)行了函數(shù)擬合，得到了有裂紋海管的極限彎曲載荷表達(dá)式。該表達(dá)式考慮了正向的初始彎曲曲率和負(fù)向的初始彎曲曲率以及在彎曲變形中海管截面的橢圓化變形。通過(guò)對(duì)該表達(dá)式中的壁厚、初始彎曲曲率、裂紋半角的研究得到的結(jié)論如下：

1) 極限彎曲載荷隨裂紋半角的增大而減小，隨壁厚的增大而增大，在裂紋為零處，極限彎曲載荷變化趨勢(shì)與無(wú)裂紋海管的極限彎曲載荷變化趨勢(shì)一致。壁厚較大的海管，極限彎曲載荷也較大，但其隨裂紋半角的增大，極限彎曲載荷變小的趨勢(shì)更快一些。對(duì)于裂紋半角較大的海管，其極限彎曲載荷隨壁厚的增大而增大，但增大趨勢(shì)趨于緩慢，故較大的壁厚有利于提高海管的極限彎曲載荷。

2) 對(duì)于有同一裂紋半角的海管，隨著初始彎曲曲率的增大，其極限彎曲載荷減小。無(wú)論是正向的初始彎曲曲率還是負(fù)向的初始彎曲曲率，都沒(méi)有改變極限彎曲載荷隨裂紋半角的增大而減小的曲線變化規(guī)律。在裂紋半角接近π/2時(shí)，極限彎曲載荷隨初始彎曲曲率變化的范圍較小，這主要是由于裂紋角越大，海管的局部越容易產(chǎn)生應(yīng)力集中，其局部越容易失穩(wěn)，使得海管的極限彎曲載荷越小。

3) 海管在進(jìn)行焊接時(shí)，將會(huì)有殘余應(yīng)力、焊縫或者各種結(jié)構(gòu)形式的裂紋等缺陷的存在。在對(duì)海管進(jìn)行多次加載-卸載的過(guò)程中，其缺陷使得海管材料力學(xué)性能發(fā)生了大的改變，遠(yuǎn)比完整海管的力學(xué)性能復(fù)雜的多，且這些缺陷在循環(huán)彎曲載荷作用下會(huì)對(duì)海管的疲勞壽命產(chǎn)生很大的影響。建議用理論、試驗(yàn)、數(shù)值模擬等多種方法結(jié)合來(lái)進(jìn)行研究。

4) 在進(jìn)行卷管鋪設(shè)時(shí)，海管的材料有多種，其結(jié)構(gòu)也較為復(fù)雜，例如，雙層管，夾心管，涂有外涂層的海管等，像這種結(jié)構(gòu)和材料都較為復(fù)雜的海管，在受到多種復(fù)雜的外載荷作用時(shí)，其變形也更為復(fù)雜，未來(lái)的研究方向可以考慮復(fù)雜結(jié)構(gòu)的海管在卷管鋪設(shè)中的彈塑性變形。

：

[1] Miller A G. Review of limit loads of structures containing defects[J]. International Jouranla of Pressure Vessels Piping, 1988,32: 197-327.

[2] Zahoor A. Ductile fracture handbook[M]. California USA: Electric Power Research Institute, 1991: 11-34.

[3] Chattopadhyay J, Dutta B K, Kushwaha H S. Experimental and analytical study of three point bend specimen and throughwall circumferentially cracked straight pipe[J]. International Journal of Pressure Vessels and Piping, 2000, 77:455-471.

[4] Chattopadhyay J ，Pavankumar T V, Dutta B K, et al. Fracture experiments on through wall cracked elbows under in-plane bending moment: Test results and theoretical/numerical analyses[J]. Engineering Fracture Mechanics, 2005, 72:1461-1497.

[5] Chattopadhyay J, Kushwaha H S, Roos E. Some recent developments on integrity assessment of pipes and elbows. Part I: Theoretical investigations[J]. International Journal of Solids and Structures, 2006, 43:2904-2931.

[6] Chattopadhyay J, Tomar A K S, Dutta B K, et al. Elastic-plastic J and COD estimation schemes for throughwall circumferentially cracked elbow under in-plane closing moment[J]. Engineering Fracture Mechanics, 2005, 72:2186-2217.

[7] Chattopadhyay J, Venkatramana W, Dutta B K, et al. Plastic collapse moment equations of throughwall axially cracked elbows subjected to in-plane bending moment[J]. Engineering Fracture Mechanics, 2008, 75:2260-2275.

[8] Chattopadhyay J, Tomar A K S. New plastic collapse moment equations of defect-free and throughwall circumferentially cracked elbows subjected to combined internal pressure and in-plane bending moment[J]. Engineering Fracture Mechanics, 2006, 73:829-854.

[9] Yahiaoui K, Moreton D N, Moffat D G. Evaluation of limit load data for cracked pipe bends under opening bending and comparisons with existing solutions[J]. International Journal of Pressure Vessels and Piping, 2002, 79:27-36.

[10] Lei Y. A review of limit load solutions for cylinders with axial cracks and development of new solutions[J]. International Journal of Pressure Vessels and Piping, 2008, 85:825-850.

[11] 帥?。鸭y管道的失效評(píng)定曲線的實(shí)例驗(yàn)證[J]．機(jī)械強(qiáng)度，2003，25(3)：251-253．

[12] 韓克江，帥健，王俊強(qiáng)，等．管道大范圍屈服斷裂評(píng)估的研究現(xiàn)狀與進(jìn)展[J]．力學(xué)進(jìn)展，2012，42(1)：29-40．

[13] 周晶，馮新，李昕．海底管線全壽命安全運(yùn)行的關(guān)鍵問(wèn)題研究[J]．工程力學(xué)，2011(增刊Ⅱ)：97-108．

[14] 金偉良，付 勇，趙冬巖，等．具有裂紋損傷的海底管道斷裂及疲勞評(píng)估[J]．海洋工程，2005，23(13)：7-16．

[15] Netto T A,Botto A. Fatigue Life Assessment of Reeled Risers[C]// Offshore Conference Technology, 2007, 1-7.

[16] Netto T A, Botto A, Lourenco M I. Fatigue performance of pre-strained pipes with girth weld defects: local deformation mechanisms under bending[J]. International journal of fatigue, 2008, 30: 1080-1091.

[17] Tkaczyk T, Dowd N P O, Nikbin K, et al. Fracture Assessment of Elastic-Plastic Steel Pipelines Subject to Multi-cycle Bending[C]//Proceedings of the Nineteenth (2009) International Offshore and Polar Engineering Conference, 2009,110-117.

[18] Tkaczyk T, Dowd N P O,Nikbin K, et al. A Non-Linear Fracture Assessment Procedure for Pipeline Materials with a Yield Plateau[C]//Proceedings of the Nineteenth (2009) International Offshore and Polar Engineering Conference, 2009, 100-109.

[19] Spence J, Findlay G E. Limit load for pipe bends under in-plane bending[C]// Proceeding of 2nd international conference on pressure vessel technology, 1973, 393-403.

[20] Calladine C R. Limit analysis of curved tubes[J].Journal of Mechanical Engineering Science, 1974, 16(2): 85-7.

[21] Goodall I W. Large deformations in plastically deforming curved tubes subjected to in-plane bending[R]. UK: CEGB RD/B/N4312, 1978.

[22] Touboul F, Ben D M, Acker D. Design criteria for piping components against plastic collapse: application to pipe bend experiments[C]//Proceeding of 6th international conference on pressure vessel technology, 1988, 73-84.

[23] 張九菊，段夢(mèng)蘭，馬建敏，等．基于深海卷管鋪設(shè)的海管橢圓度分析[J]．應(yīng)用數(shù)學(xué)和力學(xué)，2013，34(6)：552-563．