一种高韧性材料J-R阻力曲线的测试方法

摘要

本发明涉及一种高韧性材料J-R阻力曲线的测试方法，它包括以下步骤：试样加工、三点弯曲试验、裂纹长度的测定、每一载荷下的J积分计算、差值函数的建立、裂纹长度与p-V数据对的递推计算以及J-R阻力曲线的构建。本发明不采用引伸计测试裂纹嘴张开位移及卸载柔度，而是通过试验机直接采集三点弯曲试样的变形和断裂全过程的载荷-位移曲线，而后通过一系列的计算得到材料的J-R阻力曲线，从而有效地克服了高韧性材料测试中试样大变形所带来的测试仪器的限制，可测得长尺寸裂纹扩展下的J-R阻力曲线，并可用于不同温度和不同尺寸试样的J-R阻力曲线测量。

权利要求

1.一种高韧性材料J-R阻力曲线的测试方法，其特征在于，所述测试方法包括以下步骤：步骤S1，将待测试的高韧性材料加工成三点弯曲试样，对该试样加工一缺口，并预制疲劳裂纹；

步骤S2，测得在对所述试样加载和裂纹稳定扩展过程直至最后断裂的载荷p与加载点位移Ｖ的关系曲线，即p-V曲线，并测得最终断裂或停机载荷pf；

步骤S3，通过九分法计算得到所述试样断裂后断口的初始裂纹长度a0和最终裂纹扩展长度af；

步骤S4，根据初始裂纹长度a0及所述p-V曲线测得若干个载荷点pi，通过J积分基本测试法计算得到每个载荷pi下的修正应力强度因子Ki和修正J积分Ji；

步骤S5，根据初始裂纹长度a0及修正J积分Ji得到修正的裂纹长度ab(i)；将所述p-V曲线上的最终断裂或停机载荷pf的p-V数据对（pf，Vf）标称化为锚点数组（PN(f)，Vpl(f)），将最大载荷pmax之前，且不包括最大载荷pmax的所有p-V数据对（pi，Vi）标称化为标称p-V数组（PN(i)，Vpl(i)）；将标称位移Vpl(i)>0.001的标称p-V数组（PN(i)，Vpl(i)）及锚点数组（PN(f)，Vpl(f)）进行曲线拟合后建立差值函数步骤S6，根据所述初始裂纹长度a0和最终裂纹扩展长度af，将裂纹长度划分为若干个不同的长度段a0、a0+Δ、a0+2Δ、a0+3Δ……af，并对任一长度段标记为 其中，取将初始裂纹长度 代入所述差值函数 在标称位移Vpl(i)>0.002处将p-V数据对（pi，Vi）依次代入所述差值函数 当含有裂纹长度 及p-V数据对（pi，Vi）的差值函数变化率小于等于1%时，则判定此时p-V数据对（pi，Vi）即是对应初始裂纹长度 的p-V数据对同理，将裂纹长度 及p-V数据对 之后的p-V数据对（pi，Vi）依次代入所述差值函数 当差值函数变化率小于等于1%，得到对应裂纹长度 的p-V数据对依次类推，得到与裂纹长度 相对应的p-V数据对

其中，要求p-V数据对 为最终p-V数据对

（pf,Vf）；当裂纹长度 与最终p-V数据对（pf，Vf）不满足差值函数变化率小于等于1%时，调整裂纹长度 之前的与裂纹长度 对应的p-V数对 并要求裂纹长度 之前的裂纹长度与相应p-V数据对在满足差值函数变化率小于等于1%的前提下，使得最终p-V数对（pf,Vf）与裂纹长度 差值函数变化率小于等于1%；

步骤S7，通过公式 计算得到裂纹扩展量 通过阻力曲线测试法中的J

积分计算公式计算得到每个载荷pi下对应裂纹扩展量 的J积分 并以裂纹扩展量为横坐标，J积分 为纵坐标，将所有的数据对 描在坐标图上，并光滑连接在一起，即得到材料J-R阻力曲线。

2.根据权利要求1所述的高韧性材料J-R阻力曲线的测试方法，其特征在于，所述步骤S1还包括定义所述试样高度为W，试样宽度为B，试样跨距S=4W，所述缺口与疲劳裂纹的总深度，即初始裂纹长度为a0。

3.根据权利要求1所述的高韧性材料J-R阻力曲线的测试方法，其特征在于，所述步骤S2中的p-V曲线通过对所述试样进行三点弯曲加载试验测得。

