纯铁压缩的分子动力学模拟 - 图文-.doc

纯铁压缩的分子动力学模拟 摘要:本文介绍了利用lammps分子动力学模拟软件和VMD可视化处理软件模拟纯铁压缩过程的方法。得出了在压缩过程中,铁原子的排列和的运动规律。得到铁分子的运动细节,从而更加利于我们了解铁的性质。 关键词:纯铁压缩;lammps;VMD;分子动力学模拟; New Roman: Molecular dynamics simulation of the pure compression WangQiang,LiuChang (1. School of materials science and engineering, , Jiamusi University, Jiamusi 154007; Abstract:This paper introduces the simulation software and the method of VM visualization simulation software of pure compression process using lammps molecular dynamics. In the compression process, the arrangement of the iron atoms and the movement of the iron atoms are obtained. The motion details of the iron molecules are obtained, which help us to understand the properties of iron. Key words:Pure compression; lammps;VMD;molecular dynamics simulation; 1引言(四号,黑体 材料的力学性能随尺度发生变化尺度减小, 材料中缺陷存在的几率降低, 材料的强度提高同时尺度的变化可能导致材料内在变形竞争机制的改变, 例如多晶材料晶粒粒径在微米级以上时, 强度主要受位错强化机制控制, 而粒径进入纳米级后, 材料的变形主要来源于晶界滑移等机制原子尺度下, 微观效应占主导地位, 材料的理化、力学性能表现出与宏观不同、甚至相反的特性。Brenner发现金属单晶晶须拉伸强度与晶须直径呈反比,Fleck在微米级细铜丝的扭转试验中观察到尺寸效应纳米电机系统(NEMS的出现同迫切要求了解纳米尺度下材料的力学行为, 当前从实验上较难获得详细的信息, 而分子动力学模拟可以提供相关细节。 分子动力学模拟,是指对于原子核和电子所构成的多体系统,用计算机模拟原子核的运动过程,并从而计算系统的结构和性质,其中每一原子核被视为在全部其它原子核和电子所提供的经验势场作用下按牛顿定律运动。我们用c语言编写程序,VMD动画演示得到原子在压缩过程中的变化。在控制温度不变的情况下,得到了金属铁分子的动力学模拟过程。通过不断压缩,观察每次压缩过程中以及压缩后铁原子的排列,得到金属铁的运动细节,从而更加利于我们了解铁的性质。 2原理 材料的力学性能随尺度发生变化尺度减小, 材料中缺陷存在的几率降低, 材料的强度提高同时尺度的变化可能导致材料内在变形竞争机制的改变, 例如多晶材料晶粒粒径在微米级以上时, 强度主要受位错强化机制控制, 而粒径进入纳米级后, 材料的变形主要来源于晶界滑移等机制原子尺度下, 微观效应占主导地位, 材料的理化、力学性能表现出与宏观不 同、甚至相反的特性。Brenner 发现金属单晶晶须拉伸强度与晶须直径呈反比,Fleck 在微米级细铜丝的扭转试验中观察到尺寸效应纳米电机系统(NEMS的出现同迫切要求了解纳米尺度下材料的力学行为, 当前从实验上较难获得详细的信息, 而分子动力学模拟可以提供相关细节。 分子动力学模拟,是指对于原子核和电子所构成的多体系统,用计算机模拟原子核的运动过程,并从而计算系统的结构和性质,其中每一原子核被视为在全部其它原子核和电子所提供的经验势场作用下按牛顿定律运动。我们用c 语言编写程序,VMD 动画演示得到原子在压缩过程中的变化。在控制温度不变的情况下,得到了金属铁分子的动力学模拟过程。通过不断压缩,观察每次压缩过程中以及压缩后铁原子的排列,得到金属铁的运动细节,从而更加利于我们了解铁的性质。 分子动力学通过直接模拟原子的运动过程, 使我们能够详细了解模拟对象的演化发展历史分子动力学模拟的一个关键在于原子势函数的选取原子势早期一般采用简单的对势, 但对势无法正确描述弹性常数, 其结果不理想世纪年代提出的镶嵌原子法、有效介质理论更客观地反映了原子间多体作用的本质, 可得到较合理的结果.认为体系总能量为: 2/1020]12exp[]1exp[∑∑∑∑?? ??????????