使用间隙-偏移边界方法（Gap–SBM）的等几何多片耦合与任意细化和参数化

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：Isogeometric multipatch coupling with arbitrary refinement and parametrization using the Gap–Shifted Boundary Method

　　

本研究提出了一种基于间隙-偏移边界方法（Gap–SBM）的新型等几何多片耦合技术。该方法能够在元素尺寸、多项式次数、方向和参数化任意不同的非协调片之间实现高阶、鲁棒且完全嵌入的耦合。与经典等几何多片策略相比，所提方法避免了需要水密界面或匹配节点向量，并通过在间隙区域上进行积分而不引入额外自由度来保持离散系统的条件数。该耦合通过无惩罚的Nitsche公式施加。该框架还允许在嵌入边界周围插入局部细化的子片，使得能够选择性地引入h-和p-细化，同时保持一致的全局离散化。一组全面的数值实验表明，所提出的多片耦合即使在薄间隙区域或高度非协调片的情况下也能实现最优收敛和线性系统的良好条件数。所提出的方法通过二维设置中的数值示例进行了评估。

**等几何多片耦合研究解读**

**研究背景与问题**
等几何分析（Isogeometric Analysis, IGA）通过采用与计算机辅助设计（Computer-Aided Design, CAD）相同的基函数（如NURBS）统一了几何建模与数值模拟，具备精确几何表示、高阶精度和光滑性优势。然而，工业CAD中复杂几何通常由修剪后的多片NURBS曲面组成，相邻片之间存在参数化不兼容、节点向量不匹配或微小几何间隙，形成非协调（nonconforming）配置。传统的体拟合方法需要生成水密界面或匹配控制点布局，计算成本高昂且难以自动处理复杂装配体。不拟合方法（如有限元细胞法（Finite Cell Method, FCM）和浸入几何分析）虽能绕过网格协调要求，但面临小切割单元导致的系统病态、收敛率下降和不稳定问题。偏移边界方法（Shifted Boundary Method, SBM）通过将物理边界替换为背景网格边界的替代边界并利用泰勒展开施加边界条件，避免了切割单元积分，但其诺伊曼（Neumann）边界条件处理因通量泰勒展开精度有限而损失一阶收敛。为此，间隙-偏移边界方法（Gap–Shifted Boundary Method, Gap–SBM）被提出，它在保留SBM优势的同时，通过在真实边界与替代边界之间的间隙区域内进行精确积分，恢复了诺伊曼条件的最优收敛，且不引入额外自由度。本研究将Gap–SBM创新性地应用于等几何多片耦合，以解决非协调片之间鲁棒、高阶且灵活的耦合需求。

**研究内容与结论**
研究人员提出了一种基于Gap–SBM的等几何多片耦合框架，允许相邻片在网格分辨率、多项式次数和参数化上任意不同，无需水密界面或匹配节点向量。耦合通过无惩罚（penalty-free）Nitsche公式施加，避免调节稳定参数，保持离散系统条件数。框架还支持在嵌入边界周围插入局部细化的子片，实现选择性的h-细化（局部加密）和p-细化（局部升阶）。通过一系列二维数值实验（以泊松问题为模型）验证了该方法的最优收敛性、良好条件数以及局部细化的有效性。论文发表在《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》。

**主要关键技术方法**
1. **Gap–SBM**：将背景网格的替代边界（surrogate boundary）定义为不穿过物理域的网格边界，通过高阶泰勒展开（包含所有混合偏导数）将代理片（patch）的基函数延拓到真实边界与替代边界之间的间隙区域，实现间隙内的积分而不引入额外自由度。
2. **无惩罚Nitsche公式**：采用斜对称（skew-symmetric）变体，设置稳定参数α=0，避免参数调节，保证界面条件的一致性和稳定性。
3. **Coons型几何填充**：将间隙区域划分为线性Coons型三角形或四边形单元，并构造高阶高斯积分规则，确保几何简易且积分精确。
4. **步进收敛分析**：通过逐步增加各参数方向的节点跨度数进行细化，评估误差和条件数随网格尺寸的变化。

**研究结果**

**3.1 多片耦合评估**
通过制造解方法在正方形域上对比单片配置、两片体拟合配置和两片Gap–SBM嵌入配置。所有配置在L²、H¹和L^∞范数下均达到多项式次数p=1,2,3时的最优收敛率（p+1）。Gap–SBM耦合引入的额外误差极小且呈平行偏移，不影响渐近收敛行为，证实了无惩罚Nitsche界面处理是稳定且最优收敛的。

