零主元错误

主元是指在「直接法」求解器（例如高斯消元法）中的对角元素，是一个关键的数值。求解器使用这个主元来消除其他行。当主元变为零或非常小（考虑到计算机舍入误差，例如小于1e-12），就会接近零除以零的情况，导致解发散或计算无法继续。这在数值上表明刚度矩阵

最典型的是「刚体位移模式」没有被约束的情况。例如，在三维结构分析中，如果物体在任何方向都没有被固定，而是施加荷载，那么物体可以无限运动。此时刚度矩阵没有逆矩阵（奇异），主元会变为零。另一个例子是材料常数或单元设置错误导致某个方向的刚度在理论上为零的情况。

是的，会出现。在热传导分析中，如果所有边界都是绝热的，且内部没有热源，温度就无法唯一确定（只有常数不确定性）。这是热学意义上的「刚体模式」。在流体分析中（特别是压力泊松方程），如果没有指定压力参考点（例如某一点的压力值），压力就会有常数位移的不确定性，这会导致零主元问题。

直接法求解器有「主元选择」策略。例如，部分主元选择会在当前列中找到绝对值最大的分量，交换行将其设为主元。Ansys Mechanical的稀疏直接求解器（默认）高度执行此操作。但如果主元降到例如

完全正确。求解器在进行行列消除时，会报告最后保留下来的（无法消除的）具有奇异性的自由度。在这个例子中，节点12345的Z方向位移（UZ）相关的方程在所有其他方程中是独立的，其系数极小或为零。这意味着该节点的Z方向在物理上或数值上没有被约束，提供了强有力的线索。

是的，完全正确。Abaqus/Standard的「人工刚度」或某些求解器的「弱弹簧」选项，会自动向奇异自由度添加极小的刚度（例如矩阵最大对角项的

按照系统的检查清单效率最高。第一步是「边界条件」。关注错误消息中的自由度（UX、UY、UZ、ROTX等），确认那个方向的位移或旋转在所有部件中是否被适当约束。第二步是「接触条件」。检查是否发生了接触分离「抬起」，初始接触状态是否合适。第三步是「单元特性」。壳单元的面外刚度或梁单元的旋转自由度是否合适。

Ansys Mechanical有「刚体模式检查」功能。在添加了弱弹簧的状态下运行一个静力分析步骤，将位移结果的云图「用户定义」设置为非常小的范围（例如0.001 mm）。异常大的位移（刚体运动）显示的区域就是约束不足的区域。Abaqus也有类似的方法，可以施加非常小的诊断荷载来可视化位移。

常见的是各向异性材料定义错误。例如，复合材料各层方向错误，或某个方向的杨氏模量被误输入为接近零的值（如1e-10），就会导致该方向没有刚度，成为零主元的原因。此外，超弹性材料模型（Mooney-Rivlin等）在外推到未定义应变范围时，数值上也会变得不稳定。

是的，各有特点。Ansys Mechanical的默认直接求解器检测到奇异性时，会以「零主元」错误停止，并详细报告节点号和自由度。Abaqus/Standard也类似，但也使用「数值奇异（numerical singularity）」这样的消息。COMSOL Multiphysics报告「奇异矩阵」，且往往在求解器设置中有「忽略奇异性继续」的选项。但根本原因都是一样的。

现象会改变，但问题不会解决。直接法求解器（Ansys Sparse、Abaqus Direct）会早期检测到奇异性并报错停止。而前处理共轭梯度法（PCG）等迭代法求解器会继续尝试计算，但无法满足收敛标准（残差范数），最终表现为不收敛（或非常慢）。只要模型根本上约束不足，就无法得到正确的解。

MSC NASTRAN会显示「秩不足（RANK DEFICIENT）」或「奇异矩阵（SINGULAR MATRIX）」这样的致命错误消息。LS-DYNA（显式方法）由于求解包含惯性项的动力方程，只要质量矩阵正则，即使刚度矩阵奇异计算也能继续，但会产生不自然的刚体运动，结果物理上无意义。两种情况都需要重新检视边界条件。

常见陷阱是「重复节点」或「微小间隙」。两个部件看起来接触，但网格没有共享，又没有适当定义接触条件时，会在物理上保持分离。另外，使用壳单元或梁单元时，面内被约束但「面外旋转（ROTZ）」没有被约束也会导致奇异。梁的旋转自由度默认通常是释放的，需要特别注意。

非线性分析中，每个增量步、每次迭代都会重新构造切线刚度矩阵。随着分析进行，接触状态变化导致约束丧失（抬起）、材料极度软化、单元过度变形导致数值不稳定，这些都会导致矩阵在中途变奇异。特别是含摩擦的接触，约束频繁开关，容易触发这个错误。

物理上有意义的「弱弹簧」添加是最后手段。例如，在宇宙空间中自由浮动的卫星姿态控制分析，无法约束重心的平移和旋转。这种情况下，可以通过添加与质量矩阵成比例的极小阻尼（人工阻尼），或虚拟添加非常柔软的弹簧来抑制刚体模式。但必须评估这个弹簧刚度对结果的影响，例如通过固有值分析确认刚体模式的频率足够低。