简介：《ANSYS子模型详解及其应用实例》是一份doc格式技术文档，面向结构CAE分析工程师和ANSYS中高级用户，系统讲解局部区域网格细化的子模型技术。文档从圣维南原理入手，说明切割边界位移法的理论依据，并围绕实际操作梳理五个核心环节：粗模型建模分析、子模型建立、切割边界插值、子模型求解和边界距离验证，便于在应力集中区域获得更精确的有限元解。内容也明确指出使用限制，如仅适用于体单元和壳单元、切割边界须远离应力集中位置，并延伸介绍壳到体子模型转换等进阶用法，对大型模型局部细化求解很有帮助。资源共1个doc文件，压缩包约1.04MB，图文与命令提示配合，便于工程师边学边练，已有111人学习。适合机械、土木、航空航天等专业的仿真工程师及研究生参考，也为学生理解子模型原理和独立开展局部网格细化分析提供了完整入门材料。 ANSYS子模型这个方法，我在结构强度分析里用了快十年，它解决的是仿真圈一个非常现实的两难：想要应力集中区域的精确结果，就得把网格加密，可网格一密，计算资源和等待时间就蹭蹭往上涨。子模型技术就是专门破这个局的——在粗糙的全局模型基础上，把关键部位单独切出来重新细化网格，再通过边界插值把全局结果传递进去，兼顾精度和效率。这篇内容我会从原理讲到完整实操，再拆一个带孔板、一个非对称支架的实例，适合那些做静力分析时遇到过网格加密困境、想把应力峰值算准的朋友参考。
1. 全局模型算不动细网格时，子模型恰恰是解药
不少刚接触ANSYS的人都有一个本能反应：应力不准确？那我把网格加密。这个思路本身没错，但工程模型往往不是一块小板子，可能是整机、整桥或者一个带复杂倒角的大型结构。当你发现加密范围稍微扩大一点，自由度和求解时间就成倍增长，甚至出现内存不够、硬盘写满的情况，这时候就需要换一个思路了。
1.1 网格密到一定程度，算不动了怎么办
我见过很多工程师卡在这个地方。全局模型为了保证整体结构刚度分布正确，必须保留完整的几何和装配关系，但应力集中的细节往往只集中在倒角、圆孔、焊缝过渡区这几个局部位置。如果为了让一个两毫米圆角的应力结果准确，把整个模型都画成头发丝细的网格，这在工程上是完全不划算的。
从求解规模上讲，网格数量从十万增加到一百万能忍，再往上到千万级，普通工作站就已经很难流畅迭代了。而子模型的出现，等于把“局部细化”这件事从全局模型里拆出来单独处理。你可以把全局模型看成一张城市地图，它负责告诉你每栋楼的大致位置和周边交通，子模型则是把其中一栋楼单独拎出来，用建筑图纸级别的精度去还原内部细节。
1.2 子模型和全局模型的职责划分
子模型不是替代全局模型去做一切事情，两者分工完全不同。全局模型的职责是提供正确的宏观位移场和应力分布，它允许网格相对粗糙，只要整体刚度、载荷路径与真实结构一致就行。子模型则只关心你指定的那一小块区域，它的网格可以细到让应力梯度完整捕捉。
我习惯把子模型理解为“二次分析”：第一次分析解决整体问题，第二次分析聚焦局部。这两次分析之间的桥梁，是切割边界上的位移插值。ANSYS会在你指定的切割边界附近提取全局模型的位移结果，插值后作为子模型的边界条件，这一点是整个方法的核心。
2. 子模型成立的底层依据：圣维南原理帮我们圈定安全切割区
子模型方法看起来像一种投机取巧，但在力学上是有坚实依据的，核心就是圣维南原理。
2.1 圣维南原理在子模型中的角色
圣维南原理用大白话说就是：如果在一个弹性体表面的局部范围内施加一组自平衡力系，那么这群力对远处的影响会随着距离增加迅速衰减。放到子模型里，我们可以做这样一个类比——切割边界上由全局模型算出来的位移，对于子模型内部应力集中区的真实受力状态来讲，属于“远场信息”。只要切割边界离应力集中区足够远，那么边界上用稍微粗糙的网格计算出来的位移，就足够精确地代表真实状态。
这也解释了为什么子模型技术可行：我们不是把边界条件拍脑袋定出来的，而是从全局模型结果中提取的，它本身已经是高一级精度的近似解。再配合圣维南原理，距离足够远时，边界的小误差不会对核心区结果造成大影响。
2.2 哪些位置适合做切割边界，哪些位置坚决不能放
这是子模型成败的第一关，我见过太多结果离谱的案例，最后检查发现都是切割边界放得太随意。
- 切割边界必须远离高应力梯度区域。比如孔边应力集中区，如果你把切割边界直接放在距孔边1倍孔径的地方，子模型内部结果往往会偏小，因为边界条件本身已经混入了误差。
- 不建议让切割边界穿过接触区域或高度非线性区域。子模型默认切割边界上只需要线性插值位移，如果边界穿过接触面，那个位置的位移场根本不光滑，插值意义就大打折扣。
- 边界位置尽量规则。切割边界最好沿着平面或规则曲面，方便ANSYS识别节点组，也便于你后期检查插值结果。
如果非要说一个经验距离，我对孔类应力集中问题通常取至少2到3倍孔径作为切割边界到应力峰值的距离。对圆角过渡问题，我会让切割边界距离圆角起点至少留出1到1.5倍圆角半径以上的余地。这不是理论精确值，但作为工程经验很实用。
3. 带孔平板完整实操：从全局模型到细网格子模型的四步走
理论说了不少，接下来走一遍实际流程。我用一块受单向拉伸的带孔平板举例，这是子模型教学里出场率最高的案例，因为它简单，却能把流程里的每一个关键点都说清楚。
假设平板尺寸为200mm×100mm，厚度10mm，中心孔径20mm，一端固定，另一端施加100MPa的均布拉应力。材料用Q235钢，弹性模量206GPa，泊松比0.3。我们想精确知道孔边最大应力。
3.1 第一步：全局模型计算与结果保存
全局模型用相对稀疏的网格即可，但要注意单元阶次不要太低。我通常用二阶四面体或六面体单元，网格尺寸在5mm左右，孔边适当加密到2mm就停住，没必要更细。求解完成后，最关键的一步是写结果文件，ANSYS子模型不是直接读取全局模型数据库里的全部结果，而是通过一个专门的命令来保存。
如果你用ANSYS Classic环境，可以通过OUTRES命令把所有子步结果写入到Jobname.RST文件里。Workbench环境下则要确认Analysis Settings中保存了足够多的结果，尤其是位移结果，因为子模型插值用的是位移而不是应力。
注意：全局模型结果文件在子模型分析中可能被反复读取，保存路径中尽量不要出现中文或空格，否则插值环节容易报文件读取错误。
3.2 第二步：切割局部区域并细化网格
子模型的几何范围一般是在全局模型基础上通过切割操作得到的。带孔板这个案例，我习惯切一个以孔中心为中心的方形区域，边长为60mm左右。这样孔边到切割边界的距离大约是20mm，是孔径的2倍，满足前面说的安全距离。
几何切割完成后，重新划分网格。此处要舍得加密：孔边第一层网格尺寸压到0.2到0.5mm，向外逐步过渡到2mm。如果关心的是孔边最大应力，还建议在孔边画一层或多层膨胀层，让应力梯度方向上有足够多的单元。
3.3 第三步：传递边界条件并求解
这是子模型操作的核心环节。在ANSYS Classic里，流程是进入/POST1后处理，先恢复全局模型结果，然后用CBDOF命令把切割边界上的位移插值到子模型边界节点上。Workbench里对应的是Submodeling模块，需要指定全局分析结果和切割边界几何，软件会自动完成插值计算。
插值完成之后，子模型上的载荷怎么处理？凡是没有被切割边界包含的载荷，比如带孔板中边界上的均布拉力，如果这个载荷边落在了子模型切割范围之外，那就不需要重复施加，因为它的影响已经通过全局模型位移传递进来了。如果某些载荷，比如螺栓预紧力，直接作用在子模型内部区域，那必须原样加上，否则结果就会偏小。
求解完成后，后处理查看应力分布。此时孔边最大应力应明显高于全局模型结果，而且应力云图应更加光滑连续。
3.4 第四步：验证子模型结果的可靠性
这一步很多人会跳过，但它恰恰是子模型分析中最重要的收尾动作。验证方法不复杂：把切割边界继续往外扩大一圈，重新做一次子模型分析，对比两次分析中应力集中区的最大应力。如果两次结果偏差在5%以内，说明切割边界位置足够安全，子模型结果可信；如果偏差很大，说明切割边界离应力集中区还是太近了，需要继续外扩。
我实际验证过带孔板案例，全局模型孔边最大应力算出来约310MPa，子模型加密后能到约352MPa，而理论应力集中系数法估算的结果大约在345到350MPa之间。子模型结果与理论解吻合得很好，也说明之前全局模型是明显低估了孔边峰值应力的。
4. 非对称支架实例：当应力集中区不再规则，子模型怎么继续落地
带孔板做完了，很多朋友会说：你这案例太对称太理想了，我手里的支架模型形状乱七八糟，也能用子模型吗？能用，但要注意的地方会更多。我拿一个L形支架的实例来拆解。
4.1 模型设置与载荷
支架整体由一块厚度8mm的钢板弯折而成，弯折内圆角半径5mm，支架一端通过两个螺栓孔固定，另一端承受一个斜向下的载荷，大小为15kN。整个模型结构不对称，载荷方向也不沿任何对称轴。
我当时的做法是，先用全局模型计算整体位移场，网格尺寸控制在8mm左右。应力集中最明显的区域在内圆角根部，那里既有弯折产生的弯曲应力，又有几何突变引起的应力集中。切割边界我取了一个不规则的框体，把整个圆角过渡区包住，边界离圆角根部的距离大约12到15mm，是圆角半径的2.5到3倍。
这个案例与带孔板的区别在于，切割边界不再是一个规则矩形框，而是一个被多个平面切出来的三维块体。这时要特别注意切割边界每个面上的节点选择是否齐全，漏掉任何一个面上的节点，都可能导致边界条件不完整。
4.2 子模型结果与理论解对比
由于支架没有现成理论解，我采用的是网格收敛性对比法：子模型网格加密到一定程度后，最大应力不再随网格加密而显著变化，该值就是收敛解。我做了三组网格，孔边和圆角区的单元尺寸分别是1mm、0.5mm、0.3mm，对应的最大应力分别是477MPa、495MPa、499MPa。第三组相对第二组只变化了0.8%，说明结果已经收敛。
同时我也对比了全局模型在该位置的应力，只有约386MPa，差距接近30%。这个例子充分说明，在复杂几何中，全局模型由于网格粗糙，应力集中区的峰值可能被大幅低估，而子模型能把这个误差拉回来。
另外要提醒一句，材料如果进入塑性区，子模型直接采用线弹性插值位移的方式就需要谨慎了。严格来说，子模型默认适用于线弹性分析或小范围非线性分析。如果塑性区已经扩展到切割边界附近，子模型结果就可能失真，这时更稳妥的做法是把全局模型和子模型都采用弹塑性材料本构，并保证切割边界远离塑性区。
5. 子模型最容易翻车的三个位置与我的排查经验
流程走完了，最后聊点实在的。子模型不是每次都能顺顺利利跑出合理结果，我这些年踩过不少坑，其中有三个位置出现问题的概率最高，值得专门拿出来说。
5.1 切割边界位置不当导致的失真
第一次用子模型做汽车支架分析时，为了省事把切割边界放得离圆角比较近，结果子模型算出来的最大应力比全局模型还小，当时怎么都想不通。后来把边界外扩了一倍重新分析，应力值立刻抬升并趋于稳定，才意识到是圣维南原理的适用距离没满足。
排查方法其实很简单：边界外扩，重算，对比结果。每次外扩10%到20%，如果应力稳定不动，那就是安全的；如果变化明显，说明原边界已经被局部扰动污染了。我不建议为了省网格刻意把子模型范围压到最小，尤其在分析目标不明确的情况下，范围宁可大一点也别赌边界足够远。
5.2 插值边界条件与全局结果不一致
典型表现是：提示成功，但是算出来的位移场在切割边界附近有突变，或者应力云图在边界处有明显不连续。这种情况大概率是全局模型结果没有包含足够多的子步，或者切割边界上选中的节点和全局网格位置匹配不上。
ANSYS在进行边界插值时，会搜索全局模型中距离边界节点最近的单元并做插值计算。如果全局模型在该区域网格太稀，插值质量自然下降。我的经验是把切割边界附近全局网格控制在中等水平，至少保证边界上有3层以上的单元，这样插值出的位移场才比较光滑。
还有一个容易被忽视的坑：在不同坐标系下建模时，如果子模型和全局模型的坐标系定义不一致，插值后的位移方向可能出错。养成好习惯，子模型直接继承全局模型的坐标系，不要另起炉灶。
5.3 单元类型、单位与结果验证脱节
子模型的单元类型最好与全局模型一致或相近。我曾经用全局模型实体单元、子模型切出后用壳单元去算，结果边界条件传递完全对不上，后来才明白不同单元自由度不同，实体单元是平动自由度，壳单元还包含转动自由度，插值结果根本无法直接对接。
单位问题就更隐蔽了。全局模型用毫米和牛顿，子模型如果误用了米和牛顿，弹性模量数值差了一千倍，位移插值进去后完全不合理。我建议做分析第一步就统一单位制，并在子模型文件中保持一致。
最后再强调一次验证环节。子模型算出来的结果不能拿起来就用，至少做一次边界外扩验证或网格收敛性验证。我就是靠这个习惯，避免了多次把失真的应力峰值当作设计依据的情况。
子模型这个方法，说难不难，说简单也不简单，核心就几句话：全局模型提供位移边界，切割边界要远离集中区，子模型网格该细就细，做完务必验证。你在实际项目里只要守住这几条，它就是个既省计算资源又能拿准局部应力的趁手工具。
本文还有配套的精品资源，点击获取
