剪切带相遇后会怎样？聚合物管内压膨胀的局部变形机制

文章题目：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 时需要更加谨慎。来源：我的博士日记