条件数 — CAE用语解释

直接相关。在求解FEM的刚性方程 [K]{u}={f} 时，矩阵K的条件数κ(K)越大，数值误差放大越严重。κ(K) = ||K|| * ||K^-1|| 定义，1.0为理想，超过1e6时需要谨慎对待数值稳定性。许多「精度不足」「不收敛」等FEM分析问题的根源往往都与条件数有关。

FEM中条件数恶化的原因主要有三个。①网格的宽高比极端（细长单元）——刚性在一个方向上极端高。②材料常数差异很大的异质接合——钢和橡胶混合模型中刚性比会达到10^5量级。③接近不可压缩材料（橡胶、软组织）——体积刚性与偏差刚性比值很大。这些都可以通过网格和材料设置来处理，关键是先确定根本原因。

直接法（LU分解）中舍入误差被放大，解变得不精确。迭代法（CG法、GMRES）收敛变慢或发散。例如κ=1e8的系统，双精度（有效数字15～16位）的有效数字会降到7～8位。实际中常见的问题是「位移看起来收敛了但应力值异常」——这是条件数差导致解失去精度的信号。

有几种方法。最有效的是预处理（前处理）。用不完全LU分解（ILU）或代数多重网格（AMG）作为迭代法的前处理，条件数会大幅下降，反复法收敛性改善。网格方面的对策包括：保持单元宽高比在2～5以下、异材边界处单元尺寸逐步变化等。Abaqus和Ansys默认做自动缩放，但大规模模型中能否自选预处理器是选择求解器的关键。

大多数FEM求解器会在日志中输出「max/min ratio」或「condition estimate」。Abaqus会输出warning: zero pivot 或 excessive distortion，OpenFOAM的fvSolution设置中solverPerformance日志可以看残差历史。Python中可以用np.linalg.cond(K)直接计算，但大规模矩阵计算本身很重，通常用估计器（RandomSVD等）。最快的办法是解析前检查网格质量报告中的宽高比。