NAFEMS FE 固有値解析ベンチマーク — 理論解と検証
概要
先生！ 今日はNAFEMS FE 固有値の話なんですよね？ どんなものなんですか？
本記事ではNAFEMSが定める代表的な固有値ベンチマーク問題（FV32: 厚肉円筒、FV52: 片持ち梁など）について詳しく見ていこう。
ベンチマークシリーズの具体的な数値例とかあると、もっとピンとくるんですけど…
ベンチマーク問題一覧
「ベンチマーク問題一覧」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
FV32: 自由振動する厚肉円筒
次は自由振動する厚肉円筒の話ですね。どんな内容ですか？
#### 幾何形状
- 内半径: $r_i = 1.0$ m
- 外半径: $r_o = 2.0$ m
- 軸方向長さ: $L = 2.0$ m
- 完全自由境界
おお〜、自由振動する厚肉円筒の話、めちゃくちゃ面白いです！ もっと聞かせてください。
FV52: 片持ち梁の自由振動
次は片持ち梁の自由振動の話ですね。どんな内容ですか？
#### 幾何形状
- 長さ: $L = 10.0$ m
- 断面: $W \times H = 1.0 \times 2.0$ m
- 片持ち支持（$x = 0$ を固定）
先輩が「自由振動する厚肉円筒だけはちゃんとやれ」って言ってた意味が分かりました。
支配方程式
いよいよ数式ですね…！ NAFEMS FE 固有値ではどんな方程式が出てくるんですか？
自由振動の固有値問題:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
ここで $\mathbf{K}$ は剛性マトリクス、$\mathbf{M}$ は質量マトリクス、$\omega$ は固有角振動数、$\boldsymbol{\phi}$ は固有モードベクトル。
先生の説明分かりやすい！ 自由振動の固有値問題のモヤモヤが晴れました。
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか？
片持ち梁のEuler-Bernoulli理論解:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
ここで $\lambda_1 = 1.8751$, $\lambda_2 = 4.6941$, $\lambda_3 = 7.8548$ ...
なるほど…自由振動の固有値問題って一見シンプルだけど、実はすごく奥が深いんですね。
NAFEMS FV52 参照解と数値解の比較
予算も時間も限られてるんですけど、コスパ最強はどれですか？
参照固有振動数
参照固有振動数って、具体的にはどういうことですか？
| モード | 理論解 (Hz) | Ansys Mechanical | Abaqus | Nastran | 相対誤差(%) |
|---|---|---|---|---|---|
| 1次曲げ（面内） | 44.623 | 44.620 | 44.625 | 44.618 | < 0.02 |
| 2次曲げ（面内） | 109.44 | 109.42 | 109.45 | 109.41 | < 0.03 |
| 1次曲げ（面外） | 17.849 | 17.847 | 17.850 | 17.846 | < 0.02 |
| 2次曲げ（面外） | 43.78 | 43.77 | 43.78 | 43.77 | < 0.03 |
| 1次ねじり | 27.35 | 27.34 | 27.36 | 27.34 | < 0.04 |
| 1次軸振動 | 125.0 | 124.98 | 125.01 | 124.97 | < 0.03 |
メッシュ収束性（FV52 第1次曲げモード）
メッシュの良し悪しってどうやって判断するんですか？
| メッシュ密度 | 要素数 | 自由度数 | $f_1$ (Hz) | 誤差(%) |
|---|---|---|---|---|
| 粗い | 50 | 900 | 45.85 | 2.75 |
| 中程度 | 200 | 3,300 | 44.82 | 0.44 |
| 細かい | 800 | 12,600 | 44.66 | 0.08 |
| 非常に細かい | 3,200 | 49,200 | 44.63 | 0.02 |
要素タイプ別の比較（FV52 第1次モード）
結局どれを選べばいいか、判断基準を教えてもらえますか？
| 要素タイプ | 節点数 | $f_1$ (Hz) | 誤差(%) |
|---|---|---|---|
| TET4 | 12,600 | 47.52 | 6.49 |
| TET10 | 12,600 | 44.72 | 0.22 |
| HEX8 | 12,600 | 44.95 | 0.73 |
| HEX20 | 12,600 | 44.63 | 0.02 |
| BEAM2 | 300 | 44.62 | < 0.01 |
梁要素は理論解に最も近い結果を返す（Euler-Bernoulli理論と直接対応するため）。ソリッド要素ではHEX20が最も高精度。
今日はNAFEMS FE 固有値について色々教えてもらって、かなり理解が深まりました！ ありがとうございます、先生！
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
各項の物理的意味
- 保存量の時間変化項：対象とする物理量の時間的変化率を表す。定常問題では零となる。【イメージ】浴槽にお湯を張るとき、水位が時間と共に上がる——この「時間あたりの変化速度」が時間変化項。バルブを閉じて水位が一定になった状態が「定常」であり、時間変化項はゼロ。
- フラックス項（流束項）：物理量の空間的な輸送・拡散を記述する。対流と拡散の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は、NAFEMS FE 固有値を含む幅広い解析分野において、実務者の知見を最大限に活かせる環境の実現を探求しています。まだ道半ばですが、共に歩んでいただける方を募集しています。
開発パートナー登録 →