CN114004171B - 一种确定管道止裂韧性裂纹尖端张开角的方法 - Google Patents

一种确定管道止裂韧性裂纹尖端张开角的方法 Download PDF

Info

Publication number
CN114004171B
CN114004171B CN202111226026.4A CN202111226026A CN114004171B CN 114004171 B CN114004171 B CN 114004171B CN 202111226026 A CN202111226026 A CN 202111226026A CN 114004171 B CN114004171 B CN 114004171B
Authority
CN
China
Prior art keywords
crack
pipeline
pressure
natural gas
velocity curve
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
CN202111226026.4A
Other languages
English (en)
Other versions
CN114004171A (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.)
Centre Of Excellence For Advanced Materials
Original Assignee
China University of Petroleum East China
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 University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202111226026.4A priority Critical patent/CN114004171B/zh
Publication of CN114004171A publication Critical patent/CN114004171A/zh
Application granted granted Critical
Publication of CN114004171B publication Critical patent/CN114004171B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Human Resources & Organizations (AREA)
  • Algebra (AREA)
  • General Business, Economics & Management (AREA)
  • Strategic Management (AREA)
  • Primary Health Care (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Tourism & Hospitality (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

本发明涉及一种确定管道止裂韧性裂纹尖端张开角(CTOA)的方法，包括：将利用BWRS气体状态方程计算气体压力‑减压波速度曲线及利用基于CTOA的裂纹扩展速度预测模型构建天然气管道的压力‑裂纹扩展速度曲线放置于同一坐标系中，通过改变CTOA值，使气体压力‑减压波速度曲线与压力‑裂纹扩展速度曲线相切，将气体压力‑减压波速度曲线与压力‑裂纹扩展速度曲线的切点CTOA作为天然气管道的止裂韧性。相对于现有技术，本发明以CTOA作为止裂韧性指标，对于管道抵抗断裂的能力表征的更准确，对于不同钢级和管道尺寸适用性更高。

Description

一种确定管道止裂韧性裂纹尖端张开角的方法
技术领域
本发明属于天然气输送管道设计制造技术领域，具体地，涉及一种利用基于裂纹尖端张开角(CTOA)的双曲线模型确定管道止裂韧性CTOA的方法。
背景技术
随着国民经济的不断发展，天然气的需求随之日益增长，管道作为天然气输送最有效的手段，正朝着高钢级、大口径、高压力的方向发展。高钢级管道尽管强度高韧性好，但仍不可避免的在生产或运输过程中产生缺陷或裂纹，这些裂纹在管内高压作用下迅速扩展至成百甚至上千米，造成灾难-性的后果。因此，必须对管道长程裂纹扩展进行有力的控制，而止裂韧性的确定则是断裂控制的核心。
Battelle双曲线模型(BTCM)是目前工业化程度最高、应用最广泛的止裂韧性预测模型，该模型以夏比冲击韧性为指标。工程实践表明，该模型对于X50、X65等早期管线钢止裂韧性的预测是较为准确的，对于X70及以上高韧性管线钢却并不适用。这是由于夏比冲击试样较小，其预测得到的夏比冲击能量主要代表对塑性变形而非对裂纹扩展的抵抗，因此，夏比冲击韧性是否可以代表管道的真实韧性值得商榷。
由于落锤撕裂试验试样尺寸与试验工况更接近于管道裂纹扩展的真实状态，因此被认为是替代夏比冲击试验的最好方法。基于上述认识，日本高强度管线钢管(HLP)委员会以落锤撕裂能代替BTCM中的夏比冲击韧性，通过对全尺寸X70管道爆破实验数据的拟合构建了HLP模型；之后，日本Sumitomo公司对HLP模型引入管道尺寸效应修正参数，形成了HLP-Sumitomo模型，但两模型并未得到普遍认可。随着管道钢级的提升，基于夏比冲击韧性或基于落锤撕裂能的止裂控制方法预测得到的止裂韧性变得越来越不确定。这给高钢级天然气管道止裂设计造成巨大影响，也给后续管道安全运行造成巨大威胁。
多种试验数据表明，材料CTOA仅与裂纹稳定扩展阶段有关，有望克服夏比冲击韧性及落锤撕裂能的根本缺陷，成为表征天然气管道延性断裂止裂韧性的首选参数，然而当前尚未形成确定止裂韧性CTOA的相关方法。
发明内容
为克服现有技术的缺陷，本发明提供一种利用基于CTOA的双曲线模型确定管道止裂韧性CTOA的方法。
为实现上述目的，本发明的技术方案如下：
一种确定管道止裂韧性CTOA的方法，包括以下步骤：
1)根据天然气管道设计要求，确定天然气管道的规格参数和最大设计输送压力P0；所述天然气管道的规格参数包括：天然气管道外径和天然气管道壁厚；
2)采用BWRS气体状态方程计算气体压力-减压波速度曲线，输入的参数包括温度、压力、各天然气组分摩尔百分比，输出为气体压力-减压波速度曲线；
3)利用基于CTOA的裂纹扩展速度预测模型构建天然气管道的裂纹扩展阻力方程，绘制天然气管道的压力-裂纹扩展速度曲线；
4)将步骤2)获取的气体压力-减压波速度曲线及步骤3)获取的天然气管道的压力-裂纹扩展速度曲线放置于同一坐标系中，然后不断改变CTOA值，使天然气管道的气体压力-减压波速度曲线与天然气管道的压力-裂纹扩展速度曲线相切，两曲线刚好相切时的CTOA则为天然气管道在当前工况下实现自身止裂需要的止裂韧性。
步骤2)所述的BWRS状态方程计算天然气管道的气体压力-减压波速曲线，按照如下公式计算：
Figure BDA0003314355650000031
其中，P为系统压力，kPa；ρ为气相或液相的密度，kmol/m3；R为通用气体常数，J/(mol·k)；T为系统温度，K；A0，B0，C0，D0，E0，a，b，c，d，α,γ为方程特征参数。
步骤3)所述的天然气管道裂纹扩展阻力方程，按照如下公式计算：
Figure BDA0003314355650000032
Figure BDA0003314355650000033
其中，Pa为止裂压力，MPa；σf为流变应力，MPa；t为管道壁厚，mm；MT为鼓胀系数；E为弹性模量，MPa；σy为屈服应力，MPa；εl为颈缩开始时的有效塑性应变；K为材料常数；n是应变速率敏感性指数；r是管道半径，mm；Vf为裂纹扩展速度。
本发明具有以下有益效果：
(1)本发明所述的利用基于CTOA的双曲线模型确定管道止裂韧性CTOA的方法在具体操作时，通过改变管道的CTOA，使天然气管道的气体压力-减压波速曲线与压力-裂纹扩展速度曲线相切进而确定天然气管道自身止裂的临界止裂韧性，原理简单，计算精度高。
(2)该方法用于计算压力-裂纹扩展速度曲线的公式(1)和(2)由理论手段结合数值模拟得到，推导过程中未对材料及管道尺寸做任何假设，因此本发明对于任意钢级和管道尺寸均普遍适用，克服了传统方法适用性差的问题。
(3)该方法所确定的止裂韧性指标为管道CTOA，而非传统夏比冲击韧性、落锤撕裂能等参数，对于高钢级管道止裂韧性表征的更准确。一方面提高了管道止裂设计的可靠性，降低了长程裂纹扩展事故发生的可能性，另一方面可以显著节约管道生产及建设成本，应用价值高。
附图说明
图1为本发明确定管道止裂韧性CTOA的示意图；
图2为本发明中止裂韧性CTOA的计算结果图。
具体实施方式
下面结合附图对本发明做进一步详细描述：
如图1所示，本发明提供一种管道止裂韧性CTOA的确定方法，其特征在于，包括以下步骤：
1)通过BWRS状态方程获取天然气管道的气体压力—减压波速曲线，输入的参数包括：温度、压力、各天然气组分摩尔百分比，输出为气体压力-减压波曲线。其中，BWRS状态方程为：
Figure BDA0003314355650000041
其中，P为系统压力，kPa；ρ为气相或液相的密度，kmol/m3；R为通用气体常数，J/(mol·k)；T为系统温度，K；A0，B0，C0，D0，E0，a，b，c，d，α,γ为方程特征参数。
2)天然气管道压力—裂纹扩展速度曲线采用如下所示的公式1和2进行计算：
Figure BDA0003314355650000051
Figure BDA0003314355650000052
其中，Pa为止裂压力，MPa；σf为流变应力，MPa；t为管道壁厚，mm；MT为鼓胀系数；E为弹性模量，MPa；σy为屈服应力，MPa；εl为颈缩开始时的有效塑性应变；K为材料常数；n是应变速率敏感性指数；r是管道半径，mm；Vf为裂纹扩展速度。
3)将步骤1)获取的天然气管道的气体压力-减压波速曲线及步骤2)获取的天然气管道的压力—裂纹扩展速度曲线放置于同一坐标系中，然后不断改变CTOA的数值，使天然气管道的气体压力-减压波速曲线与天然气管道的压力—裂纹扩展速度曲线相切，天然气管道的气体压力-减压波速曲线与天然气管道的压力—裂纹扩展速度曲线相切时的CTOA即为天然气管道自身质量的临界CTOA。
下面以具体实例进行说明：
(1)输送气体组分如表1所示，天然气温度为12℃,通过BWRS状态方程计算天然气管道的气体压力—减压波速曲线。
表1天然气组分
天然气组分 甲烷 乙烷 丙烷
摩尔百分比 92％ 6％ 2％
(2)确定天然气管道的规格参数和输送压力P0如表2所示，根据公式1和2计算天然气管道的压力—裂纹扩展速度曲线。
表2天然气管道参数
Figure BDA0003314355650000053
Figure BDA0003314355650000061
(3)根据气体压力-减压波速曲线与压力—裂纹扩展速度曲线得到管道止裂韧性CTOA＝8.2°，如图2所示。由此可知，当管道管材测试所得CTOA大于等于8.2°时，管道即使出现裂纹，也会因自身材料性能而止裂，不会发生长程扩展。

Claims (3)

1.一种确定管道止裂韧性裂纹尖端张开角的方法，其特征在于，包括以下步骤：
1)根据天然气管道设计要求，确定天然气管道的规格参数和最大设计输送压力P0；
2)通过BWRS气体状态方程计算气体压力-减压波速度曲线；
3)利用基于CTOA的裂纹扩展速度预测模型构建天然气管道的裂纹扩展阻力方程，绘制天然气管道的压力-裂纹扩展速度曲线；
4)将步骤1)将步骤2)获取的气体压力-减压波速度曲线及步骤3)获取的天然气管道的压力-裂纹扩展速度曲线放置于同一坐标系中，然后不断改变CTOA值，使天然气管道的气体压力-减压波速度曲线与天然气管道的压力-裂纹扩展速度曲线相切，两曲线刚好相切时的CTOA则为天然气管道在当前工况下实现自身止裂需要的止裂韧性。
2.权利要求1所述的确定管道止裂韧性裂纹尖端张开角的方法，其特征在于，BWRS气体状态方程为：
Figure FDA0003314355640000011
其中，P为系统压力，kPa；ρ为气相或液相的密度，kmol/m3；R为通用气体常数，J/(mol·k)；T为系统温度，K；A0，B0，C0，D0，E0，a，b，c，d，α,γ为方程特征参数。
3.权利要求1所述的确定管道止裂韧性裂纹尖端张开角的方法，其特征在于，天然气管道的裂纹阻力方程为：
Figure FDA0003314355640000021
Figure FDA0003314355640000022
其中，Pa为止裂压力，MPa；σf为流变应力，MPa；t为管道壁厚，mm；MT为鼓胀系数；E为弹性模量，MPa；σy为屈服应力，MPa；εl为颈缩开始时的有效塑性应变；K为材料常数；n是应变速率敏感性指数；r是管道半径，mm；Vf为裂纹扩展速度。
CN202111226026.4A 2021-10-21 2021-10-21 一种确定管道止裂韧性裂纹尖端张开角的方法 Active CN114004171B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111226026.4A CN114004171B (zh) 2021-10-21 2021-10-21 一种确定管道止裂韧性裂纹尖端张开角的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111226026.4A CN114004171B (zh) 2021-10-21 2021-10-21 一种确定管道止裂韧性裂纹尖端张开角的方法

