4. CALPHAD方法原理与多元相图计算流程

CALPHAD（Calculation of Phase Diagrams）的核心思想并不神秘：它把“相图”看作热力学模型在给定温度、压力和成分条件下的稳定相平衡计算结果，而模型中的可调参数则来自实验与第一性原理数据的联合评估。对多元体系而言，CALPHAD 的真正价值在于“用低元数据外推高元行为”，因此它既是一套参数优化方法，也是一套从二元、三元到多元的逐级外推流程。本节以 Al-Zn-Mg-Cu 四元系为例，把一次完整的相图计算任务拆解为可执行步骤，并说明每一步的输入、操作与判断依据。

一、CALPHAD 的基本原理与模型层次

CALPHAD 用吉布斯自由能函数描述每个相。对置换固溶体相（如 fcc、液相），常用 Redlich-Kister 多项式表达过剩自由能：

G_m = Σ x_i G_i^0 + RT Σ x_i ln x_i + Σ_i Σ_j x_i x_j Σ_k L_ij^k (x_i - x_j)^k

其中 G_i^0 是纯组元参考态自由能，L_ij^k 是二元交互作用参数，可含温度线性项 a + bT。对三元及以上，还需加入三元交互参数 L_ijk。对化学计量化合物（如 Al₂Cu、MgZn₂），则用形成焓与熵描述。CALPHAD 参数优化的本质，就是让这些 L 参数在实验相平衡数据、热化学数据和第一性原理数据之间取得最优折中。

二、以 Al-Zn-Mg-Cu 系为例的任务分解

Al-Zn-Mg-Cu 是典型的高强铝合金体系，其相图计算涉及 fcc-Al 基体、MgZn₂、Al₂CuMg、Al₂Cu、T 相（AlZnMgCu 三元化合物）等多个相。一次完整的 CALPHAD 计算任务可分解为以下六个阶段。

阶段 1：实验数据评估与筛选。 首先收集 Al-Zn、Al-Mg、Al-Cu、Zn-Mg、Zn-Cu、Mg-Cu 六个二元系，以及 Al-Zn-Mg、Al-Zn-Cu、Al-Mg-Cu、Zn-Mg-Cu 四个三元系的相平衡数据。数据来源包括：DTA 测定的液相线温度、EPMA 测定的相成分、XRD 物相鉴定、以及量热法测定的形成焓。评估时需检查数据自洽性，例如同一合金的液相线温度与相成分是否对应同一平衡状态。对 Al-Zn-Mg 三元系，重点筛选 T 相与 MgZn₂ 的相区边界数据；对 Al-Cu-Mg 三元系，重点筛选 Al₂CuMg 的溶解度线。示例中可假设从文献中筛选出 120 组二元数据、85 组三元数据，并标注数据质量等级（A/B/C）。

阶段 2：热力学参数优化。 采用 Thermo-Calc 的 PARROT 模块或 Pandat 的 PanOptimizer，以最小二乘为目标函数，对每个二元系独立优化。以 Al-Zn 系为例，液相与 fcc 相的交互参数 L_AlZn^0 和 L_AlZn^1 需同时拟合液相线、固相线和活度数据。优化时先固定纯组元数据，再逐步释放交互参数。判断收敛的标准是：计算相图与实验数据的偏差在液相线温度上小于 5 K，在相成分上小于 2 at.%。若偏差过大，需检查是否遗漏了某个中间相或数据点存在系统误差。

阶段 3：三元系外推与验证。 将优化好的二元参数直接外推至三元系，加入三元交互参数 L_AlZnMg 等。外推后与三元实验数据对比，若 T 相相区计算值与实验偏差超过 10 K，则需调整三元参数。这一步的关键是“先外推、后修正”，而不是一开始就拟合三元数据，否则会破坏二元系的物理合理性。

阶段 4：四元系外推与计算。 将三个三元系的参数外推至 Al-Zn-Mg-Cu 四元系。此时通常不再引入四元交互参数，除非有明确的四元化合物或实验证据。计算时设定温度范围 300–700 °C，成分范围 Al 70–95 wt.%，Zn 0–10 wt.%，Mg 0–5 wt.%，Cu 0–5 wt.%。使用 Thermo-Calc 的 POLY 模块或 Pandat 的 Phase Diagram 模块，计算等温截面、垂直截面和液相面投影图。

阶段 5：结果判读与实验验证。 对计算得到的 475 °C 等温截面，观察 fcc-Al 相区、MgZn₂ 相区、Al₂CuMg 相区的位置与实验金相照片是否一致。若发现某相区计算值偏大，可回溯至三元参数调整。验证时建议选取 3–5 个关键合金成分，用 DTA 测定其液相线温度，与计算值对比，偏差应控制在 5 K 以内。

阶段 6：数据库封装与复用。 将优化后的参数按 TDB 格式封装，包含所有二元、三元参数及纯组元数据。TDB 文件需注明版本号、优化日期、数据来源和适用范围。例如：

ELEMENT AL FCC_A1 26.9815 4577.3 28.3215 !
ELEMENT ZN HCP_ZN 65.38 5657.0 41.631 !
PARAMETER L(AL,ZN;0) 298.15 10000-5*T; 6000 !

封装后的数据库可直接用于后续多元合金设计，如 7075 铝合金的凝固路径计算。

三、计算流程中的关键判断与常见问题

在参数优化阶段，最常见的问题是“过拟合”：为了拟合某个实验点而引入高阶交互参数，导致外推至多元系时出现异常相区。判断方法是检查参数的温度依赖性是否合理，例如 L 参数不应随温度剧烈变化。另一个问题是“数据冲突”：不同文献的液相线温度相差 20 K 以上，此时应优先采用实验方法更可靠、成分标定更准确的数据，并在 TDB 中标注数据来源。

在四元外推阶段，需注意“三元交互参数的传递性”。若 Al-Zn-Mg 三元系优化时引入了 L_AlZnMg，该参数在四元系中仍然有效，但需检查其是否与 Al-Zn-Cu 三元系的参数产生矛盾。若矛盾，应回到三元系重新优化，而不是在四元系中强行修正。

通过上述六个阶段，CALPHAD 方法将分散的实验数据转化为可计算的热力学数据库，再通过逐级外推获得多元相图。对 Al-Zn-Mg-Cu 系而言，这一流程不仅能给出平衡相图，还能用于计算凝固路径、析出相分数和热处理窗口，是多元合金设计的基础工具。实际执行时，建议先完成二元系的严格优化，再逐步外推，并在每一步保留完整的参数版本和验证记录，以便后续追溯和修正。