纳米金属材料单向拉伸应力应变关系的数值模拟研究.pdf

首届 FEPG 用户年会论文集 纳米金属材料单向拉伸应力应变关系的数 值模拟研究 李林安 邸玉贤 李鸿琦 （天津大学机械学院力学系 天津 300072） 摘要：本文利用 “复合有限元”方法，将分子动力学方法引入到连续体的有限元方法中，用 Morse 势函数 推导出纳米金属固体材料 （如α−Fe ，Cu，Al 等）晶界区域的本构关系，建立由晶粒、晶界、孔洞等组 成的复合模型。借助有限元源程序自动生成系统 fepg 软件，调整晶粒尺寸，研究在单向线性逐步加载过程 中，晶粒尺寸对材料非线性应力应变关系的影响，以及材料宏观弹性模量在加载过程中的变化情况。得出 结论，随着晶粒尺寸的减小，材料的非线性特征逐渐增强，材料的延伸性也随之增强。 关键词：纳米金属，分子动力学，有限元，非线性，应力应变关系 1． 引言 纳米材料又称超微细材料，是由微小颗粒一一绝大多数是晶体，其特征尺度至少在一个 方向上为纳米量级——组成的固体，其典型的晶粒尺度为 l～100nm。晶粒大小是影响传统 金属多晶材料(晶粒尺寸在微米以上量级)力学性能的重要因素。当晶粒尺寸减小到纳米尺度 时，它的硬度、强度和延展性都发生很大的变化[1]。通过大量的实验测试、计算模拟及理论 分析，证明金属纳米材料具有非常独特的力学性能及结构—性能关系。研究变形机理的最大 困难在于难以从实验获得可靠的本征力学性能参量。目前纳米材料样品制备技术上存在的问 题以及纳米材料本身的结构稳定性，使准确测量纳米材料本征力学性能数据十分困难。 计算机技术的飞速发展为材料科学的理论研究提供了良好的条件。。对材料的结构和性 能的计算机模拟的计算模型有两类：连续体模型和原子级模型。连续钵模型就是把材料看作 连续的介质，采用的方法大多是有限元方法，它的研究对象是有限小的单元。原子级模型是 把材料看作由许多单个原子的聚集体，它的研究对象则是单个原子。材料的各个宏观量则由 所有原子的统计量给出。目前，经常采用的原子级模拟方法有分子动力学方法、蒙特卡罗方 法和晶格动力学方法。 材料的应力-应变关系反映了该材料的所受应力随应变的变化关系。它在材料的性质研 究中起着重要的作用。不同的材料，一般来讲，它的应力-应变曲线特征是有所不同的。晶 粒大小是影响晶体材料力学性能的重要因素。本文提出一种 “复合有限元”思想，将分子 动力学方法引入到传统的有限元方法中，模拟研究晶粒尺寸对纳米金属材料在单向拉伸过 程中非线性应力应变关系的影响。 2.计算模拟方法 在晶粒内部原子严格排布在晶格格点上，有关实验及研究表明，纳米材料在受宏观载荷 的情况下，在晶粒内部几乎不发生位错。在接近晶粒表面，原子开始偏移晶格，发生畸变， 晶界区域原子无规则排列，材料破坏也主要发生在这里。可以假想纳米材料是由许多晶粒堆 62 首届 FEPG 用户年会论文集 积成的整体，晶粒和晶粒之间由晶界材料填充。这样，就可以把纳米材料看成是由晶界基体 和晶粒增强相组成的复合材料，当然可以考虑三叉晶界和孔洞等。由于晶粒内部一般不发生 破坏，可以认为晶粒是由线弹性均质材料组成。可以借助分子动力学原理，用原子势函数找 出晶粒之间力的相互作用关系，在晶界区域填充晶界材料，使这种晶界材料的本构关系符合 晶粒之间的相互作用关系，这样就建立起了连续体模型。 [2] 在本文中，我们采用如下的 Morse 势函数 ： φ(r ) D[exp{−2α(r −r )} −2 exp{−α(r −r )}] （1） ij ij 0 ij 0 式中， D 为原子结合能，α 为与原子有关的常数，r 为原子 i 与 j 的原子间距，r 为 ij 0 两原子间的平衡间距，对势函数求导可以得到原子间作用力函数。 将纳米晶体材料的细观结