NAFEMS LE1 楕円膜の内圧問題 — V&V結果総括レポート

いい質問だね。LE1は楕円膜に内圧をかけたときの応力を求める問題で、評価点Dにおける $\sigma_{yy}$ の参照解が 92.7 MPa だ。きみの結果の誤差は $(92.7 - 92.5)/92.7 \approx 0.22\%$ だから、十分に高精度だよ。でも大事なのは単一の結果だけじゃなくて、メッシュを細かくしたときの収束挙動を確認することなんだ。

そのとおり。メッシュを細かくするたびに解が参照解に近づいていく「単調収束」が確認できれば、解の信頼性が高いと言える。逆に細かくしても値が振動するようだと、モデルに何か問題がある可能性がある。今回は7つのソルバーで検証したから、その結果を一緒に見てみよう。

LE1は、楕円形の穴が開いた板（平面応力状態）に一様な内圧を加える問題だ。楕円の長軸・短軸比が応力集中に影響して、評価点での $\sigma_{yy}$ が解析解で正確に求まるから、FEAソルバーの精度検証に最適なんだ。例えるなら、自動車のボディに楕円形の開口部（窓やドア穴）があるとき、その周辺の応力集中を正しく計算できるかのテストだね。

十分に細かいメッシュ（48×32分割、二次要素）で計算した結果がこれだ。参照解は $\sigma_{yy} = 92.70$ MPa だよ。

そう、LE1は線形弾性の平面応力問題だから、ソルバーの実装が正しければ理論的に同じ解に収束する。ソルバー間のばらつきは最大でも0.09 MPa（0.1%未満）で、これは数値積分の精度や要素定式化のわずかな違いに起因するものだ。実務では「どのソルバーを使うか」よりも「メッシュをどう切るか」の方がはるかに重要ってことが、このデータからよく分かるだろう？

いい着眼点だ。要素タイプごとに3段階のメッシュ密度で比較した結果がこれだよ。

驚くよね。TRIA3（一次三角形要素）は「定ひずみ要素」とも呼ばれていて、要素内のひずみが一定値しか表現できない。だから応力集中のような勾配の大きい場で正確な値を出すには、膨大な数の要素が必要になる。一方、QUAD8は二次の変位関数を持っているから、応力の空間的な変化を滑らかに表現できるんだ。例えば自動車のフィレット部の応力解析で、TRIA3しか使えない状況だと結果を信用できないレベルの誤差が出ることがある。

有限要素法の誤差は、要素サイズ $h$ に対して次のように減少する：

GCI（Grid Convergence Index）は、Roacheが提案したメッシュ収束性の定量的指標だ。3段階のメッシュで得た結果からRichardson外挿を行い、「真値」に対する離散化誤差の上限を推定する。今回のQUAD8の結果ではGCI = 0.27%で、これは「メッシュをこれ以上細かくしても結果は0.27%以上は変わらない」という意味だよ。実務のV&V報告書ではGCI < 1%が一般的な合格基準として使われる。

同じ精度レベルを達成するのに必要な計算コストを比較してみよう。

そう、これがまさにポイントだ。「要素数を増やして力技で精度を上げる」よりも、「二次要素で賢くメッシュを切る」方が圧倒的にコストパフォーマンスがいい。QUAD8の粗96要素で達成できる精度を、QUAD4で実現しようとすると1,536要素以上が必要で、DOFは4倍以上になる。大規模モデルになるほど、この差は計算時間・メモリの両面で効いてくる。

ASME V&V 10やNAFEMSガイドラインに基づいて、4つの判定項目を設定した。結果はこのとおりだ。

まとめると、LE1は線形弾性の2D平面応力問題として、全ソルバーが参照解を高精度に再現できることが確認された。実務への教訓として覚えておいてほしいのは3点だ：

そう、まずは「二次要素でメッシュを切って、3段階で収束確認」を習慣にしよう。LE1のような簡単な問題で練習しておくと、複雑なモデルでもスムーズにV&Vができるようになるよ。