晶体相场模型模拟纳米孪晶材料.pdfVIP

晶体相场模型模拟纳米孪晶材料 摘要：本文以纳米孪晶材料为研究对象，在等体积变形条件下，采用晶体相场模型，采用等体积 方法，利用Matlab 工具，模拟了纳米孪晶材料的凝固形成，以及在单向拉伸下的形变过程，通过 对模拟出的微观原子排列进行分析，讨论了纳米孪晶材料微观变形机制。 关键词：纳米孪晶 晶体相场模型 单向拉伸 微观形变 引言 纳米孪晶(1—100nm)材料由于具有与传统粗晶材料 （平均晶粒尺寸 [1] d1mm）不同的力学性质而备受关注 。近10 年来，低能态共格孪晶界由于其 [2] 明显优于传统非共格界面的特殊强韧化效应而受到材料学界的广泛关注 。实验发 [3] 现 ，微米或亚微米尺度的孪晶片层强化效果并不显著，但当孪晶片层厚度薄至纳 米量级时其强化效果就显著增加。位错－共格孪晶交互作用使得纳米孪晶材料具有 非常独特的结构－力学性能关系。纳米尺度共格孪晶界面不仅能够有效阻碍位错的 运动，表现出类似于传统晶界强化的效果；同时，孪晶界附近可提供丰富的位错存 [4] [5] 储空间(不同于晶界强化) ，从而保证材料足够的加工硬化与稳定塑性变形能力 。 事实上，孪晶的微观结构特征(如孪晶片层厚度、孪晶长度、孪晶界取向、晶体取向 等)以及外在加载方式、受力状态等均会直接影响变形过程中可启动位错类型，从而 显著影响位错与孪晶界的反应结果。因此了解孪晶的形成与微观形变机制是必要的。 但是，由于纳米尺度的实验存在着很大的困难，实验上制备纳米孪晶结构不易控制， 难以获得理想的纳米孪晶层结构，目前尚无法通过实验详细地了解纳米材料的力学 性能，因此本文用晶体相场法 （Phase Field Crystal，PFC），使用Matlab 工 具对其进行了相关模拟。 1晶体相场法研究理论 1.1晶体相场模型的特点 传统的相场模型是用空间均衡的场来表征平衡态，即分别用具有不同恒定值的 序参量表示固、液两相，忽略了晶体相中固有周期对称性结构特性。因此，传统相 场模型最大的局限性在于不能体现晶体点阵的任何信息，从而不能从根本上反映出 源于原子尺度上晶格交互作用所引起的弹、塑性变形机制和相转变动力学。 要建立一种基于微观尺度的数值模拟方法以预测其宏观性质。 [7] 与传统相场模型不同，2002 年，Elder 等 人提出了基于密度泛函理论的PFC 模型，该模型可以看作是相场模型的一个更加深化的扩展形式，其构造的平衡自由 能产生了周期对称性的原子密度状态，而不是空间上均衡一致的状态。这个自由能 函数自洽地包含弹性能、各向异性和对称特性 （如弹性常数间的关系、位错的类型 和数量、小角度晶界能以及点阵位置等），把弹性和结晶学自然地结合到了一起，所 以说PFC 模型可以反映晶体由对称性所决定的所有性质。通过耗散和保守性的弛豫 动力学方程，PFC模型可以描述从原子振动到热/质传输的不同时间尺度上的原子密 PFC 度场的演变，揭示原子效应的影响。同时，通过给定合适的初始条件， 模型可 以得到有多个不同取向晶粒的多晶系统。由此可以看到PFC 方法是相场方法和分子 动力学方法的桥梁，可以方便地应用于纳米孪晶系统的研究。 1.2 晶体相场模型理论计算原理 2007 年，Elder 等人直接利用经典密度泛函理论推导出PFC 模型的简化形式及 动力学方程，由此可以建立不同类型的PFC 模型以适应不同的研究体系。 1.2.1晶体相场模型的自由能函数 [32] 根据对经典密度泛函理论 表达式进行合适的近似 ： F      