CN112992280B - 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 - Google Patents

一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 Download PDF

Info

Publication number
CN112992280B
CN112992280B CN202110279581.7A CN202110279581A CN112992280B CN 112992280 B CN112992280 B CN 112992280B CN 202110279581 A CN202110279581 A CN 202110279581A CN 112992280 B CN112992280 B CN 112992280B
Authority
CN
China
Prior art keywords
alpha
beta
model
phase
ternary alloy
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
Application number
CN202110279581.7A
Other languages
English (en)
Other versions
CN112992280A (zh
Inventor
纪飞飞
周晓刚
郑爱权
石彩华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tonggao Intelligent Manufacturing (Hainan) Co.,Ltd.
Original Assignee
Suzhou Chien Shiung Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Suzhou Chien Shiung Institute of Technology filed Critical Suzhou Chien Shiung Institute of Technology
Priority to CN202110279581.7A priority Critical patent/CN112992280B/zh
Publication of CN112992280A publication Critical patent/CN112992280A/zh
Application granted granted Critical
Publication of CN112992280B publication Critical patent/CN112992280B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like

Abstract

本发明属于材料分子动力学领域,涉及三元合金双相材料分子动力学建模及势函数拟合方法,具体公开了一种Ti‑Al‑V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法。首先,根据α+β双相钛合金高温X‑ray衍射试验将常温下的Ti‑Al‑V三元合金分为Ti‑Al为主的α相和Ti‑V为主的β相,并采用Voronoi细分法构建α和β相的多晶模型;合并α相和β相多晶模型,构建Ti‑Al‑V三元合金α+β多相多晶模型;基于所构建的多相多晶模型采用二元合金2NN‑MEAM势函数拟合原子间相互作用关系,表征Ti‑Al‑V三元合金体系原子相互作用势。

Description

一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势 函数拟合方法
技术领域
本发明属于材料分子动力学领域,涉及三元合金双相材料分子动力学建模及势函数拟合方法,具体公开了一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法。
背景技术
分子动力学是一门结合物理、数学和化学的综合技术。主要依靠粒子间相互作用立场和牛顿运动定律来模拟分子体系的运动,是以统计物理学为桥梁,通过研究物质微观性质把握宏观性能。其势函数和模型结构对模拟的准确性具有很大的影响。
对于Ti-Al-V类的三元合金相,目前主要采取的建模方法是建立Ti原子晶胞模型,然后以一定数量的Al和V原子随机替代晶胞中的部分Ti原子。所建的模型虽然能反映Ti、Al、V各原子成分比,但是无法实现Ti-Al基材料HCP结构和BCC结构并存的材料特性。此外,根据α+β双相钛合金高温X-ray衍射试验可知,在α相中主要为Ti和Al元素,V元素主要集中在β相中,且沿着α相边界分布,当高温(约995℃)时,V原子才能扩散进入α相,而通过原子随机替代的方式无法控制原子分布。因此,基于Al和V元素在不同相分布规律构建α相和β相并存的Ti-Al-V三元合金模型对提高仿真准确性十分必要。
针对Ti-Al-V三元合金体系势函数的开发,目前主要采用多体势拟合。分布对Ti-Al元素采用嵌入原子势EAM,Ti-V之间采用Ti-Ni-V三元合金改进型嵌入原子势,Al-V间采用Morse势。该拟合方式构成的Ti-Al-V三元合金势函数能够较为准确的分析高温辐照效应研究及热力学作用下相变行为,但是对金属材料微观缺陷及受力学作用下的应力应变特性、多相多晶结构相互作用等方面的研究则存在着明显的误差。因此开发能够准确描述Ti-Al-V三元合金受载荷作用下相互作用关系的势函数对研究材料的力学性能及微观机理十分重要。
发明内容
针对上述关于Ti-Al-V三元合金体系建模及势函数拟合方法存在的不足,本发明提供一种新的Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法。本发明根据α+β双相钛合金高温X-ray衍射试验中α相、β相中的元素分布规律,构建了一种α+β双相钛合金几何模型,并根据建模思路确定Ti-Al-V三元合金体系势函数拟合方法。
本发明提供的模型构建及势函数拟合方法具体包括以下步骤:
1)根据试验确定Ti-Al-V三元合金中α相和β相各自的成分构成,并采用Voronoi细分法构建α和β相的多晶模型;本发明根据α+β双相钛合金高温X-ray衍射试验将常温下的Ti-Al-V三元合金分为Ti-Al为主的α相和Ti-V为主的β相。
2)合并α相和β相,构建Ti-Al-V三元合金α+β多相多晶模型;
3)采用二元合金2NN-MEAM势函数拟合原子间相互作用关系,表征Ti-Al-V三元合金体系原子相互作用势。
在步骤1)中,α+β双相钛合金高温X-ray衍射试验结果显示:在α相中主要为Ti和Al元素,V元素主要集中在β相中,且沿着α相边界分布,当高温(约995℃)时,V原子才能扩散进入α相。因此,将常温下Ti-Al-V三元合金分模型拆分为以Ti-Al为主要成分的α相(HCP结构)和Ti-V为主要成分的β相(BCC结构)。
在步骤1)中,采用Voronoi细分法,以相同节点设置方式分别构建HCP-Ti和BCC-Ti多晶模型;计算所构建晶胞模型原子总数,并根据Ti、Al、V元素各自成分占比确定所构建模型中Al和V原子数量分别为为n1和n2;随机选取HCP-Ti晶胞中n1个Ti原子替换为Al原子,BCC-Ti晶胞中的n2个Ti原子替换为V原子,以获得HCP-TiAl和BCC-TiV多晶模型。
在步骤2)中,选取Voronoi部分节点,删除HCP-TiAl多晶模型中节点处晶粒,保留BCC-TiV节点处晶粒,并合并HCP-TiAl和BCC-TiV多晶模型,最终获得HCP+BCC结构的多晶多相Ti-Al-V三元合金晶胞模型。
在步骤3)中,采用2NN-MEAM势函数分别拟合Ti-Ti、Ti-Al、Ti-V、Al-V、Al-Al、V-V二元合金体系相互作用关系,最终表征Ti-Al-V三元合金体系各原子相互作用势。
有益效果:本发明根据α+β双相钛合金高温X-ray衍射试验结果将Ti-Al-V三元合金模型分为由Ti、Al元素组成的α相和由Ti、V元素组成的β相,将三元合金转化为两个二元合金的组合;采用2NN-MEAM势函数分别拟合Ti-Ti、Ti-Al、Ti-V、Al-V、Al-Al、V-V二元合金体系相互作用关系,表征Ti-Al-V三元合金体系各原子相互作用势,实现了双相多晶模型的构建及势函数拟合,简化了分子动力学仿真计算过程,提高了计算精度,从而减少实验次数,节省实践成本。本发明为Ti-Al-V三元合金分子动力学仿真及载荷作用下微观结构演变规律研究奠定基础,也为新开发多相多晶建模及势函数开发提供借鉴。
附图说明
图1为HCP-Ti多晶模型。
图2为BCC-Ti多晶模型。
图3为α+β双相Ti-Al-V多晶模型。
具体实施方式
为了使本发明的、技术方案及优点更加清楚明白,下面结合附图及实施例对本发明做进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例
以Ti6Al4V双相钛合金材料为例具体说明Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法。
在Ti6Al4V双相钛合金中,Ti、Al和V原子占比(wt%)分别为90%、6%和4%,则Ti、Al、V原子个数比例则为86.2%:10.2%:3.6%。
假设晶胞中总原子个数为n,Al原子个数n1=n×10.2%,V原子个数n2=n×3.6%,Ti原子个数为n3=n×86.2%。
1)根据α+β双相钛合金高温X-ray衍射试验将常温下的Ti-Al-V三元合金分为Ti-Al为主的α相和Ti-V为主的β相,并采用Voronoi细分法构建α和β相的多晶模型;
2)合并α相和β相,构建Ti-Al-V三元合金α+β多相多晶模型;
3)采用二元合金2NN-MEAM势函数拟合原子间相互作用关系,表征Ti-Al-V三元合金体系原子相互作用势。
具体地,首先根据α+β双相钛合金高温X-ray衍射试验结果可知,Ti6Al4V材料常温下α相主要由Ti和Al元素组成,β相主要由Ti和V元素组成。
采用ATOMSK建模软件构建HCP-Ti单晶模型(如图1所示)和BCC-Ti单晶模型(如图2所示),采用相同的Voronoi细分法构建HCP-Ti和BCC-Ti多晶,具体多晶晶胞尺寸及节点设置方法如下:
box 200 180 210
node 0 0 0[100][010][001]
node 40 80 60 56-83 45
node 0.8*box 0.6*box 0.9*box[11-1][112][1-10]
node 50 5 60[110][1-10][001]
node 0.75*box 0.3*box 0.45*box-31.4 28.7 90.0
node 60 100 80random。
所构建的HCP_Ti和BCC_Ti多晶模型原子数目分别为nHCP和nBCC,则总原子数n=nHCP+nBCC,分别以n1个Al原子随机替换HCP_Ti中的Ti原子,以n2个V原子替换BCC_Ti中的Ti原子,得到HCP_TiAl多晶模型和BCC_TiV多晶模型。
删除HCP_TiAl多晶中节点1和6处的晶粒,同时删除BCC_TiV多晶中的2-5节点,合并两个多晶晶胞得到Ti-Al-V三元合金α+β双相多晶模型(如图3所示)。
采用2NN-MEAM势函数分别对Ti、Al、V单质进行拟合,最终得Ti-Ti、Al-Al和V-V势函数各参数如表1所示。
表1 2NN-MEAM拟合Ti-Ti、Al-Al和V-V势函数各参数
Figure BDA0002977580230000061
类似地,采用2NN-MEAM势函数分别对Ti-Al、Ti-V和Al-V二元合金势函数拟合,各参数如表2所示。
表2 2NN-MEAM拟合Ti-Al、Ti-V和Al-V二元合金势函数各参数
Figure BDA0002977580230000062
Figure BDA0002977580230000071
注:X,Y分别代表不同元素
表3势函数拟合各参数说明
Figure BDA0002977580230000072
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进,这些改进也应视为本发明的保护范围。

