10節点四面体要素(TET10)

TET10は3次元構造解析の実務標準だ。ほぼ任意の形状に自動メッシュが生成でき、精度も高い。現代の3次元FEMの大部分はTET10で行われていると言ってよい。

TET4の4つの頂点に加えて、各辺の中点に6つの中間節点を追加した10節点の四面体。各節点に3自由度で合計30自由度。

形状関数は二次多項式：

そう。これがTET4（ひずみ一定）との決定的な違いだ。曲げによる応力勾配を1つの要素内で表現できる。

テイラー展開で考えるとわかりやすい。真の変位場を多項式で近似するとき：

収束の速さ（メッシュを細かくしたときの誤差減少率）：

だからTET10は粗いメッシュでも十分な精度が出る。TET4と同じDOF数なら、TET10のほうが遥かに正確だ。

TET10は形状関数が二次だから、B行列は1次。$B^T D B$ は2次の多項式になり、これを正確に積分するには4点Gauss積分が必要だ。

ほとんどない。TET4と違って、TET10の完全積分（4点）で体積ロッキングは通常起きない。低減積分にすると精度が落ちるだけで利点が少ない。

ただしAbaqusのC3D10M（Modified）は特殊な積分手法を使う。ホバーグラス制御付きの低減積分で、接触問題での安定性が向上している。通常のC3D10とC3D10Mの使い分けは後で詳しく説明する。

中間節点が辺上にあるため、辺を曲線にできる。CADの曲面に中間節点をスナップさせることで、円筒や球面のような曲面をTET4よりはるかに正確に近似できる。

例えば円孔の周囲をTET4でモデル化すると、円が多角形で近似される。TET10なら各辺が曲線になり、円に沿う。応力集中の評価精度に直結する。

中間節点が直線的に配置された「直辺TET10」になり、曲面の近似精度がTET4並みに落ちる。多くのプリプロセッサでは自動的にCAD面にスナップするが、設定を確認すること。

要点：

その通り。完璧ではないが（HEX20のほうが効率は良い）、メッシュ生成の手間を含めた総合コストではTET10が最も優れている。

これは実務で最も重要な使い分けの一つだ。

AbaqusのマニュアルはC3D10Mを一般推奨としている。C3D10は接触面でチェッカーボードパターン（応力の振動）が出やすいのが弱点。C3D10Mはこの問題を解決した改良版だ。

ただしC3D10Mは1つのホバーグラスモードを持つため、非常にまれにホバーグラス変形が問題になることがある。C3D10のほうが「素直な」結果が出る場面もある。

NastranのCTETRA(10)とAnsysのSOLID187はそれぞれ独自の安定化手法を内蔵しており、C3D10Mに相当する対策が組み込まれている。名前が分かれていないだけで、内部的には改良されている。

1. 中間節点のCADスナップ — 曲面上の中間節点をCAD面に投影。これを忘れると曲面の近似精度が低下

隣接要素のサイズ比が1.5倍以下が推奨。2倍を超えるとインターフェースで精度が落ちる。多くのプリプロセッサの自動メッシュは漸変率を自動制御するが、極端なサイズ差がある場合は手動調整が必要。

適応的メッシュリファインメント（AMR）は、誤差推定に基づいてメッシュを自動細分化する手法だ。TET10は四面体の分割が容易だからAMRとの相性が良い。

手順：

Ansys WorkbenchのConvergence機能が最も手軽だ。目標応力の変化率を指定すると、自動的にメッシュを細分化して収束を確認する。Abaqusでは *ADAPTIVE MESH が流体に対応しているが、構造のh-refinementは手動が多い。

要点：

1. CADインポート＋簡略化 — 不要な形状特徴を除去

着目部位の形状特徴に応じて：

そう。フィレットの応力集中を捕捉するには十分なメッシュ密度が必要だ。フィレットを省略してエッジにすると応力特異点になるから、疲労評価が必要な部位ではフィレットを必ず保持すること。

2次元断面を押し出した形状（パイプ、Hビーム等）ではsweepメッシュでHEX要素を効率的に生成できる。TET10でこういう形状をメッシュすると不必要にDOF数が増える。

