3. CALPHAD热力学模型与数据库

CALPHAD方法的核心，是把复杂多元合金的吉布斯自由能写成可计算、可参数化的数学表达式，再用实验相图、热化学数据和第一性原理结果去拟合这些表达式中的参数。模型决定“用什么函数描述相”，数据库决定“这些函数里的参数取什么值”。理解这两层关系，是后续用Thermo-Calc、Pandat或JMatPro进行高温合金设计的前提。

从理想溶液到亚正规溶液：自由能模型的递进

最基础的模型是理想溶液。它假设不同原子混合时无热效应，混合焓为零，混合熵只由组态熵决定：

G_m^ideal = Σ x_i G_i^0 + RT Σ x_i ln x_i

其中 x_i 为组元 i 的摩尔分数，G_i^0 为该组元在参考态下的摩尔吉布斯自由能。理想溶液能解释“混合使自由能下降”这一基本趋势，但对高温合金中Ni-Al、Ni-Cr、Ni-Ti等强相互作用体系，误差很大，因为它完全忽略了原子间的键合差异。

正规溶液在理想溶液基础上加入一个与成分相关的超额自由能项，且假设混合焓与温度无关、混合熵仍为理想组态熵：

G_m^ex = x_A x_B L_AB

L_AB 称为相互作用参数，单位为J/mol。若 L_AB < 0，表示A-B键比A-A、B-B键更稳定，倾向于形成有序或化合物；若 L_AB > 0，则倾向于偏聚或分层。正规溶液已能描述许多二元体系的液相和固溶体，但对非对称的相图（如液相线明显不对称）仍显不足。

亚正规溶液进一步让相互作用参数随成分变化，最常用的形式就是Redlich-Kister多项式：

G_m^ex = x_A x_B Σ_{v=0}^{n} L_AB^{(v)} (x_A - x_B)^v

其中 v 为多项式阶数，L_AB^{(v)} 为第 v 阶相互作用参数，通常写成 A + BT 的形式以引入温度依赖。Redlich-Kister展开的物理意义是：用 (x_A - x_B) 的幂次去修正非对称性，阶数越高，能拟合的相图细节越丰富。对高温合金常用的Ni基、Co基体系，液相和γ相通常取到 v=0~2 即可获得较好精度；对存在强短程有序的体系，可能需要更高阶或改用其他模型。

化合物能量模型与有序-无序模型

高温合金中大量存在金属间化合物相，如γ′-Ni₃Al（L1₂结构）、γ″-Ni₃Nb（D0₂₂结构）、η-Ni₃Ti、δ-Ni₃Nb等。这些相不能用简单的置换固溶体模型描述，需要化合物能量模型（Compound Energy Formalism, CEF）。CEF把每个相看成由若干亚点阵组成，每个亚点阵可容纳不同组元，相的自由能由各亚点阵上的组元占据分数和相互作用参数共同决定：

G_m = Σ_I Π_I y_i^I G_I^0 + RT Σ_I a_I Σ_i y_i^I ln y_i^I + G_m^ex

其中 I 为亚点阵编号，a_I 为亚点阵位数，y_i^I 为组元 i 在亚点阵 I 上的占据分数，G_I^0 为端元化合物能量。以γ′-Ni₃Al为例，通常采用两亚点阵模型 (Ni,Al)₃(Ni,Al)₁，四个端元分别为Ni₃Ni、Ni₃Al、Al₃Ni、Al₃Al，其中Ni₃Al是实际有序结构对应的端元。通过拟合有序-无序转变温度、有序度及相平衡数据，可以同时描述γ和γ′两相。

有序-无序模型则用于处理同一相在不同温度下从无序固溶体转变为有序结构的情形。典型做法是把有序相和无序相统一为一个模型，用长程有序参数η描述有序度，当η=0时退化为无序固溶体，η≠0时出现超点阵。对高温合金而言，γ/γ′共格组织正是有序-无序模型与CEF结合描述的重点：γ为无序FCC，γ′为L1₂有序相，两者在成分和结构上连续过渡，模型必须能同时给出两相的自由能和相平衡。

高温合金常用热力学数据库

模型是骨架，数据库是血肉。目前高温合金CALPHAD计算最常用的商业数据库包括：

• TCNI系列（Thermo-Calc）：覆盖Ni基高温合金的γ、γ′、γ″、δ、η、σ、μ、碳化物、硼化物等相，是学术和工业界使用最广的Ni基数据库之一。TCNI8及以上版本对Co、Cr、Mo、W、Ta、Ti、Al等组元的多元相互作用参数做了较系统的评估。
• TTNI系列（ThermoTech）：同样面向Ni基高温合金，特点是包含较多Co基和Ni-Co-Cr体系的评估，适合做Co基或Ni-Co基合金设计。
• PanNi系列（Pandat）：与PanPhaseDiagram配合使用，界面友好，适合快速计算平衡相分数、相图和凝固路径。
• JMatPro：内置Ni基、Fe基、Co基数据库，偏向工程应用，可直接输出力学性能、物理性能随温度的变化，适合做合金筛选和工艺窗口评估。

使用数据库时需要注意三点：第一，数据库有适用成分和温度范围，超出范围外推可能失真；第二，不同数据库对同一体系的参数评估来源不同，计算结果可能有差异，关键设计应交叉验证；第三，数据库中的相名称和亚点阵模型要与自己的合金体系对应，不能直接套用。

从模型到计算：一个Ni-Al二元示例

以Ni-Al二元系为例，用亚正规溶液模型描述液相和γ相，用CEF描述γ′-Ni₃Al。液相RK参数可写为 L_liq = −80000 + 20T J/mol（示例值，仅用于说明形式），γ相 L_γ = −50000 + 10T J/mol，γ′采用两亚点阵 (Ni,Al)₃(Ni,Al)₁，端元能量 G_Ni₃Al^0 设为 −40000 J/mol（示例值）。将这些参数输入Thermo-Calc，计算得到的Ni-Al相图应能再现γ′有序化温度约1385℃、共晶温度约1338℃等关键特征。若计算值与实验偏差较大，优先调整γ′端元能量和γ/γ′相互作用参数，而不是盲目增加RK阶数。

实际高温合金设计通常从二元、三元体系出发，逐步加入Cr、Co、Mo、W、Ta、Ti等组元，每加入一个组元都要检查数据库是否覆盖该组元与已有组元的相互作用。对多元体系，CALPHAD计算的核心输出包括：平衡相分数随温度的变化、γ′固溶温度、TCP相（σ、μ、P）析出倾向、凝固路径与偏析行为。这些结果直接支撑合金成分优化和热处理制度制定。

掌握模型与数据库的关系后，下一步就是学会用Thermo-Calc或Pandat实际建立计算、读取相图、判断结果合理性。模型选择是否正确、数据库是否适用，最终都要通过计算与实验的对比来检验。