一种疲劳短裂纹扩展速率预测模型的确定方法及装置与流程

1.本技术涉及疲劳裂纹扩展速率预测技术领域，尤其涉及一种疲劳短裂纹扩展速率预测模型的确定方法及装置。


背景技术：

2.随着我国经济水平及科学技术的飞速发展，对轻量化、高寿命化、大型化的海洋船舶、石油开采等工程结构件的需求日益增加。但随着结构件增加，疲劳断裂的危害性也随之增加。而小于1.5mm的疲劳短裂纹约通常占据结构件疲劳总寿命的50
‑
80％，如果采用传统paris公式则会获得非保守的疲劳寿命，带来一定的风险性。通常通过消除裂纹闭合效应获得有效应力强度因子δk
eff
计算疲劳短裂纹扩展速率。
3.然而疲劳短裂纹通常分为微观结构短裂纹(msc)及物理短裂纹(psc)，并且在msc的裂纹扩展速率会受到极大的晶界阻碍，呈现递减趋势，但目前不同材料具有不同的msc，而psc阶段的裂纹扩展速率几乎不会受到晶界阻碍，呈现增长趋势，致使疲劳短裂纹呈现一个先减后增的v形趋势，但目前尚未存在有效应力强度因子δk
eff
描述短裂纹v形扩展的模型。


技术实现要素：

4.有鉴于此，本发明提出一种基于δk
eff
来预测疲劳短裂纹扩展速率的分段模型，为检检测海洋平台结构件提供重要理论依据及工业价值，具有重要意义。
5.为解决上述技术问题，本说明书实施例是这样实现的：
6.本说明书实施例提供的一种疲劳短裂纹扩展速率预测模型的确定方法，其特征在于，所述方法包括：
7.对材料进行不同应力范围条件下疲劳短裂纹扩展实验，记录实验数据；
8.根据实验数据计算有效应力强度因子δk
eff
和疲劳短裂纹扩展速率da/dn；
9.计算微观结构短裂纹长度对应的有效应力强度因子δk
eff
‑
msc
；
10.基于δk
eff
和δk
eff
‑
msc
确定第一预测疲劳断裂纹速度模型；
11.确定第一预测疲劳断裂纹速度模型中的拟合常数，带入所述第一预测疲劳断裂纹速度模型中。
12.可选的，所述根据实验数据计算有效应力强度因子δk
eff
，具体包括：
13.根据实验数据计算应力强度因子δk；
[0014][0015]
其中δk为应力强度因子范围，y为形状因子需按照实际试样所对应标准进行计算，δσ为应力，a为表面裂纹长度；
[0016]
根据所述应力强度因子计算有效应力强度因子δk
eff
；
[0017]
[0018]
其中δk
eff
为有效应力强度因子范围，σ
op
为裂纹张开应力，σ
max
为最大应力，r为应力比。
[0019]
可选的，所述根据实验数据计算疲劳短裂纹扩展速率da/dn，具体包括：
[0020][0021]
其中，δn为循环次数增加数值，δa为表面裂纹增加长度，a
i 1
为循环次数为n
i 1
时表面裂纹长度。
[0022]
可选的，所述计算微观结构短裂纹长度对应的有效应力强度因子δk
eff
‑
msc
，具体包括：
[0023]
根据应力δσ和微观结构短裂纹长度计算微观结构短裂纹长度a
msc
所对应的应力强度因子δk
msc
；
[0024]
根据δk
msc
微观结构短裂纹长度a
msc
对应的有效应力强度因子δk
eff
‑
msc
。
[0025]
可选的，
[0026]
其中，微观结构短裂纹长度a
msc
的长度是通过da/dn
‑
a趋势图定义的，a
msc
为裂纹扩展速率递减与增加的临界点，δk
msc
为a
msc
对应的应力强度因子，y为形状因子需按照实际试样所对应标准进行计算，δσ为应力；
[0027]
可选的，
[0028]
其中，δk
eff
‑
msc
为微观结构短裂纹长度a
msc
对应的有效应力强度因子范围，σ
op
为裂纹张开应力，σ
max
为最大应力，r为应力比。
[0029]
可选的，所述基于δk
eff
和δk
eff
‑
msc
确定第一预测疲劳断裂纹速度模型，具体包括：
[0030]
基于δk
eff
和δk
eff
‑
msc
，计算不同δσ下疲劳短裂纹扩展速率：
[0031]
da/dn＝c|δk
eff
‑
δk
eff
‑
msc
|
m
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0032]
其中c和m为拟合常数，其大小随δσ数值大小增加而增加。
[0033]
可选的，所述确定第一预测疲劳断裂纹速度模型中的拟合常数，具体包括：
[0034]
将获得的不同δσ下da/dn
‑
δk
eff
数据在双对数坐标系中绘图，并对其按照公式(6)进行非线性曲线拟合，获得不同δσ下对应的拟合常数c、m，后分别在双对数坐标系中绘制c
‑
δσ和m
‑
δσ并对其线性拟合，得到公式(7)和(8)：
[0035]
c＝10
a
(δσ)
b
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0036]
m＝10
x
(δσ)
y
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0037]
其中a，b和x，y分别为拟合常数
[0038]
可选的，利用公式(7)和(8)通过线性外推法分别获得σ0条件下所对应的拟合常数和最终确定公式(9)为描述材料不同应力范围条件下疲劳短裂纹v形扩展速率模型：
[0039][0040]
其中，应力归一化系数β为δσ/σ0；σ0为流变应力，取材料屈服强度与抗拉强度之和
的平均值。
[0041]
本说明书实施例提供的一种疲劳短裂纹扩展速率预测模型的确定装置，所述装置包括：
[0042]
应力实验模块，用于对材料进行不同应力范围条件下疲劳短裂纹扩展实验，记录实验数据；
[0043]
第一计算模块，用于根据实验数据计算有效应力强度因子δk
eff
和疲劳短裂纹扩展速率da/dn；
[0044]
第二计算模块，用于计算微观结构短裂纹长度对应的有效应力强度因子δk
eff
‑
msc
；
[0045]
模型确定模块，用于基于δk
eff
和δk
eff
‑
msc
确定第一预测疲劳断裂纹速度模型；
[0046]
常数确定模块，用于确定第一预测疲劳断裂纹速度模型中的拟合常数，带入所述第一预测疲劳断裂纹速度模型中。
[0047]
本说明书实施例采用的上述至少一个技术方案能够达到以下有益效果：
[0048]
本发明基于有效应力强度因子δk
eff
及有效应力强度因子δk
eff
‑
msc
及paris公式提出分段函数精准描述某种材料疲劳短裂纹v形扩展趋势，整个测试过程精度高、效率高，应用范围广，具有广阔应用前景和经济效益，海洋工程材料性能评估提供重要理论依据及实践基础。
附图说明
[0049]
此处所说明的附图用来提供对本技术的进一步理解，构成本技术的一部分，本技术的示意性实施例及其说明用于解释本技术，并不构成对本技术的不当限定。在附图中：
[0050]
图1为本说明书实施例提供的一种疲劳短裂纹扩展速率预测模型的确定方法的流程示意图；
[0051]
图2为实施例中所述疲劳短裂纹扩展速率
‑
裂纹长度趋势图；
[0052]
图3为实施例中材料不同应力范围下实验数据及拟合曲线；
[0053]
图4为实施例中eh36钢c、m与应力范围之间的线性拟合图。
[0054]
图5为本说明书实施例提供的对应于图1的一种疲劳短裂纹扩展速率预测模型的确定装置的结构示意图。
具体实施方式
[0055]
为使本技术的目的、技术方案和优点更加清楚，下面将结合本技术具体实施例及相应的附图对本技术技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0056]
以下结合附图，详细说明本技术各实施例提供的技术方案。
[0057]
图1为本说明书实施例提供的一种方法的流程示意图。从程序角度而言，流程的执行主体可以为搭载于应用服务器的程序或应用客户端。
[0058]
如图1所示，该流程可以包括以下步骤：
[0059]
步骤110：对材料进行不同应力范围条件下疲劳短裂纹扩展实验，记录实验数据；
[0060]
步骤120：根据实验数据计算有效应力强度因子δk
eff
和疲劳短裂纹扩展速率da/dn；
[0061]
步骤130：计算微观结构短裂纹长度对应的有效应力强度因子δk
eff
‑
msc
；
[0062]
步骤140：基于δk
eff
和δk
eff
‑
msc
确定第一预测疲劳断裂纹速度模型；
[0063]
步骤150：确定第一预测疲劳断裂纹速度模型中的拟合常数，带入所述第一预测疲劳断裂纹速度模型中。
[0064]
具体的，一种疲劳短裂纹扩展速率预测模型，可以包括以下步骤：
[0065]
s1、通过公式(1)和(2)分别计算不同应力范围(δσ)下的应力强度因子范围δk和有效应力强度因子范围δk
eff
，疲劳短裂纹扩展速率da/dn通过公式(3)计算：
[0066][0067]
其中δk为应力强度因子范围，y为形状因子需按照实际试样所对应标准进行计算，δσ为应力范围，a为表面裂纹长度；
[0068][0069]
其中δk
eff
为有效应力强度因子范围，σ
op
为裂纹张开应力，σ
max
为最大应力，r为应力比；
[0070][0071]
其中δn为循环次数增加数值，δa为表面裂纹增加长度，a
i 1
为循环次数为n
i 1
时表面裂纹长度；
[0072]
s2、通过公式(4)和(5)计算微观结构短裂纹长度a
msc
所对应的有效应力强度因子范围δk
eff
‑
msc
；
[0073][0074]
其中微观结构短裂纹长度a
msc
的长度是通过da/dn
‑
a趋势图定义的，a
msc
长度为裂纹扩展速率递减与增加的临界点，δk
msc
为a
msc
对应的应力强度因子范围，y为形状因子需按照实际试样所对应标准进行计算，δσ为应力范围；
[0075][0076]
其中δk
eff
‑
msc
为微观结构短裂纹长度a
msc
对应的有效应力强度因子范围，σ
op
为裂纹张开应力，σ
max
为最大应力，r为应力比；
[0077]
s3、基于δk
eff
和δk
eff
‑
msc
，提出公式(6)计算不同δσ下疲劳短裂纹扩展速率：
[0078]
da/dn＝c|δk
eff
‑
δk
eff
‑
msc
|
m
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0079]
其中c和m为拟合常数，其大小随δσ数值大小增加而增加；
[0080]
s4、将获得的不同δσ下da/dn
‑
δk
eff
数据在双对数坐标系中绘图，并对其按照公式(6)进行非线性曲线拟合，获得不同δσ下对应的拟合常数c、m，后分别在双对数坐标系中绘制c
‑
δσ和m
‑
δσ并对其线性拟合，得到公式(7)和(8)：
[0081]
c＝10
a
(δσ)
b
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0082]
m＝10
x
(δσ)
y
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0083]
其中a，b和x，y分别为拟合常数
[0084]
s5、引入应力归一化系数β，计算材料流变应力σ0，利用公式(7)和(8)通过线性外推法分别获得σ0条件下所对应的拟合常数和最终确定公式(9)描述材料不同应力范围条件下疲劳短裂纹v形扩展速率：
[0085][0086]
其中应力归一化系数β为δσ/σ0；σ0为流变应力，取材料屈服强度与抗拉强度之和的平均值。
[0087]
经由上述的技术方案可知，与现有技术相比，本发明基于δk
eff
及δk
eff
‑
msc
及paris公式提出分段函数精准描述某种材料疲劳短裂纹v形扩展趋势，整个测试过程精度高、效率高，应用范围广，具有广阔应用前景和经济效益，海洋工程材料性能评估提供重要理论依据及实践基础。
[0088]
基于图1的方法，本说明书实施例还提供了该方法的一些具体实施方式，下面进行说明。
[0089]
实施例1
[0090]
本实施例公开了一种疲劳短裂纹扩展速率预测模型，包括以下步骤：
[0091]
s1、对eh36钢进行δσ＝160，192和106.67mpa三种不同应力范围条件下疲劳短裂纹扩展实验；
[0092]
其中疲劳短裂纹长度为0
‑
1.5mm；
[0093]
通过公式(1)计算单边缺口试样的应力强度因子δk，后通过公式(2)消除裂纹闭合效应并获得有效应力强度因子范围δk
eff
，疲劳短裂纹扩展速率da/dn通过公式(3)计算：
[0094][0095]
其中应力强度因子δk值计算方法需按单边缺口试样所计算，δσ为应力范围，a为表面裂纹长度，l为缺口深度，w为试样宽度；
[0096][0097]
其中δk
eff
为有效应力强度因子范围，σ
op
为裂纹张开应力，σ
max
为最大应力，r为应力比；
[0098][0099]
其中δn为循环范围，δa为表面裂纹增加长度，a
i 1
为循环次数为n
i 1
时表面裂纹长度；
[0100]
s2、获得如图2所示的疲劳短裂纹扩展速率
‑
裂纹长度趋势图并寻找材料裂纹扩展速率递减与增加的临界点，即疲劳微观结构短裂纹长度a
msc＝
0.2mm，后通过公式(4)和(5)计算a
msc
所对应的δk
eff
‑
msc
，提出基于δk
eff
和δk
eff
‑
msc
的分段模型公式(6)计算疲劳短裂纹速率：
[0101][0102]
其中δk
msc
为微观结构短裂纹长度对应的应力强度因子范围，δk
msc
按照单边缺口
sent试样相关标准进行计算，为所述sent试样相关标准中系数方程，δσ为应力范围，a
msc
为微观结构短裂纹长度，l为缺口深度，w为试样宽度，；
[0103][0104]
其中δk
eff
‑
msc
为微观结构短裂纹长度对应的有效应力强度因子范围，σ
op
为裂纹张开应力，σ
max
为最大应力，r为应力比；
[0105]
da/dn＝c|δk
eff
‑
δk
eff
‑
msc
|
m
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0106]
其中c和m为拟合常数，其大小随应力范围数值大小增加而增加；
[0107]
s3、通过公式(4)和(5)计算δσ＝160，192和106.67mpa三种应力范围条件下所对应的δk
eff
‑
msc
分别为14.07，17.63和11.59mpa*m
0.5
，通过公式(6)获得所述三种应力范围条件下基于δk
eff
和δk
eff
‑
msc
的分段模型及其对应的拟合常数c和m，如图3及以下公式所示：
[0108]
da/dn＝9.65*10
‑9|δk
eff
‑
14.07|
0.388
(δσ＝160mpa)
[0109]
da/dn＝1.16*10
‑8|δk
eff
‑
17.63|
0.467
(δσ＝192mpa)
[0110]
da/dn＝6.24*10
‑9|δk
eff
‑
11.59|
0.259
(δσ＝106.67mpa)
[0111]
后对不同应力范围下拟合常数c、m与应力范围大小进行线性拟合，如图4所示，得到公式(7)和(8)并通过线性外推获得eh36钢流变应力σ0＝460mpa条件下所对应和并引入应力归一化系数β，确定公式(9)描述eh36钢不同应力范围条件下疲劳短裂纹v形扩展速率：
[0112]
c＝10
‑
10.373
(δσ)
1.058
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0113]
m＝10
‑
2.621
(δσ)
1.003
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0114][0115]
其中公式(9)应力归一化系数β为δσ/σ0；σ0为流变应力并取屈服强度与抗拉强度之和的平均值。
[0116]
基于同样的思路，本说明书实施例还提供了上述方法对应的装置。图5为本说明书实施例提供的对应于图1的一种疲劳短裂纹扩展速率预测模型的确定装置的结构示意图。如图5所示，该装置可以包括：
[0117]
应力实验模块510，用于对材料进行不同应力范围条件下疲劳短裂纹扩展实验，记录实验数据；
[0118]
第一计算模块520，用于根据实验数据计算有效应力强度因子δk
eff
和疲劳短裂纹扩展速率da/dn；
[0119]
第二计算模块530，用于计算微观结构短裂纹长度对应的有效应力强度因子δk
eff
‑
msc
；
[0120]
模型确定模块540，用于基于δk
eff
和δk
eff
‑
msc
确定第一预测疲劳断裂纹速度模型；
[0121]
常数确定模块550，用于确定第一预测疲劳断裂纹速度模型中的拟合常数，带入所述第一预测疲劳断裂纹速度模型中。
[0122]
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的
包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
[0123]
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
[0124]
以上所述仅为本技术的实施例而已，并不用于限制本技术。对于本领域技术人员来说，本技术可以有各种更改和变化。凡在本技术的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本技术的权利要求范围之内。