三元相图计算

三元合金的凝固路径和相组成，靠实验逐点测定极其耗时——一个 A-B-C 三元系，成分空间是二维三角，实验点动辄上百。我这个项目客户要做一种铸造铝合金，想用 CALPHAD 方法做 三元相图计算，把液相线等温线和三元共晶点算出来，直接指导配料和工艺。

CALPHAD 的原理是先有” thermodynamic database”——把各相的 Gibbs 自由能表达成温度、成分的函数，参数来自实验和热力学优化。计算时软件（Thermo-Calc / Pandat）在给定成分和温度下，令各相化学势相等，求解平衡相组成。我用的数据库包含了 A、B、C 三个组元及其二元、三元交互作用参数。

第一步是验证数据库质量。这是最容易被跳过的致命环节。如果某些二元或三元交互参数缺失，软件会外推或报错，算出来的共晶点可能偏几十度。我的做法是先用二元子系统（A-B、B-C、A-C）分别计算，和文献实验相图对比，确认三个二元系都对，再扩展到三元。我这个项目的 B-C 二元系初版参数偏旧，共晶温度差了 15℃，换新版数据库后吻合。

第二步是算液相面投影。我在整个成分三角内扫描，计算每一点的液相线温度（液相开始凝固的温度），用 Gibbs 三角表示：颜色代表液相线温度（高温红→低温黄），等温线把三角分成若干区域，标出 L+α、L+β、α+β+γ 等相区。这张图直接告诉客户”某个配料会在什么温度开始凝固、先析出什么相”。

最关键的信息是三元共晶点：成分约在 A:B:C≈50:30:20 附近、温度约 256℃，是 L→α+β+γ 三相同时凝固的最低点。客户据此把配料定在共晶点附近，因为共晶成分流动性最好、热裂倾向最低，最适合铸造。

第三步是 Scheil 非平衡凝固。上面算的是平衡凝固，实际铸造冷速快、溶质来不及扩散，是非平衡过程。我用 Scheil 模型（假设固相内无扩散、液相完全混合）算了非平衡凝固路径，发现实际凝固区间比平衡宽了约 30℃，且会出现平衡相图里没有的偏析组织。这部分给客户的工艺意义是：实际铸造要控制冷速，否则偏离相图预测。

判读逻辑：液相面投影是”配料→凝固起点”的图，共晶点是”最优铸造成分”的锚，Scheil 曲线是”实际会偏多少”的修正。三者合起来才能指导真实工艺。

踩坑记录：第一，数据库版本。不同版本数据库对三元系参数差异大，必须交叉验证，我用了两套数据库取交集确认。第二，步长。液相面扫描步长太粗会漏掉窄相区，我用细网格重扫关键成分区。第三，亚稳相。Scheil 非平衡下可能出现亚稳相，软件默认只算平衡相，需手动开启亚稳选项。第四，外推陷阱。三元交互参数缺失时软件外推结果看起来”合理”但可能是错的，必须用二元验证兜住。

从我的工程经验看，三元相图计算 最怕”拿来数据库就点计算”。我给客户的建议是：先验证二元子系统、再用 Scheil 补非平衡、最后才定配料。更多相图计算，见 [三元相图计算](https://www.keyanxueshu.com/category/xtjs/xcfjs/)；需要定制计算，[三元相图计算](https://www.keyanxueshu.com/) 上有入口。