贺屹, 蔡力勋^*,, 陈辉^*,, 彭云强
西南交通大学力学与工程学院应用力学与结构安全四川省重点实验室,成都 610031

A SEMI ANALYTICAL METHOD TO SOLVE J-INTEGRAL FOR MODE-I CRACK COMPONENTS

HeYi, CaiLixun^*,, ChenHui^*,, PengYunqiang
Applied Mechanics and Structure Safety Key Laboratory of Sichuan Province,School of Mechanics and Engineering,Southwest Jiaotong University,Chengdu 610031,China
中图分类号:O346.1
文献标识码:A

通讯作者:通讯作者:蔡力勋,教授,博士生导师,主要研究方向：断裂与疲劳、强度理论、材料测试理论与技术. E-mail:lix_cai@263.net;陈辉,博士研究生,主要研究方向：材料测试理论与技术. E-mail：chen_hui5352@163.com通讯作者:蔡力勋,教授,博士生导师,主要研究方向：断裂与疲劳、强度理论、材料测试理论与技术. E-mail:lix_cai@263.net;陈辉,博士研究生,主要研究方向：材料测试理论与技术. E-mail：chen_hui5352@163.com
收稿日期:2018-01-24
接受日期:2018-04-4
网络出版日期:2018-06-10
版权声明:2018《力学学报》编辑部《力学学报》编辑部 所有
基金资助:国家自然科学基金资助项目(11472228).

展开

摘要
表征裂纹尖端应力应变场程度的J积分是一个定义明确、理论严密的弹塑性断裂力学基础参量. 目前J积分的计算主要是依靠塑性因子法和有限元法,但对各类裂纹构元获得J积分以及载荷-位移关系的解析公式以实现材料断裂韧性理论预测和材料测试是断裂力学的重要和困难的任务. 以J积分为参量的材料断裂测试中应用最广的是I型裂纹试样的断裂韧性测试. 本文在平面应变条件下,针对断裂韧性测试中使用的6种I型裂纹构元,基于能量等效假设,提出了J积分-载荷和载荷-位移的工程半解析统一表征方法,进而结合有限元分析的少量计算获得J积分-载荷和载荷-位移关系的半解析公式待定参数. 分析表明,6种I型裂纹构元的J积分-载荷和载荷-位移统一公式的预测结果与有限元结果吻合良好. 新提出的J积分-载荷工程半解析公式包含了材料的弹性模量、应力强度系数和应变硬化指数,能够广泛适应不同的材料,且运用该公式能够方便获取任意载荷点对应的J积分值. 应用新方法可便于获得各类I型裂纹构元的J积分-载荷和载荷-位移工程半解析公式.

关键词：J积分;塑性断裂;半解析公式;能量等效假设;I型裂纹
Abstract
The J-integral to characterize the singular level of the stress and strain field at the crack tip is definite and rigorous and is a basic parameter of elastoplastic fracture mechanics. The calculation of J-integral mainly depends on the plastic factor method and the finite element method at present. For theoretical predicting and testing of material fracture toughness, it is important and difficult to obtain analytical expressions about J-integral-load and load-displacement relations of cracked components. The most widely used test for structure integrity evaluation with J-integral is the ductile fracture toughness of type-I cracked specimens. Here, based on the Chen-Cai energy equivalence hypothesis, a unified characterization method of J-integral-load and load-displacement relation is proposed for six Mode-I cracked components which are commonly used in fracture toughness test under the plane strain condition. Then, the undetermined parameters of the engineering semi-analytical formulas of the J-integral-load and the load-displacement relations are obtained by a small amount of finite element analysis. The results show that the J-integral-load and load-displacement relation predicted by the unified semi-analytical formulas are in good agreement with those from finite element method. The engineering semi-analytical J-integral-load formula, which contains the elastic modulus, stress strength coefficient and strain hardening exponent of materials, can be widely adapted for different materials. And the J-integral value corresponding to arbitrary load points can be easily obtained by the formula. The presented novel method is convenient to establish the engineering semi-analytical formulas of J-integral-load and load-displacement relations for various type-I cracked components or specimens.

