含腐蝕缺陷的非粘結(jié)柔性管骨架層承載分析

魯中歧 崔俊國 肖文生 王魁濤 張俊

摘要：柔性管最內(nèi)層骨架層與運輸介質(zhì)直接接觸，不可避免地會有腐蝕發(fā)生，進而影響管道的承壓能力。為研究腐蝕缺陷對非粘結(jié)柔性管骨架層承壓能力的影響，結(jié)合已有研究成果，借用ANSYS Workbench軟件建立了柔性管骨架層的有限元模型，研究了腐蝕形貌尺寸及腐蝕形式對骨架層失效壓力的影響規(guī)律，并借助分析結(jié)果構(gòu)建了馬爾可夫預(yù)測模型，對管道腐蝕狀態(tài)進行了預(yù)測。研究結(jié)果表明：小范圍的局部腐蝕缺陷對骨架層的承載能力影響并不大，而均勻腐蝕深度對其影響較為劇烈，且隨腐蝕深度的增加，塑性極限載荷下降趨勢明顯加快，根據(jù)該影響規(guī)律，擬合了以腐蝕深度為自變量、骨架層塑性極限載荷為因變量的關(guān)系式，經(jīng)檢驗，該模型誤差較小，擬合數(shù)據(jù)吻合；現(xiàn)有的腐蝕狀態(tài)劃分方式對非粘結(jié)柔性管骨架層不適用，可借助馬爾可夫預(yù)測模型進行預(yù)測。所得結(jié)論可為柔性管的腐蝕程度評價及安全使用提供借鑒。

關(guān)鍵詞：非粘結(jié)柔性管；腐蝕缺陷；均勻腐蝕；骨架層；承載能力；失效壓力；馬爾可夫

中圖分類號：TE832 文獻標識碼：A DOI：10.16082／j．cnki．issn．1001-4578.2022.12.019

Bearing Capacity of Unbonded Flexible Pipe Carcass with Corrosion Defects

Lu Zhongqi' Cui Junguo' Xiao Wensheng' Wang Kuitao Zhang Jun'

（1.College of Mechanical and Electrical Engineering，China University of Petroleum（East China）;2.CNOOC Research Institu- te Co.，Lid.）

Abstract：The innermost carcass of flexible pipe directly contacts with the transported medium， and corrosion occurs inevitably，which will affect the pressure bearing capacity of the pipeline. Inorder to investigate the influ- ence of corrosion defects on the pressure bearing capacity of unbonded flexible pipe carcass， based on the existing research results，a finite element model of flexible pipe carcass was established by using the ANSYS Workbench software.The influences of corrosion morphology and size and corrosion form on the carcass failure pressure were studied.Based on the analysis resulIs， a Markov prediction model was established to predict the corrosiomn state of the pipeline.The results show that small local corrosion defects have little influence on the bearing capacity of car- cass，but uniform corrosion depths have a significant influence on the bearing capacity.With the increase of corro- sion depth， the plastic limit load decreases obviously. Taking the corrosion depth as the independent variable and the carcass plastic limit load as the dependent variable，a fitting relationship was established. The verification indi- cates that this model has small error and well agrees with the fitting data. The available classification method of cor- rosion state is not applicable to the unbonded flexible pipe carcass， and the corrosion state can be predicted by means of Markov prediction model.The conclusions provide references for the corrosion assessment of flexible pipe and its safe application.

Keywords：unbonded flexible pipe; corrosion defect; uniform corrosion; carcass; bearing capacity;fail- ure pressure;Markov

0引言

