CN110822294B - 一种含腐蚀缺陷的海底管道结构弯曲失效评估方法 - Google Patents

一种含腐蚀缺陷的海底管道结构弯曲失效评估方法 Download PDF

Info

Publication number
CN110822294B
CN110822294B CN201810917814.XA CN201810917814A CN110822294B CN 110822294 B CN110822294 B CN 110822294B CN 201810917814 A CN201810917814 A CN 201810917814A CN 110822294 B CN110822294 B CN 110822294B
Authority
CN
China
Prior art keywords
pipeline
corrosion
bending moment
plastic
plastic limit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810917814.XA
Other languages
English (en)
Other versions
CN110822294A (zh
Inventor
杨冬平
王强
丛军
郭爱洪
王伟斌
牛更奇
支景波
朱丽国
高莎莎
马文峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Technology Inspection Center of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Technology Inspection Center of Sinopec Shengli Oilfield Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Technology Inspection Center of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201810917814.XA priority Critical patent/CN110822294B/zh
Publication of CN110822294A publication Critical patent/CN110822294A/zh
Application granted granted Critical
Publication of CN110822294B publication Critical patent/CN110822294B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00—Protection or supervision of installations
    • F17D5/02—Preventing, monitoring, or locating loss
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N17/00—Investigating resistance of materials to the weather, to corrosion, or to light
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/20—Investigating strength properties of solid materials by application of mechanical stress by applying steady bending forces

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Immunology (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Environmental Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Ecology (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种含腐蚀缺陷的海底管道结构弯曲失效评估方法，具体为：检测管道最大尺寸腐蚀缺陷位置及腐蚀尺寸；根据管道检测数据、工况条件及环境条件计算管道承受的波浪及海流载荷，最终确定出被测管道所承受的最大弯矩位置和最大弯矩值；推导出含腐蚀缺陷管道的塑性极限弯矩公式，按照该公式可计算被测管道的塑性极限承载弯矩；将计算得到的最大弯矩值与塑性极限弯矩进行比较，在考虑安全系数的前提下判断被测管道是否失效。本发明的有益效果是：本发明能够对含有任意位置的小范围长轴向腐蚀缺陷的管道，进行整体弯曲强度计算评估，省去了有限元建模及分析过程，降低了计算难度，提高了计算的效率，同时又满足结果精度要求。

Description

一种含腐蚀缺陷的海底管道结构弯曲失效评估方法
技术领域
本发明涉及海底管道失效评估技术领域，特别涉及一种含腐蚀缺陷的海底管道结构弯曲失效评估方法。
背景技术
近些年来，我国对石油、天然气等资源的需求量不断上升。石油、天然气等资源的开发和利用也越来越重要。海底管道作为油气资源开发、生产以及运输的“生命线”，确保管道结构安全有着极其重要的作用。由于海洋环境复杂，海底管道常常会因淘蚀和海水冲刷而处于悬空状态。当海底管道外防腐层出现损坏时，海水会沿防腐层损坏的方向发生侵蚀，易在管道外侧任意位置形成长腐蚀缺陷(环向尺寸小而轴向长度大的腐蚀缺陷，理想状态为轴向腐蚀长度不小于
Figure GDA0001867257710000011
其中，R是管道半径，t是管道厚度)。管道腐蚀后，会造成管道局部变薄，极限承载能力降。管道的腐蚀是造成海底管道发生破坏的重要原因之一。另外，海底管道往往会遭受滑坡、河床变形等自然灾害，以及因淘蚀和海水冲刷导致悬空，此时，弯矩载荷便成为导致海底管道发生破坏的主要因素。
准确评估腐蚀管道的极限承载力或者极限弯矩(通过极限承载力可以计算极限弯矩，这两者可以互相转换)是确保管道正常工作的前提。在试验研究和理论模拟的基础上，国内外学者对腐蚀管道的极限承载能力进行了大量的研究。20世纪70年代，美国机械工程师协会提出了腐蚀管道极限内压载荷规范：ASME B31G，仍是目前最为广泛运用于管道极限承载能力评估的规范。考虑B31G在评估腐蚀管道极限内压载荷过于保守，Kiefer、Vieth对B31G中的流应力、Folias系数等参数进行了改进，进而提出了修正版的极限内压载荷规范：MB31G(0.8dl)。21世纪初期，挪威船级社(DNV)在考虑轴压和弯曲载荷的前提下，通过对海底腐蚀管道进行了一系列的数值模拟和试验研究，提出了DNV-RP-F101腐蚀管道极限承载能力的评估规范。2003年，Choi等对X65天然气腐蚀管道进行了爆破试验和数值模拟，回归出了含矩形腐蚀缺陷和椭圆腐蚀缺陷管道的极限内压公式。Zhu和Leis等应用有限应变理论和变形失稳理论，提出了适用于完好管道的塑性失效模型，并将模型扩展到腐蚀缺陷管道中。
但是，目前对于海底管道在受弯作用下塑性失效的极限承载能力的评估方法，多是将腐蚀缺陷设定在弯曲应力最大的位置，从而得到含腐蚀海底管道塑性失效临界荷载的下限。但实际海底管道表面腐蚀位置存在任意性，即腐蚀位置不一定位于弯曲应力最大的部位，因此研究包含任意位置海底管道受弯塑性极限承载能力的计算方法更具有实际应用价值，而目前要进行任意腐蚀位置缺陷的失效判定只能通过建立管道有限元模型进行模拟分析实现。
发明内容
针对传统评估方法是通过将腐蚀缺陷设定于管道弯曲最大位置，从而得出腐蚀海的管道塑性失效临界荷载下限，存在评测结果不准确、实际应用价值不高的问题，本发明提供了一种含腐蚀缺陷的海底管道结构弯曲失效评估方法，其包含任意位置海底管道受弯塑性极限承载弯矩的计算方法，而不仅限于腐蚀缺陷仅存在于管道弯曲最大位置的情况。
本发明提供了一种含腐蚀缺陷的海底管道结构弯曲失效评估方法，所述评估方法具体为：
步骤S1：选择一服役管道作为被测管道，检测被测管道最大尺寸腐蚀缺陷的位置及腐蚀尺寸；
步骤S2：根据步骤S1的检测数据，进行海浪荷载和海流载荷计算，确定被测管道所承受的弯矩大小及弯曲方向，从而确定出被测管道所承受的最大弯矩位置和最大弯矩值；
步骤S3：推导出含腐蚀缺陷管道的塑性极限弯矩公式，并依据步骤S1的检测结果按照塑性极限弯矩公式，计算出被测管道的塑性极限承载弯矩；
步骤S4：将步骤S2计算得到的最大弯矩值与步骤S3计算得到的塑性极限承载弯矩进行比较，若最大弯矩值等于或大于塑性极限承载弯矩，则被测管道失效；若最大弯矩值小于塑性极限承载弯矩，则被测管道未失效；在实际评估中需要考虑安全系数时，用许用塑性极限弯矩代替步骤S3计算得到的塑性极限弯矩进行比较，此时许用塑性极限弯矩取S3计算得到的塑性极限弯矩除以安全系数n后的值，安全系数n依据管道规范要求确定。
其中，所述步骤S1中，所述腐蚀缺陷的位置包括腐蚀坑位置角α；所述腐蚀尺寸包括腐蚀坑的腐蚀宽度2Rθ、腐蚀深度h和腐蚀长度L。
所述步骤S2中采用海浪谱进行所述海浪载荷的计算，同时要考虑海流载荷的联合作用。
其中，所述步骤S3中被测管道的塑性极限承载弯矩的计算过程，具体如下：
步骤S301：依据被测管道的腐蚀缺陷的位置及腐蚀尺寸，确定塑性中性轴的位置，并按照公式(1)和公式(2)计算得被测管道的受拉区面积AI和受压区面积AΠ：
AΙ＝2Rmtβ (1)
AΠ＝2Rmt(π-β)-2Rhθ (2)
其中：
Figure GDA0001867257710000031
为平均半径；R为管道外壁半径；t为被测管道的厚度，t远小于R；h为腐蚀深度，θ为腐蚀宽度的半角；
步骤S302：依据被检测管道的腐蚀缺陷的位置及腐蚀尺寸，能够得出塑性中性轴的角度β、腐蚀缺陷的位置角α、塑性中性轴与水平轴的夹角γ，进而按照公式(3)-(6)得出受拉区区域和受压区区域的形心坐标(yc1，zc1)和(yc2，zc2)
Figure GDA0001867257710000032
Figure GDA0001867257710000033
Figure GDA0001867257710000034
Figure GDA0001867257710000035
步骤S303：若被测管道处于塑性极限承载状态，则应满足以下两个条件：
a)：受拉区与受压区的形心连线平行于z轴，即yc1＝yc2；
b)：受拉区面积等于受压区面积，即A＝AI＝AΠ；
由以上两个条件，计算可得：
Figure GDA0001867257710000036
Figure GDA0001867257710000041
步骤S304：依据步骤S303的计算结果，进而得到被测管道在塑性极限承载状态下的极限弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβcosγ-2Rm 2hsinθcosα) (9)
当缺陷相对z轴对称时，即α＝γ＝0，其塑性极限承载弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβ-2Rm 2hsinθ) (10)
步骤S305：将被测管道的腐蚀缺陷位置及腐蚀尺寸，代入到公式(9)或公式(10)，即可得到被测管道的塑性极限承载弯矩。
为了实现上述发明目的，本发明还提供了一种海底管道塑性极限承载弯矩的计算方法，所述计算方法具体为：
步骤S1：依据被测管道的腐蚀缺陷的位置及腐蚀尺寸，确定塑性中性轴的位置，并按照公式(1)和公式(2)计算得被测管道的受拉区面积AI和受压区面积AΠ：
AΙ＝2Rmtβ (1)
AΠ＝2Rmt(π-β)-2Rhθ (2)
其中：
Figure GDA0001867257710000042
为平均半径；R为管道外壁半径；t为被测管道的厚度，t远小于R；h为腐蚀深度，θ为腐蚀宽度的半角；
步骤S2：依据被检测管道的腐蚀缺陷的位置及腐蚀尺寸，能够得出塑性中性轴的角度β、腐蚀缺陷的位置角α、塑性中性轴与水平轴的夹角γ，进而按照公式(3)-(6)得出受拉区区域和受压区区域的形心坐标(yc1，zc1)和(yc2，zc2)
Figure GDA0001867257710000043
Figure GDA0001867257710000044
Figure GDA0001867257710000045
Figure GDA0001867257710000046
步骤S3：若被测管道处于塑性极限承载状态，则应满足以下两个条件：
a)：受拉区与受压区的形心连线平行于z轴，即yc1＝yc2；
b)：受拉区面积等于受压区面积，即A＝AI＝AΠ；
由以上两个条件，计算可得：
Figure GDA0001867257710000051
Figure GDA0001867257710000052
步骤S4：依据步骤S303的计算结果，进而得到被测管道在塑性极限承载状态下的极限弯矩承载弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβcosγ-2Rm 2hsinθcosα) (9)
当缺陷相对z轴对称时，即α＝γ＝0，其塑性极限承载弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβ-2Rm 2hsinθ) (10)
步骤S5：将被测管道的腐蚀缺陷位置及腐蚀尺寸，代入到公式(9)或公式(10)，即可得到被测管道的塑性极限承载弯矩。
其中，所述步骤S1中，所述腐蚀缺陷的位置包括腐蚀坑位置角α；所述腐蚀尺寸包括腐蚀坑的腐蚀宽度2Rθ、腐蚀深度h和腐蚀长度L。
本发明的有益效果是：本发明能够针对管道任意位置的小范围长腐蚀缺陷，计算出在纯弯矩载荷下腐蚀管道的极限承载弯矩，并与管道所承受的最大弯矩位置和最大弯矩值进行对比，对含腐蚀管道进行整体强度评测，省去了含腐蚀缺陷的管状结构弯曲失效评估的有限元建模计算过程，降低了工程应用的计算难度，提高了计算的效率，同时又满足了结果精度要求。
附图说明
图1为本发明具体实施方式中含腐蚀缺陷管道的有限元网格图。
图2为本发明具体实施方式中有限元分析破坏图。
图3为本发明具体实施方式中有限元弯矩与端部转角曲线图。
图4为本发明对比试验二中试件尺寸及腐蚀缺陷尺寸。
具体实施方式
由于海洋环境复杂，海底管道常常会因淘蚀和海水冲刷而处于悬空状态，此时弯矩载荷便成为海底管道发生破坏的主要因素。当海底管道外防腐层出现损坏时，海水会沿防腐层损坏的方向发生侵蚀，很容易在管道外侧任意位置形成长腐蚀缺陷。本申请针对任意位置的小范围长腐蚀缺陷，进行了数值模拟以及有限元分析推导出了在纯弯矩载荷下的极限弯矩计算公式，并与有限元分析结果进行对比，验证该计算公式的准确性。
为能清楚说明本方案的技术特点，下面通过具体实施方式，对本方案进行阐述。
实施例1失效评估方法
本发明实施例提供了一种含腐蚀缺陷的海底管道结构弯曲失效评估方法，具体为：
步骤S1：选择一现役管道作为被测管道，检测被测管道最大尺寸腐蚀缺陷的位置及腐蚀尺寸；
步骤S2：根据被测管道的服役环境状态参数和步骤S1的检测数据，进行海浪荷载和海流载荷的合力计算，确定被测管道所承受的弯矩大小及弯曲方向，从而确定出被测管道所承受的最大弯矩位置和最大弯矩值；
步骤S3：推导出含腐蚀缺陷管道的塑性极限弯矩公式，并依据步骤S1的检测结果按照塑性极限弯矩公式，计算出被测管道的塑性极限弯矩；
步骤S4：将步骤S2计算得到的最大弯矩值与步骤S3计算得到的塑性极限承弯矩进行比较，若最大弯矩值等于或大于塑性极限承弯矩，则被测管道失效；若最大弯矩值小于塑性极限弯矩，则被测管道未失效，本案例为验证结果取安全系数为1。
其中，步骤S1中，腐蚀缺陷的位置包括腐蚀坑位置角α；腐蚀尺寸包括腐蚀坑的腐蚀宽度2Rθ、腐蚀深度h和腐蚀长度L；步骤S2中采用海浪谱进行所述海浪载荷和海流载荷的计算，海浪载荷和海流载荷的计算方法为现有技术，也不是本申请研究的重点，在此不再详细赘述。
其中，所述步骤S3中被测管道的塑性极限弯矩的计算过程，具体如下：
步骤S301：依据被测管道的腐蚀缺陷的位置及腐蚀尺寸，确定塑性中性轴的位置，并按照公式(1)和公式(2)计算得被测管道的受拉区面积AI和受压区面积AΠ：
AΙ＝2Rmtβ (1)
AΠ＝2Rmt(π-β)-2Rhθ (2)
其中：
Figure GDA0001867257710000071
为平均半径；R为管道外壁半径；t为被测管道的厚度，t远小于R；h为腐蚀深度，θ为腐蚀宽度的半角；
步骤S302：依据被检测管道的腐蚀缺陷的位置及腐蚀尺寸，能够得出塑性中性轴的角度β、腐蚀缺陷的位置角α、塑性中性轴与水平轴的夹角γ，进而按照公式(3)-(6)得出受拉区区域和受压区区域的形心坐标(yc1，zc1)和(yc2，zc2)
Figure GDA0001867257710000072
Figure GDA0001867257710000073
Figure GDA0001867257710000074
Figure GDA0001867257710000075
步骤S303：若被测管道处于塑性极限承载状态，则应满足以下两个条件：
a)：受拉区与受压区的形心连线平行于z轴，即yc1＝yc2；
b)：受拉区面积等于受压区面积，即A＝AI＝AΠ；
由以上两个条件，计算可得：
Figure GDA0001867257710000076
Figure GDA0001867257710000077
步骤S304：依据步骤S303的计算结果，进而得到被测管道在塑性极限承载状态下的极限弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβcosγ-2Rm 2hsinθcosα) (9)
当缺陷相对z轴对称时，即α＝γ＝0，其塑性极限承载弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβ-2Rm 2hsinθ) (10)
步骤S305：将被测管道的腐蚀缺陷位置及腐蚀尺寸，代入到公式(9)或公式(10)，即可得到被测管道的塑性极限弯矩。
实施例2塑性极限弯矩的计算方法
本发明实施例提供了一种海底管道塑性极限弯矩的计算方法，具体为：
步骤S1：依据被测管道的腐蚀缺陷的位置及腐蚀尺寸，确定塑性中性轴的位置，并按照公式(1)和公式(2)计算得被测管道的受拉区面积AI和受压区面积AΠ：
AΙ＝2Rmtβ (1)
AΠ＝2Rmt(π-β)-2Rhθ (2)
其中：
Figure GDA0001867257710000081
为平均半径；R为管道外壁半径；t为被测管道的厚度，t远小于R；h为腐蚀深度，θ为腐蚀宽度的半角；
步骤S2：依据被检测管道的腐蚀缺陷的位置及腐蚀尺寸，能够得出塑性中性轴的角度β、腐蚀缺陷的位置角α、塑性中性轴与水平轴的夹角γ，进而按照公式(3)-(6)得出受拉区区域和受压区区域的形心坐标(yc1，zc1)和(yc2，zc2)
Figure GDA0001867257710000082
Figure GDA0001867257710000083
Figure GDA0001867257710000084
Figure GDA0001867257710000085
步骤S3：若被测管道处于塑性极限承载状态，则应满足以下两个条件：
a)：受拉区与受压区的形心连线平行于z轴，即yc1＝yc2；
b)：受拉区面积等于受压区面积，即A＝AI＝AΠ；
由以上两个条件，计算可得：
Figure GDA0001867257710000086
Figure GDA0001867257710000087
步骤S4：依据步骤S303的计算结果，进而得到被测管道在塑性极限承载状态下的极限弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβcosγ-2Rm 2hsinθcosα) (9)
当缺陷相对z轴对称时，即α＝γ＝0，其塑性极限承载弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2tsinβ-2Rm 2hsinθ) (10)
步骤S5：将被测管道的腐蚀缺陷位置及腐蚀尺寸，代入到公式(9)或公式(10)，即可得到被测管道的塑性极限弯矩。
其中，所述步骤S1中，腐蚀缺陷的位置包括腐蚀坑位置角α；腐蚀尺寸包括腐蚀坑的腐蚀宽度2Rθ、腐蚀深度h和腐蚀长度L。
对比试验一
传统的评估方法，需要进行包含腐蚀缺陷的管道弯曲失效评估的有限元分析。本实验对同一管道分别采用有限元分析和本申请的计算公式进行对比。其中，实验用管道参数如下：材料型号为API 5L X56，直径为325mm，管道厚度为14mm，长度为1100mm；腐蚀缺陷为外腐蚀，腐蚀缺陷的深度为7mm，腐蚀的宽度尺寸大小与位置参数，如表1所示。
表1腐蚀的宽度尺寸大小与位置参数
序号 θ(°) α(°) 序号 θ(°) α(°)
1 15° 0° 10 30° 45°
2 15° 15° 11 30° 60°
3 15° 30° 12 30° 75°
4 15° 45° 13 45° 0°
5 15° 60° 14 45° 15°
6 15° 75° 15 45° 30°
7 30° 0° 16 45° 45°
8 30° 15° 17 45° 60°
9 30° 30° 18 45° 75°
本实验中在对有缺陷的海底管道的有限元分析中，采用Abaqus软件进行分析模拟。在模拟过程中，采用三维8结点线性单元对管道进行分析。管道网格划分采用不同应力梯度对应不同网格密度的方法。在缺陷处采用小尺寸单元形成精细网格，在其他区域采用尺寸较大的单元形成粗糙网格。这个方法得到一个典型的含腐蚀缺陷管道的有限元网格如图1所示。有限元模型管道参数如前所述，模型约束采用对称约束，同时在一端部设置参考点，参考点与端部截面之间采用耦合约束，弯矩通过参考点施加在管道模型上。
腐蚀管道在纯弯作用下存在着三种主要失效模式，分别为椭圆化、局部屈曲和管壁开裂；当腐蚀缺陷处于受拉状态，腐蚀深度较浅管道容易发生椭圆化失效，而腐蚀深度较深管道容易发生断裂失效，且失效模式同腐蚀长度有关；当腐蚀缺陷处于受压状态时，腐蚀深度较浅管道容易发生椭圆化失效，腐蚀深度较深管道容易发生局部屈曲，且失效模式同腐蚀长度有关。
其中，管道中腐蚀模型18的有限元破坏图如图2所示，由图2能够看出，缺陷管道在纯弯矩载荷作用下，当达到屈服应力发生破坏时，其最终的破坏形式为腐蚀缺陷部分的外凸破坏以及缺陷周围部分发生局部的屈曲破坏。腐蚀模型18的弯矩与端部转角曲线图如图3所示，由有限元弯矩与端部转角曲线即可得到有限元计算的极限弯矩值。
分别采用有限元分析和本申请中的计算公式对管道腐蚀模型1-18依次进行分析，得到腐蚀模型1-18的极限弯矩值，如表2和图3所示。
表2有限元分析与公式计算结果对比
Figure GDA0001867257710000101
Figure GDA0001867257710000111
由表2的数据能够看出，当缺陷大小不变时，随着决定缺陷位置的参数α的增大，极限承载弯矩在逐渐增大，但是增大的幅度不大。当缺陷位置不变时，随着缺陷大小的增大，极限承载弯矩有着明显的减小。由图3能够看出，有限元的计算结果曲线与公式计算结果曲线的变化趋势基本相同，两者结果的最大误差为腐蚀模型1，其误差大小为1.37％，最小误差为腐蚀模型3，其最小误差为0.05％。可见，采用本申请公式计算出的极限弯矩基本上与有限元计算的结果相同，但公式计算出的结果相对于有限元结果来说偏危险，即采用本申请的计算公式以及含有本申请计算公式的评估方法，更能准确且有效的预知管道的失效情况。
对比试验二
本试验选取单层含腐蚀缺陷的海底管道试件进行弯曲试验测试，试件尺寸及腐蚀缺陷尺寸分别如表4和图4所示。本试验对试件弯曲试验测试均在四立柱压力试验机上进行，压力机量程5000kN，管道在轴压弯曲变形过程中用位移计测量位移大小，位移计量程最大可达200mm；施加在管道上的荷载大小通过加载和控制系统的荷载盘读数。加载过程中，分阶段加载，到达加载值后，保持荷载不动，读取荷载和对应的位移值，从而使荷载和位移读数保持一致。测试记录数据如表4所示。
表3试验件管道的基本参数
Figure GDA0001867257710000112
表4试件弯曲试验测试记录数据表
Figure GDA0001867257710000113
Figure GDA0001867257710000121
通过表3和表4的数据可以知道，试件轴压极限载荷为3860kN，弯曲极限载荷为1382kN。
通过本发明的评估方法分别计算弯曲极限载荷，并进行误差计算，具体如下：
(1)承载力计算：
中性轴计算公式
Figure GDA0001867257710000122
弯矩计算公式为
Figure GDA0001867257710000132
转化为力为F＝2M/0.55＝2×378.9/0.55＝1377.8kN
(2)实验测试值为F1＝1382kN；
(3)误差计算，公式为
Figure GDA0001867257710000131
由上述内容能够看出，本申请的评估方法与实际测量基本一致，误差很小能够满足工程应用对结果精度的要求。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种含腐蚀缺陷的海底管道结构弯曲失效评估方法，其特征在于，所述评估方法具体为：
步骤S1：选择一现役管道作为被测管道，检测被测管道最大尺寸腐蚀缺陷的位置及腐蚀尺寸；
步骤S2：根据步骤S1的检测数据，进行海浪荷载和海流载荷计算，确定被测管道所承受的弯矩大小及弯曲方向，从而确定出被测管道所承受的最大弯矩位置和最大弯矩值；
步骤S3：推导出含腐蚀缺陷管道的塑性极限弯矩公式，并依据步骤S1的检测结果按照塑性极限弯矩公式，计算出被测管道的塑性极限弯矩；
步骤S4：将步骤S2计算得到的最大弯矩值与步骤S3计算得到的塑性极限弯矩进行比较，若最大弯矩值等于或大于塑性极限弯矩时，则被测管道失效；若最大弯矩值小于塑性极限弯矩，则被测管道未失效；
在实际评估中需要考虑安全系数时，用许用塑性极限弯矩代替步骤S3计算得到的塑性极限弯矩进行比较，此时许用塑性极限弯矩取步骤S3计算得到的塑性极限弯矩除以安全系数n后的值，安全系数n依据管道规范要求确定；
其中，所述步骤S3中被测管道的塑性极限弯矩的计算过程，具体如下：
步骤S301：依据被测管道的腐蚀缺陷的位置及腐蚀尺寸，确定塑性中性轴的位置，并按照公式(1)和公式(2)计算得被测管道的受拉区面积AI和受压区面积AΠ：
AΙ＝2Rmtβ (1)
AΙ＝2Rmt(π-β)-2Rhθ (2)
其中：
Figure FDA0003056440790000011
为平均半径；R为管道外壁半径；t为被测管道的厚度，t远小于R；h为腐蚀深度，θ为腐蚀宽度的半角；
步骤S302：依据被检测管道的腐蚀缺陷的位置及腐蚀尺寸，能够得出塑性中性轴的角度β、腐蚀缺陷的位置角α、塑性中性轴与水平轴的夹角γ，进而按照公式(3)-(6)得出受拉区区域和受压区区域的形心坐标(yc1，zc1)和(yc2，zc2)
Figure FDA0003056440790000021
Figure FDA0003056440790000022
Figure FDA0003056440790000023
Figure FDA0003056440790000024
步骤S303：若被测管道处于塑性极限承载状态，则应满足以下两个条件：
a)：受拉区与受压区的形心连线平行于z轴，即yc1＝yc2；
b)：受拉区面积等于受压区面积，即A＝AI＝AΠ；
由以上两个条件，计算可得：
Figure FDA0003056440790000025
Figure FDA0003056440790000026
步骤S304：依据步骤S303的计算结果，进而得到被测管道在塑性极限承载状态下的极限弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2t sinβcosγ-2Rm 2h sinθcosα) (9)
当缺陷相对z轴对称时，即α＝γ＝0，其塑性极限承载弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2t sinβ-2Rm 2h sinθ) (10)
步骤S305：将被测管道的腐蚀缺陷位置及腐蚀尺寸，代入到公式(9)或公式(10)，即可得到被测管道的塑性极限弯矩。
2.根据权利要求1所述的含腐蚀缺陷的海底管道结构弯曲失效评估方法，其特征在于，所述步骤S1中，所述腐蚀缺陷的位置包括腐蚀坑位置角。
3.根据权利要求1或2所述的含腐蚀缺陷的海底管道结构弯曲失效评估方法，其特征在于，所述步骤S1中，所述腐蚀尺寸包括腐蚀坑的腐蚀宽度、腐蚀深度和腐蚀长度。
4.根据权利要求1-3任一项所述的含腐蚀缺陷的海底管道结构弯曲失效评估方法，其特征在于，所述步骤S2中采用海浪谱进行所述的海浪载荷计算，并考虑海流载荷的联合作用。
5.一种海底管道塑性极限弯矩的计算方法，其特征在于，所述计算方法具体为：
步骤S1：依据被测管道的腐蚀缺陷的位置及腐蚀尺寸，确定塑性中性轴的位置，并按照公式(1)和公式(2)计算得被测管道的受拉区面积AI和受压区面积AΠ：
AΙ＝2Rmtβ (1)
AΙ＝2Rmt(π-β)-2Rhθ (2)
其中：
Figure FDA0003056440790000031
为平均半径；R为管道外壁半径；t为被测管道的厚度，t远小于R；h为腐蚀深度，θ为腐蚀宽度的半角；
步骤S2：依据被检测管道的腐蚀缺陷的位置及腐蚀尺寸，能够得出塑性中性轴的角度β、腐蚀缺陷的位置角α、塑性中性轴与水平轴的夹角γ，进而按照公式(3)-(6)得出受拉区区域和受压区区域的形心坐标(yc1，zc1)和(yc2，zc2)
Figure FDA0003056440790000032
Figure FDA0003056440790000033
Figure FDA0003056440790000034
Figure FDA0003056440790000035
步骤S3：若被测管道处于塑性极限承载状态，则应满足以下两个条件：
a)：受拉区与受压区的形心连线平行于z轴，即yc1＝yc2；
b)：受拉区面积等于受压区面积，即A＝AI＝AΠ；
由以上两个条件，计算可得：
Figure FDA0003056440790000036
Figure FDA0003056440790000037
步骤S4：依据步骤S303的计算结果，进而得到被测管道在塑性极限承载状态下的极限弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2t sinβcosγ-2Rm 2h sinθcosα) (9)
当缺陷相对z轴对称时，即α＝γ＝0，其塑性极限承载弯矩为：
Mp＝σyA(zc2-zc1)＝σy(4Rm 2t sinβ-2Rm 2h sinθ) (10)
步骤S5：将被测管道的腐蚀缺陷位置及腐蚀尺寸，代入到公式(9)或公式(10)，即可得到被测管道的塑性极限弯矩。
6.根据权利要求5的所述的海底管道塑性极限弯矩的计算方法，其特征在于，所述步骤S1中，所述腐蚀缺陷的位置包括腐蚀坑位置角。
7.根据权利要求5或6的所述的海底管道塑性极限弯矩的计算方法，其特征在于，所述步骤S1中，所述腐蚀尺寸包括腐蚀坑的腐蚀宽度、腐蚀深度和腐蚀长度。
CN201810917814.XA 2018-08-13 2018-08-13 一种含腐蚀缺陷的海底管道结构弯曲失效评估方法 Active CN110822294B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810917814.XA CN110822294B (zh) 2018-08-13 2018-08-13 一种含腐蚀缺陷的海底管道结构弯曲失效评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810917814.XA CN110822294B (zh) 2018-08-13 2018-08-13 一种含腐蚀缺陷的海底管道结构弯曲失效评估方法