4.根据权利要求1所述的高韧性材料J-R阻力曲线的测试方法，其特征在于，所述步骤S3包括测量所述试样断裂后断口的裂纹初始长度a01、a02、a03……a09，以及最终稳定扩展长度af1、af2、af3………af9，并按式（1）和式（2）分别计算得到初始裂纹长度a0和最终裂纹扩展长度af，

5.根据权利要求1所述的高韧性材料J-R阻力曲线的测试方法，其特征在于，所述步骤S4包括：首先，通过式（3）和式（4）计算得到修正应力强度因子Ki，此时ai取初始裂纹长度a0，然后，通过式（5）、式（6）和式（7）计算得到修正J积分Ji，在式（3）至式（7）中，取ai=a0，S为试样跨距，W为试样高度，B为试样宽度，BN为试样有效厚度，取b0=（W-a0）为试样初始韧带长度，E为材料的弹性模量，v为材料的泊松比，Atot(i)为对应载荷pi时p-V曲线下包围的总面积，tgα为p-V曲线初始弹性段对应的直线斜率，Apl(i)为对应载荷pi时p-V曲线包围的塑性部分面积，Jpl(i)为对应载荷pi时p-V曲线塑性部分的Ｊ积分，η为参数，该参数η通过取ai=a0代入式（17）中计算得到，

6.根据权利要求5所述的高韧性材料J-R阻力曲线的测试方法，其特征在于，所述步骤S5包括：首先，通过式（8）计算得到修正的裂纹长度ab(i)，

其中，Ji为修正J积分，σf为流动应力，其值取为σf=(σy+σu)/2，σy和σu分别为材料的屈服强度和抗拉强度；

其次，对最大载荷pmax之前，且不包括最大载荷pmax的载荷pi通过式（9）计算得到标称载荷PN(i)，其中，W为试样高度，B为试样宽度，η为参数，该参数η通过取ai=ab(i)代入式（17）中计算得到，再次，对最大载荷pmax之前，且不包括最大载荷pmax的载荷pi对应的位移Vi按照式（10）、式（11）和式（12）计算得到标称位移Vpl(i)，在式（10）至式（12）中，取ai=ab(i)，BN为试样有效厚度，E为材料的弹性模量，Ci为对应载荷pi时的卸载柔度；

同理，通过将最终裂纹扩展长度af替代式（9）中的ab(i)，计算得到最终载荷PN(f)；通过将最终裂纹扩展长度af替代式（10）中的ai，并根据式（11）和式（12）计算得到最终位移Vpl(f)；

然后，将标称p-V数组（PN(i)，Vpl(i)）及锚点数组（PN(f)，Vpl(f)）按式（13）进行曲线拟合，其中，c1、c2、c3、c4为相应的拟合系数；

最后，建立由式（14）表示的差值函数

7.根据权利要求6所述的高韧性材料J-R阻力曲线的测试方法，其特征在于，所述步骤S7中的J积分 通过将式（3）、式（4）、式（10）、式（11）、式（12）以及式（15）至式（18）代入式（5）中计算得到，其中：取bi-1=W-ai-1，对应 的Jpl(i)取为零；

8.根据权利要求7所述的高韧性材料J-R阻力曲线的测试方法，其特征在于，当所述试样加工有侧槽时，则BN为去除侧槽厚度的试样有效厚度，否则，BN=B。

说明书

一种高韧性材料J-R阻力曲线的测试方法

技术领域

[0001] 本发明涉及材料力学性能测试技术，尤其涉及一种高韧性材料J-R阻力曲线的测试方法。

背景技术

[0002] 在压力容器和管道的安全可靠性设计和评价中，目前国内外的研究发展趋势是按先漏后爆(LBB：leak before break)的准则进行设计和评价。LBB分析评价的目的是确证：穿过裂纹的流体泄漏在裂纹达到快速失稳扩展的极限长度前可以被检测到。因此，在LBB分析评价中，主要是用弹塑性断裂力学的方法计算裂纹扩展的稳定性，为此需要准确测定材料在较长延性裂纹扩展时的阻力曲线，即J-R阻力曲线。另一方面，在含缺陷结构的完整性评价中，为进行延性撕裂评定，也需要测定材料的J-R阻力曲线。

