冷速对液态合金Ca50Zn50快速凝固过程中微观结构演变的影响

采用分子动力学方法对六种不同冷速对原子尺寸相差较大的液态合金Ca50Znso凝固过程中微观结构演变的影响进行了模拟研究,并采用双体分布函数、Honeycutt—Andersen（HA）键型指数法、原子团类型指数法（CTIM一2）、可视化等方法进行了深入分析,结果表明：系统存在一个临界冷速,介于和5×10^11K／s与1×10^11K／s之间,在临界冷速以上（如1×10^11K／s,1×10^11K／s,1×10^11K／s和5×10^11K／s）时,系统形成以1551,1541,1431键型或二十面体基本原子团（12012000）等为主体的非晶态结构;在临界冷速以下时,系统形成以1441和1661键型或bcc基本原子团（1460800）为主体（含有少量的hcp（1200066）和fcc（12000120）基本原子团）的部分晶态结构．在非晶形成的冷速范围内,其总双体分布函数的第一峰明显分裂成与近邻分别为Zn—Zn,Ca—Zn,Ca—Ca相对应的三个次峰;且随着冷速的下降,同类原子近邻的次峰峰值升高、异类原子近邻的次峰峰值下降;Zn原子容易偏聚,随着冷速降低,二十面体的数量增多,非晶态结构也越稳定．在晶态形成的冷速范围内,Zn原子己大量偏聚形成大块bcc晶态结构,Ca原子也部分形成hcp和fcc晶态结构．     

冷速对液态金属Pb凝固过程中微观团簇结构演变影响的模拟研究

采用分子动力学方法对六种不同冷却速率对液态金属Pb凝固过程中微观团簇结构演变的影响进行了模拟跟踪研究.采用双体分布函数、Honeycutt-Andersen键型指数法、原子团类型指数法(CTIM-2)、平均配位数等方法对凝固过程中微观团簇结构的演变进行了分析.结果表明：系统存在一个形成非晶态或晶态结构的临界冷速(介于5×10¹²与1×10¹²之间),大于这个临界冷速时,系统将形成以1551,1541和1431键型为主的非晶态结构;当小于这个临界冷速时,系统将先形成以1441和1661键型或以bcc基本原子团(14 6 0 8 0 0)为主的晶态结构,并稳定存在一段时间,然后又迅速转变为以1421和1422键型为主或以fcc基本原子团(12 0 0 0 12 0)和hcp基本原子团(12 0 0 0 6 6)以一定比例并存的部分晶态结构;同时发现,冷速对系统中fcc结构和hcp结构的相对比例有明显的影响,冷速越低,fcc结构所占的比例越多,越倾向于形成完美的fcc晶态结构. 关键词： 液态金属铅 凝固过程 团簇结构演变 分子动力学模拟     

冷却速率对液态金属Cu凝固过程中微观结构演变影响的模拟研究

采用Quantum Sutton-Chen(Q-SC)多体势对液态金属Cu在四个不同冷却速率下的凝固过程进行了分子动力学模拟研究. 通过双体分布函数、键型指数、配位数、均方位移及可视化分析, 结果表明：冷却速率对液态金属Cu的微观结构演变有决定性影响. 当冷却速率为1.0×10¹⁴K/s时得到非晶态结构；当冷速分别为1.0×10¹³K/s，1.0×10¹²K/s和1.3×10¹¹K/s时，系统形成以1421键型为主体的面心立方(fcc)与六角密集(hcp)共存的混合晶体结构；且其结晶温度分别为373K，773K和873K，即冷速越慢，其结晶温度越高，结晶程度也越高；且冷速越慢，1421键型越多，混合晶体中面心立方(fcc)结构所占的比例越高. 同时发现，原子的平均配位数的变化与1551，1441，1661键型的变化密切相关, 反映出体系对称性结构的变化规律与配位数的变化有关. 在可视化分析中，进一步采用中心原子法展现出非晶态与晶体结构的2D截面，及在3D下混合晶体中两个基本原子团分别为面心立方(fcc)与六角密集(hcp)基本原子团的具体结构. 关键词： Q-SC多体势 液态金属Cu 凝固过程 分子动力学模拟     

Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟

本文用分子动力学模拟研究了Ni熔体以不同冷速凝固后微观结构的演变规律, 并通过理论计算确定出了Ni熔体凝固后获得理想非晶的临界条件. 模拟结果发现冷速小于10¹¹ K/s时, Ni 熔体凝固后形成晶态组织; 冷速在10¹¹ K/s到10^14.5 K/s之间时, Ni熔体凝固后形成由晶态结构与非晶态结构组成的混合组织. 冷速小于10¹⁰ K/s, Ni 熔体凝固后形成的晶态组织具有fcc结构; 冷速在10¹⁰ K/s到10^14.5 K/s之间时, Ni熔体凝固后组织中的晶态由fcc和hcp结构层状镶嵌排列构成. 通过分析模拟结果和计算结果, 确定出了Ni熔体凝固后形成理想非晶的临界冷速为10^14.5 K/s. 并发现Ni熔体中临界晶核(冷速等于10^14.5 K/s)和亚临界晶核(冷速大于10^14.5 K/s) 均由fcc和hcp组成层状偏聚结构, 这表明Ni熔体中生长的晶体、临界晶核和晶胚的结构是相同的. 关键词： 分子动力学模拟 晶体团簇 临界冷速 结构     

液态金属Ga急冷凝固中微观结构转变的模拟研究

采用分子动力学方法对液态金属Ga的快速凝固过程进行计算机模拟跟踪研究,运用HA键型指数法和原子团类型指数法分析了金属原子Ga的成键类型和形成的基本原子团结构.发现:与通常的液态金属(如Al)相反,随着温度的降低,二十面体及与二十面体相关的1551, 1541键型数目明显下降;与立方体(fcc,bcc),六角密集结构相关的1421,1422键型数目明显增加;而与菱面体相关的1321,1311键型的数目却显著增加,逐渐占据优势.最后形成一种新型的以菱面体结构为主、夹杂着立方体(fcc,bcc)、六角密集(hcp)等团簇结构所组成的非晶态结构.这正是非晶态金属Ga的g(r)曲线分裂的第二峰的顺序为前低后高,而与非晶态金属Al的g(r)曲线(其分裂的第二峰的顺序为前高后低)明显不同的微观结构上的物理根源.     

冷速对液态金属Na凝固过程中微观结构影响的模拟研究

采用分子动力学方法对液态金属Na在四种不同冷速下的快速凝固过程进行了模拟跟踪研究.采用双体分布函数g(r)曲线、Honeycutt-Andersen键型指数法和原子团类型指数法对凝固过程中微观结构的变化进行了分析.结果表明：冷却速率对微结构的转变有决定性影响，当冷速为1.0×10¹⁴和1.0×10¹³K/s时，系统形成以1551和1541键型或以缺陷多面体基本原子团(13 1 10 2)和二十面体基本原子团(12 0 12 0)为主体的非晶态结构；当冷速为1.0×10¹²和1.0×10¹¹K/s时，系统则形成以1441和1661键型或以体心立方基本原子团(14 6 0 8)为主体的晶态结构.同时发现：不同冷速对液态金属Na在液态和过冷态时微观结构的影响甚小；但不同冷速对其固态(非晶态利晶态)时的微观结构有显著的影响，且要在液-固转变点(分别在玻璃转变温度T_g和晶化起始温度T_c)附近或以后才能充分展现出来.根据这一特点，有可能建立另一种确定液态金属T_g和T_c的新方法.原子团类型指数法比键型指数法更有利于研究液态、非晶态等无序体系和一些晶化体系的具体结构特征. 关键词： 液态金属Na 凝固过程 分子动力学模拟 原子团类型指数法     

不同冷速对金属铜Cu凝固过程中微观团簇结构演变影响的模拟研究

采用分子动力学方法和Quantum Sutton-Chen(Q-SC)多体势对由5万个液态金属铜(Cu)原子构成的系统在三个不同冷却速率下的凝固过程中微观团簇结构转变的影响进行了模拟研究.运用双体分布函数、Honeycutt-Andersen(HA)键型指数法、原子团类型指数法(CTIM-2)和可视化分析等方法,对凝固过程中微观团簇结构的演变特性进行了分析研究.结果发现:由非晶体向晶体转变的临界速度约为1.0×1013 K/s,在此冷速下系统形成非晶体和晶体以一定比例并存的混合结构;在冷速为1.0×1014 K/s冷却时系统形成以1551、1541、1532、1431键型为主的非晶体结构,非晶转化温度约为673K;在以4.0×1012 K/s速度冷却时,系统从673K就开始结晶,并形成以1421和1422二种键型或由这二种键型构成的面心立方(FCC)(12 0 0 0 12 0)和六角立方(HCP)(12 0 0 0 6 6)基本原子团为主的晶体结构,尤其是由1421键型构成的面心立方(12 0 0 0 12 0)基本原子团在晶体生长和微观团簇结构形成过程中占主导地位.同时发现,冷速对金属Cu系统中的FCC结构和HCP结构的相对比例有显著的影响,冷速越低,FCC基本原子团以及由其构成的团簇结构越多.     