**3.2 耦合方法对比**
在狄利克雷（Dirichlet）边界条件下，比较九片体拟合、两片体拟合和两片Gap–SBM嵌入三种耦合策略。所有策略均达到最优收敛率，但两片体拟合由于裁剪基函数沿界面支撑受限而产生最大误差。Gap–SBM配置在误差常数上接近九片体拟合，同时提供最大的灵活性：界面无需矩形，参数化可独立选择，局部细化无需重划背景网格。

**3.3 诺伊曼边界附近的局部细化**
针对嵌入诺伊曼边界，传统SBM处理损失一阶收敛。通过嵌入局部细化的子片，分别实施h-细化（网格尺寸缩小4倍）和p-细化（内片从二次升至三次），显著降低了局部误差。收敛研究表明，单一单片和两片等分辨率Gap–SBM的误差曲线几乎重合，表明全局误差由偏移诺伊曼边界主导；p-细化配置降低了误差，但未完全恢复最优收敛，这归因于经典SBM对诺伊曼条件的固有限制。

**3.4 条件数分析**
刚阵条件数在Gap–SBM耦合下随网格尺寸h呈O(h^-2)增长，无指数爆炸现象，避免了小切割单元不稳定。局部h-或p-细化仅引起条件数的温和偏移，不影响渐近增长率。对于高次（p=3）基函数，条件数在粗网格下出现平台期，但整体仍保持稳定。

**讨论与结论**
本文基于Gap–SBM提出了一种新型等几何多片耦合框架，实现非协调片之间的鲁棒、高阶耦合。该框架支持任意网格分辨率、多项式次数和参数化差异，无需水密界面，耦合通过无惩罚Nitsche公式施加。数值实验证实了所有范数下的最优收敛，且界面处理不引入显著额外误差。局部细化（h-和p-）有效降低了嵌入诺伊曼边界附近的误差，p-细化尤为高效。条件数分析表明O(h^-2)增长，无病态或不稳定迹象。该方法为工业IGA工作流提供了高度灵活且数值稳健的耦合策略，特别适用于复杂多片壳、裁剪实体以及需要局部细化的场景。研究结论部分翻译如下：
**结论**：本文提出了一种基于Gap–SBM的新型高阶等几何多片框架，能够实现非匹配片之间的鲁棒且精确的耦合。该方法支持网格分辨率、多项式次数和参数化上的任意差异，无需水密界面或匹配控制点分布。耦合通过无惩罚Nitsche公式施加，避免了稳定参数调节，并保持了离散系统的条件数。数值实验证实了该方法在所有考虑范数下均达到最优收敛，且界面处理不引入超出底层离散化所固有的显著误差。该框架还针对嵌入诺伊曼边界进行了局部细化测试，表明带有h-或p-细化的嵌入子片可显著降低局部误差，其中p-细化在改善诺伊曼边界附近精度方面尤为有效，尽管经典SBM的限制阻止了完全最优收敛的恢复，但结果显示了整体求解质量的明显提升。条件数分析进一步验证了方法的稳定性：对于低和中等多项式次数，条件数遵循预期的O(h^-2)趋势，没有病态或不稳定迹象；即使在强局部细化或薄间隙区域的配置下，该方法仍保持数值稳定。插入局部细化片的能力为在复杂边界、局部非线性区域或物理特征（如边界层）附近提高精度提供了灵活机制。该框架也适用于非光滑问题（如塑性或损伤），其中低阶离散化在局部区域有利。在工业环境中，嵌入具有减少基函数支撑的辅助片为缓解交叉串扰提供了有前景的策略。此外，该方法可通过将移动物体嵌入随物体一致移动的片中简化移动浸入体的模拟。其公式自然延伸到三维，其中使用替代边界确保干净的几何交点并避免任意复杂的体积修剪。总之，所提出的Gap–SBM耦合技术为多片等几何分析提供了一个统一且高度灵活的框架，其高阶精度、鲁棒性、几何灵活性和一致条件数的组合使其成为先进工业工作流和复杂多物理应用的引人注目的解决方案。不过，本研究仅限于二维线性椭圆问题，扩展到三维主要涉及额外的算法和实现方面，需进一步研究；此外，还需要针对强非线性材料行为以及更复杂的偏微分方程类别进行有效性评估。