一种基于余氏理论的钛合金β相高温固溶组织模拟方法

一种基于余氏理论的钛合金
β
相高温固溶组织模拟方法
技术领域
1.本发明属于钛合金组织模拟领域，具体涉及一种基于余氏理论的钛合金β相高温固溶组织模拟方法。


背景技术：

2.钛及钛合金具有比强度高、耐腐蚀、生物相容性好以及优异的高温稳定性等性能，被广泛应用于航空航天、化学工业、舰船、生物医疗等领域，钛合金组织性能研究也备受关注。当前，物理冶金模型、神经元网络、元胞自动机、相场模拟等广泛应用到钢铁合金组织模拟中，这些模拟方法因其各自局限性，在指导材料成分—制备工艺—组织间的优化匹配中仍很难广泛应用，使得钛合金中组织模拟方法研究仍不深入。另外,基于余氏理论的钛合金组织模拟目前更是没有相关报道。


技术实现要素：

3.本发明的目的在于提供一种基于余氏理论的钛合金β相高温固溶组织模拟方法，可在计算机上实现组织模拟与晶粒尺寸测量，实现钛合金成分—制备工艺—组织间的优化匹配，解决了现有技术中钛合金组织性能研究存在局限性的问题。
4.本发明技术方案如下：
5.一种基于余氏理论的钛合金β相高温固溶组织模拟方法，包括如下步骤：
6.步骤1、将钛合金成分的质量百分数转化为原子百分数；
7.步骤2、利用余氏理论计算钛合金中β
‑
ti
‑
al、β
‑
ti
‑
al
‑
m、β
‑
ti
‑
m原子团簇电子结构，即结合能、最强共价键上的公用电子对数；
[0008][0009]
其中其中c
al
为al的原子分数，c
m
为不包括al在内的合金原子分数，z为合金元素种类数；p为参数，与al的含量有关，p＝1.0～2.0；ap代表各原子团簇的结合能或最强共价键上的公用电子对数；
[0010]
步骤3、利用β相高温原子团簇可能比例份额计算元胞空间nx
×
ny中β相形核点数目，计算公式为
[0011][0012]
其中，n
j
为j种原子团簇的不考虑温度的形核点数目；f
j
为j种原子团簇的可能比例份额，见公式(1)；
[0013]
利用原子团簇结合能建立不同固溶温度下各原子团簇在元胞空间nx
×
ny中的形核点数目，计算公式为
[0014][0015]
其中s
j
为j种原子团簇不同固溶温度的形核点数目；k
p
＝1000～1444；t为固溶温度，k；ec为j种原子团簇结合能；r为摩尔气体常数；
[0016]
步骤4、将步骤3计算的不同固溶温度的形核点数目按照生命游戏规则进行生长，在晶粒生长过程中，根据原子团簇结合能或最强共价键上的共用电子对数的大小确定形核点长大与消失的概率，并进行排序；
[0017]
步骤5、将公式(1)至(3)以及生命游戏规则编制成计算机程序，向计算机程序中输入步骤3所得的形核点数目，然后按照计算步长范围进行模拟，模拟出钛合金β相高温固溶组织图。
[0018]
所述的步骤4中，将形核点能长大的形核数目归入到形核点生长模块，将形核点消失的形核数目归入到形核点消失模块。
[0019]
所述的步骤5中，同时获得晶粒尺寸的统计结果。
[0020]
本发明的优点及有益效果是：
[0021]
1.本发明可以在计算机上实现钛合金不同固溶温度的组织模拟，这为优化钛合金固溶处理温度参数提供了较好的方法手段。
[0022]
2.本发明可以在计算机上开展不同计算步长的组织模拟，从而实现钛合金在不同固溶时间的组织模拟，这为优化钛合金固溶处理时间提供了较好的方法手段。
[0023]
3.本发明利用计算机可以模拟钛合金不同固溶处理工艺参数下的模拟平均晶粒尺寸，晶粒尺寸模拟结果与实际符合较好。
[0024]
4.本发明具有操作简单、计算成本低、精度高的特点。
附图说明
[0025]
图1a为tc4合金1000℃固溶0.5h时的金相显微图；
[0026]
图1b为tc4合金1000℃固溶0.5h时的形核点数目分布图；
[0027]
图1c为tc4合金1000℃固溶0.5h时的模拟的组织图；
[0028]
图1d为tc4合金1000℃固溶0.5h时的晶界处的形核点数目分布图；
[0029]
图2a为tc4合金1050℃固溶0.5h时的金相显微图；
[0030]
图2b为tc4合金1050℃固溶0.5h时的形核点数目分布图；
[0031]
图2c为tc4合金1050℃固溶0.5h时的模拟的组织图；
[0032]
图2d为tc4合金1050℃固溶0.5h时的晶界处的形核点数目分布图。
[0033]
图3a为tc4合金1100℃固溶0.5h时的金相显微图；
[0034]
图3b为tc4合金1100℃固溶0.5h时的形核点数目分布图；
[0035]
图3c为tc4合金1100℃固溶0.5h时的模拟的组织图；
[0036]
图3d为tc4合金1100℃固溶0.5h时的晶界处的形核点数目分布图。
具体实施方式
[0037]
以下结合实例对本发明进行详细说明。
[0038]
一种基于余氏理论的钛合金β相高温固溶组织模拟方法，包括如下步骤：
[0039]
步骤1、将钛合金成分的质量百分数转化为原子百分数；
[0040]
步骤2、利用余氏理论计算钛合金中β
‑
ti
‑
al、β
‑
ti
‑
al
‑
m、β
‑
ti
‑
m原子团簇电子结构，即结合能、最强共价键上的公用电子对数；钛合金β相高温固溶可能存在β
‑
ti
‑
al、β
‑
ti
‑
al
‑
m、β
‑
ti
‑
m原子团簇，这些原子团簇出现的概率与合金原子m含量及原子间竞争al原子的能力有关，原子间竞争al的能力可用合金原子m与al原子间相互作用强弱表征。根据余氏理论的自洽键距差可以计算这些原子团簇中原子的相互作用，即结合能、最强共价键上的公用电子对数；
[0041][0042]
其中其中c
al
为al的原子分数，c
m
为不包括al在内的合金原子分数，z为合金元素种类数；p为参数，与al的含量有关，p＝1.0～2.0；ap代表各原子团簇的结合能或最强共价键上的公用电子对数；
[0043]
步骤3、利用β相高温原子团簇可能比例份额计算元胞空间nx
×
ny中β相形核点数目，计算公式为
[0044][0045]
其中，n
j
为j种原子团簇的不考虑温度的形核点数目；f
j
为j种原子团簇的可能比例份额，见公式(1)；
[0046]
利用原子团簇结合能建立不同固溶温度下各原子团簇在元胞空间nx
×
ny中的形核点数目，计算公式为
[0047][0048]
其中s
j
为j种原子团簇不同固溶温度的形核点数目；k
p
＝1000～1444；t为固溶温度，k；ec为j种原子团簇结合能；r为摩尔气体常数；
[0049]
步骤4、将步骤3计算的不同固溶温度的形核点数目按照生命游戏规则进行生长，在晶粒生长过程中，根据原子团簇结合能或最强共价键上的共用电子对数的大小确定形核点长大与消失的概率，并进行排序；将形核点能长大的形核数目归入到形核点生长模块，将形核点消失的形核数目归入到形核点消失模块；
[0050]
步骤5、将公式(1)至(3)以及生命游戏规则编制成计算机程序，向计算机程序中输入步骤3所得的形核点数目，然后按照计算步长范围进行模拟，模拟出钛合金β相高温固溶组织图，同时获得晶粒尺寸的统计结果。
[0051]
应用实施例1：以tc4合金1000℃固溶0.5h的组织模拟为例。
[0052]
以tc4合金为例，tc4钛合金的主要成分见表1。
[0053]
步骤1.将tc4合金的质量百分数转换为原子分数，见表2。
[0054]
步骤2.利用余氏理论计算出不同原子团簇的结合能和最强键上共用电子对数(简称电子对数)，见表3。
[0055]
步骤3.利用公式(1)及表2、3的数据计算出不同原子团簇在β相区固溶时的可能比例份额，计算结果为β
‑
ti
‑
al
‑
fe，0.22696；β
‑
ti
‑
al，6.51453；β
‑
ti
‑
fe，0.09854；β
‑
ti
‑
al
‑
v，2.69173；β
‑
ti
‑
v，0.69858。在计算过程中，p取1.0。
[0056]
步骤4.根据公式(2)利用β相高温原子团簇可能比例份额计算元胞空间50
×
50中β相形核点数目(不考虑温度的影响)，分别为:n
β
‑
ti
‑
al
‑
fe
，11；n
β
‑
ti
‑
al
，326；n
β
‑
ti
‑
fe
，5；n
β
‑
ti
‑
al
‑
v
，135；n
β
‑
ti
‑
v
，35。
[0057]
步骤5.根据公式(3)利用表3中结合能建立1000℃固溶温度下各原子团簇在元胞空间50
×
50中的形核点数目，计算时k取1313，各原子团簇形核点数目分别为：s
β
‑
ti
‑
al
‑
fe
，11；s
β
‑
ti
‑
al
，324；s
β
‑
ti
‑
fe
，5；s
β
‑
ti
‑
al
‑
v
，134；s
β
‑
ti
‑
v
，35。
[0058]
步骤6.按照β相固溶晶粒生长规则：将步骤5计算的不同固溶温度的形核点数目按照生命游戏规则进行生长。在晶粒生长过程中，可根据原子团簇结合能的大小进行选择形核点的生长与消失。
[0059]
步骤7.将公式(1)至(3)以及生命游戏规则编制成计算机程序，向计算机程序中输入步骤5计算的形核点数目，然后按照计算步长50
‑
1400范围进行模拟；确定计算步长1250的模拟组织与实际组织比较吻合，模拟组织的模拟尺寸长为14.667元胞，实际组织尺寸为446.299μm，见图1a
‑
1d。
[0060]
实例2：以tc4合金1050℃固溶0.5h的组织模拟为例。
[0061]
按照应用实施例1中的步骤1
‑
4，计算tc4合金的原子团簇形核点数目(不考虑温度)。
[0062]
利用步骤1
‑
4计算结果及公式(3)计算tc4合金在1050℃固溶时元胞空间中50
×
50中的形核点数目，即s
β
‑
ti
‑
al
‑
fe
，11：s
β
‑
ti
‑
al
＝，312；s
β
‑
ti
‑
fe
＝，5；s
β
‑
ti
‑
al
‑
v
，129；s
β
‑
ti
‑
v
，34，其中k取1313。
[0063]
在计算机程序中输入本实例步骤2的形核点数目，计算步长为1295时模拟组织与
实际组织较符合，见图2a
‑
2d。按照14.6667：446.299比例尺，模拟组织尺寸为669.4470μm；实际组织为825.0845μm，模拟尺寸误差为11.69％。
[0064]
实例3：以tc4合金1100℃固溶0.5h的组织模拟为例。
[0065]
按照应用实施例1中的步骤1
‑
4，可以计算tc4合金的原子团簇形核点数目(不考虑温度)。
[0066]
利用步骤1
‑
4计算结果及公式(3)计算tc4合金在1100℃固溶时元胞空间中50
×
50中的形核点数目，即s
β
‑
ti
‑
al
‑
fe
，10；s
β
‑
ti
‑
al
，301；s
β
‑
ti
‑
fe
，5；s
β
‑
ti
‑
al
‑
v
，＝124；s
β
‑
ti
‑
v
，33。
[0067]
在计算机程序中输入本实例步骤2的形核点数目，计算步长1216时模拟组织与实际组织较符合，见图3a
‑
3d。按照14.6667：446.299比例尺，模拟组织尺寸为745.5205μm；实际组织为825.0845μm，模拟尺寸误差为9.64％。
[0068]
本发明还包括一种模拟方法所用设备，包括存储器和处理器，存储器用于存储计算机程序；处理器用于执行所述计算机程序时实现上述存储介质模拟方法的步骤。
[0069]
不同热处理工艺的钛合金具有不同的性能，而钛合金的性能依赖于合金成分及不同尺度的微观组织，合理调控合金微观组织是进行合金性能优化与调控的最有效手段之一。本发明的模拟方法可以在计算机上进行合金成分、热处理工艺及组织状态之间的内在关系的匹配设计，为钛合金组织性能优化与新合金设计提供方法手段。
[0070]
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制。凡是根据发明技术实质对以上实施例所作的任何简单修改、变更以及等效变化，均仍属于本发明技术方案的保护范围内。
[0071]
表1 ti
‑
6al
‑
4v主要成分(wt％)
[0072][0073]
表2 ti
‑
6al
‑
4v主要成分(at％)
[0074][0075]
表3原子团簇的结合能
[0076]