腐蝕缺陷嚴重影響了油氣管道的承載能力，進而影響其運載壽命。目前，國內(nèi)外學者針對腐蝕缺陷形貌尺寸對傳統(tǒng)剛性管道承載能力影響關(guān)系進行了深入研究，其中部分專家借助數(shù)值分析軟件建立了管道有限元模型，對矩形和球形等不同形狀的單個腐蝕缺陷進行了模擬，分析了腐蝕形貌對管道失效壓力的影響；不少學者也對腐蝕缺陷的長度、寬度及徑向深度對管道承載能力的影響規(guī)律進行了研究，得出了腐蝕缺陷深度對管道失效壓力影響最大、缺陷長度次之、缺陷寬度影響最小的結(jié)論［2—4］。還有一些研究人員對存在軸向、環(huán)向及交疊等情況的組合腐蝕缺陷對管道等效應(yīng)力的影響規(guī)律進行了探討，得出結(jié)論如下：軸向及環(huán)向存在多個腐蝕坑時，存在一個臨界距離，大于該臨界距離時，兩腐蝕缺陷不會相互影響；對于交疊分布缺陷，小缺陷對等效應(yīng)力的影響不可忽略［5—7］。

以上研究都是針對于傳統(tǒng)剛性管，而柔性管道相比剛性管道具有更大的優(yōu)勢，逐步得到廣泛應(yīng)用。柔性管最內(nèi)層骨架層與運輸介質(zhì)直接接觸，不可避免地會有腐蝕發(fā)生，進而影響管道的承壓能力，導致管道失效。然而針對該腐蝕缺陷對于管道承載能力影響規(guī)律的研究還比較少。國內(nèi)外相關(guān)專家通過數(shù)值分析的方法，探討了材料非線性、初始橢圓度、相鄰節(jié)之間的接觸方式及螺旋纏繞角度等因素對均勻荷載或?qū)较驂毫ο碌墓羌軐邮毫Φ挠绊懀贸隽瞬牧戏蔷€性對骨架層的抗壓潰能力影響最大，初始橢圓度及接觸方式影響次之，螺旋纏繞角度影響最小的結(jié)論。

上述針對柔性管骨架層的研究都在無腐蝕缺陷情況下進行，有關(guān)腐蝕缺陷對非粘結(jié)柔性管骨架層的承壓能力研究較少。為此，筆者結(jié)合已有研究成果，借用 ANSYS Workbench 軟件建立了柔性管骨 架層的有限元模型，研究了腐蝕形貌尺寸及腐蝕形式對骨架層失效壓力的影響規(guī)律，并借助分析結(jié)果對骨架層的腐蝕狀態(tài)進行了劃分。研究結(jié)果可為柔性管腐蝕程度評價提供借鑒。

1理論模型研究

圓環(huán)在均勻外壓載荷下的屈曲載荷值［8］可結(jié) 合鐵木辛柯彈性穩(wěn)定理論中圓環(huán)撓曲線微分方程推導，具體如下：

式中：q為圓環(huán)均勻外壓下的屈曲載荷，MPa；I為圓環(huán)矩形截面慣性矩，m；b為矩形截面長，m；t為矩形截面寬，m；R為圓環(huán)半徑，m；E為材料彈性模量，MPa。

對于傳統(tǒng)均質(zhì)圓環(huán)，截面為長方形，于是單位長度下的平面應(yīng)力臨界壓潰值qa為：

由于非粘結(jié)柔性管骨架層截面形狀十分復(fù)雜，可依據(jù)下式計算等效厚度t：

式中：1c為骨架層截面的最小主慣性矩，㎡；A為截面面積，㎡。

本文選取內(nèi)徑101.6mm（4in）非粘結(jié)柔性管道為研究對象，骨架層截面尺寸如圖1所示［4］

該骨架層截面面積A為56.2107m㎡，最小主慣性矩Ic為181.1725mm4。

將以上參數(shù)帶入式（4）可求得骨架層的等效厚度t為6.2191mm。圓環(huán)半徑取中徑［”］，即有：

式中：R1為非粘結(jié)柔性管最內(nèi)層半徑，mm。

帶入式（5）求得R為52.409 5mm，骨架層材料為316L不銹鋼，彈性模量E為206 GPa。將彈性模量E、等效厚度t及半徑R帶入式（3）求得單位長度下平面應(yīng)力臨界壓潰值為79.06 MPa。