Keywords：J-integral;plastic fracture;semi-analytical formula;energy equivalence hypothesis;mode I crack

-->0
PDF (4576KB)元数据多维度评价相关文章收藏文章
本文引用格式导出EndNoteRisBibtex收藏本文-->
贺屹, 蔡力勋, 陈辉, 彭云强. 求解I型裂纹构元J积分的半解析方法[J]. 力学学报, 2018, 50(3): 579-588 https://doi.org/10.6052/0459-1879-18-026
He Yi, Cai Lixun, Chen Hui, Peng Yunqiang. A SEMI ANALYTICAL METHOD TO SOLVE J-INTEGRAL FOR MODE-I CRACK COMPONENTS[J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(3): 579-588 https://doi.org/10.6052/0459-1879-18-026

引言

随着科技的发展,新材料的广泛使用,各种工程建筑和结构^[1,2]不断涌现,研究其力学行为以保障安全就显得尤为重要. 弹塑性断裂力学已经在工程结构设计中得到切实应用,航空部门提出了损伤容限设计思想,在选材方面有了专门的指导手册,美国电力研究院发展了以J积分为基础的弹塑性缺陷评估方法,称为EPRI方法^[3,4,5].J积分是一个定义明确,理论严密的应力、应变场参量,与线弹性断裂力学中的应力强度因子一样,J积分既能表征裂纹尖端区域应力、应变场的程度,又因与加载过程中的能量相关而容易通过实验来测定,是弹塑性断裂力学的基础参量^[6].
对于平面应变条件下断裂韧性 Jc测试,现已发展的较为成熟,Zhu等^[7]从线弹性断裂力学和弹塑性断裂力学的角度对金属材料的断裂韧性测试、评定和标准化进行了较为详细的技术综述,嵇醒^[8]对断裂力学中的判据作了扼要的综述,归纳出断裂力学判据中目前还没有较好解决的几个问题. 在现有断裂韧性评定方法中,最常见的是采用裂纹尖端具有高约束度的标准紧凑拉伸(compact tension,CT)试样或单边裂纹弯曲(single edged notched bending,SEB)试样进行断裂力学试验以获取J,这些试样已在若干测试标准^[9,10,11]中得到推荐,其中CT试样使用最广,已经涉及到蠕变分析^[12]、残余应力分析^[13]、复合材料^[14]断裂韧性测试等方面. 由Feddern和Macherauch^[15]提出的圆形紧凑拉伸试样(round compact tensile,RCT)是一种对圆棒材料裂韧性测试非常方便的试样,目前已运用到各种合金的疲劳裂纹扩展研究^[16,17]. 单边裂纹拉伸(single edged notched tension,SENT)试样因其约束特性与含裂纹管道材料约束性相近而在管道材料的断裂韧度评定中得到重视^[18,19],该试样对焊接材料的力学性能测试^[20]同样方便. 小尺寸构件不仅取样方便且能有效降低试验成本,因而研究非标准小试样具有重要理论意义和工程应用价值. Bao等^[21]开发了一种实际尺寸仅有一元硬币大小的含内侧边裂纹C形拉伸 (C-shaped inside edge-notched tension,CIET)小试样,实现了金属材料的疲劳裂纹扩展和断裂行为测试. 但晨等^[22]基于量纲一载荷分离理论建立了可用于测 定CIET试样延性断裂行为的规则化法,并完成了J阻力曲线测试. 双边裂纹板拉伸(double edge crack tension,DECT)试样也是一种在研究中使用较多的裂纹构元^[23].
现有的方法中,对于J积分的计算可以分为两大类：一类是按照简化模型或者工程估算方法进行计算,如EPRI法,另一类是通过有限元分析(finite element analysis,FEA)或者实验方法获取J积分值,如柔度法. EPRI手册J积分估算法的关键是将弹塑性J积分分解为弹性 Je和塑性 Jp两部分^[24,25,26],即
$J = J_{\rm e} + J_{\rm p} $(1)
其中弹性部分为
$J_{\rm e} = f_1 (a_{\rm e} )\dfrac{P^2}{{E}'} $(2)
式中, E'=E/(1-ν2)对应平面应变问题, E'=E对应平面应力问题; E为弹性模量; ν为泊松比. f1为与试件裂纹 大小、位置等因素有关的形状函数, ae为修正的有效裂纹长度. 塑性部分^[27]表示为
$J_{\rm p}= \alpha \sigma _0 \varepsilon _0 ch_1 \left( {\dfrac{P}{P_0 }}\right)^{n + 1} $(3)
式中, α为材料硬化系数, σ0和 ε0为屈服应力和屈服应变, n为硬化指数, c为试样剩余韧带宽度, P0为试样的参考屈服载荷; h1为 J积分全塑性解的函数,它与材料、结构、载荷以及裂纹尺寸有关. 通过有限元计算,可以得到 h1的值,EPRI手册针对多种试样提供了 h1随 n、试样几何尺寸变化的系列离散数据表格.
在断裂韧性试验时, 若采用ASTM E1820-15^[28],J积分的弹性分量 Je表示为
$J_{\rm e} = \dfrac{K^2}{{E}'} (4)$
式中, K为应力强度因子,塑性分量为
$J_{\rm p} = \eta \dfrac{U_{\rm p} }{B_{\rm N} b} $(5)
式中, b表示剩余韧带厚度, BN表示试样净厚度, Up表示塑性功如图1所示,可由载荷( P)和加载线位移( h)的试验曲线获得, η为无量纲塑性因子^[29,30]. η因子多由有限元法分析确定^[31],其步骤为建立有限元网格模型,设置材料属性,施加约束及载荷,提交运算,进而提取 J积分值,减去由 K因子得到的弹性部分,而根据 式(5)逆向求解. 在平面应变下, η近似与几何尺寸相关,在一般三维约束条件下, η不仅与几何尺寸相关,也与 n相关,但要获得 η的关于几何尺寸、硬化指数的统一表达式十分困难.
显示原图|下载原图ZIP|生成PPT
图1塑性功示意图
-->Fig. 1The schematic diagram of plastic work
-->

考虑Rice在1968年针对I型裂纹构元提出的定义式^[32]
$J = - \Big(\dfrac{\partial U}{B\partial a} \Big) $(6)
式中, B表示试样厚度, a为裂纹长度. 本文基于Chen-Cai能量等效假设,采用 J积分能量定义式(6),针对不同几何的I型裂纹,提出了 J积分的工程半解析表征方法,得到 J积分和载荷-位移的半解析统一公式.

1 单向加载下构元弹塑性问题的半解析方法

Chen和Cai^{[33,34,35,36]}通过建立变形域的应变能与积分中值点位置的材料RVE(representative volume element)的应变能密度与变形域有效体积等效,以及考虑Mises等效使得复杂应力状态下RVE能量密度与单轴应力状态下RVE能量密 度等效,进而建立反映连续固体的能量、载荷、位移和材料本构关系参数关系的通用模型.
对于幂硬化材料,考虑Ramberg-Osgood应力应变关系,即
$\varepsilon = \dfrac{\sigma }{E} + \Big (\dfrac{\sigma }{C_K } \Big)^N $(7)
式中, ε表示应变, σ表示应力, CK表示应力强度系数, N表示应力硬化指数,对于纯塑性材料,有
$\varepsilon _{\rm p} =\Big (\dfrac{\sigma }{C_K } \Big)^N $(8)
由Chen-Cai能量等效假设,可得塑性应变能 Up为
$U_{\rm p } = \dfrac{NC_K }{N + 1}V_{\rm P} \varepsilon _{\rm p}^{1 / N + 1} = \dfrac{NC_K V^\ast }{N + 1}\dfrac{V_{\rm P} }{V^\ast }\varepsilon _{\rm p p}^{1 / N + 1} $(9)
式中, V表示变形域有效体积, V*表示特征体积,且 V*=A*h*, A*为特征面积, h*表示特征位移. 特征体积、特征面积、特征位移旨在用于实现不同 a/W条件下载荷、位移、应变、能量的无量纲化. 对材料的塑性有效体积和塑性等效应变做如下幂律假定
$\left. \begin{array}{l} \dfrac{V_{\rm p} }{V^\ast } = k_1 (\dfrac{h_{\rm p} }{h^\ast })^{k_2 } \\ \varepsilon _{\rm p} = k_3 (\dfrac{h_{\rm p} }{h^\ast })^{k_4 } \end{array}\!\! \right\} $(10)
其中, k1和 k2分别为有效体积系数和指数, k3和 k4分别为有效应变系数和指数, hp表示加载点在加载方向移塑性位移, h*表示特征位移,将式(10)代入式(9)可得
$U_{\rm p} = \dfrac{NC_K V^\ast }{N + 1}k_1 k_3^{1 + 1 / N} \Big(\dfrac{h_{\rm p} }{h^\ast } \Big )^{k_4 / N + k_4 + k_2 } $(11)
求导可得载荷 P的半解析公式
$ P = \dfrac{\partial U_{\rm p} }{\partial h_{\rm p} } = \dfrac{NC_K V^\ast}{(N + 1)h^\ast }k_1 k_3^{1 + 1 / N} \cdot \\ \Big (\dfrac{k_4 }{N} + k_4 + k_2\Big ) \Big (\dfrac{h_{\rm p} }{h^\ast }\Big)^{k_4 / N + k_4 + k_2 - 1} $(12)
整理可得,在塑性状态下材料的应变能与载荷、位移之间的关系为
$\left.\begin{array}{ll} \dfrac{U_{\rm p} }{U_{\rm p}^\ast } =\left\{\begin{array} \Big (\dfrac{h_{\rm p} }{h^\ast } \Big )^{m_{\rm p} + 1} \\ \Big (\dfrac{P}{P_{\rm p}^{\ast } } \Big)^{1 / m_{\rm p} + 1} \end{array}\right. \\ \dfrac{P}{P_{\rm p}^{\ast } } = \Big (\dfrac{h_{\rm p} }{h^\ast } \Big )^{m_{\rm p} } \end{array} \right\} $(13)
其中
$\left.\begin{array}{l} U_{\rm p} ^{\ast } = \dfrac{NC_K V^\ast }{N + 1}k_1 k_3^{1 + 1 / N} \\ m_{\rm p} = \dfrac{k_4 }{N} + k_4 + k_2 - 1 \\ P_{\rm p}^{\ast } = \dfrac{(1 + m_{\rm p} )NC_K V^\ast }{(N + 1) h^\ast }k_1 k_3^{1 + 1 / N} \end{array}\!\! \right\} $(14)
类似地,在线弹性范围内可以得到材料的弹性能与载荷、位移之间的关系为
$\left.\begin{array}{ll} \dfrac{U_{\rm e} }{U_{\rm e}^\ast }=\left\{\begin{array} \Big (\dfrac{h_{\rm e} }{h^\ast }\Big)^2 \\ \Big (\dfrac{P }{P_{\rm e}^\ast }\Big)^2 \end{array}\right. \\ \dfrac{P }{P_{\rm e}^\ast } = \dfrac{h_{\rm e} }{h^\ast } \end{array} \right\} $(15)
式中
$ \left. \begin{array}{l} U_{\rm e} ^{\ast} = \dfrac{k_0 EV^\ast }{2} \\ P_{\rm e}^{\ast} = k_0 EA^\ast \end{array} \right\} $(16)

