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はどの程度正確ですか？
HEX20のほうが1オーダー速く収束する。メッシュが半分の粗さでも同等精度…。
だからHEXメッシュが作れる形状では、HEX20が最も効率的だ。問題はメッシュ生成の手間。この「精度効率 vs. メッシュ生成コスト」のトレードオフが、HEX20とTET10の選択の本質だ。
積分スキーム
HEX20の低減積分は2×2×2 = 8点ですか。HEX8の完全積分と同じ点数。
そう。HEX20の低減積分（C3D20R）はHEX8の完全積分と同じ8つの積分点で、はるかに高い精度を出す。C3D20Rは3次元FEMで最も効率的な要素の一つと評される。
C3D20Rにアワーグラスの問題はありますか？
HEX20の低減積分にはアワーグラスモードが1つだけ存在する。HEX8Rの12個と比べて極めて少なく、実用上ほとんど問題にならない。ただし1要素のパッチテストでは注意が必要。
いつHEX20を使うか
HEX20を使うべき場面は？
- 最高精度が必要 — ベンチマーク検証、精密応力評価
- 板厚方向の応力分布が重要 — 圧力容器の応力分類
- 接触問題の精密評価 — 接触面の圧力分布
- メッシュ生成が容易な形状 — 回転体、押し出し体
逆にHEX20が不利な場面は？
- 複雑な形状 — 自動HEXメッシュが困難。TET10のほうが現実的
- 大変形解析 — HEX20は歪みに弱い（中間節点が移動しやすい）
- 陽解法 — DOFが多く安定時間増分が小さい。HEX8Rのほうが効率的
まとめ
HEX20の理論を整理します。
要点：
- 20節点、Serendipity型、二次六面体 — 実用3D要素の最高精度
- C3D20R（低減積分）が最推奨 — 8積分点で最高効率
- 収束が速い — TET10より1オーダー速い
- メッシュ生成が課題 — 自動HEXは困難。Sweep/マップドメッシュ向き
- 大変形・陽解法には不向き — HEX8RやTET10のほうが安定
HEX20は「作れるなら最高の要素」ということですね。
まさにそう。「最高の精度」と「作れるならば」の間にメッシュ生成の壁がある。この壁を越えられるエンジニアにとって、HEX20は最強の武器だ。
NASAとNASTRAN — FEMの夜明け
今や世界中で使われている有限要素法ソルバー「NASTRAN」は、1960年代にNASAが開発しました。アポロ計画でロケットの構造解析が必要だったのです。当時のコンピュータはメモリ数KBの時代——今のスマートフォンの100万分の1以下の性能で、人類を月に送る構造計算をしていたのです。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値例：片持ち梁の先端荷重（L=1m, 断面50×100mm, 鋼材E=210GPa, P=1kN）
最大たわみ δ = PL³/(3EI) = 1000×1000³/(3×210000×4,166,667) ≈ 0.381 mm 最大応力 σ = PL×(h/2)/I ≈ 12.0 MPa（降伏応力235MPaに対して安全率19.6）
メッシュ密度を変えた収束性の確認：
ポイント：要素数を4倍にしても結果は0.5%しか変わらない→8,000要素で十分収束。これが「メッシュ収束性」の確認です。
構造解析の収束問題や計算コストに課題を感じていませんか？ — Project NovaSolverは、実務者が日々直面するこうした課題の解決を目指す研究開発プロジェクトです。
CAEの未来を、実務者と共に考える
Project NovaSolverは、20節点六面体要素(HEX20)における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →