CN105772689A - 基于分子动力学的铸造高铬合金的建模方法 - Google Patents
基于分子动力学的铸造高铬合金的建模方法 Download PDFInfo
- Publication number
- CN105772689A CN105772689A CN201610121445.4A CN201610121445A CN105772689A CN 105772689 A CN105772689 A CN 105772689A CN 201610121445 A CN201610121445 A CN 201610121445A CN 105772689 A CN105772689 A CN 105772689A
- Authority
- CN
- China
- Prior art keywords
- alloy
- chromium
- atom
- model
- unit
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22D—CASTING OF METALS; CASTING OF OTHER SUBSTANCES BY THE SAME PROCESSES OR DEVICES
- B22D21/00—Casting non-ferrous metals or metallic compounds so far as their metallurgical properties are of importance for the casting procedure; Selection of compositions therefor
- B22D21/02—Casting exceedingly oxidisable non-ferrous metals, e.g. in inert atmosphere
- B22D21/022—Casting heavy metals, with exceedingly high melting points, i.e. more than 1600 degrees C, e.g. W 3380 degrees C, Ta 3000 degrees C, Mo 2620 degrees C, Zr 1860 degrees C, Cr 1765 degrees C, V 1715 degrees C
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Heat Treatment Of Steel (AREA)
- Physical Vapour Deposition (AREA)
Abstract
本发明基于分子动力学的铸造高铬合金的建模方法属于合金材料结构设计领域,涉及一种基于分子动力学方法高铬合金的建模方法。建模方法利用分子动力学建模软件进行铸造高铬合金的组成结构的建模,先获取铸造高铬合金的特征参数,构造体心立方晶体结构的原子模型、构造面心立方晶体结构的原子模型;再通过置换程序,以铬原子置换出模型中的现有原子;然后判断置换元素的含量,如不满足含量要求,继续置换,直到铬元素满足合金中的要求;最后利用Voronoi算法对所构成的模型进行优化,直到满足合金中所有元素的含量要求。该建模方法能准确反映合金的微观组织成分和结构,对研究合金的性质和加工方法具有现实意义。
Description
技术领域
本发明属于合金材料结构设计领域,涉及一种基于分子动力学方法高铬合金的建模方法。
背景技术
随着航空航天、设备制造业、核工程等的迅猛发展,对材料的性能要求也越来越高,合金材料具有较高的强度、硬度,特别是耐高温、耐磨、耐腐蚀性能较纯金属来说具有相当大的优势,且合金材料能够按照要求改变材料的配比、生成条件得到不同性能的合金,所以越来越多的合金材料被应用于高精尖的产品中。高铬合金材料是核主泵推力轴瓦的常用材料,核主泵推力轴瓦由于长时间受到冲刷腐蚀的影响,故对材料的耐磨性和抗腐蚀性要求较高,且核工业本身对零部件的要求是要具有良好的表面完整性,以及更高的综合质量。分子动力学模拟(MolecularDynamicsSimulation)以其突出的方便于对所研究材料静态、动态结果分析的优势,应用于建立被加工材料、刀具等模型的建立,以及对切削加工过程的相互运动、应力、应变、加工过程的温度、产生的热以及温度场的分布、材料内部结构变化以及产生的晶格变形和缺陷等方面的研究,应用前景广泛,受到越来越多的研究者们的重视。
目前,国内对这种铸造高铬合金的研究主要集中在以实验为基础的合金性能及加工方法的研究上,成本花费较大、所需时间周期较长,并难以得出实验的微观作用机理对宏观性能的影响。所以需要利用分子动力学模拟仿真技术,通过计算机建立模型、算法,实现对材料微观结构与作用机理的研究。Jae-HyeokShim等发表的《AtomisticmodelingofnanosizedCrprecipitatecontributiontohardeninginanFe–Cralloy》一文在《JournalofNuclearMaterials》杂志的2009年第386-388卷,第56-59页,用刃型位错与纳米析出相之间的相互作用的方法研究了高Cr合金钢的α相的强化效应,没有合理考虑和建立Cr等元素置换形成的置换固溶体的合金模型,不能完整地体现合金的结构特征,为以后的研究设置了障碍。因此,仿真的分析结果难与实际吻合。
发明内容
本发明的目的在于解决对新型高铬合金自身性能以及加工过程中对机理研究的理论模型空缺的问题,提供一种基于分子动力学方法的铸造高铬合金的建模方法,建立这种铸造高铬合金的理论模型。保证合金模型组织元素位置分布的均匀性、准确性和建模模拟过程中的可靠性。
本发明采取如下技术方案是一种基于分子动力学方法的铸造高铬合金的建模方法,其特征在于,利用分子动力学建模软件进行铸造高铬合金的组成结构的建模,所述的铸造高铬合金的建模方法包括以下步骤:
步骤1,获取高铬合金的特征参数,包括合金中元素的种类、各种元素的体积分数、质量分数、各元素的凝固温度以及粒子间距和晶粒大小;
步骤2,基于高铬合金中组织结构的不同存在形式,其微观结构表现为不同的晶体结构;马氏体、铁素体等组织结构微观表现为体心立方的晶体结构,所以构造生长体心立方晶体结构,并对建好的晶粒进行旋转、平移晶粒中心的坐标至原始原点坐标位置,以方便后续的数据处理;
步骤3,合金中包含的奥氏体等组织结构微观表现为面心立方的结构,所以构造生长面心立方晶体结构,同样对生长的面心立方晶粒进行旋转和坐标平移处理,两次旋转的角度应该不同,以模拟更接近实际的晶体排列结构;
步骤4,统计铬的含量分布,编写置换程序,通过置换原子的程序,并对晶格常数的赋值得到铬原子,以铬原子置换模型中的现有原子,形成置换固溶体,置换过程中只置换原子参数不变的原子,已经置换过的原子不再置换;
步骤5,设定算法判据,若满足铬元素在合金中的体积分数和质量分数的要求,则结束原子置换过程,若未满足要求则重复步骤4的置换过程;
步骤6,利用Voronoi算法对以上步骤所构造的多晶模型进行优化处理,取一个多边形区域作为单元区域,对于这个单元区域,使超出单元左边界的模型平移到单元区域的右边,超出单元右边界的模型平移到单元的左边,超出单元上边界的模型平移到单元下边,超出单元下边界的模型平移到单元上边,使上下、左右构成一个循环,使其更符合多晶的晶体结构,以及合金中晶体的排列方式;
步骤7,通过赋值得到合金中其他元素的原子;重复步骤4到步骤6的过程,得到合金中其他所有元素的含量,最终得到铸造高铬合金的理论模型。
本发明的有益效果是本发明的构建高铬合金的微观结构模型的方法,能够准确生成真实合金中各元素的含量及空间排布方式,对研究该合金的自身性能和加工方法以及对加工中的切削力、温度场的分布、材料的结构变化都具有现实意义。
附图说明
图1:本发明建立铸造高铬合金模型方法的流程图;
图2:建立的合金多晶结构图;
图3:晶粒优化处理示意图。
具体实施方案
下面结合技术方案和附图的详细说明本发明的具体实施。
图1是本发明建立铸造高铬合金模型方法的流程图,首先获取铸造高铬合金模型的特征参数:元素的种类,各种元素的体积分数,质量分数,各元素的凝固温度以及粒子间距、晶粒大小和晶粒的分布情况。其中,元素种类包括钼(Mo)、铜(Cu)、镍(Ni)、铁(Fe)、锰(Mn)、铬(Cr)、钒(V)、硅(Si)、硫(S)。在进行分子动力学仿真过程中,工件和刀具模型是按照原子的组成规律来排列的。
在步骤2和步骤3中,由于铸造高铬合金中组织结构存在不同的形式,有马氏体、奥氏体、铁素体,这些不同的组织形式其微观结构表现为不同的晶体结构。其中马氏体、铁素体等组织结构微观表现为体心立方的晶体结构,奥氏体微观结构表现为面心立方的晶体结构,所以利用Fortran语言构造生长体心立方晶体结构,再构造面心立方的晶体结构,实际的合金结构中不同的晶粒具有不同的晶向,所以对建好的晶粒模型进行旋转来构造不同晶向的晶格结构。体心立方的晶格旋转α角,面心立方的晶格先旋转一个β角,再旋转α角。分子动力学中每个原子都是由空间坐标表示的,确定了坐标也就确定了原子位置,把晶粒形心的坐标平移至原始原点坐标位置,方便后续的数据处理。
高铬合金具有优良的耐磨性、硬度,由于合金中存在较多的铬碳化合物和铬固溶在奥氏体、马氏体及其转变的产物。铁原子的原子半径是126pm,铬原子的原子半径是128pm,可以满足置换固溶的条件,实际上铁与铬两元素是无限固溶的。所以在步骤4中,统计铬(Cr)的含量分布,编写置换程序,通过赋值得到铬原子,用铬原子随机置换模型中现有原子,置换过程中只置换原子参数不变的原子,置换过的原子就不再被置换。铬原子的原子半径比铁原子的原子半径大,置换后产生正晶格畸变,内能增高,阻碍位错滑移,使材料的强度、耐磨性增强。
设定算法判据,根据统计的铬元素在合金中的含量,编写程序判断是否满足要求,若满足则停止置换,若不满足则继续重复步骤4过程置换,直到满足要求为止。图2表示出建立的合金多晶结构图。
编写Voronoi算法程序,对步骤2到步骤5构建的模型进行优化处理。如图3所示,选取一个区域作为单元区域,在这个单元区域内随机选取n个点作为晶核,将这个单元区域扩展成9个相同的区域,对单元区域内的某个点与它相邻最近的其他点进行计算,并构造Voronoi多边形,由于这n个点是随机选取的,所以多边形的形状也是随机的,这是符合实际情况的。由于其他区域都是由单元区域扩展得到的,n个点的位置相同,对于处在单元区域外面的晶体部分,分割区域对其进行平移,即处在上边界面外的平移至下边界面内,下边界面外的平移至上边界面内,左边界面外的平移至右边界面内,有边界面外的平移至左边界面内,上下左右为一个循环,截取这个单元区域就是我们想要得到置换铬元素后的晶体微观结构模型,见图3。
按照步骤4到步骤6的过程对合金中其他元素:钼、铜、镍、钒、硅、硫置换模型中的铁元素,达到最后新型高铬合金的元素含量要求。建立铸造高铬合金的理论模型。
本发明提出的构建铸造高铬合金的微观结构模型的方法,能够准确生成真实合金中各元素的含量及空间排布方式,对研究该合金的自身性能和加工方法以及对加工中的切削力、温度场的分布、材料的结构变化都具有现实意义。
Claims (1)
1.一种基于分子动力学的铸造高铬合金的建模方法,其特征在于,利用分子动力学建模软件进行铸造高铬合金的组成结构的建模,所述的铸造高铬合金的建模方法包括以下步骤:
步骤1,获取新型高铬合金的特征参数,包括合金中元素的种类、各种元素的体积分数、质量分数、各元素的凝固温度以及粒子间距和晶粒大小;
步骤2,基于高铬合金中组织结构的不同存在形式,其微观结构表现为不同的晶体结构;马氏体、铁素体等组织结构微观表现为体心立方的晶体结构,所以构造生长体心立方晶体结构,并对建好的晶粒进行旋转、平移晶粒中心的坐标至原始原点坐标位置,以方便后续的数据处理;
步骤3,合金中包含的奥氏体等组织结构微观表现为面心立方的结构,所以构造生长面心立方晶体结构,同样对生长的面心立方晶粒进行旋转和坐标平移处理,两次旋转的角度应该不同以模拟更接近实际的晶体排列结构;
步骤4,通过置换原子的程序,并对晶格常数的赋值得到铬原子,以铬原子置换模型中的现有原子,形成置换固溶体,置换过程中只置换原子参数不变的原子,已经置换过的原子不再置换;
步骤5,设定算法判据,若满足铬元素在合金中的体积分数和质量分数的要求,则结束原子置换过程,若未满足要求则重复步骤4的置换过程;
步骤6,利用Voronoi算法对以上步骤所构造的多晶模型进行优化处理,取一个多边形区域作为单元区域,对于这个单元区域,使超出单元左边界的模型平移到单元区域的右边,超出单元右边界的模型平移到单元的左边,超出单元上边界的模型平移到单元下边,超出单元下边界的模型平移到单元上边,使上下、左右构成一个循环,使其更符合多晶的晶体结构,以及合金中晶体的排列方式;
步骤7,通过赋值得到合金中其他元素的原子;重复步骤4到步骤6的过程,得到合金中其他所有元素的含量;最终得到铸造高铬合金的理论模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610121445.4A CN105772689B (zh) | 2016-03-03 | 2016-03-03 | 基于分子动力学的铸造高铬合金的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610121445.4A CN105772689B (zh) | 2016-03-03 | 2016-03-03 | 基于分子动力学的铸造高铬合金的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105772689A true CN105772689A (zh) | 2016-07-20 |
CN105772689B CN105772689B (zh) | 2017-11-10 |
Family
ID=56387329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610121445.4A Expired - Fee Related CN105772689B (zh) | 2016-03-03 | 2016-03-03 | 基于分子动力学的铸造高铬合金的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105772689B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108427861A (zh) * | 2018-03-06 | 2018-08-21 | 中国工程物理研究院材料研究所 | 一种基于mpt工具包构建材料周期性多晶结构的方法 |
CN109545290A (zh) * | 2018-11-22 | 2019-03-29 | 南京航空航天大学 | 一种基于Voronoi分形技术的非晶合金自由体积检测方法 |
CN109858110A (zh) * | 2019-01-15 | 2019-06-07 | 燕山大学 | 基于分子动力学仿真的非晶合金中的缺陷表征方法 |
CN110008531A (zh) * | 2019-03-15 | 2019-07-12 | 北京应用物理与计算数学研究所 | 一种含间隙原子的无序固溶材料原子结构建模方法 |
CN112768009A (zh) * | 2020-12-25 | 2021-05-07 | 江苏方天电力技术有限公司 | 一种微尺度分析合金高温蒸汽腐蚀机理的方法 |
CN112992280A (zh) * | 2021-03-16 | 2021-06-18 | 苏州健雄职业技术学院 | 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 |
CN115372135A (zh) * | 2022-08-10 | 2022-11-22 | 国标(北京)检验认证有限公司 | 一种测定高温合金晶粒转动角度的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440361A (zh) * | 2013-07-19 | 2013-12-11 | 清华大学 | 一种等离子体刻蚀工艺中刻蚀产额的建模方法 |
CN104899393A (zh) * | 2015-06-19 | 2015-09-09 | 山东大学 | 一种离散相增强复合材料细观结构的建模方法 |
-
2016
- 2016-03-03 CN CN201610121445.4A patent/CN105772689B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440361A (zh) * | 2013-07-19 | 2013-12-11 | 清华大学 | 一种等离子体刻蚀工艺中刻蚀产额的建模方法 |
CN104899393A (zh) * | 2015-06-19 | 2015-09-09 | 山东大学 | 一种离散相增强复合材料细观结构的建模方法 |
Non-Patent Citations (2)
Title |
---|
JESB.SHERMAN ET AL.: "Crystalline Alloys of Organic Donors and Acceptors Based on TIPS-Pentacene", 《THE JOURNAL OF PHYSICAL CHEMISTRY》 * |
SHAOQING WANG.: "Atomic Structure Modeling of Multi-Principal-Element Alloys by the Principle of Maximum Entropy", 《ENTROPY》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108427861A (zh) * | 2018-03-06 | 2018-08-21 | 中国工程物理研究院材料研究所 | 一种基于mpt工具包构建材料周期性多晶结构的方法 |
CN108427861B (zh) * | 2018-03-06 | 2020-05-12 | 中国工程物理研究院材料研究所 | 一种基于matlab mpt工具包构建材料周期性多晶结构的方法 |
CN109545290A (zh) * | 2018-11-22 | 2019-03-29 | 南京航空航天大学 | 一种基于Voronoi分形技术的非晶合金自由体积检测方法 |
CN109858110A (zh) * | 2019-01-15 | 2019-06-07 | 燕山大学 | 基于分子动力学仿真的非晶合金中的缺陷表征方法 |
CN110008531A (zh) * | 2019-03-15 | 2019-07-12 | 北京应用物理与计算数学研究所 | 一种含间隙原子的无序固溶材料原子结构建模方法 |
CN112768009A (zh) * | 2020-12-25 | 2021-05-07 | 江苏方天电力技术有限公司 | 一种微尺度分析合金高温蒸汽腐蚀机理的方法 |
CN112768009B (zh) * | 2020-12-25 | 2023-09-22 | 江苏方天电力技术有限公司 | 一种微尺度分析合金高温蒸汽腐蚀机理的方法 |
CN112992280A (zh) * | 2021-03-16 | 2021-06-18 | 苏州健雄职业技术学院 | 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 |
CN112992280B (zh) * | 2021-03-16 | 2022-01-11 | 苏州健雄职业技术学院 | 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 |
CN115372135A (zh) * | 2022-08-10 | 2022-11-22 | 国标(北京)检验认证有限公司 | 一种测定高温合金晶粒转动角度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105772689B (zh) | 2017-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105772689A (zh) | 基于分子动力学的铸造高铬合金的建模方法 | |
Ding et al. | Effect of phase transition on micro-grinding-induced residual stress | |
Khoei | Computational plasticity in powder forming processes | |
Boër et al. | Process modelling of metal forming and thermomechanical treatment | |
Shin et al. | Predictive modeling capabilities from incident powder and laser to mechanical properties for laser directed energy deposition | |
CN102592035B (zh) | 一种车铣复合切削加工表面粗糙度及表面形貌仿真预测方法 | |
Chen et al. | A novel 3D surface generation model for micro milling based on homogeneous matrix transformation and dynamic regenerative effect | |
CN103995937B (zh) | 基于响应面和遗传算法的精密机床质量匹配设计方法 | |
Surmann et al. | Simulation of milling tool vibration trajectories along changing engagement conditions | |
Khajeian et al. | Shot peening effects on residual stresses redistribution of offshore wind monopile multi-pass weldments | |
Hsieh et al. | Particle swarm optimisation (PSO)-based tool path planning for 5-axis flank milling accelerated by graphics processing unit (GPU) | |
Molodov | Microstructural design of advanced engineering materials | |
CN105760578B (zh) | 一种含铰链间隙的连杆机构非概率时变可靠性综合方法 | |
Lv et al. | Application of phase-field modeling in solid-state phase transformation of steels | |
Zhang et al. | Filling path planning and polygon operations for wire arc additive manufacturing process | |
CN110543704A (zh) | 一种在反应堆堆芯流场计算中考虑局部结构影响的修正方法 | |
Fu et al. | Multiple-step incremental air-bending forming of high-strength sheet metal based on simulation analysis | |
CN105478765A (zh) | 一种基于金属3d打印球形粉末的紧密堆积的配粉方法 | |
Koirala et al. | Implementation of computer aided engineering for francis turbine development in Nepal | |
Hasan Al-Hafadhi et al. | Prediction and numerical simulation of residual stress in multi-pass pipe welds | |
Dadmal et al. | Finite element analysis of roller burnishing process | |
CN109210001B (zh) | 一种超塑成形空心扭转叶片的设计方法 | |
Dang et al. | A new incremental bending system to form curved metal plate by experiment and finite element simulation | |
BEELER et al. | Directional Dependence of the Displacement Energy Threshold for a FCC Metal | |
CN106372276A (zh) | 一种基于静刚度的机床结构件数据库建立方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171110 Termination date: 20200303 |