4. 相变动力学建模基础
在热处理组织应力耦合仿真中，相变动力学模型是连接温度场与组织场、进而驱动应力场演化的核心环节。若把温度场比作“热历史”，那么相变动力学就是“组织响应函数”：给定温度路径，它回答“转变量有多少、转变速率有多快、残余奥氏体还剩多少”。本节围绕扩散型相变（铁素体、珠光体、贝氏体）与切变型相变（马氏体）的数学描述展开，并说明等温转变与连续冷却转变在数值实现上的差异。
一、扩散型相变：从等温孕育到JMAK方程
铁素体、珠光体和贝氏体的共同特征是依靠原子扩散完成点阵重构，因此转变需要“孕育时间”，且转变量随温度呈C曲线特征。描述等温转变最经典的模型是Johnson-Mehl-Avrami-Kolmogorov（JMAK）方程：
X(t) = 1 - exp( -k(T) · t^n )其中X为转变体积分数，t为等温保持时间，n为Avrami指数（反映形核与长大机制），k(T)为与温度相关的速率常数，通常写成Arrhenius形式：
k(T) = k0 · exp( -Q / (R·T) )Q为有效激活能，R为气体常数，T为绝对温度。对铁素体转变，n常取1.0～1.5；珠光体转变n约2～3；贝氏体转变n约1.5～2.5。这些取值并非固定，需结合材料与实验C曲线标定。JMAK方程只适用于等温、单一机制主导的阶段，若同一温度下先后发生铁素体+珠光体+贝氏体，需要按“竞争转变”处理：各相独立计算虚拟转变量，再按可用奥氏体分数归一化。
二、连续冷却转变：从等温到非等温的加和法则
实际热处理多为连续冷却，温度随时间变化，不能直接套用等温JMAK。工程上常用Scheil加和法则（也称孕育期加和）判断转变开始：把连续冷却路径离散为若干等温小段，对每一段计算“消耗的孕育分数”：
Σ Δt_i / τ_s(T_i) ≥ 1 → 转变开始其中τ_s(T_i)为该温度下的孕育期。转变开始后，转变量可用“虚拟时间”法：把当前已转变分数X反解为等效等温时间t_eq，再叠加当前时间步Δt，最后用JMAK重新计算X。这样既保留了JMAK的物理意义，又能处理连续冷却。对珠光体和贝氏体，还需注意二者C曲线重叠区：同一温度下可能同时具备两种转变的驱动力，数值上通常按各自速率竞争消耗奥氏体。
三、切变型相变：马氏体的无扩散描述
马氏体转变不依赖扩散，而是通过切变与协作原子位移完成，因此没有孕育期，转变量主要取决于温度，且具有“变温转变”特征。最常用的描述是Koistinen-Marburger（K-M）方程：
X_m(T) = 1 - exp( -α · (M_s - T) ) (T < M_s)其中M_s为马氏体开始转变温度，α为材料常数（碳钢常取0.011～0.02 K⁻¹）。当T≥M_s时X_m=0。K-M方程假设冷速足够快、无扩散型竞争，且马氏体转变不可逆。若冷速不足，部分奥氏体可能先转变为铁素体/珠光体/贝氏体，剩余奥氏体才在M_s以下转变为马氏体，此时需先扣除已转变分数，再对剩余奥氏体应用K-M方程。
四、数值实现：时间步进与状态更新
在有限元或有限差分框架中，相变动力学通常以“单元级状态变量”形式嵌入热-组织耦合循环。每个时间步Δt内，先由温度场插值得到当前温度T，再按以下顺序更新：
- 判断当前温度区间：若T高于A1，奥氏体化；若处于铁素体/珠光体/贝氏体C曲线区间，按扩散型模型计算；若T<M_s，按K-M方程计算马氏体。
- 对扩散型相变，用Scheil加和判断孕育是否完成；未完成则只累积孕育分数，不更新转变量。
- 孕育完成后，用虚拟时间法更新各相转变量，并保证各相分数之和不超过1。
- 对马氏体，直接按当前温度计算X_m，并更新残余奥氏体分数。
- 将更新后的相分数传递给力学模块，用于计算相变塑性、体积变化和热应变。
下面用一张图概括本节的核心逻辑：温度路径如何分别进入扩散型与切变型两条数学通道，并最终汇合为组织分数与应力耦合输入。
五、参数标定与常见陷阱
模型能否用于仿真，关键在参数标定。JMAK的k0、Q、n通常由等温膨胀实验或C曲线反算；M_s和α可由膨胀曲线或DSC确定。实际标定时建议先固定n，再拟合k(T)，最后用连续冷却实验验证Scheil加和是否合理。常见陷阱有三类：一是把等温JMAK直接用于连续冷却而不做虚拟时间修正，导致转变开始偏晚；二是忽略残余奥氏体，使马氏体分数虚高；三是各相独立计算后未归一化，导致总分数超过1。只要在时间步进中坚持“先判断、再累积、后归一化”，就能避免大部分数值异常。
综上，扩散型相变的核心是“时间—温度—转变量”的JMAK与Scheil加和，切变型相变的核心是“温度—转变量”的K-M方程。两者在数值实现上共享同一时间步进框架，但判断逻辑与状态更新方式不同。理解这一分工，是后续组织应力耦合仿真能够稳定收敛的前提。