2 有限元模型分析

2.1 模型驗證

根據(jù)現(xiàn)有文獻，骨架層的螺旋纏繞角度對均勻外壓下臨界壓潰值影響較小，因此在研究中可以忽略不計。對骨架層承載能力影響較大的因素有材料非線性、初始橢圓度以及骨架層相鄰節(jié)之間的接觸方式”］。其中材料非線性指的是材料的彈塑特性，即當應(yīng)力超過材料的屈服極限時，應(yīng)力與應(yīng)變之間呈現(xiàn)一種非線性的增長關(guān)系，該因素對骨架層承載能力影響最大；骨架層在制造及安裝過程中，難免會有夾持等情況出現(xiàn)，從而導致其形狀改變，通常會呈現(xiàn)一種橢圓形狀。規(guī)定以橢圓的長徑與短徑之差和兩者之和的比值作為衡量其發(fā)生橢圓變形的標準，稱之為橢圓度。橢圓度越大，骨架層的承載能力越小。骨架層相鄰節(jié)之間的接觸方式對承載能力有所影響，例如綁定接觸或摩擦接觸以及摩擦接觸的摩擦因數(shù)等，同樣會影響骨架層的承載能力。本文所研究骨架層材料為316L不銹鋼，在0.2％初始橢圓度、相鄰節(jié)之間為摩擦接觸、摩擦因數(shù)為0.15的情況下分析其承載能力。

首先建立管道整體模型，如圖2a所示。在不考慮材料非線性、初始橢圓度及節(jié)間接觸的情況下，對其進行特征值屈曲分析，查看管道的壓潰形式，同時分析壓潰的臨界載荷值，與理論值進行對比，驗證理論模型的準確性。該分析中，骨架層的各節(jié)之間為綁定接觸，在管道軸向兩端限制軸向位移約束，同時打開弱彈簧開關(guān)，骨架層外表面施加1MPa的均勻載荷。骨架層及特征值屈曲分析結(jié)果如圖2b所示。

由圖2b可知，屈曲載荷系數(shù)為76.672，屈曲載荷等于屈曲載荷系數(shù)乘以施加載荷，由于施加的是單位載荷，所以屈曲載荷為76.672 MPa，即骨架層在該設(shè)定條件下的失穩(wěn)值為76.672 MPa，與理論值對比如表1所示。

由表1可知，骨架層特征值屈曲分析值與理論計算值很接近，驗證了用理論模型來計算無節(jié)間摩擦、完美圓狀態(tài)下彈性材料管道臨界失穩(wěn)值的準確性。

2.2 數(shù)值模擬分析

在工作實際中，上述理想條件比較少，僅靠理論模型不能滿足要求，可借助數(shù)值模擬來進行分析。結(jié)合實際及現(xiàn)有理論，材料彈塑性對骨架層的承載能力影響很大，對于塑性材料，當應(yīng)力超過比例極限后同樣采用上述分析顯然不合理。對于塑性分析有塑性極限載荷和塑性失穩(wěn)載荷2個參量來衡量含缺陷管道的承載能力。當考慮材料應(yīng)變硬化特性，載荷增加則變形增加，結(jié)構(gòu)發(fā)生無限制塑性變形時的載荷就稱之為塑性失穩(wěn)載荷；而假設(shè)材料為理想彈塑性材料時，管道所能承受的最大載荷稱之為塑性極限載荷，可用ASME規(guī)范中的2倍斜率準則來確定。

截止到目前，對于存在腐蝕缺陷的柔性管道承壓能力研究很少。本文研究了不同軸向、周向腐蝕長度及小范圍的徑向腐蝕深度對柔性管道骨架層承壓能力的影響。由于骨架層幾何尺寸及所受載荷的對稱性，為簡化模型，節(jié)省分析時間，只建立1／4模型進行后續(xù)分析。

骨架層材料為316L不銹鋼，在20℃溫度下，屈服強度為210 MPa，抗拉強度為490 MPa，彈性模量E為206GPa，泊松比為0.3，密度為7.98g／cm3，伸長率為40％。

在材料參數(shù)設(shè)置時，選用雙線性隨動強化模型，切線模量可通過下式進行簡化計算。

式中：T..為切線模量，MPa；Rm為抗拉強度，MPa；Y，為屈服強度，MPa；l為伸長率，無量綱。將以上參數(shù)帶入求得切線模量為700 MPa。

本文采用 ANSYS Workbench 軟件對非粘結(jié)柔 性管道的骨架層進行有限元分析，依據(jù)所研究的橢圓度要求建立有限元模型。由于骨架層截面的不規(guī)則性，采用掃掠的方式進行網(wǎng)格劃分，以骨架層的2個對稱截面分別為源面和目標面，首先對源面進行網(wǎng)格劃分，然后沿圓周向目標面掃掠，完成實體的網(wǎng)格劃分，劃分結(jié)果如圖3a所示。模型節(jié)與節(jié)之間接觸面設(shè)置為摩擦接觸，摩擦因數(shù)為0.15［15］；在模型上下兩截面施加對稱約束，限制其法向位移，兩側(cè)邊界施加軸向的位移約束，限制其軸向位移，外表面施加均勻壓力載荷，如圖3b所示。

塑性極限載荷又稱之為垮塌載荷，本文以此為參量來衡量骨架層的承載能力，以無腐蝕缺陷的骨架層分析為例，外表面施加50MPa均勻荷載，載荷步設(shè)置為1，步長為1000s，激活自動時間步，時間步類型選擇子步，初始載荷步、最小、最大載荷步分別為500、500和1000。總載荷隨載荷子步分步施加，在其逐步施加總變形發(fā)生突變，取該突變點位移及應(yīng)力云圖如圖4所示。經(jīng)估算該點施加載荷約為15.95 MPa，明顯低于特征值屈曲分析時的臨界壓潰值，表明對于塑性材料，材料非線性是分析其承載能力必不可少的考慮因素之一，同時也說明在塑性材料的比例極限之前沒有失穩(wěn)情況發(fā)生。

通過建立及改變骨架層內(nèi)表面不同尺寸的矩形凹陷，來模擬實際情況中管道不同的腐蝕缺陷，研究不同缺陷尺寸對管道失效壓力的影響，矩形缺陷形貌如圖5所示。紅色區(qū)域為腐蝕坑，XC為軸向長度方向，ZC為周向長度切線方向，YC為徑向深度方向。

所建立矩形腐蝕缺陷尺寸如下：

當軸向長度為3mm、徑向深度為0.5mm時，缺陷周向長度為17.732 5、35.465 1及53.197 6mm；當周向長度為35.4651mm、徑向深度為0.5mm時，缺陷軸向長度為3、5和7mm；當周向長度為35.4651mm、軸向長度為3mm時，缺陷徑向深度為0.3、0.5和0.7mm。

采用與無腐蝕缺陷模型同樣的步驟進行分析，根據(jù)分析結(jié)果，依據(jù)2倍彈性斜率準則［16］來確定骨架層的塑性極限載荷。在載荷—位移曲線中，斜率為彈性部分直線斜率2倍的直線與曲線的交點即為塑性極限載荷。采用該準則進行計算，結(jié)果如圖6所示。

將圖6依據(jù)2倍彈性斜率準則計算所得數(shù)據(jù)統(tǒng)計，結(jié)果如表2所示。

由上述分析可知，不同軸向、周向長度及徑向深度的腐蝕缺陷都能影響管道骨架層的失效壓力，但是影響甚微，可以忽略不計，即小范圍的局部腐蝕對骨架層承載能力的影響不大。由于骨架層與運輸介質(zhì)直接接觸，考慮其內(nèi)層緩慢的大面積均勻腐蝕對骨架層失效壓力的影響，因為骨架層內(nèi)表面與所運輸?shù)挠蜌馑旌辖橘|(zhì)充分接觸，所假設(shè)的均勻腐蝕邊界如圖7所示。