Claims (5)

1.一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建的方法,其特征在于,包括以下步骤:
1)根据试验确定Ti-Al-V三元合金中α相和β相各自的成分构成,包括元素含量和晶胞结构信息,然后根据获得的α相和β相各自的元素含量和晶胞结构信息采用Voronoi细分法分别构建α和β相的多晶模型;
2)合并α相和β相的多晶模型,构建Ti-Al-V三元合金α+β多相多晶模型;
在步骤1)中,采用Voronoi细分法,以相同节点设置方式分别构建HCP-Ti和BCC-Ti多晶模型;计算所构建晶胞模型原子总数n,并根据Ti、Al、V元素各自成分占比确定所构建模型中Al和V原子数量分别为n1和n2;随机选取HCP-Ti晶胞中n1个Ti原子替换为Al原子,BCC-Ti晶胞中的n2个Ti原子替换为V原子,获得HCP-TiAl和BCC-TiV多晶模型。
2.根据权利要求1所述的Ti-Al-V三元合金分子动力学α+β双相几何模型构建的方法,其特征在于,在步骤1)中,据α+β双相钛合金高温X-ray衍射试验将常温下的Ti-Al-V三元合金分为Ti-Al为主的α相和Ti-V为主的β相。
3.根据权利要求1所述的Ti-Al-V三元合金分子动力学α+β双相几何模型构建的方法,其特征在于,在步骤2)中,选取Voronoi部分节点,删除HCP-TiAl多晶模型中所选节点处晶粒,保留BCC-TiV多晶模型中所选节点处晶粒,并合并HCP-TiAl和BCC-TiV多晶模型,获得HCP+BCC结构的多晶多相Ti-Al-V三元合金晶胞模型。
4.一种势函数拟合方法,其特征在于,根据权利要求1所构建得到的Ti-Al-V三元合金α+β多相多晶模型,采用二元合金2NN-MEAM势函数拟合原子间相互作用关系,表征Ti-Al-V三元合金体系原子相互作用势。
5.一种势函数拟合方法,其特征在于,根据权利要求3所构建得到的HCP+BCC结构的多晶多相Ti-Al-V三元合金晶胞模型,采用2NN-MEAM势函数分别拟合Ti-Ti、Ti-Al、Ti-V、Al-V、Al-Al、V-V二元合金体系相互作用关系,最终表征Ti-Al-V三元合金体系各原子相互作用势。
CN202110279581.7A 2021-03-16 2021-03-16 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 Active CN112992280B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110279581.7A CN112992280B (zh) 2021-03-16 2021-03-16 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110279581.7A CN112992280B (zh) 2021-03-16 2021-03-16 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法