[0003] 传统测试材料J-R阻力曲线的方法主要有多试样法和单试样法。多试样法需要测试6个以上的试样，不仅浪费材料、时间及人力，而且由于取材位置的不同，通常得到的数据分散性较大。单试样法能较好的避免多试样法的缺陷，但是单试样法过于依赖测试仪器及操作技巧。目前电位降法和弹性卸载柔度法是具有代表性的两种单试样测试方法，裂纹扩展长度分别通过电流信号及试样柔度计算得到。电位降法需要准确判断裂纹起裂点，对仪器及操作技巧要求较高，往往难以使用。弹性卸载柔度法因其测试计算理论成熟，操作相对简单，已形成测试标准，并被广泛应用。

[0004] 然而，随着工程结构材料的不断改进和其韧性的提高，在用弹性卸载柔度法测试J-R阻力曲线时，发现由于该类材料的高韧性和高延性，裂尖严重钝化，并会在发生大范围屈服后才发生延性裂纹起裂。此时，由于试样宏观变形很大，受制于试验机引伸计量程和计算控制软件的限制，而不能测到裂纹扩展过程中完整的卸载柔度载荷-位移曲线，仅可以测试计算得到1mm左右的少量裂纹扩展时的J-R曲线，从而无法满足LBB分析评价的要求。此外，由于如核电部件等通常工作在较高温度下的材料，受测试仪器的限制，一般也难以用卸载柔度法测得材料在高温下的J-R阻力曲线。而且，在考虑地震动态载荷时的工程结构完整性评价中，需要测试评价材料在动态载荷下的J-R阻力曲线，因此，适用于准静态加载的卸载柔度法也难于使用。另外，工程结构件的完整性评定和LBB分析技术的研究发展趋势是考虑材料与结构的断裂拘束效应，以实现结构安全可靠性的准确评定和预测，为此需要测试不同尺寸和裂纹深度试样(即不同拘束)的J-R阻力曲线。而目前常用的单试样卸载柔度法主要适用于大尺寸深裂纹的高拘束试样的J-R阻力曲线。

[0005] 综上所述，针对目前常用的J-R阻力曲线测试方法应用于高韧性材料的局限性，迫切需要研发一种适用于该类材料J-R阻力曲线测试的新方法。

发明内容

[0006] 为了解决上述现有技术存在的问题，本发明旨在提供一种简单、准确且应用范围较宽的高韧性材料J-R阻力曲线的测试方法，以用于测试高韧性材料在不同温度和不同尺寸试样下的J-R阻力曲线，从而得到该类材料长尺寸裂纹扩展下的J-R阻力曲线，为工程结构的安全可靠性设计和评价提供准确的基础材料参数。

[0007] 本发明所述的一种高韧性材料J-R阻力曲线的测试方法，它包括以下步骤：

[0008] 步骤S1，将待测试的高韧性材料加工成三点弯曲试样，对该试样加工一缺口，并预制疲劳裂纹；

[0009] 步骤S2，测得在对所述试样加载和裂纹稳定扩展过程直至最后断裂的载荷p与加载点位移V的关系曲线，即p-V曲线，并测得最终断裂或停机载荷pf；

[0010] 步骤S3，通过九分法计算得到所述试样断裂后断口的初始裂纹长度a0和最终裂纹扩展长度af；

[0011] 步骤S4，根据初始裂纹长度a0及所述p-V曲线测得若干个载荷点pi，通过J积分基本测试法计算得到每个载荷pi下的修正应力强度因子Ki和修正J积分Ji；

[0012] 步骤S5，根据初始裂纹长度a0及修正J积分Ji得到修正的裂纹长度ab(i)；将所述p-V曲线上的最终断裂或停机载荷pf的p-V数据对(pf，Vf)标称化为锚点数组(PN(f)，Vpl(f))，将最大载荷pmax之前，且不包括最大载荷pmax的所有p-V数据对(pi，Vi)标称化为标称p-V数组(PN(i)，Vpl(i))；将标称位移Vpl(i)＞0.001的标称p-V数组(PN(i)，Vpl(i))及锚点数组(PN(f)，Vpl(f))进行曲线拟合后建立差值函数

[0013] 步骤S6，根据所述初始裂纹长度a0和最终裂纹扩展长度af，将裂纹长度划分为若干个不同的长度段a0、a0+Δ、a0+2Δ、a0+3Δ......af，并对任一长度段标记为 其中，取[0014] 将初始裂纹长度 代入所述差值函数 在标称位移Vpl(i)＞0.002处将p-V数据对(pi，Vi)依次代入所述差值函数 当含有裂纹长度 及p-V数据对(pi，Vi)的差值函数变化率小于等于1％时，则判定此时p-V数据对(pi，Vi)即是对应初始裂纹长度 的p-V数据对

