VPSC耦合元胞自动机模拟金属热压缩的再结晶行为
传统动态再结晶模型往往更关注宏观流动应力或平均晶粒尺寸，例如通过临界应变、Avrami 方程或经验型晶粒尺寸方程描述软化行为。这类方法虽然计算方便，但很难回答一个关键问题：不同取向晶粒在变形过程中如何累积位错？哪些晶粒更容易发生再结晶？新生晶粒又如何影响后续织构演化？这篇文章正是针对这一问题，提出了一种基于 VPSC 的 DDRX 模拟方法。核心流程框架如下：
模型的基本思想是：将每个晶粒看作嵌入均匀等效介质中的黏塑性夹杂，通过 VPSC 计算晶粒尺度的应力、应变率、滑移活动和取向演化。在此基础上，作者引入位错密度硬化模型，用 Kocks-Mecking 型方程描述位错累积与动态回复过程。当某个晶粒的位错密度达到临界值后，该晶粒具备动态再结晶形核条件。随后，模型通过形核概率判断是否生成新的 DRX 晶粒，并赋予新晶粒较低的初始位错密度和新的初始取向。
文章中的另一个重要部分是晶粒长大。由于 VPSC 是平均场模型，不包含真实晶粒邻接关系，因此作者采用平均场晶界迁移假设：晶粒长大驱动力来自单个晶粒位错密度与整体平均位错密度之间的差值。低位错密度的新生 DRX 晶粒倾向于长大，而高位错密度的变形母晶逐渐被消耗。通过不断更新晶粒半径和体积分数权重，模型可以同时描述 DRX 体积分数、平均晶粒尺寸和织构演化。
作者将该模型应用于 OFHC 铜的热变形模拟。结果表明，模型能够较好地预测流动应力曲线中的峰值和后续软化行为，也能反映动态再结晶导致的晶粒细化过程。更重要的是，模型能够进一步分析织构演化：塑性变形会使晶粒取向向特定稳定取向旋转，而动态再结晶新晶粒的不断形成会削弱原有织构强度，使整体织构趋于分散。这一点体现了多晶塑性模型相比传统经验型 DRX 模型的优势。
这篇文章给我们的启示主要有三点。第一，动态再结晶建模不应只停留在宏观应力软化层面，而应尽可能联系到晶粒尺度的位错密度、取向演化和体积分数变化。第二，VPSC 这类平均场晶体塑性模型虽然没有空间分辨能力，但计算效率高，适合进行多参数、多温度、多应变率条件下的组织演化预测。第三，对于 FCC 金属材料，尤其是铜、铝合金以及镍基高温合金，位错密度驱动的 DDRX 模型具有较强的可扩展性。
总体来看，这篇文章的最大贡献不是提出某一个复杂公式，而是建立了一条清晰的建模路线：以 VPSC 为力学骨架，以位错密度为内部变量，以临界位错密度控制形核，以储能差驱动晶粒长大，最终实现热变形过程中力学响应与微观组织演化的耦合预测。对于从事晶体塑性、动态再结晶和高温合金热加工模拟的研究者来说，这是一篇非常值得精读和复现的文章。
使用作者提供的完整公式和计算流程。我们可以使用最新的vpsc8.0为基础，并在此上面完成二次开发，加入再结晶影响。在加入时可以将再结晶形核和生长做成相对独立的模块，这样可以不破坏原始结构。
使用原始模型和参数，以及对应的温度(775K)和晶粒尺寸(13μm)进行简单的压缩模拟（ND方向压缩100%）。初始300个随机晶粒
考虑DRX：
初始随机织构：
开启再结晶后的热变形模拟织构：
开启再结晶后的热变形晶粒数量演化：
开启再结晶后的热变形应力应变曲线：
不考虑DRX：
热变形模拟织构：
热变形模拟应力应变曲线：
