高级检索

API X65管道深水铺设GMAW横向焊接温度场

周灿丰, 陈智, 焦向东, 罗雨, 高辉, 周珍珍

周灿丰, 陈智, 焦向东, 罗雨, 高辉, 周珍珍. API X65管道深水铺设GMAW横向焊接温度场[J]. 焊接学报, 2020, 41(9): 60-68. DOI: 10.12073/j.hjxb.20200310001
引用本文: 周灿丰, 陈智, 焦向东, 罗雨, 高辉, 周珍珍. API X65管道深水铺设GMAW横向焊接温度场[J]. 焊接学报, 2020, 41(9): 60-68. DOI: 10.12073/j.hjxb.20200310001
ZHOU Canfeng, CHEN Zhi, JIAO Xiangdong, LUO Yu, GAO Hui, ZHOU Zhenzhen. Study on temperature field of GMAW horizontal welding for deep water laying of API X65 pipe[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2020, 41(9): 60-68. DOI: 10.12073/j.hjxb.20200310001
Citation: ZHOU Canfeng, CHEN Zhi, JIAO Xiangdong, LUO Yu, GAO Hui, ZHOU Zhenzhen. Study on temperature field of GMAW horizontal welding for deep water laying of API X65 pipe[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2020, 41(9): 60-68. DOI: 10.12073/j.hjxb.20200310001

API X65管道深水铺设GMAW横向焊接温度场

基金项目: 国家自然科学基金项目(51305037);北京市自然科学基金项目(KZ201210017017).
详细信息
    作者简介:

    周灿丰,1970年出生,教授. 主要从事火下焊接与焊接自动化技术研究. Email:Canfeng@bipt.edu.cn.

  • 中图分类号: TG 402

Study on temperature field of GMAW horizontal welding for deep water laying of API X65 pipe

  • 摘要: 为了应对管道深水铺设GMAW横向焊接熔池下坠、窄坡口、无衬垫等挑战,建立了管道三维模型并进行了网格划分,利用材料软件生成了API X65温度场材料数据库,通过模拟温度场云图与实际焊缝截面的匹配校核了Goldak双椭球热源模型,管道焊接温度场分布SYSWELD软件模拟结果与实际焊接过程相符合. 建立了焊接温度场测试系统,通过背孔法埋设的热电偶测量了API X65管线钢平板焊接的热循环曲线. 测试曲线与SYSWELD仿真曲线具有相同的温度升降趋势,上板节点温度也均明显低于下板节点温度. 采用焊缝截面平均热循环曲线作为热源进行了多道焊模拟,并进行了试验. 结果表明,多道焊焊接变形是由焊接工艺参数与残余应力释放共同确定的.
    Abstract: To meet challenges of GMAW horizontal welding in deep water pipeline laying, such as molten pool falling, narrow gap groove and without backing support, a 3D model of pipe is established and meshed, temperature field material database of API X65 is generated by material software,the Goldak double ellipsoid heat source model is checked by matching the simulated temperature field cloud chart with the actual weld cross section, and temperature field distributions in pipe welding simulated by SYSWELD software are consistent with the actual welding process. The welding temperature field test system is established, and the thermal cycle curve of API X65 pipeline steel plate welding is measured by the thermocouple embedded in the back hole method. Tested curve has the same temperature rise and fall trend as the SYSWELD simulated curve, and temperature of upper plate node is significantly lower than that of lower plate node. The average thermal cycle curve of weld cross section is used as the heat source to simulate the multi pass welding, and the test is carried out. The results show that the welding deformation of multi pass welding is determined by the welding process parameters and residual stress release.
  • 选区激光熔化(selective laser melting, SLM)是激光增材技术的一种,以金属或陶瓷粉末为原料,在激光的熔化作用下快速成形出具有高致密度的零件[1-2].316L是常用的铁基材料因为出色的力学性能与耐腐蚀性等.然而,该工艺下的零件表面常极为粗糙,影响应用[3].鉴于SLM成形过程的迅速性和微小性,不少研究者开始建立激光热源与金属粉末的相互作用数值计算模型—介观尺度模型,研究SLM中工艺参数对熔道形貌特征的影响来合理调节扫描参数,以对成形件表面形貌进行优化[4].

    在介观尺度层面,基于离散元法的激光增材工艺介观数值模拟中,最早由德国Carolin Körner等人[5]于2011年基于二维玻尔兹曼模型实现,可以对粉末熔化再凝固的过程进行预测.然而,这仅是二维模型,无法对成形件的表面的三维形貌进行模拟及实际试验验证.后续的研究者主要焦距于三维化的基于离散元与流体体积法的介观数值模型的合理性.吴宇车等人[6-7]以TC4和H13为例建立SLM工艺下基于离散元法和流体体积法的三维介观数值模型,对在不同功率和速度下的熔道形貌进行研究,证实了数值模拟中加入蒸发反冲力模型的必要性;向羽等人[8]在基于离散元和流体体积法的介观数值模拟中,加入蒸汽反冲力模型的情况下进行TC4的SLM单熔道成形计算,有效得出了成形理想熔道形貌的工艺参数组合;单雨薇[9]对激光粉床增材的介观模拟中的常被采用的多种热源模型进行比较,发现表面热源模型模拟下的熔池形状相比体积热源模型更加具有合理性.叶唯娟[10]从上述被逐渐完善的介观数值模型出发,对两种不同路径的单层多道扫描进行研究,发现有效熔宽的扫描间距和Z形的扫描方式下能得到更为稳定的温度场. Hen等人[11]也在此类介观模型上对SLM下的单层多熔道下的成形件熔道形貌特征进行了研究,但未涉及到不同扫描间距对熔道形貌的影响性.

    可见,当前基于离散元和流体体积法的SLM介观尺度数值模拟,大都局限于与单熔道成形效果有关的扫描工艺参数研究,即使有单层多道的数值相关研究,对不同扫描间距下成形的熔道形貌的关注也较少.SLM的本质是层层累积叠加成形,每一层皆由复数条单熔道扫描搭接而成,不同的扫描间距对于成形件的表面形貌必定存在较大影响,客观上单熔道的数值模拟研究存在一定缺陷.文中也将基于离散元法和流体体积法,从单层双道的角度出发,对316L粉末的成形过程进行瞬态过程的数值模拟及验证,从表面的熔道形貌和熔道宽度两个方面,研究SLM工艺下激光功率、扫描速度和扫描间距3个扫描参数对成形件表面形貌特征的影响.

    数值模型的建立要符合基本的工艺过程,SLM的基本工艺过程可分为粉末铺设和激光扫描2个阶段.文中将直接建立粉末铺设完成后的粉床模型和介观层面下的瞬态的熔道成形模型.为能够合理阐明SLM工艺下的诸如热传导、蒸发、凝固、表面张力等效应,在模型建立过程中需要一定简化假设,具体如下:①熔化流动的金属为不可压缩牛顿流体;②所有不锈钢粉末在单个上均为形貌均匀的球体;③不考虑蒸发带来的质量损失且气液界面始终稳定[12-13].

    使用离散元模拟工具EDEM,建立单层堆积密度0.4,计算域为1 020 μm × 408 μm × 61.2 μm的长方体的球状单层粉末局部堆积模型,该计算域的选择满足熔道设计长宽、粉末层的最高粉末颗粒高度以及后处理过程中对于成形过程物理现象的观测.球状粉末D50为30 μm且在10 ~ 60 μm范围内呈正态分布.使用只考虑颗粒间的法向与切向作用力Hertz-Mindlin(no slip)模型作为颗粒生成过程中的接触模型[14].

    SLM工艺的热源输入是激光束,材料因激光所带的自由光子与材料的原子碰撞被加热.在数值模拟中[15-16],激光通常被定义为移动的高斯热源,自由光子的热量较低,易被反射难以穿透到材料表面以下较深的位置,因此在基于激光束的加热模式中,采用平面热源模型更具有合理性.

    $$ Q\left( r \right) = \frac{{2AP}}{{{\text{π}} {R^2}}}\exp \left( {\frac{{ - 2{r^2}}}{{{R^2}}}} \right) $$ (1)

    式中:A为粉末的热量的吸收率;P为激光功率;R表示半径为50 μm的光斑其数据源于设备说明;r为熔化区域任一点到激光中心距离.

    在SLM成形过程中,热量交互行为主要有热传导,热对流、热辐射和蒸发散热.基本传热方程为

    $$ \rho c\frac{{\partial T}}{{\partial Z}} = Q\left( r \right) - {q_{\rm{c}}} - k\nabla T - {q_{{\rm{e}}{\rm{v}}}} $$ (2)

    式中:$ \rho $为材料密度;$ c $为材料比热;$ T $为初始温度;$ k $为材料热导率;${q}_{{\rm{c}}}$为气体的对流换热;$ k\nabla T $为表面热辐射.

    在熔池的自由界面与上方保护气体间存在的热量交换边界条件为

    $$ {q_{{\rm{in}}}} = q - {h_{\rm{c}}}\left( {T - {T_0}} \right) - {\varepsilon _{\rm{r}}}{\sigma _{\rm{s}}}\left( {{T^4} - {T_0}^4} \right) - {q_{{\rm{ev}}}} $$ (3)

    式中:$ {q_{{\rm{in}}}} $为熔池自由界面上的热流变化;$ {h_{\rm{c}}} $为自然对流传热系数;$ {\varepsilon _{\rm{r}}} $为发射率;$ {\sigma _{\rm{s}}} $为斯特凡−玻尔兹曼常数;$ T $为熔池上表面温度;$ {T_0} $为环境温度.

    金属粉末熔化发生蒸发散热控制方程为

    $$ {q_{{\rm{ev}}}} = 0.82\frac{{\Delta H^* {P_0}}}{{\sqrt {2{\text{π}} MRT} }}{e^{\left( {\Delta * \tfrac{{T - {T_{{\rm{lv}}}}}}{{RT{T_{{\rm{lv}}}}}}} \right)}} $$ (4)

    式中:$ M $为材料摩尔质量;$ R $为理想气体常数;$ {P}_{0} $为环境压力;$ {T}_{\mathrm{l}\mathrm{v}} $为材料沸点;$ \Delta {H}^{\mathrm{*}} $为损失金属蒸气有效焓.

    为简化研究模型和提高数值计算可靠性,文中未考虑热力耦合问题.由于热量分布不均涉及到的热变形涉及高度非线性的热弹塑性问题,同时残余应力的试验测试精度也难以得到准确验证.

    粉末材料在热源范围内被加热至熔点以上形成熔池,熔池上部和下部的密度以及温度不同使得内部存在自然对流,其流动遵循N-S方程,在对流上主要考虑重力,表面张力以及蒸汽反冲力对熔池稳定性的影响[17].

    表面张力在熔池上主要体现为Marangoni效应,即

    $$ \left\{ \begin{gathered} {f_{\rm{s}}} = \sigma \kappa n\left| {\nabla \alpha } \right|\frac{{2\rho }}{{{\rho _{316{\rm{L}}}} + {\rho _{\text{α}} }}} \\ {f_{\rm{m}}} = \frac{{d\sigma }}{{dT}}\left[ {\nabla T - \left( {n \cdot \nabla T} \right)n} \right]\left| {\nabla \left( {1 - \alpha } \right)} \right|\frac{{2\rho }}{{{\rho _{316{\rm{L}}}} + {\rho _{\text{α}} }}} \\ \end{gathered} \right. $$ (5)

    式中:${f}_{{\rm{s}}}$为熔池表面张力;$ \sigma $为表面张力系数;n为熔池表面单位法向量;$ \kappa $表示曲率;$ \left|\nabla \alpha \right| $ 将界面力转化为体积表面力;${{f}}_{{\rm{m}}}$为熔池中Marangoni力.

    当熔池最高温度超过金属沸点时,金属在蒸发溢出同时对熔池表面产生一个反冲压力,亦称蒸汽反冲力,表达式为

    $$ {p_{\rm{r}}} = 0.54{P_0}\exp \left( {{L_{{\rm{lv}}}}\frac{{T - {T_{{\rm{lv}}}}}}{{RT{T_{{\rm{lv}}}}}}} \right) $$ (6)

    在FLOW-3D软件中,对传热与流动相关的控制方程进行求解;使用流体体积法对熔池流动及凝固行为实现跟踪.

    在SLM中,材料的部分热物参数难以通过试验手段获取,文中将材料的热物属性看作仅是与温度有关的变量,采用JMATPRO软件进行计算得到材料的密度、比热、粘度和热导率的热物理参数,如图1所示.

    图  1  316L热物参数
    Figure  1.  Theromphysicl parameters of 316L. (a) liquid viscosity; (b) specific heat; (c) thermal conductivity; (d) density

    在基于有限元的数值模拟中,需要考虑同种材料的粉末情况与实体块情况下的导热差异等物性参数.文中的数值模拟研究模型是基于离散元法的建模,不再需考虑粉体与块体材料的物理差异,表1为材料的参数设置[18].

    表  1  316L材料参数
    Table  1.  316L material parameters
    固相密度ρ/(kg·m−3) 固相温度T/K 液相温度T/K 沸点温度T/K 熔化潜热Hf /(J·kg−1) 蒸发潜热Hv /(J·kg−1) 表面张力r /(N·m−1) 吸收率A
    7 850 1 658 1 723 3 090 2.7 × 105 7.45 × 106 1.6 0.25
    下载: 导出CSV 
    | 显示表格

    试验涉及参数自变量多,为提高数值计算效率,设置三水平三因素正交试验,见表2,不同水平按照经验范围进行选定.

    表  2  正交试验设计
    Table  2.  Experimental design using orthogonal table
    参数变量A功率P/WB速度v/(m·s−1)C间距S/μm
    水平12000.650
    水平22500.970
    水平33001.290
    下载: 导出CSV 
    | 显示表格

    选区激光熔化的本质是若干条预定扫描轨迹下形成多条熔道的搭接集合,为简化数值模型的复杂度和运算时间,采用同扫描层相邻的双熔道数值模拟的特征效果来分析实际成形件上表面可能出现的熔道形貌特征.为方便后续分析,对正交试验不同组合进行编号,见表3.

    表  3  试验参数及编组
    Table  3.  Experimental parameters and grouping
    组号功率 P/W速度 v /(m·s−1)间距 S /μm编号
    12000.650A1B1C1
    22000.970A1B2C2
    32001.290A1B3C3
    42500.670A2B1C2
    52500.990A2B2C3
    62501.250A2B3C1
    73000.690A3B1C3
    83000.950A3B2C1
    93001.270A3B3C2
    下载: 导出CSV 
    | 显示表格

    图2为数值模拟的正交试验设计下的9组单层双道形貌图.图2a2d2g2h这4组参数组合下的熔道成形效果更为接近理想状态的熔道形貌,如图3所示.这4组参数下的双道形貌,搭接部分平滑,未熔合缺陷特征少且熔道平直,是良好的的表面形貌.

    图  2  数值模拟熔道形貌
    Figure  2.  Numerical simulation of melt tracks morphology. (a) group No. 1 No. A1B1C1; (b) group No. 2 No. A1B2C2; (c) group No. 3 No. A1B3C3; (d) group No. 4 No. A2B1C2; (e) group No. 5 No. A2B2C3; (f) group No. 6 No. A2B3C1; (g) group No. 7 No. A3B1C3; (h) group No. 8 No. A3B2C1; (i) group No. 9 No. A3B3C2
    图  3  理想熔道形貌示意
    Figure  3.  Ideal melt track morphology

    图2e虽然熔道整体平直,但有明显岛屿球化特征;图2f图2i参数组合,激光功率更大,虽然粉床受热源照射部分熔化效果好,但扫描速度较快,也存在明显未熔合缺陷和岛屿球化效应;图2b图2c下的参数组合成形的熔道形貌存在明显缩颈、熔道小角度扭曲和未熔合特征,相邻熔道的搭接不明显,激光功率低,扫描速度较快,粉末熔化时较小的熔池尺寸导致湿润角增大,在表面张力作用下极易发生岛屿球化形成不连续的熔道,不生成规律性纹路,最终形成缺陷特征多的熔道形貌.

    研究中所涉及的扫描参数包括扫描速度、激光功率和扫描间距.涉及参数较多且参数之间对于熔道形貌特征的影响存在一定耦合.

    为更好评估各扫描参数与熔道良好形貌特征之间的关系,引入综合参数激光线性能量密度(line energy density,LED),即激光功率与扫描速度的比值,评估激光的能量输入[19-20].

    图4所示,成形效果良好状态表面的熔道形貌特征需要较高的线能量密度和扫描间距的参数组合,这4组参数组合从左往右依次是1,8,4,7组,对应的线能量密度依次为333,333,416,500 J/m,对应扫描间距为0.05,0.05,0.07,0.09 mm.其中1组和8组在图中数据重合.在非良好熔道形貌区域中,线能量密度较低,在线能量密度较低的情况下,由于Plateau-rayleigh不稳定性,熔体趋向于从圆柱状向球状变化以减少其表面能,导致出现球化、不连续的熔道等不良的表面形貌特征的风险增加.虽然也有研究表示,当线能量密度较高时,熔池面积变大且表面会受到较大的反冲作用,凝固后的熔道表面粗糙[21].文中的数值结果中并无相关体现,这可能是受材料、堆积密度、参数设置、统计方法等因素影响下造成的结果误差.

    图  4  线能量密度(LED)与扫描间距对熔道形貌特征影响
    Figure  4.  Influence of linear energy density (LED) and scanning spacing on the morphology characteristics of the melt tracks

    图4的良好熔道形貌区域可知,扫描间距与线能量密度在所设置的数值模拟参数范围内满足某种线性关系,其表达式为

    $$ {x_1} = 0.239\;5x - 29.719 $$ (7)

    式中:$ {x_1} $为扫描间距,单位 μm;$ x $为线能量密度,单位J/m.当线能量密度和扫描间距满足该式时,可以得到较好的熔道形貌.

    对以上4组参数组合进行各参数的变化分析:7组和8组,功率一致,扫描速度上升50%,扫描间距下降44%;对比1组和4组,在速度一致情况下,第4组功率提升25%,间距提升40%;对比1组和8组,在间距一致情况下,第8组功率提升33%,速度提升50%.说明在影响熔道形貌完整度上,扫描速度、扫描间距和扫描功率是依次递减的. 同时,高的扫描速度下液相时间缩短,有利于抑制球化现象,提高致密度[22-23].

    图5所示,在所有9组计算结果中,熔道的宽度基本随着扫描间距和线能量密度的增加而增加,且相比扫描间距,线能量密度对熔宽的影响更大.在50 μm间距下,有2组数据点高度重合,说明在不同功率与速度的相同比值下的线能量密度对于熔宽影响微乎.

    图  5  不同扫描间距下的线能量密度与熔宽
    Figure  5.  Line energy density and melting width under different scanning distances

    具有良好熔道形貌的4组参数组合均在3组数据曲线的最末端,可显著看出线能量密度和扫描间距对熔宽的影响也具备线性的特征.

    为保证对数值模拟结果的可靠性,由表3参数组号进行实际试验,通过打印样件,观察成形件表面熔道形貌与非直接的测量熔道宽度的方式与数值结果进行对比,以验证数值模拟结果的有效性.

    使用易博三维公司研发的IGAM-2金属打印机和商业化316L粉末分别对9组不同的扫描参数组合进行单个试样打印,样件为长宽高为7 mm × 7 mm × 5 mm的长方体块.

    局限于设备因素,无法单独在每个试件表面进行指定熔道数量的扫描及作断面切割观测.双道的本质是多道的局部体现且不涉及成形后合金组织形貌研究,因此在结果观察时将直接观察成形件最表层的形貌特征.使用超景深三维显微系统VHX-1000E进行对最表层熔道形貌图像和熔道宽度数据进行采集,图6为超景深系统下的316L粉末图像.

    图  6  316L粉末形状轮廓图像
    Figure  6.  316L powder shape outline image

    在基于粉床的金属增材工艺中,粉末须满足粒径细小、粒度分布窄、球形度高等要求.从图6可以看出,粉末基本为球形.经激光粒度测试,该批316L不锈钢金属粉末的粒径分布为D10 = 24.7 μm,D50 = 35.5 μm,D90 = 50.7 μm,其中平面平均径D[3,2]与体积平均径D[4,3]分别为34.1,36.7 μm,说明该批粉末形状较为规则,粒径较为集中.图7为粉末的粒径分布图像,基本符合正态分布,这有利于粉层颗粒分布的均一性和粉末熔化与铺设过程中良好的流动性[24-25].

    图  7  316L粉末粒径分布
    Figure  7.  316L powder particle size distributio

    图8为试验后不同成形件在超景深显微镜下的表面形貌图像,其中各个分图的数字和字母编号对应表3中的正交设计表.

    图  8  不同参数组合下的成形件表面形貌
    Figure  8.  Forming surface topography under the different parameters combination conditions. (a) group No. 1 No. A1B1C1; (b) group No. 2 No. A1B2C2; (c) group No. 3 No. A1B3C3; (d) group No. 4 No. A2B1C2; (e) group No. 5 No. A2B2C3; (f) group No. 6 No. A2B3C1; (g) group No. 7 No. A3B1C3; (h) group No. 8 No. A3B2C1; (i) group No. 9 No. A3B3C2

    当激光功率为200 W时,在不同扫描速度和扫描间距下,成形件表面的熔道形貌如图8a ~ 8c所示.图8a与数值结果的熔道与实际结果具有较高一致性,但也存在少量球化特征,试件表面熔道并不完全平整,这可能与基板平台与刮刀的定位精度和平整度有关.图8b中熔道存在明显缩颈及熔道小角度扭曲状态,且同层熔道之间的搭接线不再明显,所对应的数值模拟结果中并未出现熔道不连续状态,可能是受限于模拟的熔道长度过短而没有展现.在图8c中,熔道特征难以辨认,大量岛屿球化和不连续缺陷.图8a ~ 8c的表面熔道形貌特征结果与数值模拟基本一致.

    当激光功率为250 W时,对应不同扫描速度和扫描间距的成形件表面形貌如图8d ~ 8f所示.相比图8a8d的表面更为平坦,在更高的激光功率下,粉末的熔化效果更好.图8e中由于扫描间距过大,相邻熔道搭接线不明显,较快的扫描速度使得熔化效果不足,熔道不连续且有弯曲与岛屿球化现象.图8f相比图8e有更为平直连续的熔道,这是更高的扫描速度下较小的扫描间距可以有效弥补前一条熔道的因扫描过程中能量密度不足而造成的未熔合特征.

    图8e图8f所对应的数值结果与试验验证结果略有出入,在数值模拟结果中,图8e中的熔道形貌相比于图8f更为平直连续且不存在明显缩颈;但在实际打印结果中,图8f中的熔道形貌略比图8e中平直连续,这说明数值模型在有在细节上有待完善的地方:一是模拟的熔道长度仅有几百微米长的区域,在范围尺度上可能缺少一定广泛性;二是数值模型的成形底面是假设平整的基板平面,实际成形时,仅有第一层是在较为平整的基板面上进行,其余层均是在上一层具有相对于基板平面存在高低起伏的面上进行粉末层铺设、扫描、熔化凝固过程.

    当激光功率为300 W时,对应不同扫描速度和扫描间距的成形件表面形貌如图8g ~ 8i所示. 图8g图8h中,在低扫描速度下(0.6 m/s和0.9 m/s),熔道平直完整,球化特征少.即使在0.6 m/s的低扫描速度下搭配较大的扫描间距,由于较高的能量密度输入,熔池面积大、粘度低,有利于熔池凝固前的铺展成形,熔道成形效果优于低扫描间距(50 μm),高扫描速度(0.9 m/s)参数组合.图8i中,成形件表面明显熔道不连续,扭曲现象严重且伴随大量岛屿球化特征,与图8c的特征基本一致,但图8i下的熔道比图8c中更为连续,这说明在高的扫描速度下,大幅提升激光功率和适当减小扫描间距对成形件表面形貌完整度的提升不太明显.

    在实际成形的熔宽方面,在共聚焦系统下采集到的较为完整的单道和双道的熔宽数据和数值结果所测量结果对比如图9所示.

    图  9  熔道宽度及误差
    Figure  9.  Width and error of the melt track

    不同线能量密度和扫描间距下,无论是数值模拟的计算结果还是试验结果,熔道宽度变化趋势基本一致,平均相对误差分别为单道15.91%,双道20.98%,由于多种客观原因比如刮刀的不平整、基板平面校准的误差、光路系统定位与聚焦平面的准确性等,试验与仿真之间存在误差是不可避免的,其中最大的误差体现为在9组参数组合中存在2组在表面的熔道形貌成形效果不佳的参数组合的缺陷特征的体现上.其中,数值结果下的表面熔道特征与实际结果误差较小的1,4,7,8组双道成形及搭接效果良好,熔道平整连续,搭接区域过渡平整,与数值模拟部分的结论吻合.

    综上,文中的数值模拟可以为SLM工艺下不同激光功率、扫描间距和扫描速度参数组合下的成形件的表面的熔道形貌特征提供有效参考.

    (1)对于成形效果良好表面的熔道形貌,在线能量密度313 ~ 500 J/m与扫描间距50 ~ 90 μm之间存在明显的线性关系,虽然这种线性关系可能是参数梯度设置不够小的结果.

    (2)扫描速度、扫描间距和激光功率对成形理想表面形貌的影响由大到小依次为扫描速度、扫描间距、扫描功率.

    (3)熔道宽度增大与减小主要受线能量密度影响.

  • 图  1   海底管道J形铺设示意与焊接工作站

    Figure  1.   Sketch map of J-lay and welding workstation. (a) J-laying of subsea pipeling; (b) weling station with double torches

    图  2   管道CAD图 (mm)

    Figure  2.   Pipe CAD diagram. (a) sizes of the narrow groove; (b) 3D of pipe

    图  3   有限元模型网格分组

    Figure  3.   Set collectors of finite elements. (a) 3D unit of metal and weld; (b) 2D unit of heat exchange surface; (c) constraint location node; (d) welding line and reference line; (e) welding start node; (f) welding start unit

    图  4   X65材料CCT曲线

    Figure  4.   CCT curve of API X65

    图  5   X65材料热物性曲线

    Figure  5.   Thermal parameter curves of API X65. (a) thermal conductivity curve; (b) specific heat capacity curve; (c) density distribution curve

    图  6   双椭球热源模型

    Figure  6.   Double ellipsoid heat source model

    图  7   模拟温度场云图与实际焊缝截面比对

    Figure  7.   Comparison between simulated temperature field cloud chart and actual weld cross section

    图  8   管道焊接温度场分布模拟结果

    Figure  8.   Simulated results of temperature field distribution in pipe welding. (a) at 1 s; (b) at 2 s ; (c) at 4 s; (d) at 8 s; (e) at 10 s; (f) at 14 s; (g) at 50 s; (h) at 2 000 s

    图  9   焊接设备

    Figure  9.   Welding equipment

    图  10   温度采集系统

    Figure  10.   Temperature acquisition system

    图  11   测温点分布图

    Figure  11.   Distribution of temperature measurement points. (a) front view; (b) top view

    图  12   板材试验节点热循环曲线

    Figure  12.   Thermal cycle curves of plate test nodes

    图  13   管道模拟节点热循环曲线

    Figure  13.   Thermal cycle curves of pipe simulated nodes

    图  14   热循环曲线

    Figure  14.   Thermal cycle curves. (a) nodes thermal cycle curves of weld cross section; (b) average thermal cycle curve of weld cross section

    图  15   多道焊焊接变形模拟结果(mm)

    Figure  15.   Simulated results of welding deformation(mm)

    图  16   横向多道焊焊接变形试验结果(mm)

    Figure  16.   Tested results of welding deformation(mm)

    图  17   焊接变形曲线

    Figure  17.   Welding deformation curves

    表  1   API X65钢的化学成分(质量分数,%)

    Table  1   Chemical composition of API X65

    CSiMnCrAlNPSFe
    0.18 0.23 1.4 0.2 0.026 0.028 0.021 0.006 余量
    下载: 导出CSV

    表  2   焊接工艺参数

    Table  2   Welding process parameters

    焊层焊接电流I/A焊接电压U/V焊接速度v/(mm·s−1)送丝速度vs/(m·min−1)焊缝熔深H/mm
    打底1703012.510.52.7
    2 ~ 3180301210.51.8
    4 ~ 61803012112
    7 ~ 8200301411.53
    下载: 导出CSV
  • [1]

    Gong Shunfeng, Zhang Tao, Wang Xipeng, et al. Numerical simulation on dynamic behaviour of deepwater J-lay systems[J]. Ocean Engineering, 2020, 196: 106771. doi: 10.1016/j.oceaneng.2019.106771

    [2]

    Kwon Y J, Nam B W, Kim N W, et al. A model test on a rigid pipe installation using J-Lay method[C]//The 28th International Ocean and Polar Engineering Conference. Sapporo, Japan, 2018, 255 − 261.

    [3]

    Mark J Kaiser. Offshore pipeline construction cost in the U. S. Gulf of Mexico[J]. Marine Policy, 2017, 82: 147 − 166. doi: 10.1016/j.marpol.2017.05.003

    [4]

    Zhong Wenjun, He Ning, Ni Wenchi. Structural analysis of J-lay tower on deepwater semi-submersible lifting and pipe laying vessel[C]//ISOPE International Journal of Offshore and Polar Engineering Conference. Rhode Island, Greecce, 2016, 676 − 683.

    [5]

    Lirola F, Pionetti F R, Agoumi J, et al. Development and qualification of an innovative and cost efficient heat traced flowline optimized for J-Laying[C]//Offshore Technology Conference. Houston, Texas, USA, 2016, 255 − 261.

    [6]

    Senthila B, Selvam Panner R. Dynamic analysis of a J-lay pipeline[J]. Procedia Engineering, 2015, 116: 730 − 737. doi: 10.1016/j.proeng.2015.08.358

    [7]

    Etienne Girault, François Lirola, Erwan Briand, et al. Lessons learned from the qualification program and J-laying campaign of 12-inch steel Catenary risers for the P55 deepwater project[C]//Offshore Technology Conference. Houston, Texas, USA, 2015, 1 − 24.

    [8]

    Heerema E P. Recent achievements and present trends in deepwater pipe-lay systems[C]//Offshore Technology Conference. Houston, Texas, USA, 2005, 1 − 4.

    [9]

    Killeen J P, Dyson K C , Reeves K D, et al. Large diameter SCR delivery challenges[C]//Offshore Technology Conference. Houston, Texas, USA, 2005, 1 − 13.

    [10]

    Van der Graaf J, Wolbers D, Boerkamp P. Field experience with the construction of large diameter SCR in deep water[C]//Offshore Technology Conference. Houston, Texas, USA, 2005.

    [11]

    Dick Wolbers, Rob Hovinga. Installation of deepwater pipelines with sled assemblies using the new J-Lay system of the DCV balder[C]//Offshore Technology Conference. Houston, Texas, USA, 2003, 1 − 12.

    [12]

    Mazar Atabaki M, Yazdian N, Ma J, et al. High power laser welding of thick steel plates in a horizontal butt joint configuration[J]. Optics & Laser Technology, 2016, 83: 1 − 12.

    [13]

    Guo Wei, Liu Qiang, Francis John A, et al. Comparison of laser welds in thick section S700 high-strength steel manufactured in flat (1G) and horizontal (2G) positions[J]. CIRP Annals - Manufacturing Technology, 2015, 64: 197 − 200. doi: 10.1016/j.cirp.2015.04.070

    [14]

    Vishvesh J B. Effect of residual magnetism on sidewall fusion in narrow gap gas metal arc welding[J]. Welding and Cutting, 2008, 7(5): 280 − 288.

    [15]

    Kang Y H, Na S J. Characteristics of welding and arc signal in narrow groove gas metal arc welding using electromagnetic arc oscillation[J]. Welding Research, 2003, 5: 93 − 99.

    [16]

    Goldak A, Bibby M, Moore J, et al. Computer modeling of heat flow in welding heat source[J]. Metallurgical and materials transactions, 1986, 17(3): 587 − 600.

    [17]

    Leblond J B, Devaux J. A new kinetic model for aniso-thermal metallurgical transformations in steels including effect of austenite grain size[J]. Acta metallurgica, 1984, 32(1): 137 − 146. doi: 10.1016/0001-6160(84)90211-6

    [18]

    Zhou Canfeng, Jiao Xiangdong, Chen Jiaqing, et al. Design of welding system applied in deepwater sub-sea pipeline laying[J]. Welding & Joining, 2010(7): 16-20.

    [19]

    Macdonald K A, Maddox S J. New guidance for fatigue design of pipeline girth welds[J]. Engineering Failure Analysis, 2003, 10(2): 177 − 197. doi: 10.1016/S1350-6307(02)00051-1

    [20]

    David Walters, Ricky Thethi. A step change application of threaded and coupled connections[J]. Offshore Pipeline Technology, 2002: 1 − 15.

  • 期刊类型引用(2)

    1. 樊炳倩,沈闲. 激光选区熔化参数对热物理过程影响数学模型分析. 激光杂志. 2025(01): 222-227 . 百度学术
    2. 袁健,何斌,张俊飞,陈国炎,张锁荣. 纳秒激光参数对316L不锈钢沟槽结构加工的影响. 金属加工(热加工). 2024(05): 61-66 . 百度学术

    其他类型引用(2)

图(17)  /  表(2)
计量
  • 文章访问数:  370
  • HTML全文浏览量:  22
  • PDF下载量:  14
  • 被引次数: 4
出版历程
  • 收稿日期:  2020-03-09
  • 网络出版日期:  2020-12-02
  • 刊出日期:  2020-09-24

目录

/

返回文章
返回