[0015] 同理，将裂纹长度 及p-V数据对 之后的p-V数据对(pi，Vi)依次代入所述差值函数 当差值函数变化率小于等于1％，得到对应裂纹长度 的p-V数据对[0016] 依次类推，得到与裂纹长度 相对应的p-V数据对
其中，要求p-V数据对 为最终p-V数据对(pf，
Vf)；当裂纹长度 与最终p-V数据对(pf，Vf)不满足差值函数变化率小于等于1％时，调整裂纹长度 之前的与裂纹长度 对应的p-V数对 并要求裂纹长度 之前的裂
纹长度与相应p-V数据对在满足差值函数变化率小于等于1％的前提下，使得最终p-V数对(pf，Vf)与裂纹长度 差值函数变化率小于等于1％；

[0017] 步骤S7，通过公式 计算得到裂纹扩展量 通过阻力曲线测试法中的J积分计算公式计算得到每个载荷pi下对应裂纹扩展量 的J积分 并以裂纹扩展量为横坐标，J积分 为纵坐标，将所有的数据对 描在坐标图上，并光滑连接在
一起，即得到材料J-R阻力曲线。

[0018] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述步骤S1还包括定义所述试样高度为W，试样宽度为B，试样跨距S＝4W，所述缺口与疲劳裂纹的总深度，即初始裂纹长度为a0。

[0019] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述步骤S2中的p-V曲线通过对所述试样进行三点弯曲加载试验测得。

[0020] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述步骤S3包括测量所述试样断裂后断口的裂纹初始长度a01、a02、a03……a09，以及最终稳定扩展长度af1、af2、af3………af9，并按式(1)和式(2)分别计算得到初始裂纹长度a0和最终裂纹扩展长度af，[0021]

[0022]

[0023] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述步骤S4包括：

[0024] 首先，通过式(3)和式(4)计算得到修正应力强度因子Ki，此时ai取初始裂纹长度a0，

[0025]

[0026]

[0027] 然后，通过式(5)、式(6)和式(7)计算得到修正J积分Ji，

[0028]

[0029]

[0030]

[0031] 在式(3)至式(7)中，取ai＝a0，S为试样跨距，W为试样高度，B为试样宽度，BN为试样有效厚度，取b0＝(W-a0)为试样初始韧带长度，E为材料的弹性模量，v为材料的泊松比，Atot(i)为对应载荷pi时p-V曲线下包围的总面积，tgα为p-V曲线初始弹性段对应的直线斜率，Apl(i)为对应载荷pi时p-V曲线包围的塑性部分面积，Jpl(i)为对应载荷pi时p-V曲线塑性部分的J积分，η为参数。

[0032] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述步骤S5包括：

[0033] 首先，通过式(8)计算得到修正的裂纹长度ab(i)，

[0034]

[0035] 其中，Ji为修正J积分，σf为流动应力，其值取为σf＝(σy+σu)/2，σy和σu分别为材料的屈服强度和抗拉强度；

[0036] 其次，对最大载荷pmax之前，且不包括最大载荷pmax的载荷pi通过式(9)计算得到标称载荷PN(i)，

[0037]

[0038] 其中，W为试样高度，B为试样宽度，η为参数；

[0039] 再次，对最大载荷pmax之前，且不包括最大载荷pmax的载荷pi对应的位移Vi按照式(10)、式(11)和式(12)计算得到标称位移Vpl(i)，

[0040]

[0041]

[0042]

[0043] 在式(10)至式(12)中，取ai＝ab(i)，BN为试样有效厚度，E为材料的弹性模量，Ci为对应载荷pi时的卸载柔度；

[0044] 同理，通过将最终裂纹扩展长度af替代式(9)中的ab(i)，计算得到最终载荷PN(f)；通过将最终裂纹扩展长度af替代式(10)中的ai，并根据式(11)和式(12)计算得到最终位移Vpl(f)；

[0045] 然后，将标称p-V数组(PN(i)，Vpl(i))及锚点数组(PN(f)，Vpl(f))按式(13)进行曲线拟合，

[0046]

[0047] 其中，c1、c2、c3、c4为相应的拟合系数；

[0048] 最后，建立由式(14)表示的差值函数

[0049]

[0050] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述步骤S7中的J积分 通过将式(3)、式(4)、式(10)、式(11)、式(12)以及式(15)至式(18)代入式(5)中计算得到，[0051]

[0052] 其中：取bi-1＝W-ai-1，对应 的Jpl(i)取为零；

[0053]

[0054]

[0055]

