第4章：TDB文件格式——热力学计算的“源代码”
说实话，刚接触Thermo-Calc那会儿，我最怕的就是TDB文件。看着那一堆符号、数字、逗号，感觉像在看天书。后来做项目做多了，才发现这东西其实挺有规律的。今天我就带你把它拆开看看。
TDB文件，全称是Thermodynamic Database文件。说白了，它就是热力学计算的“源代码”。你算相图、算热力学性质，底层调用的都是TDB文件里的数据。没有它，Thermo-Calc就是个空壳子。
核心观点：TDB文件 = 元素定义 + 相定义 + 参数定义。这三块搞懂了，你就掌握了自定义数据库的钥匙。
4.1 TDB文件的基本结构
一个标准的TDB文件，长什么样？我直接给你看个例子：
$ 这是一个简单的TDB文件示例
$ 注释行以$开头
ELEMENT CU FCC_A1 63.546 0.0 0.0!
ELEMENT NI FCC_A1 58.6934 0.0 0.0!
PHASE LIQUID:L % 1 1.0!
PHASE FCC_A1:F % 1 1.0!
PARAMETER G(LIQUID,CU;0) 298.15 +5194.0-5.849*T; 6000 N !
PARAMETER G(LIQUID,NI;0) 298.15 +11235.0-10.209*T; 6000 N !
嗯，看着有点乱是吧？别急，我们一条一条拆。
文件里以$开头的行是注释。我个人习惯在每段数据前加一段注释，说明这段数据是干什么的。比如：
$ 定义纯元素
$ 格式：ELEMENT 名称 晶格类型 原子量 H298 S298
这样过几个月回来再看，自己还能看懂。
4.2 元素定义
元素定义是TDB文件的基石。格式很简单：
ELEMENT 元素符号 参考相 原子量 298.15K焓值 298.15K熵值 !
注意最后那个感叹号，它是结束符。我刚开始写的时候经常漏掉，结果Thermo-Calc报错报得我怀疑人生。
举个例子：
ELEMENT CU FCC_A1 63.546 0.0 0.0!
ELEMENT NI FCC_A1 58.6934 0.0 0.0!
ELEMENT AL FCC_A1 26.9815 0.0 0.0!
这里FCC_A1是参考相。为什么选它？因为铜、镍、铝在室温下都是FCC结构，用这个作为参考态最方便。
我的经验：焓值和熵值通常设为0.0，除非你有特殊需求。这是行业惯例，别乱改。
4.3 相定义
相定义告诉Thermo-Calc：这个体系里有哪些相，每个相由哪些元素组成，比例是多少。
格式：
PHASE 相名称:晶格类型 % 晶格数 化学计量比!
比如：
PHASE LIQUID:L % 1 1.0!
PHASE FCC_A1:F % 1 1.0!
PHASE BCC_A2:B % 1 1.0!
这里LIQUID:L中的:L是晶格类型标识符。L代表液相，F代表FCC，B代表BCC。这些标识符不是随便写的，Thermo-Calc有自己的一套命名规则。
我曾经在项目中遇到一个问题：定义了一个新相，但计算时Thermo-Calc死活不认。查了半天，发现是晶格类型标识符写错了。嗯，这种坑踩过一次就记住了。
避坑指南：相名称不要用中文，不要用特殊符号，最好只用字母和下划线。我见过有人用空格，结果Thermo-Calc直接崩溃。
4.4 参数定义
参数定义是TDB文件里最复杂、也最核心的部分。它定义了每个相的热力学性质，比如吉布斯自由能。
最常见的参数是G参数，代表吉布斯自由能。格式：
PARAMETER G(相名称,元素组合;序号) 温度范围 表达式; 最高温度 N !
看个例子：
PARAMETER G(LIQUID,CU;0) 298.15 +5194.0-5.849*T; 6000 N !
这个参数的意思是：
- 液相中纯铜的吉布斯自由能
- 温度范围从298.15K开始
- 表达式是
+5194.0-5.849*T - 最高适用温度6000K
N表示这个参数是新的，不是对已有参数的修改
为什么用+5194.0-5.849*T这种形式？因为热力学数据通常用多项式拟合。温度T是变量，系数是实验拟合出来的。
对于二元交互作用参数，格式类似：
PARAMETER G(LIQUID,CU,NI;0) 298.15 +10000.0-2.0*T; 6000 N !
这个参数描述的是液相中Cu-Ni二元交互作用对吉布斯自由能的贡献。
关键点：参数表达式中的T是绝对温度，单位是开尔文。所有能量单位是焦耳/摩尔。别搞混了。
4.5 自定义数据库入门
了解了上面三块，你就可以开始写自己的数据库了。我建议你从最简单的二元体系开始。
比如，你想做一个Cu-Ni二元数据库：
$ Cu-Ni二元数据库
$ 元素定义
ELEMENT CU FCC_A1 63.546 0.0 0.0!
ELEMENT NI FCC_A1 58.6934 0.0 0.0!
$ 相定义
PHASE LIQUID:L % 1 1.0!
PHASE FCC_A1:F % 1 1.0!
$ 纯元素参数
PARAMETER G(LIQUID,CU;0) 298.15 +5194.0-5.849*T; 6000 N !
PARAMETER G(LIQUID,NI;0) 298.15 +11235.0-10.209*T; 6000 N !
PARAMETER G(FCC_A1,CU;0) 298.15 0.0; 6000 N !
PARAMETER G(FCC_A1,NI;0) 298.15 0.0; 6000 N !
$ 交互作用参数
PARAMETER G(LIQUID,CU,NI;0) 298.15 +10000.0-2.0*T; 6000 N !
PARAMETER G(FCC_A1,CU,NI;0) 298.15 +8000.0-1.5*T; 6000 N !
保存成Cu-Ni.TDB，然后在Thermo-Calc里用DATABASE RETRIEVAL模块加载它，就能用了。
我的建议：刚开始别贪多。先做两个元素的，跑通了再慢慢加。我见过有人一上来就想做10元数据库，结果数据打架，根本算不出来。
4.6 知识体系总览
下面这张图帮你理清TDB文件的核心逻辑：
4.7 常见错误与调试
写TDB文件时，你肯定会遇到错误。我总结几个最常见的：
| 错误类型 | 表现 | 解决方法 |
|---|---|---|
| 缺少结束符 | Thermo-Calc报"Unexpected end of file" | 检查每行末尾是否有! |
| 元素未定义 | 报"Unknown element" | 检查元素名称是否拼写正确，是否在ELEMENT行定义过 |
| 温度范围重叠 | 报"Temperature range conflict" | 检查同一参数的不同温度段是否有重叠 |
| 相名称重复 | 报"Duplicate phase definition" | 检查是否重复定义了同一个相 |
避坑指南：我曾经花了一整天调试一个TDB文件，最后发现是ELEMENT写成了ELEMNET。Thermo-Calc对拼写非常敏感，一个字母都不能错。
调试小技巧：先用DATABASE RETRIEVAL模块加载你的TDB文件，如果加载成功，说明基本语法没问题。然后试着计算一个简单的二元相图，看看结果是否合理。如果相图形状怪异，多半是参数有问题。
好了，TDB文件的基本内容就这些。你想想看，其实它就是一个结构化的文本文件，把热力学数据按照固定格式写进去而已。多写几次，你就能摸清它的脾气了。