弧形双悬臂梁试样裂纹扩展电测接线点布局研究

倪陈强 薛河 赵凌燕 崔英浩 王帅

摘要：直流电位降法是应用于高温高压水环境下裂纹扩展监测试验的主要方法，弧形双悬臂梁试样是一种应力腐蚀裂纹扩展试验的常用试样。为了寻求被测试样的最优接线点布局方案以提高测试精度和稳定性，论文利用有限元软件建立了弧形双悬臂梁试样三维有限元模型，模拟了多组可选的接线点布局方案，求解得到了试样三维电位场分布及其随裂纹变化的规律，并根据计算结果从信号灵敏度、信号幅度和测点位置误差敏感度3方面对试样电流加载点位置和电位差测量点位置的优化选取进行了分析。研究结果表明，不同位置电位差测点下的灵敏度相对大小关系与电流输入位置无关，电位差测点选在试样缺口侧时可获得较高的信号灵敏度;电位差测点与电流加载点的距离是信号幅度大小和测点位置敏感度的主要影响因素，较大的距离会导致信号幅度降低，而较小的距离会导致很高的测点位置误差敏感度。提出的寻优依据和分析方法也可应用于其他类型试样，为进一步研究其他测试精度的影响因素奠定了基础。

关键词：机械工程;裂纹扩展监测;有限元分析;直流电位降法;弧形双悬臂梁试样

中图分类号：TP 206+.1文献标志码：A

DOI：10.13800/j.cnki.xakjdxxb.2019.0218文章编号：1672-9315（2019）02-0309-08

0引言

核电压力容器及管道的应力腐蚀开裂显著降低了核电设备的设计寿命和安全性[1-5]，利用含裂纹试样进行应力腐蚀裂纹扩展试验是核电站安全评价研究工作的重要组成部分之一[6-9]。应力腐蚀裂纹扩展试验中常用的试样有3点弯曲试样[10]、紧凑拉伸试样[11-12]和弧形双悬臂梁试样等，其中，弧形双悬臂梁（Contoured Double Cantilever Beam，CDCB）试样是一种特殊设计的试样，在裂纹扩展过程中可保持应力强度因子恒定[13-15]。在应力腐蚀裂纹扩展试验中，试样被放置于模拟实际工作条件的高压釜密闭环境中，对裂纹扩展很难直接观测，只能用间接方法对裂纹长度进行测量。电位降法是可用于高压釜密闭环境中的一种常用裂纹扩展监测方法[16-17]，其分为直流电位降（Direct Current Potential Drop，DCPD）法和交流电位降（Alternating Current Potential Drop，ACPD）法，原理都是对被测金属试样施加激励电流（直流或交流），通过测量试样特定位置的电位差变化来确定裂纹长度变化量[18]。DCPD法不存在ACPD法的集肤效应，是裂纹扩展试验中实时测量试样裂纹动态扩展长度的主要手段[19]。裂纹长度与电位差的关系与试样形状有关，Johnson研究了无限长有限宽板中心裂纹试样的电位场，并给出了均匀电流密度条件下电位差与裂纹长度关系的解析解[20]，陈篪用解析法研究了3点弯曲试样的电位场和无销钉孔紧凑拉伸试样的电位场[21]。对于各种实际试样的理论标定，一般使用有限元方法给出电位场数值解，并回归得到近似理论标定公式。Ritchie利用有限元法求得了紧凑拉伸试样的电位场，并回归得到了裂纹长度-电位差标定公式。李智军利用有限元法对含半圆表面裂纹试样的电位场进行了有限元分析，为含裂纹试样三维电位场有限元求解奠定了基础[22]。相关研究表明，表征裂纹长度的电位差信号的灵敏度和幅度不但与材料的电导率有关，还与被测试样的形状、电流加载点和电位差测量点的布局密切相关[23-25]。另外，DCPD法裂纹测量中电位差信号非常微小（微伏级）[26]，再加之热电势、仪器温漂、环境电磁干扰等因素的影响，信号的准确测量面临很大困难[27-29]。因此，寻求被测试样最优接线点布局以获得最优的信号灵敏度和信号幅度，对于提高测试精度具有重要意义。文中利用有限元软件建立CDCB试样三维模型，求取不同裂纹长度下的试样三维电位场分布，从信号灵敏度，信号幅度和测点位置误差敏感度3方面对电流加载点和电位差测量点的合理选择进行研究。

