精品晶体塑性文章推荐_Abaqus_疲劳_UM_裂纹_理论

精品晶体塑性文章推荐

文章题目：
A benchmarking CPFEM study on dimensionality reduction: quantifying grain-scale strain deviations across crystal structures, morphologies, sizes, and textures

文章 DOI：10.1016/j.ijsolstr.2026.114237

推荐理由：
这篇文章主要研究了一个在 DIC+CPFEM 分析中经常遇到的问题：二维晶体塑性模型到底能够在多大程度上代表真实的三维晶粒变形？

目前很多研究会根据表面 EBSD 数据直接建立二维或准二维 CPFEM 模型，再与 DIC 测得的局部应变场进行比较。但是二维模型缺少厚度方向的晶粒结构以及三维邻域约束，因此即使表面的晶粒形貌和取向完全一致，其局部应变响应也可能与三维模型存在明显差异。

作者针对这一问题，系统比较了 FCC、BCC 和 HCP 三种晶体结构，并进一步考虑了晶粒形貌、晶粒尺寸、初始织构以及边界条件对二维—三维应变偏差的影响。相比单纯比较应变云图，这篇文章还采用 Bland–Altman 方法对对应晶粒进行逐一统计分析，因此对于做 DIC、EBSD、CPFEM以及局部应变场分析的研究比较有参考价值。

作者的理论框架：

作者采用经典的有限变形晶体塑性模型，包括变形梯度乘法分解、基于 Schmid 定律的滑移系分切应力、幂律滑移模型以及晶体塑性硬化模型，并通过 UMAT 在 Abaqus/Standard 中实现。

FCC 考虑 12 个滑移系；BCC 考虑 ({110}\langle111\rangle) 和 24 个滑移系；HCP 则考虑基面、柱面以及锥面 (c+a) 滑移。

作者分别建立了三维 RVE 和准二维 RVE。

三维模型尺寸为 100×100×100 μm，准二维模型尺寸为 100×100×5 μm。二维模型由对应三维模型表面直接截取得到，因此二者具有相同的表面晶粒形貌。

模型采用 C3D8 单元，网格尺寸约为 1 μm，三维模型约包含 100 万个单元，而准二维模型约为 5 万个单元。

作者考虑了三种典型晶粒形貌：

等轴晶、层片晶和柱状晶。

同时设置不同晶粒尺寸和不同初始织构，从而系统分析这些微观因素对于二维模型准确性的影响。

边界条件方面，三维模型采用四周期边界条件（4P），二维模型分别采用四周期边界和普通单轴拉伸边界。

因此：

3D-4P 与 2D-4P 的比较主要反映维度降低本身的影响；

而 3D-4P 与 2D-UT 的比较则同时包含维度降低和边界条件简化的影响。

这种设置可以将二维模型误差的两个主要来源进行一定程度的分离。

作者没有只根据应变云图判断二维和三维模型“像不像”，而是将二维和三维模型中的对应晶粒逐个匹配，并计算每个晶粒的应变差。

随后采用 Bland–Altman 分析方法，通过中位数以及 2.5% 和 97.5% 分位点定义 95% 一致性区间，从而定量评价二维模型相对于三维模型产生的系统偏差以及随机离散。

作者的主要结果：

首先，对于等轴晶随机取向模型，FCC、BCC 和 HCP 在二维化以后都会产生明显的晶粒尺度应变偏差。

其中 FCC 的偏差最小，BCC 次之，HCP 最大。作者给出的结果分别约为：

FCC：30%

BCC：50%

HCP：63%

说明晶体塑性各向异性越强，材料对于三维约束缺失通常越敏感。

另外，二维模型表现出一个比较明显的规律：

低应变晶粒容易被低估，而高应变晶粒容易被高估。

所以即使宏观平均应变差别不大，二维模型仍可能明显放大局部应变不均匀性。

晶粒形貌同样会明显影响二维—三维偏差。

对于 HCP，在周期边界条件下，等轴晶模型的偏差范围约为 1.20，而层片晶下降至 0.92，降低约 23%。这说明几何各向异性与晶体各向异性之间存在明显耦合，不能简单认为某一种晶体结构的二维模型误差始终固定不变。

晶粒尺寸的影响则更加复杂。

作者发现，晶粒尺寸对整体系统偏差影响不大，但会明显改变二维模型误差的离散程度，而且不同晶体结构、不同晶粒形貌得到的规律甚至可能完全相反。

例如等轴晶条件下，随着晶粒尺寸减小，BCC 和 FCC 的偏差有所增加，而 HCP 的偏差反而减小；对于层片晶，BCC 和 FCC 又表现出了相反的尺寸效应。

最后作者研究了初始织构。

结果表明，织构既可能降低二维化误差，也可能放大误差。

例如 FCC 的 Cube 织构可以明显降低二维—三维偏差，而 Copper 和 Brass 织构反而会增加偏差；HCP 的 Basal 织构则能够明显降低二维模型的应变离散。