実務的なアプローチ：

非平均化応力の不連続が最も信頼できる指標だ。

隣接要素間の応力値を比較して：

Abaqus/CAE, Ansys Workbench, HyperViewなど主要なポストプロセッサは全て対応している。Abaqus/CAEでは「Average: 75%」等の設定で平均化のしきい値を変えられる。100%は完全平均化、0%は完全非平均化。

繰り返しになるが、それだけ重要だということだ。TET10を使い、メッシュ品質を確認し、収束性を検証する。この3ステップを守れば、3次元FEMの結果は信頼できる。

基本的な精度はほぼ同じだが、付加機能に差がある。

TET10は通常、陽解法では使いにくい（安定時間増分が小さくなりがち）。C3D10Mの改良積分は陽解法での安定性も考慮されており、衝突解析でもTET10が使える。

そう。TET10は各ソルバーが最も力を入れて開発している要素の一つだ。自動メッシュユーザーの大半がTET10を使うから、ソルバーベンダーも品質向上に注力している。

自動メッシュの品質向上と新しい計算手法がメインだ。

通常のh-法（メッシュを細かくして精度向上）に対して、p-法はメッシュを変えずに要素の多項式次数を上げて精度向上する。TET10（p=2）をp=3, p=4と上げていく。

p-法の利点：

Simcenter Nastranのp-法要素（SOL 147/148）やMSC Nastranのp-要素が対応している。Ansysにもp-法オプションがある。ただしp-法は通常のh-法のワークフローと異なるため、普及は限定的だ。

大変形解析で要素が歪んだ場合に、解析の途中でTET10メッシュを再生成する手法。ALE（Arbitrary Lagrangian-Eulerian）法の構造版。

鍛造や押し出し成形のような大変形問題で有効。AbaqusのALE適応メッシュやDefromの自動リメッシュが実用化されている。

GNNやCNNを使って：

理想的にはそう。ただし現時点では「人間のメッシュ設計能力をAIが超える」段階には至っていない。研究は活発で、数年以内に実用的なAIメッシュアシスタントが商用ソフトに組み込まれる可能性は高い。

四面体にこだわらず、任意の多面体要素を使う手法も研究されている。VEM（Virtual Element Method）やPolygonal FEMがこの系統だ。TET10の「メッシュ品質に敏感」という弱点を、要素形状の自由度で解決する発想。

TET10は優秀だがトラブルがゼロではない。知っておくべき問題がある。

中間節点がCAD面からずれている、または隣接要素の中間節点が一致しない場合に問題が起きる。

症状：

対策：

多くのソルバーはヤコビアンが負の要素をエラーとして拒否する。NastranのFATALエラーやAbaqusの「Negative eigenvalue」が出る。ヤコビアン負の要素が1つでもあると解析が走らない。

TET10はTET4の2.5倍のDOF（10節点 vs 4節点）、HEX8の約3倍のDOF。大規模モデルではDOF数が問題になる。

対策：

TET10とHEX20の接続にはピラミッド要素（PYRA13/14）を挟む。ピラミッド要素は四角形面（HEX側）と三角形面（TET側）を持つ遷移要素だ。主要なプリプロセッサで自動生成できる。

C3D10の接触面安定性問題だ。4点完全積分のC3D10では接触面の圧力が振動する。

対策：

Abaqusで接触を含むTET10解析ではC3D10Mを使わない理由がない。C3D10を使って接触が不安定になるのは予測可能な問題だ。

同じDOF数ではHEX20のほうが精度が高いのは事実だ。TET10でHEX20と同等精度を得るには2〜5倍のDOFが必要。

ただしメッシュ生成時間を含めた総コストではTET10が有利なことが多い。HEXメッシュに1日かかる形状でも、TET10自動メッシュは数分で完了する。

その通り。精度効率（DOFあたりの精度）はHEXが上。ワークフロー効率（メッシュ生成＋解析＋結果評価のトータル時間）はTET10が上。どちらを重視するかはプロジェクトによる。

この2つを押さえておけば、TET10で大きなトラブルは起きない。TET10は実務で最も信頼できる3D要素だ。