1理论基础

电位降法是基于金属材料本身导电特性测量裂纹长度的方法。如图1所示，在含裂纹试样上电位场是裂纹尺寸的函数，在恒定激励电流下，当裂纹长度发生变化时试样的电位场发生变化，通过测定裂纹两侧电位差的变化就可以求出裂纹长度的变化。电位降法一般需要4根引线，与试样形成4个接线点。在图1中，激励电流从外侧两接线点（C1，C2）通过被测试样，电位差信号从内侧两接线点（P1，P2）引出测量。

2有限元模型

如图2所示，取CDCB试样公称宽度W=60 mm，取试样厚度B=16 mm，试样有限元模型使用商用有限元软件ABAQUS建立。需要注意的是裂纹长度a是从销钉孔算起的，这符合常用断裂力学试样标准中裂纹长度的规定。由于试样几何形状相对于裂纹面对称，且实际应用中电流加载点和电位差测量点一般也对称于裂纹面，所以只需建立试样1/2有限元模型，以节约计算成本和提高模型的网格质量。试样有限元网格如图3所示，单元类型为8节点线性热电耦合实体单元（DC3D8E）。为了保证计算结果准确，裂纹附近的单元尺寸被充分细化。材料属性方面，设置材料电导率为核电站压力容器结构材料中常用的304不锈钢的电导率1 388.89 S /mm.有关热传导的参数不作设置，这样ABAQUS将只进行电传导分析。

试样接线点布局包括电流输入点的选择和电位差测量点的选择2方面。考虑5种电流输入位置，如图3所示，均位于垂直于試样厚度方向（Z方向）的中间平面上，电流输入点为I1，I2，I3，I4和I5，电流输出点为这些点相对于裂纹面的对称点。在载荷模块中，向电流输入点加载集中电流，在边界条件中，设置试样韧带为零电势，模拟恒定电流通过试样。在分析步模块中，共建立13个分析步，通过在每个分析步中改变零电势面沿X方向的长度来模拟裂纹长度的增长。取初始裂纹长度a0=0.4 W，从初始裂纹开始按裂纹长度a每增加0.5 mm计算1次，通过13步计算获得所需要的全部计算结果。为了验证仿真计算的正确性，对公式（2）所适用的情况进行了建模分析，计算结果与解析解误差在0.1%以内，这说明所用的仿真计算方法具有足够精度。

电位差测量点的选择有多种可能，如图3所示，考查的电位差测点与电流输入点位于同一平面上，在此平面上取试样边界上所有单元结点以及其相对于裂纹面的对称点作为电位差测量点，共86种电位差测量点布局方案。在图3的半模型中，从试样缺口嘴部开始沿逆时针方向将所有测点依次标记为1，2，…，86，以此来表示选定的电位差测量位置。

3分析和讨论

图4给出了试样电位场的一个示例（电流从I2点输入，a/W=0.45），比较图4（a）和图4（b）可以看出，试样前表面的电位分布与垂直于Z轴的中间纵切面的电位分布差别不大。根据计算结果数据，试样在各纵切面内的电势分布在电流加载点附近差别较大，但在离电流加载点稍远处差别非常小，这说明在中间纵切面讨论电流加载点和电位差测量点是合宜的。

3.1信号灵敏度

电流从I1点和I2点输入时6个测点的电位差U与无量纲化裂纹长度a/W关系曲线如图5所示。从图中可以看出，电位差信号的量级为微伏级，各测点差别较大，电位差随裂纹长度的变化曲线非常平坦，信号灵敏度低。

图6给出了电流从I1，I4和I5输入时10个测点的电位差增量ΔU与无量纲化裂纹长度a/W关系曲线，从图中可见，ΔU-a/W曲线近似呈线性关系，测点1，9，17，23，29，33和39的曲线几乎重合。根据图3，这些测点均位于试样缺口一侧。这说明当测点位于试样缺口侧时，电位差测点位置对信号灵敏度影响不大。比较图6（a）、（b）和（c）可见，电流输入位置的不同对各测点灵敏度的相对关系几乎没有影响。

