应力集中部位裂纹自由表面处应力强度因子的计算方法技术

本发明专利技术属于高压及超高压容器的设计开发处理及失效评定领域，具体涉及一种应力集中部位裂纹自由表面处应力强度因子的计算方法。本发明专利技术包括：根据容器结构及载荷参数，对应力集中部位进行弹性应力分析；根据分析的结果，提取出垂直于裂纹所在平面的应力分布数据；把应力分布数据进行离散化，获得离散化数据，各段离散化数据分别进行多项式求解，获得相应段离散化数据的应力分布函数；进行裂纹形状系数的计算；进行裂纹自由表面处应力强度因子的计算。本发明专利技术有效提升了应力分布数据的拟合精度，进而确保了计算结果的准确性；同时，计算方法也非常高效简洁，工程计算中的实用性高。工程计算中的实用性高。工程计算中的实用性高。

【技术实现步骤摘要】
应力集中部位裂纹自由表面处应力强度因子的计算方法


[0001]本专利技术属于高压及超高压容器的设计开发处理及失效评定领域，具体涉及一种应力集中部位裂纹自由表面处应力强度因子的计算方法
。

技术介绍

[0002]高压及超高压容器的应力集中部位，是在制造及使用运行过程中最易损伤的地方，更是最易出现表面裂纹缺陷的区域，通常表现为半椭圆形或半圆形裂纹；
GB/T 34019—2017《
超高压容器
》
中明确表述的开孔处裂纹（
B
型）和盲底裂纹（
D
型），便是此类应力集中部位中出现的典型裂纹
。
在高压及超高压容器的设计开发处理及失效评定过程中，需要准确地计算出裂纹尖端的应力强度因子，因为它在含裂纹压力容器剩余强度评定和裂纹扩展剩余寿命计算中是不可或缺的关键参数；尤其对于某些特定裂纹，评定过程中更需重点考察自由表面处以及最深点处应力强度因子的值
。
因而，在实际工程应用中，找到一种简洁又能保证计算精度的应力集中部位裂纹自由表面和最深点其中之一处的应力强度因子的计算方法，显得尤为重要
。
[0003]当前，裂纹尖端的应力强度因子的计算方法主要有数学分析法
、
有限元法
、
边界配置法及光弹性法等，其中以有限元法居多
。
对于高压及超高压容器典型裂纹应力强度因子的计算，
ASME BPVC.
Ⅷ
.3
‑
2021《AlternativeRules for Construction of High Pressure Vessels》
中的非强制性附录
D
和
GB 34019—2017《
超高压容器
》
中的附录
F
中都有所提及，两个业内常见标准也都给出了
A
型裂纹（筒体内壁轴
‑
径向裂纹）的详细计算步骤，即：将应力分布数据进行拟合后，再根据拟合系数及裂纹形状等，进而计算出相应的应力强度因子值
。
但，对于开孔处裂纹（
B
型）和盲底裂纹（
D
型），标准仅是提到可以参照
A
型裂纹计算方法得出，而并未见具体计算流程
。
[0004]然而，实际高压及超高压容器应力集中部位的应力分布特点为：应力梯度变换大，且衰减由快变慢，即应力起初会随着从裂纹表面所测得距离的加大呈“快速下降”趋势，到一定的距离后再呈“平稳下降”趋势
。
因此，如果惯性地按照传统的
A
型裂纹计算方法，由于上述大幅度变化的应力梯度因素影响，其应力分布数据往往无法进行合适的拟合，或者拟合出来的曲线与实际数据相差甚远，随之导致应力集中部位自由表面处应力强度因子的计算结果必然会和实际值有较大偏差，从而无法作为评定过程中的准确判据
。
当然，此时可以换用数值分析方法对含裂纹结构进行专门的断裂力学分析，但计算过程会很繁琐且收敛难度大，造成的计算成本很高，这也是目前极少采用这类断裂力学分析的主因之一
。
因而，是否能研发出一种针对应力集中部位裂纹自由表面处应力强度因子计算方法，为本领域近年来所亟待解决的技术难题
。

技术实现思路

[0005]本专利技术的目的是克服上述现有技术的不足，提供一种应力集中部位裂纹自由表面处应力强度因子的计算方法，本专利技术有效提升了应力分布数据的拟合精度，进而确保了计
算结果的准确性；同时，计算方法也非常高效简洁，工程计算中的实用性高
。
[0006]为实现上述目的，本专利技术采用了以下技术方案：应力集中部位裂纹自由表面处应力强度因子的计算方法，其特征在于包括以下步骤：
S1.
根据容器结构及载荷参数，对应力集中部位进行弹性应力分析；
S2.
根据分析的结果，提取出垂直于裂纹所在平面的应力分布数据；
S3.
把应力分布数据进行离散化，获得离散化数据，各段离散化数据分别进行多项式求解，获得相应段离散化数据的应力分布函数；
S4.
进行裂纹形状系数的计算；
S5.
进行裂纹自由表面处应力强度因子的计算
。
[0007]优选的，步骤
S3
中，每三组彼此相邻的应力分布数据形成一段离散化数据，各相邻段离散化数据之间彼此衔接，以保证连续性；具体包括以下子步骤：
S31.
按下式求取第
i
段离散化数据的多项式系数
a
0i
、a
1i
及
a
2i
：
[0008][0009][0010]式中：为步骤
S2
中所提取的应力分布数据中第组对应的应力值，单位
MPa
；为步骤
S2
中所提取的应力分布数据中第组对应的应力值，单位
MPa
；为步骤
S2
中所提取的应力分布数据中第组对应的应力值，单位
MPa
；为步骤
S2
中所提取的应力分布数据中第组对应的深度值，单位
mm
；为步骤
S2
中所提取的应力分布数据中第组对应的深度值，单位
mm
；为步骤
S2
中所提取的应力分布数据中第组对应的深度值，单位
mm
；
S32.
以下式进行第
i
段离散化数据的应力分布函数
y
i
(x)
的求解：
[0011]式中：
x
为裂纹自由表面起所测得的距离，也即深度，为变量，，单位
mm
；
n
为离散化的总段数
。
[0012]优选的，步骤
S5
包括以下子步骤：
S51.
进行分应力强度因子
K
I,0
～
K
I,4
的计算：
[0013][0014][0015]式中：为椭圆形裂纹深度值，为常量，单位
mm
；
S52.
以下式计算裂纹自由表面处应力强度因子：式中：
p
i
为容器的内压力，单位
MPa
；
Q
为裂纹形状系数；
M
1B
～
M
4B
均为中间系数，计算公式如下：
[0016]G0～
G3均为裂纹自由表面处系数，按
GB 34019—2017《
超高压容器
》
中的表
F.2
取值，在表
F.2
所给出值之间的取插值
。
[0017]优选的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
应力集中部位裂纹自由表面处应力强度因子的计算方法，其特征在于包括以下步骤：
S1.
根据容器结构及载荷参数，对应力集中部位进行弹性应力分析；
S2.
根据分析的结果，提取出垂直于裂纹所在平面的应力分布数据；
S3.
把应力分布数据进行离散化，获得离散化数据，各段离散化数据分别进行多项式求解，获得相应段离散化数据的应力分布函数；
S4.
进行裂纹形状系数的计算；
S5.
进行裂纹自由表面处应力强度因子的计算
。2.
根据权利要求1所述的应力集中部位裂纹自由表面处应力强度因子的计算方法，其特征在于：步骤
S3
中，每三组彼此相邻的应力分布数据形成一段离散化数据，各相邻段离散化数据之间彼此衔接，以保证连续性；具体包括以下子步骤：
S31.
按下式求取第
i
段离散化数据的多项式系数
a
0i
、a
1i
及
a
2i
：：式中：为步骤
S2
中所提取的应力分布数据中第组对应的应力值，单位
MPa
；为步骤
S2
中所提取的应力分布数据中第组对应的应力值，单位
MPa
；为步骤
S2
中所提取的应力分布数据中第组对应的应力值，单位
MPa
；为步骤
S2
中所提取的应力分布数据中第组对应的深度值，单位
mm
；为步骤
S2
中所提取的应力分布数据中第组对应的深度值，单位
mm
；为步骤
S2
中所提取的应力分布数据中第组对应的深度值，单位
mm
；
S32.
以下式进行第
i
段离散化数据的应力分布函数
y
i
(x)
的求解：式中：
x
为裂纹自由表面起所测得的距离，也即深度，为变量，，单位

【专利技术属性】
技术研发人员：汪志福，范志超，周煜，秦宗川，牛铮，范海俊，戴兴旺，危书涛，朱金花，
申请(专利权)人：合肥通用机械研究院有限公司，
类型：发明
国别省市：