4. 切割边界位移插值机制
子模型法能否成立，关键不在子模型本身算得多细，而在于切割边界上的位移是否被“正确搬运”。粗模型在切割边界上只有稀疏的节点位移，子模型在该边界上却往往有更密的节点，两者网格不重合，因此必须通过插值把粗模型的位移场映射到子模型的边界节点上。理解这套插值机制，才能判断子模型结果是否可信。
一、位移边界条件的插值算法
在 Ansys Mechanical 中，子模型边界位移通常通过 Submodeling 边界条件（或 Nodal Displacement 配合外部插值）施加。其核心是：先由粗模型求解得到切割边界节点的位移结果，再以这些离散位移为“采样点”，在子模型边界节点位置重建位移场。Mechanical 内部常用的插值方式包括形函数插值与距离加权插值两类。
形函数插值（基于单元插值）：把粗模型切割边界所在的单元作为插值载体。对子模型边界上的任一目标节点 P，先判断它落在粗模型哪个边界单元内，再用该单元的形函数 Ni 组合单元节点位移：
u(P) = Σ N_i(ξ,η,ζ) · u_i
v(P) = Σ N_i(ξ,η,ζ) · v_i
w(P) = Σ N_i(ξ,η,ζ) · w_i其中 (ξ,η,ζ) 是 P 在该单元局部坐标下的参数坐标，ui、vi、wi 为粗模型单元节点在全局坐标系下的三个位移分量。这种方法精度高，能保持位移场在单元内的连续性，是首选方式。
距离加权插值（基于最近节点）：当目标节点落在粗模型单元之外，或单元形状畸变导致参数坐标求解失败时，Mechanical 会退化为按距离加权的方式，用最近的若干粗模型节点位移加权平均：
u(P) = Σ w_j · u_j / Σ w_j , w_j = 1 / d_j^pdj 为目标点到第 j 个粗模型节点的距离，p 为权重指数（常见取 2）。这种方式鲁棒但精度较低，尤其在位移梯度大的区域容易产生平滑误差。
Shape Function 即为形函数插值，若为 Distance Based 则为距离加权。也可通过求解信息（Solver Output）中的映射警告确认。二、插值精度影响因素
插值误差并非只由算法决定，更多来自“采样质量”和“几何匹配”。主要影响因素如下：
- 粗模型切割边界网格密度：采样点越稀，位移场重建越粗糙。若粗模型边界单元尺寸远大于子模型，插值相当于用低阶多项式拟合高梯度场，误差显著。
- 切割边界位置：应选在位移梯度平缓、远离应力集中和圣维南区的截面。若边界穿过高梯度区，插值误差会被直接带入子模型。
- 粗模型与子模型边界几何一致性：切割面必须几何重合。若子模型边界比粗模型切割面外扩或内缩，落在单元外的节点会触发距离加权，精度下降。
- 粗模型单元质量：严重畸变、负雅可比或高长宽比的边界单元，会使参数坐标求解不稳定，形函数插值退化。
- 位移分量与坐标系：插值在全局坐标系下进行。若子模型使用了局部坐标系或柱坐标系，需确认边界位移的参照系一致，否则会出现方向性错误。
- 插值阶次与单元类型：粗模型若为二阶单元（如 SOLID186），形函数为二次，插值精度高于一阶单元（SOLID185）。
三、节点映射检查方法
插值是否成功，不能只看求解是否收敛，必须做映射检查。推荐按以下顺序核查：
- 查看映射警告：求解后在 Solver Output 中搜索
not mapped、outside、extrapolat等关键字。若出现大量未映射节点，说明边界几何或容差有问题。 - 可视化边界位移：在子模型结果中插入
User Defined Result，表达式用U或USUM，只显示切割边界节点，检查位移是否连续、无突变或锯齿。 - 对比粗模型与子模型边界位移：在相同位置提取粗模型和子模型的位移值，逐点比较。若差异超过 5%（示例判据，需按项目精度要求调整），应加密粗模型边界网格或调整切割位置。
- 检查节点映射数量：在 Submodeling 边界条件的 Worksheet 或 Solution Information 中查看映射节点数与总边界节点数之比，理想应接近 100%。
- 做网格收敛验证：将粗模型边界网格加密一档重算，若子模型边界位移变化明显，说明原插值精度不足。
四、插值机制与检查流程示意
下图给出从粗模型位移采样、插值映射到子模型边界节点，再到映射检查判断的完整关系。图中同时标出形函数插值与距离加权插值两条路径，以及精度影响因素的作用位置。
从图中可以看到，插值算法位于粗模型与子模型之间，是位移传递的“咽喉”。精度影响因素作用于插值环节，而映射检查则是对插值结果的反向验证。若检查不通过，应回到粗模型加密边界网格或调整切割位置，而不是直接修改子模型内部网格。
五、实操建议与判据
结合上述机制，给出可执行的检查清单（示例判据，需按项目精度要求调整）：
| 检查项 | 操作位置 | 合格判据（示例） | 不合格处理 |
|---|---|---|---|
| 映射警告 | Solver Output | 无 not mapped / outside 记录 | 检查切割面几何与容差 |
| 映射节点比例 | Submodeling 详情 | ≥ 99% | 加密粗模型边界网格 |
| 边界位移连续性 | User Defined Result (USUM) | 无突变、无锯齿 | 改用形函数插值或调整边界 |
| 粗/子模型位移差异 | 同位置节点提取 | ≤ 5% | 加密粗模型或移动切割位置 |
| 网格收敛 | 粗模型加密重算 | 边界位移变化 ≤ 2% | 继续加密直至收敛 |
需要强调的是，插值精度与子模型内部网格密度无关。子模型内部再细，也只能在给定边界位移下求解；边界位移错了，内部结果就失去意义。因此，切割边界应选在位移梯度平缓、远离应力集中的截面，并保证粗模型在该处有足够密的网格。对于高梯度区域，宁可把切割边界外移，也不要让边界穿过应力集中区。
最后，坐标系一致性必须全文统一。本文及图示中位移分量均以全局坐标系为参照；若子模型使用局部坐标系，应在施加边界条件前明确转换关系，并在检查时确认方向无误。只有插值算法正确、影响因素受控、映射检查通过，子模型结果才具备工程可信度。