网格收敛·Richardson外推·GCI 返回
分析工具

网格收敛·Richardson外推·GCI

通过逐步细化网格观察解的变化,用Richardson外推估计真实解。实时确认收敛阶数和GCI。

网格设置

1.03.0
60300
0.2%5%

暂停时,拖动滑块即可即时更新结果。

网格细化与收敛实时可视化
单元数
单元尺寸 h
峰值应力
与收敛值的误差

自动生成报告

计算中...
理论·主要公式

$$f_{exact} \approx f_1 + \frac{f_1 - f_2}{r^p - 1}$$

Richardson外推(两网格):\(r=h_2/h_1\) 网格比、\(p\) 精度阶数、\(f_1\) 细网格解

$$\text{GCI} = F_s \frac{|f_2 - f_1|}{r^p - 1} \cdot \frac{1}{f_1}$$

网格收敛指标(GCI):\(F_s = 1.25\) 安全系数、GCI < 1% 判断为收敛

网格收敛·Richardson外推·GCI详解

🙋
什么是网格收敛?我听说CAE中网格越细计算结果越会改变,到底要细到什么程度呢?
🎓
网格收敛就是指:网格再细化计算结果基本不再改变的状态。比如汽车车盖变形计算,粗网格得到10mm,细网格得到12mm,再细化到12.1mm……这样的情况下,解就收敛到12mm附近了。这个模拟器中"维数"参数改为1D、2D、3D时,你会发现需要的网格数量增长方式完全不同。
🙋
原来如此!但是每次都用超细的网格计算太耗时了。能不能根据粗一点的结果来预测答案呢?
🎓
这正是Richardson外推法的作用!用3个不同粗度的网格计算,根据这3个结果推测"理论上无限细网格的答案"。在这个工具里输入粗、中、细三种网格下得到的关键数值(比如最大应力),就能计算出推估的"真实解"。实际工作中,这个推估值有时能让你省略掉最细网格的计算。
🙋
明白了!那这个推估有多可信呢?有没有数字标准来判断"已经收敛"?
🎓
有的,那就是GCI(网格收敛指数)。它表示"推估值和现在得到的解之间,如果再细化网格会改善百分之几"。工具里会根据GCI值用颜色表示:2%以下绿色、5%以下橙色、5%以上红色。实务中常见的做法是让GCI_fine降到3%以下来证明网格收敛。试着调整参数看GCI颜色的变化吧。

物理模型与主要公式

网格尺寸$h$与数值解$\phi$的关系,由理论精度阶数$p$按如下方式建模。这是Richardson外推的基础:

$$\phi = \phi_{exact}+ g h^p + \cdots$$

其中$\phi_{exact}$是真实解(连续体解),$g$是常数,$h$是代表网格尺寸。$p$是离散化格式的理论精度阶数(如中心差分为2阶),但实际问题复杂,所以用3个解计算的"表观精度阶数"代替。

从三个不同粗度的网格(细:1,中:2,粗:3)计算表观精度阶数$p$、外推解$\phi_{ext}^{21}$和GCI。

$$ p = \frac{1}{\ln(r_{21})}\left| \ln \left| \frac{\epsilon_{32}}{\epsilon_{21}}\right| + q(p) \right| $$ $$ \phi_{ext}^{21}= \phi_1 + \frac{\phi_1 - \phi_2}{r_{21}^p - 1}$$ $$ GCI_{fine}^{21}= F_s \frac{|\epsilon_{21}|}{r_{21}^p - 1}$$

$r_{21}=h_2/h_1$为细化比,$\epsilon_{21}=\phi_2-\phi_1$为解差,$F_s$为安全系数(通常1.25)。$q(p)$是当$p$未知时的修正项。工具基于这些公式,从输入的3个解自动计算$p$、外推解和GCI,评估收敛性。

常见问题

基本做法是用3段以上系统细化的网格,保持细化比$r$常数(如r=2)。例如粗网格单元数N、中为2N、细为4N。比值不均匀会导致表观精度阶数$p$计算不稳定,外推结果失真。
不一定。外推解只有在解处于渐近区(网格足够细的区域)时才可信。要确认GCI较小、表观精度阶数$p$接近理论值。渐近区外外推解可能偏离真值很远。
一般来说GCI < 5%即解不依赖网格。但因应用不同而异。关键是粗、中、细各网格的GCI应呈现逐级递减趋势。若GCI在某个阶段突然反弹,说明网格尚未进入渐近区。
可能原因:(1)网格还太粗,未进渐近区;(2)网格质量不佳(歪斜、宽高比差);(3)边界条件或物理模型有间断点;(4)数值格式未正确实现。首先应尝试进一步细化网格观察。

实际应用

结构强度分析(汽车、飞机):评估部件最大应力和变形时,必须消除网格依赖性,得到可信结果。轻量化设计中安全系数较小,GCI定量验证往往是设计决策的基础。

流体分析(CFD):翼型升力、阻力系数和热交换器换热、压降等预测,广泛用此法提升精度。复杂湍流模型时理论阶数$p$未知,靠计算的表观阶数判断收敛。

电磁场分析:天线辐射方向图、共振频率、器件寄生参数等网格敏感量,Richardson外推在保持计算成本低的同时实现高精度估计。

学术研究·论文发表:数值计算结果可信性的重要验证手段,已成为可重复研究的标配。建议文章中同时报告网格尺寸和GCI值。

常见误区与注意事项

第一点,"网格越细越接近真值"不总是对。几何奇异点(尖角、接触点)附近,精细网格也可能应力发散。此时应从"求最大应力值"改为"研究应力分布",重新定义CAE目标。

第二点,3个网格必须"系统地"细化。别随意选1mm、0.7mm、0.5mm,而要保持细化比$r$恒定(推荐$r \ge 1.3$),如1mm、0.77mm(1/1.3)、0.59mm(1/1.3²)。否则表观阶数$p$计算失稳,外推无信度。

第三点,不要被GCI迷蛊。GCI_fine < 2%不代表所有物理量都收敛。常见情况是位移收敛但应力未收敛。必须对关键输出量逐一验证收敛性。

使用指南

  1. 输入3段网格(粗、中、细)的单元数N和对应计算结果φ
  2. 从粗网格和中网格间的相对变化计算网格比r,迭代求解收敛阶数p
  3. 用Richardson外推法推估真实解φ₀,用GCI值量化误差范围并判断收敛性(GCI < 5%为充分收敛)

计算示例

2D翼型绕流分析Cp分布:粗网格15000单元(φ₁=-0.856)、中网格42000单元(φ₂=-0.924)。网格比r=(42000/15000)^(1/2)=1.673,计算得收敛阶数p≈1.8,Richardson外推得φ₀≈-0.971为真值估计。GCI=(|φ₂-φ₁|/φ₁)×100/(r^p-1)=2.1%在许可范围内,可判定无需进一步细化到超细网格(120000单元)。

实务中的注意