2 I型裂纹构元的半解析统一公式

对于弹塑性变形,在单向加载下裂纹构元的总能量可以通过弹性能和塑性能进行工程叠加获得,即
U(P,h)=Ue(P,he)+Up(P,hp)(17)
加载线位移
h(P)=he(P)+hp(P)(18)
将式(13)和式(15)代入式(18)即可得裂纹构元的载荷-位移统一公式
$\Big(\dfrac{PA_1 }{A^\ast }\Big)^{\tfrac{1}{m_{\rm p} }} + \dfrac{PA_2 }{A^\ast }= \dfrac{h}{h^\ast } $(19)
式中
$ \left. \begin{array}{l} A_1 = \dfrac{(1 + N)}{(1 + m_{\rm p} )NC_K k_1 k_3^{1 + 1 / N} } \\ A_2 = \dfrac{1}{k_0 E} \end{array} \right\} $(20)
且
$ \left.\begin{array}{l} A\ast = WB \Big(1 - \dfrac{a}{W} \Big)^m \\ h\ast = W \\ V\ast = W^2B\Big(1 - \dfrac{a}{W} \Big)^m \end{array} \right\} $(21)
式(21)中 W表示试样宽度, m表示与裂纹长度 a相关的特征体积折减系数,对于特定的裂纹构元,其值可固定为一个适当的常数. 联立式(13)、式(15)和式(17)得到弹塑性条件下应变能为
$ U = W^2B \Big(1 - \dfrac{a}{W}\Big)^m\left[ \dfrac{NC_K }{N + 1}k_1 k_3^{1 + 1 / N} \Big(\dfrac{P}{P_{\rm p}^{\ast}} \Big)^{1 / m_{\rm p} + 1} +\right. \\ \left. \dfrac{k_0 E}{2} \Big(\dfrac{P }{P_{\rm e}^\ast }\Big)^2 \right] $(22)
联立式(22)和式(6),可得裂纹构元的 J积分-载荷统一公式为
$ J = \Big (1 - \dfrac{a}{W}\Big)^{m - 1}\Big [a_1 \Big (\dfrac{P}{P_{\rm p}^\ast}\Big)^{a_2 } + a_3 \Big (\dfrac{P}{P_{\rm e}^\ast }\Big)^2\Big] $(23)
式中
$\left.\begin{array}{ll} a_1 = Wm\dfrac{NC_K }{N + 1}k_1 k_3^{1 + 1 / N} \\ a_2 = 1 + \dfrac{1}{m_{\rm p} } \\ a_3 = Wm\dfrac{k_0 E}{2} \end{array}\right\} $(24)
如式(23)所示的 J积分-载荷统一公式, k0~k4皆为式(14)和式(16)中的未知参数,可通过载荷位移曲线获取. 同时, J积分统一公式又可以表示成如下无量纲化形式
$ \dfrac{J}{J^\ast } = a_4 \Big (\dfrac{P}{P^\ast }\Big)^{a_2 } + a_5 \Big (\dfrac{P}{P^\ast }\Big)^2 $(25)
式中
$\left.\begin{array}{ll} P^\ast = (P_{\rm p}^\ast P_{\rm e}^\ast )^{0.5} \\ J^\ast = Wm \Big(1 - \dfrac{a}{W} \Big)^{m - 1}C_K \end{array} \right\} $(26)
$ \left. \begin{array}{l} A = \dfrac{P_{\rm e}^\ast }{P_{\rm p}^\ast } = \dfrac{(N + 1)Ek_0}{(1 + m_{\rm p} )NC_K k_1 k_3^{1 + 1 / N} } \\ a_4 = \dfrac{NA^{\tfrac{1}{2}(1 + \tfrac{1}{m_{\rm p} })}}{N + 1}k_1 k_3^{1 +1 / N} \\ a_5 = \dfrac{(1 + m_{\rm p} )N}{2 (N + 1)} k_1 k_3^{1 + 1 / N} \end{array} \right\} $(27)
对于 h不是加载线位移的I型裂纹构元情形,上述公式的形式依然可假定成立.

3 J积分公式参数的确定方法

采用有限元分析通过少量计算即可获得 J积分统一公式参数,其计算过程如图2所示. 采用有限元软件ANSYS14.5计算6种I型 裂纹构元的 J积分统一公式的参数,全部采用平面应变条件,模型单元采用plane183单元,泊松比 υ=0.3,应力强度系数 CK=1000MPa,线弹性计算时取弹性模量 E=210GPa. 表1给出了各裂纹构元的尺寸,现以CT试样为例说明参数求解过程.
显示原图|下载原图ZIP|生成PPT
图2参数获取流程框图
-->Fig. 2Flow chart of parameters acquisition
-->

Table 1
表 1
表 1各类裂纹构元的尺寸
Table 1The size of various cracked components


新窗口打开
有限元分析中,单元网格的划分对计算结果可能存在影响,本文对不同的裂尖网格密度条件下的载荷位移曲线进行了比较. 定义裂尖2 mm×1 mm范围网格数量为104时为一倍网格密度,在此基础上对裂尖区域网格加密,其余部分网格密度不变. 图3给出了网格密度对载荷位移曲线的影响,结果表明采用一倍网格时已经满足计算要求.
显示原图|下载原图ZIP|生成PPT
图3裂尖网格对载荷位移曲线的影响
-->Fig. 3Effect of crack tip grid on load-displacement curves
-->

