悪条件行列（ill-conditioning）
悪条件行列
先生、「ill-conditioned matrix」って何ですか？
理論と物理
悪条件行列の定義と影響
悪条件行列（ill-conditioning）って、具体的にどういう状態の行列を指すんですか？「条件が悪い」という感覚的な表現ではなく、数値的にどう判断するんですか？
良い質問だ。定量的には「条件数（condition number）」
有効数字が失われるとは、具体的にどういう現象がシミュレーション結果に現れるんですか？収束が遅くなるだけですか？
収束遅延は一症状に過ぎない。最も深刻なのは「解の不正確さ」だ。例えば、片持ち梁の先端たわみを計算したとき、理論値が10.0 mmなのに、悪条件行列によってソルバーが返す値が9.999 mmから10.001 mmの範囲に収まらず、9.5 mmや10.5 mmといった全く異なる値になることがある。反復解法では収束しないか、誤った解に「収束」してしまう。直接法でも、数値的な丸め誤差が解全体に増幅されて現れる。
条件数が大きくなる物理的な原因は何ですか？材料定数を極端に変えたときだけですか？
材料定数の極端な差は典型例だ。例えば、ゴム（ヤング率~1 MPa）と鋼（~210 GPa）を結合したモデルでは、剛性の比が
数値解法と実装
ソルバー内部での扱い
ソルバーは内部で悪条件行列をどう処理してるんですか？条件数が大きいと自動で検知してくれるんですか？
多くの商用ソルバーは何らかの検知機能を持つ。例えば、Abaqus/Standardの直接法ソルバーは、ピボット選択時にピボット値が極端に小さい要素を検出し、警告メッセージ（例: 「***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY...」）を出力する。これは実質的に悪条件の兆候だ。反復法ソルバーでは、反復回数が設定上限に達しても収束しない（「残差ノルム」が下がらない）ことで、ユーザーに問題を提示する。
直接法と反復法では、悪条件に対する頑健さに違いはありますか？
大きな違いがある。直接法（LDL^T分解やLU分解）は、ピボット選択（スケーリングと並べ替え）によってある程度の悪条件に対処できるが、根本的な条件数は変わらない。最終的な解の精度は丸め誤差に影響される。一方、前処理付き共役勾配法（PCG）のような反復法は、前処理（Preconditioner）の良し悪しが命運を分ける。例えば、不完全コレスキー分解を前処理に使うと、元の行列が悪条件だと前処理自体が不安定になることがある。Ansysのスパース直接法はデフォルトで強いスケーリングを行うが、それでも限界はある。
「スケーリング」とは具体的に何をしているんですか？単位系を変えることですか？
単位系の統一とは別の、数値的な処理だ。行列
実践ガイド
モデリング時の予防策
実際のモデリングで、悪条件行列を生まないための具体的なチェックリストはありますか？
最低限、以下の4点を確認せよ。
もし、ゴムと鋼のように物性差が大きい材料をどうしても結合しなければならない場合は、どうすればいいですか？
いくつかの実践的な対策がある。第一に、**サブモデリング**を使う。全体モデルでは剛体部分をリジッドボディとして扱い、局部の詳細な応力評価はサブモデルで行う。第二に、**マルチポイント拘束（MPC）** や**結合要素**を介して接続し、直接的な要素結合を避ける。第三に、ソルバー設定を変更する。Abaqusならば、直接法ソルバーで`PARAMETER, INREL, -2`（より強いスケーリングを適用）を試す。ただし、根本解決にはならないので、結果の妥当性は必ず別途検証すること。
解析中に悪条件の警告が出た場合、最初に疑うべきモデルの箇所はどこですか？
警告メッセージにノード番号や要素番号が出力されていれば、そこが震源地だ。出力されない場合は、以下の順で調査する。
ソフトウェア比較
各ソルバーの挙動と対策
Ansys MechanicalとAbaqus/Standardでは、悪条件行列に対するエラーメッセージやデフォルトの対策に違いはありますか？
違いは明確だ。**Abaqus/Standard**は比較的「おせっかい」で、数値的特異性（numerical singularity）の警告を頻繁に出力する。これは悪条件や剛体モードの可能性を示唆する。一方、**Ansys Mechanical**のデフォルト直接法ソルバー（Sparse Direct）は、内部で自動的にスケーリングとピボット選択を行うため、同じモデルでもAbaqusより警告が出にくい傾向がある。その代わり、収束しない、または明らかにおかしい結果で解析が終了することがある。Ansysで詳細なソルバー出力を有効に（`SOLVE`コマンドの`/STATUS`）すれば、ピボット情報などが確認できる。
COMSOL Multiphysicsのようなマルチフィジックスソフトでは、異なる物理場の結合によって悪条件が起きやすいと聞きました。本当ですか？
その通りだ。COMSOLで構造-音響連成や圧電解析を行う場合、支配方程式のマトリックスが自然と悪条件になりやすい。例えば、固体の変位（オーダー: `1e-6 m`）と音圧（オーダー: `1 Pa`）を直接結合すると、変数間のスケール差が巨大になる。COMSOLはデフォルトで**変数スケーリング**機能を持ち、各物理場の変数に内部スケーリング係数を掛けてマトリックスの条件数を改善しようとする。ユーザーは「分離解法」を選択し、物理場ごとに異なるソルバー設定（メッシュや時間ステップ）を適用することで、問題を緩和できる。
無償のオープンソースソルバー（CalculiX, Code_Aster）では、商用ソフトより悪条件の問題が顕著ですか？
ソルバーアルゴリズムの核心部分は大差ないが、**前処理**と**エラー処理**の完成度に差がある。CalculiX（CCX）のデフォルト反復法ソルバーは、悪条件に対して比較的脆弱で、収束しない場合が多い。その場合は、マニュアルで直接法ソルバー（SPOOLES）に切り替えることが推奨される。Code_Asterは非常に多くのソルバーと前処理オプションを提供しており、`LDLT`分解や`MUMPS`直接法ソルバーを選択できる。ただし、適切な選択には深い知識が必要で、商用ソフトのような自動チューニングは期待できない。いずれにせよ、モデリングの基本（単位、メッシュ、拘束）を正すことが最も重要だ。
トラブルシューティング
具体的なエラーとその対策
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%も結果が変わったら、その結果は信用できない。いずれも根本解決ではないが、リスク評価には有用だ。
関連トピック
なった
詳しく
報告