Publications (2)

Publication Number Publication Date
CN112992280A CN112992280A (zh) 2021-06-18
CN112992280B true CN112992280B (zh) 2022-01-11

Family

ID=76335775

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110279581.7A Active CN112992280B (zh) 2021-03-16 2021-03-16 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法

Country Status (1)

Country Link
CN (1) CN112992280B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113782104A (zh) * 2021-08-20 2021-12-10 西北工业大学 镍基单晶合金模型构建方法、装置、存储介质及设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103714579A (zh) * 2013-12-11 2014-04-09 燕山大学 一种建立多相多晶原子结构模型的方法
CN105658372A (zh) * 2013-09-05 2016-06-08 萨澳激光有限公司 通过激光束加工工件的方法、激光刀具、激光加工机、及机器控制
CN105772689A (zh) * 2016-03-03 2016-07-20 大连理工大学 基于分子动力学的铸造高铬合金的建模方法
CN108411156A (zh) * 2018-03-14 2018-08-17 中南大学 一种近β型高强度钛合金及其制备方法
CN110489934A (zh) * 2019-09-20 2019-11-22 长春理工大学 一种单晶及多晶钛内微裂纹扩展的分子动力学构建方法
CN111414693A (zh) * 2020-03-19 2020-07-14 中国矿业大学 一种基于分子动力学的异质结材料界面结合能测量方法
CN111785831A (zh) * 2020-06-23 2020-10-16 武汉理工大学 一种提升In2Se3相变材料多值存储特性的方法
CN112487568A (zh) * 2020-11-12 2021-03-12 中国科学院金属研究所 预测双相钛合金在不同升温速率下组织形貌演化与合金元素配分的相场模拟方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101307417B (zh) * 2007-05-14 2010-12-22 蒋铭瑞 高强度钛合金高尔夫杆头构件的制造方法及其产品
CN102750425B (zh) * 2012-07-17 2014-09-03 哈尔滨工业大学 一种焊接过程热影响区组织演变的模拟方法
CN107034371B (zh) * 2013-05-21 2020-06-19 麻省理工学院 稳定的纳米晶有序合金体系及其鉴定方法
CN111326220B (zh) * 2020-04-16 2023-08-15 重庆大学 一种高强韧锆钛基合金的设计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105658372A (zh) * 2013-09-05 2016-06-08 萨澳激光有限公司 通过激光束加工工件的方法、激光刀具、激光加工机、及机器控制
CN103714579A (zh) * 2013-12-11 2014-04-09 燕山大学 一种建立多相多晶原子结构模型的方法
CN105772689A (zh) * 2016-03-03 2016-07-20 大连理工大学 基于分子动力学的铸造高铬合金的建模方法
CN108411156A (zh) * 2018-03-14 2018-08-17 中南大学 一种近β型高强度钛合金及其制备方法
CN110489934A (zh) * 2019-09-20 2019-11-22 长春理工大学 一种单晶及多晶钛内微裂纹扩展的分子动力学构建方法
CN111414693A (zh) * 2020-03-19 2020-07-14 中国矿业大学 一种基于分子动力学的异质结材料界面结合能测量方法
CN111785831A (zh) * 2020-06-23 2020-10-16 武汉理工大学 一种提升In2Se3相变材料多值存储特性的方法
CN112487568A (zh) * 2020-11-12 2021-03-12 中国科学院金属研究所 预测双相钛合金在不同升温速率下组织形貌演化与合金元素配分的相场模拟方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"双相钛合金高温变形协调性的CPFEM研究";李学雄等;《金属学报》;20190731;第55卷(第7期);第929页左栏第1-3段和第930-931页第1.3节 *
"激光立体成形Ti-6Al-4V时效态合金显微组织结构研究";赵学平等;《内蒙古工业大学学报》;20191231;第38卷(第3期);第170页第2段 *
"镍基高温合金中γ/θ相界面性能的数值模拟";董卫平等;《稀有金属材料与工程》;20190531;第48卷(第5期);第1529页第1节 *

