Inconel625合金激光熔覆过程中显微组织演变的数值模拟.pdf

97MATERIALSFOR MECHANICAL ENGINEERINGVol.47No.7Jul.20232023年7 月第47 卷第7 期2023机械工程材料D0I:10.11973/jxgccl202307016Inconel6255合金激光熔覆过程中显微组织演变的数值模拟李继红，郭钊，李保铃，张敏，郭宇飞，苟川东（西安理工大学材料科学与工程学院，西安7 10 0 48）摘要：基于MATLAB模拟平台并结合晶粒的形核与生长理论，采用元胞自动机-有限差分方法构建了Inconel625合金激光熔覆层显微组织演变的数值计算模型，利用该模型对单道激光熔覆过程中合金的显微组织演变进行了模拟，并进行了试验验证，同时研究了激光热输入和异质形核数量对熔覆层晶粒形貌的影响。结果表明：采用建立的模型模拟得到熔覆层内部由柱状晶以及分布在柱状晶之间的亚晶粒组成，靠近熔合线区域由胞状晶组成，与试验结果相吻合；模拟和试验得到柱状晶的平均宽度分别为2.8 17，2.7 43m，相对误差仅为2.6 9%，验证了计算模型的可靠性；随着激光热输入的减小或异质形核数量的增大，熔覆层组织逐渐细化，晶粒数量增多，柱状晶宽度减小。关键词：激光熔覆；数值模拟；元胞自动机-有限差分法；Inconel625合金；显微组织中图分类号：TG401文献标志码：A文章编号：10 0 0-37 38(2 0 2 3)0 7-0 0 9 7-0 7Numerical Simulation of Microstructure Evolution of Inconel625 AlloyDuring Laser CladdingLI Jihong,GUO Zhao,LI Baoling,ZHANG Min,GUO Yufei,GOU Chuandong(School of Materials Science and Engineering,Xian University of Technology,Xian 710048,China)Abstract:On the basis of the MATLAB simulation platform and the grain nucleation and growth theory,thenumerical calculation model of microstructure evolution of laser cladding layer of Inconel625 alloy was constructedby cellular antomata-finite difference method.The microstructure evolution of the alloy during single pass lasercladding was simulated by the model and verified by experiments.The effects of laser heat input and heterogeneousnucleation quantity on the grain morphology of the cladding layer were studied.The results show that according tothe established model simulation,the inner part of the cladding layer was composed of columnar crystals and sub-grains distributed between columnar crystals,and the cladding layer near the fusion line consisted of cellularcrystals;the simulation was consistent with the test results.The average size of columnar crystals obtained bysimulation and test were 2.817,2.743 m,respectively,and the relative error was only 2.69%,which verified thereliability of the calculation model.With decreasing laser heat input or increasing heterogeneous nucleation quantity,the microstructure of the cladding layer was gradually refined,the number of grains increased,and the size of thecolumnar crystals decreased.Key words:laser cladding;numerical simulation;cellular automata-finite difference method;Inconel625alloy;microstructure0引言激光熔覆作为一种表面改性工艺，能够有效地收稿日期：2 0 2 2-0 2-15；修订日期：2 0 2 3-0 2-0 3基金项目：国家自然科学基金资助项目（519 7 42 43）作者简介：李继红（19 7 3一），男，河南洛阳人，副教授，博士改善材料表面的耐腐蚀、耐磨损和抗氧化等性能，延长工件的使用寿命 1-3。与传统的电弧堆焊相比，激光熔覆具有工艺精度高、稀释率低、熔覆层组织细密等优势；相比于喷涂方法，激光熔覆层与基体之间结合强度高，因此激光熔覆的应用前景十分广阔 4-。Inconel625合金具有润湿性好、高温自润滑特性优良等优点，广泛应用于石油化工、海洋工程等领域关98MATERIALS FOR MECHANICAL ENGINEERING李继红，等：Inconel625合金激光熔覆过程中显微组织演变的数值模拟2023机械工程材料键设备的表面改性工艺中 6-7。激光熔覆是一个高温、瞬态、复杂的物理冶金及化学冶金过程，难以采用传统的试验方法进行实时、全过程的深入研究 8-9 。计算机技术和数值模拟技术的迅速发展则为该工艺及其机制的研究提供了有力的分析工具10TAN等 11采用元胞自动机法和相场法建立了一种多尺度模型，用于研究激光熔覆过程中枝晶生长以及熔覆参数对熔覆层形状和微观结构的影响，熔覆层形貌和枝晶间距的计算结果与试验结果吻合良好。SHIMONO等 12 将多相场法与相图计算软件相结合，模拟了钛合金激光熔覆过程中柱状晶向等轴晶的转变过程。ACHARYA等 13 采用相场模型对激光熔覆过程中的二次元素偏析、枝晶尺寸、枝晶取向等微观结构进行了模拟，发现枝晶取向对表面粗糙度和扫描速度有很强的依赖性。孙道金等1采用宏微观相场模型模拟了纯镍激光熔覆过程中温度场与显微组织的演变，发现熔池内枝晶生长方向略偏于激光扫描方向。刘朝阳等 15 建立了三维模型，模拟了激光多道熔覆镍基单晶高温合金显微组织演变，模拟结果和试验结果吻合良好。目前，国内外关于激光熔覆的数值模拟都已经取得了一定的研究成果，但是尚未对Inconel625镍基合金激光熔覆过程中局部区域显微组织演变及合金元素的偏析行为进行深入的数值模拟研究。鉴于此，作者结合前人的研究成果，在MATLAB模拟平台上建立了Inconel625合金激光熔覆过程熔覆层显微组织转变的元胞自动机-有限差分（cellularautomata-finite difference,CA-FD)模型，利用所建模型对单道激光熔覆过程中熔覆层显微组织的演变进行了模拟，并进行了试验验证，同时研究了激光热输入和异质形核参数对熔覆层晶粒形貌的影响，以期为激光熔覆工艺的制定和改进提供参考。1CA-FD模型的建立1.1温度场计算模型的构建作者基于二维尺度对单道激光熔覆过程进行模拟和分析，激光熔覆过程及熔覆层横截面如图1所示。为了便于模型的构建以及计算分析，同时考虑到激光熔覆较低的稀释率，忽略基体材料化学成分的影响，模型中将Inconel625合金激光熔覆熔池视为Ni-Nb二元合金（Ni-3.76Nb合金），该合金的热物性参数如表1所示，激光熔覆工艺参数为激光激激光束保护气体空气粉末流熔覆层基体熔覆层熔覆方向容覆层横截面基体图1激光熔覆过程及熔覆层横截面示意Fig.1Schematic of laser cladding process and cross-section ofcladding layer表1Ni-Nb合金的热物性参数Table1Thermo-physical property paramters of Ni-Nb alloy参数数值密度/(kgm=3)8440热导率/(Wm-1.K-1)12.1比热容/(Jkg-1.K-1)430锯溶质质量分数/%3.76液相线温度/1350液相线斜率-2.75溶质平衡分配系数0.41Gibbs-Thompson 系数/(Km)3.65X10-7液相扩散系数/(m.s-1)3X10-9固相扩散系数/（m.s-1）310-12功率5kW、扫描速度5mms-l、光斑直径5mm。光热源具有能量集中的特点，易在时间和空间上形成梯度很大的温度场。不均匀的温度场会直接影响熔覆层的显微组织演变及形貌。为了简化计算模型，将热源的热流密度近似以高斯函数加以描述，具体表达式 16 为qr=Qexp(-Kr2)(1)式中：q，为距离热源中心距离为r处的热流密度；Q为作用于试样表面的有效激光功率；K为能量集中系数。在激光熔覆过程中，激光热源产生的热量以辐射的形式传递到基体试样表面，使得合金粉末和基体试样表面局部区域快速熔化形成熔池。热量以对流和传导的形式在熔池和基体试样内部进行传递，并且基体试样表面还与周围环境通过对流和辐射进行热量交换。假设热量在垂直于熔覆方向的二维平面内传播，忽略熔池内的对流，则二维瞬态导热方程可表示为99MATERIALSFOR MECHANICAL ENGINEERING李继红，等nconel625合金激光熔覆过程中显微组织演变的数值模拟2023机械工程材料aTaTaT入+Q(2)at式中：为熔覆层材料的密度；c为熔覆层材料的比热容；入为熔覆层材料的热导率；T为模拟区域温度；，y 为二维平面坐标；t为时间。基于课题组前期的研究成果 17 ，利用有限差分法求解式（2），即可得到模拟区域宏观单元网格的温度场分布。设定模拟区域的初始温度为室温Tr。传热问题的边界条件通常有3种形式：边界温度已知、边界热流已知和边界自由热交换。作者研究对象为室温环境下激光熔覆层显微组织的演变过程，因此采用模拟区域（熔覆层)与周围空气自由换热这种更接近于物理现实的边界条件形式 18 ，具体表达式如下：qw=h(Tb-T,)(3)式中：T为熔覆层与空气接触边界处的温度；h为熔覆层与周围空气接触的换热系数；qw为熔覆层和周围空气的热交换量。1.2元胞自动机模型的建立激光熔覆层的形核模型采用RAPPAZ等 19 提出的由统计法建立的、基于高斯分布函数的准连续形核模型。该模型既考虑了过冷度对熔覆层非均匀形核的影响，又考虑了形核过程的连续性。在某一过冷度下，晶粒形核密度的表达式为dNN(AT-T2maxmaxd(T)V2元exp2AT。(4)式中：T为过冷度；N为晶粒形核密度；Nmax为非均匀形核密度的最大值；T。为曲率过冷度；T，max为最大形核过冷度。过冷度是晶粒生长的驱动力，形核后晶粒需在一定过冷度下才能持续长大。过冷度由热过冷度T w、成分过冷度T。和曲率过冷度T。组成。由Gibbs-Thompson关系 2 0 1可得t时刻的过冷度为T(t)=T L-T(t)+m C L(t)-C -I K(t)(5)式中：TL为熔覆层材料液相线温度；T（t)为t时刻模拟区域微观单元的瞬时温度；mL为熔覆层材料液相线的斜率；C。为熔覆层材料的溶质初始浓度；CL（t）为t时刻液相中的溶质浓度；为Gibbs-Thompson系数；K（t）为t时刻的固/液界面平均曲率。界面元胞是处于固相和液相之间的元胞，也称糊状元胞。在凝固过程中界面元胞会向周围的液相元胞中排出多余溶质，时间步长t内排出的溶质浓度C(t)的表达式 2 1为DLtC(t):ZCi(t)-Cm,(t)(6)(）n式中：DL为溶质的液相扩散系数；nb表示界面元胞周围的液相元胞个数；Cn.（t）为液相元胞的溶质浓度；c为网格尺寸。在固/液界面推进过程中，界面元胞排出溶质后，元胞单元的固相率会发生变化，固相率增量f。的计算公式 2 2 为C(t)Af.=C.(O(1-k。)1+A(1+2n rand)(7)式中：k。为溶质平衡分配系数；A为扰动系数；n rand为 0,1 区间内的随机数。在固/液界面推进过程中，排到周围相邻液相元胞中的多余溶质会导致这些元胞的溶质浓度升高，此时液相元胞之间会出现较大的浓度梯度，因此还需要对液相元胞进行溶质扩散计算，其控制方程 2 3为aCi(t)a?CL(t)a?CL(t)DL(8)at1.3宏微观温度场模型的转换利用最近邻插值法将宏观温度场模型转换为微观温度场模型，其原理如图2 所示，图中：O代表宏观单元网格，P代表微观单元网格。具体插值方法 2 4 如下：88Tp=2L7T:/2L7(9)5-5式中：Tp为微观单元P的温度；T；为P周围宏观单元的温度；L；为P到周围宏观单元的距离；i为从5到8 求和的底数。6N73WE25S814(a)宏观网格(b)微观网格图2插值原理示意Fig.2SSchematic of interpolation principle:(a)macro grid and(b)micro grid可见微观单元网格的温度受到宏观单元网格温100MATERIALS FOR MECHANICAL.ENGINEERING李继红，等：I n c o n e l6 2 5合金激光熔覆过程中显微组织演变的数值模拟机械工程材料度以及距离L的影响。2模拟结果及试验验证2.1熔覆层凝固过中程晶粒的生长由图3可以发现：在凝固初期（0.35s），由于熔覆层底部和边缘具有较大的过冷度和温度梯度，因此这些区域快速形核并长出了大量胞状晶；随着凝固过程的进行（0.7 5s），熔覆层边缘区域的晶粒沿着温度下降最快的方向向熔覆层中心区域不断生长，部分胞状晶转变为柱状晶，并且由于不同晶粒生长方向存在差异，生长方向与温度梯度较大的方向平行的晶粒生长速率更快，抑制着周围其他晶粒的生长，导致晶粒间出现竞争生长的现象；随着凝固过程的进一步进行（1.15s），由于生长空间的限制，晶粒间的竞争生长更加明显，只有少数生长方向具有优势的晶粒继续长大，其余晶粒停止生长；随着时间的再推移（1.55s），熔覆层晶粒生长过程结束，液态金属完全凝固，熔覆层主要由大量柱状晶以及少量胞状晶组成，1000um1000m1000 m1 000 um(a)0.35 s(b)0.75 s(c)1.15 s(d)1.55 s图3熔覆层凝固过程中不同时刻晶粒生长的模拟结果Fig.3Simulation of crystal growth during cladding layer solidification at different times为进一步研究熔覆层的显微组织演变过程，在更小尺度下对不同时刻熔合线附近区域晶粒的形核与生长过程进行模拟。由图4可以看出：在凝固初期（0.0 0 5s），在较小的成分过冷条件下，熔合线附近产生大量晶核并生长出了大量胞状晶；随着凝固过程的进行（0.0 2 5s），由于生长空间的限制，靠近熔合线附近的胞状晶在生长到一定程度后，因与其周围的晶粒接触而停止生长，靠近液态金属区域的晶粒则会沿着温度下降的方向继续向熔覆层内部生长，同时在较大温度梯度的作用下转变成为柱状晶；随着凝固过程的继续进行（0.0 7 5s），柱状晶沿着温度梯度方向不断生长，并且出现明显的竞争生长现象，此时部分柱状晶被周围相邻的晶粒“吞并”，导致凝固过程中出现一定程度的晶粒粗化现象；随着柱状晶的进一步长大，在0.12 5s时整个模拟区域完全凝固，近熔合线位置为由胞状晶组成的细晶区域，离熔合线一定距离的内部则由大量贯穿式柱状晶以及部分分布在柱状晶之间的亚晶粒组成。经测量，完全凝固后柱状晶的宽度在16 m，平均宽度为2.817 m。10um10m10m10 m(a)0.005 s(b)0.025s(c)0.075 s(d)0.125 s图4凝固不同时刻熔合线附近区域晶粒生长模拟结果Fig.4Simulation of crystal growth near fusion line at different times during solidification2.2熔覆层中溶质元素锯分布的模拟结果在熔覆层中垂直于晶轴方向对锯元素含量进行了提取。由图5可以看出，熔覆层中的锯元素质量分数在3.7 6%上下波动，但是多处锯元素质量分数超过了10%。可知在熔覆层凝固过程中，元素在这些区域产生了一定程度的偏析，导致这些区域析101MATERIALS FOR MECHANICAL ENGINEERING李继红，等Inconel625合金激光熔覆过程中显微组织演变的数值模拟2023机械工程材料60504020100020406080100120距离/um图5熔覆层中锯元素质量分数分布的模拟结果Fig.5Simulation of Nb element mass fraction distribution incladding layer出了大量富锯相2.3模拟结果的试验验证采取与模拟相同的工艺参数（激光功率2 kW，扫描速度5mms-1，光斑直径5mm)在ISDL-3008型激光熔覆设备上进行了单道激光熔覆Inconel625合金粉末试验。基体材料为Q345钢板，熔覆材料为Inconel625合金粉末，其化学成分如表2 所示，粒径在2 510 0 m。熔覆结束后制取金相试样，对金相试样表面打磨、抛光，用由1gFeCl3和10 mLHCl组成的溶液腐蚀后，采用GX71型光学显微镜（OM)和VEGA3XMU型扫描电子显微镜（SEM)对熔覆层不同区域的显微组织进行观察，采用NanoMeasurerl.2粒径分布计算软件测量柱状晶的宽度表2 Inconel625合金粉末的化学成分Table 2Chemical composition of Inconel625 alloy powder元素NiCrMoNbCuFeCoMnSd质量分数/%63.3521.879.53.760.320.750.060.390.02由图6 可以看出，熔覆层内部由柱状晶以及分布在柱状晶之间的亚晶粒组成，靠近熔合线区域由胞状晶组成，与模拟结果相吻合。试验测得柱状晶宽度在15m，平均宽度为2.7 43m。计算可知，平均宽度的模拟结果与试验结果间的相对误差仅为2.6 9%，说明所建CA-FD模型具有较高的可靠性。但是观察发现，柱状晶晶粒取向的模拟结果与试验结果略有不同。造成这种现象的可能原因是，在模拟熔合线附近熔池边缘柱状晶生长时，假设此处的温度梯度垂直于生长界面，而实际激光熔覆试验时，熔池边缘的温度梯度受很多外界条件的影响。2.4热输入对晶粒形貌影响的模拟结果采用所建模型进一步研究激光熔覆工艺参数对10m20m(a)内部(b)熔合线附近图6试验得到熔覆层不同区域的显微组织Fig.6Microstructures of different regions of cladding layer by test:(a)inside and(b)near fusion line熔覆层晶粒形貌的影响。在激光熔覆过程中，热输人主要体现在激光功率和扫描速度上，而扫描速度在二维尺度上很难体现，因此通过固定扫描速度、改变激光功率来研究热输入对熔覆层晶粒形貌的影响。根据前期研究成果，热输人分别设为3.6，4.0，4.5kJcm-l，扫描速度固定为5mms-1。由图7可以看出，随着激光热输入的增大，熔覆层中的晶粒数量呈减少趋势，柱状晶宽度呈增大趋势。这是由于热输人的增大会导致单位时间内熔池所获得的热量增多，使得熔池温度升高；而熔池温度越高，形核就越困难，凝固形成的熔覆层内的晶粒数量越少，熔覆层组织粗化2.5异质形核数量对晶粒形貌影响的模拟结果凝固组织与晶粒的形核和生长过程密切相关。在激光熔覆工艺确定的情况下，模拟了异质形核数量（1X104，110 5,110 个mm-）对熔覆层晶粒形貌的影响。由图8 可以看出，随着异质形核数量的增大，熔池在凝固过程中产生晶核的概率增大，熔覆层中柱状晶宽度减小，晶粒数量增多，熔覆层组织更加细小均匀，这与文献 2 5 的研究结果一致。因此，在实际激光熔覆过程中，可以通过向激光熔覆粉末中添加适量的异质形核剂（例如TiC、Zr C等高熔点碳化物细化剂 2 6 或Ce2Os、La 2 O 等高熔点稀土氧化物细化剂2 7）来细化晶粒，从而提高熔覆层的力学性能。3结论（1）采用建立的CA-FD模型模拟得到Inconel625合金激光熔覆层内部由柱状晶以及分布在柱状晶之间的亚晶粒组成，靠近熔合线区域由胞状晶组成，与试验结果相吻合；模拟得到和试验得到柱状晶的平均宽度分别为2.8 17，2.7 43m，二者之间的相对误102MATERIALS FOR MECHANICAL ENGINEERING李继红，等：I n c o n e l6 2 5合金激光熔覆过程中显微组织演变的数值模拟2023机械工程材料1000m1000um1000um(a)3.6kJcm-1(b)4.0kJcm-l(c)4.5kJcm-l图7不同热输入下熔覆层晶粒形貌模拟结果Fig.7Simulation of cladding layer grain morphology under different heat inputs1 000 m1000 m1000 um(a)1x104个mm-2(b)1x105个mm-2(c)1x106个mm-2图8不同异质形核数量下熔覆层晶粒形貌模拟结果Fig.8Simulation of cladding layer grain morphology under different heterogeneous nucleation quantities差仅为2.6 9%，验证了计算模型的可靠性。（2）在激光熔覆凝固过程中，熔合线附近晶粒首先形核并以胞状晶的形式生长，随后转变为柱状晶向熔池中心生长。随着激光热输入的减小或者异质形核数量的增加，熔覆层组织细化，晶粒数量增多，柱状晶宽度减小。参考文献：1 JAVID Y.Multi-response optimization in laser claddingprocess of WC powder on Inconel 718JJ.CIRP Journal ofManufacturing Science and Technology,2020,31:406-417.2谭金花，孙荣禄，牛伟，等.Ni60/h-BN含量对激光熔覆钛基复合涂层组织及性能的影响 J.表面技术，2 0 19,48（10）：10 7-115.TAN J H,SUN R L,NIU W,et al.Effect of Ni6o/h-BNcontent on microstructures and properties of laser claddingtitanium-based composite coatings JJ.Surface Technology,2019,48(10):107-115.3LIANG CJ,WANG C L,ZHANG K X,et al.Nucleation andstrengthening mechanism of laser cladding aluminum alloy byNi-Cr-B-Si alloy powder based on rare earth control J.Journal of Materials Processing Technology,2021,294:117145.4高秋实，闫华，秦阳，等.钛合金表面激光熔覆Ti-Ni十TiN十MoS2/TiS自润滑复合涂层 J.材料研究学报，2 0 18，32（12）：921-928.GAO Q S,YAN H,QIN Y,et al.Self-lubricating wearresistant composite coating Ti-Ni+TiN+MoS2/TiS preparedon Ti-6Al-4V alloy by laser claddingJ.Chinese Journal ofMaterials Research,2018,32(12):921-928.5LIU H F,TAN C K I,WEI Y F,et al.Laser-cladding andinterface evolutions of Inconel 625 alloy on low alloy stelsubstrate upon heat and chemical treatmentsJ.Surface andCoatings Technology,2020,404:126607.6VERDI D,GARRIDO M A,MUNEZ C J,et al.Microscaleeffect of high-temperature exposition on laser cladded Inconel625-Crg C2 metal matrix compositeJJ.Journal of Alloys andCompounds,2017,695:2696-2705.7ZHANG B C,BI G J,WANG P,et al.Microstructure andmechanical properties of Inconel 625/nano-TiB,compositefabricated by LAAMLJJ.Materials&Design,2016,1ll:70-79.8张敏，黄超，郭宇飞，等.TC4合金焊接熔池微观组织演变的数值模拟与分析 J.中国有色金属学报，2 0 2 0，30（8）：18 7 6-1886.ZHANG M,HUANG C,GUO Y F,et al.Numericalsimulation and analysis of microstructure evolution of TC4alloy weld poolJ.The Chinese Journal of Nonferrous Metals,2020,30(8):1876-1886.9郭钊，周建新，沈旭，等.改进元胞自动机法数值模拟高温合金凝固过程枝晶生长行为J.机械工程材料，2 0 2 0，44（2）：6 5-72.103欢迎来稿欢迎订阅欢迎刊登广告欢迎来稿欢迎订阅欢迎刊登广告订欢欢欢迎刊迎迎阅稿MATERIALSFOR MECHANICAL ENGINEERING李继红，等：I n c o n e l6 2 5合金激光熔覆过程中显微组织演变的数值模拟2023机械工程材料GUO Z,ZHOU J X,SHEN X,et al.Numerical simulation ofdendritic growth behavior in solidification of superalloy by animproved cellular automaton method J.Materials forMechanical Engineering,2020,44(2):65-72.10党元晓，祁文军，芦丽丽。激光熔覆技术数值模拟研究现状及发展趋势 J.热加工工艺，2 0 16，45（6）：2 3-2 7.DANG Y X,QI W J,LU L L.Research status anddevelopment trend of numerical simulation of laser claddingtechnologyJJ.Hot Working Technology,2016,45(6):23-27.11TAN W D,WEN S Y,BAILEY N,et al.Multiscale modelingof transport phenomena and dendritic growth in laser claddingprocesses JJ.Metallurgical and Materials Transactions B,2011,42(6):1306-1318.12SHIMONO Y,OBA M,NOMOTO S,et al.Numericalsimulation of solidification in additive manufacturing of Tialloy by multi-phase field methodC/Proceedings of the28th Annual International Solid Freeform FabricationSymposium.S.1.J:University of Texas at Austin,2017:1048-1057.13ACHARYA R,SHARON J A,STAROSELSKY A.Prediction of microstructure in laser powder bed fusionprocessJJ.Acta Materialia,2017,124:360-371.14孙道金，刘继常，李钦栋.激光熔覆纯镍熔池底部组织生长的相场法模拟J.中国激光，2 0 13，40（4）：0 40 30 0 5.SUN D J,LIU J C,LI Q D.Phase-field method simulation ofmicrostructure evolution at the bottom of melt pool in coaxiallaser claddingJ.Chinese Journal of Lasers,2013,40(4):0403005.15刘朝阳，齐欢.镍基单晶高温合金激光多道搭接熔覆过程晶体生长行为和微观组织分布的研究 J.电加工与模具，2 0 14(1):41-46.LIU Z Y,QI H.Study of crystal growth and microstructureformation in multi-track laser powder deposition of nickel-base single crystal superalloy J.Electromachining&Mould,2014(1):41-46.16PAVELIC V,TANBAKUCHI R,UYEHARA0A.Experiment and computed temperature histories in gastungsten-arc welding of thin plates JJ.Welding Journal,1969,48(7):295-305.17李露露.基于元胞自动机法的Fe-C合金焊缝熔池凝固过程微观组织模拟 D.西安：西安理工大学，2 0 16.LI L L.Simulation on solidification microstructure of Fe-Cbinary alloy in weld pool using cellular automaton methodD.Xian:Xian University of Technology,2016.18潘小勇.流体力学与传热学M.南昌：江西高校出版社，2019:362.PAN X Y.Fluid mechanics and heat transfer M.Nanchang:Jiangxi University and Colleges Press,2019:362.19RAPPAZ M,THEVOZ P H.Solute diffusion model forequiaxed dendritic growth:Analytical solution J.ActaMetallurgica,1987,35(12):2929-2933.20TSAI D C,HWANG W S.Numerical simulation ofsolidification morphologies of Cu-0.6Cr casting alloy usingmodified cellular automaton model J.Transactions ofNonferrous Metals Society of China,2010,20(6):1072-1077.21ZHU M,STEFANESCU D.Virtual front tracking model forthe quantitative modeling of dendritic growth in solidificationof alloysJ.Acta Materialia,2007,55(5):1741-1755.22邓小虎.金属热变形及焊缝凝固过程的元胞自动机模拟 D.大连：大连理工大学，2 0 0 9.DENG X H.Modeling of metal hot deformation and weldsolidification process using cellular automaton methodD.Dalian:Dalian University of Technology,2009.23LI D M,LI R,ZHANG P W.A cellular automaton techniquefor modelling of a binary dendritic growth with convectionJJ.Applied Mathematical Modelling,2007,31(6):971-982.24陶文铨.数值传热学 MI.西安：西安交通大学出版社，19 8 8：15-66.TAO W Q.Numerical heat transfer theoryMJ.Xian:XianJiaotong University Press,1988:15-66.25齐海波，张云浩，冯校飞，等。多元合金激光增材制造凝固组织演变模拟J.焊接学报，2 0 2 0，41(5）：7 1-7 7.QI H B,ZHANG Y H,FENG X F,et al.Simulation ofsolidification microstructure evolution in laser additionmanufacturing of multicomponent alloyJJ.Transactions ofthe China Welding Institution,2020,41(5):71-77.26周京.Inconel601H镍基合金焊缝化学法晶粒细化的研究D.天津：河北工业大学，2 0 15.ZHOU J.A study on chemical method of grain refinement ofInconel601H nickel based alloy weld D.Tianjin:HebeiUniversity of Technology,2015.27师帅，李长荣，刘然，等.Inconel600合金异质形核的经验电子理论分析 J.热加工工艺，2 0 2 0，49（2 2）：6 8-7 0.SHI S,LI C R,LIU R,et al.Empi