CN111339670B - 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法 - Google Patents

一种镍基合金激光熔覆过程中枝晶生长数值模拟方法 Download PDF

Info

Publication number
CN111339670B
CN111339670B CN202010130633.XA CN202010130633A CN111339670B CN 111339670 B CN111339670 B CN 111339670B CN 202010130633 A CN202010130633 A CN 202010130633A CN 111339670 B CN111339670 B CN 111339670B
Authority
CN
China
Prior art keywords
molten pool
nickel
laser cladding
based alloy
growth
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010130633.XA
Other languages
English (en)
Other versions
CN111339670A (zh
Inventor
张敏
郭宇飞
黄超
郭钊
张立胜
王刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Technology
Original Assignee
Xian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Technology filed Critical Xian University of Technology
Priority to CN202010130633.XA priority Critical patent/CN111339670B/zh
Publication of CN111339670A publication Critical patent/CN111339670A/zh
Application granted granted Critical
Publication of CN111339670B publication Critical patent/CN111339670B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Laser Beam Processing (AREA)

Abstract

本发明公开了一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,首先定义熔池形状并建立枝晶的形核与生长模型,然后建立溶质分配与扩散模型,最后编写计算机程序,输入合金热物性参数以及各种焊接工艺参数,导入模拟计算软件,进行计算即可得到模拟结果。本模型研究周期短、节能能源,能够模拟镍基合金熔覆层凝固过程中枝晶的生长形貌、溶质浓度分布以及形核基底数、表面换热系数、热输入等参数对枝晶生长的影响,从而对提高镍基合金熔覆层性能起到一定的指导作用。

Description

一种镍基合金激光熔覆过程中枝晶生长数值模拟方法
技术领域
本发明属于金属材料激光熔覆数值模拟技术领域,具体涉及一种镍基合金激光熔覆过程中枝晶生长数值模拟方法。
背景技术
激光熔覆能够改善工件表面的耐磨性、耐蚀性和抗氧化性等性能,具有可控性好、熔覆层组织致密并且结合强度高、易于实现自动化控制等优点。因此,已经成为了材料表面局部强化的一种有效途径。激光熔覆既能够提高材料表面的性能,使之达到我们使用的要求,又能够节约大量的贵重战略金属原料,既可用于材料表面改性,还可用于材料表面的修复,因此,激光熔覆应用前景十分广阔。
激光熔覆过程中会发生各种复杂的现象,传统的实验研究方法不仅费时费力,而且不能直观的再现激光熔覆的具体过程。元胞自动机法作为一种研究凝固过程中各种复杂变化现象的手段,可以清晰地再现激光熔覆过程中产生的各种复杂现象,能够便于我们研究熔覆层的组织结构,从而提高熔覆层的性能。
发明内容
本发明的目的是提供一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,具有帮助科研人员研究镍基合金熔覆层的微观组织结构的特点。
本发明所采用的技术方案是,一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,具体按照以下步骤实施:
步骤1:定义熔池形状;
步骤2:建立枝晶形核与生长模型;
步骤3:建立溶质分配与扩散模型;
步骤4:模拟计算及结果导出。
本发明的特点还在于:
步骤1按照以下具体步骤实施:
步骤1.1、建模时,将整个模拟区域分为两个部分,上半部分熔池以外的区域定步骤2.1义为空气,下半部分熔池以外的区域定义为母材;
步骤1.2、上半部分熔池形状可以通过以下公式建立:
步骤1.3、下半部分熔池形状可以通过以下公式建立:
式中:(i,j)为任意一点的坐标;(i1,j1)为上半部分熔池顶部的坐标;(i2,j2)为下半部分熔池底部的坐标。
步骤2按照以下具体步骤实施:
步骤2.1、在某一过冷度ΔT时所形成的晶粒密度n(ΔT)如下式所示:
dn/d(ΔT)的表达式如下:
式中:nmax为非均匀形核密度的最大值;ΔTσ为标准曲率过冷度;ΔTmax为最大形核过冷度;
步骤2.2、晶核形成后在过冷度的作用下会不断生长,在生长过程中排出的溶质浓度为:
式中:Dl为液相扩散系数;Δt为步长时间;dx为网格尺寸;
步骤2.3、随着晶粒的长大,液相不断地转变为固相,固相率的增长可以用下式进行计算:
式中:k0为溶质平衡分配系数;A为扰动因子;rand()能够在[0,1]产生一个随机数。
步骤3具体按照以下步骤实施:
枝晶在生长过程中排出的溶质会导致枝晶周围液相溶质浓度升高,使液相元胞间出现较大的浓度梯度,加剧溶质的扩散,对于二维非稳态溶质扩散,采用如下控制方程:
式中:Dl、Ds分别表示液相扩散系数和固相扩散系数。
步骤4具体按照以下步骤实施:
基于步骤1~3所构建的模型进行编程,将编好的程序导入模拟软件Matlab中,输入合金的热物性参数以及各种焊接工艺参数,进行计算即可得到激光熔覆过程中枝晶生长的模拟结果。
本发明的有益效果是:
(1)提出了一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,为研究镍基合金熔覆层组织提供了一种新的研究方法;
(2)能够得到镍基合金熔覆层组织凝固过程的可视化结果,有助于更加方便的研究镍基合金激光熔覆过程中的复杂转变;
(3)本发明研究周期短、节约资源。
附图说明
图1是本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法的流程图;
图2是本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法的熔池简化示意图;
图3是本发明实施例1模拟镍基合金熔覆层在不同时间下枝晶生长的形貌图;
图4是本发明实施例2模拟镍基合金熔覆层在不同形核基底数下枝晶生长的形貌图;
图5是本发明实施例3模拟镍基合金熔覆层在不同表面换热系数下枝晶生长的形貌图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
实施例1
本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,如图1所示,具体按照以下步骤实施:
步骤1:定义熔池形状;
步骤2:建立枝晶形核与生长模型;
步骤3:建立溶质分配与扩散模型;
步骤4:模拟计算及结果导出。
步骤1按照以下具体步骤实施:
步骤1.1、建模时,将整个模拟区域分为两个部分,上半部分熔池以外的区域定步骤2.1义为空气,下半部分熔池以外的区域定义为母材,如图2所示;
步骤1.2、上半部分熔池形状可以通过以下公式建立:
步骤1.3、下半部分熔池形状可以通过以下公式建立:
式中:(i,j)为任意一点的坐标;(i1,j1)为上半部分熔池顶部的坐标;(i2,j2)为下半部分熔池底部的坐标。
步骤2按照以下具体步骤实施:
本发明采用适用范围较广的准连续形核模型。度的作用下会不断生长,在生长过程中排出的溶质浓度为:
式中:Dl为液相扩散系数;Δt为步长时间;dx为网格尺寸;
步骤2.3、随着晶粒的长大,液相不断地转变为固相,固相率的增长可以用下式进行计算:
式中:k0为溶质平衡分配系数;A为扰动因子;rand()能够在[0,1]产生一个随机数。
步骤3具体按照以下步骤实施:
枝晶在生长过程中排出的溶质会导致枝晶周围液相溶质浓度升高,使液相元胞间出现较大的浓度梯度,这必然加剧溶质的扩散,对于二维非稳态溶质扩散,采用如下控制方程:
式中:Dl、Ds分别表示液相扩散系数和固相扩散系数。
步骤4具体按照以下步骤实施:
基于步骤1~3所构建的模型进行编程,将编好的程序导入模拟软件Matlab中,输入合金的热物性参数以及各种焊接工艺参数,进行计算即可得到激光熔覆过程中枝晶生长的模拟结果。
按照以上实施步骤,通过本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,可计算得到镍基合金熔覆层在不同时间下枝晶生长的模拟结果,如图3所示。
图3是形核基底数为1*105,表面换热系数为8*10-2w/(cm2*k),凝固时间分别为0.175s、0.25s、0.4s时熔池中枝晶的生长形貌。从图3中可以看出,随着时间的增长,液相不断地转变为固相,熔池底部因散热速度快,过冷度大,晶体以等轴晶或胞状晶的形式生长,而熔池中心的晶体因过冷度小都以柱状晶形式向熔池中心生长,由图中还可以看到,枝晶间竞争生长激烈,那些生长方向与热流方向平行的晶体,更具有生长优势。
实施例2
本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,如图1所示,具体按照以下步骤实施:
步骤1:定义熔池形状;
步骤2:建立枝晶形核与生长模型;
步骤3:建立溶质分配与扩散模型;
步骤4:模拟计算及结果导出。
步骤1按照以下具体步骤实施:
步骤1.1、建模时,将整个模拟区域分为两个部分,上半部分熔池以外的区域定步骤2.1义为空气,下半部分熔池以外的区域定义为母材,如图2所示;
步骤1.2、上半部分熔池形状可以通过以下公式建立:
步骤1.3、下半部分熔池形状可以通过以下公式建立:
式中:(i,j)为任意一点的坐标;(i1,j1)为上半部分熔池顶部的坐标;(i2,j2)为下半部分熔池底部的坐标。
步骤2按照以下具体步骤实施:
本发明采用适用范围较广的准连续形核模型。度的作用下会不断生长,在生长过程中排出的溶质浓度为:
式中:Dl为液相扩散系数;Δt为步长时间;dx为网格尺寸;
步骤2.3、随着晶粒的长大,液相不断地转变为固相,固相率的增长可以用下式进行计算:
式中:k0为溶质平衡分配系数;A为扰动因子;rand()能够在[0,1]产生一个随机数。
步骤3具体按照以下步骤实施:
枝晶在生长过程中排出的溶质会导致枝晶周围液相溶质浓度升高,使液相元胞间出现较大的浓度梯度,这必然加剧溶质的扩散,对于二维非稳态溶质扩散,采用如下控制方程:
式中:Dl、Ds分别表示液相扩散系数和固相扩散系数。
步骤4具体按照以下步骤实施:
基于步骤1~3所构建的模型进行编程,将编好的程序导入模拟软件Matlab中,输入合金的热物性参数以及各种焊接工艺参数,进行计算即可得到激光熔覆过程中枝晶生长的模拟结果。
按照以上实施步骤,通过本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,可计算得到镍基合金熔覆层在不同形核基底数下枝晶生长的模拟结果,如图4所示。
图4是表面换热系数为8*10-2w/(cm2*k),凝固时间为0.4s,形核基底数分别为1*104、1*105时熔池中枝晶的生长形貌。从图4中可以看出,随着形核基底数的增大,熔池中晶粒数量明显增多,熔池底部的等轴晶与熔池中心的柱状晶均发生了细化。
实施例3
本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,如图1所示,具体按照以下步骤实施:
步骤1:定义熔池形状;
步骤2:建立枝晶形核与生长模型;
步骤3:建立溶质分配与扩散模型;
步骤4:模拟计算及结果导出。
步骤1按照以下具体步骤实施:
步骤1.1、建模时,将整个模拟区域分为两个部分,上半部分熔池以外的区域定步骤2.1义为空气,下半部分熔池以外的区域定义为母材,如图2所示;
步骤1.2、上半部分熔池形状可以通过以下公式建立:
步骤1.3、下半部分熔池形状可以通过以下公式建立:
式中:(i,j)为任意一点的坐标;(i1,j1)为上半部分熔池顶部的坐标;(i2,j2)为下半部分熔池底部的坐标。
步骤2按照以下具体步骤实施:
本发明采用适用范围较广的准连续形核模型。度的作用下会不断生长,在生长过程中排出的溶质浓度为:
式中:Dl为液相扩散系数;Δt为步长时间;dx为网格尺寸;
步骤2.3、随着晶粒的长大,液相不断地转变为固相,固相率的增长可以用下式进行计算:
式中:k0为溶质平衡分配系数;A为扰动因子;rand()能够在[0,1]产生一个随机数。
步骤3具体按照以下步骤实施:
枝晶在生长过程中排出的溶质会导致枝晶周围液相溶质浓度升高,使液相元胞间出现较大的浓度梯度,这必然加剧溶质的扩散,对于二维非稳态溶质扩散,采用如下控制方程:
式中:Dl、Ds分别表示液相扩散系数和固相扩散系数。
步骤4具体按照以下步骤实施:
基于步骤1~3所构建的模型进行编程,将编好的程序导入模拟软件Matlab中,输入合金的热物性参数以及各种焊接工艺参数,进行计算即可得到激光熔覆过程中枝晶生长的模拟结果。
按照以上实施步骤,通过本发明一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,可计算得到镍基合金熔覆层在不同表面换热系数下枝晶生长的模拟结果,如图5所示。
图5是凝固时间为0.4s,形核基底数为1*105,表面换热系数分别为8*10-2w/(cm2*k)、5*10-2w/(cm2*k)时熔池中枝晶的生长形貌。从图5中可以看出,随着表面换热系数的减小,晶体生长速度减慢,这是因为表面换热系数减小会导致熔池中的热量不利于散失,造成过冷度减小所造成的。
本发明的有益效果是:针对激光熔覆过程中,熔池中枝晶的形核与生长问题,提供了一种镍基合金激光熔覆过程中枝晶生长数值模拟方法。本发明可以模拟镍基合金凝固过程中枝晶的生长形貌和溶质组元的分布,同时还可以模拟形核基底数、表面换热系数、热输入等参数对枝晶生长的影响,从而对提高镍基合金熔覆层性能起到一定的指导作用。

Claims (2)

1.一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,其特征在于,具体按照以下步骤实施:
步骤1:定义熔池形状;
步骤2:建立枝晶形核与生长模型;
步骤3:建立溶质分配与扩散模型;
步骤4:模拟计算及结果导出;
步骤1按照以下具体步骤实施:
步骤1.1、建模时,将整个模拟区域分为两个部分,上半部分熔池以外的区域定义为空气,下半部分熔池以外的区域定义为母材;
步骤1.2、上半部分熔池形状可以通过以下公式建立:
步骤1.3、下半部分熔池形状可以通过以下公式建立:
式中:(i,j)为任意一点的坐标;(i 1,j 1)为上半部分熔池顶部的坐标;(i 2,j 2)为下半部分熔池底部的坐标;
步骤2按照以下具体步骤实施:
步骤2.1、在某一过冷度ΔT时所形成的晶粒密度n(ΔT)如下式所示:
dn/d(ΔT)的表达式如下:
式中:n max为非均匀形核密度的最大值;为标准曲率过冷度;ΔT max为最大形核过冷度;
步骤2.2、晶核形成后在过冷度的作用下会不断生长,在生长过程中排出的溶质浓度为:
式中:Dl为液相扩散系数;Δt为步长时间;dx为网格尺寸;
步骤2.3、随着晶粒的长大,液相不断地转变为固相,固相率的增长可以用下式进行计算:
式中:k0为溶质平衡分配系数;A为扰动因子;rand()能够在[0,1]产生一个随机数;
步骤3具体按照以下步骤实施:
枝晶在生长过程中排出的溶质会导致枝晶周围液相溶质浓度升高,使液相元胞间出现较大的浓度梯度,加剧溶质的扩散,对于二维非稳态溶质扩散,采用如下控制方程:
式中:DlDs分别表示液相扩散系数和固相扩散系数。
2.根据权利要求1所述的一种镍基合金激光熔覆过程中枝晶生长数值模拟方法,其特征在于,所述步骤4具体按照以下步骤实施:
基于步骤1~3所构建的模型进行编程,将编好的程序导入模拟软件Matlab中,输入合金的热物性参数以及各种焊接工艺参数,进行计算即可得到激光熔覆过程中枝晶生长的模拟结果。
CN202010130633.XA 2020-02-28 2020-02-28 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法 Active CN111339670B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010130633.XA CN111339670B (zh) 2020-02-28 2020-02-28 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010130633.XA CN111339670B (zh) 2020-02-28 2020-02-28 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法

Publications (2)

Publication Number Publication Date
CN111339670A CN111339670A (zh) 2020-06-26
CN111339670B true CN111339670B (zh) 2024-04-09

Family

ID=71184222

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010130633.XA Active CN111339670B (zh) 2020-02-28 2020-02-28 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法

Country Status (1)

Country Link
CN (1) CN111339670B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112185474B (zh) * 2020-09-07 2024-02-06 西安理工大学 一种Ti-45%Al合金定向凝固过程的数值模拟方法
CN112115634B (zh) * 2020-09-21 2021-04-06 哈尔滨理工大学 金属液单向凝固过程中晶粒组织三维数值预测方法
CN116145133B (zh) * 2023-02-06 2024-08-30 中北大学 一种镍基单晶高温合金激光外延生长方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103668461A (zh) * 2013-09-21 2014-03-26 北京工业大学 一种镍基超合金Rene80定向生长柱晶及单晶合金制备及零部件制造方法
CN104014768A (zh) * 2014-06-24 2014-09-03 哈尔滨理工大学 一种镁合金枝晶组织数值模拟的方法
CN107309543A (zh) * 2017-05-04 2017-11-03 南京航空航天大学 一种激光焊接熔池枝晶生长模拟方法
CN110321604A (zh) * 2019-06-18 2019-10-11 西安理工大学 一种三元合金凝固过程单个枝晶生长数值模拟方法
CN110489821A (zh) * 2019-07-29 2019-11-22 西安理工大学 一种镍基合金堆焊熔池枝晶生长数值模拟方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103668461A (zh) * 2013-09-21 2014-03-26 北京工业大学 一种镍基超合金Rene80定向生长柱晶及单晶合金制备及零部件制造方法
CN104014768A (zh) * 2014-06-24 2014-09-03 哈尔滨理工大学 一种镁合金枝晶组织数值模拟的方法
CN107309543A (zh) * 2017-05-04 2017-11-03 南京航空航天大学 一种激光焊接熔池枝晶生长模拟方法
CN110321604A (zh) * 2019-06-18 2019-10-11 西安理工大学 一种三元合金凝固过程单个枝晶生长数值模拟方法
CN110489821A (zh) * 2019-07-29 2019-11-22 西安理工大学 一种镍基合金堆焊熔池枝晶生长数值模拟方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
镍基合金焊接熔池凝固组织模拟;董志波;马瑞;王勇;占小红;魏艳红;;焊接学报;20101025(第10期);31-34 *

Also Published As

Publication number Publication date
CN111339670A (zh) 2020-06-26

Similar Documents

Publication Publication Date Title
CN111339670B (zh) 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法
CN110489821B (zh) 一种镍基合金堆焊熔池枝晶生长数值模拟方法
CN111261241B (zh) 一种镍基合金多道堆焊过程中枝晶生长数值模拟方法
CN110489820B (zh) 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法
CN107790717A (zh) 一种实现镍基合金晶体学织构调控的准连续激光金属3d打印方法
Liang et al. Experimental optimization of laser additive manufacturing process of single-crystal nickel-base superalloys by a statistical experiment design method
CN110968954B (zh) 一种基于元胞自动机的bga锡铅焊球凝固过程模拟方法
CN101075287B (zh) 连铸坯凝固组织生长过程的物理模拟方法及其装置
CN107309543A (zh) 一种激光焊接熔池枝晶生长模拟方法
CN110489818B (zh) 一种三元合金焊接熔池柱状枝晶生长数值模拟方法
CN110321604B (zh) 一种三元合金凝固过程单个枝晶生长数值模拟方法
Zhou et al. Effects of process parameters on microstructure and cracking susceptibility of a single crystal superalloy fabricated by directed energy deposition
CN104404418B (zh) 一种镍基高温合金的热处理方法
CN110263418B (zh) 一种体心立方合金微观偏析数值预测方法
CN112185474B (zh) 一种Ti-45%Al合金定向凝固过程的数值模拟方法
CN110344049A (zh) 一种单晶/定向凝固镍基高温合金的修复方法及其应用
CN103350219A (zh) 一种降低铸锭宏观偏析的浇铸方法
CN111370078B (zh) 一种Inconel 625合金多层多道激光熔覆过程中晶体生长数值模拟方法
CN113192565A (zh) 一种钛铝合金定向凝固过程晶粒生长的三维数值模拟方法
CN116189787A (zh) 基于相场动力学的轻量化车用镁合金焊接熔池内cet转变模拟方法
CN111611734B (zh) 预测激光熔覆层顶部转向树枝晶区的有限元计算方法
CN113948158B (zh) 一种通过相场法高效模拟枝晶生长的方法
CN113823357B (zh) 一种四元合金凝固过程等轴晶生长数值模拟方法
CN118709357A (zh) Inconel 718合金电弧熔覆过程中相转变数值模拟方法
CN116759027A (zh) 镍基合金凝固时晶粒溶质浓度分布模拟方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant