BLAS/LAPACKエラー

BLAS (Basic Linear Algebra Subprograms) は、線形代数の基本演算を標準化したライブラリ規格だ。例えば、連立一次方程式

その通りだ。Ansys MechanicalやAbaqus/Standardの直接法ソルバー（スパースソルバー）も、内部ではPARDISOやMUMPSといった並列ソルバーが使われ、それらがさらにBLAS/LAPACKに依存している。例えば、Intel Math Kernel Library (MKL) や OpenBLAS はBLAS/LAPACK規格に準拠した実装で、CPUのSIMD命令（AVX-512など）を活用して極限まで計算を高速化している。

典型的なのは、特異行列や条件数が極端に悪い（ill-conditioned）行列に対する分解処理だ。例えば、LAPACKの `dgesv` (倍精度連立一次方程式ソルバー) はLU分解を行うが、ピボット要素がマシンイプシロン（倍精度で約

条件数κは、行列Aの最大特異値と最小特異値の比

直接法（Abaqus/Standardのデフォルトソルバーなど）で顕著だ。なぜなら、LU分解のような正確な行列分解を行う過程で数値的特異性が露見するからだ。一方、反復法（Abaqus/StandardのIterative SolverやAnsysのPCGソルバー）は収束しない、または発散する形で問題が現れる。ただし、反復法の前処理にILU分解などを使う場合、その内部でBLAS/LAPACKが呼ばれエラーになることもある。

それはBLASルーチンの引数エラーだ。BLAS `dgemm` は一般行列の乗算

以下の順でシステマティックに調査する。1) ログファイルの前後を確認：エラー直前の「行列の条件数が大きすぎます」などの警告を見逃さない。2) 単位系の一貫性：特にインポートしたジオメトリやマテリアルプロパティをチェック。ヤング率2e11[Pa]とジオメトリ長さ0.001[mm]の組み合わせは桁差が大きすぎる。3) 拘束条件：全ての剛体移動モード（並進・回転）が拘束されているか。特に対称条件やローカル座標系に注意。

有効なのはサブストラクチャリング（部分モデル解析）だ。Abaqusなら*SUBSTRUCTURE GENERATE、AnsysならSuperelementを使って、疑わしいコンポーネントだけの剛性マトリクスを生成し、その特異性をテストする。あるいは、ソルバー設定を変えてみる。Ansys Mechanicalで直接法ソルバー（Sparse）が失敗したら、反復法ソルバー（PCG）に切り替えてみる。反復法で収束すれば、問題が数値的な特異性（直接法が敏感）なのか、物理的な不安定性（どちらでも失敗）なのか切り分けられる。

直接法ソルバーには「ピボット許容値」や「特異性許容値」を設定できる場合がある。Ansys Sparseソルバーの `PIVCHECK` や、Abaqusの `*CONTROLS, RESET` カードで `singularity criterion` を緩和できる。ただし、これは根本解決ではなく、数値的にゼロと見なす閾値を

各社とも異なるメッセージでラップしている。Ansys Mechanicalは「A failure occurred while calling a BLAS routine from the sparse solver」と明確に表示する。Abaqus/Standardは「THE MATRIX IS SINGULAR」や「TOO MANY ATTEMPTS MADE FOR THIS INCREMENT」が多く、内部エラーコード「-9」などが.datファイルに記録される。COMSOLは「Failed to evaluate Jacobian」や「Singular matrix」と表示し、線形ソルバーの詳細ログでBLAS/LAPACKエラーを追跡できる。

可能な場合もあるが、高度な作業だ。COMSOLや一部のオープンソースFEMコード（CalculiXなど）はリンクするBLASライブラリを変更できる。例えば、Reference BLASからIntel MKLやOpenBLASに切り替えると、異なる数値アルゴリズムやピボット戦略が使われ、エラーを回避できる「可能性」はある。しかし、根本原因がモデル側にあるなら、ライブラリを変えても別の形で失敗する。AnsysやAbaqusはプリコンパイル済みのMKLをバンドルしており、通常は変更できない。

MSC NASTRANは独自の行列演算モジュール（DMAP）を持つが、内部でBLAS相当の演算をしている。エラーは「USER FATAL MESSAGE 3060 (SSG3)」のような形で、「SINGULAR MATRIX」と出る。LS-DYNAの implicit solver も同様だ。興味深いのは、Code_Aster（Salome-Meca）のようなオープンソースソルバーで、エラーメッセージが非常に詳細で、問題の発生した具体的なLAPACKルーチン名（例：`dgetrf_`）と引数値をログにダンプしてくれる点だ。これは根本原因のデバッグには非常に有用だ。

ある。大規模モデルで直接法ソルバーを使用する場合、LU分解の作業用メモリ（ワークアレー）が不足すると、メモリ割り当てに失敗し、不正なポインタがBLASルーチンに渡されることがある。結果として「BLASルーチン呼び出し失敗」というエラーに帰着する。対策は、Ansysなら `SPARSE_MEMORY_FACTOR` を増やす、Abaqusなら `*MEMORY` 設定を見直す、あるいはソルバータイプを反復法（メモリ使用量が少ない）に変更することだ。

接触状態の変化が剛性マトリクスのランクを変えるからだ。例えば、あるステップまでは接触面がスティック状態で剛性が高かったが、次のステップでスリップに遷移し、接線方向の剛性が実質的にゼロになる。この急激な変化でマトリクスが数値的に特異になる。対策は、接触剛性（`KN`）を適切に設定し、増分ステップ `INC` を小さくして状態変化をなめらかにすること。Abaqusでは `*CONTACT CONTROLS` で安定化をオンにできる。

数値的な「ごまかし」になるが、以下の手順がある。1) モデル全体に微少な人工剛性を追加する。Abaqusの `*DASHPOT` 要素（質量比例ダンパー）を全体に定義し、係数を非常に小さく（