企业官网建设流程全解析

简介：《ANSYS子模型介绍及应用实例》是一份面向有限元分析工程师与科研人员的专业文档，主要解决局部网格细化求解与计算效率之间的平衡问题。当整体模型在应力集中区域网格密度不足时，不必重新划分整个模型，而是可以通过切割边界位移法对关注区域单独建立子模型，从而获得更精确的局部结果。文档系统介绍了子模型的理论基础，即圣维南原理及其适用前提，梳理了生成粗糙模型、创建子模型、切割边界插值、运行子模型分析、验证切割边界位置是否满足要求等关键步骤，同时说明了该方法只适用于体单元和壳单元等限制条件，并补充了壳到体子模型这一特殊应用技巧。资源以单个doc格式文档打包，总大小仅1.04MB，内容精炼而完整，适合需要学习ANSYS子模型操作流程、解决应力集中与网格细化工程问题的读者，目前已有111人学习使用。 做结构强度分析的人，十有八九都被同一个问题折磨过：算带孔板、齿轮齿根或者圆角过渡处的应力，网格粗一点，最大应力230 MPa；加密一圈，变成260 MPa；再加密，数值还在往上涨。全模型加密吧，单元数量涨得吓人，一台工作站算一个晚上都跑不完。我第一次碰到这个现象时，还以为是约束加错了，后来才知道这是应力集中区的典型网格敏感性问题。解决这个问题的经典手段之一，就是ANSYS里的子模型（Submodeling）技术。这篇文章我会把子模型的原理、适用边界、完整操作流程，以及一个能直接照着复现的带孔板拉伸应用实例全部讲透，适合正在做结构强度校核的工程师，刚接触有限元但想少走弯路的新手也能跟着操作。

1. 子模型到底是什么：从网格加密的“无底洞”说起

1.1 应力集中区的网格敏感性问题

应力集中不是软件算出来的“坏结果”，而是真实的物理现象。孔、缺口、圆角、焊缝趾部这些位置，力流被强行挤压，局部应力会远高于名义应力，通常用理论应力集中系数Kt来衡量。问题出在有限元计算上：峰值应力所在的区域应力梯度极大，网格太粗根本捕捉不到，算出来的最大应力会明显偏小，而且网格越粗偏得越离谱；加密网格后峰值逐渐上升，直到网格细到一定程度才会收敛。

我早年做齿轮齿根弯曲强度校核时就吃过这个亏：齿根圆角处用2 mm网格算出来最大应力320 MPa，加密到0.5 mm后变成了380 MPa，差了快20%。当时第一反应是模型有问题，反复检查边界条件都找不到原因。后来才明白，齿根区域的应力梯度摆在那里，网格密度决定了解的分辨率，粗网格只能给出一个“被平均过的低值”。

这不是软件bug，也不是操作失误，而是应力集中区共有的数值特性。尖角、裂纹尖端甚至会出现理论上的应力奇异，应力随网格细化不断增大不收敛，那就需要断裂力学专门处理。但工程上大量的圆角、孔边问题，应力峰值是收敛的，只是需要的网格尺度很小。如果整个结构都用这么细的网格去算，单元规模会直接爆炸，计算资源和时间成本完全不可接受。

1.2 子模型的核心思想和适用范围

子模型做的事情其实很朴素：先用较粗的网格对整个结构做一次全局分析，这一步只要求“大范围趋势对”；然后把关心的高应力局部区域单独“切”出来，用非常细的网格重建一个局部模型；最后把全局分析求出来的位移，作为边界条件施加到这个局部模型的切割边界上，重新求解一次，得到精确的局部应力。

这个思路的理论基础是圣维南原理：当我关心的是局部区域内部的应力时，只要切割边界离这个区域足够远，边界上具体载荷怎么分布影响并不大，真正起作用的是边界上的合力与合力矩。全局模型虽然网格粗，但在切割边界处的位移场是相对可靠的，拿它来驱动局部细网格模型完全站得住脚。

所以子模型特别适合三类场景：一是孔、圆角、缺口这些应力集中位置的高精度应力提取；二是大型结构里某一小块区域需要做疲劳、强度细化评估时，不需要把整个模型都建细；三是断裂分析中配合裂纹尖端网格做局部精细化。但它也有明显边界，如果结构发生大范围塑性屈服，或者局部区域的接触状态会显著影响整体刚度，又或者分析对象本身就是薄壁结构的局部屈曲问题，这时候子模型的圣维南原理前提就不成立了，硬用会得出错误结论。

