纳米铜丝尺寸效应分子动力学模拟.pdfVIP

维普资讯 第 34 卷 第 2期 力 学 学 报 Vo1．34,No．2 2002 年 3 月 ACTA MECHANICA SINICA M脚l h．姗 纳米铜丝尺寸效应的分子动力学模拟 梁海弋 王秀喜 吴恒安 王 字 (中国科学技术大学力学和机槭工程系，合I巴230026) 摘要 采用EAM 原子势函数对不同截面尺寸纳米铜丝拉伸性能进行了零温分子动力学模拟， 研究表 明截面变化对纳米丝拉仲性能有明显影响．由于表面原子弛豫，纳米丝存在表面张应 力以及与之相平衡的内部应力．这种本征应力的存在是纳米丝尺寸效应 的根源．纳米丝截面减 小，刚拉仲强度提高、屈服推迟、韧始拉仲模量的软化程度增加． 关键词 纳米铜丝。尺寸效应，拉伸性能，分子动力学 引 言 人们很早就知道材料的力学性能随尺度发生变化 LlJ．尺度减小，材料中缺陷存在的几率降 低 ，材料 的强度提高；同时尺度 的变化可能导致材料 内在变形竞争机制 的改变，例如多晶材料 晶粒粒径在微米级 以上时，强度主要受位错强化机制控制，而粒径进入纳米级后，材料的变形 主要来源于晶界滑移等机制 原子尺度下，微观效应 占主导地位，材料的理化、力学性能表现 出与宏观不同、甚至相反的特性． BrennerI3]发现金属单晶晶须拉伸强度与晶须直径呈反 比． Fleck[J‘在微米级细铜丝的扭转试验中观察到尺寸效应．纳米电机系统 (NEMS)的出现 【迫切 要求了解纳米尺度下材料的力学行为，当前从实验上较难获得详细的信息，而分子动力学模拟 可以提供相关细节． 分子动力学通过直接模拟原子的运动过程，使我们能够详细了解模拟对象 的演化发展历 史．分子动力学模拟的一个关键在于原子势函数的选取．原子势早期一般采用简单的对势，但 对势无法正确描述弹性常数，其结果不理想． 20世纪 80年代提出的镶嵌原子法 (EAM)[~I、 Finnis—Sinclair方法 、有效介质理论 (EMT)s[]以及 glue势 等原子势更客观地反映了原子 间多体作用 的本质，可得到较合理的结果． EAM等多体势主要考察了原子在局域背景电子云 密度环境下 的结合能，认为体系总能量可分解为 【 Eto~=∑ )+ ∑ ) ，(q) 其 中 为第 个原子处 的平均 电子云密度，函数 F反映了原子 {与背景电子云的相互作用， 函数 采用对势形式， r为原子 ，J间的距离． EAM 等原子势虽然比传统对势复杂，但计 29Ol-O2-o7收到第一蔷 2~}Ol-09-04收翻修改善． 1)国家教委博士点基金 资助疆 目． 维普资讯 第 2 期 梁海弋等 ：纳米铜丝尺寸效应 的分子动 力学模拟 209 算量并没有明显增加．目前采用 EAM 等多体势 的分子动力学在微机上可模拟数十万原子．在 并行机上则已达到千万原子量级． 周国辉等 _1l_模拟 了液体金属吸附促进裂尖位蜡发射的机制．丁家强等 【模拟了纳拳铁 的热力学性能． LiShen等 模拟了压结条件下 晶体铜中微裂纹的愈合过程． Baskes等 u 计算了纳米镍单晶剪切强度随尺寸变化 的规律． 本文采用分子动力学模拟了不同截面尺寸方形纳米铜丝 的拉伸过程 (模拟原子数从 1014 变化到 96ooo)，试图揭示尺寸效应对纳米丝拉伸特性的影响及其内在机制．为统一起见，除特 别注明外，本文长度单位采用铜的晶格常数 (0．361rim)，应力单位采用 GPa． 1 分子动力学模拟过程 纳米铜丝初始构型按几何方法生成，横截面 (XY面)为正方形，原子位置按理想点阵排 列．如图1(a)，X， Z坐标轴分别对应面心立方的 1【0o]，