20節点六面体要素(HEX20)
20節点六面体要素(HEX20)の理論基礎
HEX20 — 最高精度の3次元要素
HEX20はFEMの3次元要素で最も精度が高いんですか？
実用的な要素としては最高精度と言ってよい。20節点の二次六面体要素で、各節点に3自由度、合計60自由度。HEX8（24自由度）の2.5倍のDOFだが、精度は桁違いに高い。
形状関数
HEX20は8つの頂点ノードと12個の辺中点ノードを持つ。形状関数は三次曲面（serendipity）型：
頂点ノード：
辺中点ノード（例: $\xi_i = 0$）：
Serendipity型って何ですか？ Lagrange型とは違うんですか？
Lagrange型（27節点HEX）は面と内部にも節点を持つが、Serendipity型（20節点HEX）は辺上のみに中間節点を持ち、面の中心や内部に節点がない。Serendipity型はDOFが少ないのに精度はLagrange型とほぼ同等で、実用上はSerendipity型が圧倒的に多い。
HEX20の精度
HEX20はどの程度正確ですか？
TET10との比較で：
| 同じDOF数での精度 | HEX20 | TET10 |
|---|---|---|
| 変位 | 非常に高い | 高い |
| 応力 | 非常に高い | 高い |
| 効率（DOFあたり精度） | 最高 | HEX20の50〜70% |
収束速度：
- 変位: $O(h^4)$（HEX20）vs. $O(h^3)$（TET10）
- 応力: $O(h^3)$（HEX20）vs. $O(h^2)$（TET10）
HEX20のほうが1オーダー速く収束する。メッシュが半分の粗さでも同等精度…。
だからHEXメッシュが作れる形状では、HEX20が最も効率的だ。問題はメッシュ生成の手間。この「精度効率 vs. メッシュ生成コスト」のトレードオフが、HEX20とTET10の選択の本質だ。
積分スキーム
HEX20の積分：
| 積分 | Gauss点 | 特徴 |
|---|---|---|
| 完全積分（3×3×3） | 27 | 最高精度。ロッキングのリスクは低い |
| 低減積分（2×2×2） | 8 | シアロッキング回避。アワーグラスは1モード |
HEX20の低減積分は2×2×2 = 8点ですか。HEX8の完全積分と同じ点数。
そう。HEX20の低減積分（C3D20R）はHEX8の完全積分と同じ8つの積分点で、はるかに高い精度を出す。C3D20Rは3次元FEMで最も効率的な要素の一つと評される。
C3D20Rにアワーグラスの問題はありますか？
HEX20の低減積分にはアワーグラスモードが1つだけ存在する。HEX8Rの12個と比べて極めて少なく、実用上ほとんど問題にならない。ただし1要素のパッチテストでは注意が必要。
いつHEX20を使うか
HEX20を使うべき場面は？
逆にHEX20が不利な場面は？
まとめ
HEX20の理論を整理します。
要点：
- 20節点、Serendipity型、二次六面体 — 実用3D要素の最高精度
- C3D20R（低減積分）が最推奨 — 8積分点で最高効率
- 収束が速い — TET10より1オーダー速い
- メッシュ生成が課題 — 自動HEXは困難。Sweep/マップドメッシュ向き
- 大変形・陽解法には不向き — HEX8RやTET10のほうが安定
HEX20は「作れるなら最高の要素」ということですね。
まさにそう。「最高の精度」と「作れるならば」の間にメッシュ生成の壁がある。この壁を越えられるエンジニアにとって、HEX20は最強の武器だ。
二次六面体要素の理論的優位性
20節点六面体要素（Serendipity要素）は1966年にErgatoudisらが提案したアイソパラメトリック定式化の代表作だ。完全2次多項式を含む形状関数により、一次要素の8倍近い要素内での応力精度を持つ。1970年代のNASAアポロ計画以降の宇宙機構造解析で採用され、「金の標準要素」と呼ばれた時代がある。
20節点六面体要素(HEX20)の数値計算手法
C3D20R — 最効率の3D要素
C3D20R（AbaqusのHEX20低減積分）が「最も効率的」とされる理由を教えてください。
3つの理由がある：
1. 精度/DOF比が最高 — 同じDOF数でTET10の2〜3倍の精度
2. シアロッキングがない — 低減積分で回避済み
3. 体積ロッキングも軽微 — $\nu = 0.49$ 程度までは問題なし
C3D20（完全積分）ではなくC3D20R（低減積分）が推奨なんですね。
C3D20は27積分点で計算コストが高い上に、$\nu$ が大きい材料で体積ロッキングの兆候が出ることがある。C3D20Rは8積分点で計算が速く、ロッキングも少ない。実務ではC3D20Rが圧倒的に多い。
ソルバー別の要素名
| バリエーション | Abaqus | Nastran | Ansys |
|---|---|---|---|
| 完全積分 | C3D20 | CHEXA(20) | SOLID186(full) |
| 低減積分 | C3D20R | — | SOLID186(red.) |
| ハイブリッド | C3D20H, C3D20RH | — | u-P対応 |
NastranのCHEXA(20)は低減積分オプションがないんですか？
NastranのCHEXA(20)はデフォルトで2×2×2（低減）積分を使う。完全積分に切り替えるオプションはあるが、デフォルトが低減。つまりNastranのCHEXA(20) ≈ AbaqusのC3D20R。
中間節点の扱い
HEX20の中間節点は辺の中点にありますが、曲面ではCADにスナップさせるんですよね。
TET10と同じだ。曲面上の中間節点をCAD面にスナップさせることで、曲面の二次近似が可能になる。ただしHEX20特有の注意点がある：
中間節点が辺の中点から大きくずれると要素が退化する。曲率が大きい部分で中間節点をCADにスナップさせると、中間節点が「辺の中点」から離れすぎてヤコビアンが負になることがある。
どうすればいいですか？
対策：
- 曲率が大きい部分ではメッシュを細かくして、辺の長さに対する中間節点のずれを小さくする
- ヤコビアンチェックを必ず実施
- 中間節点のスナップ量に上限を設けるプリプロセッサの設定を利用
HEX20の適用限界
HEX20が苦手な問題はありますか？
大変形解析ではHEX20は不利だ。中間節点が物理的なアンカーなしに動くため、大変形で要素が歪んだときに中間節点が「おかしな位置」に飛ぶことがある。
金属成形のような大変形問題では：
HEX20は主に線形〜中程度の非線形（接触、小変形塑性）に使う。
まとめ
HEX20の数値手法、整理します。
要点：
- C3D20R（低減積分）が最推奨 — 8積分点で最高効率
- NastranのCHEXA(20)はデフォルトで低減積分 — C3D20R相当
- 中間節点のスナップ量に注意 — 曲率が大きい部分でヤコビアン負のリスク
- 大変形にはHEX8を使う — HEX20は線形〜中非線形向き
- 完全積分C3D20は特殊用途のみ — 通常はC3D20Rで十分
27点vs14点ガウス積分の比較
HEX20要素の完全積分には3×3×3=27ガウス点が必要だが、計算コストが高い。1974年にBarlow点理論により14点（reduced integration）でも精度を保てることが示され、ANSYS SOLID186では3×3×3がデフォルトだが、Abaqus C3D20Rは2×2×2（8点）の縮減積分を採用し計算時間を約40%短縮しつつ精度を維持する。
20節点六面体要素(HEX20)の実務適用
HEX20の実務適用
HEX20を実務で使う場面を教えてください。
HEX20が特に威力を発揮する場面：
圧力容器の応力分類
ASME Div. 2の応力分類には板厚方向の応力分布が重要。HEX20は板厚方向2要素で正確な応力勾配を捕捉できる（C3D20Rの場合）。TET10だと4要素以上必要。
ベンチマーク検証
FEM解析の精度検証（NAFEMSベンチマーク等）ではHEX20が基準要素として使われることが多い。「最も正確な結果」を出す参照解としてHEX20で解いておき、TET10やシェルの結果と比較する。
接触解析の精密評価
ボルト締結部やベアリングの接触圧分布を精密に評価する場合、HEX20の接触面はTET10より安定で滑らか。接触面のメッシュが粗くても比較的正確な圧力分布が得られる。
メッシュ生成戦略
HEX20メッシュを効率よく作る方法は？
| 手法 | 適用形状 | 効率 |
|---|---|---|
| Sweepメッシュ | 押し出し体（パイプ、シャフト） | 非常に高い |
| 回転メッシュ | 回転体（フランジ、ディスク） | 非常に高い |
| マッピングメッシュ | 直方体的領域 | 高い |
| Multi-zone | やや複雑な形状を分割 | 中程度 |
| 手動（HyperMesh等） | 任意形状 | 低い（手間大） |
回転体ならHEX20が最も効率的ですね。
圧力容器やタービンディスクのような回転体は、2次元断面を周方向に回転してHEX20を生成できる。このケースではHEX20の精度効率が最大限に発揮される。
TET10からHEX20への移行判断
TET10で解析した結果をHEX20で検証したい場合、どうすればいいですか？
着目部位だけHEX20でサブモデルを作る手法が実務的だ：
1. TET10の全体モデルで解析
2. 着目部位の境界変位を抽出
3. HEX20のサブモデルに境界条件として適用
4. HEX20の結果とTET10を比較
全体をHEX20にする必要はないんですね。
全体をHEX20にするのはメッシュ生成コストが膨大。サブモデリングで「必要な部分だけHEX20」がコスト効率最高のアプローチだ。
実務チェックリスト
HEX20のチェックリストをお願いします。
HEX20は「作れる形状では最強」。チェックリストも比較的シンプルですね。
HEX20自体のトラブルは少ない。問題の大半はメッシュ生成にある。メッシュさえ作れれば、HEX20は最も信頼できる3次元要素だ。
ターボチャージャーハウジング解析
自動車用ターボチャージャーのコンプレッサーハウジングは複雑な曲面形状を持つが、HEX20要素で約15万要素のモデルを構築すると、テトラ要素モデル（60万要素）と同等の応力精度を2019年のIHI技術報告で確認。計算時間は1/4に削減された。このため設計初期の最適化ループではHEX20が選択される。
20節点六面体要素(HEX20)のソフトウェア比較
HEX20の各ソルバーでの特徴
各ソルバーのHEX20要素に差はありますか？
基本的な精度は同等。差が出るのは付加機能と大規模性能。
| 機能 | Abaqus C3D20R | Nastran CHEXA(20) | Ansys SOLID186 |
|---|---|---|---|
| 低減積分 | 明示的に選択 | デフォルト | KEYOPT(2)で選択 |
| ハイブリッド | C3D20RH | — | u-P対応 |
| 温度依存 | 完全対応 | 対応 | 完全対応 |
| 大規模並列 | MPI対応 | MPI + GPU | MPI + GPU |
| 複合材 | — | — | SOLID186 Layered |
AnsysのSOLID186 Layeredって何ですか？
AnsysのSOLID186はLayeredオプションで積層構造（複合材）に対応する。通常のHEX20に積層方向と各層の材料を定義できる。複合材のソリッドモデリングに便利だ。
メッシュ生成ツールの比較
| ツール | HEX20メッシュ品質 | 特徴 |
|---|---|---|
| HyperMesh | 最高（手動制御） | HEXメッシュのプロ向け |
| Ansys Meshing Multi-zone | 良好（半自動） | GUI操作で使いやすい |
| Abaqus/CAE Sweep | 良好 | Sweepメッシュが安定 |
| Cubit/Trelis | 良好 | 研究用。高度な分割機能 |
| GMSH | 基本的 | TET→HEX変換は限定的 |
HEXメッシュ生成ではHyperMeshが最強ですか？
手動HEXメッシュの品質と柔軟性ではHyperMeshが圧倒的だ。自動車や航空宇宙の衝突モデルで数百万のHEX要素を手動で作る技能は、HyperMeshのスキルそのものだ。
選定ガイド
まとめると？
HEX20は「メッシュが作れるなら最高」という結論は変わりませんね。
そう。HEX20の要素技術自体はどのソルバーでも十分成熟している。ボトルネックは常にメッシュ生成だ。
HEX20の各ソルバー対応状況
Abaqus C3D20/C3D20R、ANSYS SOLID186、NX Nastran CHEXA（20節点）、LS-DYNA ELFORM=3はいずれも20節点六面体だが入力形式が異なる。Abaqusは節点番号の並び順がANSYSと異なり、混在モデルの変換時にHyperMeshのOptiStructブロック変換で誤配列が生じる事例が2015年頃に頻出し、サポートFAQの常連問題となった。
20節点六面体要素(HEX20)の先端研究
HEX20の先端研究
HEX20に関する最新の研究はありますか？
HEX20自体は成熟した要素だが、メッシュ生成の自動化と高次化の研究がある。
HEX27（Lagrange型）vs. HEX20（Serendipity型）
HEX27要素はHEX20と比べてどうですか？
HEX27は面中心ノード（6個）と体中心ノード（1個）を追加した完全二次の六面体。Serendipity型のHEX20より形状関数が「完全」で、数学的には望ましい。
ただし実用上の差は小さい。HEX27はDOFが27×3=81で、HEX20の60より35%多い。精度向上はわずかで、DOF増加に見合わない。だからHEX20（Serendipity型）が実務標準。
Spectral Element Method
スペクトル要素法は高次の多項式（p=4〜8）を使うFEMの一種。HEX要素に高次のGLL（Gauss-Lobatto-Legendre）点を配置する。地震波伝播のシミュレーション（SPECFEM3D等）で標準的に使われている。
HEX20のp=2を、p=8くらいまで上げるイメージですか？
そう。スペクトル要素法は対角質量マトリクスが自然に得られるため、陽解法の時間積分が効率的。地震学の大規模波動シミュレーションでは事実上唯一の選択肢だ。
IGA（等幾何解析）のHEX
等幾何解析（IGA）では、通常のLagrange/Serendipity基底の代わりにNURBS/B-spline基底を使う。IGAのHEX要素は：
- CAD形状と完全に一致（メッシュの形状近似誤差ゼロ）
- $C^1$ 以上の連続性（通常のFEMは $C^0$）
- 応力場が滑らか
応力のスムーズさはメリットが大きいですね。
通常のHEX20では要素間で応力が不連続（$C^0$）だが、IGAでは $C^1$ 以上の連続性により応力が要素間で滑らか。応力集中部の評価精度が向上する。
まとめ
HEX20の先端研究、まとめます。
HEX20は成熟した要素だが、高次化やIGAで新しい可能性が開けている。
非線形変形とHEX20の限界
HEX20要素は大変形解析で体積ロッキングが生じやすい。1984年にSimoらはF-bar法を提案し体積変形を独立変数として扱うことでロッキングを解消した。Abaqus C3D20HはHybrid formulation（混合変形・圧力）を採用し、ゴム・金属塑性の大変形でHEX20Rより10〜20%精度が高いとされる。
20節点六面体要素(HEX20)のトラブル対応
HEX20のトラブル
HEX20でもトラブルは起きますか？
HEX20は最も安定した要素の一つだが、メッシュ品質に起因するトラブルがある。
ヤコビアン負（要素の退化）
ヤコビアンが負になるのはどんなときですか？
中間節点の位置が不適切な場合だ。HEX20の中間節点は辺の中点付近にあるべきだが、曲面へのスナップで辺の1/4や3/4の位置に移動すると、ヤコビアンが負になることがある。
経験則として、中間節点が辺の長さの25%〜75%の範囲内に収まるべき。これを超えるとヤコビアンが負になるリスクが高い。
対策は？
完全積分での体積ロッキング
HEX20の完全積分（C3D20）で体積ロッキングは起きますか？
$\nu > 0.49$ 程度で起き得る。C3D20の27積分点は多すぎて体積拘束が過剰になることがある。
対策：
- C3D20R（低減積分）に切り替え — 8積分点で拘束が緩和
- C3D20RH（ハイブリッド+低減） — 非圧縮材の最安定選択
結局C3D20Rが最も安全ということですね。
そう。C3D20（完全積分）を使う理由は「低減積分のアワーグラスが心配」な場合だけだが、HEX20の低減積分ではアワーグラスは実質問題にならない。C3D20Rをデフォルトにして問題ない。
DOFが多すぎる
HEX20モデルのDOFが多くてメモリ不足になります。
HEX20は1要素あたり60自由度。HEX8（24自由度）の2.5倍。
対策：
- HEX8I（非適合モード）への切り替え — 精度はC3D20Rの70〜80%だがDOFは半分以下
- サブモデリング — 着目部位のみHEX20、全体はHEX8またはTET10
- 反復法ソルバー — 直接法よりメモリ効率が良い
「全体HEX8I + 着目部位HEX20のサブモデリング」が効率的なアプローチですね。
まさにそう。HEX20の精度が必要な部位は通常限定的だ。全体をHEX20にする必要はなく、サブモデリングで最適化するのが実務的だ。
まとめ
HEX20のトラブル対処、整理します。
HEX20のトラブルはHEX8より少ないですね。要素自体の品質が高いから。
その通り。HEX20のトラブルの大半はメッシュ生成の段階で起きる。要素の計算自体はほとんど問題を起こさない。良いメッシュさえ作れば、HEX20は最高の結果を出す。
ミッドサイド節点の品質確認
HEX20のミッドサイド節点が辺中点から25%以上ずれると要素のヤコビアンが負になり解析エラーとなる。商用メッシャーでは品質チェック指標「Jacobian Ratio」が1.0〜0.6の範囲推奨で、特にインポートSTEPデータから自動メッシュ生成した場合は確認必須だ。Hypermeshでは2000年代からこの指標を標準品質チェックに含めている。
関連トピック
なった
詳しく
報告