如图2的流程图所示,在有限元软件前处理时输入相应参数,计算得到不同 a/W条件下的载荷位移曲线,进行几何无关 处理( P/A*-h/h*)后可得参数 m,同时求得 a/W的适用范围,结果如图4所示. 对模型进行全塑性计算(文中模型均取 a/W=0.5),根据 式(13),可得 Pp*和 mp,并将其代入 式(14)即可得到 k1~k4,对模型进行线弹性计算,根据 式(15)可求得 k0.
计算出CT试样 J积分统一公式中各参数数值分别为 m=2.282, k0=0.1429, k1=2.2456, k2=0.0038, k3=0.3170, k4=1.0033,公式对 a/W的适用范围为0.45 ~0.75. 图5给出了全塑性条件 J积分公式预测载荷位移曲线与有限元结果的对比,弹塑性条件下公式预测的载荷位移曲线与有限元的比 较由图6(a)给出. 图7(a)和图7(b)给出了弹塑性条件下公式预测的 J积分与有限元结果的对比.
显示原图|下载原图ZIP|生成PPT
图4CT试样几何无关载荷位移曲线
-->Fig. 4Geometry-independent load-displacement curves of CT specimen
-->

显示原图|下载原图ZIP|生成PPT
图5预测与计算的CT全塑性载荷-位移曲线比较
-->Fig. 5Comparison of fully plastic load-displacement curves predicted by formula and those from FEA for CT specimen
-->

