NAFEMS LE1: 楕円膜の平面応力 — 理論と支配方程式
概要
先生！ 今日はNAFEMS LE1: 楕円膜の平面応力の話なんですよね？ どんなものなんですか？
NAFEMS線形弾性ベンチマークLE1。楕円膜に一様内圧を負荷し、内縁の点Dにおける法線応力を評価する標準問題。
参照解: $$ \sigma_{yy}(D) = 92.7 \text{ MPa} $$
先輩が「線形弾性ベンチマークだけはちゃんとやれ」って言ってた意味が分かりました。
問題設定
「問題設定」について教えてください！
- 形状: 楕円膜（長軸 2a = 4m, 短軸 2b = 2m）
- 材料: 等方弾性体（E = 210 GPa, ν = 0.3）
- 荷重: 一様内圧 p = 10 MPa
- 拘束: 対称条件（1/4モデル）
- 評価点: 点D（短軸上の内縁）
ここまで聞いて、線形弾性ベンチマークがなぜ重要か、やっと腹落ちしました！
支配方程式
いよいよ数式ですね…！ NAFEMS LE1: 楕円膜の平面応力ではどんな方程式が出てくるんですか？
平面応力条件下の弾性体の釣り合い方程式:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
構成方程式（平面応力）:
式にするとこう。一つずつ見ていこう。
先輩が「平面応力条件下の弾性だけはちゃんとやれ」って言ってた意味が分かりました。
理論解と数値解の比較
予算も時間も限られてるんですけど、コスパ最強はどれですか？
各ソルバーによるベンチマーク検証データ
各ソルバーによるベンチマーク検って、具体的にはどういうことですか？
| 評価項目 | 参照解 | Ansys Mechanical | Abaqus | MSC Nastran | COMSOL | 単位 |
|---|---|---|---|---|---|---|
| σ_yy (点D) | 92.7 | 92.68 | 92.71 | 92.65 | 92.72 | MPa |
| σ_xx (点D) | -10.0 | -10.01 | -9.99 | -10.02 | -10.00 | MPa |
| 最大主応力 | 92.7 | 92.69 | 92.70 | 92.66 | 92.71 | MPa |
メッシュ収束性検証
次はメッシュ収束性検証の話ですね。どんな内容ですか？
| メッシュ密度 | 要素数 | 自由度数(DOF) | σ_yy (MPa) | 相対誤差(%) |
|---|---|---|---|---|
| 非常に粗い | 24 | 168 | 85.3 | 7.98 |
| 粗い | 96 | 624 | 89.5 | 3.45 |
| 中程度 | 384 | 2,400 | 91.8 | 0.97 |
| 細かい | 1,536 | 9,408 | 92.5 | 0.22 |
| 非常に細かい | 6,144 | 37,248 | 92.68 | 0.02 |
要素タイプ別比較（中程度メッシュ）
「要素タイプ別比較」について教えてください！
| 要素タイプ | 要素名 | 節点数 | σ_yy (MPa) | 相対誤差(%) |
|---|---|---|---|---|
| QUAD4（4節点四角形） | CPS4 / PLANE182 | 384 | 88.2 | 4.85 |
| QUAD8（8節点四角形） | CPS8 / PLANE183 | 384 | 92.5 | 0.22 |
| TRIA3（3節点三角形） | CPS3 / PLANE182 | 768 | 82.1 | 11.4 |
| TRIA6（6節点三角形） | CPS6 / PLANE183 | 768 | 92.3 | 0.43 |
収束特性
次は収束特性の話ですね。どんな内容ですか？
- QUAD8（二次要素）: $O(h^4)$ の超収束を示す
- QUAD4（一次要素）: $O(h^2)$ の収束速度
- TRIA3（一次三角形）: 精度が低く、応力が過小評価される傾向
各ソルバーによるベンの具体的な数値例とかあると、もっとピンとくるんですけど…
離散化手法
連続的な式をバラバラにして解くって聞いたんですけど、具体的にはどうするんですか？
有限要素法（FEM）による離散化。弱形式への変換:
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
形状関数 $N_i$ を用いて変位場を近似:
式にするとこう。一つずつ見ていこう。
この式のイメージを教えてもらえますか？
要素剛性マトリクスの構築:
へぇ〜！ 有限要素法についてだいぶ理解が深まりました。メモメモ…📝
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか？
ここで $[B]$ はひずみ-変位マトリクス、$[D]$ は弾性マトリクス、$t$ は板厚。
へぇ〜！ 有限要素法についてだいぶ理解が深まりました。メモメモ…📝
行列解法アルゴリズム
行列解法アルゴリズムって、具体的にはどういうことですか？
| 解法 | 分類 | メモリ使用量 | 適用規模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小〜中規模 |
| Cholesky分解 | 直接法（対称正定値） | O(n²) | 小〜中規模 |
| PCG法 | 反復法 | O(n) | 大規模 |
| GMRES法 | 反復法 | O(n·m) | 大規模・非対称 |
| AMG前処理 | 前処理 | O(n) | 超大規模 |
あっ、そういうことか！ 有限要素法ってそういう仕組みだったんですね。
商用ツールにおける実装
いろんなソフトがあるんですよね？ それぞれの特徴を教えてください！
| ツール名 | 開発元/現在 | 主要ファイル形式 |
|---|---|---|
| MSC Nastran / NX Nastran | MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software） | .bdf, .dat, .f06, .op2, .pch |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| COMSOL Multiphysics | COMSOL AB | .mph |
ベンダーの系譜と製品統合の経緯
各ソフトの成り立ちって、結構ドラマチックだったりしますか？
MSC Nastran / NX Nastran
次はMSC Nastranの話ですね。どんな内容ですか？
NASA構造解析（NASTRAN）として1960年代にNASAが開発。MSC Softwareが商用化し、後にUGS（現Siemens）がNX Nastranを派生。MSCは2017年にHexagon ABに買収。
現在の所属: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）
Abaqus FEA (SIMULIA)
Abaqus FEAって、具体的にはどういうことですか？
1978年にHKS (Hibbitt, Karlsson & Sorensen) が開発。非線形解析に強み。2005年にDassault Systèmesが買収し、SIMULIAブランドに統合。
現在の所属: Dassault Systèmes SIMULIA
待って待って、構造解析ってことは、つまりこういうケースでも使えますか？
Ansys Mechanical (旧ANSYS Structural)
「Ansys Mechanical」について教えてください！
1970年にSwanson Analysis Systems Inc. (SASI) が開発。APDL（Ansys Parametric Design Language）ベース。
現在の所属: Ansys Inc.
ふむふむ…構造解析って意外と身近な現象と繋がってるんですね。
ファイル形式と相互運用性
異なるソルバー間での要素タイプの変換リスク:
- QUAD8要素の中間節点配置が異なる場合がある
- 積分スキーム（完全積分/低減積分）の既定値がソルバーにより異なる
- 材料モデルの入力パラメータの定義差異に注意
先生の説明分かりやすい！ フォーマットのモヤモヤが晴れました。
判定基準
「判定基準」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
NAFEMS推奨の合格基準: 参照解に対して相対誤差1%以内（適切なメッシュ密度において）
今日は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は検証プロセスの改善にも注力しています。
CAEの未来を、実務者と共に考える
Project NovaSolverは、NAFEMS LE1: 楕円膜の平面応力における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →