一种用于非界面拟合网格的移位内聚力区方法及其在晶体塑性中的应用

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：A shifted cohesive-zone method for non-interface-fitted meshes with applications to crystal plasticity

　　

移位内聚力区方法（SCZM）通过在一个替代界面上评估界面贡献，实现了在非界面拟合网格上对内聚力界面行为的精确施加。在这项工作中，研究人员将移位断裂框架扩展至内聚力区公式，并展示了其在涉及晶体塑性问题中的适用性。关键进展包括一个定义在替代界面上的简化弱形式、一个在MOOSE框架内与NEML2集成的实现，以及一个降低预处理成本的几何感知、PCA增强的点分类算法。

**论文解读文章**

**研究背景与问题**
在固体与流体力学中，界面条件（如材料界面、裂纹、晶界）对偏微分方程的解至关重要。在晶体塑性中，多晶或复合微结构的代表性体积单元（RVE）模拟中，界面行为（如晶界脱粘、滑移、损伤）常主导宏观响应，影响非弹性变形、表面能释放和失效。然而，传统晶体塑性模型通常假设相邻区域完美结合，只有通过显式界面模型才能捕捉界面软化与损伤。内聚力区公式（CZM）为模拟界面行为提供了物理一致的框架，但其在RVE计算中需要界面拟合网格，这在复杂三维微结构中生成高质量四边形或六面体网格时面临严重挑战（如网格扭曲、负雅可比行列式）。现有非界面拟合方法（如浸入式有限元法IFEM、扩展有限元法XFEM）需修改近似空间，导致非均匀函数空间，而移位边界法（SBM）通过将界面条件移至附近网格实体（替代界面）保持一致空间，但现有移位断裂法（SFM）主要限于单纯形网格，限制了在工程中常用四边形和六面体离散中的应用。因此，迫切需求一种能兼容非界面拟合网格、适用于历史依赖晶体塑性并支持高阶单元类型的移位类内聚力区方法。

**研究内容与结论**
研究人员提出了移位内聚力区方法（SCZM），首次将其应用于历史依赖的晶体塑性问题。通过推导简化弱形式，在替代界面上实施内聚力贡献，并引入方向修正项以处理替代界面法线与真实法线的偏差。在MOOSE框架中实现开源版本，无缝集成NEML2材料库，并开发了基于主成分分析（PCA）增强的几何感知点分类算法，显著加速预处理。通过制造解（MMS）验证确认SCZM在L²范数下呈一阶收敛（因牵引力处理缺少一阶泰勒修正），但在非光滑问题中实际精度与界面拟合网格（IFM）无异。在单界面基准中，SCZM在多种牵引-分离律（指数、双线性混合模式、三维耦合模型）下均准确再现IFM反应力-位移曲线，且方向修正项对消除局部应力畸变至关重要。在二维/三维多晶RVE中，SCZM在粗网格下即能收敛到IFM参考解，准确捕捉全局和局部应力集中、损伤演化及三叉晶界处复杂相互作用。论文发表在《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》。

**关键技术方法**
1. **移位内聚力区弱形式**：将真实界面Γ上的积分映射到替代界面Γ_h（由非界面拟合网格的单元面组成），保留一阶几何修正因子|n·n_h|，并增加方向修正项以补偿法线偏差。
2. **位移跳跃的移位近似**：通过一阶泰勒展开将Γ_h上的位移跳跃修正为真实界面处的值，即?u?_Π = ?u? + ??u?·d，其中d为距离向量。
3. **PCA增强点分类**：利用主成分分析选择射线方向，结合空间搜索结构减少交点检测次数，将复杂度从O(N_elemN_T)降至O(N_elemN_cand)。
4. **MOOSE集成与NEML2耦合**：在MOOSE框架中实现，支持任意物理模块和材料模型，通过NEML2调用历史依赖晶体塑性本构。
5. **距离函数与法线评估**：基于k-d树最近邻查询计算替代界面点到真实界面的投影距离和法线。
6. **界面拟合网格生成与投影**：对SCZM网格进行局部界面拟合手术，生成IFM用于可视化，并通过表面传输扩展进行解投影。

