CN112784424A - 一种钛合金焊接熔池晶粒生长三维数值模拟方法 - Google Patents
一种钛合金焊接熔池晶粒生长三维数值模拟方法 Download PDFInfo
- Publication number
- CN112784424A CN112784424A CN202110119796.2A CN202110119796A CN112784424A CN 112784424 A CN112784424 A CN 112784424A CN 202110119796 A CN202110119796 A CN 202110119796A CN 112784424 A CN112784424 A CN 112784424A
- Authority
- CN
- China
- Prior art keywords
- formula
- titanium alloy
- welding
- solute
- grain 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.)
- Granted
Links
- 238000003466 welding Methods 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000004088 simulation Methods 0.000 title claims abstract description 54
- 229910001069 Ti alloy Inorganic materials 0.000 title claims abstract description 50
- 230000008569 process Effects 0.000 claims abstract description 41
- 238000010899 nucleation Methods 0.000 claims abstract description 27
- 230000006911 nucleation Effects 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 230000001052 transient effect Effects 0.000 claims abstract description 8
- 238000009795 derivation Methods 0.000 claims abstract description 5
- 238000012546 transfer Methods 0.000 claims abstract description 5
- 238000004781 supercooling Methods 0.000 claims description 37
- 239000007791 liquid phase Substances 0.000 claims description 24
- 238000009792 diffusion process Methods 0.000 claims description 9
- 239000007788 liquid Substances 0.000 claims description 9
- 239000007787 solid Substances 0.000 claims description 9
- 238000007711 solidification Methods 0.000 claims description 8
- 230000008023 solidification Effects 0.000 claims description 8
- 230000001413 cellular effect Effects 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 7
- 239000007790 solid phase Substances 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 7
- 210000001787 dendrite Anatomy 0.000 abstract 1
- 230000001737 promoting effect Effects 0.000 abstract 1
- 239000013078 crystal Substances 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 239000007769 metal material Substances 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010891 electric arc Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000004021 metal welding Methods 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- 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
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- 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
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)
- Welding Or Cutting Using Electron Beams (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种钛合金焊接熔池晶粒生长三维数值模拟方法,具体步骤如下:步骤1、简化模型条件;步骤2、基于传热原理,构建焊接瞬态宏观温度场模型;步骤3、基于插值原理,构建微观温度场模型;步骤4、构建晶粒形核与生长模型;步骤5、模拟计算及结果导出。本发明建立的模型能够在三维尺度上模拟钛合金焊接熔池晶粒生长过程,以及研究焊接工艺参数对熔池大小/形状以及熔池内枝晶形貌的影响,从而对研究钛合金焊接过程中熔池微观组织演化过程起到了一定的促进作用。
Description
技术领域
本发明属于金属材料焊接工艺数值模拟技术领域,具体涉及一种钛合金焊接熔池晶粒生长三维数值模拟方法。
背景技术
随着航空航天、船舶化工、生物医疗等领域的快速发展,对于材料性能的要求也越来越严苛,在众多金属材料中,钛合金具有密度低、比强度高、耐高温性、耐腐蚀性好等优点。钛合金的连接工艺中,焊接具有连接性能好、结构刚度大、适应性强等优点,是钛合金最主要的连接工艺之一。焊缝的微观组织对其力学性能具有重要影响,焊接具有瞬时、高温、高压的特征,采用传统的实验方法难以对焊接过程中任意时刻的焊接熔池状况进行研究,即使用特种实验设备,也不可避免的存在成本过高的问题。因此,急需一种成本可控且物理背景明确的研究方法,实现对焊接过程的可视化分析。近年来,计算机技术飞速发展,采用数值模拟方法对金属材料焊接过程进行研究为相关学者提供了新的研究思路。目前该方法已经在金属材料的焊接领域取得了一些成果,其基于物理基础构建数学模型,通过计算机模拟软件实现了模拟结果的可视化。
目前,对于焊接领域的数值模拟计算多集中于二维模型的建立与分析,对于三维模型的建立与分析鲜有报道,然而二维模型分析问题时,只能对模型某个截面的特征进行研究,三维模型较二维模型来说,其模拟结果更加精确和直观,因此,建立一种钛合金焊接熔池晶粒生长三维数值模拟方法显得尤为重要。
发明内容
本发明的目的是提供一种钛合金焊接熔池晶粒生长三维数值模拟方法,解决了当前有关钛合金焊接过程缺少的晶粒生长三维数值模拟方法的问题。
本发明所采用的技术方案是,一种钛合金焊接熔池晶粒生长三维数值模拟方法,其特征在于,具体按照以下步骤实施:
步骤1、简化模型条件;
步骤2、基于传热原理,构建焊接瞬态宏观温度场模型;
步骤3、基于插值原理,构建微观温度场模型;
步骤4、构建晶粒形核与生长模型;
步骤5、模拟计算及结果导出。
本发明的特点还在于:
步骤1简化模型建立条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、模拟过程中为减少计算量,将物体的表面换热系数视为定值;
简化条件3、忽略动力学过冷,只考虑温度过冷、成分过冷和曲率过冷;
简化条件4、模拟区域划分为正方体单元,每一个单元即为一个元胞;
简化条件5、元胞邻域关系采用Moore型邻域,即八邻域。
步骤2按照以下具体步骤实施:
步骤2.1、模拟时选取高斯分布热源作为焊接热源,所述热源模型可由下式表示:
式(1)中:r为焊件上任一点到热源中心的距离/m;rh为电弧有效热源半径/m;U为焊接电压/V;I为焊接电流/A;η为焊接热效率;
步骤2.2、热源加载到焊件上后,在所述焊件的表面和内部进行传导,焊接三维瞬态热传导方程可由下式表示:
式(2)中:T为温度场函数值/℃;ρ为密度/Kg·m-3;c为容积比热容/J·(m3·℃)-1;λ为热导率/W·(m·℃)-1;qv为热源强度/W·m-3;
步骤2.3、确定导热的初始条件和边界条件;
初始条件可由下式确定:
T0=Tf (3)
式(3)中:T0为模拟区域初始温度/℃;Tf为室温20℃;
边界条件是指物体表面与周围介质的热交换情况,边界条件可由下式确定:
qw=hw(T-Tf) (4)
式(4)中:qw为体表面与周围介质的热交换量/J;hw为表面换热系数/W/(m2·℃)。
步骤3具体按照以下步骤实施:
步骤3.1、假定材料在各个方向的热导率相同,根据插值原理,将(2)式转化为:
步骤3.2、计算得到宏观单元的温度数据后,基于式(5),通过插值原理求解微观单元的温度值,实现温度场宏微观耦合过程。
步骤4具体按照以下步骤实施:
晶粒形核模型:
考虑形核的过冷度和形核的连续性,采用基于高斯分布函数的准连续形核模型,所述准连续形核模型由下式表示:
式(6)中: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)
式(8)中:TL为液相线温度/℃;T(t)为t时刻的瞬时温度/℃;mL为液相线斜率;C0为初始溶质浓度/wt.%;CL(t)为t时刻的液相溶质浓度/wt.%;Γ为Gibbs-Thompson系数;K(t)为t时刻的界面平均曲率;
界面元胞在凝固过程中会向周围的液相元胞排出多余溶质,单位时间ΔT内排出的溶质量ΔC可由下式计算:
式(9)中:DL为溶质液相扩散系数/m2/s;CL为固液界面处液相溶质浓度/wt.%;nb表示界面元胞周围的液相元胞,Cnb为其溶质浓度/wt.%。
固/液界面推进过程中,界面元胞排出ΔC的溶质之后,元胞单元固相率会发生变化,固相率增量可由下式表示:
式(10)中:k0为溶质分配系数;A为扰动因子;rand为0到1之间的随机数。
界面生长过程中,多余的溶质排到周围液相邻胞中会导致周围液相溶质浓度升高,需要对液相元胞进行溶质扩散计算,采用的控制方程可由下式表示:
式(11)中:DS为溶质固相扩散系数/m2/s。
步骤5体按照以下步骤实施:
步骤5.1:基于步骤1~4所构建的一种钛合金焊接熔池晶粒生长三维数值模拟模型进行编程;
步骤5.2:将编好的程序导入模拟软件Matlab中,输入钛合金的热物性参数和焊接工艺参数,即可得到钛合金焊接熔池晶粒生长三维数值模拟结果。
本发明的有益效果是:
(1)本发明提供一种钛合金焊接熔池晶粒生长三维数值模拟方法,解决了现有技术中存在的缺少关于钛合金焊接过程的三维数值模型问题;
(2)本发明可以在三维尺度上再现钛合金焊接过程中熔池晶粒生长情况,模拟结果更加精确和直观,从而为进一步研究钛合金焊接过程提供了新的研究方案;
(3)相对于传统的通过金相、扫描等实验手段研究钛合金的焊接过程及其焊缝微观组织的转变过程,本发明基于物理学基础构建了数学模型,通过编程利用计算机模拟平台对钛合金的焊接过程进行了三维数值模拟,节省了大量的人力物力投入,经济高效且节能环保。
附图说明
图1是本发明一种钛合金焊接熔池晶粒生长三维数值模拟方法的流程图;
图2是本发明一种钛合金焊接熔池晶粒生长三维数值模拟方法三维元胞邻域关系示意图;
图3是实施例1 TA15钛合金焊接速度为4mm/s时的焊接熔池晶粒生长三维模拟结果;
图4是实施例2 TA15钛合金焊接速度为5mm/s时的焊接熔池晶粒生长三维模拟结果;
图5是实施例3 TC4钛合金焊接电流为60A时的焊接熔池晶粒生长三维模拟结果;
图6是实施例3 TC4钛合金焊接电流为70A时的焊接熔池晶粒生长三维模拟结果。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种钛合金焊接熔池晶粒生长三维数值模拟方法,如图1所示,具体按照以下步骤实施:
步骤1、简化模型条件;
步骤2、基于传热原理,构建焊接瞬态宏观温度场模型;
步骤3、基于插值原理,构建微观温度场模型;
步骤4、构建晶粒形核与生长模型;
步骤5、模拟计算及结果导出。
步骤1、简化模型建立条件:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、模拟过程中为减少计算量,将物体的表面换热系数视为定值;
简化条件3、忽略动力学过冷,只考虑温度过冷、成分过冷和曲率过冷;
简化条件4、模拟区域划分为正方体单元,每一个单元即为一个元胞;
简化条件5、元胞邻域关系采用Moore型邻域,即八邻域。
步骤2、基于传热原理,构建焊接瞬态宏观温度场模型:
施加到焊件上的热量的大小、分布及传导影响着熔池金属的微观组织转变过程,焊接过程中,电弧直接作用于焊件表面形成加热斑点,选取高斯分布热源作为焊接热源,该热源模型可由下式表示:
式(1)中:r为焊件上任一点到热源中心的距离/m;rh为电弧有效热源半径/m;U为焊接电压/V;I为焊接电流/A;η为焊接热效率。
热源加载到焊件上后,会在焊件的表面和内部进行传导,焊接三维瞬态热传导方程可由下式表示:
式(2)中:T为温度场函数值/℃;ρ为密度/Kg·m-3;c为容积比热容/J·(m3·℃)-1;λ为热导率/W·(m·℃)-1;qv为热源强度/W·m-3。
求解导热问题之前,需要确定导热的初始条件和边界条件,本发明的模拟环境为室温,故初始条件可由下式确定:
T0=Tf (3)
式(3)中:T0为模拟区域初始温度/℃;Tf为室温20℃。
边界条件是指物体表面与周围介质的热交换情况,本发明所选边界条件规定物体表面自由换热,边界条件可由下式确定:
qw=hw(T-Tf) (4)
式(4)中:qw为体表面与周围介质的热交换量/J;hw为表面换热系数/W/(m2·℃)。
步骤3、基于插值原理,构建微观温度场模型:
温度场模型属于宏观范畴,而焊接凝固过程属于微观范畴,因此,需要进行温度场的宏微观耦合,由于所建模型中设定微观元胞为规则正方体,根据插值原理,可将(2)式转化为:
计算得到宏观单元的温度数据后,基于上式,通过插值原理求解微观单元的温度值,便可得到模拟区域内部元胞、表面元胞以及边角元胞的温度数据,从而实现温度场宏微观耦合过程。
步骤4、构建晶粒形核与生长模型:
晶粒形核模型:
考虑形核的过冷度和形核的连续性,采用基于高斯分布函数的准连续形核模型,该模型可由下式表示:
式(6)中: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)
式(8)中:TL为液相线温度/℃;T(t)为t时刻的瞬时温度/℃;mL为液相线斜率;C0为初始溶质浓度/wt.%;CL(t)为t时刻的液相溶质浓度/wt.%;Γ为Gibbs-Thompson系数;K(t)为t时刻的界面平均曲率。
界面元胞在凝固过程中会向周围的液相元胞排出多余溶质,单位时间ΔT内排出的溶质量ΔC可由下式计算:
式(9)中:DL为溶质液相扩散系数/m2/s;CL为固液界面处液相溶质浓度/wt.%;nb表示界面元胞周围的液相元胞,Cnb为其溶质浓度/wt.%。
固/液界面推进过程中,界面元胞排出ΔC的溶质之后,元胞单元固相率会发生变化,固相率增量可由下式表示:
式(10)中:k0为溶质分配系数;A为扰动因子;rand为0到1之间的随机数。
界面生长过程中,多余的溶质排到周围液相邻胞中会导致周围液相溶质浓度升高,因此还需要对液相元胞进行溶质扩散计算,采用的控制方程可由下式表示:
式(11)中:DS为溶质固相扩散系数/m2/s。
步骤5、模拟计算及结果导出:
步骤5.1:基于步骤1~4所构建的模型进行编程,将编好的程序导入模拟软件Matlab中,输入TC4钛合金和TA15钛合金的热物性参数,如表1和表2所示,以及焊接工艺参数,进行计算得到不同类型钛合金焊接熔池晶粒生长三维数值模拟结果以及结论。
表1 TC4钛合金模拟时计算所用热物性参数
表2 TA15钛合金模拟时计算所用热物性参数
实施例1
图3a、3b为焊接电流为60A时TC4钛合金焊接过程中熔池内晶粒生长的三维模拟结果。可以看出,熔池凝固过程中,晶粒在熔池的边缘形核,并以柱状晶的形式向熔池中心生长,由于熔池内温度过高,没有达到晶粒形核的条件,故熔池中心无晶核产生。
实施例2
图4a、4b为焊接电流为70A时TC4钛合金焊接过程中熔池内晶粒生长的三维模拟结果。可以看出,由于热输入较实施例3大,焊件被完全熔透,焊缝熔池截面形状呈“梯形”,且熔池上熔宽大于下熔宽,这是由于热源直接作用于焊件上表面导致的。
实施例3
图5a、5b为焊接速度为4mm/s时TA15钛合金焊接过程中熔池内晶粒生长的三维模拟结果。可以看出,随着焊接热源的移动,先前熔化的部分逐渐凝固,且熔池内依旧由粗大的柱状晶组成,焊件后半部分的熔池略大于前半部分,这是由于焊接过程的热量积累造成的。
实施例4
图6a、6b为焊接速度为5mm/s时TA15钛合金焊接过程中熔池内晶粒生长的三维模拟结果。可以看出,当焊接速度较大时,由于TA15钛合金的熔点较高,焊件前端由于热输入不足未能形成熔池,由于实施例4的焊接速度大于实施例3,导致焊接线能量减小,因此实施例4的熔池较实施例3小,焊缝完全凝固后,组织为粗大的柱状晶。
综合分析上述4个实施例可以得出,本发明可以完整再现钛合金焊接过程中熔池内晶粒生长的三维变化过程,同时还可以研究不同焊接工艺参数对熔池大小和形貌的影响,且模拟结果更加直观。
Claims (6)
1.一种钛合金焊接熔池晶粒生长三维数值模拟方法,其特征在于,具体按照以下步骤实施:
步骤1、简化模型条件;
步骤2、基于传热原理,构建焊接瞬态宏观温度场模型;
步骤3、基于插值原理,构建微观温度场模型;
步骤4、构建晶粒形核与生长模型;
步骤5、模拟计算及结果导出。
2.根据权利要求1所述的一种钛合金焊接熔池晶粒生长三维数值模拟方法,其特征在于,所述步骤1简化模型建立条件包括:
简化条件1、整个凝固过程只存在液相、固相和界面三种元胞状态;
简化条件2、模拟过程中为减少计算量,将物体的表面换热系数视为定值;
简化条件3、忽略动力学过冷,只考虑温度过冷、成分过冷和曲率过冷;
简化条件4、模拟区域划分为正方体单元,每一个单元即为一个元胞;
简化条件5、元胞邻域关系采用Moore型邻域,即八邻域。
3.根据权利要求2所述的一种钛合金焊接熔池晶粒生长三维数值模拟方法,其特征在于,所述步骤2按照以下具体步骤实施:
步骤2.1、模拟时选取高斯分布热源作为焊接热源,所述热源模型可由下式表示:
式(1)中:r为焊件上任一点到热源中心的距离/m;rh为电弧有效热源半径/m;U为焊接电压/V;I为焊接电流/A;η为焊接热效率;
步骤2.2、热源加载到焊件上后,在所述焊件的表面和内部进行传导,焊接三维瞬态热传导方程可由下式表示:
式(2)中:T为温度场函数值/℃;ρ为密度/Kg·m-3;c为容积比热容/J·(m3·℃)-1;λ为热导率/W·(m·℃)-1;qv为热源强度/W·m-3;
步骤2.3、确定导热的初始条件和边界条件;
初始条件可由下式确定:
T0=Tf (3)
式(3)中:T0为模拟区域初始温度/℃;Tf为室温20℃;
所述边界条件是指物体表面与周围介质的热交换情况,边界条件由下式确定:
qw=hw(T-Tf) (4)
式(4)中:qw为体表面与周围介质的热交换量/J;hw为表面换热系数/W/(m2·℃)。
5.根据权利要求1所述的一种钛合金焊接熔池晶粒生长三维数值模拟方法,其特征在于,所述步骤4具体按照以下步骤实施:
晶粒形核模型:
考虑形核的过冷度和形核的连续性,采用基于高斯分布函数的准连续形核模型,所述准连续形核模型由下式表示:
式(6)中: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)
式(8)中:TL为液相线温度/℃;T(t)为t时刻的瞬时温度/℃;mL为液相线斜率;C0为初始溶质浓度/wt.%;CL(t)为t时刻的液相溶质浓度/wt.%;Γ为Gibbs-Thompson系数;K(t)为t时刻的界面平均曲率;
界面元胞在凝固过程中会向周围的液相元胞排出多余溶质,单位时间ΔT内排出的溶质量ΔC可由下式计算:
式(9)中:DL为溶质液相扩散系数/m2/s;CL为固液界面处液相溶质浓度/wt.%;nb表示界面元胞周围的液相元胞,Cnb为其溶质浓度/wt.%。
固/液界面推进过程中,界面元胞排出ΔC的溶质之后,元胞单元固相率会发生变化,固相率增量可由下式表示:
式(10)中:k0为溶质分配系数;A为扰动因子;rand为0到1之间的随机数。
界面生长过程中,多余的溶质排到周围液相邻胞中会导致周围液相溶质浓度升高,需要对液相元胞进行溶质扩散计算,采用的控制方程可由下式表示:
式(11)中:DS为溶质固相扩散系数/m2/s。
6.根据权利要求1所述的一种钛合金焊接熔池晶粒生长三维数值模拟方法,其特征在于,所述步骤5体按照以下步骤实施:
步骤5.1:基于步骤1~4所构建的一种钛合金焊接熔池晶粒生长三维数值模拟模型进行编程;
步骤5.2:将编好的程序导入模拟软件Matlab中,输入钛合金的热物性参数和焊接工艺参数,即可得到钛合金焊接熔池晶粒生长三维数值模拟结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110119796.2A CN112784424B (zh) | 2021-01-28 | 2021-01-28 | 一种钛合金焊接熔池晶粒生长三维数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110119796.2A CN112784424B (zh) | 2021-01-28 | 2021-01-28 | 一种钛合金焊接熔池晶粒生长三维数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112784424A true CN112784424A (zh) | 2021-05-11 |
CN112784424B CN112784424B (zh) | 2024-04-09 |
Family
ID=75759441
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110119796.2A Active CN112784424B (zh) | 2021-01-28 | 2021-01-28 | 一种钛合金焊接熔池晶粒生长三维数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112784424B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113284570A (zh) * | 2021-05-25 | 2021-08-20 | 西安理工大学 | 一种铝合金焊接熔池微观组织的模拟方法 |
CN113770383A (zh) * | 2021-09-16 | 2021-12-10 | 南京智能高端装备产业研究院有限公司 | 一种基于晶粒形貌预测确定增材制造成形工艺参数的方法 |
CN114864017A (zh) * | 2022-04-11 | 2022-08-05 | 华南理工大学 | 一种铝合金铸件全域晶粒尺寸的预测计算方法和系统 |
CN114864012A (zh) * | 2022-04-08 | 2022-08-05 | 浙江海洋大学 | 一种预测合金焊接熔池沿凝固路径的晶粒尺寸的方法 |
CN114943159A (zh) * | 2022-07-11 | 2022-08-26 | 北京科技大学 | 一种基于分子动力学的流场作用下金属凝固的模拟方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004004034A (ja) * | 2002-04-26 | 2004-01-08 | Kobe Steel Ltd | 溶接部の材質予測方法 |
CN105320804A (zh) * | 2014-08-01 | 2016-02-10 | 通用汽车环球科技运作有限责任公司 | 用于铸造铝合金的材料性质预测器 |
CN107309543A (zh) * | 2017-05-04 | 2017-11-03 | 南京航空航天大学 | 一种激光焊接熔池枝晶生长模拟方法 |
CN110321604A (zh) * | 2019-06-18 | 2019-10-11 | 西安理工大学 | 一种三元合金凝固过程单个枝晶生长数值模拟方法 |
CN110489820A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法 |
CN110619157A (zh) * | 2019-08-30 | 2019-12-27 | 西安理工大学 | 一种钛合金焊接熔池微观组织演化的模拟方法 |
CN111261241A (zh) * | 2020-02-28 | 2020-06-09 | 西安理工大学 | 一种镍基合金多道堆焊过程中枝晶生长数值模拟方法 |
-
2021
- 2021-01-28 CN CN202110119796.2A patent/CN112784424B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004004034A (ja) * | 2002-04-26 | 2004-01-08 | Kobe Steel Ltd | 溶接部の材質予測方法 |
CN105320804A (zh) * | 2014-08-01 | 2016-02-10 | 通用汽车环球科技运作有限责任公司 | 用于铸造铝合金的材料性质预测器 |
CN107309543A (zh) * | 2017-05-04 | 2017-11-03 | 南京航空航天大学 | 一种激光焊接熔池枝晶生长模拟方法 |
CN110321604A (zh) * | 2019-06-18 | 2019-10-11 | 西安理工大学 | 一种三元合金凝固过程单个枝晶生长数值模拟方法 |
CN110489820A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法 |
CN110619157A (zh) * | 2019-08-30 | 2019-12-27 | 西安理工大学 | 一种钛合金焊接熔池微观组织演化的模拟方法 |
CN111261241A (zh) * | 2020-02-28 | 2020-06-09 | 西安理工大学 | 一种镍基合金多道堆焊过程中枝晶生长数值模拟方法 |
Non-Patent Citations (1)
Title |
---|
欧阳智勇;卢国杰;郭亮;张庆茂;: "激光焊接温度场模拟分析及熔池深度的计算", 应用激光, no. 01, 15 February 2018 (2018-02-15) * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113284570A (zh) * | 2021-05-25 | 2021-08-20 | 西安理工大学 | 一种铝合金焊接熔池微观组织的模拟方法 |
CN113284570B (zh) * | 2021-05-25 | 2024-03-19 | 西安理工大学 | 一种铝合金焊接熔池微观组织的模拟方法 |
CN113770383A (zh) * | 2021-09-16 | 2021-12-10 | 南京智能高端装备产业研究院有限公司 | 一种基于晶粒形貌预测确定增材制造成形工艺参数的方法 |
CN114864012A (zh) * | 2022-04-08 | 2022-08-05 | 浙江海洋大学 | 一种预测合金焊接熔池沿凝固路径的晶粒尺寸的方法 |
CN114864012B (zh) * | 2022-04-08 | 2024-05-24 | 浙江海洋大学 | 一种预测合金焊接熔池沿凝固路径的晶粒尺寸的方法 |
CN114864017A (zh) * | 2022-04-11 | 2022-08-05 | 华南理工大学 | 一种铝合金铸件全域晶粒尺寸的预测计算方法和系统 |
CN114864017B (zh) * | 2022-04-11 | 2024-04-12 | 华南理工大学 | 一种铝合金铸件全域晶粒尺寸的预测计算方法和系统 |
CN114943159A (zh) * | 2022-07-11 | 2022-08-26 | 北京科技大学 | 一种基于分子动力学的流场作用下金属凝固的模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112784424B (zh) | 2024-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112784424A (zh) | 一种钛合金焊接熔池晶粒生长三维数值模拟方法 | |
CN110619157B (zh) | 一种钛合金焊接熔池微观组织演化的模拟方法 | |
Li et al. | Reducing arc heat input and obtaining equiaxed grains by hot-wire method during arc additive manufacturing titanium alloy | |
CN110489821B (zh) | 一种镍基合金堆焊熔池枝晶生长数值模拟方法 | |
CN110489820B (zh) | 一种基于元胞自动机法的焊接熔池微观组织演化模拟方法 | |
CN110489818B (zh) | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 | |
CN113284570B (zh) | 一种铝合金焊接熔池微观组织的模拟方法 | |
CN108856974B (zh) | 一种超声场耦合共熔池双丝cmt电弧增材成型工艺方法 | |
Pavlyk et al. | Simulation of weld solidification microstructure and its coupling to the macroscopic heat and fluid flow modelling | |
CN108127117A (zh) | 一种针对飞行器异形舵轴整体快速成形的工艺方法 | |
Liu et al. | Study on 2219 aluminum alloy T-joint during dual laser-beam bilateral synchronous welding: Effect of the welding speed and incident beam angle on grain morphology | |
Liu et al. | External filler wire based GMA-AM process of 2219 aluminum alloy | |
CN105057642A (zh) | 铸件晶粒组织形成相关数值的模拟方法 | |
Su et al. | Mechanical property and characterization of TA1 titanium alloy sheets welded by vacuum electron beam welding | |
CN112185474B (zh) | 一种Ti-45%Al合金定向凝固过程的数值模拟方法 | |
Wu et al. | Double pulse low-frequency modulation for high-power double-wire pulsed GMAW | |
Fonda et al. | Crystallographic variability in additive manufacturing | |
CN111261241A (zh) | 一种镍基合金多道堆焊过程中枝晶生长数值模拟方法 | |
Ge et al. | Achievement of forging without canning for β-solidifying γ-TiAl alloy containing high content of niobium | |
Liu et al. | Microstructure evolution of large-scale titanium slab ingot based on CAFE method during EBCHM | |
CN111339670B (zh) | 一种镍基合金激光熔覆过程中枝晶生长数值模拟方法 | |
Mao et al. | Improve the manufacturing efficiency of steel bars by using hot-wire pulse arc additive manufacturing | |
CN113192565A (zh) | 一种钛铝合金定向凝固过程晶粒生长的三维数值模拟方法 | |
CN113130019A (zh) | 一种tc4钛合金焊接熔池固态相变的模拟方法 | |
Gao et al. | Cross-scale simulation for MnS precipitation of Fe-C alloy with cooling rate variation |
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 |