[0056] 在上述的高韧性材料J-R阻力曲线的测试方法中，当所述试样加工有侧槽时，则BN为去除侧槽厚度的试样有效厚度，否则，BN＝B。

[0057] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述式(6)中的参数η通过取ai＝a0代入式(17)中计算得到。

[0058] 在上述的高韧性材料J-R阻力曲线的测试方法中，所述式(9)中的参数η通过取ai＝ab(i)代入式(17)中计算得到。

[0059] 由于采用了上述的技术解决方案，本发明不采用引伸计测试裂纹嘴张开位移及卸载柔度，而是通过试验机直接采集三点弯曲试样的变形和断裂全过程的载荷-位移曲线，而后通过一系列的计算得到材料的J-R阻力曲线，从而有效地克服了高韧性材料测试中试样大变形所带来的测试仪器的限制，可测得长尺寸裂纹扩展下的J-R阻力曲线，并可用于不同温度和不同尺寸试样的J-R阻力曲线测量。

附图说明

[0060] 图1是本发明一种高韧性材料J-R阻力曲线的测试方法中三点弯曲试样的结构示意图；

[0061] 图2是本发明中p-V曲线及p-V数据对(pi，Vi)的示意图；

[0062] 图3是本发明实施例中测量的常温及高温下A508及316L的p-V曲线；

[0063] 图4(a)是本发明实施例中计算的常温下A508三点弯曲试样的标称p-V曲线及对应的拟合曲线的示意图；

[0064] 图4(b)是本发明实施例中计算的常温下316L三点弯曲试样的标称p-V曲线及对应的拟合曲线的示意图；

[0065] 图5(a)是本发明实施例中计算的高温340℃下A508三点弯曲试样的标称p-V曲线及对应的拟合曲线的示意图；

[0066] 图5(b)是本发明实施例中计算的高温340℃下316L三点弯曲试样的标称p-V曲线及对应的拟合曲线的示意图；

[0067] 图6(a)是本发明实施例中测得的常温下A508三点弯曲试样的J-R阻力曲线及其与卸载柔度法测量结果的对比示意图；

[0068] 图6(b)是本发明实施例中测得的常温下316L三点弯曲试样的J-R阻力曲线及其与卸载柔度法测量结果的对比示意图；

[0069] 图7(a)是本发明实施例中测得的高温340℃下A508三点弯曲试样的J-R阻力曲线；

[0070] 图7(b)是本发明实施例中测得的高温340℃下316L三点弯曲试样的J-R阻力曲线。

具体实施方式

[0071] 下面结合附图，给出本发明的较佳实施例，并予以详细描述。

[0072] 本发明，即一种高韧性材料J-R阻力曲线的测试方法，包括以下步骤：

[0073] 步骤S1，将待测试的高韧性材料加工成三点弯曲试样，对该试样加工一缺口，并预制疲劳裂纹。

[0074] 具体来说：根据测试需要，参考GB/T21143-2007标准，对待测试的高韧性材料，加工成图1所示的三点弯曲试样，并通过线切割方式加工窄的缺口，预制疲劳裂纹；定义试样高度为W，试样宽度为B，试样跨距S＝4W，缺口与疲劳裂纹的总深度，即初始裂纹长度为a0。

[0075] 步骤S2，测得在对试样加载后的变形和裂纹稳定扩展直至最后断裂的载荷p与加载点位移V的关系曲线，即p-V曲线，并测得最终断裂或停机载荷pf。

[0076] 具体来说：对图1所示的试样进行三点弯曲加载试验，即施加力F，测试记录得到如图2所示的完整的p-V曲线。

[0077] 步骤S3，通过九分法计算得到试样断裂后断口的初始裂纹长度a0和最终裂纹扩展长度af。

[0078] 具体来说：按GB/T21143-2007标准，在断裂后试样的断口上，测量断口的裂纹初始长度a01、a02、a03……a09，以及最终稳定扩展长度af1、af2、af3………af9，并按式(1)和式(2)分别计算得到初始裂纹长度a0和最终裂纹扩展长度af，

[0079]

[0080]

[0081] 步骤S4，根据初始裂纹长度a0及p-V曲线测得若干个载荷点pi，通过J积分基本测试法计算得到每个载荷pi下的修正应力强度因子Ki和修正J积分Ji。

[0082] 具体来说：首先，在图2所示的p-V曲线上，测得一系列载荷点pi；

[0083] 然后，通过式(3)和式(4)计算得到每个载荷pi下的修正应力强度因子Ki，[0084]

[0085]

[0086] 最后，通过式(5)、式(6)和式(7)计算得到每个载荷pi下的修正J积分Ji，[0087]

[0088]

[0089]

[0090] 在式(3)至式(7)中，取ai＝a0，S为试样跨距，W为试样高度，B为试样宽度，BN为试样有效厚度(如图1所示，当试样加工有侧槽时，则BN为去除侧槽厚度的试样有效厚度，否则，BN＝B)，取b0＝(W-a0)为试样初始韧带长度，E为材料的弹性模量，v为材料的泊松比，Atot(i)为对应载荷pi时p-V曲线下包围的总面积(即图2中斜线区面积)，tgα为p-V曲线初始弹性段对应的直线斜率(如图2所示)，Apl(i)为对应载荷pi时载荷位移曲线包围的塑性部分面积；Jpl(i)是对应载荷pi时p-V曲线塑性部分的J积分，η为参数(该参数η通过取ai＝a0代入下面的式(17)中计算得到)。

[0091] 步骤S5，根据初始裂纹长度a0及修正J积分Ji得到修正的裂纹长度ab(i)；将p-V曲线上的最终断裂或停机载荷pf的p-V数据对(pf，Vf)标称化为锚点数组(PN(f)，Vpl(f))，将最大载荷pmax之前(不包括最大载荷pmax)的所有p-V数据对(pi，Vi)标称化为标称p-V数组(PN(i)，Vpl(i))；将标称位移Vpl(i)＞0.001的标称p-V数组(PN(i)，Vpl(i))及锚点数组(PN(f)，Vpl(f))进行曲线拟合后建立差值函数

[0092] 具体来说：首先，通过式(8)计算得到修正的裂纹长度ab(i)，

[0093]

[0094] 其中，Ji为修正J积分(由式(3)至(5)可得)，σf为流动应力，其值取为σf＝(σy+σu)/2，σy和σu分别为材料的屈服强度和抗拉强度；

[0095] 其次，对最大载荷pmax之前的载荷pi(不包括最大载荷pmax)通过式(9)计算得到标称载荷PN(i)，

[0096]

[0097] 其中，W为试样高度，B为试样宽度，参数η通过取ai＝ab(i)代入下面的式(17)中计算得到(不同深度的裂纹计算得到的参数η不同)；

[0098] 再次，对最大载荷pmax之前的载荷pi(不包括最大载荷pmax)对应的位移Vi按照式(10)、式(11)和式(12)计算得到标称位移Vpl(i)，

[0099]

[0100]

[0101]

[0102] 在式(10)至式(12)中，取ai＝ab(i)，BN为试样有效厚度，E为材料的弹性模量，Ci为对应载荷pi时的卸载柔度；对加工有侧槽的试样，式(10)中的Be由式(11)计算；对未加工侧槽的试样，则取Be＝B；

[0103] 同理，通过将最终裂纹扩展长度af替代式(9)中ab(i)，计算得到最终载荷PN(f)；通过将最终裂纹扩展长度af替代式(10)中的ai，得到相应的卸载柔度Cf，并根据式(11)和式(12)，计算得到最终位移Vpl(f)；

[0104] 然后，将标称p-V数组(PN(i)，Vpl(i))及锚点数组(PN(f)，Vpl(f))按式(13)进行曲线拟合，

[0105]

[0106] 其中，c1、c2、c3、c4为相应的拟合系数；

[0107] 最后，建立由式(14)表示的差值函数

[0108]

[0109] 步骤S6，根据初始裂纹长度a0和最终裂纹扩展长度af，将裂纹长度划分为若干个不同的长度段a0、a0+Δ、a0+2Δ、a0+3Δ......af，并对任一长度段标记为 其中，取[0110] 将初始裂纹长度 代入差值函数 在标称位移Vpl(i)＞0.002处将p-V数据对(pi，Vi)依次代入差值函数 当含有裂纹长度 及p-V数据对(pi，Vi)的差值函数变化率小于等于1％时，则判定此时p-V数据对(pi，Vi)即是对应初始裂纹长度 的p-V数据对[0111] 同理，将裂纹长度 及p-V数据对 之后的p-V数据对(pi，Vi)依次代入所述差值函数 当差值函数变化率小于等于1％，得到对应裂纹长度 的p-V数据对[0112] 依次类推，得到与裂纹长度 相对应的p-V数据对
其中，要求p-V数据对 为最终p-V数据对(pf，
Vf)；当裂纹长度 与最终p-V数据对(pf，Vf)不满足差值函数变化率小于等于1％时，调整裂纹长度 之前的与裂纹长度 对应的p-V数对 并要求裂纹长度 之前的裂
纹长度与相应p-V数据对在满足差值函数变化率小于等于1％的前提下，使得最终p-V数对(pf，Vf)与裂纹长度 差值函数变化率小于等于1％。

