奇异矩阵错误
奇异矩阵
老师，出现了「Singular matrix」错误。
奇异矩阵错误的理论基础
奇异矩阵的数学意义
奇异矩阵错误是什么意思？「奇异」具体指什么状态？
在求解线性方程组
「机构」状态具体是什么？什么样的模型会出现这种情况？
典型情况是约束不足。例如，只建立一根梁的三维模型，两端都没有完全固定。梁在空间中对所有6个自由度（3个平移和3个旋转）都没有刚度。用Abaqus/Standard运行这个模型，会显示"NUMERICAL SINGULARITY"警告，位移会变成1.0E+30这样的非物理值。
材料性质或单元类型也会导致这个问题吗？
当然会。如果杨氏模量设置成0或极小值（比如1.0E-20 Pa），材料本身就没有刚度。另外，在壳单元中使用"膜单元"（只有弯曲刚度，没有面内刚度）也会导致面内方向变成机构，引起矩阵奇异。
奇异矩阵错误的数值计算方法
求解器如何检测奇异性
软件如何检测奇异矩阵？每次都计算行列式吗？那样计算量很大啊。
完全正确。对大规模矩阵，不会直接计算行列式。直接法求解器（如Ansys的稀疏直接求解器）在进行LU分解时，监测对角元素（主元素）。当主元素绝对值降到
迭代求解器（比如共轭梯度法）也会出现这个错误吗？
会出现，但表现形式不同。共轭梯度法（CG法）要求矩阵是正定对称的。奇异矩阵不满足正定性，所以迭代无法收敛，残差也不会减小。通常会报错"非正定矩阵"而不是直接说"奇异"。在Ansys的PCG求解器中，会显示"Pivot value less than tolerance"这样的信息。
我听说奇异矩阵有时候也能得出解。为什么？
这是因为用了伪逆矩阵。对奇异方程
奇异矩阵错误的实际应用
出现错误时的系统检查步骤
出现奇异矩阵错误时，应该先检查什么？
首先检查"边界条件的完整性"。对三维实体模型，需要确保6个自由度（3个平移和3个旋转）都被约束了。不能只用简支约束，要检查是否还有剩余的刚体运动。在Ansys Mechanical中，关闭"弱支撑"选项，隐藏的约束不足就会显露出来。
边界条件看起来没问题，还是出错了。下一步呢？
下一步检查"单元之间的耦合状态"。当不同类型单元（比如实体单元和壳单元）共享同一个节点时，旋转自由度可能无法传递，形成铰链式的机构。在Abaqus中，应该用"*TIE"约束或"*COUPLING"来正确连接，或用MPC（多点约束）来耦合相关自由度。
接触分析中经常出现这个错误。接触引起的奇异性怎么处理？
接触面在初始状态下相距太远或有初始贯穿，都会导致问题。在Ansys中，用"Adjust to Touch"功能调整初始接触。另外，接触单元的法向刚度（法向支撑）不能太小，经验上应该取基础实体单元刚度的10～100倍。
网格本身会不会是问题原因？
会的。"坏的单元形状"是常见原因。比如一次单元中内角极小的三角形壳单元（小于1度），或宽高比超过1000的细长单元，其雅可比行列式会变成奇异的，单元刚度矩阵计算不稳定。用网格质量检查工具（如"Skewness"和"Aspect Ratio"），确保在ISO标准或软件推荐范围内。
奇异矩阵错误的软件比较
各软件的错误消息和对策工具
Ansys、Abaqus、COMSOL处理奇异矩阵的方式有区别吗？
有明显区别。Ansys Mechanical显示"A singular matrix has been encountered..."，默认启用"弱支撑"，自动尝试继续计算。Abaqus/Standard则显示"THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES..."或"NUMERICAL SINGULARITY"，态度很严格，要求用户手动修正"*BOUNDARY"才能继续。COMSOL简洁地显示"Singular matrix."，经常建议自动网格自适应。
Ansys的"弱支撑"是什么？能否依赖它？
就是在机构方向上自动加入微小的支撑刚度。默认开启，能让计算继续进行，但这不是真正的解决方案。加入的虚假刚度可能影响结果精度。特别是在固有值分析或座屈分析中，绝对不能依赖弱支撑。正确做法是根据物理现实设置边界条件，不依赖虚假的弱支撑。
在Abaqus中，怎样定位奇异性的位置？
查看.msg文件（消息文件）。错误信息附近会明确指出"THE ROTATION ABOUT AXIS 1 AT NODE 1005 INSTANCE PART-1-1 IS NOT CONSTRAINED."这样的信息，告诉你哪个自由度（平移或旋转）、哪个节点、哪个部品没被约束。根据这个信息回到前处理模块修改约束。
COMSOL多物理耦合，会不会因为这个原因导致奇异性？
经常发生。比如把"结构力学"和"静电场"耦合时，结构位移有定义，但电位没有设置参考点（接地），电场控制方程就变成奇异的。每个物理接口都要独立地检查边界条件是否完整。COMSOL的"模型生成器"会在有问题的物理接口图标上显示警告标记，你可以从那里找到原因。
奇异矩阵错误的故障排除
高级案例和最后手段
所有检查都做了，非线性分析在某一步突然出现奇异矩阵错误。怎么回事？
这是"结构失稳"的表现。大变形导致支座失效、接触面分离导致约束消失、材料局部屈服且刚度急剧下降（应变软化），都会这样。Ansys的非线性求解器（Newton-Raphson法）每一增量步都要重新计算切线刚度矩阵，这里就会检测到奇异性。对策是减小增量步长，或改用弧长法（Riks法）。
用螺栓连接的多部件模型，只有螺栓孔周围出现奇异性。什么原因？
螺栓孔周围的网格可能出现"沙漏"模式（低减积分单元容易出现）。对策有三个：1) 在螺栓孔周围使用扫掠网格，改善单元形状。2) 改用二次单元或完全积分单元（如C3D20），替代低减积分单元（C3D8R）。3) 用刚体单元或梁单元模拟螺栓本体，用MPC进行连接。参考JIS B 1083等螺栓连接体分析指南。
找不到原因时，有最终的排查方法吗？
用"删除法"彻底简化模型。先在线性静力分析中重现错误，然后把单元类型统一成最简单的（3D用一阶完全积分六面体），材料用等向弹性，边界条件简化为单点完全固定加单点集中荷载。这个最小模型成功后，逐步加回部件、接触、复杂材料等。哪里出错就找到病根了。虽然费时，但最稳妥。
完全消除奇异矩阵错误有什么设计诀窍吗？
在模型建立初期就做"刚体模式检查"。许多商用软件都有这个功能，在运行前检查边界条件合理性。比如MSC Nastran的"ASET"和"OMIT"卡都包含刚体模式检查。另外，从设计上避免不切实际的约束（比如，现实中有微小间隙的地方不要用完全面接触）。记住：CAE是现实的近似，要时刻保持对真实物理的敬畏。
更详细
错误