2. 子模型分析流程与关键决策
子模型方法的核心思想是“两级求解、位移传递”:先在较粗的全局模型上求解整体响应,再把全局解在局部区域边界上的位移插值给一个独立的、网格更细的局部模型,作为其边界条件重新求解。理解这条主线,比记住步骤顺序更重要——它决定了每一步为什么必须那样做,以及哪些地方最容易出错。
一、完整流程:从全局解到局部解
一个可复现的子模型流程通常包含以下环节,其中前四步属于全局模型阶段,后四步属于子模型阶段。
- 全局模型建模与求解:建立覆盖整体结构的模型,网格可以较粗,但几何、材料、载荷与约束必须能反映真实受力路径。求解后必须确认全局解本身是可信的,否则子模型只是在一个错误的基础上做精细化。
- 确定局部关注区域:明确要精细化的对象,例如某个圆角、焊缝、螺栓孔、接触区或应力集中部位。区域范围应包含足够大的“过渡带”,使边界远离应力梯度剧烈变化的位置。
- 切割边界选取:在全局模型中划定子模型边界。边界应尽量落在应力变化平缓、位移梯度小的地方,通常取在距关注点2~3倍特征尺寸之外。
- 提取边界位移:从全局解中输出切割边界上节点的位移(对实体单元通常为三个平动自由度)。这一步是两级模型之间唯一的物理桥梁。
- 建立局部子模型:按关注区域的真实几何重新建模,可保留圆角、倒角、细节孔等全局模型中被简化的特征,并赋予更精细的网格。
- 施加驱动边界条件:把第4步提取的位移插值到子模型对应边界节点上,作为强制位移边界条件。子模型内部不再施加原有的外部载荷,除非该载荷确实作用在局部区域内。
- 子模型求解:提交计算,得到局部区域的高精度应力、应变或损伤结果。
- 结果验证与收敛判断:检查子模型边界处的应力是否与全局解一致,并做边界距离敏感性分析,确认结果不再随边界外移而明显变化。
在 Abaqus 中,第4步和第6步通常借助 *Submodel 与 *Boundary 配合驱动变量实现,也可通过 *Equation 或场变量插值方式传递。实际操作中,全局与子模型的单位制、坐标系必须严格一致,否则位移传递会整体错位。
二、驱动变量选择:位移是默认,但不是唯一
子模型的边界条件本质上是“把全局解在边界上的信息传给局部模型”。最常用、最稳健的驱动变量是位移,因为位移在有限元中是节点自由度,插值直接、物理意义明确,且对网格密度差异不敏感。对于实体单元,边界节点通常约束三个平动自由度;对于壳或梁,还需考虑转动自由度。
当分析涉及温度场、热应力或相变时,温度也应作为驱动变量传递。此时子模型边界上同时施加位移和温度,才能复现热-力耦合响应。若关注的是接触压力或螺栓预紧,则可能需要传递接触力或预紧力,但这类“力驱动”方式对边界位置更敏感,一般只在位移驱动不可行时使用。
选择驱动变量时有一个实用判据:该变量在边界上的梯度是否平缓。位移梯度平缓,插值误差小;若边界恰好穿过应力集中区,位移梯度大,插值误差会被放大,此时应外移边界或改用更精细的全局网格。
三、切割边界确定:位置比形状更重要
切割边界的位置直接决定子模型结果是否可信。基本原则是:边界必须远离应力集中区。圣维南原理告诉我们,局部细节对应力场的影响随距离衰减,通常取距关注点2~3倍特征尺寸之外,边界上的应力分布已接近整体应力场,插值误差可忽略。
边界形状可以灵活,矩形、圆形或沿几何特征线均可,但应尽量与全局网格线对齐,减少插值时的节点匹配误差。若全局与子模型网格不匹配,Abaqus 会自动进行插值,但边界节点过少会导致位移场失真。建议边界上至少保留两到三层全局单元,使位移插值有足够采样点。
四、精度与成本权衡:把计算力花在刀刃上
子模型方法的经济性来自“全局粗、局部细”的分工。全局模型网格粗,计算量小,用于捕捉整体变形和载荷路径;子模型网格细,计算量大,但只覆盖小区域,总成本远低于整体精细网格。以一个典型支架为例:全局模型约5万单元,子模型关注区约20万单元,若整体都用20万单元量级的网格,规模将放大数倍甚至一个数量级。
权衡时需关注三个参数:局部网格尺寸、边界距离、子模型区域大小。局部网格尺寸由关注特征的曲率决定,圆角处通常需要至少6~8个单元跨越圆弧;边界距离按前述2~3倍特征尺寸选取;子模型区域不宜过大,否则失去成本优势,也不宜过小,否则边界落入应力梯度区。
| 决策项 | 推荐做法 | 风险 |
|---|---|---|
| 驱动变量 | 优先位移,热分析叠加温度 | 力驱动对边界敏感 |
| 边界位置 | 距关注点2~3倍特征尺寸 | 过近导致插值误差放大 |
| 边界网格 | 保留2~3层全局单元 | 节点过少使位移场失真 |
| 局部网格 | 圆角处6~8单元跨弧 | 过粗无法捕捉应力集中 |
| 验证 | 边界外移20%~30%复算 | 不复算可能掩盖边界效应 |
五、常见误区与检查要点
第一,全局解本身不可信。若全局模型遗漏了重要载荷路径或约束不合理,子模型再精细也只是在错误基础上放大误差。第二,边界穿过应力集中区。这是最常见的失败原因,表现为子模型边界处应力与全局解明显不符。第三,单位制或坐标系不一致。全局与子模型若使用不同单位或局部坐标系,位移传递会整体错位。第四,忘记移除子模型内部的原载荷。子模型边界已由位移驱动,若再叠加外部载荷,会造成重复加载。
完成子模型求解后,建议做三项检查:边界处应力与全局解对比、边界外移敏感性、局部网格收敛性。三项都通过,才能认为局部精细化结果可信。子模型方法的价值不在于“算得更细”,而在于“在正确的位置算得更细”,这正是流程与决策要点需要被反复强调的原因。