虽然MEMS发展已经有二十多年的时间了,但粘着问题一直阻碍MEMS产品进入市场,是MEMS技术急需解决的关键问题。分子动力学模拟作为一种数字实验方法,可以对理论分析和实验观察都难以了解的现象做出一定的解释。分子动力学模拟是研究粘着问题的一种有效方法。针对MEMS中的粘着问题,本文选取Au材料,以分子动力学模拟为研究手段,对MEMS球面模型接触过程进行了模拟,综合考虑引力和斥力,并对模拟结果进行了分析。本文首先介绍了分子动力学模拟技术的基本原理,然后建立了分子动力学的仿真模型,并确定了分子动力学模拟的各种仿真参数。模拟中采用EAM势函数描述原子间的相互作用,采用Velocity-Verlet算法求解系统运动方程。本文研究了球面模型接触过程中原子构型的变化;讨论了在接触过程各个阶段中接触区域的压力分布情况;分析了接触半径随着接触力以及压头位移的变化关系;分析了接触力随着压头位移的变化关系以及接触力与接触区域原子层数的关系;分析了接触过程中系统总能量随时间的变化关系以及系统势能随着压头位移的变化关系。本文使用C++语言进行微/纳观接触问题的分子动力学研究,代码精简易懂,可移植性好,为分子动力学模拟技术的编程风格进行了改进。研究表明:压头逐渐趋近于基体表面时,基体表面在粘着力的作用下以“突跳”的方式与压头迅速接触,得到“突跳”现象;当压头反向撤离时,与压头相连的基体材料发生颈缩,得到颈缩分离现象;分离粘着过程中的接触力相对于粘着接触过程中的接触力呈现明显的滞后,这表明粘着力在分离粘着过程中对弹性基体的原子构型有较大影响;接触力的大小只与接触区域的上面少数几层原子有关;受压状态下,接触区域近处原子受压力,远处原子受拉力;接触过程中出现了晶体变形和原子迁移现象。