3. 材料热物性与相变参数获取
在热处理组织应力耦合仿真中，材料参数不是“填几个常数”那么简单。温度场求解依赖热导率、比热和密度；相变场求解依赖TTT/CCT曲线、JMAK方程参数和相变潜热；应力场求解则依赖热膨胀系数、相变塑性参数以及马氏体相变的Koistinen-Marburger模型参数。这些参数共同决定了温度—组织—应力三场耦合的计算路径，任何一个参数取值失真，都会在后续的应力结果中被放大。
本节围绕“参数从哪里来、以什么形式进入模型、如何判断可用”展开，重点说明热物性参数与相变参数的获取方法、典型数值范围以及工程处理策略。
一、热物性参数：温度场求解的输入基础
热导率λ、比热cp和密度ρ共同决定热扩散率α=λ/(ρcp)。热扩散率越大，工件内部温度趋于均匀的速度越快，温度梯度越小，热应力也越缓和。对于钢件，热导率随温度升高总体下降，但在相变区间会出现拐点；比热在相变区因潜热释放而出现等效峰值；密度则随温度升高略有下降，并在奥氏体化、马氏体转变时发生体积突变。
获取热物性参数通常有三条途径。第一，查手册或材料数据库，例如ASM、JMatPro、Thermo-Calc等，适合初步建模。第二，采用激光闪射法测热扩散率、差示扫描量热法测比热，再结合密度计算热导率，适合关键材料。第三，在仿真中采用随温度变化的表格曲线，而不是单一常数。以45钢为例，20℃时热导率约50 W/(m·K)，800℃时降至约25 W/(m·K)；比热在20℃约470 J/(kg·K)，在相变区可等效升至800 J/(kg·K)以上。这些数值属于教学示例范围，实际应以材料批次和测试条件为准。
二、相变潜热与膨胀系数：耦合相变与应力的桥梁
相变潜热是组织转变时吸收或释放的热量。在温度场中，它通常以等效比热的形式叠加到cp上，使相变区出现明显的“比热峰”。若忽略潜热，相变区温度会下降过快，组织转变分数被高估，进而影响应力演化。潜热值可通过DSC测量，也可由热力学数据库计算。例如奥氏体向铁素体转变的潜热约为几十kJ/kg量级，马氏体相变潜热相对较小，但在快速冷却中仍不可忽略。
热膨胀系数αth描述温度变化引起的体积变化，而相变本身也会引起体积变化。因此，在耦合模型中通常把“热膨胀”和“相变膨胀”分开处理：热膨胀系数用于温度引起的应变，相变膨胀则通过各相的比容差或相变应变项引入。以碳钢为例，奥氏体比容最小，马氏体比容最大，奥氏体向马氏体转变时体积膨胀约0.5%～1.0%。这个膨胀量虽然不大，但在约束条件下会产生显著的相变应力，是淬火开裂分析的关键来源。
三、TTT/CCT曲线与JMAK方程参数
TTT曲线描述等温转变，CCT曲线描述连续冷却转变。在仿真中，CCT曲线更常用，因为热处理过程多为连续冷却。CCT曲线给出不同冷速下组织转变的起始温度、终止温度和转变量，是判断最终组织组成的第一手依据。
JMAK方程（Johnson-Mehl-Avrami-Kolmogorov）用于描述扩散型相变的转变量随时间和温度的变化：
f = 1 - exp[-(k·t)^n]其中f为转变体积分数，t为时间，k为与温度相关的速率常数，n为Avrami指数。k通常用Arrhenius形式表达：k = k0·exp(-Q/RT)，Q为激活能，R为气体常数，T为绝对温度。n反映形核与长大机制，扩散型相变常见n在1～3之间。实际建模时，k0、Q、n三个参数需要根据TTT/CCT曲线反算或拟合得到，而不是直接照搬文献值。
反算的基本思路是：从CCT曲线读取某一冷速下转变开始和结束对应的温度与时间，代入JMAK方程，联立求解k和n；再用不同冷速的数据做回归，得到k0和Q。若曲线数据不足，可用等温TTT曲线先拟合等温JMAK参数，再通过Scheil叠加原理或加法规则推广到连续冷却。
四、马氏体相变Koistinen-Marburger模型参数
马氏体相变是非扩散型相变，转变量主要取决于温度，而不是时间。Koistinen-Marburger（K-M）模型是最常用的描述形式：
f_M = 1 - exp[-α·(M_s - T)] (T < M_s)其中fM为马氏体体积分数，Ms为马氏体开始转变温度，α为相变速率参数，T为当前温度。当T≥Ms时，fM=0。α通常取0.011～0.03 K-1，具体值与钢种和奥氏体化条件有关。Ms可由经验公式估算，例如Ms(℃)=539-423C-30.4Mn-17.7Ni-12.1Cr-7.5Mo-7.5Si（质量分数，%），也可由CCT曲线或膨胀法实测确定。
K-M模型的关键在于：马氏体转变对温度路径敏感，但对应变速率不敏感。因此在耦合仿真中，马氏体分数通常按温度增量逐步累加，而不是按时间步长直接计算。若冷却过程中出现回温或等温停留，已形成的马氏体不再增加，但可能发生回火分解，这需要在模型中额外考虑。
五、参数一致性检查与常见问题
参数获取完成后，必须做一致性检查。第一，检查温度单位：TTT/CCT曲线常用摄氏度和秒，而JMAK方程中的Arrhenius项必须用开尔文。第二，检查相变潜热是否已经包含在比热曲线中，避免重复叠加。第三，检查Ms与CCT曲线是否自洽：若CCT曲线显示在某一冷速下马氏体开始转变温度明显低于经验公式值，应以实测曲线为准。第四，检查热膨胀系数与相变膨胀是否重复计算：若软件已通过相变应变项处理体积变化，则热膨胀系数应只保留纯热膨胀部分。
常见问题包括：用室温热导率代替高温值，导致温度场计算偏慢；忽略相变潜热，导致相变区温度骤降；JMAK参数直接取文献值而未按本材料CCT曲线反算，导致组织分数偏差；K-M模型中α取值过大，导致马氏体转变在几度内完成，与实际不符。这些问题在耦合仿真中会表现为应力峰值位置偏移、残余应力符号错误或开裂风险误判。
六、可直接使用的参数记录表
为便于建模时逐项核对，下面给出一份参数记录表模板。表中“来源”一栏应填写实测、数据库或反算，“适用条件”应注明材料牌号、奥氏体化温度和冷速范围。
| 参数 | 符号/形式 | 示例值（45钢，教学示例） | 来源 | 适用条件 |
|---|---|---|---|---|
| 热导率 | λ(T) | 20℃: 50 W/(m·K)；800℃: 25 W/(m·K) | 手册/实测 | 45钢，退火态 |
| 比热 | cp(T) | 20℃: 470 J/(kg·K)；相变区等效峰值: 800 J/(kg·K) | DSC/数据库 | 含潜热等效 |
| 密度 | ρ(T) | 20℃: 7850 kg/m³；800℃: 7700 kg/m³ | 手册 | 45钢 |
| 热膨胀系数 | αth(T) | 20～800℃平均: 1.3×10-5 K-1 | 膨胀仪 | 纯热膨胀部分 |
| 相变潜热 | L | 奥氏体→铁素体: 约 60 kJ/kg | DSC/热力学计算 | 扩散型相变 |
| JMAK参数 | k0, Q, n | k0=1×105 s-1，Q=200 kJ/mol，n=1.5 | CCT反算 | 铁素体/珠光体转变 |
| K-M参数 | Ms, α | Ms=330℃，α=0.011 K-1 | 膨胀法/CCT | 马氏体转变 |
表中数值仅用于说明参数形式和量级，实际建模时应以本批次材料的测试或可靠数据库结果替换。完成参数整理后，建议先做一次纯温度场试算，检查相变区温度曲线是否出现合理的潜热平台；再做一次纯相变场试算，检查最终组织分数是否与CCT曲线预测一致；最后才进入三场耦合计算。这样可以把参数错误隔离在早期阶段，避免在应力结果中反复排查。