ゼロピボットエラー
ゼロピボット
先生、「Zero pivot」は特異行列と同じですか？
理論と物理
ゼロピボットの物理的意味
「ゼロピボットエラー」というエラーメッセージが出ました。ピボットって何ですか？なぜ「ゼロ」だと問題なんですか？
ピボットとは、連立一次方程式を解く「直接法」ソルバー（例えばガウスの消去法）において、対角成分を指す重要な数値です。ソルバーはこのピボットを使って他の行を消去していきます。ピボットがゼロ、あるいは計算機の丸め誤差を考慮して非常に小さい値（例えば1e-12以下）になると、ゼロ除算に近い状態となり、解が発散したり、計算が続行できなくなります。これは剛性マトリクス
剛性マトリクスが特異になる物理的な状況とは、具体的にどんな時ですか？
代表的なのは「剛体移動モード」が拘束されていない場合です。例えば、3次元の構造解析で、物体にどの方向の変位も拘束（固定）せずに荷重をかけると、物体は無限に動けてしまいます。この時、剛性マトリクスは逆行列を持たず（特異）、ピボットがゼロになります。別の例は、材料定数や要素の設定ミスで、ある方向の剛性が理論的にゼロになる場合です。
熱伝導解析や流体解析でも同じエラーが出ますか？その場合の「剛体移動」に相当するものは何ですか？
はい、出ます。熱伝導解析では、全ての境界が断熱（熱流束ゼロ）で、内部発熱もない場合、温度が一意に定まりません（定数だけの不定性）。これが熱的な「剛体モード」です。流体解析（特に圧力のポアソン方程式）では、圧力の基準点（例えばある1点の圧力値）を指定しないと、圧力が定数だけシフトできる不定性が生じ、これがゼロピボットの原因になります。
数値解法と実装
ソルバーアルゴリズムとピボット選択
ソルバーはピボットが小さいとどう処理するんですか？エラー以外の選択肢はないんですか？
直接法ソルバーには「ピボット選択」という戦略があります。例えば、部分ピボット選択では、現在の列の中で絶対値最大の成分を見つけ、行を交換してそれをピボットとします。Ansys Mechanicalのスパース直接ソルバー（デフォルト）はこれを高度に行います。それでもピボットが、例えば
メッセージに「at degree of freedom UZ of node 12345」のように具体的な節点と自由度が表示されることがあります。これはソルバーが特異性を検出した「場所」を示しているんですか？
その通りです。ソルバーが行・列の消去を進め、最後に残った（消去できなかった）特異性のある自由度を報告しています。この例では、12345番節点のZ方向変位（UZ）に関する方程式が、他の全ての方程式から独立で、係数が極端に小さいかゼロだったことを意味します。これは、その節点のZ方向が物理的・数値的に拘束されていない、という強力な手がかりになります。
「弱いばね」や「人工剛性」というオプションを聞いたことがあります。これはゼロピボットを回避する数値的処方箋ですか？
はい、その通りです。Abaqus/Standardの「人工剛性」や、いくつかのソルバーにある「弱いばね」オプションは、特異な自由度にごく小さな剛性（例えばマトリクス最大対角項の
実践ガイド
モデルチェックのワークフロー
ゼロピボットエラーが出た時、最初に何をチェックすべきですか？具体的な手順を教えてください。
体系的なチェックリストに従うのが効率的です。第一に「境界条件」。エラーメッセージの自由度（UX, UY, UZ, ROTXなど）に注目し、その方向の変位や回転が全てのパートで適切に拘束されているか確認します。第二に「接触条件」。接触が離れることで拘束が失われる「リフトオフ」が起きていないか、初期接触状態は適切か。第三に「要素特性」。シェル要素の面外剛性や、ビーム要素の回転自由度が適切か。
複雑なアセンブリモデルで、どこが拘束不足か目視で見つけるのが難しいです。効率的な検出方法はありますか？
Ansys Mechanicalでは「剛体モードチェック」機能があります。弱いばねを追加した状態で静的解析を1ステップ実行し、変位結果のコントアを「ユーザー定義」で非常に小さな範囲（例えば0.001 mm）に設定します。異常に大きな変位（剛体移動）を示す部分が、拘束不足の領域です。Abaqusでも同様に、非常に小さな荷重（診断荷重）をかけて変位を可視化する方法があります。
材料特性や実定数の設定ミスでも起きるとのことですが、具体的にはどんなミスがありますか？
よくあるのは、異方性材料の定義ミスです。例えば、複合材の層ごとの方向を間違えたり、ある方向のヤング率を誤ってゼロに近い値（1e-10など）で入力してしまうと、その方向の剛性がなくなり、ゼロピボットの原因になります。また、超弾性材料モデル（Mooney-Rivlinなど）で未定義のひずみ領域に外挿された場合、数値的に不安定になることもあります。
ソフトウェア比較
各ソルバーの挙動と対策ツール
Ansys、Abaqus、COMSOLで、ゼロピボットエラーへの対処方法やメッセージの出方は違いますか？
はい、特徴があります。Ansys Mechanicalのデフォルト直接ソルバーは、特異性を検出すると「ゼロピボット」エラーで停止し、節点番号と自由度を詳細に報告します。Abaqus/Standardも同様ですが、「数値特異（numerical singularity）」というメッセージも使います。COMSOL Multiphysicsは「特異マトリクス」と報告し、ソルバー設定で「特異性を無視して続行」するオプションを持っていることが多いです。いずれも根本原因は同じです。
ソルバーの種類（直接法 vs 反復法）を変えるとエラーは出なくなることがありますか？
現象は変わりますが、問題は解決しません。直接法ソルバー（Ansys Sparse, Abaqus Direct）は特異性を早期に検出してエラー停止します。一方、前処理付き共役勾配法（PCG）などの反復法ソルバーは、特異な問題に対しても計算を続けようとしますが、収束基準（残差ノルム）を満たせず、収束しない（または非常に遅い）という形で問題が現れます。モデルに根本的な拘束不足がある限り、正しい解は得られません。
NASTRANやLS-DYNAのような別のソルバーではどう扱われますか？
MSC NASTRANでは「ランク不足（RANK DEFICIENT）」や「特異マトリクス（SINGULAR MATRIX）」という致命的エラーメッセージとして現れます。LS-DYNA（陽解法）は、慣性項を含む動的方程式を解くため、質量マトリクスが正則であれば剛性マトリクスが特異でも計算は進みますが、不自然な剛体運動が発生し、結果が物理的に無意味になります。いずれも境界条件の見直しが必須です。
トラブルシューティング
ケース別対策と落とし穴
「全ての境界をきちんと固定したはずなのに、ゼロピボットエラーが出る」という場合、他にどんな原因が考えられますか？
よくある落とし穴が「重複節点」や「わずかな隙間」です。2つのパーツが接触しているように見えても、メッシュが共有されておらず、かつ接触条件が適切に定義されていないと、物理的に分離したままになります。また、シェル要素やビーム要素を使用する場合、面内は拘束されていても「面外回転（ROTZ）」が拘束されていないために特異になるケースがあります。ビームの回転自由度は、デフォルトで解放されていることが多いので要注意です。
非線形解析（接触や大変形）の途中でゼロピボットエラーが出ました。線形解析では出なかったのに、なぜですか？
非線形解析では、各増分ステップ、各反復で接線剛性マトリクスが再構成されます。解析が進むうちに、接触状態が変化して拘束が失われたり（リフトオフ）、材料が極端に軟化したり、要素が過度に変形して数値的に不安定になったりすることで、途中でマトリクスが特異になることがあります。特に、摩擦を含む接触では、拘束のオン/オフが頻繁に切り替わるため、このエラーが起きやすくなります。
どうしても拘束不足の部位が特定できない、または拘束を追加したくない場合の最終手段はありますか？
物理的に意味のある「弱いばね」の追加が最終手段です。例えば、宇宙空間で自由に浮遊する衛星の姿勢制御解析では、重心の並進・回転を拘束できません。その場合、剛体モードを抑制するために、質量マトリクスに比例したごく小さな減衰（人工減衰）や、非常に柔らかいばねを仮想的に追加する手法があります。しかし、そのばねの剛性が結果に与える影響を必ず評価し、例えば固有値解析で剛体モードが十分に低周波数であることを確認する必要があります。
なった
詳しく
報告