可以看到，这篇文章虽然没有提出新的晶体塑性本构模型，但是比较系统地回答了一个很实际的问题：

二维 CPFEM 到底在什么情况下能够比较可靠地代表三维模型？

作者的结果表明，二维模型的准确性不仅与晶体结构有关，还受到晶粒形貌、晶粒尺寸、织构以及边界条件的共同影响。

尤其对于局部应变集中、疲劳裂纹萌生和损伤预测等问题，由于真正关心的是少数高应变晶粒，因此直接使用二维 CPFEM 时需要更加谨慎。

来源：我的博士日记

EVPSC-TDT模型文章推荐

文章名称《A crystal plasticity model for hexagonal close packed (HCP) crystals including twinning and de-twinning mechanisms》DOI：10.1016/j.ijplas.2013.02.016 在镁合金、钛合金和锆合金等六方密排结构材料中，塑性变形远比传统立方晶体复杂。除了基面滑移、柱面滑移和锥面滑移外，变形孪生往往承担着不可替代的协调作用。尤其是在循环加载、反向加载和复杂应变路径下，已经形成的孪晶还可能逐渐收缩甚至消失，这一过程被称为去孪生。 2013年，Wang、Wu、Wang和Tomé在 International Journal of Plasticity 上发表论文 A crystal plasticity model for hexagonal close packed crystals including twinning and de-twinning mechanisms，提出了一套具有明确物理图像的孪生—去孪生模型，即TDT模型，并将其嵌入弹黏塑性自洽框架EVPSC中。文章以AZ31B镁合金板材和AZ31挤压棒材为对象，重点研究循环加载和预压缩后反向拉伸过程中的力学响应。许多早期晶体塑性模型主要面向单调加载。当某一孪生系被激活后，通常通过晶粒整体重取向，或者采用经验规则更新孪晶体积分数。这类方法能够近似描述织构演化，却很难处理卸载和反向加载。问题的关键在于：孪生并不是不可逆的“晶粒翻转”。在载荷反向后，孪晶界可以向相反方向迁移，孪晶区域逐渐缩小，母晶区域重新扩展。如果模型只保存孪生后的取向，而没有分别保留母晶和孪晶内部的应力状态，就无法准确描述去孪生，也难以解释镁合金循环滞回曲线中明显的拉压不对称、屈服平台和反向软化现象。TDT模型将整个孪生—去孪生过程划分为四类基本机制：第一类是孪晶形核，即在原始母晶内部引入一个具有较小初始体积分数的孪晶区域；第二类是孪晶长大；第三类是反向加载时的孪晶收缩；第四类是孪晶内部发生的再孪生。文章进一步把孪晶长大分为两个部分：由母晶平均应力驱动的母晶缩减，以及由孪晶平均应力驱动的孪晶扩展。与之对应，去孪生也被拆分为母晶扩展和孪晶缩减。这一划分非常重要。它意味着孪生界面的运动不再只由母晶应力或孪晶应力单独控制，而是同时考虑界面两侧不同的应力状态。母晶和孪晶被视为两个具有独立取向、应力和体积分数的弹黏塑性夹杂，共同与多晶体有效介质相互作用。在EVPSC模型中，每个晶粒都被看作嵌入等效多晶介质中的弹黏塑性夹杂。晶粒总应变率由弹性部分和塑性部分共同组成，因此模型能够保存加载、卸载及反向加载过程中的晶粒应力历史。这正是去孪生计算所需要的基础。论文结果表明，AZ31B板材在压缩过程中发生拉伸孪生，而在反向拉伸初期，去孪生可以在远低于初始孪生应力的宏观应力下启动。这并不一定意味着去孪生具有更低的临界分切应力，而主要源于预加载后母晶和孪晶内部保留的非均匀残余应力。换句话说，EVPSC框架可以从细观应力重新分配的角度，自然解释反向加载中的低屈服现象。模型还揭示了不同阶段的主导变形机制：压缩阶段主要由拉伸孪生承担塑性变形；反向拉伸初期以去孪生为主；当已有孪晶基本消失后，柱面滑移等机制开始承担更多变形；再次反向加载时，孪生重新启动。与此同时，晶体织构也随孪晶体积分数发生可逆或部分可逆的演化。这些结果说明，宏观循环曲线中的平台、拐点和拉压不对称，并不是简单的经验硬化现象，而是滑移、孪生、去孪生和晶粒残余应力共同作用的结果。基于文章的框架，我们可以在evpsc代码的基础上二次开发考虑循环过程中的twin-detwin。使用文章提供的参数与文章进行了相同的循环模拟加载，初始使用典型的基底织构：模拟循环变形结果如下可以看到和文章有良好的对应，（但由于初始未根据文章生成织构文件，而是使用MTEX完成了216个晶粒生成的，存在部分偏差）来源：我的博士日记