Simulation of solidification microstructure evolution in laser addition manufacturing of multicomponent alloy
-
摘要: 针对多元合金激光熔覆过程,基于有限元方法和元胞自动机技术建立了多元合金激光熔覆熔池三维传热传质及凝固组织形貌演变模型,通过自行开发的宏微观耦合接口程序实现了三维熔池模型和多元合金凝固组织演变模型的耦合. 针对IN718激光熔覆过程中的传热传质和凝固组织演变过程进行了模拟,研究了基板初始晶粒尺寸、异质形核及多层熔覆扫描路径对熔覆层凝固组织形貌的影响机理,并对模拟结果进行了试验验证. 结果表明,模拟结果与实际物理过程吻合较好,所开发的耦合模型能够真实反映多元合金激光熔覆过程.Abstract: For the process of laser cladding of multicomponent alloy, the model for heat and mass transfer of three-dimensional pool and solidification structure evolution were established based on the finite element method and cellular automata technology. The coupling of model for three-dimensional molten pool and model for multi-element alloy solidification structure evolution was realized by self-developed macro-micro coupling interface program. The heat and mass transfer and solidification structure evolution of IN718 laser cladding process were simulated. The effects of initial grain size, heterogeneous nucleation and scanning path of multilayer cladding on the solidification structure and morphology of cladding layer were investigated, and the simulation results were verified by experiments. The results showed that the simulation results were in good agreement with the actual physical process, and the developed coupling model truly reflected the laser cladding process of multicomponent alloys.
-
0. 序言
韧性作为评估钢材性能好坏的关键性能之一,直接影响钢铁材料设备的安全可靠使用. 大多数具有优异韧性的压力设备用钢材在经历焊接热循环后会产生一个韧性差的局部脆性区,焊件的热影响区常常是设备损坏的薄弱区[1-2] . 而解理裂纹最容易在这个区域起裂,并发生失稳扩展,进而导致断裂[3]. 研究发现[4-7],热影响区(heat affected zone,HAZ)中的粗晶热影响区(coarse grained heat affected zone,CGHAZ)会生成粗大的脆性组织. 在这个区域原始奥氏体晶粒会显著粗化,由于显微组织的粗化和基体上M-A组元析出,使得CGHAZ冲击韧性相对于母材严重下降[8-9]. 通常认为CGHAZ是HAZ中典型的最薄弱区域. 但是最薄弱区域对冲击韧性影响的详细机制尚不清楚. 文章通过实际施焊确定了HAZ的最薄弱区域, 分析了09MnNiDR埋弧焊焊接接头HAZ的组织、硬度、断口形貌以及最薄弱区域与冲击韧性之间的关系.
1. 试验方法
试验所用母材为09MnNiDR压力容器低温钢,板厚60 mm,采用埋弧焊,加工U形坡口,焊前预热温度为150 ℃,焊后经610 ℃ × 3 h退火处理. 具体的焊接工艺参数:热输入E = 20 kJ/cm,焊接电流I = 510 ~ 520 A,焊接电压U = 30 V,焊接速度v = 28 m/h. 焊接接头热影响区宏观形貌如图1所示. 采用的焊材为H09MnNiDR焊丝和SJ208DR焊剂. 母材化学成分如表1所示.
表 1 09MnNiDR钢化学成分(质量分数,%)Table 1. Chemical composition of 09MnNiDR steelC Si Mn Al Ni Nb 0.08 0.35 1.47 0.02 0.55 0.02 在板厚1/4处取冲击试样,缺口方向垂直于焊接方向(图2). 冲击试样的夏比V形缺口分别取在母材(base metal, BM),焊缝(weld metal, WM)以及热影响区各个不同的区域,CGHAZ缺口开在熔合线处,细晶热影响区(fine grained heat affected zone, FGHAZ)缺口开在距离熔合线1 mm处,临界热影响区(intercritical heat affected zone,ICHAZ)缺口开在距离熔合线1.5 mm处,亚临界热影响区(subcritical heat affected zone,SCHAZ) 缺口开在距离熔合线2 mm处. 通过带有千分尺的光学显微镜在左右移动距离为0.1 mm的范围内统计出缺口尖端粗晶热影响区宽度所占8 mm缺口尖端前沿的比例(图2c). 采用微机控制摆锤冲击试验机进行−70 ℃的低温冲击试验,Quanta450FEG型扫描电子显微镜(SEM)分析组织、断口形貌. 通过HAT-1000A数字显示显微硬度仪测量显微硬度.
2. 试验结果
2.1 焊接接头各区域的冲击吸收能量
焊接接头−70 ℃平均冲击吸收能量如表2所示. BM,SCHAZ,ICHAZ,FGHAZ的冲击吸收能量都在270 J以上,而FGHAZ表现出最高的冲击吸收能量,WM的冲击韧性表现良好,冲击吸收能量平均达到139 J. 焊接接头各区域的断口形貌如图3所示. BM,SCHAZ,ICHAZ,FGHAZ都是以大量韧窝和撕裂脊组成的韧性断裂. WM的断口上呈现解理面与韧性撕裂脊共存的现象,表现为准解理断裂. 通过将缺口开在熔合线处,得到含有不同比例CGHAZ的焊接接头试样,大量试验发现这些试样冲击吸收能量在20 ~ 320 J的范围内,在这个区域冲击吸收能量出现巨大的波动,断口形态也有显著差异,将在2.3节中单独讨论.
表 2 焊接接头−70 ℃下各区域的冲击吸收能量(J)Table 2. Impact energy of each region for 09MnNiDR welded joints at −70 ℃BM SCHAZ ICHAZ FGHAZ CGHAZ WM 276 281 274 304 20 ~ 320 139 2.2 焊接接头各区域的微观组织及硬度
焊接接头不同区域的组织如图4所示. 09MnNiDR母材组织为铁素体,且基体上弥散分布有少许碳化物,SCHAZ, ICHAZ, FGHAZ组织均为块状铁素体. 不同的是,ICHAZ组织为尺寸不均匀的块状铁素体,这可能是一部分没有发生奥氏体化的铁素体在冷却的过程中继续长大,从而ICHAZ的组织由晶粒尺寸差异较大的铁素体组成. FGHAZ为晶粒尺寸均匀且细小的块状铁素体组成. WM的一次组织为粗大的柱状晶,分布有少量碳化物,在晶界内形成细小的针状铁素体,二次组织为块状铁素体,且尺寸较为均匀. CGHAZ组织发生严重的粗化,在冷却的过程中,产生粗大的复合组织 (粒状贝氏体 + 板条贝氏体 ).
通过比较晶粒尺寸差异较大的BM,FGHAZ和CGHAZ(图5),发现BM的平均晶粒尺寸为15.7 µm,最大为38 µm;FGHAZ的平均晶粒尺寸为11.8 µm,最大为24 µm;而CGHAZ的平均晶粒尺寸为29.5 µm,最大为63 µm. CGHAZ晶粒尺寸粗化显著. 焊接接头各区域的显微维氏硬度如图6所示,SCHAZ硬度相比于母材发生少许软化,平均在155 ~ 166 HV,其它区域并未产生软化现象,而临近熔合线的CGHAZ的复合组织表现出最高的硬度值,平均为207 ~ 226 HV.
2.3 粗晶热影响区宽度占缺口尖端前沿比例与冲击吸收能量的关系
图7为缺口尖端开在熔合线处不同冲击吸收能量下的冲击试样的宏观和微观断口形貌,缺口尖端前沿的组织比例具有随机性,使得这些试样的冲击吸收能量波动巨大. 当冲击吸收能量为20 J时,测得粗晶区宽度所占比例达到整个缺口尖端前沿的96.5%,断口呈现完全的解理断裂. 随着缺口尖端前沿CGHAZ宽度逐渐减小,冲击吸收能量提高. 由于CGHAZ区域为完全脆性解理断裂特征,从图3可知焊接接头热影响区域中的其它区域的断口特征均为韧窝断裂形态,所以直接可以通过断口上脆性解理断裂的区域来测量缺口尖端前沿粗晶热影响区宽度分布的比例. 根据测量结果统计了图8所示缺口尖端前沿粗晶热影响区的分布比例与冲击吸收能量之间的关系,发现当缺口前沿无粗晶热影响区时,冲击吸收能量达到260 J以上,试样断裂形式为完全的韧性断裂,断口完全被韧性断裂区所覆盖. 粗晶热影响区占8 mm缺口尖端前沿的比例达到25%左右时,冲击吸收能量降低到202 J. 当比例达到70% ~ 80%时,冲击吸收能量分别下降到60,48,30 J左右,相比母材冲击韧性损失达到了80% ~ 90%. 总体趋势为:随着缺口尖端前沿粗晶热影响区占8 mm区域的比例的增加,冲击吸收能量逐渐降低.
图 7 缺口前沿占有不同比例粗晶热影响区的冲击试样断口的宏观和微观形貌Figure 7. Macroscopic and microscopic fracture surfaces of impact specimens with different CGHAZ factions. (a) macroscopic fracture-20 J; (b) microscopic fracture-20 J; (c) macroscopic fracture-60 J; (d) macroscopic fracture-60 J; (e) microscopic fracture-141 J; (f) macroscopic fracture-141 J; (g) macroscopic fracture-202 J; (h) microscopic fracture-202 J2.4 临界事件
Chen等人[10]提出解理断裂的3个准则,分别为裂纹形核控制的应变准则,控制裂纹钝化的三向应力度准则以及控制裂纹扩展的应力准则. 而临界事件被定义为这3个准则当中最困难的阶段,这个阶段控制解理断裂微观开裂过程. 对于起裂但终止扩展而残留在断裂试样上的微裂纹进行观察,以确定解理断裂的裂纹扩展过程中最艰难的阶段-临界事件. 并且残留裂纹通常在具有最大阻力的晶粒边界或者第二相被中断,只有完全解理的断口侧面的残留裂纹才能更准确的反应临界事件. 故可以观察残留裂纹的扩展路径去推断试样断裂过程中最难发生的阶段. 图9所示为CGHAZ比例占96.5%、最低冲击吸收能量为20 J的冲击试样断口侧面的残留裂纹.从图9可知残留裂纹有三种形态,一是穿过晶界,沿着晶界扩展;二是穿过晶粒边界,裂纹发生偏转,最终停留在晶粒内;三是直接停留在一个晶界上.以上形态也说明对于铁素体钢缺口试样,粗晶热影响区组织的晶粒边界对裂纹的扩展具有关键性的阻碍作用. 为了进一步推断临界事件,对完全解理的冲击试样断口上的解理小面半长轴进行测量如图7b所示,发现解理小面半长轴ac大约在26 ~ 35 µm左右,与最大晶粒尺寸的1/2 (31.5 µm)能较好的吻合,可以推论在铁素体钢缺口试样中解理断裂的临界事件是晶粒尺寸的裂纹穿过晶粒/晶粒边界的扩展,粗大晶粒为关键的微观组元,最粗大晶粒的尺寸表征了临界事件的尺寸.
3. 分析与讨论
3.1 焊接接头最薄弱区域的确定
试验中CGHAZ和FGHAZ最大亚晶粒尺寸分别为63和24 μm. CGHAZ亚晶粒尺寸显著增大,FGHAZ的铁素体晶粒分布均匀且细小,平均尺寸为11.8 μm. FGHAZ的平均晶粒尺寸仅为CGHAZ的1/3,FGHAZ微观组织具有更好的韧性. 因此,韧性的提高是由细小的晶粒引起的,而对于CGHAZ,由于硬度较高的粒状贝氏体和板条贝氏体组织的出现,使得硬度值平均达到了207 ~ 226 HV,大于焊接接头任何区域的硬度. 由于晶粒尺寸的粗化以及出现的复合组织表现出最差的韧性[11-13]. 因此,细小的块状铁素体表现出最优的冲击韧性,而CGHAZ中粗大的粒状贝氏体和铁素体组成的复合组织表现出最差的冲击韧性. 从显微组织确定CGHAZ是焊接接头的最薄弱区域. Chen等人[10]提出临界事件决定了解理微观开裂的过程中最薄弱的区域,对于解理微观开裂过程,裂纹起裂和扩展产生于最薄弱的组元(最薄弱的区域,即粗大的晶粒)中,而不是产生于最小晶粒的最强组织组元中,最大的晶粒提供最低的阻力,而最小的晶粒则提供最大的阻力,这意味着尽管形成临界事件的临界阶段为整体开裂过程提供了更大的困难,但是微观开裂仍然发生在临界开裂阶段中最低阻力的最薄弱区域(微观组元)上. 文中缺口试样中解理断裂的临界事件是晶粒尺寸的裂纹穿过晶粒/晶粒边界的扩展,即最困难的阶段. 而粗大的复合组织只在CGHAZ发现,其它亚区域并不存在. CGHAZ中由复合组织组成的粗大晶粒为关键的微观组元(薄弱组元),最粗大晶粒的尺寸称为临界事件的尺寸. 因而确定了临界事件,也就认定了裂纹扩展的最困难阶段以及最弱的阻力来自于哪个微观组元,从而也就确定了试样的最薄弱区域. 图7a中所示为CGHAZ比例占96.5%、最低冲击吸收能量为20 J的冲击试样其断裂为完全的解理断裂形态,这也进一步说明09MnNiDR中CGHAZ是最薄弱的区域.
3.2 粗晶热影响区比例对冲击韧性的影响
取样位置的微小变化会造成冲击韧性的显著差异,随着缺口尖端前沿的CGHAZ宽度所占比例增加,冲击吸收能量逐渐降低,直至达到20 J. 表明缺口尖端的微观组织所占比例对冲击韧性有很大的影响,冲击吸收能量从260 J左右变化到20 J,断裂模式从韧性断裂变为解理断裂,其中最主要原因是在缺口尖端前沿8 mm区域上存在一定比例的CGHAZ. 这表明除了由粗大的粒状贝氏体,板条贝氏体和块状铁素体组织组成的CGHAZ外,其它区域的组织在试验温度下表现出更好的韧性. 而起裂源处的解理小面半长轴尺寸和CGHAZ最大晶粒尺寸的1/2相一致,表明CGHAZ中粗大晶粒是起裂位置. 因此随着CGHAZ比例的增加,冲击韧性逐渐降低. 在实际焊接条件下,虽然无法获得完整CGHAZ的冲击韧性. 但是,根据图8中CGHAZ分布比例与冲击韧性的关系,可以通过多项式拟合获得CGHAZ完整的冲击韧性.
综上所述,焊接接头的冲击韧性的差异主要来自缺口尖端前沿CGHAZ的比例. 而晶粒尺寸以及微观组织决定了CGHAZ的强韧性. 因此CGHAZ在缺口尖端前沿的分布也就影响着试样整个焊接接头的冲击韧性. 文中揭示了09MnNiDR焊接接头最薄弱区域以及最薄弱区域与冲击韧性之间的关系,对实际焊接热影响区冲击韧性的评估建立了一定的理论根据.
4. 结论
(1) 揭示了09MnNiDR低温钢焊接接头在实际施焊热输入为20 kJ/cm的情况下,最薄弱区域为CGHAZ,CGHAZ粗大的复合组织 (粒状贝氏体 + 板条贝氏体 ) 使得冲击吸收能量最低在20 J左右,相对于母材冲击吸收能量降低约92.7%,粗大晶粒的出现以及复合组织的转变是韧性恶化的重要因素.
(2) 随着缺口尖端前沿粗晶热影响区的比例越大,冲击韧性越低. CGHAZ比例达到70% ~ 80% 时冲击吸收能量分别降低到60,48,30 J,相比母材韧性损失达到80% ~ 90%,充分体现出薄弱环节对接头冲击韧性的损害有巨大的影响.
-
表 1 熔覆层和热源参数
Table 1 Deposited layer and heat source parameters
表 2 IN718热物性参数
Table 2 IN718 thermophysical parameters
密度ρ/(kg·m−3) 液相温度Tl/K 固相温度Ts/K 导热系数k/(W·m−1·K−1) 熔化潜热Lf/(kJ·kg−1) 热熔cp/(J·kg·K−1) 7 734 1 609 1 533 0.56 + 2.9 × 10−2T – 7 × 10−6T 290 435 参考文献[9] 参考文献[9] 参考文献[9] 参考文献[17] 参考文献[9] 参考文献[9] 粘度μ/(Pa·s) α2/(m·s−1·K−2) α3/(m·s−1·K−3) 热膨胀系数β/K−1 表面张力系数 ${\frac{ { {\rm{d} }\gamma } }{ { {\rm{d} }T} }}$ /(N·m−1·K−1)5.31 × 10−3 1.58 × 10−5 2.29 × 10−6 1.3 × 10−5 2.5 × 10−5 参考文献[9] 参考文献[11] 参考文献[11] 参考文献[9] 参考文献[9] -
[1] 黄安国, 余圣甫, 李志远. 焊缝金属凝固组织元胞自动机模拟[J]. 焊接学报, 2008, 29(4): 45 − 48. doi: 10.3321/j.issn:0253-360X.2008.04.012 Huang Anguo, Yu Shengfu, Li Zhiyuan. Simulation on weld metal solidifying microstructure with cellular automaton[J]. Transactions of the China Welding Institution, 2008, 29(4): 45 − 48. doi: 10.3321/j.issn:0253-360X.2008.04.012
[2] 马瑞, 董志波, 魏艳红, 等. 镍基合金焊缝凝固组织演变过程模拟和仿真[J]. 焊接学报, 2010, 31(7): 43 − 46. Ma Rui, Dong Zhibo, Wei Yanhong, et al. Simulation of solidification microstructure evolution in molten pool of nickel base alloy[J]. Transactions of the China Welding Institution, 2010, 31(7): 43 − 46.
[3] Gandin C, Rappaz M. A 3D cellular automaton algorithm for the prediction of dendritic grain growth[J]. Acta Materialia, 1997, 45(5): 2187 − 2195. doi: 10.1016/S1359-6454(96)00303-5
[4] Yin H, Felicelli S. Dendrite growth simulation during solidification in the LENS process[J]. Acta Materialia, 2010, 58(4): 1455 − 1465. doi: 10.1016/j.actamat.2009.10.053
[5] 魏雷, 林鑫, 王猛, 等. 激光立体成形中熔池凝固微观组织的元胞自动机模拟[J]. 物理学报, 2015, 64(01): 356 − 363. Wei Lei, Lin Xin, Wang Meng, et al. Cellular automaton simulation of the molten pool of laser solid forming process[J]. Acta Physica Sinica, 2015, 64(01): 356 − 363.
[6] 王志坚, 王宗园, 宋鸿武, 等. Ti-6Al-4V激光快速成形熔池凝固过程研究[J]. 机械设计与制造, 2017, 01(08): 103 − 105. doi: 10.3969/j.issn.1001-3997.2017.08.029 Wang Zhijian, Wang Zongyuan, Song Hongwu, et al. Study on solidification process of the molten pool of laser rapid forming of Ti-6Al-4V[J]. Machinery Design & Manufacture, 2017, 01(08): 103 − 105. doi: 10.3969/j.issn.1001-3997.2017.08.029
[7] 张敏, 周玉兰, 薛覃, 等. Ti-45Al合金焊接熔池凝固过程数值模拟[J]. 焊接学报, 2018, 39(3): 6 − 10. Zhang Min, Zhou Yulan, Xue Tan, et al. Numerical simulation of solidification process of Ti-45Al alloy weld pool[J]. Transactions of the China Welding Institution, 2018, 39(3): 6 − 10.
[8] Knapp G, Mukherjee T, Zuback J, et al. Building blocks for a digital twin of additive manufacturing[J]. Acta Materialia, 2017, 135: 390 − 399. doi: 10.1016/j.actamat.2017.06.039
[9] Pottlacher G, Hosaeus H, Kaschnitz E, et al. Thermophysical properties of solid and liquid Inconel 718 Alloy[J]. Scandinavian Journal of Metallurgy, 2002, 31: 161–168.
[10] Parimi L, Ravi G, Clark D, et al. Microstructural and texture development in direct laser fabricated IN718[J]. Materials Characterization, 2014, 89: 102 − 111. doi: 10.1016/j.matchar.2013.12.012
[11] Panwisawas C, Qiu C, Anderson M, et al. Mesoscale modelling of selective laser melting: thermal fluid dynamics and microstructural evolution[J]. Computational Materials Science, 2017, 126: 479 − 490. doi: 10.1016/j.commatsci.2016.10.011
[12] Kurz W, Giovanola B, Trivedi R. Theory of microstructural development during rapid solidification[J]. Acta Metallurgica, 1986, 34(5): 823 − 830. doi: 10.1016/0001-6160(86)90056-8
[13] Rappaz M, Gandin C. Probabilistic modelling of microstructure formation in solidification processes[J]. Acta Metallurgica et Materialia, 1993, 41(2): 345 − 360. doi: 10.1016/0956-7151(93)90065-Z
[14] Gandin C, Rappaz M. A coupled finite element-cellular automaton model for the prediction of dendritic grain structures in solidification processes[J]. Acta Metallurgica et Materialia, 1994, 42(7): 2233 − 2246. doi: 10.1016/0956-7151(94)90302-6
[15] Gandin C, Desbiolles G, Rappaz M, et al. A three-dimensional cellular automaton-finite element model for the prediction of solidification grain structures[J]. Metallurgical and Materials Transactions A, 1999, 30(12): 3153 − 3165. doi: 10.1007/s11661-999-0226-2
[16] Zinoviev A, Zinovieva O, Ploshikhin V, et al. Evolution of grain structure during laser additive manufacutring. simulation by a cellular automata method[J]. Materials & Design, 2016, 106: 321 − 329.
[17] Wei H, Knapp G, Mukherjee T, et al. Three-dimensional grain growth during multi-layer printing of a nickel-based alloy Inconel 718[J]. Additive Manufacturing, 2019, 25: 448 − 459. doi: 10.1016/j.addma.2018.11.028
-
期刊类型引用(6)
1. 刘金浩,李家辰,张亮亮,吴宝生,李鹏,董红刚. Al-Mg合金搅拌摩擦焊接头的组织演变及腐蚀性能. 机械制造文摘(焊接分册). 2025(02): 13-22+36 . 百度学术
2. 金聪聪,黄立兵,黄文彬,王大力,马月婷,李鹏,董红刚. 新型铝合金MIG焊接头微观组织与力学性能. 焊接学报. 2024(07): 74-82 . 本站查看
3. 薛新莲,臧华平,王洪建,李顺方. 铝合金板激光切割实验工艺研究. 激光杂志. 2024(09): 183-187 . 百度学术
4. 陈晓,闫成旗,代俊辉. 基于IRSM的7075铝合金焊接接头抗拉强度可靠性优化设计. 焊接. 2024(10): 32-41 . 百度学术
5. 刘金浩,李家辰,张亮亮,吴宝生,李鹏,董红刚. Al-Mg合金搅拌摩擦焊接头的组织演变及腐蚀性能. 焊接学报. 2024(10): 8-18 . 本站查看
6. 邓伟,罗坤,雷基林,宋仲模,徐远志,张勇. 基于灰色关联性分析的铝合金缸体铸造工艺参数优化. 特种铸造及有色合金. 2024(11): 1511-1517 . 百度学术
其他类型引用(6)