需要提醒的是，子模型和另一个常见概念——子结构——不是一回事。子结构是把一组单元凝聚成超单元，反复求解时提高效率；子模型则是局部重新建模加密，通过切割边界位移从全局结果中提取边界条件。两者一个“缩”一个“切”，解决的问题完全不同。

2. 子模型技术原理解析：切割边界位移是怎么“搬”过去的

2.1 为什么ANSYS选择“位移”而不是“力”做边界

子模型切割边界上有两种“搬运数据”的思路：一种是从全局结果提取切割边界上的位移，把它作为位移约束施加到子模型上；另一种是提取节点力或面力，把它作为载荷施加到子模型上。ANSYS经典界面里的子模型流程默认走的是位移路径，这个选择背后是有道理的。

位移边界受圣维南原理保护。切割边界距离关注区域足够远时，边界上位移的局部细节误差会随距离迅速衰减，不会污染内部的应力结果。而力的路径对边界网格质量极其敏感，全局网格粗，切割边界上的力分布本身就不平滑，转成载荷以后很容易在边界上引起人为的应力扰动。简单说，位移是“积分量”，比力更平滑、更稳定。

另外从数值实现角度看，节点位移插值非常成熟，子模型切割边界节点不需要和全局模型节点位置重合，ANSYS会自动用周围全局单元的结果做插值，计算效率高且稳定。如果去提取力，还要考虑等效节点力转换、面积积分等一系列麻烦事，误差来源更多。所以我个人在做子模型时，从来都是老老实实用位移边界，这也是官方推荐的默认做法。

2.2 子模型操作的三段式流程

子模型在实际操作中就是一个标准的三段式流程。第一步是全局模型求解，这一步要求粗网格趋势正确，完整施加真实结构的载荷与约束，求解后保留结果文件。第二步是建立子模型，把需要细化的局部几何切出来重新划分网格，施加除了切割边界位移之外的其余边界条件，比如对称约束、固定约束等。

第三步是关键，进入求解器后执行CBDOF命令，让ANSYS从全局结果文件中读取位移场，插值到子模型的切割边界节点上，自动生成位移约束，然后求解。需要特别说明的是，ANSYS约定要把切割边界上的所有节点定义成一个名为CUT的组件，CBDOF命令默认就是从CUT组件包含的节点上读取并施加位移。这个命名习惯我一直沿用，省得每次翻帮助文档。

命令流骨架大致是下面这样：

! 全局模型分析完成后保存结果 /SOLU SOLVE FINISH ! ===== 建立子模型，以下为前处理阶段 ===== /PREP7 ! ... 创建局部几何、划分细网格 ... ! 选择切割边界节点并定义为CUT组件 NSEL,S,LOC,X,0.04 ! 示例：右侧边界 NSEL,A,LOC,Y,0.04 ! 示例：上侧边界 CM,CUT,NODE ALLSEL,ALL ! ===== 进入子模型求解 ===== /SOLU ! 施加局部模型自身的其他约束，例如对称约束 ! 从全局结果文件中读取CUT组件节点位移作为边界条件 CBDOF SOLVE FINISH

这个流程看起来简单，但有几个细节必须重视。全局模型和子模型必须采用同一套单位制，如果全局模型长度用米、力用牛顿，子模型长度用毫米、力用牛顿，读入的位移边界就完全对不上。还要保证结果文件一致，CBDOF读取的是当前工作目录下的结果文件，中间切来切去时不要搞混Jobname。CBDOF命令的具体菜单位置不同版本略有差异，以帮助文件为准。

2.3 经典界面和Workbench怎么选

子模型技术最早是在ANSYS经典界面（APDL）里发展起来的，命令流很透明，每一步干了什么清清楚楚。直到今天，很多老工程师做子模型还是习惯在经典界面里操作，就是因为流程固定、可控、不容易出幺蛾子。

Workbench里也能做子模型，但说实话，操作路径在不同版本之间差异比较大，有的版本可以通过导入切割边界位移实现，有的版本需要配合External Data模块，用起来并不算顺滑。我的个人习惯是：几何建模用SpaceClaim或Workbench DM，网格划分也先在Workbench里控制好，最后把模型导出成CDB文件，到经典界面里统一执行子模型命令流。这样既能享受Workbench的建模便利，又能保住APDL对子模型流程的精细控制，还能方便积累一套可复用的命令流模板。

3. 应用实例：带孔板拉伸应力集中子模型全程复现

3.1 模型参数与分析目标

为了把子模型的操作讲透，我选一个经典到不能再经典的例子：中心带圆孔方形薄板承受单向拉伸。这个例子在机械设计手册里能查到理论应力集中系数，方便验证子模型结果的正确性，不会出现“算完也不知道对不对”的情况。

模型参数如下：

理论最大应力大约就是100×2.8=280 MPa，后面所有数值都拿这个做基准。为了减小计算规模，利用结构和载荷的对称性，取四分之一模型计算，在对称面上施加对称约束。这个操作本身也是工程中节省计算资源的重要习惯。

3.2 全局模型：先求“趋势对”

全局模型用六面体单元，单元尺寸取5 mm，厚度方向分两层。四分之一模型尺寸是100 mm×100 mm×5 mm，算下来单元不到一千个，求解只用几秒钟。边界条件上，左侧面设UX=0，下侧面设UY=0，模拟对称边界条件；上端面施加100 MPa的均匀拉应力。

这里有个小经验：如果直接用Pressure施加载荷，要特别注意压力方向与单元面法向的关系。为了避免方向搞错，我在上端面专门建了一层表面效应单元Surf154，再施加PRES=100 MPa，方向设置为向外拉伸。这样既保证了载荷均匀，也避免了单点加载带来的局部应力扰动。

全局模型求解后直接看最大应力，孔边最大应力算出来大约215 MPa，比理论值280 MPa低了约23%。这就是前面说的网格敏感性问题，粗网格把孔边高梯度应力平均掉了。但从整体云图分布看，高应力区位置正确，力流走向合理，作为子模型的“趋势底图”完全够用。

3.3 子模型：把关心区域切出来细算

子模型区域取孔中心周围80 mm×80 mm×5 mm的局部块，也就是四分之一模型里x从0到40 mm、y从0到40 mm的区域。这块区域的边界距离孔边约30 mm，大约是孔径的1.5倍，符合切割边界的安全经验值。

网格方面，孔边单元尺寸取0.8 mm，向外过渡到2 mm，厚度方向保持足够层数，采用SOLID186高阶六面体单元。这样划分下来，子模型单元规模大约三万多，求解时间几十秒。网格过渡要平滑，不要在切割边界附近突然变粗，否则边界插值位移的精确度会打折扣。

子模型上需要施加两类边界条件。第一类是模型自身的对称约束，左侧面UX=0，下侧面UY=0。第二类是切割边界位移，通过CBDOF命令从全局结果文件中读取。这里要特别强调一个新手容易犯的错：子模型是从整体模型内部切出来的，并没有穿过上端面载荷施加面，所以千万不要在子模型上再额外施加一遍100 MPa的拉力。全局载荷的效果已经通过切割边界位移传递进来了，重复施加载荷会导致结果严重偏大。

还有一个细节，如果切割边界节点恰好落在对称面上，ANSYS在加载时会跳过已经约束的自由度。为了保险，我习惯先把对称约束写好，再执行CBDOF。如果求解时弹出约束冲突相关警告，就去检查是否有节点被重复定义。

3.4 结果对比：子模型精度到底提升了多少

子模型求解完成，可以直接查看孔边最大应力。用0.8 mm网格算出来的结果大约是270 MPa，和理论值280 MPa相比，误差约3%。这个精度对于工程强度校核来说已经非常理想了。

注意看对比：子模型只用局部三万多单元，就达到了全模型加密到1 mm的精度水平，而求解时间缩短了一个数量级。如果把全局模型加密到0.5 mm、0.2 mm，计算成本还要翻几倍，这就是子模型在工程分析里真正的价值：用最小的代价拿到关键区域的精确结果。

做完子模型之后，我一直保留一个验证步骤：把子模型切割边界上的一圈节点提取出来，和全局模型在同一位置的应力值做对比。如果两者偏差在5%以内，说明切割边界取得足够远、边界条件传递没有明显失真；如果偏差很大，说明切割边界离应力集中区太近，需要把子模型范围再往外扩大。这个验证动作虽然简单，但能避免大量“假收敛”的判断失误。

4. 常见问题与排查技巧实录

4.1 结果对不上？按顺序查这四件事

子模型结果和预期差距大时，我先按固定顺序排查。第一是单位体系，全局模型和子模型必须完全一致，这是最基础也最容易忽略的问题。第二是切割边界位置，边界离高应力区太近，圣维南原理失效，结果必然不对。第三是局部网格质量，切割边界附近的网格不能太疏，更不能出现畸形单元，否则位移插值精度受影响。