为了进一步比较所有测点的灵敏度相对关系，图7给出了裂纹长度增量Δa=0.05 W和Δa=0.1 W时5种电流输入位置下所有测点的电位差增量ΔU。从图中可以看出，在测点1（试样缺口嘴处）ΔU最大，随后单调降低。因此，当电位差测点位于试样缺口嘴处时，信号灵敏度最大。从图7还可以看出，电流从I1，I2，I3，I4输入时的曲线几乎重合，而从I5输入时ΔU显著降低。根据图3，电流输入位置I1，I2，I3，I4均位于试样缺口一侧。这说明电流输入位置位于试样缺口侧时，其对信号灵敏度的绝对大小几乎没有影响，而当电流输入点远离试样缺口侧时，各测点灵敏度显著降低。

3.2信号幅度

由于表征裂纹长度的电位差信号非常微小，信号幅度的大小也是选取接线点布局的重要依据。图8给出了电位差测点位于缺口嘴处时5种不同电流输入位置下的U-a/W曲线。从图中可以看出，从I1点输入电流可以获得最大的信号幅度，这时电位差测点和电流加载点重合。为了比较试样边界上所有测点的信号幅度，图9给出了a/W=0.45时5种不同电流输入位置下的电位差与测量位置关系曲线。从图中可见，电流输入点I1，I2和I3对应的信号幅度较大，而I4和I5对应的信号幅度较小。可见为了获得足够的信号幅度，电流输入点应选在试样缺口左侧边界或上侧直线边界部分。

当电流输入位置一定时，信号幅度最大值均出现在电位差测点与电流加载点重合时，这时可以只使用一对导线同时完成电流加载和电位差测量，但实际中通常不采用这种方案，这是因为在这种情况下，电位差信号对导线与试样之间的接触电阻非常敏感，导致测试误差大大增加[30]。从图9还可以看出，信号幅度的最小值与电流输入位置关系不大，均出现在试样右侧。

3.3电位差测量点位置误差影响

在DCPD技术中，导线通常由人工操作点焊机焊接到被测试样上，焊点位置难免会存在误差。电位差测点的位置偏差可分为沿试样边界（X-Y平面内）方向和沿厚度方向（Z方向）2个方面。

在沿试样边界方向上，从图9可以看出，随着电位差测点远离电流输入点，信号幅度快速衰减，然后趋于平缓。因此，当电位差测量点位于电流输入点附近时，电位差随测点位置变化较大，而当电位差测量点远离电流输入点时，测点电位差受测点位置变化影响不大。图10给出了a/W=0.45时所有测点电位差对测点位置的敏感度S（单位毫米的电位差变化量），从图中可见，只要电位差测点距离电流输入点较远（约0.12 W以上），位置误差的影响可以忽略。根据图9和图10也可以推断，焊点大小对试样电位场的影响是存在的，但在距离焊点较远的地方影响可以忽略。

在沿厚度方向方面，图11给出了电流从I1和I2输入时试样中间平面和前表面上沿试样边界各测点的电位差。可见，在靠近电流输入点处，中间平面上和前表面上各点的电位差相差较大，但在稍微远离电流输入点处（距离约0.12 W以上），二者相差非常微小。因此，只要电位差测量点与电流输入点距离较远，在厚度方向上电位差测点位置误差的影响可以忽略。

4结论

1）DCPD法裂纹扩展监测中电位差测点位置和电流输入位置均对信号灵敏度有影响，当电流输入位置位于试样缺口侧时，电位差信号的灵敏度大小主要取决于电位差测点位置。试样边界上各电位差测点的灵敏度相对大小关系与电流输入位置无关，电位差测点选在试样缺口侧时可获得较高的信号灵敏度，选在缺口嘴处时信号灵敏度最大;

2）电流输入点位置的变化对试样电位场分布有相当的影响。当电流输入点与电位差测点重合时，可以获得最大的信号幅度，但是，这种情况下电位差信号对测点位置误差非常敏感。为了避免测点位置误差的影响，电位差测点应选在距离电流输入点较远处，具体距离可用有限元法确定;

3）利用ABAQUS软件建立了CDCB试样的三维有限元模型，对试样的接线点布局方案进行了寻优分析，不但给出了合理的试样接线点布局，提出的尋优依据和分析方法也可应用于其他类型试样，为进一步研究环境温度波动、试样局部应变、裂纹闭合效应等因素对测试精度的影响奠定了基础。

参考文献（References）：

[1]Xue H，Sato Y，Shoji T.Quantitative estimation of the growth of environmentally assisted cracks at flaws in light water reactor components[J].Journal of Pressure Vessel Technology，2009，131（1）： 41-49.

[2]Andresen P L，Morra M M.Stress corrosion cracking of stainless steels and nickel alloys in hightemperature water[J].Corrosion，2008，64（1）：15-29.

[3]郭瑞，薛河，崔英浩.安全端焊接残余应力对裂纹尖端力学参量的影响[J].西安科技大学学报，2018，38（3）：479-483.

GUO Rui，XUE He，CUI Yinghao.Influence of residual stress in safety end on mechanical parameters at crack tip[J].Journal of Xi’an University of Science and Technology，2018，38（3）：479-483.

[4]薛河，崔英浩，赵凌燕，等.压水堆一回路环境中304不锈钢的蠕变特性分析[J].西安科技大学学报，2018，38（1）：156-161.

XUE He，CUI Yinghao，ZHAO Lingyan，et al.Creep characteristics analysis of 304 stainless steel in pressured water reactor primary circuit[J].Journal of Xi’an University of Science and Technology，2018，38（1）：156-161.

[5]薛河，庄泽城，曹婷，等.结构材料维氏硬度与屈服应力的关系分析[J].西安科技大学学报，2017，37（2）：274-279.

XUE He，ZHUANG Zecheng，CAO Ting，et al.Relationship analysis between Vickers hardness and yield stress of structure material[J].Journal of Xi’an University of Science and Technology，2017，37（2）：274-279.

[6]李永强，薛河.核电关键结构材料应力腐蚀裂纹裂尖微观力学特性分析[J].西安科技大学学报，2016，36（3）：380-384.

LI Yongqiang，XUE He.Micromechanical state at SCC tip in nuclear key structure materials[J].Journal of Xi’an University of Science and Technology，2016，36（3）：380-384.

[7]Sato Y，Xue H，Takeda Y，et al.Development of a stress corrosion cracking test methodology using tubeshaped specimens[J].Journal of Testing and Evaluation，2007，35（3）：254-258.

[8]但体纯，吕战鹏，王俭秋，等.690合金在高温水中的应力腐蚀裂纹扩展行为[J].金属学报，2010，46（10）：1267-1274.

DAN Tichun，LV Zhanpeng，WANG Jianqiu，et al.Crack growth behavior for stress corrosion cracking of 690 alloy in high temperature water[J].Acta Metallurgica Sinica，2010，46（10）： 1267-1274.

[9]Jin W K，Choi M R，Yun J K.Effect of loading rate on the fracture behavior of nuclear Piping materials under cyclic loading conditions[J].Nuclear Engineering & Technology，2016，48（6）：1376-1386.

[10]胡梦，陈凯，张乐福.基于直流电压降法的三点弯曲试样疲劳裂纹扩展速率测量方法[J].上海交通大学学报，2015，49（12）：1778-1784.

HU Meng，CHEN Kai，ZHANG Lefu.Measurement on fatigue crack growth rate of threepoint bending specimens based on direct current potential drop method[J].Journal of Shanghai Jiaotong University，2015，49（12）：1778-1784.

[11]杜東海，陈凯，余论，等.直流电压降法测量核电结构材料在空气中的疲劳裂纹扩展速率[J].原子能科学技术，2014，48（8）：1386-1391.

DU Donghai.，CHEN Kai，YU Lun，et al.Measurement of fatigue crack growth rate of reactor structural material in air based on DCPD method[J].Atomic Energy Science & Technology，2014，48（8）：1386-1391.

[12]ASTM E399-09，Standard test method for linearelastic planestrain fracture toughness KIc of metallic materials[S].West Conshohocken： ASTM International，2009.

[13]Lu Z，Shoji T，Meng F，et al.Effects of water chemistry and loading conditions on stress corrosion cracking of coldrolled 316NG stainless steel in high temperature water[J].Corrosion Science，2011，53（1）：247-262.

[14]Lu Z，Shoji T，Takeda Y，et al.Transient and steady state crack growth kinetics for stress corrosion cracking of a cold worked 316L stainless steel in oxygenated pure water at different temperatures[J].Corrosion Science，2008，50（2）：561-575.

