基于Matlab的晶粒模拟与再结晶演变程序

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. 根据用户选择是否引入第二相颗粒，调用相应的函数addgb_phase_mex或addgb_mex生成含有晶界的元胞数组ori_gb。 9. 保存元胞数组ori和含有晶界的元胞数组ori_gb。 第二段程序是与第一段程序相对应的，用于再结晶演变模拟。程序的主要流程如下： 1. 载入初始组织元胞矩阵ori。 2. 初始化元胞数组ori、晶界数组GB、位错密度数组D、元胞结晶状态数组State和颜色数组C。 3. 进行Nstep次再结晶演变： - 调用DislocationMatrix_mex函数更新位错密度数组D。 - 调用Nucleation_mex函数进行再结晶形核。 - 调用showgb_mex函数将颜色数组C转换为灰度值数组。 - 调用DynamicRX_mex函数使再结晶晶粒长大。 - 调用DislocationDRX_mex函数更新再结晶晶粒的位错密度。 - 调用addgb_mex函数生成含有晶界的元胞数组ori_gb。 4. 统计再结晶更新一次后总的晶粒个数和每个晶粒的尺寸。 5. 保存数据。 这段程序涉及到的知识点包括：随机数生成、数组操作、条件语句、循环语句、函数调用、文件操作等。其中，Mex文件是用于加快运算速度的二进制文件。程序中还涉及到晶粒均匀化、晶粒长大、晶界迁移、位错密度演变等概念和算法。

（基于原始代码解读）

一、程序核心目标与实现框架

本程序通过元胞自动机（CA）方法，在MATLAB环境中模拟奥氏体动态再结晶及晶粒正常长大过程。核心依托曲率驱动机制（控制晶粒边界迁移方向）和热激活机制（调控形核与生长速率），通过模块化设计实现从初始微观组织构建到再结晶演化的全流程模拟，支持参数可调与结果可视化，忠实还原材料热加工过程中微观组织的动态变化。

二、核心模块功能解析

（一）初始微观组织生成模块（`InitialMicrostructure.m`）

该模块用于构建模拟所需的初始晶粒结构，为后续再结晶过程提供基础拓扑模型，具体流程如下：

1. 参数输入：通过命令行交互获取关键参数，包括元胞空间尺寸（Nx、Ny）、初始形核点数量（控制初始晶粒总数）、随机数种子（保证模拟可重复性）、晶粒均匀化迭代步数（优化初始晶粒形态）。
2. 形核与生长：
   - 在空白元胞数组中随机生成指定数量的形核点，每个形核点分配1-180的随机整数作为晶体学取向标识；
   - 基于4邻域规则（neibo4mex函数）驱动形核点向周围空白元胞生长，直至填满整个元胞空间；
   - 通过8邻域迭代（neibo8mex、coarse_mex函数）对初始晶粒进行均匀化处理，减少形态畸变。
3. 第二相颗粒引入：支持用户选择是否添加第二相颗粒，若启用，需输入颗粒数量和半径，程序通过addphase函数在元胞数组中标记第二相区域（标识为-1）。
4. 晶界识别与保存：通过addgbmex或addgbphasemex函数识别晶界（相邻元胞取向不同的位置），生成晶界数组（晶界为1，晶内为0），最终将初始取向矩阵（ori）和晶界矩阵（origb）保存为.mat文件。

（二）动态再结晶模拟模块（`DRXphase.m` 与 `Drx.m`）

该模块为程序核心，模拟动态再结晶的形核、晶粒长大及位错密度演变过程，两者功能一致（Drx.m额外支持拓扑变形开关），具体逻辑如下：

1. 初始化设置：
   - 加载初始组织（ori）和参数配置（Parameters.m），定义关键状态数组：
   -D：位错密度数组（母相初始值1e+10 m⁻²，再结晶晶粒1e-10 m⁻²）；
   -State：元胞状态标识（0=未再结晶，>0=已再结晶）；
   -Grainnumber：晶粒编号（用于统计单个晶粒尺寸）；
   -DRXFraction：再结晶分数（已再结晶元胞占比）。

1. 迭代模拟过程：通过Nstep步循环实现动态演化，每步包含4个核心过程：
   -位错密度演变：调用DislocationMatrixmex函数，基于热激活机制计算母相与再结晶晶粒的位错密度变化，公式为：
   $$D = D + \Delta\varepsilon \times (k1\sqrt{D} - k2D)$$
   （$\Delta\varepsilon$为每步应变量，$k1$、$k2$为硬化/软化系数）。
   -再结晶形核：调用Nucleationmex函数，当母相与位错密度达到临界值（CriticalDislocation）且满足形核概率（Pnuc）时，在晶界处生成再结晶核心，重置位错密度并分配新晶粒编号。
   -再结晶晶粒长大：调用DynamicRXmex函数，基于曲率驱动机制，通过4邻域信息计算晶界迁移速率（与位错密度差相关），判断元胞是否被再结晶晶粒占据并更新状态。
   -拓扑变形（可选）：Drx.m中通过topomex函数，按应变进度调整元胞数组尺寸（新尺寸=初始尺寸×exp(-应变)），模拟实际变形的拓扑变化。

