相场法模拟增材制造及焊接过程中显微组织的研究进展.docx

金属凝固过程是一个复杂的过程，涉及到高温、组织相变以及熔体与基体材料之间的相互影响。随着计算机技术及数值模型的快速发展，通过数值模拟方法研究增材制造以及焊接熔池的凝固过程成为可能。

近年来，学者们通过数值模拟方法积极探索凝固过程显微组织的演变规律，以实现对材料（零件）力学性能和物理性能的预测，获取工艺调控凝固组织的理论依据，并建立工艺参数与组织演变的关系。

目前，对凝固过程中显微组织进行数值模拟的常用方法有确定性方法、蒙特卡洛法、元胞自动机法和相场法。

增材制造（AM）是一种利用计算机辅助设计逐层堆积材料的零件成形技术，具有周期短、可成形复杂结构零件、力学性能优异等特点，广泛用于航空航天、汽车船舶、武器装备等领域高端装备的制造。增材制造过程中熔池的凝固行为影响诸如溶质偏析、裂纹、气孔等缺陷的形成，同时也会影响熔池组织的尺寸和形态，最终决定零件的性能。

通过传统试验方法能够获得工艺参数对熔池组织、气孔、裂纹等的影响规律，实现优化工艺、改善构件质量的目的。然而，大量的试验不仅耗时耗力，而且由于增材制造过程中熔池体积小、凝固速率快，采用试验方法难以对熔池内部凝固过程进行观测，无法获得完整凝固组织的形成过程，只能解决宏观层面的问题。

有效控制显微组织演变，进而提高材料性能是增材制造技术的发展方向之一。同样，焊接熔池的凝固过程类似于激光增材制造熔池的凝固过程。

常用的显微组织模拟方法

确定性方法

是以晶核和生长物理模型为起点，利用确定的微观运动方程描述系统状态，以经典凝固动力学理论为基础的一种方法。

确定性模型指在给定时刻和一定体积内，晶粒的形核密度和生长情况均为确定的函数，但这些函数需要通过试验才能得到。

由此可知，确定性方法具有真实的材料组织演变的物理基础，并且使用该方法能够构建出符合系统物理本质的数学模型且能够精确预测特定系统的组织特征。

元胞自动机法

即CA法，是在20世纪40年代后期提出的，这种数值算法常用于处理和描述复杂系统在空间和时间上的演化规律。

CA法的基本思想如下：将凝固过程所需时间分解为相同的时间周期（时间步），将所需计算的凝固区域分解为一定尺寸的元胞（计算单元），每个元胞与其相邻元胞相互作用，并赋予每个元胞温度以及溶质浓度等物理量；定义一套演化规则，将被赋予元胞的物理量作为变量并且放之整个计算域皆成立，其中的每个元胞可转变为多种状态；在达到一定条件后，状态之间可以互相转换，在单个时间步内，以演变规则为基础对元胞各种物理量进行计算，再判断元胞状态，每个时间步完成一次循环，最终在一段时间内实现物理过程的模拟。

蒙特卡罗法

即MC法，是随机统计数学中的一种方法，最早应用于复杂的数学问题。蒙特卡罗法基于某种抽样或者试验得到对特定事件的期望值。

有关学者将其引入到组织演变仿真中，其原因是MC法具有良好的物理基础。

MC法的基本思路如下：将模拟区域（这里只表明空间）离散为按照规律进行分布的格点并且对每个格点赋予一个随机数Q，Q所代表的物理意义是这个格点的晶粒生长取向；如果相邻格点间的随机数Q相等，则表示两个格点属于同一个晶粒，反之则属于不同的晶粒，并形成界面；界面存在一定的界面能，在界面能的驱动下界面实现迁移，即进行晶粒生长过程中的界面迁移。

MC模型只适用于固态晶粒的生长，例如再结晶和第二相沉淀。

相场法

即PF法，以Ginzburg-Landau理论为基础，通过微分方程反映相变过程扩散、有序化势及热力学驱动力的综合作用，是模拟相变过程组织演变的理想方法。

在扩散界面处，相场参数φ连续地从-1变化到1；相对应地，自由能密度函数f(φ)呈现出先增大后减小的趋势。当实时温度T与材料熔点Tm相同时，固相和液相的自由能密度相等；当TTm时，固相的自由能密度小于液相。根据实际物理过程，可以建立其对应的自由能密度函数。相场法能够在界面平衡的情况下描述新相从母相中生长的复杂过程以及凝固过程中固/液界面的演变情况。

相场法最大的优势在于避免了复杂的固/液界面追踪，通过外部场与相场的结合使微观与宏观模型相结合，从而可直接模拟宏观相场作用下的显微组织演变，并且可以定量地研究热力学扰动、组织各向异性以及不同材料参数对显微组织演变的影响。相场法在模拟熔池显微组织生长过程方面具有广阔的应用前景。

优缺点比较

下表列出了确定性方法、MC法、CA法、相场法的优缺点。可知相对于MC法，相场法的物理意义更加明确，并且与CA法相比，其模型的合理性不需要真实物理现象来决定；确定性方法的模拟结果与试验结果吻合较差。因此，相场法是模拟凝固过程中显微组织演变的理想方法。

方法

优点

缺点

相场法

相场与多种物理场耦合能够准确地模拟凝固过程中晶粒形貌的变化，还可以模拟材料凝固过程中固/液界面的形态和运动情况

部分参数的确定方法受到质疑；计算区域有限制，计算效率较低并且