CN113571139A - 一种金属晶界建模方法、记录媒体及系统 - Google Patents
一种金属晶界建模方法、记录媒体及系统 Download PDFInfo
- Publication number
- CN113571139A CN113571139A CN202110833111.0A CN202110833111A CN113571139A CN 113571139 A CN113571139 A CN 113571139A CN 202110833111 A CN202110833111 A CN 202110833111A CN 113571139 A CN113571139 A CN 113571139A
- Authority
- CN
- China
- Prior art keywords
- grain boundary
- model
- metal
- surface model
- metal grain
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 229910052751 metal Inorganic materials 0.000 title claims abstract description 31
- 239000002184 metal Substances 0.000 title claims abstract description 31
- 239000013078 crystal Substances 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims abstract description 11
- 230000007547 defect Effects 0.000 claims abstract description 8
- 238000013459 approach Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 230000002040 relaxant effect Effects 0.000 claims description 2
- 239000000463 material Substances 0.000 abstract description 18
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 6
- 230000001052 transient effect Effects 0.000 abstract description 3
- 230000000704 physical effect Effects 0.000 abstract description 2
- 238000012827 research and development Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000003860 storage Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 102100021164 Vasodilator-stimulated phosphoprotein Human genes 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005610 quantum mechanics Effects 0.000 description 3
- 238000005204 segregation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004883 computer application Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229910052725 zinc Inorganic materials 0.000 description 1
Images
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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
本发明属于材料设计技术领域,特别涉及一种金属晶界建模方法,包括以下步骤:计算金属的晶格参数;构建原子取向表面模型;通过弛豫表面构型,获取原子配位信息,并进而确定合适的表面模型层数;通过对称镜面操作构建初始晶界模型;通过设置应力、缺陷和滑移等操作,改变晶界构型,在此基础上优化晶界模型,获得最终晶界结构信息,并据此进行材料模拟获取晶界的物理性质。通过该方法,可以迅速有效地获取不同晶界的结构信息,适用于材料的研发和设计,具有很强的应用价值。本发明还提供了一种存储有金属晶界建模程序的非暂态可读记录媒体及包含该媒体的装置,通过处理电路可以调用该程序,以执行上述方法。
Description
技术领域
本发明属于材料设计技术领域,公开了一种金属晶界建模方法、存储有能执行该方法程序的记录媒体及系统。
背景技术
材料在使用过程中大部分都是以多晶状态存在的。多晶体中相邻的晶体间由于原子各有不同的取向,从而形成晶界。晶界可认为是一种晶体缺陷。它结构复杂,原子排列不规则,同时不可避免地存在着杂质和合金添加剂等偏聚元素,有时候会因此导致材料出现晶界脆化现象,严重影响材料的物理性能。
晶界性质的研究长期以来一直是材料科学工程中最重要的领域之一,具有十分重要的理论和实践意义。为了改进材料性能,提高其使役寿命,大量的科研工作开始针对材料晶界的性质展开研究,如晶界滑移、晶界偏析以及脆裂等。伴随着计算机应用技术的不断提升,直接用计算机模拟材料的晶界结构成为一种趋势。它能够直接从原子层次构建模型,准确地描述出材料的结构与性能的对应关系。尤其是精度更高的量子力学第一原理计算方法,它不使用除基本物理常数和原子量以外的任何可调的经验和拟合参数,直接对体系的薛定谔方程进行求解,从而得到体系的总能和电子结构,据此实现对材料的性能的预测分析。
在晶界的模拟研究过程中,获取精确结构的位错结构对于性能预测的准确性至关重要。其次,合适的晶胞尺寸也能节约计算资源,加速研发过程的实施。由于晶界结构的错综复杂,无论是实验技术,还是材料模拟仿真技术对于准确获取晶界的结构信息都是非常困难的,因此设计一种方法可以快速有效地获取不同金属晶界处的原子结构信息对材料性质的表征具有重大的实用价值。
发明内容
本发明的目的在于提供一种金属晶界建模方法,采用该筛选方法可以简单快捷地获取不同金属不同晶体取向的晶界结构信息,模拟结果可与实验测量进行比较,该方法可大幅度降晶界模拟操作的复杂性,准确表征出晶界的性质,对于通过晶界调控提升材料的物理性能有着重要的促进作用。
为实现上述目的,本发明采用以下技术方案:
一种金属晶界建模方法,包括以下步骤:
(1)计算金属的晶格参数获取原胞结构信息;
(2)根据所述原胞晶体结构,构建原子取向的表面模型;
(3)弛豫所述表面模型,获取所述表面模型不同层数的原子配位信息,确定所述表面模型的层数n;
(4)根据所述层数n,重构表面模型,通过对称镜面操作,去除重复原子层,获取初始晶界模型,晶界模型层数为2n-1;
(5)优化晶界模型。
优选的,所述步骤(2)中,运用三维建模软件Materials Studio或编程方法获取周期性重复的所述表面模型。
优选的,所述步骤(5)中优化是指通过设置应力、缺陷和滑移等操作,使所述晶界模型接近于金属晶界的真实状况。
在本发明的方法中,晶界结构错综复杂,通常选取实验中观测到的某种取向的稳定晶界建模,所构建的晶界为孪晶。构建过程中需保证拓展原胞每层至少含有4个原子,这样获取的原子配位信息受晶胞尺寸的影响较小,模型才会更加准确。对于表面模型,一个晶粒上下都是真空层;对于晶界模型,与一个晶粒连接的是另一种晶体取向的晶粒。基于表面模型建立晶界模型,更加方便有效。另外,通过设置应力、缺陷和滑移等操作,使得体系与真实更加接近。
本发明的另一方案在于提供一种非暂态可读记录媒体,用以存储包含多个指令的一个或多个程序,其特征在于,当执行指令时,将致使处理电路执行上述的一种金属晶界建模方法。
本发明的又一方案在于提供一种金属晶界建模系统,包括处理电路及与其电性耦接的存储器,其特征在于,所述存储器配置储存至少一程序,所述程序包含多个指令,所述处理电路运行所述程序,能执行上述的一种金属晶界建模方法。
本发明的优点在于:
本发明建立一种金属晶界建模方法,通过该方法可获取稳定合理的晶界结构信息。以往的实验和模拟研究并没有提出完整有效的设计方案。本发明设计基于表面模型设计出一种快速有效的异质界面建模方法,其很大程度简化了晶界体系模拟研究的复杂性。应用该方法及包括该方法程序的系统可以合理地构建出文献中实验和计算报道的晶界以及未知的其它晶界,为材料的研发和设计节约投入成本。
附图说明
图1为Al金属的(210)表面的原子排布示意图;
包括侧视(a)和俯视(b)两部分;其中大球和小球分别代表不同平面的原子,球上的数字为原子层数示意。
图2为Al∑5(210)[001]晶界的侧视图;
其中虚线上下为不同排列的晶粒。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行描述,所描述的实施例是本发明一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域技术人员在没有做出创新劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
以Al∑5(210)[001]孪晶为例,一种金属晶界建模方法包括如下步骤:
2、基于面心立方的Al金属的晶体结构,利用三维建模软件Materials Studio构建(210)面堆垛的金属表面结构,如图1所示。该原胞为(2×2)的拓展结构,每层含有4个原子,共计60个原子。原子层厚度设置为15层。上下表面各厚的真空层,总计真空层厚度为
3、利用计算软件VASP弛豫Al(210)表面模型。获取不同层数的原子配位信息,结果如表1所示。表面的1-3层配位数低于体相的12个,与体相情况相差甚远;表面4-5层近邻原子配位数达到12,但是近邻距离与体相的仍然有差距。表面6-8层基与体相的配位情况基本一致,尤其是第8层的原子,它的12个近邻原子距离为几乎跟体相一样。
据此,可以认为8层的表面模拟最合适,6-7层的表面模型较为合适。
4、基于以上分析结果,确定最合适的表面模型厚度参数为n=8。据此利用Materials Studio软件构建8层(210)面堆垛的铝金属表面模型,获得晶粒1,如图2所示。对晶粒1进行对称镜面操作获得晶粒2。在此结构上去除重复的1层原子,从而获得厚度为17层的初始晶界模型。
5、将结构导入量子力学第一原理模拟计算软件VASP。通过优化上述晶界模型,获得没有其它缺陷的晶界结构,从而准确表征出Al∑5(210)[001]孪晶晶界周围原子的近邻配位信息。
6、根据实际需要,对于1-8层不同位置的原子,设置空位或者掺杂元素,如Mg和Zn等,构建有缺陷偏聚的晶界构型,重新获得含有其它缺陷体系的晶界原子的构型和近邻配位信息;另外,可以通过调节最外层(靠近真空层)原子的位置并固定该层原子,使得界面体系承受拉应力和压应力,获得不同应力条件下晶界原子的构型和近邻配位信息;最后,通过晶粒1和晶粒2的相对移动操作,可以获得滑移后晶界原子的构型和近邻配位信息。
利用该方法可以快速有效地获取不同条件下Al金属晶界处的原子结构信息,对后期Al金属材料原子层次性质的预测有着重要的意义。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机、可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
将上述方法步骤汇编成程序再存储于硬盘或其他非暂态存储介质就构成了本发明的“一种非暂态可读记录媒体”技术方案;而将该存储介质与计算机处理器电连接,通过程序实施金属晶界建模方法,则构成本发明的“一种金属晶界建模系统”技术方案。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种金属晶界建模方法,其特征在于,包括以下步骤:
(1)计算金属的晶格参数获取原胞结构信息;
(2)根据所述原胞晶体结构,构建原子取向的表面模型;
(3)弛豫所述表面模型,获取所述表面模型不同层数的原子配位信息,确定所述表面模型的层数n;
(4)根据所述层数n,重构表面模型,通过对称镜面操作,去除重复原子层,获取初始晶界模型,晶界模型层数为2n-1;
(5)优化晶界模型。
2.根据权利要求1所述的一种金属晶界建模方法,其特征在于,所述步骤(2)中,运用三维建模软件Materials Studio或编程方法获取周期性重复的所述表面模型。
3.根据权利要求1所述的一种金属晶界建模方法,其特征在于,所述步骤(5)中优化是指通过设置应力、缺陷和滑移等操作,使所述晶界模型接近于金属晶界的真实状况。
4.一种非暂态可读记录媒体,用以存储包含多个指令的一个或多个程序,其特征在于,当执行指令时,将致使处理电路执行权利要求1-3中任一项所述的一种金属晶界建模方法。
5.一种金属晶界建模系统,包括处理电路及与其电性耦接的存储器,其特征在于,所述存储器配置储存至少一程序,所述程序包含多个指令,所述处理电路运行所述程序,能执行权利要求1-3中任一项所述的一种金属晶界建模方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110833111.0A CN113571139B (zh) | 2021-09-22 | 2021-09-22 | 一种金属晶界建模方法、记录媒体及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110833111.0A CN113571139B (zh) | 2021-09-22 | 2021-09-22 | 一种金属晶界建模方法、记录媒体及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113571139A true CN113571139A (zh) | 2021-10-29 |
CN113571139B CN113571139B (zh) | 2024-03-12 |
Family
ID=78166465
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110833111.0A Active CN113571139B (zh) | 2021-09-22 | 2021-09-22 | 一种金属晶界建模方法、记录媒体及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113571139B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109086507A (zh) * | 2018-07-24 | 2018-12-25 | 北京航空航天大学 | 一种评价含复合点缺陷镍基合金变形能力的方法 |
CN109101780A (zh) * | 2018-07-24 | 2018-12-28 | 北京航空航天大学 | 一种评价稀土原子与空位在镍基合金双相界面处相互作用的方法 |
CN109187337A (zh) * | 2018-09-10 | 2019-01-11 | 南京工业职业技术学院 | 一种筛选强韧性FeAl晶界的方法 |
CN109830268A (zh) * | 2019-02-16 | 2019-05-31 | 辽宁石油化工大学 | 一种倾扭混合型晶界原子结构的建模方法 |
CN111354422A (zh) * | 2018-12-24 | 2020-06-30 | 有研工程技术研究院有限公司 | 一种微量溶质元素晶界偏聚的筛选方法 |
CN112053745A (zh) * | 2020-07-27 | 2020-12-08 | 西南交通大学 | 获取双晶晶界能的方法 |
CN112071371A (zh) * | 2020-08-28 | 2020-12-11 | 武汉大学 | 一种用于物质和材料的计算机仿真模拟方法及装置 |
US20210118530A1 (en) * | 2019-05-27 | 2021-04-22 | Beijing University Of Technology | Multi-scale method for simulating mechanical behaviors of multiphase composite materials |
US20210217501A1 (en) * | 2020-01-10 | 2021-07-15 | Fujitsu Limited | Learning device and learning method |
CN113393912A (zh) * | 2021-06-23 | 2021-09-14 | 湖南大学 | 一种双相结构高熵合金的强度预测方法 |
-
2021
- 2021-09-22 CN CN202110833111.0A patent/CN113571139B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109086507A (zh) * | 2018-07-24 | 2018-12-25 | 北京航空航天大学 | 一种评价含复合点缺陷镍基合金变形能力的方法 |
CN109101780A (zh) * | 2018-07-24 | 2018-12-28 | 北京航空航天大学 | 一种评价稀土原子与空位在镍基合金双相界面处相互作用的方法 |
CN109187337A (zh) * | 2018-09-10 | 2019-01-11 | 南京工业职业技术学院 | 一种筛选强韧性FeAl晶界的方法 |
CN111354422A (zh) * | 2018-12-24 | 2020-06-30 | 有研工程技术研究院有限公司 | 一种微量溶质元素晶界偏聚的筛选方法 |
CN109830268A (zh) * | 2019-02-16 | 2019-05-31 | 辽宁石油化工大学 | 一种倾扭混合型晶界原子结构的建模方法 |
US20210118530A1 (en) * | 2019-05-27 | 2021-04-22 | Beijing University Of Technology | Multi-scale method for simulating mechanical behaviors of multiphase composite materials |
US20210217501A1 (en) * | 2020-01-10 | 2021-07-15 | Fujitsu Limited | Learning device and learning method |
CN112053745A (zh) * | 2020-07-27 | 2020-12-08 | 西南交通大学 | 获取双晶晶界能的方法 |
CN112071371A (zh) * | 2020-08-28 | 2020-12-11 | 武汉大学 | 一种用于物质和材料的计算机仿真模拟方法及装置 |
CN113393912A (zh) * | 2021-06-23 | 2021-09-14 | 湖南大学 | 一种双相结构高熵合金的强度预测方法 |
Non-Patent Citations (4)
Title |
---|
杨志勇;李明哲;郝义磊;刘凯俐;司慧玲;刘冬梅: "黄铜晶间腐蚀机理的实验及模拟研究", 表面技术, vol. 47, no. 8, pages 1 - 3 * |
熊志华;刘继东;雷敏生;: "铜薄膜中空位的计算机模拟", 南昌大学学报(理科版), no. 02 * |
马帅;李拥华;高裕博;: "Ca对氧化铝晶界处氧空位扩散的活化机理", 上海大学学报(自然科学版), no. 04 * |
马文;祝文军;张亚林;陈开果;邓小良;经福谦;: "纳米多晶金属样本构建的分子动力学模拟研究", 物理学报, no. 07 * |
Also Published As
Publication number | Publication date |
---|---|
CN113571139B (zh) | 2024-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fritzen et al. | Computational homogenization of elasto-plastic porous metals | |
Papadakis et al. | Numerical computation of component shape distortion manufactured by selective laser melting | |
KR20200023164A (ko) | 원자 요소들 및 그 합금 물질들의 구조와 성질들을 예측하기 위한 시스템들 및 방법들 | |
Barrett et al. | An automated procedure for geometry creation and finite element mesh generation: Application to explicit grain structure models and machining distortion | |
Prakash et al. | A multiscale simulation framework of the accumulative roll bonding process accounting for texture evolution | |
Tucker et al. | Tail departure of log-normal grain size distributions in synthetic three-dimensional microstructures | |
Sakane et al. | Parallel-GPU-accelerated adaptive mesh refinement for three-dimensional phase-field simulation of dendritic growth during solidification of binary alloy | |
CN113642666B (zh) | 一种基于样本扩充和筛选的主动增强软测量方法 | |
CN108038333B (zh) | 一种大型盘铣刀平面铣削的柔性加工表面误差的预测方法 | |
March-Rico et al. | Defect cluster configurations and mobilities in α-Zr: A comparison of the BMD19 and M07 interatomic potentials | |
Duchêne et al. | Numerical simulation of a pyramid steel sheet formed by single point incremental forming using solid-shell finite elements | |
CN113571139B (zh) | 一种金属晶界建模方法、记录媒体及系统 | |
Henrich et al. | DRAGen–A deep learning supported RVE generator framework for complex microstructure models | |
Klippel | Constitutive equations for simulation of metal cutting with meshless methods on GPU | |
Saunders et al. | Influence of Grain Size and Shape of Mechanical Properties of Metal AM Materials | |
CN116525038A (zh) | 一种水凝胶最优合成方法、系统、设备及介质 | |
CN111354422A (zh) | 一种微量溶质元素晶界偏聚的筛选方法 | |
Chen et al. | Constitutive modeling of slip, twinning and detwinning for mg alloy and inhomogeneous evolution of microstructure | |
Perzyński | The effect of model size and boundary conditions on the representativeness of digital material representation simulations of ferritic-pearlitic sample compression | |
Li et al. | A method for object reconstruction based on point-cloud data via 3D scanning | |
Tokarev et al. | Multi-level direct crystal plasticity model of polycrystalline specimen: representative volume determination | |
Dmitriev et al. | The features of structure transformation caused by nano-burnishing process | |
CN108062433A (zh) | 基于增材再制造点云模型的梯度曲面分层方法 | |
Barrett et al. | The Deformation Gradient of Interfacial Defects on Twin-like Interfaces | |
Bugas et al. | Network plasticity: reduced-order modeling of deformation-driven microstructure evolution |
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 |