项目背景滚压表面强化技术是提高零件疲劳寿命和耐磨性的关键工艺其通过引入梯度塑性变形在表层实现晶粒细化和残余压应力从而极大地增强零件的机械性能。尽管该技术在工程实践中应用广泛但在滚压过程中梯度应变场驱动的位错演化与连续动态再结晶cDRX晶粒细化之间的耦合机制仍不明确这也限制了工艺参数的优化设计。为了填补这一研究空白我们通过数值模拟深入探讨了这一过程中的关键机制。研究目标本研究旨在通过数值模拟揭示滚压表面强化过程中梯度应变场如何驱动位错演化与晶粒细化。我们将重点分析不同滚压力下的位错演化与cDRX晶粒细化行为并为工艺参数优化提供理论依据。项目简介我们采用42CrMo合金结构钢作为研究对象构建了一个多尺度数值模拟框架结合Hertz接触理论、Kocks-Mecking位错模型和元胞自动机CA方法模拟了不同滚压力下的位错演化与cDRX晶粒细化过程。通过这一方法揭示了梯度应变场对晶粒细化的关键驱动作用。模拟框架滚压载荷场模型基于Hertz接触理论计算不同滚压力1000~4000N下的梯度应变、应变率和温度场分布。位错密度演化采用Kocks-Mecking位错模型描述位错密度随应变的变化并引入温度与应变率的耦合修正建立了错配角的累积模型。晶粒细化模拟基于元胞自动机CA方法模拟cDRX晶粒细化过程综合考虑位错驱动的晶粒分裂与温度驱动的晶界迁移。模拟结果滚压力的影响随着滚压力的增加表层的等效塑性应变和位错密度显著升高梯度分布呈指数衰减。增大滚压力不仅促进位错密度的增殖还加剧了表层的塑性变形。cDRX的激活当亚晶错配角超过高角晶界阈值且局部应变满足临界条件时cDRX过程被激活表层晶粒发生显著细化。细化层的厚度与晶粒尺寸随着滚压力从1000N增至4000N细化层厚度和表层晶粒细化程度逐渐增加。表层晶粒尺寸显著减小而深层保持粗晶形成明显的梯度晶粒结构。位错密度对晶粒长大的影响高位错密度对晶粒长大有抑制作用有助于维持细晶组织的稳定性。结果展示EBSD风格晶粒图通过EBSD风格的晶粒图可以直观展示不同滚压力下的梯度晶粒结构。表层晶粒逐渐细化深层保持粗晶清晰呈现细晶与粗晶的过渡。晶粒尺寸与深度的关系随着滚压力的增加表层晶粒尺寸逐渐减小而深层晶粒保持较大尺寸。位错密度与深度分布位错密度随着滚压力的增大显著增加且呈现出半对数关系。等效塑性应变与深度的分布随着滚压力的增加表层的塑性应变显著升高体现了滚压过程的塑性变形特征。细化层厚度与晶粒尺寸统计柱状图柱状图展示了细化层的厚度与表层晶粒尺寸的变化趋势进一步揭示了不同滚压力下的晶粒演化规律。系统架构本项目采用模块化MATLAB程序结构具体模块包括calculate_rolling_field用于计算梯度应变、应变率及温度场基于Hertz接触理论。dislocation_evolution实现位错密度的演化计算基于Kocks-Mecking模型加入错配角的累积计算。CA_grain_evolution基于元胞自动机方法模拟cDRX晶粒细化过程综合考虑位错驱动的晶粒分裂与晶界迁移。快速开始在MATLAB中进入项目目录cd D:\rolling_induced_gradient_grain_refinement。执行main启动完整模拟程序将自动遍历四种滚压力工况1000~4000N依次计算载荷场、位错演化和晶粒演化。结果图像和数据文件将保存在results/目录。环境要求MATLAB版本R2024b及以上。工具箱无需额外工具箱依赖MATLAB内置函数如movmean、hsv2rgb、containers.Map等。模拟验证运行run_quick_test进行7组模块验证结果自动保存到results/目录。查看结果结果图像存放在results/figures/目录包括EBSD风格图、晶粒尺寸与深度曲线、位错密度曲线、应变分布图等。数据存放在results/data/目录包括simulation_results.mat和results_summary.csv。项目亮点梯度晶粒细化的深入理解通过数值模拟揭示了滚压过程中的微观机制特别是位错演化与cDRX的耦合过程。多材料适用性该模拟框架不仅适用于42CrMo钢还可扩展至其他材料如Ti6Al4V、304SS等。可视化与数据分析通过EBSD风格图直观展示晶粒结构结合统计数据提供定量分析结果。项目总结本研究通过数值模拟揭示了滚压表面强化过程中梯度应变场如何驱动位错演化与晶粒细化的机制。该方法为滚压工艺参数优化提供了理论支持同时也为提高零件的疲劳寿命和耐磨性提供了有力的技术依据。通过这一数值模拟框架我们成功地实现了对滚压过程微观机制的深入理解并为实际应用中的工艺设计提供了宝贵的指导。这个项目为如何通过数值模拟方法揭示滚压表面强化过程中的细节提供了全新的视角期待更多的工程领域能够借鉴这一成果进一步优化零件制造工艺。