1. 子模型方法总览

在大型装配体或复杂零部件的有限元分析中，工程师经常面临一个两难：整体模型必须足够大才能反映真实的载荷传递路径与边界约束，但整体网格又不可能处处细化到能捕捉圆角、螺纹、焊缝或接触区的应力集中。若强行在整体模型上全局加密，网格规模、求解时间与内存需求会迅速失控；若只做粗糙整体分析，局部峰值应力又会被严重低估。子模型技术（Submodeling）正是为解决这一矛盾而生的层级式分析方法：先用较粗的整体模型获得全局位移场，再把该位移场作为边界条件"驱动"一个从整体模型中切割出来的局部精细模型，从而以可控的计算代价获得局部高精度结果。

基本原理：圣维南原理支撑的"切割—驱动"逻辑

子模型方法成立的理论基础是圣维南原理（Saint-Venant's Principle）：当一组自平衡力系作用于弹性体的一小部分区域时，其影响在远离作用区的部位迅速衰减。这意味着，只要局部子模型的切割边界距离应力集中区足够远，边界上的位移分布就主要由整体结构的整体变形决定，而对局部细节（如圆角半径、倒角形式）不敏感。因此，把整体模型在切割边界上的位移解插值传递给子模型，就能在局部区域复现真实的变形环境。

由此形成两个模型的明确分工：整体模型（Global/Coarse Model）负责承载真实约束、真实载荷与全局刚度分布，输出的是位移场；子模型（Local/Submodel）只保留局部几何细节，其边界条件不再是原始约束，而是从整体模型插值得到的位移。子模型不关心载荷如何施加，只关心边界如何变形——这是理解整个技术的关键。

层级关系与实现路径

在 Ansys Mechanical 中，子模型分析通常以"两个独立分析系统"的形式组织：一个静力学分析系统承载整体模型，另一个静力学分析系统承载子模型，二者通过数据链接传递位移。整体模型求解完成后，在子模型的切割边界面上插入 Imported Load（导入载荷）对象，将类型设为 Displacement（位移），并指定数据来源为整体分析的结果文件。Mechanical 会自动完成边界节点上的位移插值，无需手工编写插值程序。

实现过程中有几个必须遵守的工程约定。第一，切割边界应选在应力梯度平缓、远离关注区的截面，一般建议距离应力集中区至少 2～3 倍壁厚或特征尺寸。第二，子模型在切割边界上的网格可以比整体模型更细，但边界几何必须与整体模型对应位置一致，否则插值会引入几何误差。第三，子模型内部不得再施加与整体模型重复的约束，否则会与导入位移冲突，导致求解失败或结果失真。第四，若整体模型使用了实体单元而子模型使用壳单元，或二者单元阶次不同，插值精度会下降，应尽量保持单元类型与阶次一致。

典型工程价值与适用场景

子模型方法的价值集中体现在"用两级网格换取全局与局部兼得"这一策略上。以某支架焊接结构为例：整体模型包含多个板件与焊缝连接，网格尺寸约 8 mm，单元数约 30 万；关注区为焊缝趾部，实际需要 0.5 mm 级网格才能捕捉应力梯度。若整体加密到 0.5 mm，单元数将膨胀到数千万级，普通工作站难以承受。采用子模型后，整体模型仍以 8 mm 网格求解，子模型只截取焊缝附近 120 mm × 80 mm × 40 mm 的区域，以 0.5 mm 网格离散，单元数约 40 万，总计算量下降一个数量级以上，而焊缝趾部的峰值应力与热点应力分布则能真实反映。

子模型特别适合以下场景：一是应力集中区远离载荷与约束施加位置，满足圣维南原理的适用条件；二是局部几何细节（圆角、倒角、孔、螺纹、焊缝）在整体模型中已被简化甚至忽略，需要在子模型中恢复；三是需要做疲劳寿命评估，必须获得热点应力或缺口应力；四是接触非线性只发生在局部，整体模型用绑定接触近似、子模型用真实接触定义。反之，若关注区本身就是载荷或约束施加区，或整体模型已经足够精细，则子模型带来的收益有限。

与整体模型的关键差异

需要强调的是，子模型结果的可信度高度依赖切割边界的位置与插值质量。工程实践中应做边界敏感性检查：把切割边界向外推移一定距离重新计算，若局部峰值应力变化在 5% 以内，说明边界位置合理；若变化显著，则需继续外移。此外，子模型与整体模型的单元类型、材料本构应保持一致，避免因刚度定义差异引入虚假的位移不匹配。掌握这些要点后，子模型就能成为复杂结构局部精细化分析中最稳定、最经济的手段之一。