本文采用 Abaqus 有限元软件计算分析了AP1000 核电安全端焊缝区不同内表面三维周向裂纹(a/t=0.1~0.9和a/c=0.2~1.0)前沿应力强度因子K的分布，为进一步研究分析焊缝区应力腐蚀(SCC)裂纹的扩展形态和寿命预测提供了准确的K分布，为AP1000 核电安全端焊接接头在未来服役中的安全评价和寿命预测奠定基础。
核电的快速发展对解决全球能源危机、减少有害气体排放和应对气候变化的挑战做出积极的贡献。然而核电设备的安全是核电发展的生命线，大规模发展核电，需要全面加强核电设备的设计、制造、建设、运行等环节的安全管理和质量管理。
核电压力容器接管安全端异种金属焊接接头是一种特殊的焊接结构，属于安全重点关注部位。该接头制造难度大，涉及异种金属的焊接极易导致焊接缺陷，且复杂的焊接热循环可引起低韧性材料组织和较高的焊接残余应力会促使应力腐蚀裂纹(SCC)的扩展，最终导致一回路安全端的失效。因此在安全端结构的完整性评价和寿命预测中，需要准确计算SCC裂纹的扩展速率、扩展形态和寿命等。SCC裂纹的扩展速率和形态主要受裂纹前沿应力强度因子K的控制。影响应力强度因子的因素主要包括：安全端所受的载荷、残余应力、裂纹形态、安全端焊接结构的几何和材料等”。因此，本文采用有限元准确计算了安全端焊接接头中裂纹前沿的应力强度因子K分布以有效纳入结构几何、焊接残余应力和接头区复杂材料的影响，为AP1000 核电设备未来服役运行中安全端焊接接头中裂纹扩展的准确分析、安全评价和寿命预测定基础。
有限元模型
安全端几何结构和材料
接管安全端异种钢焊接接头是一种特殊而且复杂的焊接结构，是核反应堆一回路压力边界是安全重点关注部位。图1(a)是反应堆压力容器及安全端模型，由接管嘴与安全端通过异种金属的焊接而构成，其左端与承压设备相连，右端与一回路主管道相连。因此，承压容器、接管安全端与一回路主管道共同构成一回路压力边界，其设计工作压力为15.5MPa左右，温度为288-345℃。各部分材料如图1(b)所示：包括管嘴(A-A508)；焊层(B-A1loy82)；焊缝(C-A1loy182)；管子(D-316L)，括号中为各部分典型材料牌号。
图(a)
图(b)
图1核反应压力容器及管道整体结构图(a)和各部分材料图(b)
对于半椭圆形裂纹可以用两个参数表征周向裂纹尺寸：裂纹的相对深度a/t和缺陷的相对长度a/c，如图2所示。a为裂纹深度，c为裂纹半长，t为接管壁厚。裂纹穿透时a/t=1。
图3是安全端内表面周向未穿透半椭圆裂纹计算不同点K时的表征方法。其中a是裂纹深度，c是裂纹半长，d是半椭圆形内表面裂纹前沿夹角(从内表面点0度开始，到最深点90度方向结束)。夹角用弧度来进行表征，0”对应的弧度为0rad，90°对应的弧度为1.57rad，中间的值按比例变化。
图2 未穿透周向缺陷表征示意图
图3 安全端内表面周向未穿透裂纹的表征
用Abaqus有限元软件进行分析计算。由于结构对称性，选取安全端结构的一半建立三维有限元模型。设置核电一回路管道工作内压 15.5MPa；在安全端右侧与一回路管道连接处利用参考点耦合并施加弯矩(本文考虑弯矩为 2000kN.m)。考虑到安全端及管子尺寸较大，故加入了重力，方法是输入材料的密度属性与重力加速度后即获得相应的重力值。经过单位换算之后输入材料密度为：
图4 焊接残余应力分布
图5 整体网格划分
图6 裂纹前沿网格
计算结果及分析
根据断裂力学知识：应力强度因子K=o(πa)"，在Abaqus软件中可以定义使其直接输出应力强度因子K值。前面已经提到，半椭圆裂纹可以用相对裂纹深度a/t和相对裂纹长度a/c来进行表征。本文主要建立了弯矩为 2000 kNXm作用下 a/t=0.1，0.2，0.3，0.4，0.5，0.6，0.7，0.8，0.9这九组模型，对应每一个裂纹相对深度又分别取a/c=0.2，0.4，0.6，0.8，1 这五个值所以一共为 45 个有限元模型。对于同一载荷下任意深度和长度的裂纹前沿K分布可以通过线性插值的方法得到。
图7不同尺寸(a/t，a/c)裂纹在弯矩 M=2000kN.m的K分布
从上面一系列裂纹前沿的应力强度因子K分布可以看到，并非是裂纹最深点的K值最大。应力强度因子K的分布受到了焊接残余应力的影响，当裂纹比较深的时候K分布和焊接残余应力的分布非常接近。
结论
裂纹前沿应力强度因子K值会直接影响 SCC裂纹扩展的速率，两者存在一个函数的关系。本文用 Abaqus 软件计算出核电安全端三维周向裂纹前沿的K分布对研究核电设备未来服役运行中安全端焊接接头中裂纹扩展的准确分析、安全评价和寿命预测定了基础。返回搜狐，查看更多