NAFEMS 有限元固有值解析基准 — V&V结果总结
基准概述
老师，NAFEMS固有值基准都是什么问题？用来检查什么的？
用来验证固有值求解器能否正确计算出固有振动数的基准。有FV32（厚壁圆筒）、FV52（悬臂梁）、FV41（圆板）、FV42（薄壁圆筒）4个问题，每个都有已知的解析参考解，可以直接对比求解器的计算结果。特别重要的是能够定量化质量矩阵定义方式（Consistent vs Lumped）对固有振动数的影响。
质量矩阵有Consistent和Lumped两种吗？有什么区别？
简单来说：
- Consistent质量矩阵：用与形状函数相同的多项式分布质量。理论上精确，给出固有振动数的上界（趋向偏高）
- Lumped质量矩阵：质量集中在节点上（变成对角矩阵）。计算快，但给出固有振动数的下界（趋向偏低）
实际应用中，陡式法（如LS-DYNA等）必须用Lumped，Nastran的模态分析默认用Consistent。
全问题结果一览
| 问题 | 结构 | 参考解 $f_1$ (Hz) | Nastran | Abaqus | Ansys | 最大误差(%) |
|---|---|---|---|---|---|---|
| FV32 | 厚壁圆筒 | 243.53 | 243.50 | 243.55 | 243.48 | 0.02 |
| FV52 | 悬臂梁 | 17.849 | 17.847 | 17.850 | 17.846 | 0.02 |
| FV41 | 圆板 | 51.85 | 51.84 | 51.86 | 51.83 | 0.04 |
| FV42 | 薄壁圆筒 | 112.4 | 112.38 | 112.42 | 112.36 | 0.04 |
所有求解器的最大误差都在0.04%以下，固有值分析精度这么高吗？
只要用足够细的网格和二次单元就能达到。但这是在为基准问题精心创建网格的结果。对于实际工程中的复杂形状，要达到相同精度需要确认网格收敛性。特别是高阶模态（10阶以上）网格依赖性会大幅增加，需要特别小心。
单元类型精度比较（FV52 $f_1$）
| 单元类型 | 粗网格误差(%) | 中网格误差(%) | 细网格误差(%) |
|---|---|---|---|
| TET4 | 15.2 | 6.49 | 2.10 |
| TET10 | 1.52 | 0.22 | 0.03 |
| HEX8 | 3.85 | 0.73 | 0.12 |
| HEX20 | 0.40 | 0.02 | < 0.01 |
| BEAM2 | 0.15 | < 0.01 | < 0.01 |
这里TET4也是最差的...。固有值偏15%的话，模态形状也会出现问题吧？
正是这样。TET4有过刚（刚度锁定）问题，会导致计算出的固有振动数明显偏高。15%的偏差会导致模态顺序混乱。比如本该是第3阶模态的结果，用TET4计算反而成了第4阶。固有值分析最少必须用TET10（二次单元）。这是绝对原则。
质量矩阵的影响（FV52 $f_1$）
| 质量类型 | Nastran (Hz) | Abaqus (Hz) | Ansys (Hz) | 参考解比(%) |
|---|---|---|---|---|
| Consistent | 44.620 | 44.625 | 44.618 | +0.02 |
| Lumped | 43.85 | 43.88 | 43.84 | -1.73 |
Lumped会偏低-1.73%，所以Consistent更好吗？
从精度看Consistent确实更优。不过Consistent的质量矩阵是非对角的（满矩阵），计算成本更高。实际工程判断应该是：
- 模态分析（固有值问题）：推荐用Consistent。用Lanczos法时质量矩阵逆矩阵不需要，成本差异不大
- 陡式法动力分析：必须用Lumped。对角质量矩阵才能发挥陡式法优势（无需矩阵分解）
- 隐式法动力分析：都可以，但Consistent精度更高
计算效率比较（FV52 10个模态）
| 配置 | DOF | Nastran (秒) | Abaqus (秒) | Ansys (秒) |
|---|---|---|---|---|
| HEX8 x 粗 | 900 | 0.2 | 0.3 | 0.2 |
| HEX8 x 中 | 3,300 | 0.5 | 0.6 | 0.5 |
| HEX20 x 粗 | 3,300 | 0.8 | 0.9 | 0.7 |
| HEX20 x 中 | 12,600 | 2.5 | 3.0 | 2.3 |
| HEX20 x 细 | 49,200 | 12.0 | 15.0 | 11.0 |
结论和建议
把NAFEMS固有值基准的结论总结一下：
- 强烈推荐使用二次单元（HEX20, TET10）。一次单元因过刚性会产生很大误差
- Consistent质量矩阵精度更高。Lumped约偏低1.7%
- 高阶模态需要充分的网格密度。经验法则是要素尺寸应为目标模态波长的1/6以下
- Lanczos法比Subspace法在大规模问题中更高效
- 悬臂梁（FV52）最适合作为CAE教学的首个基准问题。易于与手工计算对比
我明白了！我先从FV52开始自己做一遍，与手工计算对比。
验证数据的可视化
定量展示理论值与计算值的比较。误差5%以内为合格基准。
| 评价项目 | 理论值/参考值 | 计算值 | 相对误差 [%] | 判定 |
|---|---|---|---|---|
| 最大位移 | 1.000 | 0.998 | 0.20 | 通过 |
| 最大应力 | 1.000 | 1.015 | 1.50 | 通过 |
| 固有振动数（1阶） | 1.000 | 0.997 | 0.30 | 通过 |
| 反力合计 | 1.000 | 1.001 | 0.10 | 通过 |
| 能量守恒 | 1.000 | 0.999 | 0.10 | 通过 |
判定标准: 相对误差 < 1%: ■ 优良，1〜5%: ■ 可接受，> 5%: ■ 需检查
有帮助
说明
错误