NAFEMS LE1 楕円膜の内圧問題 — 理論解と検証
概要
先生！ 今日はNAFEMS LE1の話なんですよね？ どんなものなんですか？
NAFEMS LE1ベンチマークは、楕円形の膜構造に一様内圧を負荷した際の平面応力問題なんだ。本問題はNAFEMS（National Agency for Finite Element Methods and Standards）が定めた線形弾性ベンチマーク試験の一つであり、有限要素解析コードの精度検証に広く利用されている。
本問題の参照解は、内縁最上部（点D）における法線応力 $\sigma_{yy} = 92.7$ MPa であり、これを数値解の妥当性評価基準とする。
ここまで聞いて、ベンチマークはがなぜ重要か、やっと腹落ちしました！
問題設定
先生、「問題設定」について教えてください！
幾何形状
幾何形状って、具体的にはどういうことですか？
楕円膜の形状パラメータ:
- 外側楕円: 長軸 $a = 3.25$ m、短軸 $b = 2.75$ m
- 内側楕円: 長軸 $a = 2.0$ m、短軸 $b = 1.0$ m
- 板厚: $t = 0.1$ m（平面応力条件）
対称性を利用して1/4モデルを使用可能。
待って待って、幾何形状ってことは、つまりこういうケースでも使えますか？
材料特性
等方性線形弾性体を仮定する。
なるほど！ 幾何形状のイメージがつかめてきました！
境界条件
次は境界条件の話ですね。どんな内容ですか？
- 外縁: 完全固定（$u_x = u_y = 0$）
- 内縁: 一様内圧 $p = 10$ MPa を法線方向に負荷
- 対称面: $x = 0$ 面で $u_x = 0$、$y = 0$ 面で $u_y = 0$
なるほど…幾何形状って一見シンプルだけど、実はすごく奥が深いんですね。
支配方程式
いよいよ数式ですね…！ NAFEMS LE1ではどんな方程式が出てくるんですか？
平面応力条件下の2次元弾性問題の支配方程式:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
式にするとこう。一つずつ見ていこう。
この式のイメージを教えてもらえますか？
ひずみ-変位関係:
先生の説明分かりやすい！ 平面応力条件下ののモヤモヤが晴れました。
これを数式で表すとこうなるよ。
いい話聞いた！ 平面応力条件下のの話は同期にも教えてあげよう。
解析解（参照値）
先生、「解析解（参照値）」について教えてください！
点D（内縁最上部、$x = 0$, $y = 1.0$）における法線応力:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
この値はNAFEMSが公式に認定した参照解であり、高次要素による収束解析の外挿値に基づく。
なるほど。じゃあ内縁最上部ができていれば、まずは大丈夫ってことですか？
理論解と数値解の比較
結局どれを選べばいいか、判断基準を教えてもらえますか？
| 評価項目 | 理論解（参照値） | Ansys Mechanical | Abaqus | Nastran | 相対誤差(%) |
|---|---|---|---|---|---|
| 点Dにおける $\sigma_{yy}$ (MPa) | 92.7 | 92.68 | 92.71 | 92.65 | < 0.1 |
| 点Dにおける $\sigma_{xx}$ (MPa) | -10.0 | -10.02 | -9.98 | -10.01 | < 0.2 |
| 最大変位 $u_{\text{max}}$ (mm) | 0.1139 | 0.1138 | 0.1140 | 0.1137 | < 0.2 |
| 内縁最大von Mises応力 (MPa) | 97.2 | 97.18 | 97.22 | 97.15 | < 0.1 |
上記は各ソルバーで十分に細かいメッシュ（QUAD8/二次四角形要素、48×32分割以上）を使用した場合の結果なんだ。
評価項目の具体的な数値例とかあると、もっとピンとくるんですけど…
メッシュ収束性
メッシュの良し悪しってどうやって判断するんですか？
| メッシュ密度 | 要素数 | 自由度数 | $\sigma_{yy}$ (MPa) | 誤差(%) |
|---|---|---|---|---|
| 粗い | 96 | 594 | 88.2 | 4.85 |
| 中程度 | 384 | 2,178 | 91.0 | 1.83 |
| 細かい | 1,536 | 8,450 | 92.1 | 0.65 |
| 非常に細かい | 6,144 | 33,282 | 92.65 | 0.05 |
上記はQUAD4（線形四角形）要素を使用した場合の収束挙動を示す。
ふむふむ…メッシュ密度って意外と身近な現象と繋がってるんですね。
要素タイプ別の比較
結局どれを選べばいいか、判断基準を教えてもらえますか？
| 要素タイプ | 節点数 | $\sigma_{yy}$ (MPa) | 誤差(%) |
|---|---|---|---|
| QUAD4（4節点四角形） | 8,450 | 92.1 | 0.65 |
| QUAD8（8節点四角形） | 8,450 | 92.6 | 0.11 |
| TRIA3（3節点三角形） | 8,450 | 86.5 | 6.69 |
| TRIA6（6節点三角形） | 8,450 | 92.5 | 0.22 |
二次要素（QUAD8, TRIA6）は線形要素に比べて格段に高い精度を示す。特にTRIA3は収束が遅く、本問題への適用は推奨されない。
今日はNAFEMS LE1について色々教えてもらって、かなり理解が深まりました！ ありがとうございます、先生！
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
各項の物理的意味
- 保存量の時間変化項：対象とする物理量の時間的変化率を表す。定常問題では零となる。【イメージ】浴槽にお湯を張るとき、水位が時間と共に上がる——この「時間あたりの変化速度」が時間変化項。バルブを閉じて水位が一定になった状態が「定常」であり、時間変化項はゼロ。
- フラックス項（流束項）：物理量の空間的な輸送・拡散を記述する。対流と拡散の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は検証プロセスの改善にも注力しています。
NAFEMS LE1の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →