金属固液间相变热力学及微结构分子动力学模拟-材料加工工程专业论文.docxVIP

山东大学博士学位论文一金属固液问相变热力学及徽结构分子动力学模拟 山东大学博士学位论文一金属固液问相变热力学及徽结构分子动力学模拟：=============；=；=======；======2=# 摘 要 本文采用常温常压分子动力学模拟的方法，利用EAM、TB作用 势，对金属的熔化结晶现象、原子团簇的熔化过程、液态金属的中程 有序结构、液固结构相关性及熔体热力学性质等进行了研究，并采用x 射线衍射方法对模拟结果进行实验验证，以确保模拟结果的可靠性。 甜Cu加热一冷却过程的模拟结果分析发现，由EAM作用势计算得 到的cu的熔化温度为1520K，结晶温度为1010K，而cu的熔点为 1356K，可见加热和冷却过程存在很大的温度滞后现象。液态cu在冷 却过程中，当温度高于1130K时，冷却过程对V．T关系曲线无影响， 即在该温度以上，不论前面的过程、初始条件以及冷却速率如何，V—T 关系是可以再现的。说明体系处于一种平衡的液态。当温度低于1130K 以后，v．T关系曲线与冷却速率有关，模拟中观察到当冷速为4 x 10t3K]s时，形成非晶，非晶形成温度约800K；而冷速为4×10¨K／s 时，形成晶体，结晶温度为1010K，结晶温度与冷却速度有很大关系， 冷却速度越低，则结晶温度就越高。斗， 对无自由表面的大块固态Ni和有自由表面的Ni原子团簇的熔化 过程进行了分子动力学模拟实验。(模拟结果表明，无自由表面的大块 固态Ni，其熔化温度要高于实验测得的Ni的熔点；而有自由表面的 Ni原子团簇，其熔化温度要低于Ni的熔点。原子团簇的熔化温度与团 簇中原子的个数有关，原子个数越多，其对应熔点越高。分析认为存 在两种不同的熔化机理：在无自由表面的条件下，熔化过程受均质形 核机理控制，液相形核需较大的过热度，使得熔化温度远高于Ni的熔 点；对团簇的熔化，由于自由表面有较高的能量，随温度升高，液相 核心首先在表面形成，充当内层的非自发形核核心，表层首先形成一 液态层，这样整个团簇的熔化由表层向中心区域推进，使得熔化温度 低于Ni的熔点，即熔化过程受非均质形核机理的控制。斗 摘 摘 要 建立了晶液两层构型方法的晶体生长分子趁应堂搓赵程序。／并用 该程序描述了Au及Ni50A150的晶体生长过程，发现在这种构型]淦属 及合金的生长都按粗糙界面的生长方式进行，液相原子连续、垂直地 向界面添加，使整个品面迅速向液相中推移，结晶很快结束。 在对液态Ni20A180进行x射线衍射和MD模拟研究中发现，液态 Ni：。A180在结构因子上存在明显的预峰，标志着体系存在一定程度的中 程有序结构。实验和模拟所得到的结构因子曲线在误差允许范围内很 好地吻合，模拟所得到的配位数和原子间距与Maret的结果在误差范 围内也很一致，说明模拟结果的可靠性。通过对MD模拟结果的分析 发现，液态体系中存在着大量的Frank多面体结构，它是由一定数量的 1551、1441和1661键对组成。预峰所对应的结构可描述为Frank多面 体通过共有原子方式连接而成的原子团簇，且每个团簇的边缘原子最 多只能成为两个Frank多面体的共有原子，其尺寸介于5．20A之问。共 对液态Ni50A150合金分别在两种不同的冷却条件下原子团簇的演 变过程、相变热力学过程进行了模拟。(两种不同的冷却速率得到两种 不同的固态组织：晶态和非晶态，团簇的种类和数量在不同的冷却条 件下也发生了不同的变化。在4 x 1013K／s的冷速下，偶分布函数的第 一峰发生了劈裂，表征体系非晶的形成；在4×10¨K／s的冷速下，偶 分布函数在低温下明显显示出晶态峰的特征，表明体系形成了晶体。 液态金属中存在由键对构成的多面体团簇。在4×10BK／s的冷速下， 标征体系非晶结构的五重对称键型155l和1541所形成的小团簇数量 在非晶体系中占主导位置，同时还存在少量的1431、1441和1661等 其它键型形成的小团簇，这样，体系中含有大量的二十面体团簇结构； 在4×1011K／s的冷速下，标征面心立方结构的1421键型所形成的小团 簇数量在晶体中占第一位，同时还存在一定量的1422、1441和1661 等其它类型的键对所形成的团簇，体系最终形成面心立方结构。一方 面，液态金属中的某些微观结构信息被保留到固态当中，另一方面液 态和固态在微观结构上存在着较大的差异，体现出液态金属冷却过程 中的结构相关性和变异性。液态Ni，oAl50在4 x 1013K／s的冷速下，其体 ：：‘=======口=====；================ 山东大学博士学位论文一金属固液间相变热力学及徽结构分子动力学模拟 积、能量与温度呈线性关系；当冷却速度为4×lo“列s时，其体积、 能量与温度的曲线在某一温度范围出现突变，该突变温度即结晶