Root crack propagation angle of load-carrying fillet weld T-joints
-
摘要: 针对Q355B低合金钢T形接头承载角焊缝根部疲劳失效的扩展路径问题,提出了一种基于受力分析计算的等效应力强度因子法(KEQ法)预测根部裂纹的扩展角度. 经有限元模拟验证,其求解应力强度因子的最大误差小于5%. 与基于有限元分析计算的最大周向应力法(MCS法)和有效结构应力法(ETS法)相比较,并结合3种不同应力水平下弯曲疲劳试验结果发现,基于等效应力强度因子法、最大周向应力法与有效结构应力法求解的裂纹扩展角度分别为25.6°,25.9°和32.2°,相较于实际疲劳试验的根部裂纹扩展角度24°,误差分别为6.67%,7.92%和34.17%. 结果表明,基于KEQ法求解裂纹扩展角度准确度最高,更适于预测承受弯曲疲劳载荷下T形接头角焊缝根部裂纹的扩展角度.Abstract: Aimed at the root fatigue failure propagation path problem of T-joints in Q355B low alloy steel, an equivalent stress intensity factor (KEQ) method based on force analysis and calculation was proposed to predict the root crack propagation angle. The maximum error of solving the stress intensity factor is less than 5% verified by finite element simulation. Compared with the maximum circumferential stress method (MCS method) and the effective structural stress method (ETS method) based on finite element analysis, and combined with the bending fatigue test results of three different stress levels, it is found that: The crack propagation angles calculated by equivalent stress intensity factor method, maximum circumferential stress method and effective structural stress method are 25.6°, 25.9° and 32.2°, respectively. Compared with the root crack propagation angle of 24° in the actual fatigue test, the errors are 6.67%, 7.92% and 34.17%, respectively. The results show that KEQ method has the highest accuracy in solving crack propagation angle, which is more suitable for predicting the root crack propagation angle of T-joint fillet weld under bending fatigue loading.
-
0. 序言
Q355B低合金钢凭借其强度高、塑性和韧性好等特点,被广泛应用于桥梁、船舶、轨道列车转向架等关键结构[1]. 在这些结构中,主要承载结构部件均经焊接而成,并在服役过程中承受着复杂的载荷条件. 而焊接接头在交变载荷条件下往往发生疲劳失效,继而成为整个结构中最薄弱的环节[2]. 因此,评估焊接结构的疲劳寿命至关重要,而焊接接头的应力状态以及相应的裂纹扩展路径更是成为了学者们关注的焦点. 评估焊接接头疲劳寿命最常用的方法有名义应力法、结构应力法等,这些方法对焊根失效问题均存在一定的局限性,准确评估疲劳寿命困难较大. 这是由于当焊根处起裂失效时,根部存在严重的应力集中,裂纹的扩展路径为焊缝区,疲劳寿命主要为焊缝区裂纹的扩展寿命,且扩展角度不同,扩展寿命不同. 当焊根处的宏观裂纹被监测到时,剩余疲劳寿命往往已不足20%[3-4],因此,准确获得根部裂纹扩展角度是评估焊根失效疲劳寿命的关键,也是结构设计中最重要的环节.
针对以上问题,Wang等人[5]采用最大周向应力准则(MCS法)模拟裂纹扩展,即裂纹垂直于最大周向应力方向从裂纹尖端开始沿径向扩展,故可以将周向应力(σθθ)达到最大时的角度作为裂纹扩展的角度. 王苹等人[3]在Xing等人[6]的基础上,针对十字接头承载角焊缝提出了有效结构应力法(ETS法). 该方法认为不同的扩展角度对应于不同的有效结构应力,其中最大有效结构应力所对应的角度即为裂纹扩展的角度. 上述两种方法均基于有限元数值模拟计算,其中MCS法具有网格敏感性,不同的网格密度有不同的结果,获得结果的精度依赖于划分网格的精密程度;相较之下,ETS法具有网格不敏感性的优点.
针对T形接头承载角焊缝根部疲劳失效问题,提出了一种新的预测裂纹扩展角度的方法——等效应力强度因子法(KEQ法). 该方法只需简单的受力分析计算即可得到根部裂纹扩展的角度,无须进行有限元模拟,可有效地避免建模时网格精度造成的误差. 为验证KEQ法的有效性和精准性,以Q355B低合金钢液压挖掘机动臂和斗杆结构中的T形接头为研究对象,从实际受力情况出发,进行了弯曲疲劳试验,测量了此接头焊根处裂纹的扩展角度. 建立了有限元模型,对KEQ法进行了验证,并与MCS法、ETS法进行了对比分析. 该研究为求解应力强度因子提供了简便途径,同时为预测T形接头角焊缝根部裂纹扩展角度问题提供了新的方法.
1. 试验方法
1.1 试验材料
文中选用工程机械领域常用的Q355B低合金钢,采用CO2气体保护焊制备T形接头,所用焊丝为ER50-6 实芯焊丝,直径为1.6 mm. 焊接工艺参数如表1所示.母材和焊丝的化学成分如表2所 示.T形接头疲劳试件取样及接头具体尺寸如图1所示. 疲劳试件经线切割、铣削加工,无加工痕迹残留,以确保疲劳试件与实际失效情况相符,由焊缝根部起裂.
表 1 焊接工艺参数Table 1. Welding processing parameters电弧电压
U/V焊接电流
I/A气体流量
Q/(L·min−1)焊接速度
v/(mm·s−1)22 200 20 5 表 2 母材及焊丝的化学成分(质量分数,%)Table 2. Chemical compositions of base metal and welding wire材料 C Si Mn S P Fe Q355B 0.20 0.25 0.60 0.04 0.02 余量 ER50-6 0.12 1.05 1.60 0.03 0.01 余量 1.2 疲劳试验
弯曲疲劳试验采用GPS300型高频疲劳试验机,其最大载荷为300 kN,试件装夹如图2所示. 弯曲疲劳试验应力比R设定为0.5,选用3个应力范围Δσ进行测试,分别为80,115 和150 MPa,每个应力范围测试6个试件.
1.3 焊接接头根部裂纹扩展角度的测量
疲劳试件失效断裂后,垂直于焊缝切取55 mm × 50 mm × 20 mm的金相试样,使用砂纸依次打磨至1000号后机械抛光至表面光亮,再使用4%硝酸+96%酒精(体积分数)的试剂进行腐蚀. 处理后的试样用于焊根处剖面的观察以及裂纹扩展角度的测量.
2. 等效应力强度因子法
裂纹尖端应力场存在奇异性,具有无法直接通过应力来反映裂纹尖端的受力特性[7]. 断裂力学理论认为应力强度因子(K)反映了裂纹尖端附近应力强弱程度,裂纹会沿着K最大的方向扩展,可通过计算最大应力强度因子间接得到裂纹扩展角度.
计算K有多种方法,如:借助数学分析的复变函数法及权函数法[8],近似计算的边界配置法[9]及有限元法,利用试验标定的柔度标定法[10]等.
对于从根部起裂的T形接头,应力强度因子手册中没有直接的公式计算K,根据BS 7910:2019 Guide to methods for assessing the acceptability of flaws in metallic structures可采用有限长板条公式进行推导,经有限元模拟验证,使用此方法计算出的解析解误差较小.
图3为有限长板条受力示意图. 根据应力强度因子手册,宽度为b的有限长板条,有一长度为a的单边裂纹,受单向均匀拉伸作用,如图3a所示,其应力强度因子KⅠ为
$$ {K}_{{\rm{I}}}={F}_{1}{\sigma }_{1}\sqrt{{\text{π}} a} $$ (1) 式中:σ1为板所受拉应力;F1为有关a和b的函数,即
$$ \begin{split} {F_1} =& 1.12 - 0.231\left( {\frac{\,a\,}{\,b\,}} \right) + 10.55{\left( {\frac{\,a\,}{\,b\,}} \right)^2}-\\ &21.72{\left( {\frac{\,a\,}{\,b\,}} \right)^3}{\text{ + }}30.39{\left( {\frac{\,a\,}{\,b\,}} \right)^4} \end{split} $$ (2) 作用在裂纹长度为a、宽度为b的有限长板条的单位厚度弯曲力矩M,如图3b所示,其应力强度因子KⅠ为
$$ {K}_{{\rm{I}}}={F}_{2}{\sigma }_{2}\sqrt{{\text{π}} a} $$ (3) $$ {\sigma }_{2}\text=\frac{6M}{{b}^{2}} $$ (4) 式中:σ2为板所受弯曲应力;F2为有关a和b的函数,即
$$ \begin{split} {F_2} = &1.122 - 1.40\left( {\frac{\,a\,}{\,b\,}} \right) + 7.33{\left( {\frac{\,a\,}{\,b\,}} \right)^2}-\\ &13.08{\left( {\frac{\,a\,}{\,b\,}} \right)^3}{\text{ + }}14.0{\left( {\frac{\,a\,}{\,b\,}} \right)^4} \end{split}$$ (5) 在含裂纹长度为a的有限长板条裂纹嘴上作用一对单位厚度的集中切向力q,如图3c所示,其应力强度因子KⅡ为
$$ {K}_{{\text{Ⅱ}}}={F}_{3}\frac{2}{\sqrt{{\text{π}} a}}q $$ (6) 式中:F3为有关a和b的函数,即
$$ {F_3} = \dfrac{{1.3 - 0.65\left( {\dfrac{\,a\,}{\,b\,}} \right) + 0.37{{\left( {\dfrac{\,a\,}{\,b\,}} \right)}^2} - 0.28{{\left( {\dfrac{\,a\,}{\,b\,}} \right)}^3}}}{{\sqrt {1 - \dfrac{\,a\,}{\,b\,}} }} $$ (7) 许多承受单轴载荷的结构部件包含随机定向的缺陷或裂纹,这些缺陷/裂纹经历混合模式的应力状态. 混合模式疲劳裂纹的一个特征是它们通常以非自相似的方式扩展[11],即当受到混合模式加载条件时,裂纹会改变扩展方向. 对于研究中承受弯曲载荷F的T形接头,焊根处未焊透可被视为初始裂纹缺陷,当承受弯曲疲劳载荷时,裂纹将沿着一定的角度进行扩展,定义扩展角度为θ. 因采用有限长板条公式进行推导,故规定等效平板的上下平面垂直于裂纹扩展面. 根据力的平移定理以及参考应力强度因子手册,将外力F平移至等效板边会产生附加弯矩M1, 将外力F分解为轴向力Fsinθ,切向力Fcosθ,并将切向力Fcosθ平移至裂纹尖端产生附加弯矩M2,受力分析如图4所示.
$$ {\sigma _1}{\text{ = }}\frac{{F\sin \theta }}{{w{t'}}} $$ (8) $$ Q_1 = \frac{{F\cos \theta }}{{{w}}} $$ (9) $$ M = {M_1} + {M_2} = F{{\text{d}}_1} + Fd\cos \theta $$ (10) $$ {\sigma _2}{\text{ = }}\frac{{6M}}{{{{w}}{{{t}}}^{'2}}} $$ (11) 式中:σ1为轴向力F sinθ产生的轴向应力; σ2为弯矩M产生的弯曲应力;Q1为单位厚度下的剪切力;d为力作用的距离,140 mm;w为取样长度,200 mm; d1为F移动的水平距离;t′为等效的有限长板条宽度;d1与t′的大小都与θ有关.
焊缝尺寸放大图如图5所示,将未焊透长度a视为初始裂纹的长度,a0为裂纹面上的等效初始裂纹长度. t′由等效初始裂纹长度a0与焊喉长度L两部分组成.
$$ {a_0} = a\cos \theta = ({t_2} - p)\cos \theta $$ (12) $$ L = \frac{{\left( {p + {{{s}}_1}} \right)\sin \alpha }}{{\sin \left( {\alpha + \theta } \right)}} $$ (13) $$ {t}'={a}_{0}+L $$ (14) $$ \alpha {\text{ = }}\arctan \left( {\frac{{{{{s}}_2}}}{{{s_1}}}} \right) $$ (15) 式中:α为焊脚角度.将上述等效关系,如图6所示,带入式(1) ~ 式(15),求出KⅠ, KⅡ.
KEQ法包含Ⅰ型与Ⅱ型两种应力强度因子,在工程结构中,Ⅰ型裂纹是引起构件疲劳失效的主导因素,Ⅰ型应力强度因子KⅠ对裂纹扩展的影响要大于Ⅱ型应力强度因子KⅡ,所以要考虑等效应力强度因子对裂纹扩展的影响,Zhang等人[12]提出的等效应力强度因子为
$$ {K}_{\text{eq}}={\left({K}_{{\text{Ⅰ}}}^{4}+8{K}_{{\text{Ⅱ}}}^{4}\right)}^{0.25} $$ (16) KEQ法中包含扩展角度θ,因此求出0° ~ 90° 中Keq的最大值,即可获得裂纹扩展角度.
3. 结果分析与讨论
3.1 KEQ法有限元模拟验证
通过式(16)计算出KEQ最大值为25.6°. 为验证KEQ法计算等效应力强度因子的准确性,对T形接头采用ABAQUS进行有限元模拟. 模型中单元为二次单元,采用二维平面应变分析,裂纹扩展角度设置为25.6°,应力范围设置为80 MPa,如图7所示. 建立不同裂纹长度的模型,分别计算Keq,并与由KEQ 法求得的Keq进行对比分析,如表3所示. 从表3可以看出,采用KEQ法求出的应力强度因子与有限元模拟的最大误差小于4%,验证了KEQ法等效过程的正确性.
表 3 不同裂纹长度下Keq对比分析Table 3. Analysis of Keq under different crack lengths等效初始
裂纹长度
a/mm应力强度因子
K/(MPa·mm−1/2)误差e(%) KEQ 法 模拟验证 a0 + 0.1 193.09 194.48 −0.93 a0 + 0.2 195.27 199.48 −2.38 a0 + 0.5 201.81 197.36 2.00 a0 + 1 212.72 207.26 2.43 a0 + 2 234.87 226.18 3.74 a0 + 4 282.54 283.16 −0.25 a0 + 6 339.12 337.24 0.55 a0 + 10 511.88 515.18 −0.64 为了进一步验证KEQ法的准确性,分别建立不同的未熔透深度a以及不同扩展角度θ的模型,计算Keq,并与KEQ法求得的Keq进行对比分析,如表4所示. 从表4可以看出,采用KEQ法求出的Keq与有限元模拟的最大误差小于5%,进一步验证了KEQ法等效过程的正确性.
表 4 a0 + 4 mm裂纹长度下Keq对比分析Table 4. Analysis of Keq at crack length of a0 + 4 mm未焊透长度
a/mm裂纹扩展
角度θ/(°)应力强度因子
K/(MPa·mm−1/2)误差e(%) KEQ 法 模拟验证 0.1 25.6 174.83 181.74 3.95 1 25.6 193.50 199.70 3.20 3 25.6 239.62 229.72 4.31 5 25.6 282.54 283.16 −0.25 8 25.6 361.90 360.81 −0.30 5 10 256.92 270.16 −4.90 5 20 282.87 278.35 −1.60 5 60 207.11 205.10 −0.97 3.2 疲劳试验及评估结果
图8为不同应力水平下裂纹扩展角度测量结果.从图8可知,Δσ = 80 MPa时,扩展角度为24°;Δσ = 115 MPa时为23.6°;Δσ = 150 MPa时为24.2°. 3个应力水平下裂纹扩展角度相差较小,均约为24°.
采用上文所述的3种方法分别计算裂纹扩展角度,在弹性范围内,其不受应力水平大小的影响,裂纹扩展角度均相同. 图9为Δσ = 80 MPa时不同测试方法下裂纹扩展角度曲线. 从图9可以看出,3种方法均呈现先上升后下降的趋势,即分别存在确定的角度使得Keq,σθθ,RETS达到最大值,对比结果如表5所示. KEQ法计算出Keq的最大值对应的角度为25.6°,采用最大周向应力法计算出的裂纹扩展角度为25.9°,采用有效结构应力法计算的裂纹扩展角度为32.2°,与实际裂纹扩展角度对比,KEQ法与MCS法误差都较小,但MCS法需要进行有限元模拟,具有网格敏感性,需要划分较为精密的网格. 因ETS法误差较大,故其适用于十字接头承载角焊缝,对于承受弯曲载荷的T形接头适用性具有局限性. 综上所述,KEQ法在无需进行有限元模拟的基础上仍具有较高的精度,适合预测承受弯曲疲劳载荷下的T形接头焊根裂纹扩展角度.
表 5 不同测试方法下裂纹扩展角度的对比Table 5. Comparison of crack propagation angle under different test methods测试方法 裂纹扩展角度θ/(°) 误差e(%) KEQ法 25.6 6.67 MCS法 25.9 7.92 ETS法 32.2 34.17 试验测试 24.0 — 4. 结论
(1)分别建立不同裂纹长度、未熔透深度、裂纹扩展角度的有限元模型计算Keq,与基于受力分析计算的KEQ法进行对比验证,最大误差小于5%.
(2) KEQ法预测T形接头承受弯曲载荷时角焊缝根部失效裂纹扩展角度,计算结果为25.6°,与疲劳试验结果相比较,误差为6.67%.
(3) MCS法与ETS法求解出的裂纹扩展角度分别为25.9°和32.2°,相较于疲劳试验结果,误差为7.92%,34.17%.
(4)经比较,KEQ法预测T形接头承载角焊缝根部裂纹扩展角度因只需要简单的力学分析,无需进行有限元模拟,最为简便,且误差较小.
-
表 1 焊接工艺参数
Table 1 Welding processing parameters
电弧电压
U/V焊接电流
I/A气体流量
Q/(L·min−1)焊接速度
v/(mm·s−1)22 200 20 5 表 2 母材及焊丝的化学成分(质量分数,%)
Table 2 Chemical compositions of base metal and welding wire
材料 C Si Mn S P Fe Q355B 0.20 0.25 0.60 0.04 0.02 余量 ER50-6 0.12 1.05 1.60 0.03 0.01 余量 表 3 不同裂纹长度下Keq对比分析
Table 3 Analysis of Keq under different crack lengths
等效初始
裂纹长度
a/mm应力强度因子
K/(MPa·mm−1/2)误差e(%) KEQ 法 模拟验证 a0 + 0.1 193.09 194.48 −0.93 a0 + 0.2 195.27 199.48 −2.38 a0 + 0.5 201.81 197.36 2.00 a0 + 1 212.72 207.26 2.43 a0 + 2 234.87 226.18 3.74 a0 + 4 282.54 283.16 −0.25 a0 + 6 339.12 337.24 0.55 a0 + 10 511.88 515.18 −0.64 表 4 a0 + 4 mm裂纹长度下Keq对比分析
Table 4 Analysis of Keq at crack length of a0 + 4 mm
未焊透长度
a/mm裂纹扩展
角度θ/(°)应力强度因子
K/(MPa·mm−1/2)误差e(%) KEQ 法 模拟验证 0.1 25.6 174.83 181.74 3.95 1 25.6 193.50 199.70 3.20 3 25.6 239.62 229.72 4.31 5 25.6 282.54 283.16 −0.25 8 25.6 361.90 360.81 −0.30 5 10 256.92 270.16 −4.90 5 20 282.87 278.35 −1.60 5 60 207.11 205.10 −0.97 表 5 不同测试方法下裂纹扩展角度的对比
Table 5 Comparison of crack propagation angle under different test methods
测试方法 裂纹扩展角度θ/(°) 误差e(%) KEQ法 25.6 6.67 MCS法 25.9 7.92 ETS法 32.2 34.17 试验测试 24.0 — -
[1] 严春妍, 易思, 张浩, 等. S355钢激光-MIG复合焊接头显微组织和残余应力[J]. 焊接学报, 2020, 41(6): 12 − 18. doi: 10.12073/j.hjxb.20191014001 Yan Chunyan, Yi Si, Zhang Hao, et al. Microstructure and residual stress of S355 steel composite welded by laser and MIG[J]. Transactions of the China Welding Institution, 2020, 41(6): 12 − 18. doi: 10.12073/j.hjxb.20191014001
[2] Wang Ping, Pei Xianjun, Dong Pingsha, et al. Analysis of weld root fatigue cracking in load-carrying high-strength aluminum alloy cruciform joints[J]. International Journal of Fatigue, 2020, 139: 105735. doi: 10.1016/j.ijfatigue.2020.105735
[3] 王苹, 米莉艳, 于谊飞, 等. 7N01铝合金十字接头抗疲劳设计[J]. 焊接学报, 2019, 40(10): 20 − 24. Wang Ping, Mi Liyan, Yu Yifei, et al. Anti-fatigue design of 7N01 aluminum alloy cross joint[J]. Transactions of the China Welding Institution, 2019, 40(10): 20 − 24.
[4] 刘永, 王苹, 马然, 等. 铝合金非承载十字接头疲劳特性[J]. 焊接学报, 2016, 37(8): 83 − 87. Liu Yong, Wang Ping, Ma Ran, et al. Fatigue property of aluminum non-load bearing cruciform joint[J]. Transactions of the China Welding Institution, 2016, 37(8): 83 − 87.
[5] Wang Cungen, Wang Shuhong, Chen Guangqi, et al. Implementation of a J-integral based maximum circumferential tensile stress theory in DDA for simulating crack propagation[J]. Engineering Fracture Mechanics, 2021, 246: 107621. doi: 10.1016/j.engfracmech.2021.107621
[6] Xing Shizhu, Dong Pingsha, Alina Threstha. Analysis of fatigue failure mode transition in load-carrying fillet-welded connections[J]. Marine Structures, 2016, 46: 102 − 126. doi: 10.1016/j.marstruc.2016.01.001
[7] Seung Jae Kim, Kyung Dong Bae, Yun Jae Kim, et al. The effect of complex crack shape on crack driving forces and crack tip stress fields: Experiment and FE analysis[J]. International Journal of Pressure Vessels and Piping, 2020, 185: 104135. doi: 10.1016/j.ijpvp.2020.104135
[8] Kou K P, Cao J L, Yang Y, et al. Weight function method for stress intensity factors of semi-elliptical surface cracks on functionally graded plates subjected to non-uniform stresses[J]. Materials, 2020, 13(14): 3155. doi: 10.3390/ma13143155
[9] Cheng Jiaxing, Sun Bin, Wang Mengyuan, et al. Analysis of III crack in a finite plate of functionally graded piezoelectric/piezomagnetic materials using boundary collocation method[J]. Archive of Applied Mechanics, 2019, 89(2): 231 − 243. doi: 10.1007/s00419-018-1462-y
[10] Alireza Akhavan Safar, Mazaher Salamat talab, Arsalan Ajdani, et al. Mode II fracture energy characterization of brittle adhesives using compliance calibration method[J]. Fatigue & Fracture of Engineering Materials & Structures, 2020, 43(9): 1928 − 1937.
[11] Silva A L L, Jesus A M P de, Xavier J, et al. Combined analytical-numerical methodologies for the evaluation of mixed-mode (I+II) fatigue crack growth rates in structural steels[J]. Engineering Fracture Mechanics, 2017, 185: 124 − 138. doi: 10.1016/j.engfracmech.2017.04.016
[12] Zhang Hui, Ali Fatemi. Short fatigue crack growth behavior under mixed-mode loading[J]. International Journal of Fracture, 2010, 165(1): 1 − 19. doi: 10.1007/s10704-010-9497-2
-
期刊类型引用(3)
1. 张明军,李晨希,邹江林,程波,张健,仝永刚,胡永乐,陈根余. AZ31B镁合金功率调制环形光斑光纤激光焊接试验研究. 机械工程学报. 2025(02): 151-161 . 百度学术
2. 刘坤,李洁,王浩,简思捷. 镁合金焊接凝固裂纹敏感性评价及晶间液相回填规律分析. 焊接学报. 2023(09): 9-15+129 . 本站查看
3. 焦婧,黄金鑫,张志凯. 论带复杂油路类镁合金铸件的清理方法. 世界有色金属. 2022(23): 175-177 . 百度学术
其他类型引用(2)