动态再结晶CA法模拟基础之初始晶粒生长规则

初始晶粒的生长规则：

热加工过程中，首先需要把钢加热到奥氏体状态，然后在适当的热工艺条件下获得所期望的性能和组织。前人在做动态再结晶的模拟时通常由实验获得初始晶粒的大小，然后按照CA法在实验所获得的结果上进行动态再结晶的模拟。这种处理方法尽管简便，但是没有考虑到材料本身的奥氏体化过程，因此在本算例中把动态再结晶之前的奥氏体化过程简化为一个晶粒正常长大的过程进行模拟。

奥氏体形成过程主要由：奥氏体形核、奥氏体长大、残余渗碳体溶解和奥氏体成分均匀化组成。因此本算例所采用的初始晶粒生长规则[1-4]为：

1.晶粒长大是一个热激活过程，位于晶界时处的原子需要克服一定的能量势垒才能跃迁到新的状态。在同一个CA步中按照概率 确定元胞发生转变的难易程度。

其中C是常数，可以通过T=Tm时转化概率为1，即此时P1进行取值。Tm为材料融化温度，Tac1为奥氏体化温度，R是气体常数，Qb是晶界迁移激活能。以上参数可因材料的不同根据论文进行选取。

由于晶粒长大时遵循以下规律：首先，弯曲晶界趋向平直，即晶界向曲率中心方向移动，以减少表面积，降低表面能；其次三个晶粒的晶界夹角不等于120度时，晶界总是向角度较小的晶粒方向移动，促使三个夹角趋向于120度。基于以上的规则考虑晶界曲率对晶粒长大的影响指定2-4规则。

2.按照8邻居（Moore）类型，如果中心元胞C5(如图1)周围有连续5个或5个以上的元胞为相同的状态，则下一个时间步，中心元胞状态变得与他们相同。

图1
3.如果C2、C4、C6、C8中任意三个状态相同，则在下一个CA模拟时间增量步中，中心元胞C5变得与他们状态相同。

4.如果C1、C3、C7、C9中任意三个状态相同，则在下一个CA模拟时间增量步中，中心元胞C5变得与他们状态相同。（如果3、4同时满足，则由于规则的顺序执行性默认执行后者，即规则4）。
如果以上的规则均不满足，则按照晶界能最小原理（晶界能越小越稳定）按照规则5执行。公式为

其中J是晶界能的度量可取为1， δ为克罗耐克求和符号， k是中心元胞C5的第k邻居， R为邻居总数，在8邻居规则中取8，4邻居则取4。
5.如果以上的任何规则均不满足，则从C5周围的元胞中任意选取一个元胞进行转变，并应用式（2）根据式（3）计算转变后的系统晶界能改变。

参考文献：

[1]. Geiger J, Roósz A, Barkóczy P. Simulation of grain coarsening in two dimensions by cellular-automaton[J]. Acta Materialia, 2001, 49(4): 623-629.

[2]. 陈飞. 热锻非连续变形过程微观组织演变的元胞自动机模拟[D]. 上海交通大学, 2012.

[3]. 关小军, 焦宪友, 周家娟, etc. 单一晶粒长大过程的元胞自动机模拟 %J 中国有色金属学报[J]. 2007, (05): 699-703.

[4]. 花福安, 杨院生, 郭大勇, etc. 基于曲率驱动机制的晶粒生长元胞自动机模型 %J 金属学报[J]. 2004, (11): 1210-1214.