Also Published As

Publication number Publication date
CN112992280A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
Van De Walle et al. The alloy theoretic automated toolkit: A user guide
Craievich et al. Structural instabilities of excited phases
CN112992280B (zh) 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法
CN107665274B (zh) 一种设计低弹性模量钛合金的方法
Yu et al. A two-stage predicting model for γ′ solvus temperature of L12-strengthened Co-base superalloys based on machine learning
Blum et al. Prediction of ordered structures in the bcc binary systems of Mo, Nb, Ta, and W from first-principles search of approximately 3,000,000 possible configurations
Li et al. Anticontrol of chaos for discrete TS fuzzy systems
Vela et al. High-throughput exploration of the WMoVTaNbAl refractory multi-principal-element alloys under multiple-property constraints
Goiri et al. Recursive alloy Hamiltonian construction and its application to the Ni-Al-Cr system
Ishida et al. MD simulation of martensitic transformations in TiNi alloys with MEAM
Yang et al. Hydride corrosion kinetics on metallic surface: a multiphase-field modeling
Göke AnyMOD–A graph-based framework for energy system modelling with high levels of renewables and sector integration
Ghazisaeidi Alloy thermodynamics via the Multi-cell Monte Carlo (MC) 2 method
Lu et al. An accurate interatomic potential for the TiAlNb ternary alloy developed by deep neural network learning method
Samolyuk et al. Thermodynamic approach to the stability of multi-phase systems: application to the Y2O3–Fe system
CN114974474A (zh) 结构材料的多尺度链条式集成设计方法及数据库系统
CN114169220A (zh) 基于可相变Fe50Mn30Cr10Co10高熵合金相变过程的原子级模拟方法
Luo et al. First principles calculation of composition dependence tracer and interdiffusion with phase change in γ/γ′ superalloy: A case study of Ir/Ir3Nb
Kulesza Modified, threshold-based circulation type classification for Central Europe, on the basis of Lityński’s classification
Ladd et al. High‐pressure mechanical instability in close‐packed Hooke’s‐law crystals
Calvo et al. On the convergence of global thermal properties of clusters extracted from simulations
Abyaneh Kinetics of two-phase electrocrystallization processes: I. Competitive growth of right-circular cones
Gadre et al. Reduced first-order density matrices and ‘‘exchange-only’’correlation factors for some closed-shell atomic systems
Nandipati et al. Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in a-Zr Grain Boundaries
Morrison Understanding the Interaction between Grain Boundaries and Precipitates in Ni-Al Using Molecular Dynamics

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
TR01 Transfer of patent right

Effective date of registration: 20240415

Address after: A10100, 4th Floor A125-88, Joint Inspection Building, Haikou Free Trade Zone, Hainan Province

Patentee after: Tonggao Intelligent Manufacturing (Hainan) Co.,Ltd.

Country or region after: China

Address before: 215411 Jianxiong Road, Taicang Science and Education New Town, Suzhou City, Jiangsu Province

Patentee before: SUZHOU CHIEN-SHIUNG INSTITUTE OF TECHNOLOGY

Country or region before: China