[0113] 具体来说：在初始裂纹长度a0和最终裂纹扩展长度af之间，将裂纹长度划分为若干个不同的长度段，划分的数量依所测J-R阻力曲线的精度而定，划分的长度段数量越多，测量计算的J-R阻力曲线精度越高，但计算工作量则较大。如划分长度分别为a0，a0+Δ，a0+2Δ，a0+3Δ......af，则对任一长度可记为 之后将裂纹长度 代入式(10)计算得到卸载柔度 将初始裂纹长度 代入式(14)，在起点标称位移Vpl(i)＞0.002处将p-V数据对(pi，Vi)依次代入式(14)，使得满足含有裂纹长度 及p-V数据对(pi，Vi)的差值函数变化率小于等于1％时(即对应载荷-位移点的 )，认定此时p-V数据对
(pi，Vi)即是对应初始裂纹长度 的p-V数据对 之后将 及 之后的
p-V数据对(pi，Vi)代入式(14)，并判断差值函数变化率小于等于1％，得到对应裂纹长度的p-V数据对 依次类推，可以分别得到裂纹长度 相对应的p-V
数据对

[0114] 步骤S7，通过公式 即式(19)，计算得到裂纹扩展量 通过阻力曲线测试法中的J积分计算公式计算得到每个载荷pi下对应裂纹扩展量 的J积分 并以裂纹扩展量 为横坐标，J积分 为纵坐标，将所有的数据对 描在坐标图上，并
光滑连接在一起，即得到材料的J-R阻力曲线。

[0115] 具体来说：J积分 通过将式(3)、式(4)、式(10)、式(11)、式(12)以及式(15)至式(18)代入式(5)中计算得到，

[0116]

[0117] 其中：b(i-1)＝W-a(i-1)，对应 的Jpl(i)取为零；

[0118]

[0119]

[0120]

[0121] 实施例

[0122] 某核电设备设计和制造公司需要测试高韧性压水堆核反应压力容器铁素体钢A508和一回路主管道奥氏体不锈钢316L材料在室温及340℃高温下的J-R阻力曲线。通过圆棒拉伸试验测得的A508及316L在室温及340℃高温下的拉伸力学性能数据见表1：

[0123] 表1材料的拉伸性能数据

[0124]

[0125] 下面按本发明的方法测试上述4个试样的J-R阻力曲线。

[0126] (1)试样加工

[0127] 将A508和316L两种材料加工成如图1所示的三点弯曲试样；试样高度W＝32mm，试样厚度B＝16mm，三点弯加载跨距S＝4W＝128mm，切割缺口和预制疲劳裂纹总深度为16mm左右；两种材料各加工两个试样。

[0128] (2)三点弯曲试验

[0129] 对两种材料的试样分别在常温和340℃的核电设计操作温度下进行三点弯曲试验；试验在Instron万能试验机上进行，以0.5mm/min的速度加载，测试记录整个加载过程中的载荷p与加载点位移V的关系曲线，即p-V曲线。4个试样所测得的p-V曲线如图3所示。

[0130] (3)裂纹长度的测定

[0131] 试验完成后，对试样进行发蓝处理，而后锯或压断。按相关标准的方法，在断裂后试样的断口上，采用九分法在工具显微镜下测量裂纹初始长度a01、a02、a03......a09，以及最终稳定扩展长度af1、af2、af3.........af9，并按式(1)和式(2)计算得到初始裂纹长度a0及最终裂纹扩展长度af，4个试样的初始及最终裂纹长度列于表2：

[0132] 表2裂纹长度测定结果

[0133]

[0134] (4)每一载荷下的J积分计算

[0135] 提取图3中4个试样的p-V数据，用表2中测量的初始裂纹长度a0替代式(3)和式(4)中的ai，得到每个载荷pi下的修正应力强度因子Ki；而后用式(5)至式(7)计算得到每个载荷pi下的修正Ji；其中，由于4个试样都未加工侧槽，取B＝BN；η、γ分别按式(17)、式(18)的计算值取为1.9和0.9。

[0136] (5)差值函数的建立

