纳米晶铜单向拉伸变形的分子动力学模拟.pdfVIP

第 34卷 第 1期 力 学 学 报 Vo】_34．NO1 2002 年 1 月 A TA MECHANICA SINICA Jan．2002 纳米晶铜单向拉伸变形的分子动力学模拟 文玉华 避 刘日武 周承恩 中国科学院力学研究所1睦 性力学国家重点实验室，北京 100080 摘要 纳米材料是由尺度在 i~100nNl的微小颗粒组成的体系，由于它具有独特的性能而备受 关注．本文简要地回顾了分子动力学在纳米材料研究中的应用，并运用它模拟了平均晶粒尺寸 从 l 79～5 38llIIl的纳米晶件的力学性质 模拟结果显示；随着晶粒尺寸的碱小，系统与晶粒内 部的原子平均能量升高，而晶界上则有所下降 纳米晶件的弹性模量要小于普通多晶件，并随 着晶粒尺寸的减小而碱小；纳米晶铜的强度随着晶粒的减小而减小，显示了反常的Hall-Perch 效应；纳米晶体的塑性变形主要是通过晶界滑移与运动，以及晶粒的转动来实现的；位错运动 起着次要的、有限的作用；在较大的应变下 (约大于 5％)，位错运动开始起作用；这种作用随 着晶粒尺寸的增加而愈加明显． 关键词 纳米晶件，力学性质，分子动力学 纳米晶体材料由于具有与粗晶材料不同的力学性质而备受人们的关注．当晶粒尺寸减小到 纳米尺度时，它的硬度、强度和延展性都发生很大的变化 【 ．已有文献指出，尺寸效应可能改 变了材料的变形机制_41．常规晶体材料的塑性变形主要是通过位错的运动来实现的，界面阻碍 位错的运动从而阻碍了变形的进行，引起材料的强化．然而，当晶粒的尺度减小到纳米尺度，界 面的体积百分数急剧增加，它可能对纳米材料的变形起着积极的推动作用．目前，对纳米材料 的变形机制提出了几种机理，如空隙相的存在，大量的三叉界，位错运动的减少，晶界滑移， 晶界运动，以及晶粒旋转等等，均还没有取得一致的结论 【 ．尽管目前已有许多实验结果， 但对纳米晶体材料的弹性与塑性变形仍畦少清晰、直观的认识．这部分是由于重复实验条件和 控制关键的参数如密度和晶界结构的困难．所以目前对纳米材料的实验研究大多还是局限于对 它的弹性模量，硬度的测量，以及蠕变测试 [s~ioj． 计算机技术的飞速发展为材料科学的理论研究提供了良好的条件．一般来说，对材料结构 和性能的计算机模拟由两部分组成： (1)材料自身的模型化；它的结构要受两个因素制约，即 结构应尽量接近实验观察到的形态和受到计算机内存和计算时间的限制． (2)对实验现察到 的物理性质及有代表性的特征进行模拟计算．对材料的结构和性能的计算机模拟的计算模型有 两类：连续体模型和原子级模型．连续体模型就是把材料看作连续的介质，采用的方法大多是 有限元方法，它的研究对象是有限小的单元．原子级模型是把材料看作由许多单个原子的聚集 体，它的研究对象则是单个原子．材料的各个宏观量则由所有原子的统计量给出．目前，经常 采用的原子级模拟方法有分子动力学方法、蒙特卡罗方法和晶格动力学方法． 分子动力学计算机模拟作为研究复杂的凝聚态系统的有力工具，既能得到原子的运动轨 迹，还能像傲实验一样做各种观察，特别是许多与原子有关的微观细节，在实际实验中无法获 2000-08-28收到第一稿． 2ooi~7-06收到修改稿． 力 学 学 报 2002年 第 34卷 得，而在计算机模拟中可以方便地得到．这种优点使分子动力学在材料科学中显得非常有吸引 力．由于纳米材料的晶粒尺寸在纳米量绒，使得运用分子动力学计算机模拟纳米材料的性质成 为可能．由于分子动力学研究不受样品制备和测试技术的限制，因此分子动力学模拟能够找出 反映纳米晶体微观结构和力学性质之间的内禀本质，有助于人们对纳米材料的微观结构与性能 之间关系的理解．在纳米材料的分子动力学模拟中，不仅需要选择正确的势函数，而且要选择 合适的晶粒数和晶粒尺寸，即所要模拟的总原子数．模拟纳米材料的总原子数一般在 10 个以 上，考虑到计算量及结果的准确性，日前采用大都是多体势函数．有关分子动力学模拟在纳米 材料中的应用，可参见文献『l11．但是分子动力学模拟要受到空间尺度(模拟的