摘要：本文建立了纯锌在霍普金森冲击加载下的动态再结晶元胞自动机（CA）模型， 系统研究了不连续动态再结晶（DDRX）、连续动态再结晶（CDRX）和孪晶诱 导动态再结晶（TDRX）三种机制的微观演化过程。模型采用两相位错密度方 法，分别跟踪未再结晶区域和再结晶区域的位错密度演化，成功捕捉了峰后 再硬化现象。基于EBSD表征的初始组织，模拟了粗晶、预变形亚晶和孪晶组 织在高应变率（3200 s⁻¹）变形下的再结晶行为。
项目概览
项目简介
动态再结晶是金属在高应变率变形过程中的重要组织演化机制，直接影响材 料的力学性能和加工工艺。纯锌作为典型的密排六方金属，在霍普金森冲击 加载下呈现复杂的再结晶行为，包括不连续动态再结晶（DDRX）、连续动态 再结晶（CDRX）和孪晶诱导动态再结晶（TDRX）三种主要机制。本文建立了 基于元胞自动机方法的微观组织演化模型，系统研究了纯锌在3200 s⁻¹应变 率下三种再结晶机制的物理过程和演化规律。
模型采用两相位错密度方法，分别跟踪未再结晶区域和再结晶区域的位错演 化，成功捕捉了实验中观察到的峰后再硬化现象。基于电子背散射衍射 （EBSD）表征的真实初始组织，模拟了粗晶、预变形亚晶和孪晶组织的微观 演化过程。结果表明，三种机制的再结晶启动应变存在显著差异：TDRX最早 （ε=0.07），DDRX次之（ε=0.10），CDRX最晚（ε=0.27），该顺序与各机 制的临界位错密度完全对应。DDRX和TDRX均表现出明显的峰后再硬化 （4-6 MPa），验证了两相位错密度模型的有效性。CDRX呈现连续的亚晶旋转 转化特征，低角晶界逐渐转化为高角晶界。
针对模型参数标定问题，提出了两阶段标定算法，将峰前硬化参数与峰后再结晶参数解耦优化，采用多起点全局优化策略实现了100%的收敛率。CDRX模型标定精度达到优秀级别（NRMSE=0.81%），DDRX和TDRX模型分别达到良好（NRMSE=5.77%）和可接受（NRMSE=10.54%）水平。
本研究为理解纯锌及其他HCP金属在高应变率下的动态再结晶行为提供了有效的数值工具，对优化热机械加工工艺和预测组织演化具有重要意义。
系统架构
本系统采用模块化架构，由核心计算层、分析标定层和可视化层组成。核心层实现基于两相位错密度模型的CA演化算法，支持DDRX/CDRX/TDRX三种机制；分析层提供两阶段自动标定算法，实现参数解耦优化；可视化层自动生成组织演化、应力应变曲线等多类图表。系统采用EBSD真实初始组织，通过主控脚本统一调度，结果结构化保存，具备100%可重复性。代码总量3651行，测试覆盖率70%，文档完整，适合学术研究和二次开发。
图1 系统架构图
技术创新
创新点1：两相位错密度CA模型
提出了分别跟踪未再结晶区域和再结晶区域位错密度演化的两相CA模型， 成功捕捉了实验中观察到的峰后再硬化现象（4-6 MPa）。
创新点2：三种DRX机制的系统对比研究
建立了DDRX、CDRX、TDRX三种机制的统一CA建模框架，揭示了再结晶 启动顺序（TDRX < DDRX < CDRX）与临界位错密度的完美对应关系。
创新点3：两阶段标定算法
提出了峰前硬化参数与峰后再结晶参数解耦的两阶段标定算法，实现了100%的标定收敛率，CDRX模型精度达到优秀级别（NRMSE=0.81%）。
创新点4：基于EBSD真实初始组织的CA模拟法
采用EBSD表征的真实微观组织作为CA模拟的初始条件，针对粗晶、预变形亚晶和含孪晶三种组织建立了专门的初始化方法。
快速开始
在MATLAB命令窗口运行 main_run_all 即可一键执行三种再结晶机制的模拟，结果自动保存在outputs/matlab_results/目录。
环境要求
硬件要求任意能运行MATLAB的计算机（8GB内存推荐），软件要求 MATLAB R2018b或更高版本（包含基础工具箱），存储空间项目70MB + 输出约 30MB，操作系统Windows、Linux或macOS均可。
运行展示
运行main_run_all.m
1.粗晶DDRX模型（coarse/目录，7张图）
图2 粗晶模型组织演化图
图3 粗晶模型应力-应变曲线图
图4 粗晶模型位错密度演化
图5 粗晶模型再结晶体积分数
图6 粗晶模型平均晶粒尺寸演化
图7 粗晶模型最终状态多场分布
图8 粗晶界附近形核位点专用可视化
2.亚晶CDRX模型（subgrain/目录，7张图）
图9 亚晶模型组织演化图
图10 亚晶模型应力-应变曲线
图11 亚晶模型位错密度演化
图12 亚晶模型再结晶体积分数
图13 亚晶模型平均晶粒尺寸演化
图14 亚晶模型最终状态多场分布
图15 亚晶模型低角/高角晶界演化专用可视化
3.孪晶TDRX模型（twin/目录，7张图）
图16 孪晶模型组织演化图
图17 孪晶模型应力-应变曲线
图18 孪晶模型位错密度演化
图19 孪晶模型再结晶体积分数
图20 孪晶模型平均晶粒尺寸演化
图21 孪晶模型最终状态多场分布
图22 孪晶界附近形核位点专用可视化
4.三模型对比图（目录，3张）
图24 三种模型应力-应变曲线对比
图25 三种模型再结晶分数对比
图23 三种模型平均晶粒尺寸对比
项目资源
配套文档
点击查看：基于MATLAB低采样率ISAR成像的快速稀疏重建算法研究 注意：免费提供！
配套文件
包括完整的项目源代码、演示视频、运行截图，开箱即用。
项目信息
作者信息
作者：Bob (张家梁)
项目编号：SD-1-2z-M
原创声明：本项目为原创作品
联系方式
开源协议
本项目采用AGPL-3.0开源协议，允许个人和组织自由使用、修改和分发代码，但基于本项目的衍生作品必须同样开源，且用于提供网络服务时需向用户提供完整源代码。本项目仅供学习研究使用，作者不对使用本项目产生的任何后果承担责任，使用者应遵守当地法律法规，合理合法使用本项目。如本项目对您的研究或工作有所帮助，欢迎引用并注明出处。
评论(0)