20節点六面体要素(HEX20)

カテゴリ: 構造解析 | 統合版 2026-04-06
CAE visualization for hex20 element theory - technical simulation diagram
20節点六面体要素(HEX20)

理論と物理

HEX20 — 最高精度の3次元要素

🧑‍🎓

HEX20はFEMの3次元要素で最も精度が高いんですか?


🎓

実用的な要素としては最高精度と言ってよい。20節点の二次六面体要素で、各節点に3自由度、合計60自由度。HEX8(24自由度)の2.5倍のDOFだが、精度は桁違いに高い。


形状関数

🎓

HEX20は8つの頂点ノードと12個の辺中点ノードを持つ。形状関数は三次曲面(serendipity)型:


頂点ノード:

$$ N_i = \frac{1}{8}(1+\xi_i\xi)(1+\eta_i\eta)(1+\zeta_i\zeta)(\xi_i\xi+\eta_i\eta+\zeta_i\zeta-2) $$

辺中点ノード(例: $\xi_i = 0$):

$$ N_i = \frac{1}{4}(1-\xi^2)(1+\eta_i\eta)(1+\zeta_i\zeta) $$

🧑‍🎓

Serendipity型って何ですか? Lagrange型とは違うんですか?


🎓

Lagrange型(27節点HEX)は面と内部にも節点を持つが、Serendipity型(20節点HEX)は辺上のみに中間節点を持ち、面の中心や内部に節点がない。Serendipity型はDOFが少ないのに精度はLagrange型とほぼ同等で、実用上はSerendipity型が圧倒的に多い。


HEX20の精度

🧑‍🎓

HEX20はどの程度正確ですか?


🎓

TET10との比較で:


