病态矩阵（ill-conditioning）
病态矩阵
老师，"ill-conditioned matrix"是什么？
病态矩阵（ill-conditioning）的理论基础
病态矩阵的定义和影响
病态矩阵（ill-conditioning）具体指什么样的矩阵状态呢？不是笼统的"条件差"的感觉，而是数值上如何判断呢？
很好的问题。在定量上，我们用"条件数（condition number）"
有效数字丢失在仿真结果中会表现为什么现象呢？只是收敛变慢吗？
收敛缓慢只是一个症状，最严重的是"解的不准确性"。例如，计算悬臂梁顶端挠度时，理论值是10.0 mm，但由于病态矩阵，求解器返回的值可能不在9.999 mm到10.001 mm范围内，而是9.5 mm或10.5 mm这样完全错误的值。对于迭代求解法，要么无法收敛，要么"收敛"到错误的解。即使是直接法，数值舍入误差也会被放大到整个解中。
条件数变大的物理原因是什么？只有材料常数差异极端时才会出现吗？
材料常数的极端差异是典型例子。比如，将橡胶（杨氏模量~1 MPa）和钢（~210 GPa）结合，刚度比可达
病态矩阵（ill-conditioning）的数值计算方法
求解器内部的处理
求解器内部如何处理病态矩阵？条件数大时是否会自动检测？
大多数商用求解器都具备某种检测功能。例如，Abaqus/Standard的直接法求解器在进行主元选择时，会检测到主元值极小的单元，并输出警告信息（如"***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY..."），这实际上是病态的征兆。对于迭代法求解器，当迭代次数达到设定上限仍未收敛（"残差范数"不下降）时，会向用户提示问题。
直接法和迭代法对病态的健壮性有区别吗？
区别很大。直接法（LDL^T分解或LU分解）可以通过主元选择（缩放和重排）在一定程度上应对病态，但根本上条件数并未改变。最终解的精度受舍入误差影响。而带前处理的共轭梯度法（PCG）等迭代法则取决于前处理（Preconditioner）的质量。例如，使用不完全Cholesky分解作前处理，如果原矩阵病态，前处理本身可能变得不稳定。Ansys的稀疏直接法默认执行强缩放，但也有其限制。
"缩放"具体是在做什么？是改变单位吗？
这是与单位统一不同的数值处理。矩阵
病态矩阵（ill-conditioning）的实际应用
建模时的预防措施
实际建模中是否有具体的检查清单来避免病态矩阵？
至少应确认以下4点。
如果必须连接橡胶和钢这样物性差异大的材料，应该怎么办？
有几种实用对策。第一，使用"子模型方法"。在整体模型中将刚体部分视为刚体，局部的详细应力评估在子模型中进行。第二，通过"多点约束（MPC）"或"连接单元"进行连接，避免单元直接结合。第三，更改求解器设置。在Abaqus中，尝试直接法求解器加参数`PARAMETER, INREL, -2`（应用更强的缩放）。但这不是根本解决方案，结果的合理性必须通过其他方式验证。
分析过程中出现病态警告时，最先应该检查模型的哪个部分？
如果警告信息输出了节点号或单元号，那就是震源所在地。如果没有输出，按以下顺序调查。
病态矩阵（ill-conditioning）的软件对比
各求解器的行为和对策
Ansys Mechanical和Abaqus/Standard在处理病态矩阵的错误信息和默认对策上有区别吗？
区别很明显。**Abaqus/Standard**相对"爱管闲事"，会频繁输出数值特异性（numerical singularity）警告，提示可能存在病态或刚体模式。而**Ansys Mechanical**的默认直接法求解器（Sparse Direct）内部会自动进行缩放和主元选择，所以对同一模型的警告频率较低。反面来说，可能会出现无法收敛或明显错误的结果而后结束分析。在Ansys中启用详细求解器输出（`SOLVE`命令的`/STATUS`）可以查看主元信息等细节。
听说COMSOL Multiphysics这样的多物理场软件，由于物理场耦合容易产生病态。这是真的吗？
完全正确。在COMSOL中进行结构-声学耦合或压电分析时，支配方程的矩阵容易自然产生病态。例如，将固体位移（量级：`1e-6 m`）与声压（量级：`1 Pa`）直接耦合时，变量间的尺度差异巨大。COMSOL默认具备**变量缩放**功能，对各物理场的变量应用内部缩放系数，以改善矩阵条件数。用户可以选择"分离求解"，对不同物理场应用不同的求解器设置（网格或时间步长），从而缓解问题。
无费开源求解器（CalculiX、Code_Aster）与商用软件相比，病态问题是否更明显？
求解器算法的核心部分差异不大，但**前处理**和**错误处理**的完成度有区别。CalculiX（CCX）的默认迭代法求解器对病态的容错能力较弱，常常无法收敛。此时推荐手动切换到直接法求解器（SPOOLES）。Code_Aster提供了大量求解器和前处理选项，可选择`LDLT`分解或`MUMPS`直接法求解器。但选择恰当的配置需要深度知识，不能指望商用软件那样的自动调优。无论如何，修正建模基础（单位、网格、约束）是最重要的。
病态矩阵（ill-conditioning）的故障排除
具体的错误和对应对策
Abaqus中出现"***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE..."警告。我已经约束了刚体模式。下一步应该做什么？
首先，重点调查警告信息中输出的**节点号**周围区域。检查该节点相连单元的网格质量（变形程度）。然后确认该节点受到的**荷载或接触**。如果只有该节点承受大的集中荷载而支撑较弱，可能产生极大的局部位移，被判为数值特异。对策包括：1) 将集中荷载分散为面压力，2) 细化并均匀化该节点周围网格，3) 在该自由度方向加入极小刚度的弹簧（`1e-6`量级）以增强数值稳定性。
Ansys的非线性分析（有接触）无法收敛，出现"pivot value"相关错误。这是病态造成的吗？
几乎可以肯定是病态或其极限情况"特异"。非线性分析的每个增量步都会更新切线刚度矩阵，接触状态急剧变化（分离→接触）时，矩阵条件数瞬间恶化。Ansys的"pivot value"错误意味着直接法求解器检测到负值或极小的主元。对策如下。
没时间重新建网格。对现有略显病态的模型，有什么简便方法验证结果的可信度吗？
有两种"简便"验证方法。第一，**求解器切换测试**。例如，在Abaqus中从默认直接法切换到迭代法（`ITERATIVE`）重新计算。若两者结果（反力、最大位移、最大应力）在小数点后2~3位一致，说明病态影响较小。第二，**扰动测试**。将杨氏模量或荷载值微调0.1%后重新分析。健全的非病态模型，结果变化应与输入变化成比例（保持线性性）。如果0.1%的输入改变导致10%的结果变化，结果就不可信。两者都不是根本解决，但对风险评估有用。
帮助
更详细
错误