高级检索

覆材厚度对不锈钢复合板焊接接头应力应变影响的数值模拟分析

冯玉兰, 吴志生, 孙智宇

冯玉兰, 吴志生, 孙智宇. 覆材厚度对不锈钢复合板焊接接头应力应变影响的数值模拟分析[J]. 焊接学报, 2024, 45(1): 73-82. DOI: 10.12073/j.hjxb.20230606001
引用本文: 冯玉兰, 吴志生, 孙智宇. 覆材厚度对不锈钢复合板焊接接头应力应变影响的数值模拟分析[J]. 焊接学报, 2024, 45(1): 73-82. DOI: 10.12073/j.hjxb.20230606001
FENG Yulan, WU Zhisheng, SUN Zhiyu. Numerical simulation of the influence of thickness of cladding material on stress and strain of welded joint of stainless steel composite plate[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2024, 45(1): 73-82. DOI: 10.12073/j.hjxb.20230606001
Citation: FENG Yulan, WU Zhisheng, SUN Zhiyu. Numerical simulation of the influence of thickness of cladding material on stress and strain of welded joint of stainless steel composite plate[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2024, 45(1): 73-82. DOI: 10.12073/j.hjxb.20230606001

覆材厚度对不锈钢复合板焊接接头应力应变影响的数值模拟分析

基金项目: 国家重点研发计划项目(2018YFA0707305);山西省重点研发计划项目(202102050201001);山西省重点研发计划项目(202101120401008);山西省青年科学基金项目(201901D211291);海安太原理工大学研究院2023开放研发项目(2023HA-TYUTKFYF012);太原科技大学研究生教育创新项目(BY 2022005)
详细信息
    作者简介:

    冯玉兰,博士研究生,高级工程师;主要研究方向焊接工艺及其自动控制;Email: 799621880@qq.com

    通讯作者:

    吴志生,博士,教授,博士研究生导师;Email: zswu1963@tyust.edu.cn

  • 中图分类号: TG 404

Numerical simulation of the influence of thickness of cladding material on stress and strain of welded joint of stainless steel composite plate

  • 摘要:

    采用钨极气体保护焊(gas tungsten arc welding,GTAW)制备了基材厚度为13 mm,覆材厚度分别为0.3,1.0,2.0 mm的304/Q355复合板焊接接头. 采用数值模拟法,分析了覆材厚度对304/Q355复合板焊接接头应力应变的影响,揭示了不同覆材厚度的不锈钢复合板焊接过程中残余应力及应变的演变特征. 结果表明,随着覆材厚度的增加,304/Q355复合板焊接接头截面区最大残余应力逐渐减小,残余应力区宽度逐渐增大,在基层焊缝向过渡层焊缝焊接过程中,应力集中区逐渐向过渡层焊缝区收缩,最高应力集中区均位于过渡层焊缝区,最大应力接近材料屈服强度. 焊后变形量随覆材厚度的增加而增大,覆材厚度2.0 mm的复合板变形量比覆材厚度0.3 mm的复合板变形量大约2 mm,变形最大的方向为焊缝厚度方向.

    Abstract:

    GTAW technology has been used to fabricate welded connections on composite plates consisting of 304 and Q355 materials. Base material thickness was 13 mm, while cladding material thicknesses ranged from 0.3 mm to 2.0 mm. The influence of different thicknesses of cladding material on the stress and strain of 304/Q355 composite plate welded joint was analyzed by numerical simulation method, and the evolution characteristics of residual stress and strain during welding of stainless steel composite plate with different thicknesses of cladding material were revealed. The results show that the X-ray diffraction test results are in good agreement with the numerical simulation results. With the increase of the thickness of the cladding material, the maximum residual stress value of the welded joint section of 304/Q355 composite plate gradually decreases, and the width of the residual stress zone gradually increases. During the welding process from the base weld to the transition weld, the stress concentration zone gradually shrinks to the transition weld zone. Finally, the highest stress concentration area is located in the weld zone of the transition layer, and the maximum stress value is close to the yield strength of the material. In addition, the deformation after welding increases with the increase of the thickness of the cladding material, and the deformation of the composite plate with the thickness of the cladding material 2.0 mm is about 2 mm higher than that of the composite plate with the thickness of the cladding material 0.3 mm, and the direction of the largest deformation is the direction of the weld thickness.

  • 随着现代工业不断推进,金属板材应用工况愈发复杂,单一金属板材已无法满足多元化的服役要求,发展趋势从单一性能向着综合性能方向发展[1-3], 双金属复合材料是通过各种复合工艺将热物理、力学等性能不同的金属结合,兼具两种金属的优异性能[4-5],既满足工业上对结构与性能需求,也可以削减材料成本,因而在汽车、石油化工、能源等领域得到广泛的使用. 不锈钢复合板具有耐腐蚀、高强度、性价比高等优点,目前被广泛应用于装备制造领域[6],焊接是重大装备的关键加工技术,因此成为金属复合材料应用的关键制造技术. 焊缝金属由基层焊缝、覆层焊缝、过渡层焊缝组成[7],3层焊缝金属间化学成分、热物理性能差异较大,易产生焊接应力与应变,应力应变是影响焊接接头质量的重要因素,焊接区域附近过大的残余应力及应变,通常是焊接接头失效的直接原因[8],而复合板覆材厚度与其焊缝成分及应力、应变有关系,因此研究复合板覆材厚度对其焊接接头应力与应变的影响,具有重要工程实际意义[9-11].

    采用304/Q355复合板为研究对象,采用钨极气体保护焊制备了基材厚度为13 mm,覆材厚度分别为0.3,1.0,2.0 mm的复合板焊接接头. 利用数值模拟与试验相结合的手段研究了复合板焊接接头表面残余应力,分析了不同覆材厚度对304/Q355复合板焊接接头应力应变的影响,揭示了不同覆材厚度的不锈钢复合板焊接过程中残余应力及应变的演变特征,研究结果为不锈钢复合板焊接接头应力应变的控制和优化提供理论参考.

    焊接过程分基层、过渡层、覆层3层焊接,基层焊缝填充材料为ER50-6,过渡层焊缝填充材料为ER309,覆层焊缝填充材料为ER308, 母材及焊接填充材料化学成分见表1,焊接工艺参数见表2. 试验选用304/Q355双金属复合板作为模型材料,基材Q355厚度均为13 mm,覆材304厚度分别为0.3,1.0,2.0 mm(A组模型,B组模型,C组模型),试板尺寸为300 mm × 200 mm,焊接坡口形式V形,如图1所示.

    表  1  母材和焊材的化学成分(质量分数,%)
    Table  1.  Chemical composition of base material and welding materials
    材料CSiMnPSNiCrMoCuN
    3040.081.002.000.0450.038.0018.000.1
    Q3550.20.551.600.0250.0250.300.300.40
    ER3080.081.000.500.0400.037.0022.000.750.75
    ER3090.151.000.50.0400.0312.0018.000.750.75
    ER50-60.1050.9750.1630.0150.0130.0160.06
    下载: 导出CSV 
    | 显示表格
    表  2  焊接工艺参数
    Table  2.  welding parameters
    焊缝道次焊接电流
    I/A
    电弧电压
    U/V
    焊接速度
    v/(mm·s−1)
    适用模型
    ① ~ ③100 ~ 1408 ~ 102A/B/C
    ④ ~ ⑤15012 ~ 152A/B/C
    ⑥ ~ ⑦1358 ~ 112C
    下载: 导出CSV 
    | 显示表格
    图  1  不同覆材厚度的坡口形式及焊道布置 (mm)
    Figure  1.  Groove form and pass arrangement. (a) 0.3 mm;(b) 1.0 mm;(c) 2.0 mm

    利用ANSYS软件,基于有限元理论建立304/Q355不锈钢复合板模型,按照上述焊接工艺进行模拟,根据焊接时所用材料和模型确定材料性能参数、热边界条件和约束条件,模拟温度场并将温度载荷通过间接耦合的方式传输到应力场,从而研究应力应变场的规律.

    建立的有限元模型如图2所示,在设计网格划分时,网格密度与计算精确度呈现正相关,网格越密集,其结果越精确, 在保证计算结果精度前提下,对远离焊缝区域网格划分相对稀疏. 利用Hypermesh软件将模型网格划分成焊缝处较密,母材处相对稀疏的形式,二者之间利用笛卡尔网格划分方式进行疏密过渡, 由于需要计算不同场的结果,不同网格适应的不同物理场自由度.在ANSYS软件中,计算温度场使用 Plane55单元或Solid70六面体单元,应力场使用Plane55和Solid45六面体单元,图2为3种模型的网格划分,上述网格划分的单元总数分别为32 200、28 840和27 580, 在3个有限元计算模型中使用了三点自由约束的形式,如图3所示.

    图  2  网格划分模型及局部放大
    Figure  2.  Meshing model and local magnificatio. (a) A; (b) B; (c) C
    图  3  约束条件
    Figure  3.  Constraints condition

    焊接过程是典型的非线性瞬态分析,在进行温度场分析需要的参数有热导率(thermal conductivity)、密度(density)、比热容(specific heat),进行应力场分析则需要的参数有屈服强度(yield strength)、热膨胀系数(thermal expansivity)、泊松比(poisson's ratio)和弹性模量(young's modulus). 由于部分材料在高温时处于熔化状态,其高温下的物理性能较难测量,对于304奥氏体不锈钢需要考虑加工硬化和退火软化[12-14]. 通过JMatProV7.0软件计算材料性能参数,304和Q355合金随温度变化的热物理性能和力学性能参数如图4所示,初始温度为25 ℃,泊松比设置为恒定值0.3,切线模量设为杨氏模量的1/10.

    图  4  304和Q355力学—热物理性能参数
    Figure  4.  Mechanical — thermophysical performance parameter of Q355 and 304. (a) mechanical 304; (b) thermal physics 304; (c) mechanical Q355; (d) thermal physics Q355

    采用Goldak双椭球体热源来模拟各单条焊道的焊接过程,热流分布包括前后两个部分

    $$ q_{\rm{f}}(x, y, z)=\dfrac{6 \sqrt{3} f_{\rm{f}} Q}{a_{\rm{f}} b c \pi \sqrt{\pi}} \exp \left(-\dfrac{3 x^2}{a_{\rm{f}}{ }^2}-\dfrac{3 y^2}{b^2}-\dfrac{3 z^2}{c^2}\right) $$ (1)
    $$ q_{\rm{r}}(x, y, z)=\dfrac{6 \sqrt{3} f_{\rm{r}} Q}{a_{\rm{r}} b c \pi \sqrt{\pi}} \exp \left(-\dfrac{3 x^2}{a_{\rm{r}}{ }^2}-\dfrac{3 y^2}{b^2}-\dfrac{3 z^2}{c^2}\right)$$ (2)

    式中:xyz为热源中心为原点的横、纵、竖坐标值;bc为形状参数;afar为前,后半轴圆长轴;Q为有效热输入;fffr分别是熔池前部和后部总输入功率的能量分布系数为1.4和0.6. 双椭球热源模型如图5所示[15].

    图  5  双椭球热源模型
    Figure  5.  Double ellipsoidal heat source model

    焊接过程中,复合板内部的热量传递以热传导为主,而复合板表面与周围环境介质的热量传递则以热辐射和热对流为主,其中,在焊缝周围热辐射损失占主要地位,而远离熔池的母材区域热对流则占主导地位[16-17].

    (1)热对流系数:指的是物体间相互接触,表面温度传递一种能力,一般而言,焊接填充材料与母材的热对流可采用考虑下式,即[18]

    $$ {q}_{\rm{c}}=-{h}_{\rm{c}}(T_{\rm{s}}-T_0) $$ (3)

    式中:hc为热对流系数为15 W/(m2·℃);Ts为复合板的表面温度;T0为环境温度为 25 ℃.

    (2)热辐射系数:高于−273.15 ℃的物体均对外产生辐射,随着温度升高辐射能量也增强,表示辐射能量的能力就是热辐射系数.

    $$ {h}_{\rm{r}}=\sigma \varepsilon {F}\left({T}^2+{T}_r^2\right)\left({T}+{T}_r\right)$$ (4)
    $$ {{{q}}_{\rm{r}}} = - \sigma \varepsilon {{F}}\left\{ {{{\left( {{{{T}}_s} + 273} \right)}^4} - {{\left( {{{{T}}_0} + 273} \right)}^4}} \right\} $$ (5)

    式中:qr为热流密度;σ为斯蒂芬—玻尔兹曼常数(5.64 × 10−14 J/mm2∙s);ε为辐射发射率 0.85;F为配置因子;hr为对流系数;Tr为辐射温度;T为瞬时温度;Ts为复合板表面温度;T0为环境温度.

    基于热弹塑性理论研究焊接应力场需遵循Von Mises屈服准则和流动准则, 将温度载荷导入应力场中,并以热应变的形式来计算应力应变,一般情况下,每个材料节点的总应变,即

    $$ \varepsilon =\varepsilon^{\text{ε}} + {\varepsilon ^{\rm{p}}} + {\varepsilon ^{{\rm{th}}}} $$ (6)

    式中:弹性应变$ \varepsilon ^\text{ε}$遵循Hookean定律;塑性应变$ \varepsilon ^{\rm{p}}$为各向同性硬化模型;εth为热应变.

    根据Hookean定律,弹塑性应力应变关系为

    $$ \{\mathrm{d} \sigma\}=[{D}]\{\mathrm{d} \varepsilon\}-\{{C}\} \mathrm{d} T$$ (7)

    式中:[D]为弹塑性矩阵;{C}为温度相关向量;{dσ}为内应力增量;{dε}为塑性应变增量;dT为单位时间.

    根据各节点和单元的状态,将集中节点力与位移载荷为总刚度矩阵与总载荷量,平衡方程式为

    $$ {K}\{\mathrm{d} \sigma\}=\{\mathrm{d}F\} $$ (8)

    式中:K为总刚度矩阵;{dF}为载荷总量;{dσ}为位移载荷矩阵.

    总载荷量表达式为

    $$ \{\mathrm{d}F\}=\sum\left(\{\mathrm{d}F\}^e + \{\mathrm{d}R\}^e\right) $$ (9)

    式中:$\{\mathrm{d}R\}^{{e}}$为等效节点力的增量,当不存在外部载荷条件时,每个单元内节点力相互平行,即$ \Sigma\{\mathrm{d}F\}^e=0$$ \{\mathrm{d}F\}=\Sigma\{\mathrm{d}R\}^e$$\{\mathrm{d}F\}^{{e}}$为等效节点载荷增量;{dF}为载荷总量.

    基于上述理论计算焊接残余应力,将模型整体划分网格,通过各节点温度,将温度载荷导入应力场, 这种方法因不考虑材料的相变应力,具有计算精度高,收敛性好的优点,既简化了模型,也可以得出相对精确的结果. 通过式(10)能够求出节点的位移增量$ \{\mathrm{d} \delta\}$,得到应变—位移方程如下

    $$ \{\mathrm{d} \varepsilon\}^{{e}}=[{B}]\{\mathrm{d} \delta\}^{{e}} $$ (10)

    式中:{dε}为等效节点塑性应变增量;[B]网格的个数矩阵;{dδ}e为节点的位移增量,把式(10)代入式(7)中计算出应力增量,来求解其应力应变过程.

    为验证所建立的有限元模型和采用的热源模型的准确性,采用覆材厚度为1.0 mm的304/Q355复合板温度场进行验证,图6为焊缝截面形貌和多焊道热源模型叠加对比,及盖面截面形貌和单道热源模型对比,可见熔池形貌与实际焊缝相近,说明经过校核的双椭球热源模型是合理准确的.

    图  6  模拟熔池与实际熔池形貌对比
    Figure  6.  Comparison of simulated melt pool with outual melt pool morphology (a) superimposed comparison of macroscopic morphology of weld section and multi-pass welding; (b) comparison between simulation results of cover welding and macroscopic morphology of weld section

    图7为覆材厚度为0.3,1.0,2.0 mm的304/Q355复合板的残余应力分布.可知,304侧横向残余应力最大值位于焊趾处,呈对称分布,覆材厚度从0.3 mm增大到2.0 mm时应力变化不大,纵向残余应力在焊缝及热影响区应力均较高. 覆材304侧横向应力远低于纵向应力,说明焊接过程,接头受到的横向拘束比纵向拘束小, 这是由于在垂直于焊缝方向上,焊缝和热影响区在加热冷却后,受到周围温度较低金属的拘束,在局部区域产生很大的横向及纵向残余应力,而在垂直于焊缝方向(纵向),由于覆材304的室温屈服强度比Q355低约131%,而其线膨胀系数却是Q355的170%,在加热过程中具有较强的塑性变形趋势,但在冷却过程会受到Q355材料的强烈约束,覆材304和基材Q355之间巨大的热物理性能的差异,使焊接加热冷却后,覆材304侧纵向残余应力较大. 而304侧横向未受到异种材料物理性能差异的影响,因此表现出的应力较低,但热影响区两侧应力较高,这是由于热影响区在加热冷却后受到周围温度较低金属的拘束作用产生. 304侧纵向残余应力随覆材厚度的增加,焊缝两端应力减小,这可能有两方面原因:①焊缝两端为自由端,在焊接冷却收缩过程中受到的约束减小;②随着覆材厚度的增加,受到基材和覆材热物理性能差异的影响减小. Q355表面最大纵向和横向残余应力均分布在焊缝及其热影响区,这是由于坡口形式为V形,坡口根部位于Q355侧,Q355侧焊缝较窄,因此焊缝残余应力宽度相对304侧较窄,焊缝横向残余应力高于纵向残余应力,说明Q355侧焊缝受到的横向拘束比纵向拘束大.

    图  7  不同覆材厚度复合板残余应力
    Figure  7.  Residual stress distribution of composite plate with different thickneos of cladding materical. (a) 0.3 mm;(b) 1.0 mm;(c) 2.0 mm

    在焊缝中间垂直焊缝方向上,作一截面观察复合板焊接过程,每一道焊缝截面的残余应力分布如图8~图10所示.覆材厚度0.3 mm的304/Q355复合板接头横截面见图8,可见随着焊接道次的增加,等效应力分布范围越来越宽, 在焊接过程中,近焊缝区由于受热影响大,产生的热应力较大, 远离焊缝的母材受热影响较小,热应力较低. 在冷却过程中,受热部位冷却收缩大,所以在焊缝及近缝区应力较大;覆材厚度1.0 mm的复合板焊接横截面见图9,其中第一道焊缝、第二道焊缝、第三道焊缝焊后最大等效应力约为340 MPa,第四道焊缝、第五道焊缝焊后最大等效应力约为355 MPa,与材料屈服强度相接近,同时,不同道次的焊缝横截面应力最大均集中在焊缝中间部位;覆材厚度2.0 mm的304/Q355不锈钢复合板焊缝横截面应力见图10,由于304和Q355两种材料的物理和化学性能不同,在焊接加热的过程中,其冷却和受热速率不同,在两种材料的交界处存在着应力不连续的情况,可以看到,每道焊缝横截面的最大应力均在材料的屈服强度附近. 总体观察,从基层焊缝到覆层焊缝焊接过程,焊缝截面应力集中区逐渐向过渡层焊缝区收缩.

    图  8  覆材厚度0.3 mm的304/Q355不锈钢复合板焊缝横截面应力云图(MPa)
    Figure  8.  Stress cloud diagram of weld section of 304/Q355 stainless steel composite plate with 0.3 mm thickness ofcladding material
    图  9  覆材厚度1.0 mm的304/Q355不锈钢复合板焊缝横截面应力云图(MPa)
    Figure  9.  Stress cloud diagram of weld section of 304/Q355 stainless steel composite plate with 1.0 mm thickness ofcladding material
    图  10  覆材厚度2.0 mm的304/Q355不锈钢复合板焊缝横截面应力云图(MPa)
    Figure  10.  Stress cloud diagram of weld section of 304/Q355 stainless steel composite plate with 2.0 mm thickness ofcladding material

    此外,观察复合板焊接最终状态的截面残余应力分布特征可以发现,3种不同覆材厚度的复合板焊接接头的残余应力分布特征类似,最大残余应力均集中在过渡层焊缝内,这与其他研究报道的坡口形式对双金属复合板多层多道焊接头残余应力分布特征一致[19-22],主要是由于碳钢和不锈钢热膨胀系数差别大,304不锈钢热膨胀系数是Q355低合金高强钢的170%,同时由于多层多道焊,中间部位的温度传递主要通过金属间的热传导,能通过空气对流等方式的区域少,由于这些原因,该区域等效应力较高,而高应力集中对于复合板的安全服役存在巨大的隐患,另外观察最后一道焊缝焊接完成后,全厚度焊缝截面的应力分布,可知随着覆材厚度的增加,复合材料界面位置的应力逐渐降低.

    焊接残余变形一般可分为纵向收缩变形、横向收缩变形、弯曲变形、角变形、波浪变形、错边变形、扭曲变形7类. 在焊接过程中,焊接变形使得部分残余应力得到释放,但变形会影响到工件的尺寸精度,所以研究焊接变形很有必要[23-24].

    图11是不同覆材厚度下的复合板焊接残余变形模拟结果.可以看出,焊接残余变形最大的区域位于焊缝及其热影响区处, 其中覆材厚度为0.3 mm的复合板变形区间为0 ~ 4.67 mm,覆材厚度为1.0 mm的复合板变形区间为0 ~ 6.4 mm,覆材厚度为2.0 mm的复合板变形区间为0 ~ 7.02 mm,随着覆层厚度的增大,焊接残余变形也随之增大.

    图  11  不同覆材厚度复合板总变形对比
    Figure  11.  Comparison of total deformation of composite plates with different thickness of cladding material. (a) 0.3 mm;(b) 1.0 mm; (c) 2.0 mm

    图12xyz方向的焊后变形云图,可见z方向变形相对较小,但随着覆材厚度的增大,z方向的变形量增大,这里横向收缩变形指x方向变形量和y方向变形量相互叠加后的变形量,覆材厚度为0.3,1.0,2.0 mm复合板焊接后在x方向的变形分别为−0.6 ~ 0.6,−0.8 ~ 0.8,−0.9 ~ 0.9 mm,由于约束形式是自由约束,所以在x方向变形基本呈现对称式分布,可见覆材厚度对x方向的焊接残余变形影响较小;在y方向的变形分别为−4.66 ~ 0.0,−6.37 ~ 0.12,−7.0 ~ 0.06 mm,可见覆材厚度对y方向的焊接残余变形影响较大;在z方向的变形分别为−0.59 ~ 0.04,−0.66 ~ 0.04,−0.62 ~ 0.10 mm,覆材厚度对z方向的焊接残余变形影响较小.

    图  12  不同覆材厚度复合板 xyz 方向变形云图
    Figure  12.  Deformation cloud image in xyz direction of composite plate with different thickness of cladding materical. (a) x; (b) y; (c) z

    图13图14分别为不同覆材厚度的Q355及304表面垂直焊缝方向横向残余应力分布曲线. 从图13可以看到,Q355表面垂直焊缝方向横向残余应力峰值均低于或近似于屈服强度,焊缝区域呈现拉应力,焊缝或焊趾处应力出现峰值. 随着离焊缝中心线距离越远,应力降低,在板边缘处应力值趋于零. 其中覆材厚度为 0.3 mm的复合板应力峰值为301 MPa,覆材厚度为1.0 mm的复合板应力峰值为321 MPa,覆材厚度为2.0 mm的复合板应力峰值为345 MPa,随着覆材厚度的增加,Q355表面的应力也随之增大,其应力峰值接近材料屈服强度,这主要是因为材料自身性能的限制,以及多道焊的过程中,低合金高强钢表面为先焊区域,前道焊缝对后道焊缝有预热的作用,从而消除了材料的加工硬化的影响,使应力峰值在一个合理区间. 复合板两侧应力值趋于零,主要是由于横向和纵向收缩的叠加,而板边缘外侧应力得到释放, 从而应力值趋于零.

    图  13  Q355 表面垂直焊缝方向横向残余应力分布
    Figure  13.  Transverse residual stress distribution in the vertical weld direction of Q355 side
    图  14  304 表面垂直焊缝方向横向残余应力分布
    Figure  14.  Transverse residual stress distribution in the vertical weld direction of 304 side

    图14可以看出,横向残余应力为驼峰状,驼峰宽度与覆层表面焊缝宽度接近, 应力最大位于焊趾处,主要是由于母材和焊接填充材料的冷却过程不一致,从而形成强约束. 随着离焊缝中心距离越远,应力降低,在垂直于焊缝的模型边缘处应力趋近于零, 在覆层厚度变化过程中,0.3 mm和1.0 mm的热影响区宽度变化不明显,2.0 mm的热影响区宽度相对较宽,这是因为同温度下304不锈钢的热导率远小于Q355的热导率(常温下Q355的热导率约为304的3倍),覆材厚度2.0 mm模型的焊接道次增多,热输入有明显增大. 此外,应力最大均出现在驼峰处,其中覆材厚度为0.3 mm的复合板应力峰值为316 MPa,覆材厚度为1.0 mm的复合板应力峰值为338 MPa,覆材厚度为2.0 mm 的复合板应力峰值为356 MPa,随着覆材厚度的增大,304表面横向残余应力也随之增大.

    采用无损检测X射线衍射法,对覆材厚度2.0 mm的304/Q355复合板焊后表面残余应力进行测量,与模拟值比较分析. 图15为测量焊接残余应力的取点位置,图16为实测值与模拟值的对比图. 通过对横向残余应力和纵向残余应力的实测值与模拟值比较可知,实测值呈现中间区域为拉应力,两边为压应力的特征. 对于横向残余应力,在热影响区的测量值分布在模拟曲线附近,说明实测值与模拟结果有较好的一致性. 在焊缝区,测量值大部分点大于模拟值,其中最大相差95 MPa,主要是由于在工业实际中焊缝区的材料屈服强度稍高于母材,而在计算时对焊接填充的材料属性做了简化处理导致的. 对于纵向残余应力,其最大误差为11.6%,最小误差为1.12%,其中最大差值也出现在焊缝区,但就整体趋势而言,实测值和模拟值拟合较好,说明模拟值具有一定的参考价值.

    图  15  测量焊接残余应力的取点位置 (mm)
    Figure  15.  Location of the point for measuring welding residual stress
    图  16  横向和纵向残余应力实测值与模拟值的对比
    Figure  16.  Comparison of measured and simulated residual stress in horizontal and rertical. (a) horizontal; (b) rertical

    采用数值模拟与试验手段相结合的方法,通过对比覆材厚度1.0 mm的304/Q355复合板焊后焊缝截面形貌和多焊道热源模型叠加,及盖面截面形貌和单道热源模型,模拟熔池形貌与实际焊缝形貌相近,验证了经过校核的双椭球热源模型的准确性.

    随着覆材厚度的增加,304/Q355复合板焊接接头截面区最大残余应力逐渐减小,残余应力区宽度逐渐增大,复合材料截面位置的应力逐渐减小,在基层焊缝向过渡层焊缝焊接过程中,应力集中区逐渐向过渡层焊缝区收缩,最后最高应力集中区均位于过渡层焊缝区,最大应力接近材料的屈服强度.

    随着覆材厚度的增大,焊接残余变形逐渐增大,焊接残余变形最大的区域位于焊缝及其热影响区,覆材厚度0.3 mm的复合板变形量为0 ~ 4.67 mm,覆材厚度1.0 mm的复合板变形量为0 ~ 6.4 mm,覆材厚度为2.0 mm的复合板变形量为0 ~ 7.02 mm,随着覆材厚度的增加,Q355及304表面的应力随之增大,其应力峰值接近材料屈服强度.

  • 图  1   不同覆材厚度的坡口形式及焊道布置 (mm)

    Figure  1.   Groove form and pass arrangement. (a) 0.3 mm;(b) 1.0 mm;(c) 2.0 mm

    图  2   网格划分模型及局部放大

    Figure  2.   Meshing model and local magnificatio. (a) A; (b) B; (c) C

    图  3   约束条件

    Figure  3.   Constraints condition

    图  4   304和Q355力学—热物理性能参数

    Figure  4.   Mechanical — thermophysical performance parameter of Q355 and 304. (a) mechanical 304; (b) thermal physics 304; (c) mechanical Q355; (d) thermal physics Q355

    图  5   双椭球热源模型

    Figure  5.   Double ellipsoidal heat source model

    图  6   模拟熔池与实际熔池形貌对比

    Figure  6.   Comparison of simulated melt pool with outual melt pool morphology (a) superimposed comparison of macroscopic morphology of weld section and multi-pass welding; (b) comparison between simulation results of cover welding and macroscopic morphology of weld section

    图  7   不同覆材厚度复合板残余应力

    Figure  7.   Residual stress distribution of composite plate with different thickneos of cladding materical. (a) 0.3 mm;(b) 1.0 mm;(c) 2.0 mm

    图  8   覆材厚度0.3 mm的304/Q355不锈钢复合板焊缝横截面应力云图(MPa)

    Figure  8.   Stress cloud diagram of weld section of 304/Q355 stainless steel composite plate with 0.3 mm thickness ofcladding material

    图  9   覆材厚度1.0 mm的304/Q355不锈钢复合板焊缝横截面应力云图(MPa)

    Figure  9.   Stress cloud diagram of weld section of 304/Q355 stainless steel composite plate with 1.0 mm thickness ofcladding material

    图  10   覆材厚度2.0 mm的304/Q355不锈钢复合板焊缝横截面应力云图(MPa)

    Figure  10.   Stress cloud diagram of weld section of 304/Q355 stainless steel composite plate with 2.0 mm thickness ofcladding material

    图  11   不同覆材厚度复合板总变形对比

    Figure  11.   Comparison of total deformation of composite plates with different thickness of cladding material. (a) 0.3 mm;(b) 1.0 mm; (c) 2.0 mm

    图  12   不同覆材厚度复合板 xyz 方向变形云图

    Figure  12.   Deformation cloud image in xyz direction of composite plate with different thickness of cladding materical. (a) x; (b) y; (c) z

    图  13   Q355 表面垂直焊缝方向横向残余应力分布

    Figure  13.   Transverse residual stress distribution in the vertical weld direction of Q355 side

    图  14   304 表面垂直焊缝方向横向残余应力分布

    Figure  14.   Transverse residual stress distribution in the vertical weld direction of 304 side

    图  15   测量焊接残余应力的取点位置 (mm)

    Figure  15.   Location of the point for measuring welding residual stress

    图  16   横向和纵向残余应力实测值与模拟值的对比

    Figure  16.   Comparison of measured and simulated residual stress in horizontal and rertical. (a) horizontal; (b) rertical

    表  1   母材和焊材的化学成分(质量分数,%)

    Table  1   Chemical composition of base material and welding materials

    材料CSiMnPSNiCrMoCuN
    3040.081.002.000.0450.038.0018.000.1
    Q3550.20.551.600.0250.0250.300.300.40
    ER3080.081.000.500.0400.037.0022.000.750.75
    ER3090.151.000.50.0400.0312.0018.000.750.75
    ER50-60.1050.9750.1630.0150.0130.0160.06
    下载: 导出CSV

    表  2   焊接工艺参数

    Table  2   welding parameters

    焊缝道次焊接电流
    I/A
    电弧电压
    U/V
    焊接速度
    v/(mm·s−1)
    适用模型
    ① ~ ③100 ~ 1408 ~ 102A/B/C
    ④ ~ ⑤15012 ~ 152A/B/C
    ⑥ ~ ⑦1358 ~ 112C
    下载: 导出CSV
  • [1] 张婷, 许浩, 李仲杰. 层状金属复合材料的发展历程及现状[J]. 工程科学报, 2021, 43(1): 67 − 75.

    Zhang Ting, Xu Hao, Li Zhongjie. Development and current situation of layered metal composites[J]. Journal of Engineering Science, 2021, 43(1): 67 − 75.

    [2] 刘晓涛, 张延安, 崔建忠. 层状金属复合材料生产工艺及其新进展[J]. 材料导报, 2002(7): 41 − 43.

    Liu Xiaotao, Zhang Yanan, Cui Jianzhong. Production technology and new progress of layered metal composites[J]. Materials Reviews, 2002(7): 41 − 43.

    [3]

    Zhu Min, Wu Wei, Qian Weifang, et al. A brief review on welding of stainless steel clad plates: issues and future perspectives[J]. The International Journal of Advanced Manufacturing Technology, 2021, 110(1−2): 49 − 59.

    [4] 班慧勇, 梅镱潇, 石永久. 不锈钢复合钢材钢结构研究进展[J]. 工程力学, 2021, 38(6): 1 − 23.

    Ban Huiyong, Mei Yixiao, Shi Yongjiu. Research progress of stainless steel composite steel structure[J]. Engineering Mechanics, 2021, 38(6): 1 − 23.

    [5]

    Gan Shiming, Liu Huaying, Zhai Zhiping, et al. A review of welding residual stress test methods[J]. China Welding, 2022, 31(2): 45 − 55.

    [6] 叶丰, 彭星煜, 张丽江, 等. Incoloy 825/L360QS双金属复合管焊接残余应力的数值模拟[J]. 焊接, 2023(9): 48 − 53.

    Ye Feng, Peng Xingyu, Zhang Lijiang, et al. Numerical simulation of welding residual stress in Incoloy 825/L360QS bimetallic composite pipe[J]. Welding& Joining, 2023(9): 48 − 53.

    [7] 王莅. 304/Q345R 复合板焊接接头残余应力数值分析[D]. 太原: 太原科技大学, 2020.

    Wang Li. Numerical analysis of residual stress of welded joint of 304/Q345R composite plate[D]. Taiyuan: Taiyuan University of Science and Technology, 2020.

    [8]

    Han Tao, Gu Shiwei, Xu Liang, et al. Study on stress and deformation of keyhole gas tungsten arc-welded joints [J]. China Welding, 2020, 29(1): 17-25.

    [9] 袁嘉欣, 邵飞, 白林越, 等. 基于试验和数值模拟的TC1/1060/6061爆炸焊接复合板界面分析[J]. 焊接学报, 2023, 44(9): 81 − 87.

    Yuan Jiaxin, Shao Fei, Bai Linyue, et al. Interface analysis of TC1/1060/6061 explosive welded composite plate based on test and numerical simulation[J]. Transactions of the China Welding Institution, 2023, 44(9): 81 − 87.

    [10] 王铭. C276/A36 复合板激光焊接工艺及接头组织性能研究[D]. 镇江: 江苏科技大学, 2022.

    Wang Ming. Study on laser welding process and joint microstructure properties of C276/A36 composite plate [D]. Zhenjiang: Jiangsu University of Science and Technology, 2022.

    [11]

    Suo Li, Jian Li, Guangai Sun, et al. Modeling of welding residual stress in a dissimilar metal butt-welded joint between P92 ferritic steel and SUS304 austenitic stainless steel[J]. Journal of Materials Research and Technology, 2023, 23: 4938 − 4954. doi: 10.1016/j.jmrt.2023.02.123

    [12]

    Luo Wenze, Hu Long, Deng Dean. Numerical simulation and development of efficient calculation method for residual stress of SUS316 saddle tube-pipe joint[J]. Acta Metallurgica Sinica, 2022, 58(10): 1334 − 1348.

    [13]

    Deng Dean, Kyoshima Shoichi. Influence of annealing temperature on calculation accuracy of welding residual stress in a SUS304 stainless steel joint[J]. Acta Metallurgica Sinica, 2014, 50(5): 626 − 632.

    [14] 周壮壮, 陈放. 基于刻槽结构的复合杆爆炸焊接数值模拟[J]. 焊接学报, 2023, 44(8): 41 − 48. doi: 10.12073/j.hjxb.20220930001

    Zhou Zhuangzhuang, Chen Fang. Numerical simulation of explosive welding of composite rod based on notched structure[J]. Transactions of the China Welding Institution, 2023, 44(8): 41 − 48. doi: 10.12073/j.hjxb.20220930001

    [15]

    Wen Lianglu, Jing Longsun, Han Su, et al. Research on residual stress field and damage distribution for thick plate welded joints[J]. Hunan Daxue Xuebao/Journal of Hunan University Natural Sciences, 2014, 41(11): 24 − 31.

    [16]

    Mehran Ghafouri, Antti Ahola, Joseph Ahn, et al. Numerical and experimental inrestigations on the melding residual streoses and distortions of the short fillet melds in hish strength steel plates [J]. Engineering structwes, 2022, 260: 114269.

    [17]

    An Gyubaek, Park Jeongung, Lim Woontaek, et al. Characteristics of welding residual stress distribution in dissimilar weld joints[J]. Metals, 2022, 12(3): 405. doi: 10.3390/met12030405

    [18] 朱敏, 郑乔, 吴巍, 等. 坡口形式对双金属复合板多层多道焊接头残余应力演变的影响[J]. 机械工程学报, 2019, 58(10): 51 − 58.

    Zhu Min, Zheng Qiao, Wu Wei, et al. Influence of groove form on the evolution of residual stress in multi-pass welded joint of bimetal composite plate[J]. Journal of Mechanical Engineering, 2019, 58(10): 51 − 58.

    [19] 周广涛, 胡庆睿, 刘彪, 等. 随焊高速气流场辅助高强铝合金薄板焊接应力演变及控制变形机理[J]. 焊接学报, 2023, 44(2): 32 − 39.

    Zhou Guangtao, Hu Qingrui, Liu Biao, et al. Welding stress evolution and deformation control mechanism of high strength aluminum alloy sheet assisted by welding high speed gas field[J]. Transactions of the China Welding Institution, 2023, 44(2): 32 − 39.

    [20] 贾登峰. 304/Q345R过渡层焊缝对温度场和应力场影响的数值模拟研究[D]. 太原: 太原科技大学, 2018.

    Jia Dengfeng. Numerical simulation of effect of 304/Q345R transition layer weld on temperature field and stress field [D]. Taiyuan: Taiyuan University of Science and Technology, 2018.

    [21]

    Wei Yinhua, Dong Zhibin, Liu Bin, et al. Stress distributions of welding joints in titanium-steel composite pressure vessel under working conditions[J]. Sicence and Technology of Welding and Joining, 2011, 16(8): 709 − 716.

    [22]

    Hu Xiaodong, Yang Yicheng, Song Ming. Experimental and numerical investigations on the thermomechanical behavior of 304 stainless steel/Q345R composite plate weld joint[J]. Materials, 2019, 12: 3489.

    [23] 彭刚勇. 激光选区熔化成形钛合金温度场和应力场数值模拟[D]. 武汉: 华中科技大学, 2018.

    Peng Gangyong. Numerical simulation of temperature field and stress field of laser selective melting forming titanium alloy [D]. Wuhan: Huazhong University of Science and Technology, 2018.

    [24]

    Jiang Wenchang, Wang Binying, Gong Jianming, et al. Finite element analysis of the effect of welding heat input and layer number on residual stress in repair welds for a stainless steel clad plate[J]. Materials & Design, 2011, 32(5): 2851 − 2857.

  • 期刊类型引用(2)

    1. 李小兵,束长荣,陆立婷,陈程,郑传波,镇凡,麻晗. 不锈钢复合板的焊接工艺研究现状及进展. 中国冶金. 2025(01): 15-31 . 百度学术
    2. 顾传宇,严英杰,刘亚东,邓军,朱强,江秀臣. 有载分接开关真空管操动机构力学特性动态仿真. 上海理工大学学报. 2024(05): 475-483 . 百度学术

    其他类型引用(1)

图(16)  /  表(2)
计量
  • 文章访问数:  153
  • HTML全文浏览量:  12
  • PDF下载量:  40
  • 被引次数: 3
出版历程
  • 收稿日期:  2023-06-05
  • 网络出版日期:  2023-11-12
  • 刊出日期:  2024-01-30

目录

/

返回文章
返回