ThermalCalc计算相图：CALPHAD数据库驱动下的多元相平衡与Scheil凝固路径

ThermalCalc计算相图——准确说是Thermo-Calc（用户常把Thermo-Calc简称为ThermalCalc）——是CALPHAD社区应用最广的相图计算软件。CALPHAD方法的核心思想是把各相的Gibbs自由能写成温度、成分和压力的函数，通过数据库中的模型参数逐相描述，再用自由能最小化求解给定条件下的平衡相组成。听起来很系统，但计算出的相图是否可靠，几乎完全取决于后端数据库的质量。

数据库选择的代价

Thermo-Calc支持多种热力学数据库：TCFE（钢和铁合金）、TCNI（镍基高温合金）、TCTI（钛合金）、SSOL（焊料合金）、TCHEA（高熵合金）等。每个数据库对同一元素和相的Gibbs自由能描述可能不一致——因为模型参数的优化是在特定成分范围和相的集合内进行的。

团队在Inconel 718合金的相图计算中遇到过数据库选择问题。TCNI10（镍基合金数据库v10）和TCFE10（铁基数据库v10）都包含Inconel 718的所有主元素（Ni, Cr, Fe, Nb, Mo, Ti, Al），但两个数据库对γ’相（Ni₃(Al,Ti)）的自由能描述有明显差异——TCNI10预测的γ’溶线温度约960°C，TCFE10预测约1010°C，差了50°C。差距的来源是TCNI10对γ/γ’两相区做过专门优化（镍基合金的核心需求），而TCFE10的重点在铁基奥氏体/铁素体平衡。对于镍基高温合金，选TCNI是正确的——领域专用数据库在它的设定成分范围内精度更高。

亚稳相和平衡相：凝固路径的关键分歧

平衡相图告诉你在给定温度下热力学上最稳定的相组成。但实际凝固过程（尤其是冷却速率>1 K/s）远离平衡——Scheil凝固模型假设固相内扩散可以忽略、液相成分完全均匀，给出的是非平衡凝固路径。

在Inconel 718的Scheil计算中，Laves相（(Fe,Ni)₂(Nb,Mo)六方结构）在平衡相图中到约850°C就消失了（转化为δ相Ni₃Nb），但在Scheil路径中Laves相残留到固相线以下，最后在约650°C被冻结在凝固组织中。Laves相的残留对合金的力学性能有直接危害——它消耗了本应形成γ”强化相的Nb。这个现象是平衡计算完全看不到的。

这就引出了一个实用规则：ThermalCalc计算相图时，平衡模型用于判断相稳定性和固溶热处理窗口，Scheil模型用于预测实际的凝固微观组织和元素偏析。两个模型互补使用，不能只做一个。

计算参数对结果的敏感性

相图计算中的3个技术参数需要仔细设定。第一：温度步长在相边界附近要细化。在大块液相或单相区步长可以设20-50°C，但在液相线、固相线和固相转变附近，步长应缩到2-5°C，否则可能跳过窄的相区。团队在TCFE10计算AISI 316L不锈钢的相图时，初算用了10°C步长——漏掉了M₂₃C₆碳化物在720-740°C之间仅20°C宽的析出窗口。改为2°C步长后这个窄相区才浮现出来。

第二：成分空间的网格密度。对于等温截面图（isothermal section），三元系通常用吉布斯三角形上的50×50网格做相分布——这个精度可以发现宽度≥2 at.%的相区，但对于<1 at.%宽度的窄相（如金属间化合物线），50×50可能完全漏掉。需要做局部加密——在可疑区域重新跑200×200甚至更高的网格。

第三：初始平衡起始点的选择。Thermo-Calc用全局自由能最小化算法，但在复杂体系中可能有多个局域极小值。做等温截面时从退火温度（如1000°C的单相区）平衡开始，逐步降温，每个温度步用上一温度的平衡结果作为初始猜测——这条”温度路径”有助于避开局域极小和虚假相转变。

计算结果的实验验证

ThermalCalc计算相图的结果必须用实验数据做标定。DSC（差示扫描量热）和高温XRD是两种互补的验证手段。DSC给出相转变温度——Inconel 718的γ’溶线温度DSC测定为958±5°C，TCNI10计算结果960°C，偏差2°C。但DSC只能告诉你”有相转变”，说不了”什么相变什么”——这是高温XRD或原位中子衍射的领地。

对于工程合金设计中的相图计算问题，核心策略是多数据库交叉校验+Scheil补充非平衡路径+关键节点实验验证。更多关于CALPHAD方法和多元合金计算的案例，站上有详细复盘。

更多内容请访问 https://www.keyanxueshu.com/