分子动力学模拟纯镁熔体凝固过程中结构演化

采用分子动力学模拟方法模拟了在周期性边界条件下由500个原子构成的液态Mg模型系统的凝固过程,分别考察了在5×10^14 K/s、5×10^13 K/s、1×10^13 K/s 、1×10^12 K/s的冷却速率下液态Mg熔体的凝固行为。模拟结果很好地重现了实验值。模拟中原子间作用势采用FS势,结构分析采用径向分布函数、均方位移、系统总能量分析、H-A键对分析技术等方法。结果表明,当冷却速率为5×10^14 K/s时,系统形成以1541键对为主的非晶态结构;当冷却速率分别为5×10^13 K/s、1×10^13 K/s、1×10^12 K/s时,系统形成以1421、1422键对为主的hcp晶态结构;另外,在快速冷却形成非晶的过程中,大部分bcc结构被保留下来,而在慢冷形成晶态的过程中,大部分bcc结构最终演化形成了hcp结构。     

冷速对Cu10Ag90合金凝固结构的影响

通过分子动力学对液态Cu10Ag90合金在四种冷速条件下进行快速凝固模拟。结果显示,1 × 1010和1 × 1011 K/s下系统的平均原子能量分别在750 K和650 K发生突变,冷速越低最终平均原子能量越低;1 × 1012和1 × 1013 K/s的双体分布函数第二峰出现分裂,表明结构处于非晶态。从1 × 1011K/s开始出现尖锐小峰,表明此冷速开始出现晶化现象,1 × 1010 K/s下分裂的峰更加尖锐明显,说明体系形成结晶度较高的晶体结构;1 × 1010 和1 × 1011 K/s下系统凝固后晶体结构含量由高到低分别为fcc, hcp, bcc。冷速越低晶体结构数目越多,系统的有序度更高,结构熵越低。     

液态Mg7Zn3合金快速凝固过程中微观结构演变机理的模拟研究

采用分子动力学方法对液态Mg₇Zn₃合金的快速凝固过程进行了模拟研究,并采用双体分布函数、Honeycutt-Andersen键型指数法、原子团类型指数法(CTIM)以及遗传跟踪等方法对凝固过程中团簇结构的形成、演变特性进行了分析.结果表明:在以冷速为1×10¹² K ·s^-1的凝固过程中, Zn-Zn原子之间的相互成键的概率明显增加,形成以1551,1541,1431键型为主的非晶态结构.其特征键型1551随温度变化 关键词： 液态MgZn合金 快速凝固过程 微观结构演变 分子动力学模拟     

A simulation study of microstructure evolution during solidification process of liquid metal Ni

A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method (EAM) model potential. The bond-type index method of Honeycutt--Andersen (HA) and a new cluster-type index method (CTIM-2) have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution: below the crystallization temperature $T_{\rm c}$, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of $2.0\times10^{13}$\,K\,$\cdot$\,s$^{-1}$ and $1.0\times10^{12}$\,K\,$\cdot$\,s$^{-1}$, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp (expressed by (12 0 0 0 6 6) in CTIM-2) and the fcc (12 0 0 0 12 0) basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc (12 0 0 0 12 0) basic clusters consisting of 1421 bond-type, and their crystallization temperatures $T_{\rm c}$ would be 1073 and 1173\,K, respectively.     

Microstructural evolution and martensitic transformation mechanisms during solidification processes of liquid metal Pb

