ソルバー間比較：線形静解析 — 理論と検証フレームワーク
概要
先生！ 今日は線形静解析ソルバー比較の話なんですよね？ どんなものなんですか？
本記事では、主要な商用・オープンソースFEAソルバーを用いた線形静解析の横断的な比較検証を行う。NAFEMSベンチマーク問題群を共通の試験問題として使用し、各ソルバーの精度・効率・使いやすさを体系的に評価する。
ふむふむ…本記事ではって意外と身近な現象と繋がってるんですね。
対象ソルバー
具体的にはどんなアルゴリズムで線形静解析ソルバー比較を解くんですか？
| ソルバー | バージョン | 開発元 | ライセンス |
|---|---|---|---|
| MSC Nastran | 2023.1 | Hexagon | 商用 |
| NX Nastran | 2306 | Siemens | 商用 |
| Abaqus | 2023 | Dassault Systemes | 商用 |
| Ansys Mechanical | 2023 R2 | Ansys Inc. | 商用 |
| COMSOL | 6.1 | COMSOL AB | 商用 |
| Code_Aster | 15.4 | EDF | OSS (GPL) |
| CalculiX | 2.21 | — | OSS (GPL) |
試験問題一覧
総合比較結果
予算も時間も限られてるんですけど、コスパ最強はどれですか？
NAFEMS LE1（楕円膜、QUAD8、48x32メッシュ）
次は楕円膜の話ですね。どんな内容ですか？
| ソルバー | $\sigma_{yy}$ (MPa) | 誤差(%) | 計算時間(秒) | メモリ(MB) |
|---|---|---|---|---|
| 参照解 | 92.70 | — | — | — |
| MSC Nastran | 92.68 | 0.02 | 0.8 | 25 |
| NX Nastran | 92.69 | 0.01 | 0.7 | 24 |
| Abaqus | 92.71 | 0.01 | 1.2 | 32 |
| Ansys Mechanical | 92.68 | 0.02 | 1.0 | 28 |
| COMSOL | 92.72 | 0.02 | 1.5 | 45 |
| Code_Aster | 92.66 | 0.04 | 1.8 | 38 |
| CalculiX | 92.63 | 0.08 | 2.2 | 42 |
NAFEMS LE10（厚板、HEX20、20x10x6メッシュ）
メッシュって、具体的にはどういうことですか？
| ソルバー | $\sigma_{yy}$ (MPa) | 誤差(%) | 計算時間(秒) | メモリ(MB) |
|---|---|---|---|---|
| 参照解 | 5.38 | — | — | — |
| MSC Nastran | 5.377 | 0.06 | 2.3 | 85 |
| NX Nastran | 5.378 | 0.04 | 2.1 | 82 |
| Abaqus | 5.382 | 0.04 | 2.8 | 95 |
| Ansys Mechanical | 5.377 | 0.06 | 2.5 | 88 |
| COMSOL | 5.385 | 0.09 | 3.2 | 120 |
| Code_Aster | 5.374 | 0.11 | 3.5 | 105 |
| CalculiX | 5.371 | 0.17 | 4.0 | 110 |
NAFEMS LE11（熱応力、CAX8、中メッシュ）
「熱応力」について教えてください！
| ソルバー | $\sigma_{yy}$ (MPa) | 誤差(%) | 計算時間(秒) |
|---|---|---|---|
| 参照解 | -105.0 | — | — |
| MSC Nastran | -104.88 | 0.11 | 0.8 |
| NX Nastran | -104.90 | 0.10 | 0.7 |
| Abaqus | -105.05 | 0.05 | 1.0 |
| Ansys Mechanical | -104.92 | 0.08 | 0.9 |
| COMSOL | -105.08 | 0.08 | 1.2 |
| Code_Aster | -104.82 | 0.17 | 1.5 |
| CalculiX | -104.75 | 0.24 | 1.8 |
支配方程式
数式は苦手なんですけど…線形静解析ソルバー比較の式の「意味」を教えてもらえますか？
線形静解析の基本方程式:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
ここで $\mathbf{K}$ は全体剛性マトリクス、$\mathbf{u}$ は節点変位ベクトル、$\mathbf{F}$ は外力ベクトル。
応力の算出:
待って待って、線形静解析の基本方程ってことは、つまりこういうケースでも使えますか？
これを数式で表すとこうなるよ。
いやぁ、線形静解析ソルバー比較って奥が深いですね… でも先生の説明のおかげでだいぶ整理できました！
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
各項の物理的意味
- 保存量の時間変化項：対象とする物理量の時間的変化率を表す。定常問題では零となる。【イメージ】浴槽にお湯を張るとき、水位が時間と共に上がる——この「時間あたりの変化速度」が時間変化項。バルブを閉じて水位が一定になった状態が「定常」であり、時間変化項はゼロ。
- フラックス項（流束項）：物理量の空間的な輸送・拡散を記述する。対流と拡散の2種類に大別される。【イメージ】対流は「川の流れがボートを運ぶ」ように流れに乗って物が運ばれること。拡散は「インクが静止した水中で自然に広がる」ように濃度差で物が移動すること。この2つの輸送メカニズムの競合が多くの物理現象を支配する。
- ソース項（生成・消滅項）：物理量の局所的な生成または消滅を表す外力・反応項。【イメージ】部屋の中でヒーターをつけると、その場所に熱エネルギーが「生成」される。化学反応で燃料が消費されると質量が「消滅」する。外部から系に注入される物理量を表す項。
仮定条件と適用限界
- 連続体仮定が成立する空間スケールであること
- 材料・流体の構成則（応力-歪み関係、ニュートン流体則等）が適用範囲内であること
- 境界条件が物理的に妥当かつ数学的に適切に定義されていること
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 代表長さ $L$ | m | CADモデルの単位系と一致させること |
| 代表時間 $t$ | s | 過渡解析の時間刻みはCFL条件・物理的時定数を考慮 |
検証データの視覚化
理論値と計算値の比較を定量的に示す。誤差5%以内を合格基準とする。
| 評価項目 | 理論値/参照値 | 計算値 | 相対誤差 [%] | 判定 |
|---|---|---|---|---|
| 最大変位 | 1.000 | 0.998 | 0.20 | PASS |
| 最大応力 | 1.000 | 1.015 | 1.50 | PASS |
| 固有振動数（1次） | 1.000 | 0.997 | 0.30 | PASS |
| 反力合計 | 1.000 | 1.001 | 0.10 | PASS |
| エネルギー保存 | 1.000 | 0.999 | 0.10 | PASS |
判定基準: 相対誤差 < 1%: ■ 優良、1〜5%: ■ 許容、> 5%: ■ 要検討
V&V検証の効率化は、シミュレーションの信頼性を支える基盤です。 — Project NovaSolverは検証プロセスの改善にも注力しています。
線形静解析ソルバー比較の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →