金属钛拉伸的分子动力学研究.pdfVIP

· 4 · 材 料 开 发 与 应 用 2011年 6月 文章编号：1003—1545(2011)03-0004-04 金属钛拉伸的分子动力学研究 张 强 (南车四方车辆有限公司，山东 青岛 266331) 摘 要：本文利用分子动力学的研究方法，采用了钛的嵌入式原子势，建立了沿 [0001]晶向和 [0l10]方向拉 伸模型，结果表明：两种方向的拉伸均包含弹性变形阶段、屈服阶段、颈缩阶段、断裂阶段。沿 [0001]方向拉 伸时，滑移系少，取向偏离软取向。变形时屈服强度为3．55GPa，屈服应变为0．063，断裂时的应变达到0．55。 沿EOl10]方向拉伸时滑移系多，取向接近软取向。变形的屈服强度为2．4GPa，屈服应变为0．067，断裂时的 应变达到 1．2，沿[0l10]的拉伸变形更容易进行，表现出的塑性更好。 关键词：塑性变形；分子动力学；金属钛；拉伸 中图分类号：TG111．5 文献标识码：A 金属钛是二十一世纪最重要的结构材料之 在刚开始进行驰豫时，体系的势能上下波动较大， 【13 一 ， 由于其具有熔点高、硬度大、可塑性强、密 随着驰豫的进行，当驰豫时间到达 10ps时，其势能 度小、比强度高、耐高温、无毒、无磁性且抗腐 达到稳定，图2中总能量的变化与势能相似 ，刚开 蚀 等一系列优点而越来越广泛地应用于各领 始进行驰豫时，其能量波动较大，随着驰豫时间的 域。分子动力学方法是一种计算机模拟实验方 增长，到达 10ps时其总能量也达到平衡。 法 j，已经广泛应用于纳米材料的塑性变形研 究，本文采用分子动力学的方法研究金属钛拉伸 的塑性变形行为。 1 金属钛拉伸模型 1．1 嵌入式原子势 嵌入式原子势(EAM)是 Daw和Baskes|4根 据Hohen—kohn定理建立的势函数模型，具有运 算量小、运算效率高、适合金属间的作用等优点。 图1 驰豫时势能变化 嵌入式 (EAM)原子势的总能量 E可以表 示为： E=÷∑~~ij(rij)+∑Fi(pi) (1) 。 ji(jcj) i 其中，咖表示距离为 的原子 i和原子 之间 的对势能，Fi表示电子密度为 的原子 i嵌入能。 1．2 金属钛拉伸模型的建立 在模拟进行前，首先对 HCP型金属钛拉伸模 型进行驰豫，共驰豫 50ps，每次驰豫时间步长为 0．005ps，共驰豫 10000步，在模拟时温度保持在 300K。驰豫使体系达到平衡，从图1中可以看出， 图2 驰豫时总能量的变化 收稿 日期 ：2010—12—17 第26卷第3期 张 强：金属钛拉伸的分子动力学研究 ·5 · 由于金属钛存在各向异性，沿不同晶向拉伸 中可以看出，将金属钛的拉伸应力应变 曲线分为 其力学性能不同，本文 中建立两种拉伸模型，如 弹性变形阶段、屈服阶段、应力降低、断裂阶段。 表 1所示。 (1)弹性变形阶段 (a—b)：在此阶段，应力 表 1 沿两种不同晶向的拉伸模型 和应变呈正 比关系，随着应变的增大应力呈线性 增大，这一阶段的变形符合胡克定律，这时候原 子已经偏离的平衡位置 ，在此阶段只有原子的晶