[0137] 考虑钝化效应，按式(8)得到修正后的裂纹长度ab(i)，其中的η按式(17)计算，其值为1.9。将修正的裂纹长度ab(i)替代式(10)中的ai，得到对应载荷pi时的卸载柔度Ci。对图3中的4个试样的p-V曲线，对最大载荷pmax之前(不包括最大载荷pmax)的所有p-V数据对(pi，Vi)中的载荷pi和位移Vi分别按式(9)和式(12)进行标称化，计算得到标称载荷PN(i)和标称位移Vpl(i)。根据最终裂纹扩展长度af，按照公式(10)得到相应的卸载柔度Cf，标称将p-V曲线上的最终断裂或停机载荷pf的p-V数据对(pf，Vf)，按照式(9)、式(12)得到锚点数组(PN(f)，Vpl(f)。对从标称位移Vpl(i)＞0.001到最大载荷点pmax(不包括最大载荷点pmax)范围内的标称p-V数组(PN(i)，Vpl(i))及锚点数组(PN(i)，Vpl(i))按照公式(13)进行曲线拟合。常温及高温340℃下三点弯曲试样所求得的拟合函数曲线见图4(a)、(b)及图5(a)、(b)，拟合的函数公式见表3；之后按照式(14)建立差值函数 见表4；

[0138] 表3由式(13)拟合的函数公式

[0139]

[0140] 表4测试试样尺寸对应的差值函数

[0141]

[0142]

[0143] (6)裂纹长度与p-V数据对的递推计算

[0144] 根据初始裂纹长度a0和最终裂纹扩展长度af，将裂纹长度划分为若干个不同的长度段a0、a0+Δ、a0+2Δ、a0+3Δ......af，并对任一长度段标记为 其中，取将裂纹长度 代入式(10)计算得到卸载柔度 将初始裂纹长度 (即a0)代入表4
中对应的差值函数 在起点标称位移Vpl(i)＞0.002处将p-V数据对(pi，Vi)依次代入表
4中对应的差值函数 使得满足含有裂纹长度 及p-V数据对(pi，Vi)的差值函数变化率小于等于1％时(即 )，认定此时p-V数据对(pi，Vi)即是对应初始裂纹长度
的p-V数据对 之后将 (即a0+Δ)及 之后的p-V数据对代入表4
中对应的差值函数 并判断差值函数变化率小于等于1％，得到对应裂纹长度 的p-V数据对 同理递推计算可以得到与裂纹长度 相对应的p-V数据对
其中要求 为最终p-V数据对(pf，
Vf)。当 与最终p-V数据对(pf，Vf)不满足差值函数变化率小于等于1％时，调整 之前对应p-V数据对 要求 之前裂纹长度与相应p-V数据对在满足差值函数变化
率小于等于1％的前提下，使得最终p-V数据对(pf，Vf)与 差值函数变化率小于等于1％。
以常温A508材料为例，计算得到的裂纹长度 及相对应的p-V数据对 列于下面
的表5中。

[0145] (7)J-R阻力曲线的构建

[0146] 由递推计算得到裂纹长度 与p-V数据对 后，对应裂纹扩展长度 的J积分 通过式(3)至式(5)、式(10)至式(12)以及式(15)至式(18)计算得到，且这些公式中的ai均用 代替，对应 的Jpl(0)取为零(Jpl(i)由式(15)和式(16)依次迭代计算得到)；用式(19)计算裂纹扩展长度 以常温A508材料为例，计算得到的J积分 和裂纹扩展长度 列于表5：

[0147] 表5常温A508材料裂纹长度与p-V数据对及J积分的计算结果

[0148]

[0149] 以J积分 为纵坐标，以对应的裂纹扩展量 为横坐标，即将所有计算的数据对描在坐标图上，并光滑连接，即得到J-R阻力曲线。

[0150] 图6(a)、(b)分别为用本发明测得的常温下A508和316L的J-R阻力曲线，为了比较，用卸载柔度法测得的J-R阻力曲线也描在图中。由图可见，本发明测量的J-R阻力曲线与卸载柔度法一致，说明本方法是可靠的；然而，卸载柔度法仅得到了1mm以下的J-R阻力曲线，而本发明测得了长尺寸裂纹扩展的J-R阻力曲线。图7(a)、(b)分别为用本发明测得的340℃高温下的A508和316L长尺寸裂纹扩展时的J-R阻力曲线。

[0151] 以上所述的，仅为本发明的较佳实施例，并非用以限定本发明的范围，本发明的上述实施例还可以做出各种变化。即凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰，皆落入本发明专利的权利要求保护范围。本发明未详尽描述的均为常规技术内容。