《武汉工程大学学报》 2012年12期
50-53
出版日期：2013-01-11
ISSN:1674-2869
CN:42-1779/TQ
含裂纹法兰接管的应力强度因子数值计算
0引言法兰连接是炼油装置中连接设备与管路的常见元件.由于炼油过程中腐蚀性介质的影响，接管与法兰连接处环焊缝高应力区容易发生开裂，给管路运行造成重大安全隐患.对于这种情况往往只能采取临时补焊或整体打包等方式进行处理，如图1、2所示.为了确保装置的长周期安全运行，有必要对含裂纹的法兰接管情况进行断裂力学参数分析，为装置的安全评定提供参考依据.图1补焊处理的接管法兰
Fig.1The welded nozzle flange图2整体打包处理的接管法兰
Fig.2The overall packaged nozzle flange下面利用有限元方法对法兰接管组合结构环焊缝处的应力强度因子进行分析.1含裂纹的法兰接管有限元分析1.1概述裂纹应力强度因子的计算方法有三大类：解析法、数值法和实验方法.解析法只能计算简单裂纹，对于复杂结构的裂纹难以获得精确解析解\[13\].实验方法是能解决一些较复杂的裂纹问题，但是实验法在裂纹模型预制和裂纹前沿多点数据测量上还存在不足\[4\].数值方法中的三维有限元法以其强大的模拟和数值计算功能已经被广泛用于应力强度因子的计算\[510\].设r为裂纹尖端附近某点到裂纹尖端的距离，因裂纹尖端应力的各个分量都与r成反比，当r趋近于0时，在裂纹尖端，应力的各个分量都会趋于无穷大，如图3所示，这种特性就叫做应力奇异性.图3裂纹尖端应力奇异性
Fig.3Stress singularity in the model of crack tip由于裂纹尖端应力奇异性，为获得较高的计算精度，采用比较成熟的14节点法.该方法采用14奇异单元，该单元是使围绕裂纹尖端的第一圈等参单元各边的中间节点移到离裂纹尖端的14处\[11\].第12期何家胜，等：含裂纹法兰接管的应力强度因子数值计算
武汉工程大学学报第34卷
1.2有限元模型的建立研究对象为脱硫系统换热器401/2壳程上方的法兰及接管，如图4所示.法兰上方焊接弯头接管，在环焊缝处产生了角度为α的穿透性裂纹，如图5所示.该法兰为PN2.5级带颈对焊管法兰，弯头接管内半径r为75 mm，壁厚t为4.5 mm，弯曲半径RM为225 mm ,弯头另一端接长度为1 000 mm的直管，材料都为20号钢.
图4法兰弯头接管结构
Fig.4The Structure of nozzle flange图5法兰弯头接管裂纹横断面
Fig.5The crack cross-section of nozzle flange1.2.1建立裂纹尖端奇异单元对于三维裂纹，其奇异单元不能自动生成，必须采用手动建立.因此建模时采用自下而上的建模方法，先建立节点，再由节点生成单元.在裂纹前沿沿其长度方向平分为9段，即形成10层节点，每一层节点是以裂纹尖端节点为圆心，形成10圈同心圆节点，如图6所示，通过连接各个节点生成裂纹前沿单元，所有的单元都是8节点的Solid45单元.然后将围绕裂纹前沿的第一圈单元转化成奇异单元.
图6裂纹尖端奇异单元
Fig.6The singular element in the model of crack tip1.2.2法兰接管整体模型的建立采用实体建模，建立法兰接管实体模型，并划分网格.接着采用节点耦合的方法，将裂纹前沿最外圈节点与整体模型裂纹前沿连接处的节点进行耦合，得到含裂纹法兰接管实体模型的有限元模型，如图7所示.图7裂纹整体有限元模型
Fig.7The finite element model of crack1.3施加边界条件在上述建立的12模型断裂模型中施加位移边界条件和力边界条件后，就可以计算出裂纹前沿各点的应力强度因子.力边界条件：在法兰与接管内表面施加内压P内=1.1 MPa ,同时在直接管的右端面，沿接管轴向即x轴正方向，施加由内压所产生的轴向载荷P轴=9.167 MPa.位移边界条件：在法兰的下端面对其位移施加全约束，即对其施加x、y、z三个方向的位移约束.具体边界条件施加情况，如图8所示.图8接管法兰边界条件
Fig.9The boundary conditions of nozzle flange1.4计算应力强度因子施加边界条件求解后得到裂纹前沿各处的应力强度因子，裂纹前沿共有10个节点，如图9所示.图9有限元模型中裂纹前沿节点示意图
Fig.9The nodes at the front of crack in
the finite model依次求得各处节点Ⅰ型裂纹的应力强度因子KI，总共计算了7组不同裂纹角度的数据.结果如表1所示.
表1裂纹前沿各节点应力强度因子KI
Table 1Stress intensity factor at the front of crack, KI
MPa·mm1/2
节点号α/(°)454035302520151140.13132.88125.47117.64109.1399.3487.622147.03138.62130.12121.28111.84101.2388.783145.20136.25127.27118.05108.3497.6285.294142.35132.96123.60114.09104.2193.4981.405138.97129.17119.48109.7399.7289.0677.276135.20125.03115.04105.0794.9884.4172.977131.10120.58110.31100.1690.0179.5768.528126.59115.76105.2494.9484.7874.5263.939121.50110.4299.7289.3479.2369.2359.1910109.5098.8488.6078.7769.3460.2351.41从表1中可知，Ⅰ型裂纹的应力强度因子KI的最大值在裂纹前沿靠近法兰接管的外表面处.KI结果分布如图10所示.图10裂纹前沿各节点应力强度因子KI
Fig.10Stress intensity factor at the front of crack ,KI
从图10中可以看出，裂纹前沿的KI值基本呈倒“U型”分布，内表面处的应力强度因子值要比外表面的小得多；应力强度因子的最大值在裂纹前沿靠近法兰接管的外表面处.根据上面的结果可以推测，当裂纹要继续扩展时，裂纹前沿靠近法兰接管外表面处的KI最先达到材料的断裂韧度KIC，即裂纹在靠近法兰接管外表面处较先扩展.应力强度因子是判断裂纹是否扩展的主要参量，而在目前又无法得到复杂裂纹前沿应力强度因子解析解的情况下，通过数值模拟的方法求出复杂裂纹前沿各点应力强度因子的近似解，再结合平面应变断裂韧度KIC的变化规律及裂纹扩展判据KI≥KIC，就能定量的判断裂纹扩展情况.该结论对含缺陷压力容器及管道的安全评定具有重要的数据价值.2结语a. 通过节点建模法及在裂纹前沿设置应力奇异单元建立了含穿透裂纹法兰接管的有限元分析模型.b. 通过有限元计算得到了不同裂纹长度的裂纹前沿各处的应力强度因子KI的数值解，根据计算得到的数据确定了裂纹前沿应力强度因子的分布情况并确定了最大应力强度因子的位置.c. 上述计算得到的含穿透裂纹法兰接管的应力强度因子数值解可用于含缺陷压力容器及管道的安全评定.