CN113284570B - 一种铝合金焊接熔池微观组织的模拟方法 - Google Patents
一种铝合金焊接熔池微观组织的模拟方法 Download PDFInfo
- Publication number
- CN113284570B CN113284570B CN202110571832.9A CN202110571832A CN113284570B CN 113284570 B CN113284570 B CN 113284570B CN 202110571832 A CN202110571832 A CN 202110571832A CN 113284570 B CN113284570 B CN 113284570B
- Authority
- CN
- China
- Prior art keywords
- aluminum alloy
- supercooling
- welding
- temperature
- nucleation
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000003466 welding Methods 0.000 title claims abstract description 66
- 229910000838 Al alloy Inorganic materials 0.000 title claims abstract description 56
- 238000004088 simulation Methods 0.000 title claims abstract description 43
- 230000008569 process Effects 0.000 claims abstract description 47
- 230000006911 nucleation Effects 0.000 claims abstract description 28
- 238000010899 nucleation Methods 0.000 claims abstract description 28
- 238000007711 solidification Methods 0.000 claims abstract description 11
- 230000008023 solidification Effects 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000012546 transfer Methods 0.000 claims abstract description 4
- 238000004781 supercooling Methods 0.000 claims description 41
- 239000007791 liquid phase Substances 0.000 claims description 22
- 239000007790 solid phase Substances 0.000 claims description 9
- 239000007788 liquid Substances 0.000 claims description 7
- 238000009792 diffusion process Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 238000005192 partition Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 239000002904 solvent Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- 239000013078 crystal Substances 0.000 abstract description 14
- 230000009466 transformation Effects 0.000 abstract description 6
- 238000001816 cooling Methods 0.000 abstract description 4
- 210000004027 cell Anatomy 0.000 description 30
- 238000011160 research Methods 0.000 description 5
- 239000000126 substance Substances 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 229910001385 heavy metal Inorganic materials 0.000 description 2
- 229910045601 alloy Inorganic materials 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 210000001787 dendrite Anatomy 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Arc Welding In General (AREA)
Abstract
本发明公开了一种铝合金焊接熔池微观组织的模拟方法,具体步骤如下:步骤1、简化模型条件;步骤2、基于传热原理,构建宏观温度场模型;步骤3、基于插值原理,构建微观温度场模型;步骤4、基于凝固理论,构建晶粒形核与生长模型;步骤5、模拟计算及结果导出。本模型能够模拟铝合金焊接过程中晶粒的形核与生长过程,熔池内部发生的CET(柱状晶‑等轴晶)转变过程,同时还能够模拟焊接工艺参数以及冷却条件对熔池内部晶粒形貌的影响,从而为研究铝合金的焊接过程以及优化铝合金的焊接工艺提供了一种全新的方法。
Description
技术领域
本发明属于金属材料焊接工艺数值模拟技术领域,具体涉及一种铝合金焊接熔池微观组织的模拟方法。
背景技术
近年来,随着全球工业化进程的不断发展,环境污染和能源短缺问题日益突出,在满足使用要求的条件下,如果能够采用轻金属代替重金属,就可以使结构轻量化,减少重金属污染。铝合金因其质量轻、密度小、强度高、以及加工性能好等优点被广泛应用于航空航天、交通运输以及民用建筑等领域。铝合金的加工制造技术主要以铸造、冲压和焊接最为常见,在实际焊接中,铝合金的焊接方式主要是熔化焊,熔化焊过程中输入到焊件上的热量以及热量的分布对焊接接头产生重要影响,直接影响接头的微观组织和力学性能,传统的铝合金焊接过程研究常采用实验法,这种方法存在费时费力、经济效益低的缺点,而且焊接过程的热量具有瞬时、集中、动态的特征,传统的方法难以实时观察焊接熔池的温度场变化,进而无法实时观察熔池内晶粒组织的演变过程。
随着计算机技术的发展,数值模拟方法的出现为研究铝合金的焊接过程提供了一种新的思路,其中,元胞自动机法因其具有明确的物理基础以及先进的概率性思想,在各个学科中都得到了长足的发展。可见,采用元胞自动机法模拟铝合金焊接过程中微观组织演变情况是切实可行的,这将为研究铝合金的焊接过程提供一种全新的研究方法,因此,建立一种铝合金焊接熔池微观组织的模拟方法显得尤为重要。
发明内容
本发明的目的是提供一种铝合金焊接熔池微观组织的模拟方法,解决了当前研究中存在有关铝合金焊接过程中缺少的数值模拟方法的问题。
本发明所采用的技术方案是,一种铝合金焊接熔池微观组织的模拟方法,具体按照以下步骤实施:
步骤1、简化模型条件;
步骤2、基于传热原理,构建宏观温度场模型;
步骤3、基于插值原理,构建微观温度场模型;
步骤4、基于凝固理论,构建晶粒形核与生长模型;
步骤5、模拟计算及结果导出。
本发明的特点还在于:
步骤1简化模型条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、将焊接熔池简化为规则的半圆形状;
简化条件3、模拟过程中忽略动力学过冷,只考虑温度过冷、成分过冷和曲率过冷;
简化条件4、模拟区域划分为正方形单元,每一个单元即为一个元胞;
简化条件5、元胞邻域关系采用V.Neumann型邻域,即四邻域。
步骤2具体按照以下步骤实施:
模拟时选取高斯分布热源作为焊接热源,该热源模型由下式表示:
式中:r为焊件上任一点到热源中心的距离/m;rh为电弧有效热源半径/m;U为焊接电压/V;I为焊接电流/A;η为焊接热效率。
热源加载到焊件上后,在焊件的表面和内部进行传导,热传导方程由下式表示:
式中:T为温度场函数值/K;ρ为密度/Kg·m-3;Cρ为容积比热容/J·(m3·K)-1;λ为热导率/W·(m·K)-1。
确定导热的初始条件和边界条件,初始条件由下式确定:
T(x,y)=Tf (3);
式中:Tf为模拟区域初始温度/K,Tf为室温293K。
边界条件由下式确定:
qw=hw(T-Tf) (4);
式中:T为试件表面温度/K;Tf为周围介质温度/K;hw为换热系数/W/(m2*K)。
步骤3具体按照以下步骤实施:
假定材料在各个方向的热导率相同,根据插值原理,将(2)式转化为:
式中:Δx为网格尺寸;Δt为时间步长;(i,j)为微观元胞的坐标;Ti,j为微观元胞(i,j)的当前温度值/K;T′i,j为微观元胞(i,j)的下一时刻温度值/K。
步骤4具体按照以下步骤实施:
晶粒形核模型:
考虑形核的过冷度和形核的连续性,采用基于高斯分布函数的准连续形核模型,该模型由下式表示:
式中:Nmax为非均匀形核密度的最大值/m3;ΔTθ为曲率过冷度/℃;ΔTmax为最大形核过冷度/℃。
晶粒生长模型:
过冷度是晶粒生长的推动力,过冷度ΔT包括三个部分,由热过冷ΔTw、成分过冷ΔTc和曲率过冷ΔTθ组成,ΔTθ由下式表示:
ΔT=ΔTw+ΔTc+ΔTθ (7);
由吉布斯汤普森关系可得t时刻的过冷度为:
ΔT(t)=TL-T(t)+mL(CL(t)-C0)-ΓK(t) (8);
式中:TL为液相线温度/K;T(t)为t时刻的瞬时温度/K;mL为液相线斜率;C0为初始溶质浓度/wt.%;CL(t)为t时刻的液相溶质浓度/wt.%;Γ为Gibbs-Thompson系数;K(t)为t时刻的界面平均曲率;
界面元胞在凝固过程中会向周围的液相元胞排出多余溶质,单位时间Δt内排出的溶质量ΔC由下式计算:
式中:DL为溶质液相扩散系数/m2/s;CL为固液界面处液相溶质浓度/wt.%;nb表示界面元胞周围的液相元胞,Cnb为其溶质浓度/wt.%;
固/液界面推进过程中,界面元胞排出ΔC的溶质之后,元胞单元固相率会发生变化,固相率增量由下式表示:
式中:k0为溶质分配系数;A为扰动因子;rand为0到1之间的随机数;
采用如下公式(11)对液相元胞进行溶质扩散计算:
步骤5具体按照以下步骤实施:
步骤5.1:基于步骤1~4所构建的一种铝合金焊接熔池微观组织的数值模型进行编程;
步骤5.2:将编好的程序导入模拟软件Matlab中,输入相关铝合金的热物性参数和焊接工艺参数,即得到铝合金焊接熔池微观组织的数值模拟结果。
本发明的有益效果是:
(1)本发明提供一种铝合金焊接熔池微观组织的模拟方法,解决了现有技术中存在缺少的关于铝合金焊接过程的数值模型问题;
(2)本发明可以再现铝合金不同工艺参数下,焊接过程中任意时刻熔池内晶粒生长情况,模拟结果更加精确和直观,从而为进一步研究铝合金焊接过程提供了新的研究方案;
(3)相对于传统的通过金相、扫描的实验手段研究铝合金的焊接过程及其焊缝微观组织的转变过程,本发明基于物理学基础构建了数学模型,通过编程利用计算机模拟平台对铝合金的焊接过程进行了数值模拟,节省了研究成本,经济高效且节能环保。
附图说明
图1是本发明一种铝合金焊接熔池微观组织的模拟方法的流程图;
图2是本发明一种铝合金焊接熔池微观组织的模拟方法元胞邻域关系示意图;
图3(a)、(b)是本发明一种铝合金焊接熔池微观组织的模拟方法实施例1 1060铝合金焊接熔池边缘晶粒的形核与生长过程模拟结果;
图4(a)、(b)是本发明一种铝合金焊接熔池微观组织的模拟方法实施例2 1060铝合金焊接熔池中心发生的CET(柱状晶-等轴晶)转变的模拟结果;
图5(a)、(b)是本发明一种铝合金焊接熔池微观组织的模拟方法实施例3 5A06铝合金不同焊接工艺参数(热输入)下熔池内微观组织的模拟结果;
图6(a)、(b)是本发明一种铝合金焊接熔池微观组织的模拟方法实施例4 5A06铝合金不同冷却工艺参数(扰动振幅)下熔池内微观组织的模拟结果。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
如图1所示,为本发明一种铝合金焊接熔池微观组织的模拟方法的流程图,具体按照以下步骤实施:
步骤1、简化模型建立条件:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、将焊接熔池简化为规则的半圆形状;
简化条件3、模拟过程中忽略动力学过冷,只考虑温度过冷、成分过冷和曲率过冷;
简化条件4、模拟区域划分为正方形单元,每一个单元即为一个元胞;
简化条件5、元胞邻域关系采用V.Neumann型邻域,即四邻域,如图2所示。
步骤2、基于传热原理,构建宏观温度场模型:
焊接时,电弧热源作用在焊件上的区域称为加热斑点,加热斑点上的热流密度分布近似的用高斯函数来描述,如下式所示:
式中:r为焊件上任一点到热源中心的距离/m;rh为电弧有效热源半径/m;U为焊接电压/V;I为焊接电流/A;η为焊接热效率。
焊接热量作用在焊件上后会在焊件的表面和内部传导,二维瞬态导热方程如下式所示:
式中:T为温度场函数值/K;ρ为密度/Kg·m-3;Cρ为容积比热容/J·(m3·K)-1;λ为热导率/W·(m·K)-1。
求解导热问题之前,需要确定导热的初始条件和边界条件,初始条件可由下式确定:
T(x,y)=Tf (3);
式中:Tf为模拟区域初始温度/K;Tf为室温293K。
边界条件是指物体表面与周围介质的热交换情况,边界条件可由下式确定:
qw=hw(T-Tf) (4);
式中:T为试件表面温度/K;Tf为周围介质温度/K;hw为换热系数/W/(m2*K)。
步骤3、基于插值原理,构建微观温度场模型:
温度场模型属于宏观范畴,而晶粒生长属于微观范畴,因此,需要进行温度场的宏微观耦合,由于所建模型中设定微观元胞为规则正方形,且模拟区域内部无其他热源,根据插值原理,可将(2)式转化为:
式中:Δx为网格尺寸;Δt为时间步长;(i,j)为微观元胞的坐标;Ti,j为微观元胞(i,j)的当前温度值/K;T′i,j为微观元胞(i,j)的下一时刻温度值/K。
计算得到宏观单元的温度数据后,基于上式,通过插值原理求解微观单元的温度值,便可得到模拟区域内部元胞、表面元胞的温度数据,从而实现温度场宏微观耦合过程。
步骤4、基于凝固理论,构建晶粒形核与生长模型:
晶粒形核模型:
考虑形核的过冷度和形核的连续性,采用基于高斯分布函数的准连续形核模型,该模型可由下式表示:
式中:Nmax为非均匀形核密度的最大值/m3;ΔTθ为曲率过冷度/K;ΔTmax为最大形核过冷度/K。
晶粒生长模型:
根据热力学定律,过冷度是晶粒生长的推动力,晶粒形核后需在一定过冷度作用下才能持续生长,过冷度ΔT包括三个部分,由热过冷ΔTw、成分过冷ΔTc和曲率过冷ΔTθ组成,ΔTθ可由下式表示:
ΔT=ΔTw+ΔTc+ΔTθ (7);
由吉布斯汤普森关系可得t时刻的过冷度为:
ΔT(t)=TL-T(t)+mL(CL(t)-C0)-ΓK(t) (8);
式中:TL为液相线温度/K;T(t)为t时刻的瞬时温度/K;mL为液相线斜率;C0为初始溶质浓度/wt.%;CL(t)为t时刻的液相溶质浓度/wt.%;Γ为Gibbs-Thompson系数;K(t)为t时刻的界面平均曲率。
界面元胞在凝固过程中会向周围的液相元胞排出多余溶质,单位时间Δt内排出的溶质量ΔC可由下式计算:
式中:DL为溶质液相扩散系数/m2/s;CL为固液界面处液相溶质浓度/wt.%;nb表示界面元胞周围的液相元胞,Cnb为其溶质浓度/wt.%。
固/液界面推进过程中,界面元胞排出ΔC的溶质之后,元胞单元固相率会发生变化,固相率增量可由下式表示:
式中:k0为溶质分配系数;A为扰动因子;rand为0到1之间的随机数。
界面生长过程中,多余的溶质排到周围液相邻胞中会导致周围液相溶质浓度升高,因此还需要对液相元胞进行溶质扩散计算,采用的控制方程可由下式表示:
步骤5、模拟计算及结果导出:
步骤5.1:基于步骤1~4所构建的模型进行编程,将编好的程序导入模拟软件Matlab中,输入相关参数,本发明中分别模拟了1060铝合金和5A06铝合金的焊接过程,合金的化学成分和热物性参数如表1至表4所示,进行计算得1060铝合金、5A06铝合金焊接熔池微观组织转变的数值模拟结果以及结论。
表1 1060铝合金的化学成分(质量百分比,wt%)
表2 1060铝合金模拟时计算所用热物性参数
表3 5A06铝合金的化学成分(质量百分比,wt%)
表4 5A06铝合金模拟时计算所用热物性参数
下面通过实施例对模拟结果进行分析。
实施例1
图3a、3b为1060铝合金焊接熔池边缘晶粒的形核与生长过程。可以看出,熔池开始凝固时,晶粒首先在熔池边缘形核,由于本发明采用的是连续型形核模型,熔池边缘在过冷度的作用下随机产生了大量晶核,如图3a所示。随着时间的推移,晶粒开始由熔池边缘向熔池中心生长,由于此时温度梯度很大,晶粒以柱状晶的形式快速向熔池中心生长,熔池中心由于温度过高无晶核产生,如图3b所示。
实施例2
图4a、4b为1060铝合金焊接熔池中心发生的CET(柱状晶-等轴晶)转变模拟结果。可以看出,随着时间的推移和熔池凝固过程的进行,熔池内部的温度梯度逐渐减小,熔池中心的温度达到了等轴晶形核和生长的要求,熔池内部发生了CET(柱状晶-等轴晶)转变,靠近熔池中心区域柱状晶和等轴晶呈现竞争生长的趋势,熔池完全凝固后组织由边缘的柱状晶和中心的等轴晶组成。
实施例3
图5a、5b为5A06铝合金不同焊接工艺参数(热输入)下熔池内微观组织的模拟结果。可以看出,随着焊接热输入的减小,焊接熔池轮廓会相应的减小,这样会导致晶粒生长的空间缩小,晶粒呈现细化的趋势,从图5a和图5b的对比可以看出,熔池中的柱状晶和等轴晶都发生了明显的细化。因此,在实际铝合金焊接过程中,应在保证焊接构件连接质量的条件下尽量减小焊接热输入以达到细化晶粒组织的目的。
实施例4
图6a、6b为5A06铝合金不同冷却工艺参数(扰动振幅)下熔池内微观组织的模拟结果。可以看出,当增大扰动振幅时,熔池内的晶粒取向发生了细微的变化,扰动振幅的增大促进了二次枝晶的生长,进而使得熔池内部的晶粒产生了细化的趋势。
综合分析上述4个实施例可以得出,本发明可以再现铝合金焊接过程中熔池内晶粒的生长过程,同时还可以研究不同焊接工艺参数和冷却工艺参数对熔池微观组织形貌的影响。
Claims (1)
1.一种铝合金焊接熔池微观组织的模拟方法,其特征在于,具体按照以下步骤实施:
步骤1、简化模型条件;
所述步骤1简化模型建立条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、将焊接熔池简化为规则的半圆形状;
简化条件3、模拟过程中忽略动力学过冷,只考虑温度过冷、成分过冷和曲率过冷;
简化条件4、模拟区域划分为正方形单元,每一个单元即为一个元胞;
简化条件5、元胞邻域关系采用V.Neumann型邻域,即四邻域;
步骤2、基于传热原理,构建宏观温度场模型;
所述步骤2具体按照以下步骤实施:
模拟时选取高斯分布热源作为焊接热源,该热源模型由下式表示:
式中:r为焊件上任一点到热源中心的距离/m;rh为电弧有效热源半径/m;U为焊接电压/V;I为焊接电流/A;η为焊接热效率;
热源加载到焊件上后,在焊件的表面和内部进行传导,热传导方程由下式表示:
式中:T为温度场函数值/K;ρ为密度/Kg·m-3;Cρ为容积比热容/J·(m3·K)-1;λ为热导率/W·(m·K)-1;
确定导热的初始条件和边界条件,初始条件由下式确定:
T(x,y)=Tf (3);
式中:Tf为模拟区域初始温度/K,Tf为室温293K;
边界条件由下式确定:
qw=hw(T-Tf) (4);
式中:T为试件表面温度/K;Tf为周围介质温度/K;hw为换热系数/W/(m2*K);
步骤3、基于插值原理,构建微观温度场模型;
所述步骤3具体按照以下步骤实施:
假定材料在各个方向的热导率相同,根据插值原理,将(2)式转化为:
式中:Δx为网格尺寸;Δt为时间步长;(i,j)为微观元胞的坐标;Ti,j为微观元胞(i,j)的当前温度值/K;T′i,j为微观元胞(i,j)的下一时刻温度值/K;
步骤4、基于凝固理论,构建晶粒形核与生长模型;
所述步骤4具体按照以下步骤实施:
晶粒形核模型:
考虑形核的过冷度和形核的连续性,采用基于高斯分布函数的准连续形核模型,该模型由下式表示:
式中:Nmax为非均匀形核密度的最大值/m3;ΔTθ为曲率过冷度/℃;ΔTmax为最大形核过冷度/℃;
晶粒生长模型:
过冷度是晶粒生长的推动力,过冷度ΔT包括三个部分,由热过冷ΔTw、成分过冷ΔTc和曲率过冷ΔTθ组成,ΔTθ由下式表示:
ΔT=ΔTw+ΔTc+ΔTθ (7);
由吉布斯汤普森关系可得t时刻的过冷度为:
ΔT(t)=TL-T(t)+mL(CL(t)-C0)-ΓK(t) (8);
式中:TL为液相线温度/K;T(t)为t时刻的瞬时温度/K;mL为液相线斜率;C0为初始溶质浓度/wt.%;CL(t)为t时刻的液相溶质浓度/wt.%;Γ为Gibbs-Thompson系数;K(t)为t时刻的界面平均曲率;
界面元胞在凝固过程中会向周围的液相元胞排出多余溶质,单位时间Δt内排出的溶质量ΔC由下式计算:
式中:DL为溶质液相扩散系数/m2/s;CL为固液界面处液相溶质浓度/wt.%;nb表示界面元胞周围的液相元胞,Cnb为其溶质浓度/wt.%;
固/液界面推进过程中,界面元胞排出ΔC的溶质之后,元胞单元固相率会发生变化,固相率增量由下式表示:
式中:k0为溶质分配系数;A为扰动因子;rand为0到1之间的随机数;
采用如下公式(11)对液相元胞进行溶质扩散计算:
步骤5、模拟计算及结果导出;
所述步骤5体按照以下步骤实施:
步骤5.1:基于步骤1~4所构建的一种铝合金焊接熔池微观组织的数值模型进行编程;
步骤5.2:将编好的程序导入模拟软件Matlab中,输入相关铝合金的热物性参数和焊接工艺参数,即得到铝合金焊接熔池微观组织的数值模拟结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110571832.9A CN113284570B (zh) | 2021-05-25 | 2021-05-25 | 一种铝合金焊接熔池微观组织的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110571832.9A CN113284570B (zh) | 2021-05-25 | 2021-05-25 | 一种铝合金焊接熔池微观组织的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113284570A CN113284570A (zh) | 2021-08-20 |
CN113284570B true CN113284570B (zh) | 2024-03-19 |
Family
ID=77281384
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110571832.9A Active CN113284570B (zh) | 2021-05-25 | 2021-05-25 | 一种铝合金焊接熔池微观组织的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113284570B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114864012B (zh) * | 2022-04-08 | 2024-05-24 | 浙江海洋大学 | 一种预测合金焊接熔池沿凝固路径的晶粒尺寸的方法 |
CN114864017B (zh) * | 2022-04-11 | 2024-04-12 | 华南理工大学 | 一种铝合金铸件全域晶粒尺寸的预测计算方法和系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489818A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
CN110489821A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种镍基合金堆焊熔池枝晶生长数值模拟方法 |
CN110489820A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法 |
CN110619157A (zh) * | 2019-08-30 | 2019-12-27 | 西安理工大学 | 一种钛合金焊接熔池微观组织演化的模拟方法 |
CN111261241A (zh) * | 2020-02-28 | 2020-06-09 | 西安理工大学 | 一种镍基合金多道堆焊过程中枝晶生长数值模拟方法 |
CN111259513A (zh) * | 2019-09-25 | 2020-06-09 | 昆明理工大学 | 基于ca理论对铝合金副车架性能优化的方法 |
CN112185474A (zh) * | 2020-09-07 | 2021-01-05 | 西安理工大学 | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 |
CN112784424A (zh) * | 2021-01-28 | 2021-05-11 | 西安理工大学 | 一种钛合金焊接熔池晶粒生长三维数值模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8655476B2 (en) * | 2011-03-09 | 2014-02-18 | GM Global Technology Operations LLC | Systems and methods for computationally developing manufacturable and durable cast components |
TWI589373B (zh) * | 2015-11-20 | 2017-07-01 | Metal Ind Res & Dev Ct | 用於連續鑄造製程之金屬凝固微觀組織模擬預測方法 |
-
2021
- 2021-05-25 CN CN202110571832.9A patent/CN113284570B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489818A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
CN110489821A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种镍基合金堆焊熔池枝晶生长数值模拟方法 |
CN110489820A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法 |
CN110619157A (zh) * | 2019-08-30 | 2019-12-27 | 西安理工大学 | 一种钛合金焊接熔池微观组织演化的模拟方法 |
CN111259513A (zh) * | 2019-09-25 | 2020-06-09 | 昆明理工大学 | 基于ca理论对铝合金副车架性能优化的方法 |
CN111261241A (zh) * | 2020-02-28 | 2020-06-09 | 西安理工大学 | 一种镍基合金多道堆焊过程中枝晶生长数值模拟方法 |
CN112185474A (zh) * | 2020-09-07 | 2021-01-05 | 西安理工大学 | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 |
CN112784424A (zh) * | 2021-01-28 | 2021-05-11 | 西安理工大学 | 一种钛合金焊接熔池晶粒生长三维数值模拟方法 |
Non-Patent Citations (8)
Title |
---|
TC4合金焊接熔池微观组织演变的数值模拟与分析;张敏;黄超;郭宇飞;张立胜;李继红;;中国有色金属学报(第08期);全文 * |
Ti-45Al合金焊接熔池凝固过程数值模拟;张敏;周玉兰;薛覃;李继红;;焊接学报(第03期);全文 * |
ULCB钢焊接熔池凝固过程数值模拟;张敏;惠媛媛;周玉兰;李继红;;电焊机(第07期);全文 * |
元胞自动机法模拟铝合金三维枝晶生长;许林,郭洪民,杨湘杰;铸造(第06期);全文 * |
基于CAFE模拟双辊连铸纯铝凝固微观组织;陈守东;陈敬超;吕连灏;;材料科学与工艺(第03期);全文 * |
基于CA-FE的双辊连铸纯铝凝固组织模拟;陈守东;陈敬超;吕连灏;;材料工程;20121020(第10期);全文 * |
焊接熔池凝固过程组织演变模拟;李玉斌;蒙大桥;刘柯钊;谢志强;;焊接学报(第04期);全文 * |
镍基合金TIG焊接熔池及热影响区组织模拟;刘仁培;陈莉莉;魏艳红;;焊接学报(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113284570A (zh) | 2021-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112784424B (zh) | 一种钛合金焊接熔池晶粒生长三维数值模拟方法 | |
CN113284570B (zh) | 一种铝合金焊接熔池微观组织的模拟方法 | |
CN110489820B (zh) | 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法 | |
Roehling et al. | Modulating laser intensity profile ellipticity for microstructural control during metal additive manufacturing | |
CN110489821B (zh) | 一种镍基合金堆焊熔池枝晶生长数值模拟方法 | |
CN110489818B (zh) | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 | |
Xie et al. | Investigation on the Laves phase formation during laser cladding of IN718 alloy by CA-FE | |
CN110619157B (zh) | 一种钛合金焊接熔池微观组织演化的模拟方法 | |
CN110321604B (zh) | 一种三元合金凝固过程单个枝晶生长数值模拟方法 | |
Lingda et al. | A phase-field simulation of columnar-to-equiaxed transition in the entire laser welding molten pool | |
Wang et al. | Effects of scanning speed on microstructure in laser surface-melted single crystal superalloy and theoretical analysis | |
CN105868468B (zh) | 一种基于元胞自动机的新型邻居捕获方法 | |
CN111261241B (zh) | 一种镍基合金多道堆焊过程中枝晶生长数值模拟方法 | |
CN112185474B (zh) | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 | |
Fonda et al. | Crystallographic variability in additive manufacturing | |
Yao et al. | Effect of steel strip feeding on the columnar-equiaxed solidification in a large continuous casting round bloom | |
CN111640472A (zh) | 镍基单晶高温合金的固溶热处理模拟方法 | |
CN111339670B (zh) | 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法 | |
CN116186920A (zh) | 一种改善车用6系铝合金激光焊对接接头下塌及接头软化的工艺方法 | |
CN116230142B (zh) | 一种铝合金凝固动力学过程的介观尺度预测方法 | |
Zhang et al. | Research on 316 stainless steel low-power pulsed laser-induced arc additive manufacturing by different deposition routes | |
CN113130019B (zh) | 一种tc4钛合金焊接熔池固态相变的模拟方法 | |
CN105740537A (zh) | 一种焊接过程气孔形成与演变的二维计算机模拟方法 | |
CN113823357B (zh) | 一种四元合金凝固过程等轴晶生长数值模拟方法 | |
CN116702505A (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 |