Altair（HyperWorks）有限元仿真教程【网格密度和结果收敛性】

网格密度和结果收敛性密切相关，但影响两者相互关系的因素存在争议。为了满足规定的时间和精度要求，用户需要在建模时间、模型准确度、计算时间以及人力成本之间做出权衡。

从数值精确性的角度来说，准确的网格是指网格细化以后，仿真结果没有出现显著的差别。尽管这个概念听起来比较简单，但很多因素是需要考虑的。网格细化必须精确的描述待分析的对象。

通过单元分割细化网格的方法可能会产生一些问题，除非新建的节点足够接近原始几何。网格细化时，要注意原始单元的使用假设。比如，当壳单元细化到不再有效时，就要用实体单元代替。

确定网格密度时，需要遵循以下原则。

几何细节要求

确定获得预期结果需要描述的最小几何特征。

一个尖锐特征可能会导致应力集中，但该特征可能不在零件的承重区域。要抓住这个细节则需要在整体结构分析后进行一次局部结构分析。

设计细节可用性

检查设计数据的可用性。如果设计数据不完整或使用有限元分析的方法帮助设计，最好使用简化模型，但要注意不能过度简化模型以致忽略了某些研究项目。

设计完成后，进行最终验证时，尽可能考虑详尽的细节并在时间允许条件下细化网格（包括结果精度检查时的网格细化迭代）。

与以前工作对比

如果要将当前仿真结果与相同或类似零件的分析结果进行对比，可以考虑使用类似网格密度的模型。如果考察不同模型的相关性，需注意检查不合理的边界条件、载荷工况、低质量的网格以及不恰当的网格密度等。如果发现测试缺陷，则需要重新建立测试标准。只有在完全理解别人工作的情况下，才可考虑使用他们的模型。

预期变形

确定可能出现的变形以及捕捉这些变形所需的网格密度。预估最大变形和弯曲位置，观察节点密度与变形模式之间的相互关系。注意以下变形模式：

收敛性检查

如果计划进行收敛性检查，至少需要进行一次网格细化。仿真结果中，如果相邻单元应力差别较大，则该区域可能没有很好的捕捉到应力变化。此时，建议进行一些网格细化工作。

一些程序可根据内置标准或用户自定义标准自动计算收敛性。另外，还有一些程序可以自动细化模型中非收敛区域的网格。这些功能均与载荷工况相关，不同的载荷工况将产生不同的网格模型。

P单元细化办法是通过增加每个单元或非收敛区域单元的复杂性实现的。所有自动化的功能均需要“近似完美”的模型，否则收敛性检查会将重点集中到与分析无关的小细节或满足一般收敛条件的单元、边界条件及载荷上。建模时，即使不打算进行收敛性检查，也尽量创建可以收敛的模型。例如，下图的平面结构使用壳单元模拟可获得较好的收敛性，而用“较为精确的”实体单元模拟则结果无法收敛。

因为壳单元在载荷方向不会产生应力，可获得与壳结构理论一致的收敛结果。相反，实体单元尝试计算“点载荷”，使得应力随着单元尺寸的减小而不断增大，以致结果不收敛。

变形或应力

多数情况下，如果只关心结构的变形和刚度结果，使用少量单元模拟结构即可。如果只关心结构变形情况，那么可以考虑使用更少量的单元。例如，分析矩形截面悬臂梁的变形时，一个梁单元或一个壳单元就可以了；而分析简支梁边界处的应力变化时，则需要大量的壳单元或实体单元来模拟。如果一个区域中出现了较大的应力变化，如框架一角存在网状结构，为获得合理的变形结果，需要使用大量的单元来模拟该区域。

关注领域

确定关注领域。如果能够提前预测出现问题的位置，则需要仔细对待这些位置。如果时间紧迫，可以细 化该区域网格，其他地方使用粗糙网格近似模拟。如果没有相关经验并且没有一点线索，那么可以先使用均匀网格，然后再细化应力集中区域网格。根据工程经验，典型的应力集中区域出现在以下位置：

