一种预测高钢级管道环焊缝极限应变能力方法及装置制造方法及图纸

本发明专利技术提供了一种预测高钢级管道环焊缝极限应变能力的方法及装置，其特征在于，所述方法包括：建立近场动力学

【技术实现步骤摘要】
一种预测高钢级管道环焊缝极限应变能力方法及装置


[0001]本专利技术涉及弹塑性金属材料
，尤其涉及一种预测高钢级管道环焊缝极限应变能力的方法及装置。

技术介绍

[0002]为了提高天然气输量及管道输送的经济性，管道工程正朝着高钢级、大口径、高压输送的方向发展。随着X80级以上高钢级管道的应用，以及在役老旧管道越来越接近失效的高发期，环焊缝开裂成为管道失效的主要因素，焊接过程导致环焊缝区域材料性能劣化、产生焊接裂纹等缺陷，当管道穿越地质灾害多发地区时，在位移作用下引起局部应变集中而开裂，制约了管道的应变能力。如何评价环焊缝的极限应变能力成为管道断裂控制中的新问题。当前针对管道应变能力的经验公式大多基于有限元计算，但传统有限元方法求解裂纹扩展不连续问题时无法有效解决裂纹尖端奇异性，且需额外引入复杂的外部失效准则判断断裂发生，计算效率低，难度大，影响了应变能力的评估。需提出一种适用于高钢级管道环焊缝裂纹扩展失效分析的数值计算方法，这也是高钢级天然气管道面临的亟待解决的科学问题。
[0003]近场动力学方法(Peridynamic，PD)，是近二十年新兴的一种非局部、无网格粒子、积分形式的计算方法。该方法把材料模型离散成带有物理信息的物质点，使模型的本构关系在介质不连续处仍然有效。该理论认为每个材料点与其近场域内其他材料点之间存在相互作用。PD理论最显著的优点是：避免了基于连续理论的传统模型在求解不连续问题时的奇异性和复杂性，且损伤和断裂是作为模拟的自然结果被直接包含在本构关系中的。因此，随着模拟过程的进行，产生裂纹萌生、扩展现象，不需要任何的额外准则，可以模拟裂纹萌生和扩展的全过程。PD方法发展至今已在岩石、复合纤维脆性和微弹性材料损伤断裂领域得到了应用，基于键(bond
‑
based theory)伸长的失效准则认为：粒子间的伸长量和力密度的关系呈线性，当变形后的伸长量达到极限伸长量s0时，则判断为失效，此过程不可逆。然而对于大多数金属来说，损伤常伴随着塑性变形，与弹性变形不同，塑性变形依赖于加载历史，外力和变形不再服从线性关系，如图1，弹塑性本构模型中的失效准则不能再简单的使用最终状态的伸长量来判断。
[0004]现有技术将PD方法与有限元方法(Finite element method，FEM)进行耦合，并在有限元软件Abaqus中运行。耦合模型中将两个材料点之间的键看成是杆单元，材料点与其近场域内所有点之间将生成大量的杆单元，导致输入文件过大，影响运算效率；键基近场动力学模型将二维和三维结构的泊松比取为固定值1/3和1/4，限制了材料的本构关系。近场动力学方法避免了有限元方法在求解裂纹扩展等不连续问题的弊端，能够模拟裂纹扩展的真实过程，但其应用领域更多是岩石和纤维等脆性材料，对于金属材料尤其是高钢级长输管道的开裂失效鲜有研究近场动力学方法在计算效率以及弹塑性损伤失效准则方面仍存在不足，制约了其在弹塑性金属材料裂纹扩展分析中的应用。

技术实现思路

[0005]鉴于上述问题，提出了本专利技术以便提供一种克服上述问题或者至少部分地解决上述问题的一种预测高钢级管道环焊缝极限应变能力的方法及装置。
[0006]本专利技术的一个方面，提供了一种预测高钢级管道环焊缝极限应变能力的方法，所述方法包括：
[0007]建立近场动力学
‑
有限元耦合模型；
[0008]施加载荷边界条件；
[0009]基于有限元软件ABAQUS中有限元用户动态子程序模块，利用用户自定义子程序实现提出的近场动力学
‑
有限元耦合模型，对所述近场动力学
‑
有限元耦合模型进行计算，预测高钢级管道环焊缝极限应变能力。
[0010]进一步地，所述建立近场动力学
‑
有限元耦合模型，包括：根据节点、单元和材料属性，建立PD弹塑性本构模型；
[0011]配置失效准则。
[0012]进一步地，所述建立PD弹塑性本构模型，包括：屈服面表达式如下：
[0013][0014][0015]式中，为i点的等效应力，为屈服应力，σ
y0
为初始屈服应力，K为等向强化模量；
[0016]表达式如下：
[0017][0018][0019][0020]式中，μ为剪切模量，δ为近场域半径，为材料点i和j之间键的弹性伸长量，|ξ|为材料点i和j之间的距离，V
(j)
为材料点j的体积，a、a
μ
、a
k
和均为近场动力学参数，E为弹性模量，v为泊松比，N为近场域节点数，h为二维结构的厚度，
[0021]θ
(i)
为体积膨胀，表达式如下：
[0022][0023]式中，和Λ
(i)(j)
均为近场动力学参数，
[0024]s
(i)(j)
为材料点i和j之间的伸长，表达式如下：
[0025][0026][0027][0028]式中，y
(j)
为材料点j的坐标，y
(i)
为材料点i的坐标，M<ξ>为变形矢量方向；等效塑性拉伸表达式如下：
[0029][0030][0031][0032]式中，为前一分析步等效塑性拉伸应变，为等效塑性拉伸应变增量，A0为系数，为塑性拉伸应变增量。
[0033]进一步地，所述失效准则按以下方式进行配置：
[0034]将消除两个粒子间作用所需的极限能量释放率作为判别粒子间作用是否失效的标准，表达式如下：
[0035][0036]式中，为两个粒子间的能量释放率，g
c
为消除两个粒子间作用所需的极限能量释放率，判断两个粒子间作用消失后，引入函数ψ
(i)(j)
判断两个粒子间作用失效与否，函数ψ
(i)(j)
取值0或1，表达式如下：
[0037][0038]g
c
计算式如下：
[0039][0040]式中，Δx为PD模型单元尺寸，h为二维结构厚度，G
c
为材料的极限能量释放率，N
c
为
近场动力学参数，w
(i)(j)
为材料点i和j之间的能量密度，由力密度和键伸长量关系曲线所包围的面积；
[0041]材料点x
i
的局部损伤值为：
[0042][0043][0044]式中，H
i
为材料点i的近场域范围。
[0045]进一步地，所述施加载荷边界条件，包括：
[0046]在原始域R的边界定义一个非零体积域R
c
，所述非零体积域宽度与近场域尺寸δ相同，外部载荷以体积力的形式施加到虚拟层R
c
的材料点。
[0047]进一步地，所述近场动力学
‑
有限元耦合模型受到的外部载荷的表达式为：
[0048][0049本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种预测高钢级管道环焊缝极限应变能力方法，其特征在于，所述方法包括：建立近场动力学
‑
有限元耦合模型；施加载荷边界条件；基于有限元软件ABAQUS中有限元用户动态子程序模块，利用用户自定义子程序实现提出的近场动力学
‑
有限元耦合模型，对所述近场动力学
‑
有限元耦合模型进行计算，预测高钢级管道环焊缝极限应变能力。2.根据权利要求1所述的方法，其特征在于，所述建立近场动力学
‑
有限元耦合模型，包括：根据节点、单元和材料属性，建立PD弹塑性本构模型；配置失效准则。3.根据权利要求2所述的方法，其特征在于，所述建立PD弹塑性本构模型，包括：屈服面表达式如下：表达式如下：式中，为i点的等效应力，为屈服应力，σ
y0
为初始屈服应力，K为等向强化模量；表达式如下：表达式如下：表达式如下：式中，μ为剪切模量，δ为近场域半径，为材料点i和j之间键的弹性伸长量，|ξ|为材料点i和j之间的距离，V
(j)
为材料点j的体积，a、a
μ
、a
k
和均为近场动力学参数，E为弹性模量，v为泊松比，N为近场域节点数，h为二维结构的厚度，θ
(i)
为体积膨胀，表达式如下：式中，和Λ
(i)(j)
均为近场动力学参数，s
(i)(j)
为材料点i和j之间的伸长，表达式如下：
式中，y
(j)
为材料点j的坐标，y
(i)
为材料点i的坐标，M<ξ>为变形矢量方向；等效塑性拉伸表达式如下：表达式如下：表达式如下：式中，为前一分析步等效塑性拉伸应变，为等效塑性拉伸应变增量，A0为系数，为塑性拉伸应变增量。4.根据权利要求2所述的方法，其特征在于，所述失效准则按以下方式进行配置：将消除两个粒子间作用所需的极限能量释放率作为判别粒子间作用是否失效的标准，表达式如下：式中，为两个粒子间的能量释放率，g
c
为消除两个粒子间作用所需的极限能量释放率，判断两个粒子间作用消失后，引入函数ψ
(i)(j)
判断两个粒子间作用失效与否，函数ψ
(i)(j)
取值0或1，表达式如下：g
c
计算式如下：式中，Δx为PD模型单元尺寸，h为二维结构厚度，G
c
为材料的极限能量释放率，N
c
为近场动力学参数，w
(i)(j)
为材料点i和j之间的能量密度，由力密度和键伸长量关系曲线所包围的面积；材料点x
i
的局部损伤值为：
式中，H
i
为材料点i的近场域范围。5.根据权利要求1所述的方法，其特征在于，所述施加载荷边界条件，包括：在原始域R的边界定义一个非零体积域R
c
，所述非零体积域宽度与近场域尺寸δ相同，外部载荷以体积力的形式施加到虚拟层R
c
的材料点。6.根据权利要求5所述的方法，其特征在于，所述近场动力学
‑
有限元耦合模型受到的外部...

【专利技术属性】
技术研发人员：田骁，裴晶晶，季淮君，周凡，王泺，王晋，路航，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：