CN109255139B - 一种高温管道中表面多裂纹合并方法 - Google Patents
一种高温管道中表面多裂纹合并方法 Download PDFInfo
- Publication number
- CN109255139B CN109255139B CN201710574812.0A CN201710574812A CN109255139B CN 109255139 B CN109255139 B CN 109255139B CN 201710574812 A CN201710574812 A CN 201710574812A CN 109255139 B CN109255139 B CN 109255139B
- Authority
- CN
- China
- Prior art keywords
- crack
- cracks
- depth
- flat plate
- length
- 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
Links
Images
Classifications
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种高温管道中含表面多裂纹干涉程度分析和合并方法,本发明提出高温下多裂纹干涉程度取决于裂纹深度和相邻裂纹之间的距离,采用有限元计算的方法获得干涉程度计算公式的参数。在实际管道评估中根据服役温度下多裂纹干涉程度计算公式,评估多裂纹干涉程度,并根据干涉程度进行多裂纹复合。本发明的有益效果是:提出了高温管道含表面裂纹的多裂纹干涉程度计算公式、裂纹合并方法,本方法完善了现有评估标准的不足,可以提高预测的精度。
Description
技术领域
本发明涉及含裂纹高温管道寿命评估技术领域,具体涉及含多裂纹高温管道中多裂纹干涉程度分析和复合方法,提高剩余寿命预测和可靠性评估的精度,用于高温管道运行维护和管理。
背景技术
高温环境下服役的蒸汽管道,在生产和服役的过程中难免会产生各种缺陷或裂纹。蠕变裂纹的萌生和扩展是含缺陷高温构件的一个主要失效机制,并导致其在设计寿命之前失效。而对于含裂纹部件不能简单判断为不安全和不可靠。通常基于合乎使用原则,按照断裂力学的理论,进行严格的理论分析对存在的裂纹做出科学的评定,以判断裂纹是否危害部件安全可靠性,以及对裂纹的发展、后果等做出判断,并判断含裂纹的高温承压设备在规定的使用工况下可不可以继续安全使用。若裂纹的存在不会对安全生产构成威胁,则可以被允许存在;若裂纹的存在暂时不会对安全生产构成威胁,但进一步发展可能会有潜在威胁,则可以被允许在监控下使用;若裂纹降低使用要求后不会构成威胁,则考虑降低要求使用;已经对安全生产构成威胁的裂纹,必须采取措施,返修或者停止使用。合乎使用原则不仅保障了高温焊接结构的安全性,还考虑到经济性,能够获得非常可观的经济效益。目前通用的含裂纹高温部件评定程序和标准,如美国的ASME/API RP 579,法国的RCC-MR(Appendix 16)、英国的R5和BS 7910、我国的GB/T 19624-2004等标准和规范都是基于此原则。
高温部件中发现的裂纹通常都不是独立存在的,往往是两个或多个裂纹(Multiple cracks)共存。在实际工程结构中,多个共面裂纹最较为常见;而非共面裂纹也通常规则化为共面裂纹进行考虑。但现行含裂纹高温结构的评估标准和方法,对多个缺陷的处理方法比较保守。一方面现行方法大多通过弹塑性分析建立的;另一方面,只是根据多裂纹间的几何关系简单判定是否复合的判定,即当多裂纹的距离关系满足一定条件时进行复合,没有考虑裂纹深度等影响,会导致不必要的维修和更换。
发明内容
本发明目的旨在改进现有含裂纹高温管道寿命预测和评估技术中表面多裂纹处理方法,降低预测的保守性,有助于高效、安全、合理的实现含缺陷高温管道的可靠性评估。
本发明所提出的高温管道表面多裂纹的干涉效应和复合方法包括以下步骤:
(1)、高温管道表面多裂纹干涉效应计算方法,适用温度范围是超过材料0.4Tm且小于Tm,其中Tm是材料的熔点,单位K,选取公式模型为:
其中,Mcreep代表蠕变下多裂纹的干涉程度,决定高温管道中多裂纹是否需要复合;a合并前为裂纹原始深度,即相邻两个裂纹深度值较大的一个,c合并前为裂纹原始半长,即邻两个裂纹半长值较大的一个,t是管道厚度,s是两个相邻裂纹的距离,其中a合并前、c合并前、t、s的单位均为mm。
式中,A1,A2,A3,A4为4个待定无量纲参数。
参数的确定方法来源于数据拟合得出的经验公式,数据由有限元模拟得到,包括以下步骤:
S1:建立平板有限元模型,在平板中插入预制表面椭圆形裂纹,如图1所示,裂纹半宽为c;裂纹最深处的深度为a,裂纹深度a=0.6c;平板长度L为20c,宽度W为10c,平板厚度t=5a,裂纹尖端最小单元尺寸选为a/50mm。在平板不含裂纹的自由端施加任意恒定应力,计算裂纹前沿表面点(两个裂纹前沿表面点的位置如图1所示)处的断裂参量C*分布,记为(当材料处于稳态蠕变情况时,蠕变裂尖断裂参量可用C*积分表示,C*的计算公式如下:
参见文献:(Yun-Jae Kim,Jin-Su Kim,Nam-Su Huh,Young-Jin Kim.EngineeringC-integral estimates for generalised creep behavior and finite elementvalidation.International Journal of Pressure Vessels and Piping,2002,79:427-443).
S2:建立平板有限元模型,在平板中插入两个等大预制表面椭圆裂纹,如图2所示。两个裂纹尺寸相同,裂纹半宽均为c;裂纹最深处的深度均为a,裂纹深度a=0.6c;平板长度L为20c,宽度W为10c,平板厚度t=5a,裂纹尖端最小单元尺寸选为a/50mm。两个裂纹的距离s分别选为0.5c,1c,2c和3c。在平板不含裂纹的自由端施加任意恒定应力,计算不同距离下两个裂纹靠近位置即表面点的断裂参量C*分布,记为两个裂纹靠近位置表面点如图2所示。
S5:采用非线性数据拟合的方法(最小二乘数法),获得蠕变多裂纹干涉程度计算公式的4个参量A1,A2,A3,A4。
(2)、高温管道表面缺陷的规则化,高温管道材料与步骤1中有限元模型所用材料形同。
依据高温管道无损检验的结果,确定缺陷或裂纹的尺寸和准确位置。考虑无损检验的精度,对探测尺寸到裂纹尺寸乘上一个安全系数。安全系数根据现行标准中对无损探伤要求进行。根据检测缺陷的长度和深度,把表面缺陷规则化为椭圆形表面缺陷,规则化时取缺陷最深位置的深度为规则化裂纹的深度,取裂纹最长处的长度为规则化裂纹的长度,a代表规则化裂纹的深度,c代表规则化裂纹的半长。如图3所示,对每个检测裂纹都按照上述步骤处理,分别标记为ai,ci,其中i代表裂纹的个数和次序,为1,2,3,4……。测量管道厚度t。
3、高温管道多裂纹复合程序。
S6:对于检测到多个裂纹,相邻的裂纹分别进行多裂纹干涉程度分析,分析时均采用裂纹原始尺寸。
a合并前=max(a1,a2)
c合并前=max(c1,c2)
其中a合并前为裂纹原始深度,c合并前为裂纹原始半长,a1,a2为相邻两个裂纹的半长,c1,c2为相邻两个裂纹的半宽。
S7:若Mcreep≥20.5,按照下面方法计算合并后裂纹尺寸:
a合并后=max(a1,a2)
2c合并后=2c1+2c2+s1
其中a合并后相邻两个裂纹合并后的深度,c合并后为相邻两个裂纹合并后的半宽,s1为相邻两个裂纹的间距。
S8:若存在两条以上裂纹,两两相邻的裂纹重复S6和S7步骤进行多裂纹干涉程度分析和裂纹复合。
优选的,所述适用温度范围是超过材料0.4Tm且小于0.8Tm。
优选的,所述安全系数为1.5。
优选的,所述无损检验通过超声探伤或磁粉探伤完成。
与现有技术相比,本发明的有益效果是:
本发明公开了一种高温管道中含表面多裂纹干涉程度分析和合并方法,本发明提出高温下多裂纹干涉程度取决于裂纹深度和相邻裂纹之间的距离,采用有限元计算的方法获得干涉程度计算公式的参数。在实际管道评估中根据服役温度下多裂纹干涉程度计算公式,评估多裂纹干涉程度,并根据干涉程度进行多裂纹复合。本发明的有益效果是:提出了高温管道含表面裂纹的多裂纹干涉程度计算公式、裂纹合并方法,本方法完善了现有评估标准的不足,可以提高预测的精度。
附图说明
图1所示为本发明的含单裂纹平板的示意图。
图2所示为本发明的含等大表面双裂纹平板的示意图。
图3所示为本发明中表面缺陷规则化为椭圆形裂纹示意图。
图4所示为本发明中实施案例含多裂纹管道示意图。
具体实施方式
下面结合具体实例进一步说明本发明的技术方案。
本发明的含多裂纹高温管道中多裂纹复合及重定义方法:
本实施例选取高等级耐热钢P92钢,使用温度为923K,管道规格是外径388mm,内径288mm,厚度t′为50mm,经无损检测发现三个表面裂纹,考虑安全系数后,尺寸分别为a1′=25mm,c1′=40mm,s1′=10mm,a2′=20mm,c1′=20mm,s2′=50mm,a3′=15mm,c3′=30mm,如图4所示。
(1)材料属性
弹性模量E(MPa) | 泊松比 | 蠕变指数A(MPa<sup>-n</sup>/h) | 蠕变硬化指数n |
125000 | 0.3 | 3.5085E-31 | 12.72 |
(2)P92钢多裂纹干涉效应公式计算。
S1:建立平板有限元模型,在平板中插入预制表面椭圆形裂纹,裂纹半宽为c′=10mm;裂纹深度为a′=0.6c′=6mm;平板长度L′为20c′=200mm,宽度W′为10c′=100mm,平板厚度t′=5a′=30mm,裂纹尖端最小单元尺寸选为a′/50mm=0.012mm。在平板不含裂纹的自由端施加应力70MPa,计算裂纹前沿的分布,如表3所示。
S2:建立平板有限元模型,在平板中插入两个等大预制表面椭圆裂纹。裂纹半宽为c′=10mm;裂纹深度为a′=0.6c′=6mm,;平板长度L′为20c′=200mm,宽度W′为10c′=100mm,平板厚度t′=5a=30mm,裂纹尖端最小单元尺寸选为a′/50mm=0.012mm。两个裂纹的距离s′分别选为0.5c′=5mm,1c′=10mm,2c′=20mm和3c′=30mm。在平板不含裂纹的自由端施加应力70MPa,计算不同距离下两个裂纹靠近位置分布,如表3所示。
表3多裂纹干涉因子计算条件及结果
S5:采用非线性拟合方式,依据计算结果,获得P92钢蠕变下多裂纹干涉效应计算公式如下:
(3)P92钢高温管道多裂纹复合
尺寸分别为a1′=25mm,c1′=40mm,s1′=10mm,a2′=20mm,c2′=20mm,s2′=60mm,a3′=15mm,c3′=30mm,如图所示
由于检测到三个裂纹,按顺序两两裂纹依次进行多裂纹干涉程度分析,
a合并前′=max(a1′,a2′)=25mm
c合并前′=max(c1′,c2′)=40mm
s1′/c合并前′=0.25
根据公式计算Mcreep′=2.07>20.5,则两个裂纹需要合并,合并后裂纹尺寸为:
a合并后1′=max(a1′,a2′)=25mm
2c合并后1′=2c1′+2c2′+s1′=80+40+10=130mm
c合并后1′=65mm
第二个和第三个裂纹的多裂纹干涉效应分析
a合并前′=max(a2′,a3′)=20mm
c合并前′=max(c2′,c3′)=30mm
s2′/c合并前′=50/30=1.67
根据公式计算Mcreep′=1.29<20.5,则两个裂纹不需要合并
综上所检测三条裂纹,根据本方法分析后,按照两条裂纹分别进行后续的寿命和可靠性评估,两条裂纹尺寸为a1新′=25mm,c1新′=65mm,a2新′=15mm,c2新′=30mm。
以上所述仅是本发明的优选实施方式,应当指出的是,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.一种高温管道中表面多裂纹的复合方法,其特征在于:包括以下步骤:
(1)、计算高温管道表面多裂纹干涉效应,适用温度范围为超过材料0.4Tm且小于Tm,其中Tm是材料的熔点,单位K,选取公式模型为:
其中,Mcreep代表蠕变下多裂纹的干涉程度,决定高温管道中多裂纹是否需要复合;a合并前为裂纹原始深度,即相邻两个裂纹深度值较大的一个,c合并前为裂纹原始半长,即邻两个裂纹半长值较大的一个,δ是管道厚度,s是两个相邻裂纹的距离,其中a合并前、c合并前、δ、s的单位均为mm,A1,A2,A3,A4为4个待定无量纲参数;
参数的确定方法来源于数据拟合得出的经验公式,数据由有限元模拟得到,包括以下步骤:
S1:建立平板有限元模型,在平板中插入预制表面椭圆形裂纹,裂纹半长为c;裂纹深度为a,裂纹深度a=0.6c;平板长度L为20c,宽度W为10c,平板厚度t=5a,裂纹尖端最小单元尺寸选为a/50mm,在平板不含裂纹的自由端施加任意恒定应力,计算裂纹前沿表面点处的断裂参量C*分布,记为
S2:建立平板有限元模型,在平板中插入两个等大预制表面椭圆裂纹,两个裂纹尺寸相同,裂纹半长均为c;裂纹深度均为a,裂纹深度a=0.6c;平板长度L为20c,宽度W为10c,平板厚度t=5a,裂纹尖端最小单元尺寸选为a/50mm,两个裂纹的距离s分别选为0.5c,1c,2c和3c,在平板不含裂纹的自由端施加任意恒定应力,计算不同距离下两个裂纹靠近位置即表面点的断裂参量C*分布,记为
S5:采用非线性数据拟合的方法,获得蠕变多裂纹干涉程度计算公式的4个参量A1,A2,A3,A4;
(2)、高温管道表面缺陷的规则化,其中高温管道材料与步骤(1)中有限元模型所用材料相同;
依据高温管道无损检验的结果,确定缺陷或裂纹的尺寸和准确位置,考虑无损检验的精度,对探测到裂纹尺寸乘上一个安全系数以将裂纹尺寸放大,安全系数根据现行标准中对无损探伤要求进行,根据检测缺陷的长度和深度,把表面缺陷规则化为椭圆形表面缺陷,规则化时取缺陷最深位置的深度为规则化裂纹的深度,取裂纹最长处的长度为规则化裂纹的长度,ai代表规则化裂纹的深度,ci代表规则化裂纹的半长其中i代表裂纹的个数和次序,为1,2,3,4……,测量管道厚度t;
(3)、高温管道多裂纹复合程序:
S6:对于检测到多个裂纹,相邻的裂纹分别进行多裂纹干涉程度分析,分析时均采用裂纹原始尺寸:
a合并前=max(a1,a2)
c合并前=max(c1,c2)
其中a合并前为裂纹原始深度,c合并前为裂纹原始半长,a1,a2为相邻两个裂纹的深度,c1,c2为相邻两个裂纹的半长;
S7:若Mcreep≥20.5,按照下面方法计算合并后裂纹尺寸:
a合并后=max(a1,a2)
2c合并后=2c1+2c2+s1
其中a合并后相邻两个裂纹合并后的深度,c合并后为相邻两个裂纹合并后的半长,s1为相邻两个裂纹的间距;
S8:若存在两条以上裂纹,两两相邻的裂纹重复S6和S7步骤进行多裂纹干涉程度分析和裂纹复合。
2.根据权利要求1所述的一种高温管道中表面多裂纹的复合方法,其特征在于:所述步骤(1)中适用温度范围是超过材料0.4Tm且小于0.8Tm。
3.根据权利要求1所述的一种高温管道中表面多裂纹的复合方法,其特征在于:所述安全系数为1.5。
4.根据权利要求1所述的一种高温管道中表面多裂纹的复合方法,其特征在于:所述无损检验通过超声探伤或磁粉探伤完成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710574812.0A CN109255139B (zh) | 2017-07-14 | 2017-07-14 | 一种高温管道中表面多裂纹合并方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710574812.0A CN109255139B (zh) | 2017-07-14 | 2017-07-14 | 一种高温管道中表面多裂纹合并方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109255139A CN109255139A (zh) | 2019-01-22 |
CN109255139B true CN109255139B (zh) | 2022-11-08 |
Family
ID=65051771
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710574812.0A Active CN109255139B (zh) | 2017-07-14 | 2017-07-14 | 一种高温管道中表面多裂纹合并方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109255139B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859601B (zh) * | 2019-04-26 | 2024-03-19 | 天津大学 | 一种含有环向埋藏裂纹v型坡口焊缝海底冶金复合管gui插件自动建模方法 |
CN111539142B (zh) * | 2020-04-20 | 2023-12-22 | 中车青岛四方机车车辆股份有限公司 | 一种管道裂纹扩展驱动力计算方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101162234A (zh) * | 2007-11-23 | 2008-04-16 | 华东理工大学 | 汽轮机高温部件的剩余寿命预测方法 |
CN104268383A (zh) * | 2014-09-17 | 2015-01-07 | 合肥通用机械研究院 | 一种含裂纹类缺陷高温压力管道的安全评定方法 |
CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
-
2017
- 2017-07-14 CN CN201710574812.0A patent/CN109255139B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101162234A (zh) * | 2007-11-23 | 2008-04-16 | 华东理工大学 | 汽轮机高温部件的剩余寿命预测方法 |
CN104268383A (zh) * | 2014-09-17 | 2015-01-07 | 合肥通用机械研究院 | 一种含裂纹类缺陷高温压力管道的安全评定方法 |
CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
Non-Patent Citations (2)
Title |
---|
"Experimental investigation of creep crack growth behavior in nickel base superally by constant displacement loading method at elevated temperature";Rui Bao etc.;《Materials Science & Engineering A》;20160413;全文 * |
"核电焊接材料晶间腐蚀评价与双裂纹断裂合并行为研究";邱智伟;《中国优秀硕士学位论文全文数据库 基础科学》;20130615(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109255139A (zh) | 2019-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6350771B1 (ja) | 圧潰強度予測方法 | |
CN109255139B (zh) | 一种高温管道中表面多裂纹合并方法 | |
Muránsky et al. | Evaluation of a self-equilibrium cutting strategy for the contour method of residual stress measurement | |
Zou et al. | Fatigue inspection and maintenance optimization: A comparison of information value, life cycle cost and reliability based approaches | |
JP2009174859A (ja) | 機械部品の余寿命評価方法 | |
Forsyth | Structural health monitoring and probability of detection estimation | |
Gao et al. | Critical strain based ductile damage criterion and its application to mechanical damage in pipelines | |
JP6225087B2 (ja) | 材料評価のための前処理方法 | |
JP4616778B2 (ja) | 高強度鋼溶接部の寿命評価方法 | |
JP2003207489A (ja) | 金属材料の損傷評価方法及び装置 | |
Nahal et al. | Failure probability assessment for pipelines under the corrosion effect | |
JP6077512B2 (ja) | 渦電流検査法を用いた伝熱管腐食量予測方法 | |
Mahajan et al. | Numerical characterization of compact pipe specimen for stretch zone width assessment | |
Song et al. | Residual stresses in weld repairs and mitigation by design | |
CN103984845B (zh) | 一种输气管道完整性评价周期的计算方法及系统 | |
Salama et al. | Characterization of the uncertainties in the inspection results of ultrasonic intelligent pigs | |
EP3109615B1 (en) | Ductile fracture evaluation method and device | |
CN113688546A (zh) | 一种金属焊缝役致缺陷接收表的制定方法 | |
Iwamatsu et al. | Evaluation procedure of limit load for non-aligned multiple flaws | |
KR20110034986A (ko) | 증기발생기 전열관 외경축 균열에 대한 모터구동 회전 탐촉자 와전류 탐상검사방법 | |
Mirsayar | Three dimensional investigation of mode I stress intensity factor variations in crack front using finite element method | |
Yoshida et al. | Effect of thickness measurement procedure on stress analysis of pipes with local metal loss | |
RU2014109580A (ru) | Способ определения срока службы трубопровода | |
Li | Acoustic emission monitoring and critical failure identification of bridge cable damage | |
RU2603884C1 (ru) | Способ определения допускаемого внутреннего давления для участка линейной части магистрального трубопровода со стресс-коррозионной трещиной |
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 |