焊接工艺因其经济便利等特点被广泛应用于各类工程结构的连接[1-2]。然而,焊接结构往往存在初始缺陷,导致疲劳问题显著[3]。目前,针对焊接节点疲劳寿命的研究方法主要有两种:一种是基于S-N曲线的传统疲劳分析方法[4-5];另一种是基于裂纹扩展速率公式的断裂力学方法[6-7]。扩展有限元方法在断裂力学方法的基础上,凭借无需进行网格重划分和预测裂纹扩展路径的优点,成为疲劳分析的主要方法之一。
由于焊接过程中加热和冷却过程不均匀,焊接接头会产生残余应力和不均匀的微观结构[8]。而焊接残余应力往往是导致焊接结构疲劳破坏的主要原因[9-10],因此,对残余应力的研究一直是热点。Lee等[11]建立了基于循环塑性本构方程的三维热弹塑性有限元模型,评估循环荷载下的应力重分布情况,预测残余应力的周期性松弛。王春生等[12]在断裂力学数值模型中引入焊接残余应力场,分析其对正交异性钢桥面板疲劳裂纹静、动态扩展行为的影响。Jiang等[13]基于对残余应力的试验测量,分析了循环荷载下残余应力重分布的机理和对疲劳性能的影响,提出了焊接接头疲劳预测模型。实际工程中,为提高焊接结构的疲劳性能,往往会对焊接结构进行消除残余应力处理。焊后热处理工艺作为一种有效消除残余应力的方法被广泛应用于各类焊接结构[14]。许多学者对焊后热处理工艺消除残余应力开展了数值计算和试验研究。Venkata等[15]以P91钢对接接头为研究对象,建立数值模型模拟了不同热处理参数对对接接头残余应力的消除情况并进行试验验证。Mitra等[16]通过有限元和试验方法确定了温度和时间对应力释放的影响,并建立了热处理参数与残余应力的关系函数,以达到不同温度和时间组合实现相同的应力释放。逯世杰等[17]模拟了焊接过程及焊后热处理的应力场,并研究了蠕变效应在焊接过程中的作用,分析了焊后热处理消除残余应力的机理。Pokorny′等[18]提出了在外部荷载与残余应力条件下车轴的疲劳寿命评估方法,并分析了热处理对结构的影响。现有对热处理后疲劳性能的研究多以圆形结构为主,T型节点结构相对复杂,热处理操作困难,对T型节点的热处理和疲劳性能研究相对较少。因此,有必要建立T型焊接结构疲劳分析模型,研究热处理工艺对焊接结构疲劳性能的影响。
采用直接耦合热力分析法,建立T型焊接节点的热弹塑性数值模型,进行T型焊接节点的残余应力分析;用热处理方法进行消除残余应力计算;基于线弹性断裂力学理论和扩展有限元方法,建立疲劳分析有限元模型并引入残余应力场,分别计算在不考虑残余应力、考虑残余应力及热处理消除残余应力情况下的疲劳寿命。
参考文献[19]的试验数据建立T型焊接接头疲劳有限元模型。T型接头由一块400 mm×100 mm×16 mm和一块200 mm×100 mm×15 mm钢板焊接而成,焊缝为全熔透焊缝。T型焊接节点尺寸如图 1所示,焊缝简化为45°坡脚形式,建模为直角三棱柱。
利用直接耦合的方式模拟焊接行为,即采用温度-位移耦合分析步,单元类型采用八节点热耦合六面体单元(C3D8T)。热分析采用双椭球热源进行加热,热源模型如图 2所示。对应的热源模型数学公式表示为
式中:q为热流密度, J/(m2·s);f为热源移动函数;Q为单位时间热量, J;a、b、c为热源形状参数,取值分别为6、6、2.8、2.8 mm。
式中:η为焊接热效率;U为焊接电压,V;I为焊接电流,A。焊接电压取280 V,电流取25 A,热效率取0.8。热源焊接速度根据实际情况取为10 mm/s。采用dflux子程序控制热源在有限元模型上的移动。为简化计算,假定焊缝金属与母材具有相同的热物理和力学性能参数。焊接分析中钢材热物理与热力学参数主要参考文献[17],如图 3和图 4所示。
图 5左侧图片为第1道焊缝焊接完成时焊缝区域的温度场,右侧为试验的熔池形貌。假设以图中1 420 ℃为熔点,焊缝区域和少部分母材处于熔化状态,对比文献[19]试验方法得到的焊接热影响区的熔池形貌,红色区域与试验结果较为接近。由此可推断,焊接有限元温度场与实际焊接结果大致相同,间接验证了有限元方法的合理性。
以焊接温度场计算结果作为热荷载条件进行焊接应力计算。焊接时,由于焊接过程中的局部和瞬时加热以及整体冷却导致不均等的热传导和膨胀,冷却后接头局部区域发生弹塑性变形,产生较大的残余应力。图 6中纵向残余应力与文献[17]的数值模拟和实测结果较为接近,说明残余应力数值模拟较为合理。
图 7(a)显示了焊接模拟获得的横向焊接残余应力分布。焊缝起始和终止位置处于残余压应力状态,最大峰值压应力为130 MPa;焊缝中间区域处于残余拉应力状态,最大峰值拉应力为140.9 MPa。右方焊接残余应力分布范围略小于左方,这主要是因为,虽然焊接边界条件相同,但是由于焊接起始位置位于右方,导致约束条件和散热条件的不同,从而导致横向残余应力分布略有差异。
图 7(b)显示了焊接模拟获得的纵向焊接残余应力分布。焊缝起始和终止位置处于残余压应力状态,最大峰值压应力为91.3 MPa;焊缝中间区域处于残余拉应力状态,最大峰值拉应力为292.5 MPa。
以计算得到的焊接残余应力作为初始条件进行焊后热处理过程的模拟计算。焊接并冷却至室温后,以20 ℃/min的速度将焊接节点加热至600 ℃,再保温30 min,最后进行炉冷,热处理曲线如图 8所示。应力的释放主要通过两个机制,首先,由于计算过程使用了随温度变化的材料性能参数,故随着焊后热处理过程中温度的升高,材料的强度逐渐降低;其次,在热处理阶段考虑材料蠕变引起的塑性变形会导致应力松弛,且此阶段的蠕变服从Field & Field蠕变方程
式中:εcr为等效蠕变应变;t为时间,s;σ为等效应力;a、b、c取值参考文献[17],a=-457+3.9T-0.01T2、b=11-0.09T、c=289-22.4T+0.006T2,T为与温度有关的常数。
为定量研究经过焊接后热处理过程(PWHT)焊接残余应力的消除情况,提取热处理前后沿焊缝轴线上的横向和纵向残余应力进行对比分析。图 9中残余应力与文献[17]的数值模拟和实测结果较为接近,说明热处理数值模拟较为合理。对比图 7与图 9,残余应力分布形式没有发生变化。经过焊接后热处理,拉伸残余应力的峰值大幅下降,但压缩残余应力的峰值变化不大,说明焊后热处理可以很好地消除由焊接产生的拉伸残余应力。横向残余应力的峰值由热处理前的140.9 MPa下降到59.1 MPa,下降了58.1%。纵向残余应力的峰值由热处理前的292.5 MPa下降到98.5 MPa,下降了66.3%。
Belytschko等[20]于1999年首次提出扩展有限元(XFEM)理论。该理论克服了有限元裂纹扩展对更新网格的需求。有限元位移公式引入富集函数,识别裂纹表面位移的跳跃和裂纹尖端的奇异性。
式中:u(x)为位移矢量;NI(x)为节点形函数;uI为节点位移矢量;H(x)为Heaviside函数;aI为考虑位移跳跃的增强节点自由度;Fj(x)为裂纹尖端函数;bIj为描述裂纹尖端的奇异性的节点自由度。对应的增强节点如图 10所示。
Heaviside函数为
式中:x为裂纹面任意点;x*为裂纹上最接近x的点;n为垂直于裂纹x*向外的单位矢量。
奇异函数Fj(x)为
式中:r、θ为裂纹尖端局部坐标系极坐标。
Paris和Erdogan于1963年提出用于表达裂纹扩展的Paris公式
式中:C和m为材料参数;a为裂纹长度;N为循环次数;ΔK=(Kmax+Kmin)/2,Kmax和Kmin为裂纹尖端最大和最小应力强度因子。
Paris公式可以表示为
式中:c1和c2为材料参数;ΔG为应变能释放率幅值,当裂纹尖端满足Gth < ΔG < Gpl时,裂纹开始扩展,Gth为能量释放率阈值,Gpl为能量释放率上限值。
T型节点焊接有限元模型几何参数与疲劳有限元模型相同。在存在较大空隙的情况下,空隙形状(即长宽比)对整体疲劳寿命的影响很小[21],故将初始裂纹模拟为a=1 mm的半圆形表面缺陷。通常,裂纹仅在腹板周围的两条焊缝中的一条上萌生并扩展,而另一侧的焊缝保持完整。因此,仅在其中一条焊缝上引入半圆形初始裂纹。
为节约计算量,对焊接接头进行局部加密处理。焊趾周围的单元尺寸为1 mm。远离焊趾的区域使用2、4 mm的单元尺寸(图 11),单元总数为101 512。焊接钢的Paris常数为m=2.88和C=6.77×10-13 MPa·m1/2 [22],并考虑了超强匹配接头对裂纹扩展率的影响[23]。其他材料参数如表 1所示。T型节点腹板上方施加循环荷载,荷载为35 kN,应力比R=-1,频率为3 Hz。边界条件为两端固定(图 12)。
以裂纹深度达到钢板厚度(16 mm)的一半作为疲劳失效标准[19]。图 13比较了试验结果[19]和XFEM计算结果。从图中可看到,两次试验平均失效循环次数为12 995,根据XFEM预测得出的失效循环次数为13 673(误差为+5%),XFEM得到的疲劳裂纹扩展预测结果与试验数据基本一致。试验数据和有限元结果略有差距,主要原因为:1)材料参数来源于贯穿裂纹扩展试验,而焊接对材料参数有一定影响;2)在焊趾中间采用半圆形形状的初始裂纹,实际初始裂纹可能并不对称,也可能并非从焊趾中间萌生,或者有多个初始裂纹萌生并融合成一条裂纹。
将焊接过程产生的残余应力、热处理后的残余应力作为初始条件,模拟得到不考虑残余应力、考虑焊接残余应力、热处理消除残余应力的T型焊接节点疲劳寿命如图 14所示。三者的失效循环次数分别为13 673、12 078、12 812。考虑焊接残余应力的T型焊接节点相较于不考虑残余应力的模型疲劳寿命减少了11.7%,其疲劳扩展速率即裂纹长度与循环次数的比值相较于不考虑残余应力的模型更大,说明焊接残余应力加快了疲劳裂纹的扩展。热处理后T型节点疲劳寿命相比不进行热处理的T型节点疲劳寿命增加了6.1%,说明热处理对疲劳性能影响较小。其主要原因可能是热处理不仅会降低结构的残余应力,而且对材料的性能也会产生相应的影响,材料性能的改变也会对结构的疲劳性能产生影响。残余应力的降低固然会导致疲劳强度的提高,但材料性能发生了改变,导致疲劳性能改变并不明显。裂纹扩展前期速度慢,后期速度变快,这可能是因为热处理虽然减小了焊缝中心区的高拉伸应力,前期疲劳扩展速率减小,但压缩应力也随之减小,后期疲劳裂纹更容易扩展。
图 15展示了考虑残余应力的裂纹动态扩展情况。考虑残余应力和热处理的有限元模型疲劳寿命虽无显著变化,但和不考虑残余应力的模型相比,其疲劳扩展路径有较大区别。疲劳裂纹由不考虑残余应力时的中间向两端均匀发展变为一端扩展快而另一端扩展慢。其疲劳裂纹扩展路径与试验[19]得到的裂纹断面较为接近。产生该现象的原因为焊接残余应力的分布不均匀导致疲劳裂纹一端更易发展。
考虑焊接残余应力及热处理对疲劳寿命的影响,建立T型焊接节点热弹塑性有限元模型和扩展有限元模型,分别进行残余应力、热处理和疲劳分析。研究结果表明:
1) 未考虑残余应力时,T型焊接节点虽然疲劳寿命与试验较为接近,但疲劳裂纹向焊缝两端均匀扩展,裂纹扩展形式与试验不符。考虑残余应力后,疲劳裂纹变为一端扩展快而另一端扩展慢,扩展路径与试验更为接近。
2) 经过600 ℃、30 min热处理后,T型焊接节点焊接残余应力得到充分消除,纵向残余应力峰值降幅为66.3%,横向残余应力峰值降幅为58.1%。
3) 热处理减小了焊缝中心区的高拉伸应力,前期疲劳扩展速率减小,但压缩应力也随之减小,后期疲劳裂纹更容易扩展。热处理后,T型节点疲劳寿命相比不进行热处理的T型节点疲劳寿命增加了6.1%。