CN110489818A - 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 - Google Patents
一种三元合金焊接熔池柱状枝晶生长数值模拟方法 Download PDFInfo
- Publication number
- CN110489818A CN110489818A CN201910689735.2A CN201910689735A CN110489818A CN 110489818 A CN110489818 A CN 110489818A CN 201910689735 A CN201910689735 A CN 201910689735A CN 110489818 A CN110489818 A CN 110489818A
- Authority
- CN
- China
- Prior art keywords
- cellular
- partalloy
- ternary alloy
- interface
- welding pool
- 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.)
- Granted
Links
Landscapes
- Manufacture And Refinement Of Metals (AREA)
- Arc Welding In General (AREA)
Abstract
本发明公开了一种三元合金焊接熔池柱状枝晶生长数值模拟方法,具体步骤如下:首先对焊接熔池凝固条件进行简化,然后建立枝晶的形核与生长模型,最后编写计算机程序,输入合金热物性参数以及各种焊接工艺参数,进行计算即可得到模拟结果。本模型能够模拟三元合金焊接熔池柱状枝晶的形核与生长、溶质浓度的分布,以及研究温度梯度、扰动振幅、各向异性强度等焊接工艺参数对枝晶生长的影响,从而对研究三元合金焊接熔池凝固过程中柱状枝晶的形核与生长起到一定的促进作用。
Description
技术领域
本发明属于金属材料焊接工艺数值模拟技术领域,具体涉及一种三元合金焊接熔池柱状枝晶生长数值模拟方法。
背景技术
焊接熔池在凝固过程中常常以柱状枝晶的形式生长,因此,长期以来对枝晶生长规律的研究受到国内外专家学者的广泛关注。然而影响焊接熔池微观组织的因素很多,仅仅通过实验的方法研究结晶过程不仅耗时、耗力和耗资,而且很难直接观察到枝晶生长的整个过程。因此,通过数值模拟技术研究凝固规律,预测显微组织具有巨大的研究潜力,近几十年来,通过材料领域科研工作者的共同努力,建立了各种各样的微观组织模拟方法,其中最著名的有相场法(PF)和元胞自动机法(CA)。虽然元胞自动机法运用于材料领域相对于相场法较晚,但其以高效和较强的工程应用能力而发展迅猛。
目前,常见的焊接熔池凝固组织的CA模型主要适用于二元合金,但在实际应用中大多数情况下使用的是三元合金,而对于三元合金焊接熔池凝固组织数值模拟的报道则较少,因此,建立一种三元合金焊接熔池柱状枝晶生长CA模型就尤为重要。
发明内容
本发明的目的是提供一种三元合金焊接熔池柱状枝晶生长数值模拟方法,解决了现有技术中存在缺少适用于三元合金焊接熔池柱状枝晶生长CA模型的问题。
本发明所采用的技术方案是,一种三元合金焊接熔池柱状枝晶生长数值模拟方法,具体按照以下步骤实施:
步骤1、简化三元合金焊接熔池凝固条件;
步骤2、建立三元合金焊接熔池柱状枝晶形核模型;
步骤3、定义捕获规则;
步骤4、建立三元合金焊接熔池柱状枝晶生长模型;
步骤5、模拟计算及结果导出。
本发明的特点还在于:
步骤1简化模型条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、元胞邻域关系采用V.Neumann型邻域,即四邻域;
简化条件3、三元合金凝固过程中设定溶质组元为B和C,忽略溶质之间的互扩散;
简化条件4、焊接熔池近似为对称分布,只考虑一半焊接熔池的凝固;
简化条件5、将熔池理想化为四分之一圆弧,圆弧内部定义为液态金属,外部定义为母材;
简化条件6、熔池中任意元胞到熔池中心的距离R由简化条件5可以得出:
式中:(i,j)为任意元胞的坐标,(s,s)为熔池中心的坐标。
步骤2具体按照以下步骤实施:
假定在时间步长Δt内,过冷度的增量为δ(ΔT),则单位体积(v)内形核数新增加的量δn(ΔT)为:
δn(ΔT)=n(ΔT+δT)-n(ΔT)
新晶粒形核的位置通过每一步长内形核概率确定,所述每一步长内形核概率的数值为:
pn=δn/n=δn(ΔT)VC
式中:pn为每一步长内形核概率,n为网格的单元总个数,Vc是网格单元的体积。
步骤3具体按照以下步骤实施:
步骤3.1、将模拟区域划分为正方形网格,每一个网格即为一个元胞;
步骤3.2、凝固开始时,基于步骤2的枝晶形核模型,熔池内会自动形成晶核元胞,晶核元胞周围的邻居元胞则被捕捉成为界面元胞,熔池内其余元胞均处于液相状态,即液相元胞;
步骤3.3、选定一个的晶核元胞,对该晶核元胞周围的界面元胞进行固相分数求解并判定,若界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞;
步骤3.4、对步骤2.3所得新的固相元胞周围的界面元胞进行固相分数求解并判定,界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞,以此类推,直至所有液相元胞转变为固相元胞。
步骤4具体按照以下步骤实施:
在稳态生长过程,枝晶在总过冷度的作用下开始凝固生长,总过冷度可由下式表示:
ΔT(tn)=Tl-T(tn)+ml1×(Cl1(tn)-C1)+ml2×(Cl2(tn)-C2)-Γ(θ)×k(tn)
式中:Tl为液相线温度;T(tn)为tn时刻液态金属的温度;ml1、ml2分别为溶质组元B、C的液相线斜率;Cl1(tn)、Cl2(tn)分别为溶质组元B、C在tn时刻的液相溶质浓度;C1、C2分别为溶质组元B、C的初始溶质浓度;Γ(θ)为Gibbs-Thompson系数;k(tn)为tn时刻的界面曲率;
枝晶在总过冷度的作用下不断凝固生长,当总过冷度给定后,界面液相溶质浓度可由下式得到:
在界面胞凝固时,溶质会向周围的液相排出,时间Δt内,穿过固/液界面面积Δx的溶质含量可表示为:
式中:Dli为i组元的液相扩散系数;jie为界面元胞邻近的液相元胞;Cjie为jie的溶质浓度;
金属凝固过程中伴随着液相份数的不断减少和固相份数的不断增加,在Δt时间内固相份数的增量可以有由下式:
凝固过程中,界面的推进速度与元胞固相份数的增量成正比,表示为:
在时间tn内,某一目标界面单元网格j内的固相份数为:
式中:N为迭代的次数。
步骤5具体按照以下步骤实施:
步骤5.1:将基于步骤1~4所构建的三元合金焊接熔池柱状枝晶生长模型导入模拟软件Matlab中;
步骤5.2:输入合金的热物性参数以及各种焊接工艺参数,进行计算得到三元合金焊接熔池中柱状枝晶生长的模拟结果。
本发明的优点在于:
(1)对于目前缺少三元合金凝固过程的数值模拟问题,提出了一种三元合金焊接熔池柱状枝晶生长的元胞自动机模型;
(2)本发明可以模拟凝固过程的动态变化,为进一步研究凝固机理提供了一种新的研究方法;
(3)相对于传统的试验方法,本发明即节省时间和资源,又环保节能。
附图说明
图1是本发明一种三元合金焊接熔池柱状枝晶生长数值模拟方法的流程图;
图2是本发明一种三元合金焊接熔池柱状枝晶生长数值模拟方法的四邻域元胞关系示意图;
图3是本发明一种三元合金焊接熔池柱状枝晶生长数值模拟方法的焊接熔池简化示意图;
图4a是实施例1Fe-0.8%C-0.3%Si三元合金C元素的固相溶质浓度分布状态和枝晶生长形貌图;
图4b是实施例1Fe-0.8%C-0.3%Si三元合金Si元素的固相溶质浓度分布状态和枝晶生长形貌图;
图5a是实施例1Fe-0.8%C-0.3%Si三元合金温度梯度为1*104K/m时,C元素的的固相溶质浓度分布状态和枝晶生长形貌图;
图5b是实施例1Fe-0.8%C-0.3%Si三元合金温度梯度为1*104K/m时,Si元素的的固相溶质浓度分布状态和枝晶生长形貌图;
图6a是实施例1Fe-0.8%C-0.3%Si三元合金扰动振幅为1时,C元素的固相溶质浓度分布状态和枝晶生长形貌图;
图6b是实施例1Fe-0.8%C-0.3%Si三元合金扰动振幅为1时,Si元素的固相溶质浓度分布状态和枝晶生长形貌图;
图7a是实施例1Fe-0.8%C-0.3%Si三元合金各向异性强度为0时,C元素的固相溶质浓度分布状态和枝晶生长形貌图;
图7b是实施例1Fe-0.8%C-0.3%Si三元合金各向异性强度为0时,Si元素的固相溶质浓度分布状态和枝晶生长形貌图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
实施例1
一种三元合金焊接熔池柱状枝晶生长数值模拟方法,如图1所示,具体按照以下步骤实施:
步骤1、简化三元合金焊接熔池凝固条件;
步骤2、建立三元合金焊接熔池柱状枝晶形核模型;
步骤3、定义捕获规则;
步骤4、建立三元合金焊接熔池柱状枝晶生长模型;
步骤5、模拟计算及结果导出。
步骤1简化模型条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、元胞邻域关系采用V.Neumann型邻域,即四邻域,如图2所示;
简化条件3、三元合金凝固过程中设定两溶质组元分别为B和C,忽略溶质之间的互扩散;
简化条件4、焊接熔池近似为对称分布,为了节省计算时间,只考虑一半焊接熔池的凝固;
简化条件5、将熔池理想化为四分之一圆弧,圆弧内部定义为液态金属,外部定义为母材,如图3所示;
简化条件6、熔池中任意元胞到熔池中心的距离R由简化条件5可以得出:
式中:(i,j)为任意元胞的坐标,(s,s)为熔池中心的坐标。
步骤2具体按照以下步骤实施:
在金属凝固的形核过程中为了实现形核在时间上的连续性,通常将过冷度的变化与时间联系在一起。在数值计算中,假定在时间步长Δt内,过冷度的增量为δ(ΔT),则单位体积(v)内形核数新增加的量δn(ΔT)为:
δn(ΔT)=n(ΔT+δT)-n(ΔT)
单位时间内所形成的晶核数δn为δn(ΔT)与液态金属整体体积V的乘积,而新晶粒形核的位置通过每一步长内形核概率确定,所述每一步长内形核概率的数值为:
pn=δn/n=δn(ΔT)VC
式中:pn为每一步长内形核概率,n为网格的单元总个数,Vc是网格单元的体积。
步骤3具体按照以下步骤实施:
步骤3.1、将模拟区域划分为正方形网格,每一个网格即为一个元胞;
步骤3.2、凝固开始时,基于步骤2的枝晶形核模型,熔池内会自动形成晶核元胞,晶核元胞周围的邻居元胞则被捕捉成为界面元胞,熔池内其余元胞均处于液相状态,即液相元胞;
步骤3.3、选定一个的晶核元胞,对该晶核元胞周围的界面元胞进行固相分数求解并判定,若界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞;
步骤3.4、对步骤2.3所得新的固相元胞周围的界面元胞进行固相分数求解并判定,界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞,以此类推,直至所有液相元胞转变为固相元胞。
步骤4具体按照以下步骤实施:
在稳态生长过程,枝晶在总过冷度的作用下开始凝固生长。枝晶外围的总过冷度主要包扩曲率过冷度、成分过冷度、热过冷度,总过冷度可由下式表示:
ΔT(tn)=Tl-T(tn)+ml1×(Cl1(tn)-C1)+ml2×(Cl2(tn)-C2)-Γ(θ)×k(tn)
式中:Tl为液相线温度;T(tn)为tn时刻液态金属的温度;ml1、ml2分别为溶质组元B、C的液相线斜率;Cl1(tn)、Cl2(tn)分别为溶质组元B、C在tn时刻的液相溶质浓度;C1、C2分别为溶质组元B、C的初始溶质浓度;Γ(θ)为Gibbs-Thompson系数;k(tn)为tn时刻的界面曲率;
枝晶在总过冷度的作用下不断凝固生长,当总过冷度给定后,界面液相溶质浓度可由下式得到:
在界面胞凝固时,溶质会向周围的液相排出,时间Δt内,穿过固/液界面面积Δx的溶质含量可表示为:
式中:Dli为i组元的液相扩散系数;jie为界面元胞邻近的液相元胞;Cjie为jie的溶质浓度;
金属凝固过程中伴随着液相份数的不断减少和固相份数的不断增加,在Δt时间内固相份数的增量可以有由下式:
凝固过程中,界面的推进速度与元胞固相份数的增量成正比,表示为:
在时间tn内,某一目标界面单元网格j内的固相份数为:
式中:N为迭代的次数。
步骤5具体按照以下步骤实施:
步骤5.1:将基于步骤1~4所构建的三元合金焊接熔池柱状枝晶生长模型编写成计算机程序导入模拟软件Matlab中;
步骤5.2:输入合金的热物性参数以及各种焊接工艺参数,如表1所示,进行计算得到Fe-0.8%C-0.3%Si三元合金焊接熔池中柱状枝晶生长的模拟结果,模拟结果如图4a~7b所示。
表1
表1为Fe-0.8%C-0.3%Si三元合金模拟时计算所用热物性参数。
图4a、4b是时间步长为6000,温度梯度为100K/m,扰动振幅为0.1,各向异性强度为0.3时,Fe-0.8%C-0.3%Si三元合金两组元的固相溶质浓度分布状态和枝晶生长形貌。
从图4a、4b中可以看出,柱状晶底部存在二次枝晶,而柱状晶顶部没有二次枝晶,这是由于枝晶在生长过程中存在着激烈的竞争生长所造成的。并且所有柱状晶都向着熔池中心生长。C组元的溶质浓度富集程度大于Si组元的溶质浓度富集程度。
图5a、5b是时间步长为6000,温度梯度为1*104K/m,扰动振幅为0.1,各向异性强度为0.3时,Fe-0.8%C-0.3%Si三元合金两组元的固相溶质浓度分布状态和枝晶生长形貌。
图5a、5b与图4a、4b对比可以看到,其他条件不变,随着温度梯度的增大,枝晶生长速度变慢,这是由于温度梯度增大,过冷度减小所造成的。
图6a、6b是时间步长为6000,温度梯度为100K/m,扰动振幅为1,各向异性强度为0.3时,Fe-0.8%C-0.3%Si三元合金两组元的固相溶质浓度分布状态和枝晶生长形貌。
图6a、6b与图4a、4b对比可以看到,其他条件不变,随着扰动振幅的增大,枝晶生长速度也有所减小,二次枝晶明显粗化,二次枝晶数量增大。
图7a、7b是时间步长为6000,温度梯度为100K/m,扰动振幅为0.1,各向异性强度为0时,Fe-0.8%C-0.3%Si三元合金两组元的固相溶质浓度分布状态和枝晶生长形貌。
图7a、7b与图4a、4b对比可以看到,其他条件不变,随着各向异性强度的减小,枝晶生长速度明显减小,一次枝晶生长方向多元化,二次枝晶显著增大,并且有三次枝晶形成。
实施例2
一种三元合金焊接熔池柱状枝晶生长数值模拟方法,如图1所示,具体按照以下步骤实施:
步骤1、简化一种三元合金焊接熔池凝固条件;
步骤2、建立一种三元合金焊接熔池柱状枝晶形核模型;
步骤3、定义捕获规则;
步骤4、建立三元合金焊接熔池柱状枝晶生长模型;
步骤5、模拟计算及结果导出。
步骤1简化模型条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、元胞邻域关系采用V.Neumann型邻域,即四邻域,如图2所示;
简化条件3、三元合金凝固过程中设定两溶质组元分别为B和C,忽略溶质之间的互扩散;
简化条件4、焊接熔池近似为对称分布,为了节省计算时间,只考虑一半焊接熔池的凝固;
简化条件5、将熔池理想化为四分之一圆弧,圆弧内部定义为液态金属,外部定义为母材,如图3所示;
简化条件6、熔池中任意元胞到熔池中心的距离R由简化条件5可以得出:
式中:(i,j)为任意元胞的坐标,(s,s)为熔池中心的坐标。
步骤2具体按照以下步骤实施:
在金属凝固的形核过程中为了实现形核在时间上的连续性,通常将过冷度的变化与时间联系在一起。在数值计算中,假定在时间步长Δt内,过冷度的增量为δ(ΔT),则单位体积(v)内形核数新增加的量δn(ΔT)为:
δn(ΔT)=n(ΔT+δT)-n(ΔT)
单位时间内所形成的晶核数δn为δn(ΔT)与液态金属整体体积V的乘积,而新晶粒形核的位置通过每一步长内形核概率确定,所述每一步长内形核概率的数值为:
pn=δn/n=δn(ΔT)VC
式中:pn为每一步长内形核概率,n为网格的单元总个数,Vc是网格单元的体积。
步骤3具体按照以下步骤实施:
步骤3.1、将模拟区域划分为正方形网格,每一个网格即为一个元胞;
步骤3.2、凝固开始时,基于步骤2的枝晶形核模型,熔池内会自动形成晶核元胞,晶核元胞周围的邻居元胞则被捕捉成为界面元胞,熔池内其余元胞均处于液相状态,即液相元胞;
步骤3.3、选定一个的晶核元胞,对该晶核元胞周围的界面元胞进行固相分数求解并判定,若界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞;
步骤3.4、对步骤2.3所得新的固相元胞周围的界面元胞进行固相分数求解并判定,界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞,以此类推,直至所有液相元胞转变为固相元胞。
步骤4具体按照以下步骤实施:
在稳态生长过程,枝晶在总过冷度的作用下开始凝固生长。枝晶外围的总过冷度主要包扩曲率过冷度、成分过冷度、热过冷度,总过冷度可由下式表示:
ΔT(tn)=Tl-T(tn)+ml1×(Cl1(tn)-C1)+ml2×(Cl2(tn)-C2)-Γ(θ)×k(tn)
式中:Tl为液相线温度;T(tn)为tn时刻液态金属的温度;ml1、ml2分别为溶质组元B、C的液相线斜率;Cl1(tn)、Cl2(tn)分别为溶质组元B、C在tn时刻的液相溶质浓度;C1、C2分别为溶质组元B、C的初始溶质浓度;Γ(θ)为Gibbs-Thompson系数;k(tn)为tn时刻的界面曲率;
枝晶在总过冷度的作用下不断凝固生长,当总过冷度给定后,界面液相溶质浓度可由下式得到:
在界面胞凝固时,溶质会向周围的液相排出,时间Δt内,穿过固/液界面面积Δx的溶质含量可表示为:
式中:Dli为i组元的液相扩散系数;jie为界面元胞邻近的液相元胞;Cjie为jie的溶质浓度;
金属凝固过程中伴随着液相份数的不断减少和固相份数的不断增加,在Δt时间内固相份数的增量可以有由下式:
凝固过程中,界面的推进速度与元胞固相份数的增量成正比,表示为:
在时间tn内,某一目标界面单元网格j内的固相份数为:
式中:N为迭代的次数。
步骤5具体按照以下步骤实施:
步骤5.1:将基于步骤1~4所构建的三元合金焊接熔池柱状枝晶生长模型编写成计算机程序导入模拟软件Matlab中;
步骤5.2:输入合金的热物性参数以及各种焊接工艺参数,如表2所示,进行计算得到Fe-0.5%C-0.02%P三元合金焊接熔池中柱状枝晶生长的模拟结果。
表2
表2为Fe-0.5%C-0.02%P三元合金模拟时计算所用热物性参数。
实施例3
一种三元合金焊接熔池柱状枝晶生长数值模拟方法,如图1所示,具体按照以下步骤实施:
步骤1、简化三元合金焊接熔池凝固条件;
步骤2、建立三元合金焊接熔池柱状枝晶形核模型;
步骤3、定义捕获规则;
步骤4、建立三元合金焊接熔池柱状枝晶生长模型;
步骤5、模拟计算及结果导出。
步骤1简化模型条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、元胞邻域关系采用V.Neumann型邻域,即四邻域,如图2所示;
简化条件3、三元合金凝固过程中设定两溶质组元分别为B和C,忽略溶质之间的互扩散;
简化条件4、焊接熔池近似为对称分布,为了节省计算时间,只考虑一半焊接熔池的凝固;
简化条件5、将熔池理想化为四分之一圆弧,圆弧内部定义为液态金属,外部定义为母材,如图3所示;
简化条件6、熔池中任意元胞到熔池中心的距离R由简化条件5可以得出:
式中:(i,j)为任意元胞的坐标,(s,s)为熔池中心的坐标。
步骤2具体按照以下步骤实施:
在金属凝固的形核过程中为了实现形核在时间上的连续性,通常将过冷度的变化与时间联系在一起。在数值计算中,假定在时间步长Δt内,过冷度的增量为δ(ΔT),则单位体积(v)内形核数新增加的量δn(ΔT)为:
δn(ΔT)=n(ΔT+δT)-n(ΔT)
单位时间内所形成的晶核数δn为δn(ΔT)与液态金属整体体积V的乘积,而新晶粒形核的位置通过每一步长内形核概率确定,所述每一步长内形核概率的数值为:
pn=δn/n=δn(ΔT)VC
式中:pn为每一步长内形核概率,n为网格的单元总个数,Vc是网格单元的体积。
步骤3具体按照以下步骤实施:
步骤3.1、将模拟区域划分为正方形网格,每一个网格即为一个元胞;
步骤3.2、凝固开始时,基于步骤2的枝晶形核模型,熔池内会自动形成晶核元胞,晶核元胞周围的邻居元胞则被捕捉成为界面元胞,熔池内其余元胞均处于液相状态,即液相元胞;
步骤3.3、选定一个的晶核元胞,对该晶核元胞周围的界面元胞进行固相分数求解并判定,若界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞;
步骤3.4、对步骤2.3所得新的固相元胞周围的界面元胞进行固相分数求解并判定,界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞,以此类推,直至所有液相元胞转变为固相元胞。
步骤4具体按照以下步骤实施:
在稳态生长过程,枝晶在总过冷度的作用下开始凝固生长。枝晶外围的总过冷度主要包扩曲率过冷度、成分过冷度、热过冷度,总过冷度可由下式表示:
ΔT(tn)=Tl-T(tn)+ml1×(Cl1(tn)-C1)+ml2×(Cl2(tn)-C2)-Γ(θ)×k(tn)
式中:Tl为液相线温度;T(tn)为tn时刻液态金属的温度;ml1、ml2分别为溶质组元B、C的液相线斜率;Cl1(tn)、Cl2(tn)分别为溶质组元B、C在tn时刻的液相溶质浓度;C1、C2分别为溶质组元B、C的初始溶质浓度;Γ(θ)为Gibbs-Thompson系数;k(tn)为tn时刻的界面曲率;
枝晶在总过冷度的作用下不断凝固生长,当总过冷度给定后,界面液相溶质浓度可由下式得到:
在界面胞凝固时,溶质会向周围的液相排出,时间Δt内,穿过固/液界面面积Δx的溶质含量可表示为:
式中:Dli为i组元的液相扩散系数;jie为界面元胞邻近的液相元胞;Cjie为jie的溶质浓度;
金属凝固过程中伴随着液相份数的不断减少和固相份数的不断增加,在Δt时间内固相份数的增量可以有由下式:
凝固过程中,界面的推进速度与元胞固相份数的增量成正比,表示为:
在时间tn内,某一目标界面单元网格j内的固相份数为:
式中:N为迭代的次数。
步骤5具体按照以下步骤实施:
步骤5.1:将基于步骤1~4所构建的三元合金焊接熔池柱状枝晶生长模型编写成计算机程序导入模拟软件Matlab中;
步骤5.2:输入合金的热物性参数以及各种焊接工艺参数,进行计算得到Ti-6%Al-4%V三元合金焊接熔池中柱状枝晶生长的模拟结果。
表4
表4为Ti-6%Al-4%V三元合金模拟时计算所用热物性参数。
本发明的有益效果是:本发明可以模拟三元合金焊接熔池凝固过程中柱状枝晶的形核与生长、溶质组元的分布状态,同时还可以模拟温度梯度、扰动振幅、各向异性强度等因素对枝晶生长的影响,从而对研究三元合金焊接熔池凝固过程中柱状枝晶的形核与生长起到一定的促进作用。
Claims (6)
1.一种三元合金焊接熔池柱状枝晶生长数值模拟方法,其特征在于,具体按照以下步骤实施:
步骤1、简化三元合金焊接熔池凝固条件;
步骤2、建立三元合金焊接熔池柱状枝晶形核模型;
步骤3、定义捕获规则;
步骤4、建立三元合金焊接熔池柱状枝晶生长模型;
步骤5、模拟计算及结果导出。
2.根据权利要求1所述的一种三元合金焊接熔池柱状枝晶生长数值模拟方法,其特征在于,所述步骤1简化模型条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、元胞邻域关系采用V.Neumann型邻域,即四邻域;
简化条件3、三元合金凝固过程中设定溶质组元为B和C,忽略溶质之间的互扩散;
简化条件4、焊接熔池近似为对称分布,为了节省计算时间,只考虑一半焊接熔池的凝固;
简化条件5、将熔池理想化为四分之一圆弧,圆弧内部定义为液态金属,外部定义为母材;
简化条件6、熔池中任意元胞到熔池中心的距离R由简化条件5可以得出:
式中:(i,j)为任意元胞的坐标,(s,s)为熔池中心的坐标。
3.根据权利要求1所述的一种三元合金焊接熔池柱状枝晶生长数值模拟方法,其特征在于,步骤2具体按照以下步骤实施:
假定在时间步长Δt内,过冷度的增量为δ(ΔT),则单位体积(v)内形核数新增加的量δn(ΔT)为:
δn(ΔT)=n(ΔT+δT)-n(ΔT)
新晶粒形核的位置通过每一步长内形核概率确定,所述每一步长内形核概率的数值为:
pn=δn/n=δn(ΔT)VC
式中:pn为每一步长内形核概率,n为网格的单元总个数,Vc是网格单元的体积。
4.根据权利要求1所述的一种三元合金焊接熔池柱状枝晶生长数值模拟方法,其特征在于,步骤3具体按照以下步骤实施:
步骤3.1、将模拟区域划分为正方形网格,每一个网格即为一个元胞;
步骤3.2、凝固开始时,基于步骤2的枝晶形核模型,熔池内会自动形成晶核元胞,晶核元胞周围的邻居元胞则被捕捉成为界面元胞,熔池内其余元胞均处于液相状态,即液相元胞;
步骤3.3、选定一个的晶核元胞,对该晶核元胞周围的界面元胞进行固相分数求解并判定,若界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞;
步骤3.4、对步骤2.3所得新的固相元胞周围的界面元胞进行固相分数求解并判定,界面元胞的固相分数大于1,则该界面元胞转变为固相元胞,新转变的固相元胞周围的液相元胞则被捕捉成为新的界面元胞,以此类推,直至所有液相元胞转变为固相元胞。
5.根据权利要求1所述的一种三元合金焊接熔池柱状枝晶生长数值模拟方法,其特征在于,步骤4具体按照以下步骤实施:
在稳态生长过程,枝晶在总过冷度的作用下开始凝固生长,所述总过冷度可由下式表示:
ΔT(tn)=Tl-T(tn)+ml1×(Cl1(tn)-C1)+ml2×(Cl2(tn)-C2)-Γ(θ)×k(tn)
式中:Tl为液相线温度;T(tn)为tn时刻液态金属的温度;ml1、ml2分别为溶质组元B、C的液相线斜率;Cl1(tn)、Cl2(tn)分别为溶质组元B、C在tn时刻的液相溶质浓度;C1、C2分别为溶质组元B、C的初始溶质浓度;Γ(θ)为Gibbs-Thompson系数;k(tn)为tn时刻的界面曲率;
枝晶在总过冷度的作用下不断凝固生长,当总过冷度给定后,界面液相溶质浓度可由下式得到:
在界面胞凝固时,溶质会向周围的液相排出,时间Δt内,穿过固/液界面面积Δx的溶质含量可表示为:
式中:Dli为i组元的液相扩散系数;jie为界面元胞邻近的液相元胞;Cjie为jie的溶质浓度;
金属凝固过程中伴随着液相份数的不断减少和固相份数的不断增加,在Δt时间内固相份数的增量可以有由下式:
凝固过程中,界面的推进速度与元胞固相份数的增量成正比,表示为:
在时间tn内,某一目标界面单元网格j内的固相份数为:
式中:N为迭代的次数。
6.权利要求1所述的一种三元合金焊接熔池柱状枝晶生长数值模拟方法,其特征在于,所述步骤5体按照以下步骤实施:
步骤5.1:将基于步骤1~4所构建的三元合金焊接熔池柱状枝晶生长模型导入模拟软件Matlab中;
步骤5.2:输入合金的热物性参数以及各种焊接工艺参数,进行计算得到三元合金焊接熔池中柱状枝晶生长的模拟结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910689735.2A CN110489818B (zh) | 2019-07-29 | 2019-07-29 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910689735.2A CN110489818B (zh) | 2019-07-29 | 2019-07-29 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110489818A true CN110489818A (zh) | 2019-11-22 |
CN110489818B CN110489818B (zh) | 2022-10-21 |
Family
ID=68548550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910689735.2A Active CN110489818B (zh) | 2019-07-29 | 2019-07-29 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110489818B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259585A (zh) * | 2020-01-14 | 2020-06-09 | 大连交通大学 | 一种焊缝晶粒形核生长的仿真方法及系统 |
CN111681714A (zh) * | 2020-07-02 | 2020-09-18 | 兰州大学 | 一种在定向凝固包晶合金中生长非典型三次枝晶的方法 |
CN112185474A (zh) * | 2020-09-07 | 2021-01-05 | 西安理工大学 | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 |
CN113139254A (zh) * | 2021-05-14 | 2021-07-20 | 苏州大学 | 描述6xxx系铝合金熔化焊熔池中柱状晶形貌的方法 |
CN113284570A (zh) * | 2021-05-25 | 2021-08-20 | 西安理工大学 | 一种铝合金焊接熔池微观组织的模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750425A (zh) * | 2012-07-17 | 2012-10-24 | 哈尔滨工业大学 | 一种焊接过程热影响区组织演变的模拟方法 |
US20170147723A1 (en) * | 2015-11-20 | 2017-05-25 | Metal Industries Research & Development Centre | Method of simulatively predicting a metal solidification microstructure for a continuous casting process |
CN107309543A (zh) * | 2017-05-04 | 2017-11-03 | 南京航空航天大学 | 一种激光焊接熔池枝晶生长模拟方法 |
-
2019
- 2019-07-29 CN CN201910689735.2A patent/CN110489818B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750425A (zh) * | 2012-07-17 | 2012-10-24 | 哈尔滨工业大学 | 一种焊接过程热影响区组织演变的模拟方法 |
US20170147723A1 (en) * | 2015-11-20 | 2017-05-25 | Metal Industries Research & Development Centre | Method of simulatively predicting a metal solidification microstructure for a continuous casting process |
CN107309543A (zh) * | 2017-05-04 | 2017-11-03 | 南京航空航天大学 | 一种激光焊接熔池枝晶生长模拟方法 |
Non-Patent Citations (3)
Title |
---|
张敏等: "Fe-C合金焊接熔池凝固过程CET转变的数值模拟", 《焊接学报》 * |
张敏等: "Ti-Al合金凝固过程中的枝晶生长研究", 《兵器材料科学与工程》 * |
魏艳红等: "基于CA模拟焊缝凝固过程枝晶生长的分析", 《焊接学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259585A (zh) * | 2020-01-14 | 2020-06-09 | 大连交通大学 | 一种焊缝晶粒形核生长的仿真方法及系统 |
CN111681714A (zh) * | 2020-07-02 | 2020-09-18 | 兰州大学 | 一种在定向凝固包晶合金中生长非典型三次枝晶的方法 |
CN111681714B (zh) * | 2020-07-02 | 2023-06-20 | 兰州大学 | 一种在定向凝固包晶合金中生长非典型三次枝晶的方法 |
CN112185474A (zh) * | 2020-09-07 | 2021-01-05 | 西安理工大学 | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 |
CN112185474B (zh) * | 2020-09-07 | 2024-02-06 | 西安理工大学 | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 |
CN113139254A (zh) * | 2021-05-14 | 2021-07-20 | 苏州大学 | 描述6xxx系铝合金熔化焊熔池中柱状晶形貌的方法 |
CN113139254B (zh) * | 2021-05-14 | 2023-01-17 | 苏州大学 | 描述6xxx系铝合金熔化焊熔池中柱状晶形貌的方法 |
CN113284570A (zh) * | 2021-05-25 | 2021-08-20 | 西安理工大学 | 一种铝合金焊接熔池微观组织的模拟方法 |
CN113284570B (zh) * | 2021-05-25 | 2024-03-19 | 西安理工大学 | 一种铝合金焊接熔池微观组织的模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110489818B (zh) | 2022-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110489818A (zh) | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 | |
CN102289542B (zh) | 一种铸件宏观偏析数值模拟的方法 | |
CN106407623B (zh) | 基于相场法的瞬态下焊接过程微观组织演变模拟方法 | |
CN110489821B (zh) | 一种镍基合金堆焊熔池枝晶生长数值模拟方法 | |
CN104014768B (zh) | 一种镁合金枝晶组织数值模拟的方法 | |
CN107309543A (zh) | 一种激光焊接熔池枝晶生长模拟方法 | |
CN110968954B (zh) | 一种基于元胞自动机的bga锡铅焊球凝固过程模拟方法 | |
CN110245449B (zh) | 一种镁合金铸造件成分不均匀性数值预测方法 | |
CN110489820A (zh) | 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法 | |
CN110263418B (zh) | 一种体心立方合金微观偏析数值预测方法 | |
Liu et al. | A cellular automaton-lattice Boltzmann method for modeling growth and settlement of the dendrites for Al-4.7% Cu solidification | |
CN103729511A (zh) | 一种复杂结构铸件铸造过程中成分偏析程度的预测方法 | |
CN111339670B (zh) | 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法 | |
Xu et al. | Multiscale modeling and simulation of directional solidification process of turbine blade casting with MCA method | |
CN116230142A (zh) | 一种铝合金凝固动力学过程的介观尺度预测方法 | |
CN101767189A (zh) | 一种钢锭中固相移动的模拟方法 | |
Wang et al. | Numerical modeling of equiaxed crystal growth in solidification of binary alloys using a lattice Boltzmann-finite volume scheme | |
CN113284570B (zh) | 一种铝合金焊接熔池微观组织的模拟方法 | |
CN112185474B (zh) | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 | |
Guo et al. | Competitive growth of diverging columnar grains during directional solidification: A three-dimensional phase-field study | |
Tsai et al. | A three dimensional cellular automaton model for the prediction of solidification morphologies of brass alloy by horizontal continuous casting and its experimental verification | |
Wang et al. | Numerical simulation on dendritic growth of Al–Cu alloy under convection based on the cellular automaton lattice Boltzmann method | |
CN113192565A (zh) | 一种钛铝合金定向凝固过程晶粒生长的三维数值模拟方法 | |
CN110970095B (zh) | 一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法 | |
ZHU et al. | Numerical simulation of recalescence of 3-dimensional isothermal solidification for binary alloy using phase-field approach |
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 |