一种钛合金焊接熔池晶粒生长三维数值模拟方法技术

本发明专利技术公开了一种钛合金焊接熔池晶粒生长三维数值模拟方法，具体步骤如下：步骤1、简化模型条件；步骤2、基于传热原理，构建焊接瞬态宏观温度场模型；步骤3、基于插值原理，构建微观温度场模型；步骤4、构建晶粒形核与生长模型；步骤5、模拟计算及结果导出。本发明专利技术建立的模型能够在三维尺度上模拟钛合金焊接熔池晶粒生长过程，以及研究焊接工艺参数对熔池大小/形状以及熔池内枝晶形貌的影响，从而对研究钛合金焊接过程中熔池微观组织演化过程起到了一定的促进作用。

【技术实现步骤摘要】
一种钛合金焊接熔池晶粒生长三维数值模拟方法
本专利技术属于金属材料焊接工艺数值模拟
，具体涉及一种钛合金焊接熔池晶粒生长三维数值模拟方法。
技术介绍
随着航空航天、船舶化工、生物医疗等领域的快速发展，对于材料性能的要求也越来越严苛，在众多金属材料中，钛合金具有密度低、比强度高、耐高温性、耐腐蚀性好等优点。钛合金的连接工艺中，焊接具有连接性能好、结构刚度大、适应性强等优点，是钛合金最主要的连接工艺之一。焊缝的微观组织对其力学性能具有重要影响，焊接具有瞬时、高温、高压的特征，采用传统的实验方法难以对焊接过程中任意时刻的焊接熔池状况进行研究，即使用特种实验设备，也不可避免的存在成本过高的问题。因此，急需一种成本可控且物理背景明确的研究方法，实现对焊接过程的可视化分析。近年来，计算机技术飞速发展，采用数值模拟方法对金属材料焊接过程进行研究为相关学者提供了新的研究思路。目前该方法已经在金属材料的焊接领域取得了一些成果，其基于物理基础构建数学模型，通过计算机模拟软件实现了模拟结果的可视化。目前，对于焊接领域的数值模拟计算多集中于二维模型的建立与分析，对于三维模型的建立与分析鲜有报道，然而二维模型分析问题时，只能对模型某个截面的特征进行研究，三维模型较二维模型来说，其模拟结果更加精确和直观，因此，建立一种钛合金焊接熔池晶粒生长三维数值模拟方法显得尤为重要。
技术实现思路
本专利技术的目的是提供一种钛合金焊接熔池晶粒生长三维数值模拟方法，解决了当前有关钛合金焊接过程缺少的晶粒生长三维数值模拟方法的问题。本专利技术所采用的技术方案是，一种钛合金焊接熔池晶粒生长三维数值模拟方法，其特征在于，具体按照以下步骤实施：步骤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)式转化为：式(5)中：(i,j,z)为微观元胞的三维坐标；为微观元胞p(i,j,z)的当前温度值/℃；为微观元胞p(i,j,z)的下一时刻温度值/℃；步骤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是实施例1TA15钛合金焊接本文档来自技高网...

【技术保护点】
1.一种钛合金焊接熔池晶粒生长三维数值模拟方法，其特征在于，具体按照以下步骤实施：/n步骤1、简化模型条件；/n步骤2、基于传热原理，构建焊接瞬态宏观温度场模型；/n步骤3、基于插值原理，构建微观温度场模型；/n步骤4、构建晶粒形核与生长模型；/n步骤5、模拟计算及结果导出。/n

【技术特征摘要】
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·℃)。


4.根据权利要求3所述的一种钛合金焊接熔池晶粒生长三维数值模拟方法，其特征在于，所述步骤3具体按照以下步骤实施：
步骤3.1、假定材料在各个方向的热导率相同，根据插值原理，将(2)式转化为：



式(5)中：(i,j,z)为微观元胞的三维坐标；为微观元胞p(i,j,z)的当前温度值/℃；为微观元胞p(i,j,z)...

【专利技术属性】
技术研发人员：李继红，郭钊，郭宇飞，李保铃，苟川东，张敏，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：陕西;61