通过对 Bauschinger 效应中背应力演化机制的实验与模拟分析，可以揭示短程效应和长程效应如何共同作用于材料的塑性行为，尤其是背应力。基于这些发现，本节能够进一步探讨如何通过数值模拟来预测和分析这些现象。在多晶材料中，循环载荷下的残余应力累积与晶粒间的异质变形密切相关。这种微尺度的非协调变形不仅导致宏观屈服面的偏移，还在晶界区域形成了复杂的能量存储与耗散过程。
图 3.16 ABAQUS 软件 Python 二次开发平台多尺度有限元流程图
为进一步量化晶粒数量对残余应力及其储能分配（ESR）的影响，本节研究结合 Abaqus Python 二次开发平台，通过 Python 脚本驱动参数化建模与批量计算，建立 RVE 进行多场耦合拟合与优化。基于 Abaqus 二次开发平台，集成 UMAT[83]子程序，旨在实现钛合金多晶集合体循环软化行为的晶粒尺度模拟。在进行晶体塑性有限元模拟时，结合 ABAQUS 软件与 Python 二次开发可以显著提升模拟分析的效率和灵活性。通过 Python 脚本，可以实现模型的自动化处理，例如生成多晶几何模型、定义材料属性和加载条件等。尽管 ABAQUS 本身没有内置显示晶界的功能，但通过 Python 二次开发可以创建插件来显示晶界，从而使多晶塑性分析更加直观。此外，使用 Voronoi 模型可以生成随机取向的晶粒，而通过 Python 脚本可以实现可视化以及多尺度晶体塑性有限元模型的构建（如图3.16）。
图 3.17 示意性多晶显示通过三个晶粒的路径，沿着这三个晶粒确定应力、取向和边界条件。
图 3.18 不同材料构成的多晶模型 RVE：(a上)双晶粒 RVE；(b下)十晶粒 RVE。
Szczepinski[82]的双晶模型是用于解释金属材料中由于残余应力而产生的残余微应力储能的一个重要工具。本文使用 Abaqus Python 脚本二次开发平台基于Chaboche 理想弹塑性本构来构建双晶（图 3.18(a)）与十晶粒的 RVE 模型（图3.18(b)），用于复现说明并优化 Szczepinski 双晶模型，同时验证该理论在多晶条件下的普适性。其次是对该模型的边界条件与循环载荷的施加（图 3.19），为了显著提高计算稳定性以及仿真结果的合理性进行以下的设置：首先，在 X=0、Y=0、Z=0 三个基准面分别固定对应方向位移（如 X=0 面仅限制 U1 自由度），基于对称约束的准则，这样既可以有效消除刚体位移又可以保留试样在非约束方向的自然变形（X 轴），同时也可以合理避免传统全约束导致的应力失真；随后采用了动态生成的严格单调递增的时间序列幅值曲线（每次循环周期为 1s），同时增设了时间戳（ε=1e-9）用于确保数据点的唯一性，同时有效规避传统循环加载引起的数据振荡问题，从而在准确控制每周期内“拉伸-卸载-压缩-卸载”四个阶段的加载路径。
图 3.19 RVE 载荷及边界条件施加示意图
2 晶粒数量对残余应力及储能演化的影响
以下是对这个模型的详细解释和相关公式通过构建由异质屈服强度晶粒(≠ )组成的理想化单元（图 3.20(a)），揭示了多晶材料在循环加载中残余应力与能量演化的微观机制。
图 3.20 双晶理论模型与有限元结果：(a) Szczepinski 理论模型；(b) Szczepinski 有限元迟滞回线；(c) Mughrabi 理论模型；(d) Mughrabi 有限元迟滞回线
图 3.21 Mughrabi 双晶循环应力云图:(a)Mises 应力云图；(b)S11 应力云图。
图 3.22 十晶粒 RVE 等效塑性应变 PEEQ 演化云图
图 3.23 Mughrabi 双晶模型的残余应力储能演化曲线：(a)随时间演化曲线；(b)随应变演化曲线
图 3.26 多晶 RVE 循环加载应力-应变曲线：(a上)每个晶粒的迟滞回线；(b下)循环硬化演化曲线
图 3.30 RVE 裂纹萌生及扩展示意图：(a)-(d)表示不同循环阶段的裂纹扩展情况
【注意：推文仅为技术分享，不售卖模型代码，如需代做和辅导请添加下面老师，答疑免费】
