NAFEMS LE11 厚肉シェルベンチマーク — V&V結果総括
問題定義
先生、NAFEMS LE11ってどんなベンチマーク問題ですか？
LE11は厚肉の複合形状体（円柱+テーパー+球）に温度荷重を与える問題だ。形状が解析的に記述可能で、温度場も既知の関数で与えられるため、正確な参照解が存在する。ポイントは、厚肉構造をシェル要素でモデル化した場合の精度を評価することにある。薄肉シェル要素では板厚方向の応力分布を捉えきれないため、厚肉シェル要素（3Dソリッド要素）との比較が重要なんだ。
つまり「シェル要素を使っていいかどうかの判断基準」を提供するベンチマークですね？
そういうことだ。実務でよくある「この構造、シェルでいける？ ソリッドにすべき？」という判断にエビデンスを与えてくれる。温度荷重の場合は板厚方向の温度勾配が重要だから、シェル要素の限界が特に顕著に出る。
要素タイプ別の精度比較
どの要素タイプが一番精度が良いんですか？
検証結果を見ればはっきりする。細メッシュでの参照解との誤差を見てくれ：
| 要素タイプ | 粗メッシュ誤差(%) | 中メッシュ誤差(%) | 細メッシュ誤差(%) | 推奨 |
|---|---|---|---|---|
| TET4（1次四面体） | 22.5 | 12.1 | 5.8 | 非推奨 |
| TET10（2次四面体） | 3.2 | 0.76 | 0.18 | 細メッシュ以上 |
| HEX8（1次六面体） | 8.5 | 2.8 | 1.14 | 非常に細かい以上 |
| HEX20（2次六面体） | 1.2 | 0.25 | 0.08 | 中メッシュ以上 |
| CAX4（1次軸対称） | 6.5 | 1.43 | 0.35 | 細メッシュ以上 |
| CAX8（2次軸対称） | 0.45 | 0.11 | 0.03 | 中メッシュ以上 |
TET4が圧倒的にダメですね...粗メッシュで22.5%も誤差があるとは。
TET4は線形の形状関数しか持たないから、温度勾配による曲げ変形を表現できない。要素内で応力が一定（定ひずみ要素）だから、本来は滑らかに変化する応力場をカクカクの階段状にしか近似できないんだ。この問題のように温度勾配が重要なケースでは特に悲惨な結果になる。温度荷重問題にTET4は使ってはいけない、これは鉄則だ。
メッシュ収束性の評価
2次要素なら粗いメッシュでもいいですか？
2次要素は収束が速いけど、粗すぎるとまだ誤差は残る。LE11の結果からわかる重要な法則がある：2次要素は収束次数が約2（誤差がメッシュサイズの2乗に比例して減少）、一方1次要素は収束次数が約1。つまり2次要素はメッシュを半分にすると誤差が1/4になるけど、1次要素は半分にしかならない。だから同じ精度を得るのに必要な要素数が桁違いに少なくて済む。
GCI（Grid Convergence Index）評価
Richardson外挿に基づくGCI評価で、離散化誤差の上界を推定した：
| パラメータ | HEX20 | TET10 | CAX8 |
|---|---|---|---|
| 観測収束次数 $p$ | 2.05 | 1.88 | 2.12 |
| $\text{GCI}_{\text{fine}}$ (%) | 0.18 | 0.52 | 0.08 |
GCIが全部1%以下ということは、結果を信頼していいレベルですね？
その通り。ASME V&V 20-2009の基準では、GCI < 1%であれば「十分に収束している」と判定できる。このベンチマークでは2次要素を使えば中メッシュでも実用上問題ない精度が得られることが実証された。
計算コスト比較
| モデルタイプ | 要素数 | DOF | 計算時間(秒) | メモリ(MB) | 誤差(%) |
|---|---|---|---|---|---|
| 3D HEX20（細） | 3,200 | 52,800 | 15.2 | 450 | 0.08 |
| 3D TET10（細） | 8,500 | 52,800 | 18.5 | 520 | 0.76 |
| 軸対称 CAX8（細） | 800 | 7,200 | 0.8 | 15 | 0.03 |
軸対称モデルの効率がすごい...計算時間19分の1で精度は最高ですか？
軸対称性を活用できる形状なら、絶対にCAX8を使うべきだ。3Dソリッドの20分の1の計算コストで、より高精度な結果が得られる。これは大規模モデルのデバッグ時にも有効で、まず軸対称でパラメータを確認してから3Dに展開する、という戦略は実務で非常に有効だよ。
結論と推奨事項
| 判定項目 | 基準 | 結果 | 判定 |
|---|---|---|---|
| メッシュ収束性 | 単調収束（3水準以上） | 4水準で確認 | 合格 |
| 参照解との一致 | 誤差 < 2% | 全ソルバー < 0.25% | 合格 |
| GCI | < 1% | 0.08〜0.52% | 合格 |
| 3D/軸対称の整合性 | 差 < 0.5% | 0.04% | 合格 |
LE11ベンチマークから得られる実務上の教訓をまとめると：
- 軸対称モデルが使える場合は、CAX8要素が圧倒的に効率的
- 3Dモデルが必要な場合は、HEX20要素を推奨
- TET4要素は温度荷重問題に不適切（精度不足）
- 温度場の付与方法として、解析関数またはユーザーサブルーチンを推奨
- 厚肉構造にシェル要素を使う場合は、板厚方向の積分点数を十分に取る（5点以上推奨）
すごく実践的な知見ですね。自分の解析でも要素タイプの選定をもっと慎重にやります！
検証データの視覚化
理論値と計算値の比較を定量的に示す。誤差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%: ■ 要検討
なった
詳しく
報告