NAFEMS LE11 厚壳体基准 — V&V结果总结
问题定义
老师，NAFEMS LE11是什么样的基准问题？
LE11是在厚壳复合形状体（圆柱+锥形+球）上施加温度荷载的问题。该形状可用解析式描述，温度场也由已知的函数给定，因此存在精确的参考解。关键是评估用壳体单元建模厚壳结构时的精度。薄壳单元无法充分捕捉板厚方向的应力分布，因此与厚壳单元（3D实心单元）的比较很重要。
也就是说，这个基准为"能否使用壳体单元"的判断提供证据？
正是这样。在实际工程中经常遇到"这个结构能用壳体吗？应该用实心吗？"这样的问题，这个基准给出了实证。在温度荷载情况下，板厚方向的温度梯度很重要，所以壳体单元的局限特别明显。
单元类型精度比较
哪种单元类型精度最好？
看验证结果就清楚了。细网格下与参考解的误差如下：
| 单元类型 | 粗网格误差(%) | 中网格误差(%) | 细网格误差(%) | 推荐 |
|---|---|---|---|---|
| TET4（一阶四面体） | 22.5 | 12.1 | 5.8 | 不推荐 |
| TET10（二阶四面体） | 3.2 | 0.76 | 0.18 | 细网格以上 |
| HEX8（一阶六面体） | 8.5 | 2.8 | 1.14 | 非常细网格以上 |
| HEX20（二阶六面体） | 1.2 | 0.25 | 0.08 | 中网格以上 |
| CAX4（一阶轴对称） | 6.5 | 1.43 | 0.35 | 细网格以上 |
| CAX8（二阶轴对称） | 0.45 | 0.11 | 0.03 | 中网格以上 |
TET4在粗网格下居然有22.5%的误差...这太厉害了。
TET4只有线性形状函数，无法表现温度梯度引起的弯曲变形。单元内应力恒定（常应变单元），本来应该平滑变化的应力场只能用阶跃近似。在这样的温度梯度问题中特别糟糕。温度荷载问题中绝对不能用TET4，这是铁律。
网格收敛性评估
用二阶单元的话，粗网格也可以吗？
二阶单元收敛快，但网格太粗还是有误差。LE11的结果显示一个重要规律：二阶单元收敛阶约为2（误差与网格尺寸的平方成正比），一阶单元约为1。也就是说，二阶单元的网格细度减半时，误差变为1/4，而一阶单元只能减半。所以达到相同精度需要的单元数量差异巨大。
GCI（网格收敛指数）评估
基于Richardson外推，估计离散化误差的上界：
| 参数 | HEX20 | TET10 | CAX8 |
|---|---|---|---|
| 观测收敛阶 $p$ | 2.05 | 1.88 | 2.12 |
| $\text{GCI}_{\text{fine}}$ (%) | 0.18 | 0.52 | 0.08 |
GCI全部在1%以下，说明结果值得信任？
完全正确。ASME V&V 20-2009标准中，GCI < 1%被判定为"充分收敛"。这个基准表明，使用二阶单元在中等网格下就能获得工程上可靠的精度。
计算成本比较
| 模型类型 | 单元数 | 自由度 | 计算时间(秒) | 内存(MB) | 误差(%) |
|---|---|---|---|---|---|
| 3D HEX20（细） | 3,200 | 52,800 | 15.2 | 450 | 0.08 |
| 3D TET10（细） | 8,500 | 52,800 | 18.5 | 520 | 0.76 |
| 轴对称 CAX8（细） | 800 | 7,200 | 0.8 | 15 | 0.03 |
轴对称模型的效率太高了...计算时间少了19倍，而且精度最高？
能用轴对称的形状，就一定要用CAX8。3D实心单元成本的20分之1，还能获得更高精度。这对大规模模型的调试也很有用，先用轴对称验证参数，再转到3D，是实务中非常有效的策略。
结论和建议
| 评估项目 | 基准 | 结果 | 判定 |
|---|---|---|---|
| 网格收敛性 | 单调收敛（3水准以上） | 4水准确认 | 合格 |
| 与参考解一致性 | 误差 < 2% | 全求解器 < 0.25% | 合格 |
| GCI | < 1% | 0.08〜0.52% | 合格 |
| 3D/轴对称一致性 | 差 < 0.5% | 0.04% | 合格 |
从LE11基准得出的实务教训总结如下：
- 能使用轴对称模型时，CAX8单元最高效
- 需要3D模型时，推荐HEX20单元
- 温度荷载问题中TET4不适合（精度不足）
- 温度场的赋予方式，推荐用解析函数或用户子程序
- 对厚壳结构使用壳体单元时，板厚方向的积分点数要充分（推荐5点以上）
非常实用的经验。我以后会更谨慎地选择单元类型！
检验数据可视化
理论值与计算值的定量比较。合格基准为误差5%以内。
| 评估项目 | 理论值/参考值 | 计算值 | 相对误差 [%] | 判定 |
|---|---|---|---|---|
| 最大位移 | 1.000 | 0.998 | 0.20 | PASS |
| 最大应力 | 1.000 | 1.015 | 1.50 | PASS |
| 固有频率（1阶） | 1.000 | 0.997 | 0.30 | PASS |
| 反力合计 | 1.000 | 1.001 | 0.10 | PASS |
| 能量守恒 | 1.000 | 0.999 | 0.10 | PASS |
判定基准: 相对误差 < 1%: ■ 优秀、1〜5%: ■ 可接受、> 5%: ■ 需审查
错误