验证 — CAE术语解释

这在ASME V&V 20标准中有明确定义。验证(Verification)是"是否正确求解方程"，确认(Validation)是"是否求解的是正确的方程"。例如，在弹性分析中，输入杨氏模量为200GPa时，求解器能否按照胡克定律

主要有两种方法。一种是代码验证，由软件制造商（如Ansys、Siemens等）进行，用来查找程序中的错误。另一种是计算验证，由用户进行，检查自己的网格、求解设置等是否合理。其中最具代表性的是网格依赖性调查。例如，在梁的挠度分析中，将单元尺寸从10mm、5mm、2.5mm逐步细化，检查结果是否收敛到某个值。如果不收敛，可能需要更细的网格，或重新调整特异点附近的网格。

完全同意，必须用定量指标。常用的是相对误差。对于感兴趣的量（如最大应力），比较网格i和i+1的结果，计算

很好的观点。主要误差包括(1)离散化误差、(2)舍入误差、(3)迭代误差。网格细化会减少(1)，但增加元素数时(2)会累积，(3)的迭代次数可能增加。例如，在Ansys Mechanical进行非线性分析时，将力的收敛标准从默认的0.5%改为0.1%会减少迭代误差，但计算时间会增加。理解这些权衡很关键。

残差反映了控制方程（如Navier-Stokes方程）的"满足程度"，属于迭代误差的指标。但即使残差降到10^-6，如果离散化误差很大，整体精度仍然很低。实践中，应该同时关注残差收敛和监测点的物理量（如翼的阻力系数Cd）是否达到稳定状态。Fluent和STAR-CCM+都有这种监测值历史绘图功能。

离散化方案决定了离散化误差的"阶"。一阶迎风的误差与网格尺寸h成正比(O(h))，二阶迎风与h²成正比(O(h²))。也就是说，网格细化为原来的1/2时，一阶方案的误差约为1/2，二阶方案约为1/4。在网格收敛性研究中，高阶方案通常以更少的网格数达到高精度。但在解的梯度陡峭的区域，可能会发振，这时需要加数值粘性等措施。

首先，明确感兴趣的输出量(Quantity of Interest: QoI)。不是"想知道强度"，而是要具体化为"螺栓连接部第一主应力的最大值"这样的表述。其次，选择对该QoI影响最大的参数（通常是网格尺寸），计划参数研究。例如，以5mm为基础网格尺寸，再做1.5倍、0.75倍、0.5倍三个变工况，共4个算例，将结果变化绘制成图表。

这正是关键所在——局部细分化。先用第一遍计算结果，找出应力集中部、流速或温度梯度大的区域。使用Abaqus/CAE的"网格控制"或Ansys Meshing的"Face Sizing"，只在感兴趣区域把网格尺寸分别改为1/2、1/4等。对于理论上的奇异点（如孔边缘），无论多细网格应力都会发散，此时需要改用应力集中系数等别的评价指标。

根据行业不同，实践方式多种多样。汽车工业中多数遵循内部检查清单。至少应记录这些项目：使用的软件及版本（如Abaqus 2023）、单元类型和积分点数（如C3D8R：8节点六面体低减积分）、网格统计（单元数、节点数）、收敛判据和结果（如米塞斯应力最大值相对变化≤2%）、最后采用的网格方案的依据。这些对可追溯性至关重要。

有的，近年来越来越先进。Ansys Workbench配备"Convergence"工具，对指定QoI自动执行h细分法或p细分法（提高单元阶数），并用图表显示结果收敛情况。COMSOL Multiphysics可设置"网格序列"研究步骤，能在一系列不同网格上连续运行参数扫描式的分析。Abaqus/CAE本身没有太多自动工具，但用户通常编写Python脚本来改变网格尺寸并批量运行。

非常重要。低质量的单元会显著增加离散化误差，有时甚至导致求解器发散。各软件都有强大的检查功能。例如，Siemens NX Nastran的"Element Quality Check"能一并评估Jacobian比、翘曲量、宽度比等，并将超过阈值的单元用红色标记。一般目标是：宽度比≤10，倾斜角≤60度。对于难以避免的低质量单元所在区域，需要进行敏感性分析（将其结果与良好网格区域比较）来评估影响。

主要两点。第一，"代码验证"信息有限。商业软件（如Ansys、Dassault）投入巨资进行基准测试和代码验证，但开源软件的这部分历史往往不清楚。第二，默认设置的可信度。OpenFOAM非常灵活，但也意味着用户要自己正确选择离散化方案和求解器设置。因此，必须用已知的精确解或标准问题（如NASA的后流圆柱问题）来自我验证求解器设置。

这可能是"网格锁定"或"数值不稳定"的迹象。网格锁定发生在完全积分的一阶单元（C3D8）处理接近不可压缩材料（泊松比接近0.5）时，单元会表现出过度刚性。对策是使用低减积分单元(C3D8R)或混合单元。数值不稳定在CFD中用高阶迎风差分且网格过细时会出现。此时需要调整求解器的松弛系数，或改用略偏一阶的方案。

这是个问题，说明粗网格的空间分解能不足以捕捉物理现象。在非线性问题（如接触、材料非线性）中特别明显。例如，Abaqus接触分析中，粗网格会导致接触面"跳过"，产生过大贯入并发散。从验证的角度，收敛的细网格结果被视为"更准确"，但如果其计算成本不可承受，可以选择中等网格在确保稳定收敛的前提下，作为采用方案，并在最终报告中说明理由。

这是常见的重要现象。位移是全局行为，相对粗糙的网格就能收敛；应力是导数（位移的梯度），尤其在应力集中部需要非常细的网格。判定标准是哪个QoI对设计最关键。如果目标是强度设计，优先考虑最大应力的收敛性。此时不是全局细化，而是局部细化应力集中区。最终报告要写明："位移在5mm网格收敛，应力在2mm网格收敛，故应力评估采用2mm网格"。