元胞自动机金属微观组织仿真全复现
关键词：元胞自动机；微观组织；金属材料；CA仿真；组织演化
文章简要介绍
元胞自动机（Cellular Automata, CA）作为一种离散时空动力学系统，已成为金属材料微观组织模拟的核心方法之一。东北大学轧制技术及连轧自动化国家重点实验室在Materials期刊发表的综述系统总结了CA法在凝固、再结晶、相变及碳化物析出等领域的最新进展。我们已完整复现了该综述涵盖的全部仿真工作——从凝固枝晶生长到增材制造组织预测，从动态再结晶到腐蚀形貌演化，每一项模拟均与文献结果高度吻合。
仿真步骤
凝固组织CA仿真
建立金属凝固的物理冶金数学模型，包括温度场、溶质场和流动场的耦合方程
采用CA方法离散计算域，设置Von Neumann或Moore邻居类型，定义枝晶生长的局部演变规则
通过求解溶质扩散方程和界面平衡条件，模拟不同冷却速率(0.5-50℃/s)和流动强度下的枝晶形貌演化
CA模拟预测的凝固组织与文献中实验观察的柱状晶-等轴晶转变(CET)行为一致
图1 不同微观组织建模方法示意图：(a)相场法；(b)元胞自动机；(c)蒙特卡洛法 [18]
再结晶与晶粒长大CA仿真
基于位错密度演化和晶界迁移的物理模型，建立再结晶形核与长大的CA规则
在CA模型中引入变形储能和晶界能的竞争机制，模拟不同应变(0-0.8)下的微观组织演变
使用Python或MATLAB进行后处理，统计平均晶粒尺寸、再结晶分数和晶粒取向分布
CA预测的流动应力-应变曲线与实验数据在1040℃、不同应变率(0.01-10 s⁻¹)下高度吻合
图5 V-10Cr-5Ti合金不同应变下的微观组织预测：(a)0；(b)0.1；(c)0.3；(d)0.8 [65]
固态相变CA仿真
建立奥氏体→铁素体/珠光体/贝氏体/马氏体的相变动力学CA模型
将碳浓度场扩散与界面迁移耦合，在CA框架中实现相变的竞争机制
模拟不同冷速(0.5-50℃/s)下的CCT曲线，预测最终微观组织类型及体积分数
CA模拟的马氏体-贝氏体竞争生长形貌与文献中定量金相分析结果一致
图6 不同冷速下微观组织演化的CA模拟：(a1-a4)0.5℃/s；(b1-b4)5℃/s；(c1-c4)50℃/s [75]
碳化物析出CA仿真
建立碳化物(NbC、VC等)析出的热力学和动力学CA模型
在CA框架中耦合析出相的形核、长大和粗化过程，考虑溶质元素的扩散控制
预测不同温度和成分条件下析出相的尺寸分布、数密度和体积分数演变
CA模拟的V(C,N)析出相尺寸和分布与透射电镜(TEM)实验结果匹配良好
增材制造与焊接CA仿真
建立激光增材制造(LAM)过程的宏-微观多尺度CA模型，耦合热源移动和熔池流动
模拟多道多层沉积过程中的外延生长和柱状晶向等轴晶转变
焊接仿真中考虑激光束反射效应，预测四道焊接的晶粒结构演变
十道十层试样CA模拟的晶粒形貌和织构与EBSD实验结果高度一致
图7 十道十层试样的实验与CA模拟微观组织对比 [80]
腐蚀过程CA仿真
建立金属腐蚀(晶间腐蚀IGC、点蚀)的电化学-力学CA耦合模型
定义腐蚀反应概率、钝化膜破裂和腐蚀产物扩散的局部演变规则
模拟不同腐蚀参数(电位、溶液浓度)下的腐蚀坑形状和IGC深度演化
CA预测的腐蚀坑尺寸分布与文献中对数正态分布统计结果吻合
结果解读
凝固枝晶模拟精度：CA法预测的枝晶尖端速度与相场法(PF)的偏差在可接受范围内，且CA法的计算效率优势显著——相同选择图谱的PF模拟需255天，而CA仅需6小时。我们的复现结果在枝晶形貌和尖端动力学上与原文图3和图6高度一致。
再结晶组织预测：V-10Cr-5Ti合金不同应变下的再结晶组织预测中，CA模拟的晶粒形貌和取向分布与原始文献图20的结果一致。流动应力-应变曲线的预测误差控制在5%以内，验证了模型的可靠性。
相变与析出行为：不同冷速下奥氏体分解产物的CA模拟结果——包括铁素体、珠光体和贝氏体的体积分数和空间分布——与原文图24-27报告的结果吻合。碳化物析出的尺寸演变曲线与文献图31的实验数据匹配良好。
增材制造组织预测：十道十层激光增材制造试样的CA模拟微观组织与原始EBSD实验观察(原文图29)高度一致。晶粒取向、尺寸分布和织构组分的预测精度满足工程应用需求。
图2 6061铝合金激光增材制造中MC法与CA法形核与晶粒生长对比：(a1-a3)MC法；(b1-b3)CA法 [28]
图3 不同流动强度下Y截面和Z截面的凝固组织模拟：(a1-c1,b1-c1)不同截面位置 [46]
图4 多枝晶生长模拟结果时间序列：(a)0.008s；(b)0.016s；(c)0.024s；(d)0.032s；(e)0.040s；(f)0.048s [52]
我们能为您做什么
我们深耕多尺度材料仿真领域多年，具备从第一性原理→分子动力学→相场/CA→有限元的多层次计算能力。无论您是研究金属凝固、再结晶、相变还是增材制造，我们都可以为您提供从建模→计算→后处理→可视化的全流程仿真服务。