分別以均勻腐蝕深度δ為0.1、0.2、······、0.8mm的管道骨架層進行承載能力分析，同樣運用2倍彈性斜率準則來確定管道的塑性極限載荷，計算結(jié)果如圖8所示。

將圖8計算結(jié)果整理，結(jié)果如表3所示，其中Δp，為相鄰兩均勻腐蝕深度之間所對應(yīng)的失效壓力差值。

由表3及圖8可知，隨著腐蝕深度的增加，骨架層失效壓力下降速度整體上逐漸增加，呈現(xiàn)為拋物線趨勢，對其進行非線性的擬合，如圖9所示。擬合公式為：

式中：y為骨架層失效壓力，x為均勻腐蝕深度。

擬合數(shù)據(jù)及誤差如表4所示，經(jīng)檢驗，該模型決定系數(shù)為0.98，平均絕對誤差為0.1398，均方誤差為0.0414，誤差較小，擬合數(shù)據(jù)符合良好。

3腐蝕預(yù)測模型建立

海底管道下一時刻的腐蝕狀態(tài)只與當前時刻有關(guān)。馬爾可夫性質(zhì)如下［17—18］：一個隨機過程在給定現(xiàn)在及過去所有狀態(tài)情況下，其下一時刻所處狀態(tài)的條件概率分布僅取決于當前狀態(tài)。這與管道的腐蝕狀態(tài)轉(zhuǎn)移性質(zhì)相符合，因此可借用離散的馬爾可夫過程，即馬爾可夫鏈模型來對管道的腐蝕狀態(tài)進行預(yù)測。馬氏鏈表示如下：

建立馬爾可夫預(yù)測模型的關(guān)鍵在于構(gòu)建狀態(tài)轉(zhuǎn)移矩陣：

P表示系統(tǒng)在t時刻處于狀態(tài)i，下一時刻（t＋l）處于狀態(tài)j的概率；n為系統(tǒng)所有可能的狀態(tài)個數(shù)。預(yù)測模型的建立流程如圖10所示。

對于傳統(tǒng)剛性管道，根據(jù)SY／T 6151—2009《鋼制管道管體腐蝕損傷評價方法》，以管道的最大腐蝕深度為指標，將不同腐蝕程度的管道劃分為“輕、中、重、嚴重及穿孔”5個等級［19］，依次對應(yīng)5種狀態(tài)，計算各狀態(tài)之間的轉(zhuǎn)移概率，可構(gòu)建馬爾可夫狀態(tài)轉(zhuǎn)移矩陣，實現(xiàn)對管道所處腐蝕狀態(tài)的預(yù)測。而對于非粘結(jié)柔性管道骨架層復(fù)雜的幾何形狀，上述狀態(tài)劃分方式明顯不適用。依據(jù)不同腐蝕深度的骨架層承載能力分析，結(jié)合管道所處工況，可進行管道腐蝕狀態(tài)的劃分。以某管道2020—2021年工況為例，該管道輸送油氣水3項物質(zhì)，其中輸氣量（標準狀態(tài)）為77x104m3／d，輸油量、輸水量都為4.62m3／d，內(nèi)壓范圍為3.10～4.27 MPa，外壓范圍為12～15 MPa，出、入口溫度分別為25.0和39.2℃。

管道不同管段的骨架層承受不同的內(nèi)壓與外壓，當某管段內(nèi)壓最小、外壓最大時，骨架層承受最大的壓力載荷，為11.30 MPa；當某管段內(nèi)壓最大、外壓最小時，骨架層承受最小的壓力載荷，為7.73 MPa，因此該管道可能承受的壓力載荷范圍為7.73～11.3 MPa。根據(jù)式（7）計算可得，該失效壓力區(qū)間所對應(yīng)的骨架層均勻腐蝕深度為0.609～0.854mm，即當均勻腐蝕深度小于0.5mm時，管道處于安全狀態(tài)，定義腐蝕狀態(tài)為輕度腐蝕；當腐蝕深度處于0.5～0.6mm之間時，雖然管道失效壓力大于工況所施加的載荷，但該階段的失效壓力隨腐蝕深度變化較快，有即將失效的風險，定義腐蝕狀態(tài)為中度腐蝕；當均勻腐蝕深度處于0.60～0.85mm之間時，管道中某些位置骨架層失效壓力小于工況所施加的壓力載荷，故定義腐蝕狀態(tài)為重度腐蝕；當繼續(xù)腐蝕，腐蝕深度大于0.85mm時，管道骨架層失效壓力小于工況所施加的最小載荷，管道完全失效，定義該階段腐蝕狀態(tài)為失效。即該工況腐蝕狀態(tài)劃分如下：①輕度腐蝕，均勻腐蝕深度小于0.5mm；②中度腐蝕，均勻腐蝕深度在0.5～0.6mm之間；③重度腐蝕，均勻腐蝕深度在0.60～0.85mm之間；④失效，均勻腐蝕深度大于0.85mm。

腐蝕過程中，若不采取措施，劣化狀態(tài)會越發(fā)嚴重，隨著時間的延長，或停留在當前狀態(tài)，或越來越差，因此狀態(tài)④，即失效狀態(tài)只能停留至該狀態(tài)，不能向其他狀態(tài)轉(zhuǎn)移，其值為1。根據(jù)上述劃分的4種腐蝕狀態(tài)，可得馬爾可夫狀態(tài)轉(zhuǎn)移矩陣如下：

式（10）中，P（i，j＝1，2，3，4）為狀態(tài)轉(zhuǎn)移概率，表示當前時刻為i狀態(tài)，下一時刻為j狀態(tài)的概率。已知初始狀態(tài)為T（ω）＝［a，b，c，d］，其中，a、b、c、d分別為初始時刻管道處于各狀態(tài)的概率，則n年后管道所處腐蝕狀態(tài)為T（m）＝T（o）xP（）。至此，完成管道所處腐蝕狀態(tài)的預(yù)測。

4 結(jié)論

基于ANSYS Workbench軟件對非粘結(jié)柔性管 道骨架層進行了特征值屈曲分析，并與理論值作對比。另外，在不考慮螺旋纏繞角度、初始橢圓度為0.2％、考慮材料彈塑性及節(jié)間摩擦接觸條件下，對腐蝕缺陷形貌及尺寸對骨架層承載能力的影響規(guī)律進行了分析，并借助分析結(jié)果構(gòu)建了馬爾可夫預(yù)測模型，得出如下主要結(jié)論：

（1）結(jié)合鐵木辛柯彈性穩(wěn)定理論中的圓環(huán)撓曲線微分方程，根據(jù)單位長度彎曲剛度相等將骨架層等效成截面為矩形的圓環(huán)，計算其臨界失穩(wěn)值；建立骨架層有限元模型并對其進行特征值屈曲分析，將分析結(jié)果與理論計算結(jié)果進行對比，誤差為3.02％，驗證了理論模型對于無節(jié)間摩擦、完美圓狀態(tài)下彈性材料管道的承載能力計算的可行性。

（2）對于塑性材料，采用理論模型或特征值屈曲方式來進行分析顯然不合理，本文以塑性極限載荷為參量來衡量骨架層的承載能力。通過分析矩形腐蝕缺陷的周向、軸向長度及徑向深度對骨架層承載能力的影響，發(fā)現(xiàn)小范圍的局部腐蝕缺陷對其影響不是很大，可以忽略，考慮到骨架層內(nèi)層與運輸介質(zhì)直接接觸，分析了骨架層內(nèi)層均勻腐蝕對承載能力的影響，分析結(jié)果顯示該種腐蝕影響較大。