（三）结果统计与可视化模块

贯穿模拟全过程，实现组织演化的实时展示与量化数据输出：

1. 实时可视化：通过image函数绘制当前组织，colormap(colorcube)区分不同取向晶粒，晶界显示为黑色；pause控制动画速度，直观观察再结晶区域扩展。
2. 量化结果计算：
   - 再结晶特征：统计再结晶分数（Fraction）、晶粒数量（NumberOfDRXGrain）、平均面积及直径（MeanAreaOfDRXGrain、MeanDiameterOfDRXGrain）；
   - 晶粒尺寸分布：通过bwlabel和regionprops函数标记晶粒并计算面积、等效直径，生成分布数据（RCA结构体）；
   - 力学参数：计算平均位错密度（AverDislocation）和应力（$\sigma=0.5\times G\times b\times \sqrt{\bar{D}}$，$G$为剪切模量，$b$为伯格斯矢量）。
3. 结果保存：模拟结束后，将所有统计数据保存为data.mat，通过plotexportpdf.m导出组织图像（支持标注与坐标轴隐藏）。

（四）参数配置与辅助函数（`Parameters.m` 及 SourceFunction文件夹）

1. 参数配置（Parameters.m）：集中定义物理参数，包括变形温度（T）、应变速率（StrainRate）、总应变（Strain）、元胞实际尺寸（L0）、临界位错密度（CriticalDislocation）、形核概率（Pnuc）等，支持用户根据材料特性调整。
2. 辅助函数：提供15个核心支撑函数，按功能分为：
   - 邻域与生长：neibo4.m（4邻域）、neibo8.m（8邻域）、growth.m（晶粒生长）等；
   - 晶界与第二相：addgb.m（晶界识别）、addphase.m（第二相生成）等；
   - 统计与可视化：avegrainsize.m（平均晶粒尺寸）、microstructureplot.m（组织绘图）等。

三、程序运行流程

1. 运行InitialMicrostructure.m，输入参数生成初始组织（ori.mat、ori_gb.mat）；
2. 调整Parameters.m中的物理参数（如温度、应变速率）；
3. 运行DRXphase.m（无拓扑变形）或Drx.m（有拓扑变形）启动模拟，实时观察组织演化；
4. 模拟结束后，通过data.mat查看量化结果，或导出图像分析。

四、程序核心特性

1. 物理机制忠实性：严格基于曲率驱动与热激活机制，模拟过程贴合动态再结晶的物理规律；
2. 参数可调性：核心参数集中配置，支持不同材料、工艺条件下的模拟需求；
3. 结果完整性：同步输出可视化组织与量化数据（再结晶分数、晶粒尺寸分布等），直接支撑分析结论。

本程序旨在通过元胞自动机方法，精准还原奥氏体动态再结晶与晶粒长大的微观过程，为材料热加工工艺优化提供数值模拟工具。

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. 根据用户选择是否引入第二相颗粒，调用相应的函数addgb_phase_mex或addgb_mex生成含有晶界的元胞数组ori_gb。 9. 保存元胞数组ori和含有晶界的元胞数组ori_gb。 第二段程序是与第一段程序相对应的，用于再结晶演变模拟。程序的主要流程如下： 1. 载入初始组织元胞矩阵ori。 2. 初始化元胞数组ori、晶界数组GB、位错密度数组D、元胞结晶状态数组State和颜色数组C。 3. 进行Nstep次再结晶演变： - 调用DislocationMatrix_mex函数更新位错密度数组D。 - 调用Nucleation_mex函数进行再结晶形核。 - 调用showgb_mex函数将颜色数组C转换为灰度值数组。 - 调用DynamicRX_mex函数使再结晶晶粒长大。 - 调用DislocationDRX_mex函数更新再结晶晶粒的位错密度。 - 调用addgb_mex函数生成含有晶界的元胞数组ori_gb。 4. 统计再结晶更新一次后总的晶粒个数和每个晶粒的尺寸。 5. 保存数据。 这段程序涉及到的知识点包括：随机数生成、数组操作、条件语句、循环语句、函数调用、文件操作等。其中，Mex文件是用于加快运算速度的二进制文件。程序中还涉及到晶粒均匀化、晶粒长大、晶界迁移、位错密度演变等概念和算法。