等几何间隙平移边界方法（Gap–SBM）

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：The isogeometric Gap-Shifted Boundary Method

　　

在二维和三维中的高阶等几何间隙平移边界方法（Gap–SBM）公式。对于浸入Dirichlet和Neumann边界条件的最优收敛性。没有额外的自由度且最优条件数缩放。计算机辅助设计（CAD）和立体光刻（STL）衍生几何的自动分析。

本研究针对等几何分析（IGA）中复杂几何的网格生成困难及浸入方法的小切割单元问题，提出了一种高阶等几何间隙平移边界方法（Gap–SBM）。研究背景在于，传统等几何分析虽能精确表示几何，但构造符合复杂边界的网格耗时且自动化困难；浸入方法（如有限细胞法FCM）虽解耦网格与几何，但小切割单元导致病态系统矩阵和收敛率下降。平移边界方法（SBM）通过替换为替代边界解决了小切割单元问题，但对Neumann边界条件收敛次优。为此，研究人员引入Gap–SBM，通过显式重构和积分真实边界与替代边界之间的间隙区域，恢复最优精度。论文发表在《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》。研究人员得出以下主要结论：Gap–SBM在二维和三维线性弹性问题中，对浸入Dirichlet和Neumann边界条件均达到最优收敛率，且条件数呈最优的h^?2缩放，无额外自由度。该方法对替代边界位置敏感性低，能自动处理CAD和STL衍生几何，计算开销可控。重要意义在于，为等几何分析提供了一种鲁棒、精确且条件数良好的浸入框架，避免了传统浸入方法的病态问题，并拓展了高阶格式在复杂几何中的应用。

研究人员为开展研究用到的主要关键技术方法包括：①二维间隙区域采用Coons patch参数化构建曲线四边形单元，基于层次二分法进行内部细分以处理小尺度特征；②三维间隙区域采用局部金字塔-四面体拓扑分解，包括Type-1（面连接）、Type-2（边连接）和Type-3（顶点连接）单元，并引入曲线重建仅当雅可比行列式正定；③边界条件处理策略包括直接施加、基于层次二分法的插值转移以及增强SBM外推（使用完整多元泰勒展开）；④采用无罚参数的Nitsche方法和对称或斜对称变体，结合间断伽辽金框架处理间隙单元界面不连续性；⑤数值积分采用高斯求积，曲线单元使用张量积求积。

**4.1 精确几何上的方法比较**
通过菱形孔等精确几何测试，Gap–SBM对Dirichlet边界条件在所有范数中达到最优收敛率，与体适形参考解一致；对Neumann边界条件，标准SBM出现次优收敛和强振荡，而Gap–SBM恢复最优率且对替代边界位置高度鲁棒。

**4.2 曲线嵌入几何上的边界条件处理比较**
在同心圆环测试中，比较了直接施加、插值转移和增强SBM外推三种策略。对于p=1，插值与SBM外推显著提高精度；p=2时几何近似收敛快，三种策略差异可忽略；p=3时，SBM外推和插值在细网格上提升精度，但插值在极细网格上因固定插值点数量而收敛减缓。

**4.3 内部细分对复杂嵌入几何的影响**
在包含小尺度特征的复杂几何中，即使背景网格远大于特征尺寸，结合层次内部细分（最多100层）和增强SBM外推，所有边界条件处理策略均恢复最优收敛率，误差接近体适形解。细分后几何误差指标（面积误差和最大边界距离）显著下降。

**4.4 计算成本比较**
在圆形基准测试中，Gap–SBM的预处理成本近似与DOF^0.5成正比，组装成本线性增长，求解成本超线性增长，额外开销在细网格下可忽略。在复杂几何中，Gap–SBM（无/有细分）的组装和求解时间接近标准SBM，且随细化趋于一致。

**4.5 条件数的影响**
在圆形几何中，未稳定化IBRA条件数出现严重切割位置依赖振荡（峰值达10³⁶），而SBM和Gap–SBM均呈最优(1/h)²增长，且Gap–SBM的振荡更小，常数更低。

**4.6 三维数值实验**
4.6.1 球体基准：精确NURBS球体测试中，Gap–SBM对p=1,2,3均达到相对L²范数p+1和H¹范数p的收敛率，条件数呈h^?2缩放且常数低于标准SBM。峰值内存随曲线重建增加（p=2从2.9 GB到6.3 GB，p=3从8.4 GB到18.8 GB）。
4.6.2 斯坦福兔子STL-to-分析工作流：线性间隙重建下，p=1和p=2的Gap–SBM在细网格上恢复预期收敛率，标准SBM次优。几何构造额外开销与标准SBM同量级，组装时间线性增长。p=1解在细网格达到与体适形有限元解（约1.1×10⁶自由度）相当的精度。

**讨论与结论**
讨论部分指出，Gap–SBM通过显式间隙积分和增强SBM外推，在二维和三维线性弹性问题中实现了最优收敛和良好条件数，且无额外自由度。主要局限性在于三维高阶曲线间隙单元需大量求积点，导致内存和组装成本增加。未来工作将聚焦于更紧凑的求积规则（如求积压缩和矩拟合）以及并行预处理和组装。研究结论部分翻译如下：本文提出了Gap–SBM的高阶扩展，旨在结合浸入和平移边界方法的灵活性与高阶离散所需的精度、鲁棒性和数值特性。该方法解决了二维和三维线性弹性问题，通过Coons patch参数化（二维）和局部金字塔-四面体分解（三维）实现间隙积分。数值实验表明，Gap–SBM对Dirichlet和Neumann边界条件均恢复最优收敛率，对替代边界位置敏感性低，条件数呈最优(1/h)²缩放且无自由度增加。方法适用于CAD和STL衍生几何，计算开销可控。当前主要限制是三维曲线间隙单元的高求积点开销，未来将改进求积策略并扩展到非线性问题。