同理,本文还得到了SEB和SENT等其他5种试样 J积分统一公式的未知参数,表2给出了其参数数值,图6给出了SEB、SENT等试样 载荷位移曲线公式预测与有限元结果对比,图7给出了各试样构元在 CK=100MPa ~2 000 MPa, N=2.5~10时统一公式对J积分的预测结果与有限元对比.
Table 2
表 2
表 2各类裂纹构元的J积分统一公式参数
Table 2Unified formula parameters of J-integral for various cracked components


新窗口打开
显示原图|下载原图ZIP|生成PPT
图6式(19)预测的弹塑性条件下载荷位移曲线与有限元比较
-->Fig.6Comparison of elastoplastic load-displacement curves predicted by Eq.(19) with those from FEA
-->

显示原图|下载原图ZIP|生成PPT
图7式(23)预测的J积分与有限元结果比较
-->Fig. 7Comparison of J-integral predicted by Eq.(23) with those from FEA
-->

3 结 论

(1) 基于Chen-Cai能量等效假设,结合 J积分能量定义式,提出了一种解析求解I型裂纹构元 J积分的方法, 这种新方法 旨在针对一系列I型裂纹标准试样和非标准试样可以较易实现获得 J积分半解析的表达式,以便用于测试和断裂问题的理性分析.
(2)对于6种I型裂纹构元,通过有限元的计算给出了 J积分-载荷和载荷-位移统一公式的参数;依靠 J积分统一公式预测的 J积分结果相较有限元结果都吻合较好.
The authors have declared that no competing interests exist.

---

