行列分解エラー

いい質問だ。単なる計算エラーではなく、多くの場合、モデルに物理的に「解けない」問題があることを示している。例えば、構造解析で剛性マトリクス

代表的な原因は三つだ。第一は、前述の拘束不足。第二は、材料定数の設定ミス。例えば、ヤング率を誤って

出る。物理は違えど、離散化して連立一次方程式

その通り、分解アルゴリズムの実行中だ。直接法ソルバー（Abaqus/Standardのデフォルトなど）は、通常、ピボット選択付きのLU分解を行う。この過程で、対角項（ピボット）の絶対値が、事前に定義された許容値より小さくなるとエラーを発する。例えば、MKL PARDISOソルバーでは、デフォルトのピボット許容値は

一時的な回避策に過ぎず、根本解決ではない。確かにAbaqusでは`*CONTROLS, RESET`カードでピボット許容値を変更できるが、これは数値的に特異な行列を無理やり分解させる行為だ。結果として、解の精度が著しく低下したり、数値的に不安定で物理的に無意味な結果（例えば、一部の節点が異常変位）を生むリスクが高い。まずはモデルの拘束や材料定数を見直すべきだ。

出る場合と出ない場合がある。反復法は行列を直接分解しないので「zero pivot」エラーそのものは出ない。代わりに、係数マトリクスが正定値でない場合、共役勾配法（CG）は収束せずに失敗する。あるいは、条件数が極端に大きい（ill-conditioned）場合、反復が進まないか、発散する。Ansys Mechanicalでスパース直接ソルバーの代わりにPCGソルバーを使った場合、特異モデルでは「Solver pivot warning」ではなく「Solution not converged」という形で問題が表面化する。

まずはエラーメッセージと警告メッセージを詳細に読むこと。Abaqusの`.msg`ファイルやAnsysの`.out`ファイルには、エラーが発生した節点番号や要素番号が記載されていることが多い。次に、以下のチェックリストを順に実施する：1) 全てのパーツに適切な拘束（支持）が定義されているか。2) 接触定義に隙間や重複はないか。3) 材料定数（ヤング率、密度、熱伝導率など）にゼロや非現実的な小さな値が入っていないか。4) メッシュに歪んだ要素（アスペクト比1000以上など）がないか。

それは非常に有益な情報だ。節点12345の第3自由度（多くの場合、Z方向の並進移動）に関連する剛性が、数値的にゼロであることを意味する。具体的には、その節点がZ方向に全く拘束されていない、または、その節点に接続されている全ての要素がZ方向の剛性を持たない（例えば、XY平面内のみのトラス要素）可能性が高い。その節点周辺の拘束条件と要素タイプを重点的に確認せよ。

接触設定が原因の可能性が高い。特に、「摩擦」を定義した接触面で、初期接触状態が「開いている」場合、接線方向の剛性が定義されず、特異性を生むことがある。まずは摩擦係数をゼロに設定して試す。また、Ansysの「Weak Springs」やAbaqusの「Stabilization」機能を試す前に、接触ペアの初期隙間を`*CONTACT INTERFERENCE`（Abaqus）や`CNOF`（Ansys）で適切に調整することを検討する。これらは数値的安定化であり、物理的には「仮想的なばね」を追加する行為だからだ。

エラーの本質は同じだが、メッセージと内蔵ツールに違いがある。Abaqus/Standardは「***ERROR: SOLVER PROBLEM. NUMERICAL SINGULARITY」や「ZERO PIVOT」と明確に表示し、`.msg`ファイルで問題の節点・自由度を特定できる。一方、Ansys Mechanicalは「A zero pivot was encountered...」と出し、`/SOLU`コマンドで`PIVCHECK, ON`として詳細を出力させる必要がある。対策面では、Ansysには「Inertia Relief」や「Weak Springs」オプションがGUIから簡単に適用できるが、Abaqusでは`*RIGID BODY`や`*COUPLING`を用いて剛体運動を制御する方法が一般的だ。

COMSOL Multiphysicsは、デフォルトで「定数」や「変数」を用いて自動的に未定義の自由度（圧力基準点など）を処理しようとするため、初歩的な特異性エラーは起こりにくい設計だ。しかし、ユーザーが物理場インターフェースをカスタマイズして矛盾した境界条件を設定すれば、同様に「特異マトリクス」エラーが発生する。LS-DYNAのような陽解法ソルバーは、質量マトリクス

ソルバーアルゴリズムの堅牢性という点では、商用ソフトの高度にチューニングされたMKLやPARDISOライブラリに比べると、オープンソースソルバーは厳しい条件で数値的に不安定になる可能性はある。しかし、エラーの根本原因はモデル自体にある。CalculiX（CCX）も「zero pivot」エラーを報告する。むしろ、商用ソフトが内部で自動的に行う「数値的安定化」をオープンソースでは手動で設定する必要があるため、ユーザーが問題の本質（拘束不足）に早く気づけるという側面もある。

以下のような複合的な原因が考えられる。第一に、「多重拘束」。一つの節点に対して、`*COUPLING`（Abaqus）や`CP`コマンド（Ansys）による拘束と、通常の境界条件（`*BOUNDARY`）が競合し、過剰拘束ではなく矛盾した拘束を生み、数値的特異性を引き起こすことがある。第二に、「ソルバーの数値スケーリング問題」。単位系が混在（長さはmm、力はN）したモデルでは、剛性マトリクスの係数のオーダーが

最も確実な方法は、一貫した単位系を使用することだ。SI単位系（m, kg, s, N）に統一するか、または広く使われるmm単位系（mm, tonne, s, N）に統一する。もしモデル修正が難しい場合は、ソルバーのスケーリングオプションを利用する。Abaqusでは`*CONTROLS, PARAMETERS=FIELD`でスケーリング係数を調整できる。Ansysの直接ソルバー（SPARSE）は自動スケーリングがデフォルトで有効だが、`EQSLV, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,