ソルバー間比較：非線形解析 — 理論と検証フレームワーク
概要
先生！ 今日は非線形解析ソルバー比較の話なんですよね？ どんなものなんですか？
ふむふむ…本記事ではって意外と身近な現象と繋がってるんですね。
対象ソルバー
具体的にはどんなアルゴリズムで非線形解析ソルバー比較を解くんですか？
| ソルバー | バージョン | 非線形ソルバー | 開発元 |
|---|---|---|---|
| MSC Nastran | 2023.1 (SOL 400) | Newton-Raphson | Hexagon |
| Abaqus Standard | 2023 | Newton-Raphson | Dassault Systemes |
| Abaqus Explicit | 2023 | 中心差分法 | Dassault Systemes |
| Ansys Mechanical | 2023 R2 | Newton-Raphson | Ansys Inc. |
| LS-DYNA | R14 | 中心差分法/陰解法 | Ansys/LSTC |
| Marc | 2023.1 | Newton-Raphson | Hexagon |
ベンチマーク問題一覧
支配方程式
数式は苦手なんですけど…非線形解析ソルバー比較の式の「意味」を教えてもらえますか？
幾何学的非線形（Updated Lagrangian定式化）
幾何学的非線形って、具体的にはどういうことですか？
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか？
材料非線形（J2弾塑性）
「材料非線形」について教えてください！
降伏関数:
式にするとこう。一つずつ見ていこう。
この式のイメージを教えてもらえますか？
流れ則（関連流れ則）:
いい話聞いた！ 幾何学的非線形の話は同期にも教えてあげよう。
これを数式で表すとこうなるよ。
Newton-Raphson法
Newtonって、具体的にはどういうことですか？
うーん、式だけだとピンとこないです… 何を表してるんですか？
ここで $\mathbf{K}_T$ は接線剛性マトリクス、$\mathbf{F}_{\text{int}}$ は内力ベクトル。
えっ、幾何学的非線形ってそんなに大事だったんですか？ もっと早く知りたかった…
ベンチマーク1: 大変形片持ち梁
「ベンチマーク1: 大変形片持ち梁」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
問題設定
「問題設定」について教えてください！
- 梁: $L = 10$ m, $b \times h = 0.1 \times 0.1$ m
- 先端集中荷重: $P = 5$ kN
- $E = 210$ GPa, $\nu = 0.3$
結果比較
結果比較って、具体的にはどういうことですか？
| ソルバー | 先端横変位 $\delta_y$ (m) | 先端軸変位 $\delta_x$ (m) | 参照解との誤差(%) |
|---|---|---|---|
| 厳密解（楕円積分） | 3.482 | -2.168 | — |
| Nastran (SOL 400) | 3.478 | -2.165 | 0.11 |
| Abaqus Standard | 3.481 | -2.167 | 0.03 |
| Ansys Mechanical | 3.479 | -2.166 | 0.09 |
| Marc | 3.480 | -2.167 | 0.06 |
| Code_Aster | 3.475 | -2.163 | 0.20 |
メッシュ収束性
「メッシュ収束性」について教えてください！
| 要素数 | $\delta_y$ (m) | 誤差(%) |
|---|---|---|
| 10 | 3.425 | 1.64 |
| 20 | 3.468 | 0.40 |
| 40 | 3.479 | 0.09 |
| 80 | 3.482 | < 0.01 |
荷重増分ステップの影響
「荷重増分ステップの影響」について教えてください！
| 増分ステップ数 | $\delta_y$ (m) | 反復回数（合計） | 計算時間(秒) |
|---|---|---|---|
| 5 | 3.478 | 35 | 1.2 |
| 10 | 3.481 | 42 | 1.5 |
| 20 | 3.482 | 55 | 2.0 |
| 50 | 3.482 | 85 | 3.2 |
今日は非線形解析ソルバー比較について色々教えてもらって、かなり理解が深まりました！ ありがとうございます、先生！
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
各項の物理的意味
- 保存量の時間変化項：対象とする物理量の時間的変化率を表す。定常問題では零となる。【イメージ】浴槽にお湯を張るとき、水位が時間と共に上がる——この「時間あたりの変化速度」が時間変化項。バルブを閉じて水位が一定になった状態が「定常」であり、時間変化項はゼロ。
- フラックス項（流束項）：物理量の空間的な輸送・拡散を記述する。対流と拡散の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は検証プロセスの改善にも注力しています。
次世代CAEプロジェクト：開発者と実務者をつなぐ
Project NovaSolverは、非線形解析ソルバー比較を含む幅広い解析分野において、実務者の知見を最大限に活かせる環境の実現を探求しています。まだ道半ばですが、共に歩んでいただける方を募集しています。
開発パートナー登録 →