Publications (2)

Publication Number Publication Date
CN114004171A CN114004171A (zh) 2022-02-01
CN114004171B true CN114004171B (zh) 2022-08-23

Family

ID=79923366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111226026.4A Active CN114004171B (zh) 2021-10-21 2021-10-21 一种确定管道止裂韧性裂纹尖端张开角的方法

Country Status (1)

Country Link
CN (1) CN114004171B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462696A (zh) * 2014-12-11 2015-03-25 中国石油天然气集团公司 一种天然气管道止裂预测的方法
CN109187202A (zh) * 2018-07-24 2019-01-11 中国石油天然气集团有限公司 利用天然气管道全尺寸爆破试验计算管道止裂韧性的方法
CN113190789A (zh) * 2021-04-29 2021-07-30 中国石油天然气集团有限公司 一种输气管道裂纹扩展速度计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462696A (zh) * 2014-12-11 2015-03-25 中国石油天然气集团公司 一种天然气管道止裂预测的方法
CN109187202A (zh) * 2018-07-24 2019-01-11 中国石油天然气集团有限公司 利用天然气管道全尺寸爆破试验计算管道止裂韧性的方法
CN113190789A (zh) * 2021-04-29 2021-07-30 中国石油天然气集团有限公司 一种输气管道裂纹扩展速度计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高钢级输气管道裂纹延性扩展及止裂模拟技术研究;贺云婷;《焊管》;20181028(第10期);全文 *

Also Published As

Publication number Publication date
CN114004171A (zh) 2022-02-01

Similar Documents

Publication Publication Date Title
CN105117536A (zh) 一种rpv含裂纹类缺陷的简化弹塑性断裂力学分析方法
CN112948909B (zh) 各向同性双管混凝土柱承载力计算方法及系统
CN105043905B (zh) 一种基于稳态流变速率切线的岩石长期强度参数确定方法
CN108710721A (zh) 一种基于失效评定图的未爆先漏分析方法
CN103926136B (zh) 管线钢平面应变断裂韧性和安全临界壁厚的确定方法
CN104992004B (zh) 一种管道止裂预测用临界单位面积损伤应变能的确定方法
CN108562495B (zh) 一种完好钢制管道极限内压的计算方法
CN114004171B (zh) 一种确定管道止裂韧性裂纹尖端张开角的方法
Wang et al. Confined buckle propagation of offshore pipelines under external pressure
CN113190789B (zh) 一种输气管道裂纹扩展速度计算方法
CN103164624B (zh) 获得平行隧道轴线的匀质地下管线状态数据的方法
CN110110386B (zh) 一种用于高钢级输气管道的钢套筒止裂器设计方法
CN114896879B (zh) 一种粘弹性管道直接水锤压力的计算方法
Amaro et al. CTOA testing of pipe line steels using MDCB specimens
Amano et al. Evaluation of leak/rupture behavior for axially part-through-wall notched high-strength line pipes
CN110516269B (zh) 一种用于海底管道工程临界评估的优化参考应变方法
Jia et al. A modified engineering critical assessment method for deeply-embedded cracks in metallic pipelines subjected to large plastic strain
Selker et al. Elastic Anisotropy in Linepipe Materials
CN110566755B (zh) 一种用于x100输气管道玻璃纤维复合材料止裂器的设计方法
Hashemi et al. A single specimen CTOA test method for evaluating the crack tip opening angle in gas pipeline steels
CN106844964A (zh) 一种压力管道输水系统非恒定流模型的优化方法
CN106769531A (zh) 一种低硬度p91管件持久曲线外推函数的建立方法
CN108645715B (zh) 由无缝试件确定混凝土岩石断裂韧度的方法
Hosoe et al. Underwater burst tests for evaluating unstable ductile crack arrestability in offshore pipelines
Deng et al. Study on the burst pressure of support cylinder in expandable liner hanger based on matching between yield strength and fracture toughness

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20240801

Address after: 523808 room 902-905, building 2, No.4 Qinyuan Road, Songshanhu Park, Dongguan City, Guangdong Province

Patentee after: CENTRE OF EXCELLENCE FOR ADVANCED MATERIALS

Country or region after: China

Address before: 266580 No. 66 Changjiang West Road, Huangdao District, Qingdao, Shandong.

Patentee before: CHINA University OF PETROLEUM (EAST CHINA)

Country or region before: China

TR01 Transfer of patent right