FEA有限元分析的误差来源诊断与网格无关性验证

误差不是越少越好，是要知道误差在哪

FEA有限元分析本质上是一种近似方法——连续体被离散为有限个单元，精确解被近似解替代。追求”零误差”既不可能也无必要，工程需要的是知道误差有多大、在哪个方向、是否可接受。

本项目中遇到过两种极端态度：一种认为”软件算出来的就是对的”，不做任何验证；另一种追求极致精度，网格加密到上百万单元，计算几小时却不知道结果是否更可靠。正确做法是建立系统化的误差诊断流程，有的放矢地控制关键误差源。关于有限元分析的完整方法体系，站内有系列文章。

四层误差分类

边界条件误差的量级最大且最容易被忽视——一个错误的固定约束可以使局部应力偏差100%以上。材料数据误差虽然不可完全消除，但可以通过灵敏度分析确定哪些参数影响最大。

网格无关性验证

标准流程

Richardson外推法

假设离散化误差具有渐近形式：

f(h) = f_exact + C·h^p

其中h是网格尺寸，p是收敛阶，C是常数。用三套网格的结果f₁, f₂, f₃（对应h₁>h₂>h₃）可以求解：

• 收敛阶：p = ln((f₂-f₁)/(f₃-f₂)) / ln(r)，其中r=h₁/h₂=h₂/h₃
• 精确解估计：f_exact ≈ f₃ + (f₃-f₂)/(r^p – 1)
• 细网格误差：e₃ ≈ |f₃-f_exact| / |f_exact|

实际案例

某支架结构应力分析，三套网格结果：

计算：

• r = 8/4 = 2（等比缩小）
• p = ln((168-142)/(176-168)) / ln(2) = ln(3.25)/ln(2) = 1.70
• f_exact ≈ 176 + (176-168)/(2^1.70 – 1) = 176 + 8/2.25 = 179.6 MPa
• 细网格误差 = |176-179.6|/179.6 = 2.0%

收敛阶p=1.70，接近理论值2（线性单元的理论收敛阶），说明网格质量正常。细网格误差2%，满足工程5%要求。

收敛曲线的异常形态

假收敛：当网格加密时结果看似收敛，但收敛方向错误。常见于过约束情况——约束越强，结果越”稳定”，但越不物理。区分方法：检查约束处的反力是否符合平衡条件。

能量范数误差估计

ANSYS和Abaqus内置了基于Zienkiewicz-Zhu的能量范数误差估计：

原理

对每个单元计算”恢复应力”（smoothed stress）和”计算应力”（FEA stress）之差，用差值的能量范数作为误差度量：

η_e = ||σ_recovered – σ_computed||_energy

全局误差比：η = √(Ση²_e) / √(Σ||σ_computed||²)

判定标准

局部误差图

ANSYS中的”Error”结果项显示每个单元的误差分布。本项目利用这个信息做自适应加密——在error>15%的单元上自动加密，迭代直到全局误差<5%。

模型简化误差

常见简化及其误差

Saint-Venant原理的验证：集中力影响范围约为载荷特征尺寸的2-3倍。在影响范围之外，集中力和分布力的结果差异<5%。因此只要关注区域距载荷足够远，集中力简化是合理的。

材料灵敏度分析

当材料参数有不确定性时，灵敏度分析可以确定哪些参数对结果影响最大：

本项目在某复合材料结构分析中做灵敏度分析：

E₁的灵敏度最高，说明这个参数的精度最需要保证。而E₂即使偏差10%对应力影响不到1%，可以适当放宽精度要求。

验证与确认（V&V）

区分验证和确认

验证方法

确认方法

误差报告模板

本项目在每份FEA分析报告中包含标准误差评估表：

关于科研学术网中FEA验证的更多案例，站内有详细报告。

FEA有限元分析的误差诊断不是可选步骤，而是分析交付的必要组成部分。一个没有误差评估的分析报告是不完整的——分析者不仅要给出数值结果，还要给出”这个结果有多可信”的量化评估。本项目坚持的准则是：每个分析结论都附带误差范围，让决策者知道安全裕度有多少是真实的、多少是仿真不确定性的。