CN111539142A - 一种管道裂纹扩展驱动力计算方法及*** - Google Patents
一种管道裂纹扩展驱动力计算方法及*** Download PDFInfo
- Publication number
- CN111539142A CN111539142A CN202010313384.8A CN202010313384A CN111539142A CN 111539142 A CN111539142 A CN 111539142A CN 202010313384 A CN202010313384 A CN 202010313384A CN 111539142 A CN111539142 A CN 111539142A
- Authority
- CN
- China
- Prior art keywords
- calculating
- driving force
- dimensionless coefficient
- finite element
- crack
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000004364 calculation method Methods 0.000 claims abstract description 69
- 239000000463 material Substances 0.000 claims abstract description 24
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 230000001419 dependent effect Effects 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 8
- 238000013075 data extraction Methods 0.000 claims description 8
- 238000007405 data analysis Methods 0.000 claims description 4
- 238000005482 strain hardening Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 23
- 238000004891 communication Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- LPLLVINFLBSFRP-UHFFFAOYSA-N 2-methylamino-1-phenylpropan-1-one Chemical compound CNC(C)C(=O)C1=CC=CC=C1 LPLLVINFLBSFRP-UHFFFAOYSA-N 0.000 description 3
- 241000682719 Adina Species 0.000 description 3
- 241000132539 Cosmos Species 0.000 description 3
- 235000005956 Cosmos caudatus Nutrition 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000001983 electron spin resonance imaging Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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]
-
- 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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明实施例提供一种管道裂纹扩展驱动力计算方法及***,其中,管道裂纹扩展驱动力计算方法,包括:在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;边界条件包括对称约束和固定约束;对有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在有限元计算后提取数据;通过有限元模型模拟得到的数据,计算得到第一无量纲系数;第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;根据第一无量纲系数,计算得到裂纹扩展驱动力。简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
Description
技术领域
本发明涉及断裂力学领域,尤其涉及一种管道裂纹扩展驱动力计算方法。
背景技术
管道在服役过程中受复杂地质条件(如不连续多年冻土带、地震、滑坡和地面沉降)、安装施工和实际运营等因素的影响,会产生1-3%的名义轴向应变。如果管道无法承受安装和服役时环境所产生的应变,就会发生泄漏或破裂。这不仅会造成巨大的经济损失,还将造成环境污染、生态破坏等灾难性后果,严重影响社会平稳运行。为了保证管道的结构完整性,需要依据基于断裂力学的BS7910和DNVGL-RP-F108等国际标准,对其进行工程临界评估(ECA)。在对承受塑性变形管道进行工程临界评估时,针对特定的材料和构件形式,需要精确测定其裂纹扩展驱动力。
迄今为止,失效评估图(FAD)是应用最为广泛的含缺陷结构断裂评估方法,其已被BS7910,R6,FITNET和API579等缺陷评定标准所采纳。失效评估图法有参考应力法和参考应变法两种,考虑断裂和失稳两种失效模式。目前,国际上使用较为广泛的几种基于应力的裂纹扩展驱动力(J积分)计算方法为:美国电力研究所(Electric power researchinstitute,EPRI)法、参考应力法和改进的参考应力法等。
上述方法的精度值极大地依赖于极限载荷解的选取,但在计算特定材料和构件形式的极限载荷解时,存在结果不够准确(导致断裂评估结果不够保守或过于保守)、计算极限载荷解的中间变量物理意义和理论基础不明确、未考虑裂纹尺寸和材料性能对极限载荷解的影响等问题。
因此,如何提供一种管道裂纹扩展驱动力计算方法及***,提出一个理论基础和物理意义明确、能精准计算裂纹扩展驱动力的简化参考应力法,对承受塑性变形管道进行断裂评估,成为亟待解决的问题。
发明内容
针对现有技术中的缺陷,本发明实施例提供一种管道裂纹扩展驱动力计算方法及***。
第一方面,本发明实施例提供一种管道裂纹扩展驱动力计算方法,包括,
在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;
通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
可选的,所述管道裂纹扩展驱动力计算方法,
所述通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数,具体包括:
通过所述有限元模型模拟以获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,计算得到第一无量纲系数;
其中,c为裂纹半长度;D为管道外径;a为裂纹深度;t为管道壁厚。
可选的,所述管道裂纹扩展驱动力计算方法,
所述通过所述有限元模型模拟以获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,计算得到第一无量纲系数,具体包括:
获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据;
通过参数化分析,构建第一无量纲系数λ的表达式:
使用Levenberg-Marquardt法对所述第一无量纲系数λ的表达式进行多项式拟合,得到第一无量纲系数λ的表达式中各拟合系数值;
将所述各拟合系数值代入λ的表达式中,得到所述第一无量纲系数λ的计算公式;
根据所述第一无量纲系数的计算公式,计算得到第一无量纲系数。
可选的,所述管道裂纹扩展驱动力计算方法,
所述第一无量纲系数λ的表达式中各拟合系数值为:
β0=1.1002,β1=0.0684,β2=0.0082,β3=0.1865,β4=2.6559,
β5=0.0289,β6=1.6237,β7=0.0003,β8=0.4599,β9=0.5851,
β10=0.4302;
所述第一无量纲系数λ的计算公式为:
可选的,所述管道裂纹扩展驱动力计算方法,
所述根据所述第一无量纲系数,计算得到裂纹扩展驱动力,具体包括:
根据所述第一无量纲系数λ,计算参考应力σref具体公式为:
σref=σ/λ;
其中,σ是远端施加应力;
根据所述参考应力σref,计算得到对应于σref的参考应变εref,进一步计算得到裂纹扩展驱动力;
参考应变εref计算公式为:
裂纹扩展驱动力计算公式为:
其中,α是硬化常数,σ0为屈服强度,ε0为屈服应变,n为应***化指数;
可选的,所述管道裂纹扩展驱动力计算方法,
所述第一无量纲系数λ的计算公式适用于任何屈服强度σ0、10≤D/t≤30、0.1≤a/t≤0.5、0.05≤2c/πD≤0.20、10≤n≤25以及0.015≤εn≤0.04的情形;其中,εn为名义应变。
可选的,所述管道裂纹扩展驱动力计算方法,
所述有限元模型采用ABAQUS软件或ANSYS软件构建。
第二方面,本发明实施例提供一种管道裂纹扩展驱动力计算***,包括:
模型设置模块,用于在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
数据提取模块,与模型设置模块相连,用于对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;
参数计算模块,与数据提取模块相连,用于通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
数据分析模块,与参数计算模块相连,用于根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
第三方面,本发明实施例提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上所述管道裂纹扩展驱动力计算方法的各个步骤。
第四方面,本发明实施例提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如上所述管道裂纹扩展驱动力计算方法的各个步骤。
本发明实施例提供一种管道裂纹扩展驱动力计算方法及***,通过构建一个取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的管道裂纹扩展驱动力计算方法流程图;
图2为本发明实施例提供的含外表面椭圆形裂纹的管道有限元模型示意图;
图3为本发明实施例提供的管道横截面示意图;
图4为本发明实施例提供的管道裂纹扩展驱动力计算***结构示意图;
图5为本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例提供的管道裂纹扩展驱动力计算方法流程图,如图1所示,该方法包括:
步骤S1,在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
步骤S2,对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;
步骤S3,通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
步骤S4,根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
具体的,步骤S1,通过有限元分析软件材料属性模块在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并通过有限元分析软件载荷模块设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
图2为本发明实施例提供的含外表面椭圆形裂纹的管道有限元模型示意图,如图2所示,所述管道本体受到远端均布拉伸载荷,在管道的厚度方向开有外表面椭圆形裂纹,裂纹面垂直于管道轴向,即承受轴向应变。
步骤S2,对所述有限元模型划分网格,有限元分析软件交互模块设置裂纹和输出参量,进行有限元计算,并通过有限元分析软件处理模块在所述有限元计算后提取数据;
步骤S3,通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
步骤S4,根据所述第一无量纲系数,以及相关裂纹扩展驱动力计算所需参数,计算得到裂纹扩展驱动力。
在上述实施例的基础上,步骤S1,在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件之前,还包括通过有限元分析软件构建所述管道的有限元模型;需要说明的是,目前流行的有限元分析软件主要有NASTRAN、ADINA、ANSYS、ABAQUS、MARC、MAGSOFT、COSMOS等,在本发明实施例中,不对构建限元模型的有限元分析软件进行限定,具体使用的有限元分析软件可根据实际情况进行选择。
本发明实施例提供一种管道裂纹扩展驱动力计算方法,通过构建一个取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
基于上述实施例,可选的,在上述管道裂纹扩展驱动力计算方法中,所述通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数,具体包括:
通过所述有限元模型模拟以获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,计算得到第一无量纲系数;
其中,c为裂纹半长度;D为管道外径;a为裂纹深度;t为管道壁厚。
具体的,由于所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数,因此通过所述有限元模型模拟以获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,通过大量的数据数学拟合得到第一无量纲系数λ与裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n的经验公式,进一步计算得到第一无量纲系数。
图3为本发明实施例提供的管道横截面示意图,如图3所示,其中,c为裂纹半长度,单位为mm;D为管道外径,单位为mm;a为裂纹深度,单位为mm;t为管道壁厚,单位为mm;r为管道半径,单位为mm。
本发明实施例提供一种管道裂纹扩展驱动力计算方法,通过构建一个与裂纹长度比、裂纹深度比、管道径厚比和硬化指数有关,仅取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
基于上述实施例,可选的,在上述管道裂纹扩展驱动力计算方法中,所述通过所述有限元模型模拟以获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,计算得到第一无量纲系数,具体包括:
获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据;
通过参数化分析,构建第一无量纲系数λ的表达式:
使用Levenberg-Marquardt法对所述第一无量纲系数λ的表达式进行多项式拟合,得到第一无量纲系数λ的表达式中各拟合系数值;
将所述各拟合系数值代入λ的表达式中,得到所述第一无量纲系数λ的计算公式;
根据所述第一无量纲系数的计算公式,计算得到第一无量纲系数。
具体的,获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,如下表所示:
对上述数据进行参数化分析,通过参数化分析,构建第一无量纲系数λ的表达式:
所述第一无量纲系数λ的表达式为变量是裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n的函数表达式。
使用Levenberg-Marquardt法对所述第一无量纲系数λ的表达式进行多项式拟合,得到第一无量纲系数λ的表达式中各拟合系数值;
将所述各拟合系数值代入λ的表达式中,得到所述第一无量纲系数λ的计算公式;
计算第一无量纲系数时,只需得到管道的裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n,就可将上述数据代入所述第一无量纲系数的计算公式,计算得到第一无量纲系数。
本发明实施例提供一种管道裂纹扩展驱动力计算方法,通过构建一个与裂纹长度比、裂纹深度比、管道径厚比和硬化指数有关,仅取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
基于上述实施例,可选的,在上述管道裂纹扩展驱动力计算方法中,所述第一无量纲系数λ的表达式中各拟合系数值为:
β0=1.1002,β1=0.0684,β2=0.0082,β3=0.1865,β4=2.6559,β5=0.0289,β6=1.6237,β7=0.0003,β8=0.4599,β9=0.5851,β10=0.4302;
所述第一无量纲系数λ的计算公式为:
具体的,根据上述表格内的不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,以及第一无量纲系数λ的表达式。使用Levenberg-Marquardt法对所述第一无量纲系数λ的表达式进行多项式拟合,得到第一无量纲系数λ的表达式中各拟合系数值,其中,β0=1.1002,β1=0.0684,β2=0.0082,β3=0.1865,β4=2.6559,β5=0.0289,β6=1.6237,β7=0.0003,β8=0.4599,β9=0.5851,β10=0.4302。
将上述各拟合系数值代入λ的表达式中,得到所述第一无量纲系数λ的计算公式为:
本发明实施例提供一种管道裂纹扩展驱动力计算方法,通过构建一个与裂纹长度比、裂纹深度比、管道径厚比和硬化指数有关,仅取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
基于上述实施例,可选的,在上述管道裂纹扩展驱动力计算方法中,所述根据所述第一无量纲系数,计算得到裂纹扩展驱动力,具体包括:
根据所述第一无量纲系数λ,计算参考应力σref,具体公式为:
σref=σ/λ;
其中,σ是远端施加应力;
根据所述参考应力σref,计算得到对应于σref的参考应变εref,进一步计算得到裂纹扩展驱动力;
参考应变εref计算公式为:
裂纹扩展驱动力计算公式为:
其中,α是硬化常数,σ0为屈服强度,ε0为屈服应变,n为应***化指数;
具体的,获取管道的裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n,并将上述数据代入上述第一无量纲系数λ的计算公式,计算得到第一无量纲系数λ。
根据所述第一无量纲系数λ,得到修正的极限载荷解PL的计算公式:
PL=2πrtσλ;
其中,r为管道半径;σ表示远端施加应力,它是物体外部施加载荷的大小与受载荷面积之比,单位为MPa。
通过上述计算参考应力σref的公式计算得到参考应力σref,单位为MPa;
其中,硬化常数α,杨氏模量E和I型应力强度因子KI的数值可以查阅标准“BS7910.Guide to methods for assessing the acceptability of flaws in metallicstructures”附录中的M4.1部分得到;其余计算使用数据均通过所述有限元模型模拟以获取。
本发明实施例提供一种管道裂纹扩展驱动力计算方法,通过构建一个与裂纹长度比、裂纹深度比、管道径厚比和硬化指数有关,仅取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
基于上述实施例,可选的,在上述管道裂纹扩展驱动力计算方法中,所述第一无量纲系数λ的计算公式适用于任何屈服强度σ0、10≤D/t≤30、0.1≤a/t≤0.5、0.05≤2c/πD≤0.20、10≤n≤25以及0.015≤εn≤0.04的情形;其中,εn为名义应变。
具体的,由于在本实施例中,用于进行参数化分析,通过参数化分析,构建第一无量纲系数λ的表达式;以及使用Levenberg-Marquardt法对所述第一无量纲系数λ的表达式进行多项式拟合,得到第一无量纲系数λ的表达式中各拟合系数值;使用的表格中的不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,均属于10≤D/t≤30、0.1≤a/t≤0.5、0.05≤2c/πD≤0.20、10≤n≤25的范围内。
因此,上述第一无量纲系数λ的计算公式仅适用于任何屈服强度σ0、10≤D/t≤30、0.1≤a/t≤0.5、0.05≤2c/πD≤0.20、10≤n≤25以及0.015≤εn≤0.04的情形;其中,εn为名义应变。
在实际情况中,也可更改上述范围,重新使用不同的数据进行参数化分析,构建第一无量纲系数λ的表达式,并对所述第一无量纲系数λ的表达式进行多项式拟合,得到第一无量纲系数λ的表达式中各拟合系数值,进而得到第一无量纲系数λ的计算公式。
本发明实施例提供一种管道裂纹扩展驱动力计算方法,通过构建一个与裂纹长度比、裂纹深度比、管道径厚比和硬化指数有关,仅取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
基于上述实施例,可选的,在上述管道裂纹扩展驱动力计算方法中,所述有限元模型采用ABAQUS软件或ANSYS软件构建。
具体的,目前流行的有限元分析软件主要有NASTRAN、ADINA、ANSYS、ABAQUS、MARC、MAGSOFT、COSMOS等,在本发明实施例中,优选采用ABAQUS软件或ANSYS软件构建限元模型。
本发明实施例提供一种管道裂纹扩展驱动力计算方法,通过构建一个与裂纹长度比、裂纹深度比、管道径厚比和硬化指数有关,仅取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
下面结合具体实例对本发明技术方案做进一步详述,管道本体受到远端均布拉伸载荷,在管道的厚度方向开有外表面椭圆形裂纹,裂纹面垂直于管道轴向,即承受轴向应变。
选取X80管线钢,以管道外径D=355.6mm,裂纹深度a=3.02mm,裂纹半长度c=100mm,管道壁厚t=15.9mm,屈服强度σ0=498MPa,硬化指数n=17.8,杨氏模量E=207000Gpa,材料泊松比v=0.3的含椭圆形外表面轴向裂纹的管道作为研究对象,假设此管道承受的远端施加应力σ=400MPa。
管道裂纹扩展驱动力计算方法,具体包括:
确定不同裂纹长2c/πD、裂纹深度a/t、管道径厚比D/t等参数:
其中,2c/πD=0.179,a/t=0.190,D/t=22.365;
将上述参数代入第一无量纲系数λ的计算公式中,计算得到第一无量纲系数λ具体数值;
根据第一无量纲系数λ,以及计算参考应力σref的公式σref=σ/λ,计算得到参考应力σref=σ/λ=400/0.886=451.572MPa。
因此,当σ=400MPa时的裂纹扩展驱动力:
图4为本发明实施例提供的管道裂纹扩展驱动力计算***结构示意图,如图4所示,管道裂纹扩展驱动力计算***,包括:
模型设置模块410,用于在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
数据提取模块420,与模型设置模块410相连,用于对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;
参数计算模块430,与数据提取模块420相连,用于通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
数据分析模块440,与参数计算模块430相连,用于根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
具体的,模型设置模块410,用于通过有限元分析软件材料属性模块在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并通过有限元分析软件载荷模块设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
图2为本发明实施例提供的含外表面椭圆形裂纹的管道有限元模型示意图,如图2所示,所述管道本体受到远端均布拉伸载荷,在管道的厚度方向开有外表面椭圆形裂纹,裂纹面垂直于管道轴向,即承受轴向应变。
数据提取模块420,与模型设置模块410相连,用于对所述有限元模型划分网格,有限元分析软件交互模块设置裂纹和输出参量,进行有限元计算,并通过有限元分析软件处理模块在所述有限元计算后提取数据;
参数计算模块430,与数据提取模块420相连,用于通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
数据分析模块440,与参数计算模块430相连,用于根据所述第一无量纲系数,以及相关裂纹扩展驱动力计算所需参数,计算得到裂纹扩展驱动力。
在上述实施例的基础上,在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件之前,还包括通过有限元分析软件构建所述管道的有限元模型;需要说明的是,目前流行的有限元分析软件主要有NASTRAN、ADINA、ANSYS、ABAQUS、MARC、MAGSOFT、COSMOS等,在本发明实施例中,不对构建限元模型的有限元分析软件进行限定,具体使用的有限元分析软件可根据实际情况进行选择。
本发明实施例提供的管道裂纹扩展驱动力计算***用于执行在上述管道裂纹扩展驱动力计算方法中,其具体的实施方式与方法实施方式一致,此处不再赘述。
本发明实施例提供一种管道裂纹扩展驱动力计算***,通过构建一个取决于裂纹尺寸和材料性能的第一无量纲系数,将含裂纹结构件的参考应力直接与远端施加应力相关联,简化了裂纹扩展驱动力计算方法,并能够精准地计算裂纹扩展驱动力,方便了对承受塑性变形管道进行断裂评估,提高了评估结果的准确性。
图5为本发明实施例提供的电子设备的实体结构示意图,如图5所示,所述电子设备可以包括:处理器(processor)510、通信接口(Communications Interface)520、存储器(memory)530和通信总线540,其中,处理器510,通信接口520,存储器530通过通信总线540完成相互间的通信。处理器510可以调用存储器530中的逻辑指令,以执行上述管道裂纹扩展驱动力计算方法的各个步骤。例如:在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
此外,上述的存储器530中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面,本发明实施例还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行,用以实现执行上述各实施例提供的管道裂纹扩展驱动力计算方法。例如包括:在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种管道裂纹扩展驱动力计算方法,其特征在于,包括,
在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;
通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
2.根据权利要求1所述的管道裂纹扩展驱动力计算方法,其特征在于,
所述通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数,具体包括:
通过所述有限元模型模拟以获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,计算得到第一无量纲系数;
其中,c为裂纹半长度;D为管道外径;a为裂纹深度;t为管道壁厚。
3.根据权利要求2所述的管道裂纹扩展驱动力计算方法,其特征在于,
所述通过所述有限元模型模拟以获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据,计算得到第一无量纲系数,具体包括:
获取不同裂纹长度比2c/πD、裂纹深度比a/t、管道径厚比D/t和硬化指数n情况下的第一无量纲系数λ数据;
通过参数化分析,构建第一无量纲系数λ的表达式:
使用Levenberg-Marquardt法对所述第一无量纲系数λ的表达式进行多项式拟合,得到第一无量纲系数λ的表达式中各拟合系数值;
将所述各拟合系数值代入第一无量纲系数λ的表达式中,得到所述第一无量纲系数λ的计算公式;
根据所述第一无量纲系数的计算公式,计算得到第一无量纲系数。
6.根据权利要求4所述的管道裂纹扩展驱动力计算方法,其特征在于,
所述第一无量纲系数λ的计算公式适用于任何屈服强度σ0、10≤D/t≤30、0.1≤a/t≤0.5、0.05≤2c/πD≤0.20、10≤n≤25以及0.015≤εn≤0.04的情形;其中,εn为名义应变。
7.根据权利要求1所述的管道裂纹扩展驱动力计算方法,其特征在于,
所述有限元模型采用ABAQUS软件或ANSYS软件构建。
8.一种管道裂纹扩展驱动力计算***,其特征在于,包括:
模型设置模块,用于在预先建立的管道的有限元模型中设置相应的弹塑性本构关系,并设置拉伸载荷及边界条件;所述边界条件包括对称约束和固定约束;
数据提取模块,与模型设置模块相连,用于对所述有限元模型划分网格,设置裂纹和输出参量,进行有限元计算,并在所述有限元计算后提取数据;
参数计算模块,与数据提取模块相连,通过所述有限元模型模拟得到的数据,计算得到第一无量纲系数;所述第一无量纲系数取决于裂纹尺寸和材料性能,是与修正的极限载荷解相关的参数;
数据分析模块,与参数计算模块相连,用于根据所述第一无量纲系数,计算得到裂纹扩展驱动力。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至7任一项所述管道裂纹扩展驱动力计算方法的步骤。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现如权利要求1至7任一项所述管道裂纹扩展驱动力计算方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010313384.8A CN111539142B (zh) | 2020-04-20 | 2020-04-20 | 一种管道裂纹扩展驱动力计算方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010313384.8A CN111539142B (zh) | 2020-04-20 | 2020-04-20 | 一种管道裂纹扩展驱动力计算方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111539142A true CN111539142A (zh) | 2020-08-14 |
CN111539142B CN111539142B (zh) | 2023-12-22 |
Family
ID=71979165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010313384.8A Active CN111539142B (zh) | 2020-04-20 | 2020-04-20 | 一种管道裂纹扩展驱动力计算方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111539142B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111914372A (zh) * | 2020-08-17 | 2020-11-10 | 北京三一智造科技有限公司 | 裂纹扩展深度计算方法、装置及电子设备 |
CN113297770A (zh) * | 2021-06-04 | 2021-08-24 | 中国核动力研究设计院 | 核电厂压力管道泄漏参数获取方法及泄漏模型和应用 |
CN113378424A (zh) * | 2021-04-29 | 2021-09-10 | 中国石油天然气集团有限公司 | 一种评估管道环焊接头承受轴向变形能力的方法及*** |
CN116663146B (zh) * | 2023-05-30 | 2023-11-17 | 西安理工大学 | 非圆形管道沿程阻力的计算方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016045024A1 (zh) * | 2014-09-25 | 2016-03-31 | 华东理工大学 | 一种高温环境下的结构材料断裂韧性测量与确定方法 |
CN108710721A (zh) * | 2018-04-03 | 2018-10-26 | 苏州热工研究院有限公司 | 一种基于失效评定图的未爆先漏分析方法 |
CN108920792A (zh) * | 2018-06-21 | 2018-11-30 | 北京工业大学 | 一种基于小裂纹扩展的搅拌摩擦焊构件疲劳寿命预测方法 |
CN109255139A (zh) * | 2017-07-14 | 2019-01-22 | 天津大学 | 一种高温管道中表面多裂纹合并方法 |
CN109759737A (zh) * | 2019-01-29 | 2019-05-17 | 中国电力科学研究院有限公司 | 含气孔缺陷的gil铝合金焊接管道安全评定方法 |
CN110516269A (zh) * | 2018-05-22 | 2019-11-29 | 天津大学 | 一种用于海底管道工程临界评估的优化参考应变方法 |
WO2020042753A1 (zh) * | 2018-08-30 | 2020-03-05 | 长沙理工大学 | 季节性腐蚀与疲劳耦合作用下钢筋混凝土桥梁寿命预测方法 |
CN110895641A (zh) * | 2018-09-07 | 2020-03-20 | 天津大学 | 一种含有环向表面裂纹v型坡口焊缝海底冶金复合管工程临界评估方法 |
-
2020
- 2020-04-20 CN CN202010313384.8A patent/CN111539142B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016045024A1 (zh) * | 2014-09-25 | 2016-03-31 | 华东理工大学 | 一种高温环境下的结构材料断裂韧性测量与确定方法 |
CN109255139A (zh) * | 2017-07-14 | 2019-01-22 | 天津大学 | 一种高温管道中表面多裂纹合并方法 |
CN108710721A (zh) * | 2018-04-03 | 2018-10-26 | 苏州热工研究院有限公司 | 一种基于失效评定图的未爆先漏分析方法 |
CN110516269A (zh) * | 2018-05-22 | 2019-11-29 | 天津大学 | 一种用于海底管道工程临界评估的优化参考应变方法 |
CN108920792A (zh) * | 2018-06-21 | 2018-11-30 | 北京工业大学 | 一种基于小裂纹扩展的搅拌摩擦焊构件疲劳寿命预测方法 |
WO2020042753A1 (zh) * | 2018-08-30 | 2020-03-05 | 长沙理工大学 | 季节性腐蚀与疲劳耦合作用下钢筋混凝土桥梁寿命预测方法 |
CN110895641A (zh) * | 2018-09-07 | 2020-03-20 | 天津大学 | 一种含有环向表面裂纹v型坡口焊缝海底冶金复合管工程临界评估方法 |
CN109759737A (zh) * | 2019-01-29 | 2019-05-17 | 中国电力科学研究院有限公司 | 含气孔缺陷的gil铝合金焊接管道安全评定方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111914372A (zh) * | 2020-08-17 | 2020-11-10 | 北京三一智造科技有限公司 | 裂纹扩展深度计算方法、装置及电子设备 |
CN111914372B (zh) * | 2020-08-17 | 2023-10-31 | 北京三一智造科技有限公司 | 裂纹扩展深度计算方法、装置及电子设备 |
CN113378424A (zh) * | 2021-04-29 | 2021-09-10 | 中国石油天然气集团有限公司 | 一种评估管道环焊接头承受轴向变形能力的方法及*** |
CN113297770A (zh) * | 2021-06-04 | 2021-08-24 | 中国核动力研究设计院 | 核电厂压力管道泄漏参数获取方法及泄漏模型和应用 |
CN116663146B (zh) * | 2023-05-30 | 2023-11-17 | 西安理工大学 | 非圆形管道沿程阻力的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111539142B (zh) | 2023-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111539142A (zh) | 一种管道裂纹扩展驱动力计算方法及*** | |
Zhao et al. | Tensile strength of Q345 steel with random pitting corrosion based on numerical analysis | |
CN110795879B (zh) | 一种连续的钢管混凝土构件压弯承载力计算方法 | |
Yan et al. | Ductile fracture simulation of constructional steels based on yield-to-fracture stress–strain relationship and micromechanism-based fracture criterion | |
CN109214020B (zh) | 一种储罐弹塑性象足屈曲临界载荷获取方法及装置 | |
Sellier et al. | Weakest link and localisation WL2: a method to conciliate probabilistic and energetic scale effects in numerical models | |
Wang et al. | An Experimental‐Numerical Combined Method to Determine the True Constitutive Relation of Tensile Specimens after Necking | |
CN111539144A (zh) | 一种含裂纹结构件的断裂韧性计算方法及*** | |
Krishan et al. | Strength calculation of short concrete-filled steel tube columns | |
Zheng et al. | Arch-dam crack deformation monitoring hybrid model based on XFEM | |
CN114398805B (zh) | 一种水岩耦合作用下裂隙岩石蠕变模型构建方法及*** | |
JP2007163392A (ja) | 鋼管の局部座屈性能評価方法、鋼管の材質設計方法、鋼管 | |
CN107066728B (zh) | 一种钛合金潜水器耐压球壳极限承载力估算方法 | |
JP7292755B2 (ja) | 弾塑性解析方法及び弾塑性解析プログラム | |
Jeong et al. | Fully plastic crack opening analyses of complex-cracked pipes for Ramberg-Osgood materials | |
CN111157338B (zh) | 一种金属哑铃型试样大应变范围压缩硬化曲线的测量方法 | |
CN111177842B (zh) | 一种盾构隧道管片接触面的动力特性确定方法、装置及存储介质 | |
Young et al. | Advances in COD Equations—Circumferential Through-Wall Cracks | |
Halicka et al. | Designing of cylindrical concrete tanks with regard to buckling and second order effects | |
CN111797553B (zh) | 复杂荷载条件下腐蚀管道的极限内压预测方法及装置 | |
Cho et al. | Engineering estimates of crack opening displacement for non-idealized circumferential through-wall cracks in pipe | |
CN115221581B (zh) | 一种基于不同承载力指标的剪力墙损伤参数确定方法 | |
CN116384201A (zh) | 橡胶本构模型参数的确定方法、装置、设备及存储介质 | |
Fu et al. | Calculation and improvement of complex mechanical properties of the mega arch dam under multiple loads based on finite element methods | |
CN115493787A (zh) | 带随机分布脱空缺陷的钢管混凝土结构冲击性能模拟方法 |
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 |