Publications (2)

Publication Number Publication Date
CN110822294A CN110822294A (zh) 2020-02-21
CN110822294B true CN110822294B (zh) 2021-07-13

Family

ID=69547107

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810917814.XA Active CN110822294B (zh) 2018-08-13 2018-08-13 一种含腐蚀缺陷的海底管道结构弯曲失效评估方法

Country Status (1)

Country Link
CN (1) CN110822294B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111797553B (zh) * 2020-06-28 2022-09-30 中国石油大学(北京) 复杂荷载条件下腐蚀管道的极限内压预测方法及装置
CN112052609B (zh) * 2020-08-14 2022-09-16 大连理工大学 一种不规则缺陷管道内压作用下的失效压力计算方法
CN112052595B (zh) * 2020-09-08 2022-10-21 大连理工大学 一种腐蚀钢制管道外压临界弹塑性屈曲压力计算方法
CN112052594B (zh) * 2020-09-08 2022-10-21 大连理工大学 一种含双腐蚀缺陷钢制管道临界弹塑性屈曲压力计算方法
CN112052616B (zh) * 2020-09-08 2022-10-21 大连理工大学 一种含任意腐蚀缺陷钢制管道外压临界弹塑性屈曲压力计算方法
CN112525692B (zh) * 2020-11-30 2022-12-02 中国石油天然气集团有限公司 一种埋地含缺陷管道承受极限载荷评估方法
CN112697608B (zh) * 2020-12-10 2022-09-13 太原理工大学 一种钢构件双向压弯下全截面塑性抗弯承载力的判定方法
CN113742957B (zh) * 2021-07-22 2023-08-18 海洋石油工程股份有限公司 一种基于弹性屈曲失效模式的海底管道壁厚设计方法
CN113792447B (zh) * 2021-07-30 2023-08-22 海洋石油工程股份有限公司 一种基于点荷载作用的大口径海底管道局部屈曲设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6000277A (en) * 1998-06-05 1999-12-14 Southwest Research Institute Method for the prediction of rupture in corroded transmission pipes subjected to combined loads
CN101726456A (zh) * 2008-10-15 2010-06-09 中国石油天然气集团公司 一种含腐蚀缺陷注汽管道补偿器弯管剩余强度评价方法
CN102636367A (zh) * 2012-04-23 2012-08-15 浙江大学 一种模拟风力及海流载荷的多自由度动力加载装置
CN103995957A (zh) * 2014-04-18 2014-08-20 华东理工大学 针对蠕变损伤及体积型缺陷的承压结构的安全评定方法
CN105403344A (zh) * 2015-12-16 2016-03-16 浙江大学 管道实时应力的获取方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6000277A (en) * 1998-06-05 1999-12-14 Southwest Research Institute Method for the prediction of rupture in corroded transmission pipes subjected to combined loads
CN101726456A (zh) * 2008-10-15 2010-06-09 中国石油天然气集团公司 一种含腐蚀缺陷注汽管道补偿器弯管剩余强度评价方法
CN102636367A (zh) * 2012-04-23 2012-08-15 浙江大学 一种模拟风力及海流载荷的多自由度动力加载装置
CN103995957A (zh) * 2014-04-18 2014-08-20 华东理工大学 针对蠕变损伤及体积型缺陷的承压结构的安全评定方法
CN105403344A (zh) * 2015-12-16 2016-03-16 浙江大学 管道实时应力的获取方法