第四是载荷重复施加。很多人子模型算出来应力偏大，一查就是自己在子模型上按照全局模型的习惯又施加了一遍外力载荷。记住一点：切割边界位移已经包含全局载荷的效应，子模型只需要补充它自身特有的约束和局部载荷。

我把几个高频问题整理成一张速查表，方便排查：

4.2 切割边界到底切在哪里才安全

切割边界位置的选择，是子模型分析中最需要经验判断的环节。太保守，子模型区域太大，细网格单元数量上去，计算优势就没了；太激进，边界贴近应力集中区，边界位移本身就不准确，局部结果被污染。

我常用的经验法则：切割边界到关注区域的距离，不要小于1.5到2倍特征尺寸。所谓特征尺寸，就是引起应力集中的那个几何量，孔径、圆角半径、齿高都可以算。比如孔边应力集中，边界离孔边至少1.5倍孔径；齿轮齿根分析，边界离齿根圆角底部的距离至少1.5到2倍齿高。这只是起步值，实际还是以边界应力一致性验证为准。

验证方法我刚才已经提过，就是把子模型切割边界上的应力结果和全局模型对应位置的应力画在同一张曲线图上。如果两条曲线贴合良好，说明边界位置足够安全；如果偏差明显，果断把子模型区域往外扩大一圈再算。这个验证环节多花十分钟，换来的是结果可信度的大幅提升。

4.3 单元类型与网格细节的取舍

子模型的精度高度依赖网格质量，特别是孔边、圆角这类区域。三维模型优先选择六面体高阶单元SOLID186，其次选择高阶四面体SOLID187，尽量避免低阶四面体SOLID285这类单元，它们太“硬”，弯曲和应力梯度捕捉能力差。曾经有一个同事图省事，子模型全部用低阶四面体，算出来的齿根应力比六面体结果低了差不多15%，白折腾了两天。

孔边或者圆角切点位置，网格要保证足够的层数，我习惯在最大应力可能出现的路径上至少布置5层以上单元。从关注区域向外围过渡时，单元尺寸逐步放大，不要出现相邻单元尺寸突变超过3倍的情况。切割边界附近的网格也不能太粗，毕竟这里的位移还要从全局结果插值过来，边界网格本身的分辨率会影响插值精度。

对于对称模型，还要注意对称面与切割边界重叠时，节点约束和边界位移可能出现重复。我处理时会把对称约束先写好，然后执行CBDOF，ANSYS一般会正确处理已约束自由度。如果求解过程出现约束冲突的警告，不要无视，回到模型里检查对应节点的约束状态。

4.4 哪些情况不适合用子模型

子模型不是万能药，有些场景硬套会得出完全错误的结果。首先是大范围材料非线性问题，如果结构发生大范围塑性流动，局部应力重分布明显，全局模型粗网格对塑性区和位移场的描述本身就不可靠，切割边界位移带过来的误差会淹没局部细节。这种情况应该用完整的细化模型，或者使用自适应网格方法。

其次是接触问题，尤其是有滑移的接触。接触状态的改变会显著影响局部刚度与力流路径，切割边界位移无法准确传递这种强非线性信息，子模型算出来的接触压力往往与真实情况偏差很大。这时候更适合的做法是对接触区域单独建立精细化接触模型，或者用子结构技术配合接触迭代。

还有一类情况是结构刚度主要由局部细部特征决定，比如薄壁结构的局部屈曲问题。这类问题的失稳行为对网格和边界条件极其敏感，切割边界位移已经隐含了“全局模型判断该区域不会失稳”的假设，子模型自然算不出屈曲。断裂力学中的裂纹尖端分析虽然可以使用子模型配合奇异单元，但需要额外的特殊处理，不建议新手直接套用常规子模型流程。

最后再说两句

做的子模型越多，我越觉得这个技术不是简单的“步骤背诵”，而是对计算力学基本概念的一次综合运用。圣维南原理、位移插值、网格收敛性、边界条件合理性，每一个环节都在考验对“近似”的理解：全局模型追求趋势正确，子模型追求局部精确，两者配合，才是工程分析的常态思维。

我自己的习惯是，无论子模型算得多顺，最后都留一步做验证，把切割边界一圈节点的应力从子模型和全局模型里分别提出来，叠加画成曲线。这个动作多花五分钟，却是我判断结果能不能交付的最重要依据。子模型不是万能的，但在线弹性局部精化这个场景下，它依然是我最常用的高效工具。希望这篇实例能帮你把这个工具用得顺手。

本文还有配套的精品资源，点击获取