CA法模拟动态再结晶,晶粒正常长大,利用元胞自动机生成拓扑晶粒模型,参数可调。
元胞胞自动机模拟动态再结晶母相晶粒生成。
本程序基于曲率驱动机制以及热激活机制,matlab编写,本程序模拟奥氏体晶粒正常长大过程。
程序均有注释,仅作学习交流使用
这段程序主要是用于生成晶粒取向随机分布的均匀化晶粒组织。程序的主要流程如下:
1. 用户输入元胞空间大小(Nx和Ny)、形核点数目(numnucl)和随机数种子(myseed),以及晶粒均匀化长大步数(nstep)。
2. 设置随机数种子。
3. 初始化元胞数组ori,并在其中随机生成numnucl个形核点,形核点的取向值分布范围为1-180。
4. 初始化当前元胞数组的邻居数组ori2、ori4、ori6、ori8。
5. 当元胞数组中存在空白元胞时,进行以下操作:
- 调用neibo4_mex函数生成当前元胞数组的4邻居数组ori2、ori4、ori6、ori8。
- 调用growth_mex函数对当前元胞数组进行晶粒均匀化长大。
6. 用户选择是否引入第二相颗粒,如果是,则输入第二相颗粒数目(number_phase)和尺寸(radius),并调用addphase函数将第二相颗粒添加到元胞数组ori中。
7. 进行nstep次晶粒均匀化长大:
- 调用neibo8_mex函数生成当前元胞数组的8邻居数组ori1、ori2、ori3、ori4、ori6、ori7、ori8、ori9。
- 调用coarse_mex函数根据晶界能减小趋势使晶粒长大。
8. 根据用户选择是否引入第二相颗粒,调用相应的函数a