（3）分析均勻腐蝕深度對管道失效壓力的影響規(guī)律發(fā)現(xiàn)，失效壓力曲線呈拋物線形式，對其進行非線性擬合，得到均勻腐蝕深度與骨架層失效壓力的關(guān)系式，經(jīng)檢驗，該關(guān)系式計算所得結(jié)果與仿真結(jié)果吻合較好。

（4）在管道腐蝕狀態(tài)預(yù)測時，由于現(xiàn)有的腐蝕狀態(tài)劃分方式對非粘結(jié)柔性管骨架層的不適用性，基于分析結(jié)果，建立了針對柔性管的腐蝕狀態(tài)劃分新方式，并且建立了馬爾可夫預(yù)測模型進行預(yù)測。

參考文獻

[1]羅懿.基于ANSYS Workbench有限元法的外腐蝕管 道失效壓力研究［J］．石油化工設(shè)備，2019，48（2）：7-11.

LUO Y. Study on the external corrosion failure pressure of pipeline based on ANSYS Workbench finite element method [J].Petro-Chemical Equipment，2019，48 （2）：7-11.

[2]張旭昀，孫磊，劉天宇，等.基于ANSYS Work—bench有限元法含內(nèi)腐蝕缺陷套管應(yīng)力模擬研究［J］．化工機械，2015，42（1）：97—100．

ZHANG X Y，SUN L，LIU T Y，et al. Stress analysis of corroded oil casing with internal defects base on AN- SYS Workbench [J].Chemical Engineering & Machin- ery，2015，42（1）：97-100.

[3]杜捷.基于ANSYS的油氣腐蝕管道的受力分析［J］．質(zhì)量技術(shù)監(jiān)督研究，2018（5）：33—36，40．

DU J. Stress analysis of oil and gas corroded pipeline based on ANSYS [J].Quality and Technical Supervi- sion Research，2018（5）：33-36，40.

[4]楊洋，楊茜，賀美玲，等.利用ANSYS進行油氣管道外腐蝕剩余強度分析［J］.內(nèi)蒙古石油化工，2015，41（10）：19-22.

YANG Y，YANG Q，HE M L， et al. Residual strength analysis of extemal corrosion of oil and gas pipelines by ANSYS [J].Inner Mongolia Petrochemi- cal Industry，2015，41（10）：19-22.

[5]崔銘偉，王媛媛，賀杰，等.相互作用腐蝕管道剩余強度評價方法對比研究［J］．表面技術(shù)，2016，45（8）：56-62.

CUI M W，WANG Y Y，HE J，et al. Comparison study on assessment methods for residual strength of in- teraction corroded pipeline [J].Surface Technology， 2016，45（8）：56-62.

[6]劉洲，雍歧衛(wèi)，全琪.基于ANSYS腐蝕管道組合缺陷剩余強度評價［J］．自動化與儀器儀表，2016（2）：59-61.

LIU Z， YONG Q W， QUAN Q. Residual strength eval- uation of combined defects of corroded pipeline based on ANSYS [J]. Automation & Instrumentation，2016 （2）：59-61.

[7]王戰(zhàn)輝，馬向榮，白雪，等.油氣管道組合缺陷等效應(yīng)力［J］．化工科技，2020，28（1）：8—11．

WANG Z H，MA X R， BAI X， et al. Equivalent stress of combined defects in oil and gas pipelines [J]. Science&Technology in Chemical Industry，2020，28 （1）：8-11.

[8]李偉民.海洋非粘結(jié)柔性管壓潰理論分析及數(shù)值模擬研究［D］．青島：中國海洋大學，2015．

LI W M. Theoretical analysis and numerical simulation on collapse for marine unbonded flexible pipe [D]. Qingdao： Ocean University of China，2015.

[9]徐孝軒，佘振興，黃玫，等.柔性立管內(nèi)部骨架層徑向承載能力計算研究［J］．石油機械，2019，47（10）：33-37.

XU X X， SHE Z X， HUANG M，et al. Calculation of radial bearing capacity of internal carcass layer of flexi- ble riser [J].China Petroleum Machinery，2019，47 （10）：33-37.

[10]TANG M G，LU Q Z，YAN J，et al. Buckling col-lapse study for the carcass layer of flexible pipes using a strain energy equivalence method [J].Ocean Engi- neering，2016，111：209-217.

[11]CUAMATZI-MELENDEZ R， CASTILLO-HERNANDEZ O，VAZQUEZ-HERN?NDEZ A O， et al. Finite ele- ment and theoretical analyses of bisymmetric collapses in flexible risers for deepwaters developments [J].O- cean Engineering，2017，140：195-208.

[12]傅一欽，余楊，余建星，等.考慮層間接觸的非黏結(jié)柔性立管骨架層有限元數(shù)值模擬研究［J］.中國造船，2019，60（4）：100—107．

FU YQ，YU Y，YU J X，et al. Numerical finite el- ement simulation of unbonded flexible riser skeleton with consideration of Layer-to-Layer contact [J].Shipbuilding of China，2019，60（4）：100-107.

[13]徐孝軒，李文博，程子云，等.非黏結(jié)柔性立管內(nèi)部骨架層失效研究進展［J］．石油機械，2018，46（8）：48-54，59.

XU XX，LIWB，CHENG Z Y， et al. Advances on failure of carcass layer in un-bonded flexible riser [J]. China Petroleum Machinery， 2018， 46（8）：48- 54，59.

[14]NETO A G， DE ARRUDA MARTINS C. A compara- tive wet collapse buckling study for the carcass layer of flexible pipes [J].Joumal of Offshore Mechanics and Arctic Engineering，2012，134（3）：031701.

[15]任少飛，唐文勇，薛鴻祥.防摩擦層對非粘結(jié)柔性立管截面力學性能影響［J］.華中科技大學學報（自然科學版），2015，43（12）：60—64．

REN S F， TANG W Y， XUE H X. Impact of anti- friction layer on the cross-sectional mechanical proper- ties of unbonded flexible risers [J].Joumal of Hua- zhong University of Science and Technology（Natural Science Edition），2015，43（12）：60-64.

[16]陳孫藝.屈服極限和流變應(yīng)力以及塑性極限載荷的確定方法綜述［C］／／第三屆全國管道技術(shù)學術(shù)會議壓力管道技術(shù)研究進展精選集.烏魯木齊：機械工業(yè)出版社，2006：66—72．

CHEN S Y. Method to obtain the yield limit and the flow stress and the plastic limit load of structure[C]// Selected Collection of Research Progress of Pressure Pipeline Technology at the Third National Pipeline Technology Academic Conference.Urumqi：China Machine Press，2006：66-72.

[17]VALOR A， CALEYO F， ALFONSO L， et al. Mark- ov chain models for the stochastic modeling of pitting corrosion [J].Mathematical Problems in Engineer- ing，2013，2013：108386.

[18]何成剛.馬爾科夫模型預(yù)測方法的研究及其應(yīng)用［D］．合肥：安徽大學，2011．

HE C G.The research of Markov model forecasing method and application [D].Hefei：Anhui Universi- ty，2011.

[19]袁賡.油氣管道的腐蝕及預(yù)測研究［D］.大連：大連理工大學，2011．

YUAN G.Research of corrosion and its prediction for oil-gas steel pipelines [D].Dalian：Dalian Universi- ty of Technology，2011.

第一作者簡介：魯中歧，生于1997年，2020年畢業(yè)于青島理工大學機械制造及其自動化專業(yè)，現(xiàn)為在讀碩士研究生，研究方向為石油裝備與工程。地址：（266580）山東省青島市。E-mail：luzq777＠163．com。

通信作者：崔俊國，E-mail：cuijunguo＠upc.edu.cn。

收稿日期：2022—07—19

（本文編輯 王剛慶）