The microstructural evolution and the martensitic transformation (bcc–hcp and bcc–fcc) mechanisms during the solidification process of liquid metal Pb were studied by molecular dynamics simulation. Results indicate that, with the decrease of temperature, the system undergoes two phase transitions: from the liquid state into a metastable bcc phase first and then from the bcc phase into a coexisting crystal structure of hcp and fcc phases. Moreover, the complicated martensitic transformation processes are clearly observed by cluster type index method (CTIM) and the tracing method. The two transformation mechanisms are very analogous at the atomic level; the essential difference between them is that, in the bcc–hcp transformation, two adjacent layers shift in opposite directions, whereas in the bcc–fcc transformation, the top layer and bottom layer shift in opposite directions relative to the middle layer. The specific mechanisms for the bcc–hcp and bcc–fcc transformations are confirmed to correspond to the revised Burgers mechanism and Bain mechanism, respectively.     

液态金属In凝固过程中微观结构转变的模拟研究

采用分子动力学方法对液态金属In的快速凝固过程进行计算机模拟跟踪研究.运用HA键型指 数法和原子成团类型指数法分析了金属原子In的成键类型和形成的原子团簇结构.发现：与通常的液态金属(如Al)相反，随着温度的降低，二十面体及与二十面体相关的1 551 键越来越少；与四方体，六角密堆积相关的1421，1422和1431键数目总和变化很小；而与菱 面体相关的132l，1311，1301和1201的数目却随着温度的降低而显著增加，逐渐占据优势 .最后形成一种新型的以菱面体结构为主、夹杂着立方体(fcc，bcc) 关键词： 液态金属In 微结构转变 团簇结构 分子动力学 计算机模拟     

Structural evolution in the crystallization of rapid cooling silver melt

The structural evolution in a rapid cooling process of silver melt has been investigated at different scales by adopting several analysis methods. The results testify Ostwald’s rule of stages and Frank conjecture upon icosahedron with many specific details. In particular, the cluster-scale analysis by a recent developed method called LSCA (the Largest Standard Cluster Analysis) clarified the complex structural evolution occurred in crystallization: different kinds of local clusters (such as ico-like (ico is the abbreviation of icosahedron), ico-bcc like (bcc, body-centred cubic), bcc, bcc-like structures) in turn have their maximal numbers as temperature decreases. And in a rather wide temperature range the icosahedral short-range order (ISRO) demonstrates a saturated stage (where the amount of ico-like structures keeps stable) that breeds metastable bcc clusters. As the precursor of crystallization, after reaching the maximal number bcc clusters finally decrease, resulting in the final solid being a mixture mainly composed of fcc/hcp (face-centred cubic and hexagonal-closed packed) clusters and to a less degree, bcc clusters. This detailed geometric picture for crystallization of liquid metal is believed to be useful to improve the fundamental understanding of liquid–solid phase transition.     

单轴应变驱动铁bcc—hcp相转变的微观模拟

用分子动力学方法模拟了沿〈001〉晶向应变加载和卸载情况下单晶铁中体心立方(bcc)与六方密排(hcp)结构的相互转变,分析了相变的可逆性和微结构演化特征.微观应力的变化显示样品具有超弹性性质,而温度变化表明在相变和逆相变过程中均出现放热现象.相变起始于爆发式均匀形核,晶核由块状颗粒迅速生长为沿{011}晶面的片状分层结构; 而卸载逆相变则从形核开始就呈现片状形态,且相界面晶面指数与加载相变完全一致,表现出形态记忆效应.在两hcp晶核生长的交界面易形成面心立方(fcc)堆垛层错. fcc通过在hcp晶粒内     

冷速对液态金属Ga凝固过程中微观结构演变影响的模拟研究

采用分子动力学方法对液态金属Ga凝固过程中不同冷却速率对微观结构演变的影响进行了模拟跟踪研究. 运用HA键型指数法和原子成团类型指数法(CTIM)分析了金属原子Ga的成键类型和形成的基本原子团结构. 结果发现，冷却速率对凝固过程中的微观结构起着非常重要的作用. 在以1.0×10¹⁴，1.0×10¹³，1.0×10¹²K/s的速率冷却时，系统形成以与1311，1301键型相关的菱面体结构为主，夹杂着立方体、六角密集等其他团簇结构所构成的非晶态结构；在以1.0×10¹¹K/s的速率冷却时，系统明显发生结晶，其结晶转变温度T_c约为198K，且冷却速率越慢，结晶转变温度T_c越高，形成以与1421键型相关的斜方晶体(经可视化分析确认)为主要构型的晶态结构. 这将为研究液态金属的结晶转变过程提供一种新方法. 关键词： 液态金属Ga 凝固过程 微结构转变 分子动力学模拟