参考文献文献选项 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1]	胡海岩, 赵永辉, 黄锐. 飞机结构气动弹性分析与控制研究 . 力学学报, 2016, 48(1): 1-27 [本文引用: 1] (Hu Haiyan, Zhao Yonghui, Huang Rui.Studies on aeroelastic analysis and control of aircraft structures .Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(1): 1-27 (in Chinese)) [本文引用: 1]
[2]	柳占立, 庄茁, 孟庆国等. 页岩气高效开采的力学问题与挑战 . 力学学报, 2017, 49(3): 507-516 [本文引用: 1] (Liu Zhanli, Zhuang Zhuo, Meng Qingguo, et al.Problems and challenges of mechanics in shale gas efficient exploitation .Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(3): 507-516 (in Chinese)) [本文引用: 1]
[3]	Kumar V, German MD, Shih CF.Engineering Approach for Elastic-plastic Fracture Analysis. General Electric Co., Schenectady, New York: Corporate Research and Development Dept., 1981 [本文引用: 1]
[4]	Kumar V, German MDWWW.Development on Elastic-plastic Fracture Analysis (EPRI NP-3607). Palo Alto, California: Research & Development Center of EPRI, 1984 [本文引用: 1]
[5]	Zahoor A.Ductile Fracture Handbook. vols 1-3, EPRI NP-6301-D/N14- 1, Novetech Corp . and Electric Power Research Institute (Palo Alto, CA), 1989, 1990, 1991 [本文引用: 1]
[6]	Rice JR.A path independent integral and the approximate analysis of strain concentration by notches and cracks .Journal of Applied Mechanics, 1968, 35(2): 379-386 [本文引用: 1]
[7]	Zhu XK, Joyce JA.Review of fracture toughness (G, K, J, CTOD, CTOA) testing and standardization .Engineering Fracture Mechanics, 2012, 85: 1-46 [本文引用: 1]
[8]	嵇醒. 断裂力学判据的评述 . 力学学报, 2016, 48(4): 741-753 [本文引用: 1] (Ji Xing.A critical review on criteria of fracture mechanics .Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(4): 741-753 (in Chinese)) [本文引用: 1]
[9]	GB/T 21143-2014. 金属材料准静态断裂韧度的统一试验方法. 北京:中国标准出版社,2015 [本文引用: 1] (GB/T 21143-2014. Metallic materials-unified method of test for determination of quasistatic fracture toughness. Beijing:Standards Press of China, 2015 (in Chinese)) [本文引用: 1]
[10]	ASTM E399-09. Standard test methods for linear-elastic plane-strain fracture toughness Klc of metallic materials . Annual Book of ASTM Standards: Vol.3.01. PhiladelPhia, PA: Ameriean Society for Testing and Materials, 2009 [本文引用: 1]
[11]	ASTM E1820-09. Standard test methods for measurement of fracture toughness// Annual Book of ASTM Standards: Vol.3.01. PhiladelPhia, PA: Ameriean Society for Testing and Materials, 2009 [本文引用: 1]
[12]	Sun PJ, Wang GZ, Xuan FZ, et al.Three-dimensional numerical analysis of out-of-plane creep crack-tip constraint in compact tension specimens .International Journal of Pressure Vessels and Piping, 2012, 96: 78-89 [本文引用: 1]
[13]	Traore Y, Paddea S, Bouchard PJ, et al.Measurement of the residual stress tensor in a compact tension weld specimen .Experimental Mechanics, 2013, 53(4): 605-618 [本文引用: 1]
[14]	Quino G, El Yagoubi J, Lubineau G.Characterizing the toughness of an epoxy resin after wet aging using compact tension specimens with non-uniform moisture content .Polymer Degradation and Stability, 2014, 109: 319-326 [本文引用: 1]
[15]	Feddern G, Macherauch E.New specimen geometry for KIC measurements. Z .Metallkunde, 1973, 64(12): 882-884 [本文引用: 1]
[16]	Szkodo M, Bień A. Influence of laser processing of the low alloy medium carbon structural steel on the development of the fatigue crack .Surface and Coatings Technology, 2016, 296: 117-123 [本文引用: 1]
[17]	Klysz S, Lisiecki J, Leski A, et al.Least squares method modification applied to the NASGRO equation .Journal of Theoretical and Applied Mechanics, 2013, 51(1): 3-13 [本文引用: 1]
[18]	Chiesa M, Nyhus B, Skallerud B, et al.Efficient fracture assessment of pipelines. A constraint-corrected SENT specimen approach .Engineering Fracture Mechanics, 2001, 68(5): 527-547 [本文引用: 1]
[19]	Kingklang S, Uthaisangsuk V.Plastic deformation and fracture behavior of X65 pipeline steel: Experiments and modeling .Engineering Fracture Mechanics, 2018, 191: 82-101 [本文引用: 1]
[20]	Hertelé S, De Waele W, Verstraete M, et al.J-integral analysis of heterogeneous mismatched girth welds in clamped single-edge notched tension specimens .International Journal of Pressure Vessels and Piping, 2014, 119: 95-107 [本文引用: 1]
[21]	Bao C, Cai LX, Dan C.Estimation of fatigue crack growth behavior for small-sized C-shaped inside edge-notched tension (CIET) specimen using compliance technique .International Journal of Fatigue, 2015, 81: 202-212 [本文引用: 1]
[22]	但晨, 蔡力勋, 包陈等. 用于断裂韧度测试的C形环小试样的规则化方法与应用 . 机械工程学报, 2015, 51(14): 54-65 [本文引用: 1] (Dan Chen, Cai Lixun, Bao Chen, et al.Normalization method used to determine fracture toughness with C-ring small sized specimen and its application .Journal of Mechanical Engineering, 2015,51(14):54-65 (in Chinese)) [本文引用: 1]
[23]	Treifi M, Oyadiji SO, Tsang DKL.Computations of the stress intensity factors of double-edge and centre V-notched plates under tension and anti-plane shear by the fractal-like finite element method .Engineering Fracture Mechanics, 2009, 76(13): 2091-2108 [本文引用: 1]
[24]	Bucci RJ, Paris PC, Landes JD, et al.J integral estimation procedures //Fracture Toughness: Part II. ASTM International, 1972 [本文引用: 1]
[25]	Ernst HA, Paris PC, Landes JD.Estimations on J-integral and tearing modulus T from a single specimen test record //Fracture Mechanics. ASTM International, 1981 [本文引用: 1]
[26]	Zahoor A.Fracture of circumferentially cracked pipes .Journal of Pressure Vessel Technology, 1986, 108(4): 529-531 [本文引用: 1]
[27]	Zahoor A.Evaluation of J-integral estimation scheme for flawed throughwall pipes .Nuclear Engineering and Design, 1987, 100(1): 1-9 [本文引用: 1]
[28]	Goldman NL, Hutchinson JW.Fully plastic crack problems: the center-cracked strip under plane strain .International Journal of Solids and Structures, 1975, 11(5): 575-591 [本文引用: 1]
[29]	ASTM E1820-15: International. Standard test method for measurement of fracture toughness . ASTM International, 2015 [本文引用: 1]
[30]	Merkle JG, Corten HT.A J-integral analysis for the compact specimen, considering axial force as well as bending effects . Paper no. 74, PVP 33. American Society for Testing and Materials, 1974 [本文引用: 1]
[31]	包陈, 蔡力勋, 但晨等.考虑裂尖约束效应的小尺寸 CIET试样延性断裂行为 . 机械工程学报, 2017, 53(2): 34-44 [本文引用: 1] (Bao Chen, Cai Lixun, Dan Chen, et al.Ductile fracture behavior of small-sized CIET specimen considering crack front constraint effect .Journal of Mechanical Engineering, 2017, 53(2): 34-44 (in Chinese)) [本文引用: 1]
[32]	Rice JR.A path independent integral and the approximate analysis of strain concentration by notches and cracks .Journal of Applied Mechanics, 1968, 35(2): 379-386 [本文引用: 1]
[33]	Chen H, Cai LX.Theoretical model for predicting uniaxial stress-strain relation by dual conical indentation based on equivalent energy principle .Acta Materialia, 2016, 121: 181-189 [本文引用: 1]
[34]	Chen H, Cai LX.Unified elastoplastic model based on strain energy equivalence principle .Applied Mathematical Modelling, 2017, 52: 664-671 [本文引用: 1]
[35]	Chen H, Cai LX.Unified ring-compression model for determining tensile properties of tubular materials .Materials Today Communications, 2017, 13: 210-220 [本文引用: 1]
[36]	Chen H, Cai LX.Theoretical conversions of different hardness and tensile strength for ductile materials based on stress-strain curves .Metallurgical and Materials Transactions A, 2018, 49(4): 1090-1101 [本文引用: 1]