Also Published As

Publication number Publication date
CN110822294A (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN110822294B (zh) 一种含腐蚀缺陷的海底管道结构弯曲失效评估方法
Wang et al. Strain-based design-advances in prediction methods of tensile strain capacity
CN110765505B (zh) 一种含表面划伤复合凹陷油气管道极限内压的预测方法
Besel et al. Corrosion assessment method validation for high-grade line pipe
Chauhan et al. Developments in methods for assessing the remaining strength of corroded pipelines
de Souza et al. Burst tests on pipeline containing long real corrosion defects
Kashyzadeh et al. Assessment of oil storage tanks performance containing cracks and cavities
CN108647441B (zh) 一种联合载荷作用下受损管道剩余弯矩极限载荷计算方法
Zhu et al. Corrosion assessment models for predicting remaining strength of corroded thick-walled pipelines
Kaveh et al. Structural integrity of offshore pipelines considering buckling and fracture limit-states
Li et al. Research on remaining bearing capacity evaluation method for corroded pipelines with complex shaped defects
Lu et al. A finite element-based analysis approach for computing the remaining strength of the pressure equipment with a local thin area defect
Kihiu et al. Stress characterization of autofrettaged thick-walled cylinders
Ameh et al. Dent and gouge defects assessment: A case study of gas pipeline
Zhao et al. Petroleum Science
Sun et al. Study on Evaluation Method of Failure Pressure for Pipeline with Axially Adjacent Defects
Yamaguchi Investigation of burst pressure in pipes with square wall thinning by using FEA and API579 FFS-1
Mokhtari et al. Numerical analysis of pit-to-crack transition under corrosion fatigue using a stochastic pit generation algorithm
CN113591348B (zh) 一种火电厂在役汽水管道焊缝三维应力计算方法
Gandomkar et al. The Crack Resolution Required to Ensure Fatigue Life of a Rescue Submarine’s Transfer Skirt
WU et al. Crack Evaluation of an Ultra-High-Pressure Chamber Used for Deep-Sea Environment Simulation
Widiyanto et al. Numerical analysis of stiffened offshore pipe subjected to environmental loading: A study case using external pressure
Dent et al. New axially loaded full ring test method for assessment of susceptibility of girth welds and parent pipe to sour service cracking
Chen et al. Burst capacity solutions for submarine pipeline with long corrosion defects
Kemp et al. Using Finite Element Analysis to Prioritize ILI Calls for Combined Features: Dents in Bends

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Yang Dongping

Inventor after: Ma Wenfeng

Inventor after: Wang Qiang

Inventor after: Cong Jun

Inventor after: Guo Aihong

Inventor after: Wang Weibin

Inventor after: Niu Gengqi

Inventor after: Zhi Jingbo

Inventor after: Zhu Liguo

Inventor after: Gao Shasha

Inventor before: Yang Dongping

Inventor before: Ma Wenfeng

Inventor before: Wang Qiang

Inventor before: Guo Aihong

Inventor before: Cong Jun

Inventor before: Wang Weibin

Inventor before: Niu Gengqi

Inventor before: Zhi Jingbo

Inventor before: Zhu Liguo

Inventor before: Gao Shasha

GR01 Patent grant
GR01 Patent grant