FEA 网格基准：碟形弹簧分析的单元尺寸与密度准则

1. 网格基准的目的

有限元分析 (FEA) 的精度高度依赖于网格质量。网格基准 是指为保证计算结果（位移、应力、频率等）达到工程可接受精度而必须遵守的最小网格尺寸、单元层数和分布密度。对于碟形弹簧这类薄壁弯曲结构，网格过粗会导致“剪切锁死”、弯曲刚度过大、应力严重低估等问题。

网格基准通常来源于三个方面： - 理论要求：基于单元形函数插值能力，弯曲问题厚度方向至少需要多少层单元； - 经验法则：根据前人分析和试验对比总结出的推荐网格密度； - 收敛性验证：通过加密网格，比较结果变化，当变化小于目标容差时，确认网格足够。

2. 厚度方向网格数基准（弯曲主导问题）

碟形弹簧以弯曲变形为主，截面应力沿厚度线性分布（纯弹性）。为了准确捕捉上下表面的应力极值，单元在厚度方向的节点数必须能够描述线性变化的弯曲正应力。

• 线性完全积分单元（如 C3D8）：弯曲时会产生“剪切锁死”，需要足够多层数才能缓解。推荐 最少 4 层。
• 线性缩减积分单元（如 C3D8R）：弯曲性能较好，沙漏控制得当，不少于 3 层。
• 二次单元（如 C3D20、C3D15）：每个节点具有更高阶插值，厚度方向 至少 2 层 即可较准确描述弯曲应力。
• 二次四面体单元（C3D10）：在厚度方向也需要 不少于 2 层，且总单元数较多。

碟簧厚度方向推荐：对于厚度 $t$，采用线性六面体时厚度方向应划分 4~6 层，二次六面体时 2~3 层，以保证 OM 点和 I 点应力误差 < 3%。

3. 圆周方向网格数基准

碟形弹簧是轴对称体，但承受均布轴向载荷时变形仍为轴对称。圆周方向离散不足会使内外缘的几何描述粗糙，导致应力失真。

• 全模型：推荐圆周方向不少于 48 个单元（即每 7.5° 一个单元）。
• 对称模型（如取 1/4 或 1/8 扇区并施加对称边界）：扇区内圆周单元数应不少于 12~16，等效全周不少于 48 个单元。
• 如果仅关心力‑位移曲线：圆周方向单元数可稍少，如 36 个，但应力评估仍需 48 以上。

4. 基于应力梯度的网格尺寸估算

在应力梯度大的区域（如内外缘），单元尺寸应小于应力变化特征长度。工程上常用如下经验关系估算局部单元尺寸：

其中 $L_{grad} = \left| \frac{\sigma}{\partial \sigma/\partial x} \right|$ 为应力变化的特征长度。可先通过一个较粗的网格预分析，找出应力梯度最大的区域，然后局部细化，使细化区的单元尺寸满足上述关系。

对于碟形弹簧，内外缘倒角或圆角附近应力梯度极高，应进行局部网格细化，使圆角处至少有 4~6 个单元 覆盖半径。

5. 接触区域的网格要求

若分析中考虑碟簧之间的接触（并联/串联叠合），则接触面上的网格必须足够精细，以正确计算接触压力。

• 接触区域的单元尺寸应小于预期接触半宽的 1/3。
• 对于面-面接触，主面和从面的网格密度宜相近，避免刚度悬殊导致的穿透。
• 厚度方向仍应保证 ≥ 4 层，因为接触只发生在表面，但弯曲应力的准确传递要求足够的厚度层数。

6. 网格收敛性判断基准

确定网格是否满足精度要求，应进行至少三种不同密度的网格计算（如粗、中、细），比较关键结果量：

• 力值（反力）：要求误差 < 1%；
• OM 点最大压应力：要求误差 < 3%；
• I 点拉应力：要求误差 < 5%。

只有当连续两个密度之间的变化低于上述阈值时，才认为网格收敛。通常中等网格（厚度 4 层，圆周 48 单元）即可满足大多数要求。

7. 显式动力学网格基准

对于显式冲击分析，网格尺寸除了满足精度要求，还直接影响临界时间步长 $\Delta t_{crit}$。单元越小，$\Delta t_{crit}$ 越小，计算成本急剧上升。因此需要平衡精度与效率：

• 最小单元尺寸 $L_{min}$ 应不小于产生可靠接触和应力分布所需的尺寸，但也不宜过细，防止时间步长过小。
• 经验：碟簧厚度方向单元尺寸不应低于 $t/5$，否则时间步长可能降至 $10^{-8}$ 秒以下，导致计算无法在合理时间内完成。
• 可通过质量缩放适当增大时间步长，但须保证动能变化 < 2%，不显著影响动态响应。

8. 子模型网格基准

在使用子模型法时，全局粗网格可放松要求（但仍需能合理传递位移），而子模型内网格必须精细。子模型切割边界应远离应力集中区至少 2~3 倍厚度，网格在该区域内可逐步过渡。子模型内部，重点区域应满足前述弯曲网格要求。

总结：碟形弹簧 FEA 网格基准以厚度方向 4 层（线性单元）、圆周 48 单元为核心，结合应力梯度局部细化、接触网格匹配和收敛性验证，确保分析结果在工程可接受误差范围内。合理的网格基准是获得可靠载荷、应力和寿命预测的前提。