纳米金属丝拉伸破坏及其应变率效应.pdfVIP

第27卷专辑 固体力学学报 VoL27S．Issue 2006年12月 ACTAMECHANlCASOLlDAS1N1CA December2006 纳米金属丝拉伸破坏及其应变率效应。 黄 丹 章 青 卓家寿 (河海大学土木工程学院工程力学系．南京，210098) 擒耍 建立了适于研究纳米金属快速变形破坏过程的分子动力学模型，并对不同应变率工况下不同截面 尺寸单晶镍纳米丝的零温单向拉伸破坏过程进行了分子动力学模拟．模拟得到各种纳米镍丝的应力-应变曲线、屈 服应变、屈服强度、断裂强度和初始弹性模量，提出了纳米金属丝快速变形力学性能的应变率效应预测公式并加以 验证．计算表明金属纳米丝的屈服应变与尺寸和应变率无关。屈服强度、断裂强度和弹性模量与应变率呈对数关 系． 关t调 纳米金属，力学性船，分子动力学。应变率效应 O 引言 形机制的影响． 宏观材料的强度会随着应变率的改变而变化， 1原理与方法 高速冲击载荷会导致材料的强度迅速增大[1]．在纳 本文采用经典分子动力学方法研究多原子系统 米尺度上直接从原子问作用来研究材料的力学性能 的运动过程，以面心立方晶格镍单晶为代表研究金 随应变率变化的规律和变形机制，既可以“从头开 属纳米丝的单向拉伸变形破坏过程及其应变率效 始”研究材料力学性能的应变率效应问题，也可以 应．镍晶体的原子间作用采用Voter等[7]根据修正 揭示纳米材料特定的力学行为和性能规律．Hu等嘲 的镶嵌原子法提出的镍的多体势描述，势函数的截 应用分子动力学方法模拟了不同的边界条件、晶向 断半径取为4倍最小原子距离，进行一次截断．系统 和变形速度情况下口．铁单晶的单向静态拉伸变形 过程，试图寻找决定纳米单晶体力学性能和强度的 行等温控制．纳米镍丝模型如图1，笛卡尔坐标系的 关键因素lIkedac3]等曾模拟了金属纳米丝中高应变 率导致的结构非晶化，BranicioE幻等研究了室温下将模型z，y方向设为自由表面，采用 纳米镍丝的拉伸强度随应变率增加而提高的现象以 及高应变率下的大变形问题；王秀喜等[5．．】给出了 面无约束，允许侧向自由独立变形．沿z向施加周期 纳米铜丝拉伸变形过程中应变率效应的敏感区、不 性边界条件，使模型呈一维纳米丝结构．本文按照图 敏感区和突变区以及各区对应的变形机制．但目前 1的几何构型共建立了4组尺寸不同的晶体模型，计 尚较少见从纳米尺度对材料力学性能的应变率效应 算元胞大小分别为(口o=0．352nnl，为镍的理想晶 建立量化公式． 格常数)： 本文建立了适于研究纳米金属快速变形破坏过 程的分子动力学模型和方法．基于镶嵌原子势描述 ；5ao)； 金属原子间作用，采用分子动力学方法研究了不同 尺寸金属纳米丝在不同加载速率下的单向拉伸变形 破坏过程，分析了应变率对金属纳米丝拉伸力学性 能的影响，提出了纳米丝屈服强度、断裂强度和初始 在求解多原子系统动力学方程的过程中，采用 弹性模量随应变率变化的预测公式，通过不同尺寸 Verlet算法的速度形式[1幻进行迭代计算，通过优化 模型加以验证，并讨论了不同应变率对纳米材料变 计算和稳定性测试本文将时间步长选取为4fs．在 和河海大学引进人才科研启动基金资助． 2005《*27# ～潞 日i*4d#*g《g 分子动力学模拟过程中，先进行温度初始化和长时 间初拍驰豫．然后在t001]方向施加拉伸应变，加载 过程中通过调整每次加载所经历的时间步长．共进