同じDOF数での精度HEX20TET10
変位非常に高い高い
応力非常に高い高い
効率(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が不利な場面は?


    🎓
    • 複雑な形状 — 自動HEXメッシュが困難。TET10のほうが現実的
    • 大変形解析 — HEX20は歪みに弱い(中間節点が移動しやすい)
    • 陽解法 — DOFが多く安定時間増分が小さい。HEX8Rのほうが効率的

    • まとめ

      🧑‍🎓

      HEX20の理論を整理します。


      🎓

      要点:


      • 20節点、Serendipity型、二次六面体 — 実用3D要素の最高精度
      • C3D20R(低減積分)が最推奨 — 8積分点で最高効率
      • 収束が速い — TET10より1オーダー速い
      • メッシュ生成が課題 — 自動HEXは困難。Sweep/マップドメッシュ向き
      • 大変形・陽解法には不向き — HEX8RやTET10のほうが安定

      🧑‍🎓

      HEX20は「作れるなら最高の要素」ということですね。


      🎓

      まさにそう。「最高の精度」と「作れるならば」の間にメッシュ生成の壁がある。この壁を越えられるエンジニアにとって、HEX20は最強の武器だ。


      Coffee Break よもやま話

      二次六面体要素の理論的優位性

      20節点六面体要素(Serendipity要素)は1966年にErgatoudisらが提案したアイソパラメトリック定式化の代表作だ。完全2次多項式を含む形状関数により、一次要素の8倍近い要素内での応力精度を持つ。1970年代のNASAアポロ計画以降の宇宙機構造解析で採用され、「金の標準要素」と呼ばれた時代がある。

      各項の物理的意味
      • 慣性項(質量項):$\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で統一

      数値解法と実装

      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が圧倒的に多い。


      ソルバー別の要素名

      バリエーションAbaqusNastranAnsys
      完全積分C3D20CHEXA(20)SOLID186(full)
      低減積分C3D20RSOLID186(red.)
      ハイブリッドC3D20H, C3D20RHu-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で十分

      Coffee Break よもやま話

      27点vs14点ガウス積分の比較

      HEX20要素の完全積分には3×3×3=27ガウス点が必要だが、計算コストが高い。1974年にBarlow点理論により14点(reduced integration)でも精度を保てることが示され、ANSYS SOLID186では3×3×3がデフォルトだが、Abaqus C3D20Rは2×2×2(8点)の縮減積分を採用し計算時間を約40%短縮しつつ精度を維持する。

      線形要素(1次要素)

      節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意(低減積分やB-bar法で緩和)。

      2次要素(中間節点付き)

      曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨:応力評価が重要な場合。

      完全積分 vs 低減積分

      完全積分:過剰拘束(ロッキング)のリスク。低減積分:アワーグラスモード(零エネルギーモード)のリスク。適材適所で選択。

      アダプティブメッシュ

      誤差指標(ZZ推定量等)に基づく自動細分化。応力集中部の精度を効率的に向上。h法(要素分割)とp法(次数増加)がある。

      ニュートン・ラフソン法

      非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。

      修正ニュートン・ラフソン法

      接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。

      収束判定基準

      力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$(一般に $\epsilon = 10^{-3}$〜$10^{-6}$)。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$

      荷重増分法

      全荷重を一度に負荷せず、小刻みに増加させる。弧長法(Riks法)は荷重-変位関係の極値点を越えて追跡可能。

      直接法 vs 反復法のたとえ

      直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す(直接法)より、見当をつけて開き、前後に調整する(反復法)方が効率的なのと同じ原理。

      メッシュの次数と精度の関係

      1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。

      実践ガイド

      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のチェックリストをお願いします。


      🎓
      • [ ] C3D20R(低減積分)を使用しているか
      • [ ] 中間節点がCAD曲面にスナップしているか
      • [ ] ヤコビアンが全要素で正か
      • [ ] 曲率が大きい部分でメッシュが十分細かいか
      • [ ] 板厚方向に最低2要素あるか
      • [ ] 大変形の場合、HEX8への切り替えを検討したか

      • 🧑‍🎓

        HEX20は「作れる形状では最強」。チェックリストも比較的シンプルですね。


        🎓

        HEX20自体のトラブルは少ない。問題の大半はメッシュ生成にある。メッシュさえ作れれば、HEX20は最も信頼できる3次元要素だ。


        Coffee Break よもやま話

        ターボチャージャーハウジング解析

        自動車用ターボチャージャーのコンプレッサーハウジングは複雑な曲面形状を持つが、HEX20要素で約15万要素のモデルを構築すると、テトラ要素モデル(60万要素)と同等の応力精度を2019年のIHI技術報告で確認。計算時間は1/4に削減された。このため設計初期の最適化ループではHEX20が選択される。

        解析フローのたとえ

        解析の流れは、実は料理とそっくりです。まず材料を買い出し(CADモデルの準備)、下ごしらえをして(メッシュ生成)、火にかけて(ソルバー実行)、最後に盛り付ける(後処理で可視化)。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう? 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。

        初心者が陥りやすい落とし穴

        あなたはメッシュ収束性を確認していますか? 「計算が回った=結果が正しい」と思っていませんか? これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。

        境界条件の考え方

        境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら? どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。

        ソフトウェア比較

        HEX20の各ソルバーでの特徴

        🧑‍🎓

        各ソルバーのHEX20要素に差はありますか?


        🎓

        基本的な精度は同等。差が出るのは付加機能と大規模性能。


        機能Abaqus C3D20RNastran CHEXA(20)Ansys SOLID186
        低減積分明示的に選択デフォルトKEYOPT(2)で選択
        ハイブリッドC3D20RHu-P対応
        温度依存完全対応対応完全対応
        大規模並列MPI対応MPI + GPUMPI + 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のスキルそのものだ。


        選定ガイド

        🧑‍🎓

        まとめると?


        🎓
        • 精密応力評価(圧力容器等) → C3D20R + Sweepメッシュ
        • ベンチマーク検証 → C3D20R(参照解として)
        • 大規模HEXモデルNastran CHEXA(20)(スケーラビリティ)
        • 複合材のソリッドモデル → Ansys SOLID186 Layered
        • HEXメッシュの手動作成HyperMesh

        • 🧑‍🎓

          HEX20は「メッシュが作れるなら最高」という結論は変わりませんね。


          🎓

          そう。HEX20の要素技術自体はどのソルバーでも十分成熟している。ボトルネックは常にメッシュ生成だ。


          Coffee Break よもやま話

          HEX20の各ソルバー対応状況

          Abaqus C3D20/C3D20R、ANSYS SOLID186、NX Nastran CHEXA(20節点)、LS-DYNA ELFORM=3はいずれも20節点六面体だが入力形式が異なる。Abaqusは節点番号の並び順がANSYSと異なり、混在モデルの変換時にHyperMeshのOptiStructブロック変換で誤配列が生じる事例が2015年頃に頻出し、サポートFAQの常連問題となった。

          選定で最も重要な3つの問い

          • 「何を解くか」:20節点六面体要素(HEX20)に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
          • 「誰が使うか」:初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車(GUI)とMT車(スクリプト)の違いに似ている。
          • 「どこまで拡張するか」:将来の解析規模拡大(HPC対応)、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。

          先端技術

          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の先端研究、まとめます。


          🎓
          • HEX27 vs. HEX20 — 理論的にはHEX27が優位だが、DOF効率でHEX20が実用的
          • Spectral Element — 高次HEXで波動伝播。地震学の標準
          • IGA — NURBS基底で形状近似誤差ゼロ。応力の連続性が向上

          • HEX20は成熟した要素だが、高次化やIGAで新しい可能性が開けている。


            Coffee Break よもやま話

            非線形変形とHEX20の限界

            HEX20要素は大変形解析で体積ロッキングが生じやすい。1984年にSimoらはF-bar法を提案し体積変形を独立変数として扱うことでロッキングを解消した。Abaqus C3D20HはHybrid formulation(混合変形・圧力)を採用し、ゴム・金属塑性の大変形でHEX20Rより10〜20%精度が高いとされる。

            トラブルシューティング

            HEX20のトラブル

            🧑‍🎓

            HEX20でもトラブルは起きますか?


            🎓

            HEX20は最も安定した要素の一つだが、メッシュ品質に起因するトラブルがある。


            ヤコビアン負(要素の退化)

            🧑‍🎓

            ヤコビアンが負になるのはどんなときですか?


            🎓

            中間節点の位置が不適切な場合だ。HEX20の中間節点は辺の中点付近にあるべきだが、曲面へのスナップで辺の1/4や3/4の位置に移動すると、ヤコビアンが負になることがある。


            🎓

            経験則として、中間節点が辺の長さの25%〜75%の範囲内に収まるべき。これを超えるとヤコビアンが負になるリスクが高い。


            🧑‍🎓

            対策は?


            🎓
            • 曲率が大きい部分でメッシュを細かくする
            • プリプロセッサの「mid-node projection limit」を設定
            • ヤコビアンが負の要素を手動で修正

            • 完全積分での体積ロッキング

              🧑‍🎓

              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 + 着目部位HEX20のサブモデリング」が効率的なアプローチですね。


              🎓

              まさにそう。HEX20の精度が必要な部位は通常限定的だ。全体をHEX20にする必要はなく、サブモデリングで最適化するのが実務的だ。


              まとめ

              🧑‍🎓

              HEX20のトラブル対処、整理します。


              🎓
              • ヤコビアン負 → 中間節点のスナップ量を制限。メッシュを細かく
              • 体積ロッキング → C3D20R(低減積分)に切り替え
              • DOF過大サブモデリング。全体はHEX8I、局所のみHEX20
              • C3D20Rがデフォルト — 完全積分C3D20を使う必要はほぼない

              • 🧑‍🎓

                HEX20のトラブルはHEX8より少ないですね。要素自体の品質が高いから。


                🎓

                その通り。HEX20のトラブルの大半はメッシュ生成の段階で起きる。要素の計算自体はほとんど問題を起こさない。良いメッシュさえ作れば、HEX20は最高の結果を出す。


                Coffee Break よもやま話

                ミッドサイド節点の品質確認

                HEX20のミッドサイド節点が辺中点から25%以上ずれると要素のヤコビアンが負になり解析エラーとなる。商用メッシャーでは品質チェック指標「Jacobian Ratio」が1.0〜0.6の範囲推奨で、特にインポートSTEPデータから自動メッシュ生成した場合は確認必須だ。Hypermeshでは2000年代からこの指標を標準品質チェックに含めている。

                「解析が合わない」と思ったら

                1. まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
                2. 最小再現ケースを作る——20節点六面体要素(HEX20)の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
                3. 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
                4. 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
                関連シミュレーター

                この分野のインタラクティブシミュレーターで理論を体感しよう

                シミュレーター一覧

                関連する分野

                熱解析製造プロセス解析V&V・品質保証
                この記事の評価
                ご回答ありがとうございます!
                参考に
                なった
                もっと
                詳しく
                誤りを
                報告
                参考になった
                0
                もっと詳しく
                0
                誤りを報告
                0
                Written by NovaSolver Contributors
                Anonymous Engineers & AI — サイトマップ
                プロフィールを見る