[15]Xue H，Li Z，Lu Z，et al.The effect of a single tensile overload on stress corrosion cracking growth of stainless steel in a light water reactor environment[J].Nuclear Engineering & Design，2011，241（3）：731-738.

[16]Kim D S，Ahn S B，Lee K S，et al.The development of crack measurement system using the direct current potential drop method for use in the hot cell[R].Ibaraki： Japan atomic energy research institute，1999.

[17]ern I.The use of DCPD method for measurement of growth of cracks in large components at normal and elevated temperatures[J].Engineering Fracture Mechanics，2004，71（6）：837-848.

[18]Ritchie R O，Bathe K J.On the calibration of the electrical potential technique for monitoring crack growth using finite element methods[J].International Journal of Fracture，1979，15（1）：47-55.

[19]許君.两种结构钢的环境促进裂纹扩展行为与裂纹监测技术[D].北京：机械科学研究总院，2010.

XU Jun.Environmentallyassisted cracking behavior of two structural steels and the monitoring technology[D].Beijing：China Academy of Machinery Science and Technology，2010.

[20]Johnson H H.Calibrating the electric potential method for studying slow crack growth[J].Materials Research And Standards，1965，5（9）： 442-445.

[21]陈凯，杜东海，陆辉，等.用直流电压降法研究316LN不锈钢的疲劳裂纹扩展行为[J].机械工程材料，2016，40（2）：7-10，36.

CHEN Kai，DU Donghai，LU Hui，et al.Research on fatigue crack growth behavior of 316LN stainless steel by direct current potential drop method[J].Materials for Mechanical Engineering，2016，40（2）：7-10，36.

[22]李智军，薛河.直流电位降裂纹测深仪的数值标定[J].西安科技大学学报，2012，32（1）：116-120.

LI Zhijun，XUE He.Numerical calibration of crack monitor based on direct current potential drop[J].Journal of Xi’an University of Science and Technology，2012，32（1）：116-120.

[23]Hicks M A，Pickard A C.A comparison of theoretical and experimental methods of calibrating the electrical potential drop technique for crack length determination[J].International Journal of Fracture，1982，20（2）：91-101.

[24]王威强，李培宁，琚定一.缺口裂纹起裂和扩展检测的直流电位法[J].实验力学，1992，7（1）：51-59.

WANG Weiqiang，LI Peining，JU Dingyi.Detection of notch crack initiation and propagation by DC potential drop method[J].Journal of Experimental Mechanics，1992，7（1）：51-59.

[25]路明旭，肖欣然，姜伟之.一种开裂紧凑试样的电位分析[J].北京航空航天大学学报，1985 （1）：110-117.

LU Mingxu XIAO Xinran JIANG Weizhi.The electrical potential analysis of cracked compact specimens[J].Journal of Beijing Insititute of Aeronautics and Astronautics，1985 （1）：110-117.

[26]Merah N.Detecting and measuring flaws using electric potential techniques[J].Journal of Quality in Maintenance Engineering，2003，9（2）：160-175.

[27]吕宝桐，路民旭.金属材料低温疲劳裂纹起始寿命的微机辅助电位法测试[J].实验力学，1991，6（3）：272-278.

LU Baotong，LU Minxu.An experimental study on low temperature fatigue crack initiation of metals from notches by a microcomputer aided potential drop measurement system[J].Journal of Experimental Mechanics，1991，6（3）：272-278.

[28]劉延利，路民旭.一种新的微机辅助疲劳裂纹扩展综合测试系统[J].北京航空航天大学学报，2000，26（2）：125-128.

LIU Yanli，LU Minxu.Set of new microcomputer aided system for comprehensively measuring FCG[J].Journal of Beijing University of Aeronautics and Astronautics，2000，26（2）：125-128.

[29]张文忠，张伟国.核电站结构部件环境敏感断裂的系统监测和寿命预测[J].腐蚀与防护，1998，19（5）：199-201.

ZHANG Wenzhong，ZHANG Weiguo.System monitoring and life prediction of environmentally assisted cracking of nuclear reactor components[J].Corrosion & Protection，1998，19（5）：199-201.

[30]Sposito G.Advances in potential drop techniques for nondestructive testing[D].London： Imperial College London，2009.