**研究结果**
**4.1 点分类算法性能**：通过对比暴力全交叉扫描，PCA增强方法在边界分段数N_T从100增至1600时保持近似恒定运行时间，而暴力法线性增长；在机翼旋转测试中，PCA方法始终优于固定方向法，对几何方向不敏感。
**4.2 基于制造解的线性各向同性弹性验证**：在平面应变下构造含线性内聚力关系（t_coh = -g）的精确解，SCZM在Delaunay三角网格上呈一阶L²收敛（因P1单元牵引力为零阶外推），但实际精度与IFM无可区分差异；当牵引力梯度为零时（E⁺=E^-），收敛率升至二阶。
**4.3 单界面各向同性弹性基准**：在[0,1]²域界面x=0.5处施加指数型TSL，SCZM反应力-位移曲线与IFM高度吻合，能量释放误差在加载步细化下呈二阶收敛（因路径积分误差主控）；旋转界面测试中，SCZM的von Mises应力沿对角线准确，而缺少SCZM或无方向修正项时出现系统偏移；不同TSL（双线性、三维耦合）下均一致匹配IFM。
**4.4 单界面晶体塑性基准**：应用有限应变单晶粘塑性模型（Voce硬化），倾斜界面，SCZM（含方向修正）的反应力-时间曲线与IFM几乎重合，无方向修正时偏差显著；|t_x|云图显示方向修正项消除了替代界面上的应力畸变。
**4.5 二维多晶RVE**：5晶粒RVE，晶界采用双线性混合模式TSL，SCZM在六层网格细化下反应力-位移和最大损伤演化收敛至IFM；变形场和牵引力幅度沿对角线分布与IFM一致；应力集中区（三叉晶界）的峰值和位置准确再现。
**4.6 三维多晶RVE**：15晶粒三维域，SCZM在非界面拟合六面体网格上准确复制IFM的位移场、反应力-最大损伤曲线；投射到真实界面的损伤场分布与IFM高度一致，仅局部区域因高阶修正缺失出现微小差异。

**总结与结论**
研究表明，SCZM通过替代界面和方向修正，在非界面拟合网格上精确实现内聚力界面行为，避免界面拟合网格生成难题，且对历史依赖晶体塑性具有鲁棒性。实验结论部分翻译如下：
“移位内聚力区方法（SCZM）通过在替代界面上评估界面贡献，实现了在非界面拟合网格上对内聚力界面行为的精确施加。在这项工作中，研究人员将移位断裂框架扩展至内聚力区公式，并展示了其在涉及晶体塑性问题中的适用性。关键进展包括一个定义在替代界面上的简化弱形式、一个在MOOSE框架内与NEML2集成的实现，以及一个降低预处理成本的几何感知、PCA增强的点分类算法。所提出的公式通过基于制造解的研究得到验证，证实了与替代界面上牵引力处理相关的一阶收敛性，而标准的非界面拟合方法未能收敛。在一系列基准问题中，包括各向同性弹性、多种牵引-分离律和历史依赖晶体塑性，SCZM以高保真度再现了界面拟合结果。我们强调，尽管牵引型界面条件在牵引力场中未采用泰勒修正（对比之下位移跳跃经过了泰勒展开），因此该方法在L²范数下仅为一阶精确，但SCZM结果在实际中与界面拟合（IFM）参考几乎不可区分。这是因为所考虑的界面和晶体塑性问题本质上非光滑，即使是符合性有限元也无法达到最优精度，因此SCZM的名义阶数降低几乎没有实际影响。结果还凸显了方向修正项的重要性，它对于准确解析局部应力场并避免在未对准界面附近产生振荡伪影是必不可少的。在二维和三维多晶RVE上的应用表明，SCZM无需界面拟合离散化即可捕捉复杂微结构中的界面响应。未来工作方向包括结合Gap-SBM框架实现高阶一致性，耦合自适应或八叉树离散化提高大规模三维模拟效率，以及扩展到多物理场界面问题（如热-力耦合）和演化界面。”