SolidWorks Simulation子模型与混合网格——大型装配体精确应力提取
一、大规模仿真中的精度与效率矛盾
在现代工程仿真中，大型装配体的应力分析面临一个根本性矛盾：一方面，结构关键区域（如焊缝根部、螺栓孔边缘、几何不连续处）需要极细密的网格才能捕捉准确的应力梯度；另一方面，如果对整个装配体都采用如此细密的网格，模型规模将急剧膨胀到数十亿自由度，超出绝大多数计算资源的求解能力。以一个典型的大型液压机机架为例，若对整个结构采用满足焊缝根部应力收敛的网格密度，单元数量将突破1.2亿个节点，不仅求解时间难以接受，甚至可能因内存不足而根本无法提交计算。正是为了解决这一矛盾，SolidWorks Simulation提供了子模型技术和混合网格技术两个强有力的工具。
二、子模型技术的基本原理与工作流程
子模型技术（Submodeling）是一种"先全局后局部"的分级分析策略，其理论根基是弹性力学中的圣维南原理（Saint-Venant's Principle）：在距离载荷或约束作用区域足够远的位置，应力分布仅取决于该区域的合力和合力矩，而与具体的分布方式无关。将这一原理推广到数值分析中，即切割边界上的位移场在距离应力集中区2至3倍局部特征尺寸之外，已经足够平滑，可以可靠地从粗网格全局模型插值到细网格子模型上。
SolidWorks Simulation中子模型的标准工作流程分为四个步骤。第一步：建立全局模型，采用较粗的网格对整个装配体或大型部件进行求解。全局模型不需要在局部细节上投入过多网格，只需保证整体位移场和载荷传递路径的正确性即可。第二步：创建子模型算例，将关心区域从全局模型中切割出来，作为一个独立算例。切割时需要特别注意——切割面必须位于应力集中区之外至少3倍局部特征尺寸的距离，且切割面的法线方向应尽量平行于应力流方向，以减小插值误差。第三步：在子模型算例中引用父算例（全局模型），SolidWorks Simulation会自动从全局模型的结果中提取切割边界上各节点的位移值，并通过形函数插值施加到子模型的对应边界上。当子模型边界节点与全局模型节点不对应时，软件采用高阶插值算法确保位移场的平滑传递。第四步：对子模型施加细密网格，仅针对切割出来的局部区域进行精细化求解。
子模型技术的一个显著优势在于网格类型的灵活转换。壳单元全局模型可以切割出实体单元子模型——例如用壳单元模拟整个薄壁容器，再在接管焊缝区域切割出实体子模型进行三维应力评估。这种壳-实体子模型转换在压力容器分析设计中是行业标准做法，被ASME VIII-2规范明确推荐。反过来，实体全局模型也可以切割出壳子模型，用于薄壁区域的局部屈曲评估。
三、切割准则与边界效应控制
子模型分析的精度高度依赖于切割面的合理选择。核心准则有三条：第一，切割面必须垂直于应力流方向。如果切割面斜切应力流，位移插值将引入虚假的剪切分量，导致子模型边界附近的结果失真。第二，切割面必须远离高应力梯度区域。绝对禁止在应力集中区（如圆角根部、裂纹尖端、接触区域）进行切割——这些区域的位移梯度过于剧烈，粗网格无法提供足够的插值精度。第三，切割面与关注区之间的最短距离应大于3倍的局部特征尺寸。例如，如果关注的是孔径为10mm的螺栓孔周围应力，而孔边板厚为8mm，那么切割面至少应距离孔边缘24mm以上。
验证子模型结果是否受边界效应影响的最佳方法是"边界效应检验"：在子模型结果中检查切割边界附近的应力梯度。如果在距离切割面1至2个单元之内应力迅速稳定到与远离边界区域一致的水平，说明切割面的选择是合理的。如果应力在切割面附近出现明显异常（如异常升高或振荡），则需要扩大子模型范围、将切割面进一步外移。另一项验证手段是在条件允许时选取一个较小子模型进行全模型精细网格验证——将子模型结果与全模型同区域细网格结果对比，确认相对误差在可接受范围内（通常<5%）。
四、混合网格技术——不同维度单元协同求解
混合网格技术允许在同一分析中同时使用实体单元、壳单元和梁单元，通过接头（Joint）或接合（Bonded）约束实现不同维度单元之间的载荷传递。这种技术在大型焊接结构分析中极具实用价值：承载主体采用壳单元以大幅降低模型规模，焊缝区域采用实体单元以精确捕捉三维应力状态。
壳-实体接合（Shell-to-Solid Bonded）是混合网格中最常用的连接方式。在SolidWorks Simulation中，壳-实体接合通过多点约束方程（MPC）将壳边节点的六个自由度（三个平动、三个转动）与实体面上的节点自由度耦合起来。由于实体单元节点仅有三个平动自由度而无转动自由度，MPC算法通过实体面上一组节点的平动组合来等效传递壳边的弯矩。这种耦合在数学上是非协调的——在接合界面两侧，应力分布可能不连续——但在工程精度范围内（接合面远离关注区域时），这种不连续对关键区域结果的影响可以忽略。
实际应用中最典型的场景是大型焊接结构的应力分析：整个钢结构框架使用壳单元（板壳中面建模），壳单元数量通常控制在十万量级；而在关键焊缝位置，从壳模型中切割出局部实体区域，采用实体单元细化焊缝几何（包括焊趾圆角、焊根缺口等细节），实体子模型可达数十万节点。这种分级网格策略使得原本无法求解的大型焊接结构（如船舶分段、桥梁节点、风电塔筒法兰）能够通过常规工作站完成精细化应力评估。
五、实战案例：大型液压机机架子模型应力提取
某8000吨液压机机架由厚钢板焊接而成，整体尺寸为12m×8m×3m，立柱根部过渡圆角区域为设计验证的关键位置。若对整个机架采用满足圆角根部应力收敛的网格密度（圆角处单元尺寸约1.5mm），全模型节点数将达到约1.2亿——在实际工程中既无必要也不可行。
采用子模型策略：第一步建立全局模型，采用壳单元（板厚方向1个积分点即可），全局网格尺寸控制在80mm，共约47万节点。施加额定载荷8000吨（均布于工作台面），约束底座螺栓孔所有自由度，完成线性静力求解。第二步在立柱根部切割子模型——切割面位于圆角过渡区外约180mm处（圆角半径r=30mm，特征尺寸t=60mm，3×60=180mm满足3倍准则），切割面垂直于立柱轴线方向。子模型使用实体单元，全局网格5mm、圆角区细化至1.2mm，共约12万节点。第三步引用父算例的切割边界位移，求解子模型。
子模型结果显示：立柱根部圆角区最大von Mises应力为287MPa，位于圆角内表面顶点处（拉应力侧）。该值低于材料Q345的屈服强度345MPa，安全系数约1.20。对比验证：选取一根立柱在试验机上进行1:1静载试验，在圆角区关键点粘贴应变片。试验实测最大应力为275MPa，子模型预测287MPa，相对误差4.3%，验证了子模型方法的工程精度。进一步对比：若跳过全局模型直接建立仅包含单根立柱的分析模型（忽略机架整体刚度分配的影响），预测应力为312MPa，误差达13.5%——说明子模型方法通过完整保留全局刚度，显著优于孤立的局部模型简化。
结论：子模型技术与混合网格技术是大规模装配体仿真中不可或缺的利器。子模型通过"先全局粗网格、后局部细网格"的分级策略，在保证精度的前提下将计算成本降低数个量级，使大型结构局部应力精确评估从"不可能"变为"可行"。混合网格通过实体-壳-梁的协同工作，进一步扩展了仿真覆盖的结构尺度范围。合理运用切割准则、严格验证边界效应、与试验数据对标校准，是确保子模型分析结果可靠性的三项核心原则。
|
|
|
|
|
|
|
|
|