• 小孔

• 载荷作用点

• 边界定义点

• 连接点

• 狭窄特征

• 突变特征

单元类型

网格密度与使用的单元类型有很大关系，因此，针对不同的分析选择合适的单元类型是非常重要的。选择单元类型时，考虑该结构是否可以使用二维单元模拟，或者需要使用完全三维实体单元才能准确模拟结构的应力状态、几何以及局部变形情况。

例如，进行刚度和动力学分析时，发动机用壳和梁单元描述就可以了。但如果应力也作为一个考察项目，该部件就变成了三维实体，上述建模方式就无法获得合适的应力结果。发动机建模时，实体单元模型和壳、梁混合模型所需的时间差别很大。选择单元类型时应注意一下几点：

• 尽量不要混合使用三角形和四边形单元，如果必须使用三角形单元，应将其放置在非重点区域。

• 验证三角形单元是否适用某个分析时，应使用全三角形单元模型对比，不要对比四边形单元模型与细化

的三角形单元模型。

• 尽可能使用一致尺寸的单元。

• 尚未完全理解分析对象前，不要使用壳和实体混合单元或梁和壳/实体混合单元。

• 尽量用实体单元模拟实体结构

• 用壳单元模拟薄板结构 (厚度 < 10~20 倍边长)

• 了解所有选用单元的应用假设

• 不要使用退化单元

• 不要将四面体单元与其他类型单元混合使用。如果必须使用，应将其放在结构外表面位置。

• 如果创建六面体单元需要太多时间和精力，可以考虑使用四面体单元（有时候这样可以节省几周的时间）。

• 如果创建壳单元或者实体单元需要太多时间和精力，可以考虑使用梁单元。

使用一个壳单元模拟加固件是不合理的刚性梁网状加固件梁单元广泛用于模拟梁结构、紧固件及连接件等。然而，梁单元的复杂性可能成为各类单元中最易用错的单元。正确使用梁单元，需要注意以下几点：

• 剪切中心

• 翘曲约束

• 长深比

• 剪切变形

• 边界处的复杂应力状态

• 分析结果可视化

复杂截面梁的应力结果只有通过经典梁理论验证后才能使用。这里，圆形截面梁是个例外。因为它容易建模，而且分析结果与理论结果匹配良好（局部加载位置除外）。

经过几十年的发展，我们已经可以利用梁单元模拟各种复杂几何结构（通常为了节省计算时间），但是很少人能够理解其中的复杂性。因此，选用这类单元时，应做好迎接复杂状况挑战的准备。

时间约束

确定网格密度时，时间约束是一个重要的考察因素。一个折中的方案是平衡人力耗时和计算机耗时。通常，使用网格自动划分工具创建大模型可以节省大量的时间。另外，不断增加的计算机性能也在逐步降低使用特殊功能缩减单元规模的需求。

因而，人力资源逐步集中到提高生产力和市场竞争力方面。选择建模功能和单元类型的同时还需要考虑外在因素，以获得最佳的投入产出比。一旦确定了网格划分方案，在实施之前，建议与他人讨论该方案的可行性。

确定模型复杂性时，另外一个需要考虑的因素是模型的“保质期”。如果模型可以很长时间使用，为了减少后续工作量，应该寻求一种简化的模型修改和定义方式。这是一个常常被忽略的问题，正确认识这一问题，可以节省很多时间。

网格密度和结果收敛性总结

1. 确定获得预期结果所需考虑的最小几何细节

2. 评估设计细节可用性

3. 检查当前项目的分析要求并与先前工作做比较

4. 预期结构变形以及该结果对网格的要求

5. 计划进行收敛性检查

6. 如果关注应力结果，应选择一个较为详细的建模方案。考虑使用一个两阶段的解决方案，详细定义第二个模型，以降低解决方案的广度和复杂度。

7. 确定关注领域

8. 选择合适的单元类型。使用混合网格时需要特别注意。

9. 考虑时间约束