4節点四面体要素(TET4)
4節点四面体要素(TET4)の理論基礎
TET4要素とは
先生、TET4って3次元解析の最も基本的な要素ですか？
最も基本的ではあるが、実務で使うべきではない要素だ。これを先に言っておく。理由は後で詳しく説明するが、まず理論を理解しよう。
いきなり「使うな」宣言ですか…。
TET4の理論を学ぶことは重要だ。FEMの基礎を理解するには最適な要素だからね。ただし実務で精度のある結果が欲しいならTET10を使うべきだ。
形状関数
TET4は4つの節点を持つ四面体で、各節点に3自由度（$u, v, w$）。合計12自由度。
形状関数は体積座標（barycentric coordinates）$L_1, L_2, L_3, L_4$ で表される：
変位の補間：
形状関数が線形（1次）だから、変位も要素内で線形に変化するんですね。
そう。変位が線形ということは、ひずみが要素内で一定。これがTET4の致命的な弱点だ。
定ひずみの問題
ひずみが一定だと何が問題なんですか？
応力が急変する部位（応力集中部、曲げ応力の勾配がある部位）で、応力の変化を全く表現できない。
例えば曲げを受ける梁をTET4でモデル化すると：
- 理論上は断面の上面が引張、下面が圧縮で線形分布
- TET4では各要素内が一定応力だから、階段状の粗い近似にしかならない
- 正確な曲げ応力を出すには極めて細かいメッシュが必要
TET10（二次要素）なら要素内でひずみが線形に変化するから、曲げも表現できる。
その通り。TET10は10節点で各辺に中間節点を持ち、二次の形状関数を使う。要素内でひずみが線形変化するから、曲げの応力分布を1要素の中で表現できる。
TET4の剛性マトリクス
TET4の要素剛性マトリクスは12×12で、B行列が定数行列だから数値積分なしに閉じた形で計算できる：
$V_e$ は要素体積、$[B]$ はひずみ-変位行列（定数）、$[D]$ は構成則マトリクス。
数値積分が不要なのはメリットですね。
計算は速い。ただし精度が低いから、結局メッシュを細かくする必要があり、DOF数が増えて計算時間のメリットが相殺される。
定量的な精度比較
TET4とTET10の精度差を数値で示してもらえますか？
片持ち梁（$L/h = 10$）の先端たわみを理論値と比較：
| 要素タイプ | メッシュ | DOF | たわみ誤差 |
|---|---|---|---|
| TET4 | 粗い（100要素） | 200 | -40%（過小） |
| TET4 | 中程度（1,000要素） | 1,500 | -15% |
| TET4 | 細かい（10,000要素） | 12,000 | -5% |
| TET10 | 粗い（100要素） | 1,200 | -2% |
| TET10 | 中程度（1,000要素） | 9,000 | -0.3% |
応力精度ではさらに差が広がる。TET4は応力が要素内一定だから、応力集中のピーク値はメッシュをどんなに細かくしても理論値に近づくのが遅い。
TET4を使うべき場面
TET4を使う正当な理由はありますか？
ほとんどない。強いて言えば：
- 陽解法（explicit）の衝撃解析 — TET4は質量マトリクスが単純で、陽解法の安定時間増分に有利な場合がある（LS-DYNAの一部ケース）
- 流体メッシュとの結合 — CFDコードがTET4ベースの場合、構造側もTET4で整合させることがある
- プリプロセッサの制約 — 一部の古い自動メッシュ生成器がTET4しか出力できない
ただしいずれも「やむを得ない場合」であり、TET10が使えるなら常にTET10を使うべきだ。
まとめ
TET4の理論を整理します。
要点：
- 4節点、線形形状関数、要素内ひずみ一定 — FEMの最も基本的な3D要素
- 精度が低い — 曲げや応力集中を表現するにはTET10の5〜10倍のDOFが必要
- 実務では使わない — TET10が代替として常に優れる
- 理論学習には最適 — FEMの原理を理解するための教科書的な要素
- 「TET4は使うな」がFEMの鉄則 — これを知っているだけで解析品質が劇的に上がる
「使わない」理由を理解するために学ぶ、というのが面白いですね。
TET4の限界を知ることはFEMの精度を理解する上で非常に重要だ。なぜTET4がダメなのか説明できるエンジニアは、FEMを正しく使えるエンジニアだ。
TET4の線形形状関数と定ひずみ特性
4節点四面体要素（TET4）は体積座標L1〜L4を用いた線形形状関数を持ち、要素内でひずみが一定（Constant Strain Tetrahedron, CST）になる。1943年にクーラント（R. Courant）が三角形要素の概念を提示し、1956年にターナーらが3次元四面体へ拡張した。最もシンプルな3Dソリッド要素として実装が容易だが、ひずみが一定なため曲げや応力集中の再現に要素数が膨大に必要となる欠点がある。
4節点四面体要素(TET4)の数値計算手法
TET4の実装詳細
TET4の実装は他の要素と比べてシンプルですか？
最もシンプルだ。B行列が定数なので数値積分が不要。FEMプログラミングの入門として最適だ。
B行列の計算
4節点の座標 $(x_i, y_i, z_i)$ から体積座標の勾配を計算し、B行列を構成する。TET4のB行列は6×12の定数行列：
ここで $b_i, c_i, d_i$ は節点座標から計算される定数。$V$ は四面体の体積。
$V$ はどう計算しますか？
4節点の座標から：
$V > 0$ なら節点の順序が正しい（右手系）。$V < 0$ なら要素が裏返っている。
ソルバー別の要素名
| ソルバー | 要素名 | 備考 |
|---|---|---|
| Nastran | CTETRA（4節点版） | PSOLIDプロパティ |
| Abaqus | C3D4 | 完全積分（積分点1つ） |
| Ansys | SOLID285 / SOLID185（退化） | SOLID285はu-p形式 |
| LS-DYNA | *ELEMENT_SOLID（TET4） | 陽解法で使用 |
AbaqusのC3D4は積分点が1つ…ということは要素中心で1点だけ評価？
そう。B行列が定数だからどこで評価しても同じ結果。1点積分で厳密。これはTET4の数学的な美しさだが、実用上は精度の低さに直結する。
TET4の改良版
TET4を改良した要素はありますか？
いくつかの改良が試みられている：
平均化ひずみ法（MINI要素）
TET4の中心に追加のバブル関数を導入し、要素内の変位場を豊かにする。非圧縮問題での体積ロッキング対策。AbaqusのSOLID285（Ansys）はこの系統。
S-TET4（Smoothed TET4）
Enhanced Assumed Strain (EAS)
TET4に追加のひずみモードを導入。ロッキングを回避しつつ精度を向上。一部の研究コードで実装。
改良TET4はTET10の代替になりますか？
精度はTET4より大幅に向上するが、TET10には及ばないことが多い。改良TET4の主な利点は既存のTET4メッシュを活かせること。メッシュを作り直す手間がない。
TET4メッシュの品質管理
TET4メッシュの品質指標は？
| 指標 | 理想値 | 許容範囲 |
|---|---|---|
| アスペクト比 | 1.0 | < 5.0 |
| 体積比（理想体積/実体積） | 1.0 | > 0.1 |
| ヤコビアン | 正 | 必ず正（負は要素反転） |
| 最小内角 | 70.5° | > 10° |
最小内角が10°以下の「ペシャンコ」な四面体は精度が出ないんですね。
TET4は要素形状への敏感性が高い。歪んだ要素は精度が著しく低下する。TET10のほうが要素形状の歪みに対して頑健だ。これもTET10を推奨する理由の一つ。
まとめ
TET4の実装詳細、整理します。
要点：
- B行列が定数 — 数値積分不要。実装が最もシンプル
- $V < 0$ で要素反転 — メッシュ品質の基本チェック
- 改良TET4（S-TET4, EAS等） — 精度向上するがTET10には及ばない
- 要素形状への敏感性が高い — 歪んだTET4は精度が著しく低下
- FEMプログラミングの入門に最適 — 実装して学ぶべき要素
TET4の閉形式積分
TET4の剛性マトリクスは体積座標の積分公式∫L1^a L2^b L3^c L4^d dV = (a!b!c!d!)/(a+b+c+d+3)! × 6Vを使って閉形式で解ける。この公式によりガウス積分なしに厳密な剛性マトリクスを計算できるため計算効率が高い。要素体積V=(1/6)|[x2-x1,y2-y1,z2-z1;...;x4-x1,y4-y1,z4-z1]|の行列式で求まり、Vが負になると要素頂点の定義順序が逆（Left-hand）であることを示す。
4節点四面体要素(TET4)の実務適用
TET4の実務的な扱い
TET4は使うなと言いつつ、実務で遭遇することはありますか？
ある。以下のケースでTET4を目にすることがある：
1. レガシーモデル — 古い解析で作られたTET4メッシュが残っている
2. 自動メッシュの不注意 — プリプロセッサのデフォルトがTET4になっていた
3. 陽解法の大規模モデル — LS-DYNAで計算コスト削減のためTET4を使用
4. CFD連成 — 流体側のメッシュがTET4ベースで、構造側も合わせた
TET4メッシュに遭遇したときの対処
既存のTET4メッシュがあるとき、どうすればいいですか？
3つの選択肢：
1. TET10に変換 — 各辺に中間節点を追加。多くのプリプロセッサで自動変換可能
2. メッシュを作り直す — TET10で新規メッシュ生成
3. そのまま使うが結果を慎重に解釈 — 応力値は信用せず、力の流れの確認のみ
TET4→TET10の変換は簡単ですか？
要素の辺の中点に節点を追加するだけなので、ほとんどのプリプロセッサで1クリックで変換できる。曲面上の中間節点は曲面にスナップさせる設定が必要（そうしないと中間節点がCAD面から外れる）。
Abaqusでは、C3D4メッシュからC3D10Mへの変換が *NGEN コマンドで可能。HyperMeshではelem typeの変更で一括変換できる。
TET4の結果をどう使うか
TET4で解析してしまった結果は全く使えないんですか？
使い方次第だ：
- 変位のトレンド — 定性的にはおおむね正しい。力の流れの確認に使える
- 反力 — 正確。要素タイプに依存しない（平衡条件）
- 応力のピーク値 — 信用できない。特に応力集中部は過小評価
- 平均応力（広い領域の平均） — おおむね正しい
「応力のピーク値は信用できない」がTET4の最大の問題ですね。
設計で最も重要なのが応力のピーク値（応力集中、最大von Mises応力）だから、TET4はまさに「最も重要な結果が出せない」要素なんだ。
ベンチマーク問題での比較
TET4とTET10の差を体感するためのベンチマーク問題は？
最も分かりやすいのは片持ち梁の先端たわみだ。
手順：
1. $L = 100$ mm, $h = 10$ mm, $b = 10$ mm の片持ち梁
2. 先端に集中荷重 $P = 100$ N
3. 理論解: $\delta = PL^3/(3EI) = 100 \times 100^3 / (3 \times 200000 \times 833.3) = 0.02$ mm
4. TET4とTET10で同じメッシュサイズで解析して比較
これなら手軽に試せますね。
FEMを始めたら最初にやるべき演習だ。TET4で解いて「なぜ理論値と合わないのか」を考え、TET10に変えて「なるほど」と理解する。この体験がFEMの精度に対する感覚を育てる。
実務チェックリスト
TET4に関するチェックリストをお願いします。
「TET4でないことを確認する」が最初のステップ。これだけで解析品質が保証される。
初心者がFEMで最も多く犯すミスが「TET4を使ってしまうこと」だ。このチェックを習慣化するだけで、FEM解析のレベルが1段上がる。
TET4の適用限界と使用場面
TET4は曲げ変形の精度が低いため、曲げが支配的な薄肉構造には使えない。実務上はメッシュ接続の「接着剤要素」として、HEX要素ゾーンとTET10ゾーンの境界に挟む用途で使われる場合がある。溶接熱影響部（HAZ）の再結晶組織像のように形状が複雑すぎてHEXメッシュを作れない部位の「最後の手段」として機能するが、応力精度を確認するには必ず隣接精密要素との比較検証が必要だ。
4節点四面体要素(TET4)のソフトウェア比較
各ソルバーのTET4の扱い
各ソルバーではTET4をどう位置づけていますか？
どのソルバーでも「非推奨」なんですね。
Abaqusのマニュアルには明確に「C3D4 should generally not be used as a stress/displacement element」と書いてある。Nastranも同様。ソルバーベンダー自身がTET4の使用を推奨していない。
メッシュ生成ソフトのデフォルト設定
自動メッシュでTET4がデフォルトになっているソフトはありますか？
注意が必要だ：
- GMSH — デフォルトは TET4 (order 1)。
Mesh.ElementOrder = 2でTET10に変更 - Ansys Meshing — デフォルトは要素の次数「Program Controlled」で通常TET10
- Abaqus/CAE — TET10 (C3D10M) がデフォルト
- HyperMesh — 要素タイプを手動で指定
GMSHのデフォルトがTET4なのは罠ですね。
オープンソースのGMSH + CalculiXの組み合わせで解析する初心者がTET4のまま解析して「精度が出ない」と悩むケースが多い。必ず ElementOrder = 2 を設定すること。
選定ガイド
TET4に関する選定ガイドは？
「TET4は使わない」が唯一の選定ガイドですね（笑）。
そう。TET4に関する最良のアドバイスは「TET10を使え」だ。これ以上シンプルなガイドラインはない。
TET4の主要ソルバー実装と推奨事項
TET4はNastranのCTETRA(4)、AbaqusのC3D4、AnsysのSOLID285（ハイブリッド線形四面体）として実装されている。ほぼすべての商用コードがTET4を実装しているが、構造応力解析への使用には警告が付く。AnsysのMechanical2023ではMeshing設定でTET4をデフォルト無効化し、常にTET10を推奨する設定に変更された。LS-DYNAでは明示的衝突解析でTET4を使う場合、付加質量法（ELFORM=10）で砂時計を防ぐ特殊定式がある。
4節点四面体要素(TET4)の先端研究
TET4の改良研究
TET4を改良する研究はまだ続いているんですか？
活発だ。TET4メッシュは自動生成が容易だから、もしTET4の精度をTET10並みに上げられれば、メッシュ生成の大幅な効率化になる。
Smoothed FEM (S-FEM)
S-FEMとは？
シンガポール国立大学のLiu G.R.教授のグループが開発した手法だ。通常のFEMではひずみを要素内で計算するが、S-FEMではひずみを要素間の「平滑化ドメイン」で計算する。
TET4に適用した場合：
- ES-FEM (Edge-based S-FEM) — 辺をベースにした平滑化。TET4でTET10並みの精度
- NS-FEM (Node-based S-FEM) — 節点をベースにした平滑化。上界推定（変位が大きめ）
- FS-FEM (Face-based S-FEM) — 面をベースにした平滑化
ES-FEMでTET4がTET10並みになるなら、実用的ですね。
理論的にはそう。ただし商用ソルバーへの実装が進んでいない。研究コード（GEOFEM等）では使えるが、Abaqus/Nastranのような主要ソルバーには組み込まれていない。
Virtual Element Method (VEM)
仮想要素法（VEM）は任意多角形/多面体要素を使えるFEMの拡張だ。TET4の代わりに、Voronoi分割による多面体要素を使うことで、メッシュの自由度を高める。
四面体にこだわらなくていいんですか？
VEMでは五角形や六角形の面を持つ多面体要素も使える。形状の自由度が高いため、メッシュ生成が容易になる。研究段階だが将来的にはTET4/TET10の代替になる可能性がある。
Immersed Boundary法
Immersed Boundary FEMはCAD形状を「背景メッシュに埋め込む」手法だ。背景メッシュは単純なHEX8やTET4の格子で、CAD形状は背景メッシュを切断する。切断された要素には特殊な積分手法を適用する。
メッシュ生成そのものを不要にする発想ですね。
そう。FCM（Finite Cell Method）やCutFEMがこの系統。CADから直接FEM解析に移行する「メッシュレスCAE」のビジョンを体現する手法で、活発に研究されている。
まとめ
TET4の先端研究、まとめます。
TET4の「精度が低い」問題は、要素技術の革新で解決される可能性がある。ただし現時点ではTET10を使うのが最も確実だ。
TET4の安定化と体積ロッキング
TET4は完全非圧縮（ν=0.5）材料では体積ロッキングが顕著で、変位が実解の1/100以下になることがある。1980年代にHugheらが提案したF-bar法（体積平均F-bar）を線形TET4に適用すると体積ロッキングを大幅に軽減できる。最新のSelectively Reduced Integration（SRI）TET4はロッキングなしで四面体メッシュの速度を維持できるとして、生体力学（骨格筋シミュレーション）コードで実用化されている。
4節点四面体要素(TET4)のトラブル対応
TET4のトラブル
TET4を使ってしまった場合のトラブル対処を教えてください。
TET4最大のトラブルは「使ってしまったこと」自体だが、具体的な症状と対処を見ていこう。
応力が理論値の半分以下
曲げ問題で最大応力が理論値の50%しか出ません。
TET4の典型的症状だ。定ひずみ要素は曲げの応力勾配を表現できないため、最大応力を大幅に過小評価する。
対策：
1. TET10に変換 — 最善の策。精度が劇的に改善
2. メッシュを5倍に細分化 — TET4のままでも精度は上がるが、計算コストが膨大
3. 結果を使わない — TET4の応力は設計判断に使えないと割り切る
応力コンターがガタガタ
応力コンターが要素ごとにバラバラの色で、滑らかではありません。
TET4は要素内でひずみ一定だから、隣接要素間で応力が大きくジャンプする。非平均化（unaveraged）の応力コンターを見るとチェッカーボードパターンになる。
平均化（averaged）すれば滑らかになりますか？
見た目は滑らかになるが、平均化で情報が失われる。特に応力集中のピーク値は平均化で過小評価される。「滑らかなコンター = 正確」ではない。非平均化コンターの不連続が大きい場所はメッシュが粗い証拠だ。
変位は合うのに応力が合わない
たわみは理論値に近いのに、応力が大きくずれます。
FEMは変位法だから、変位は比較的早く収束するが、応力（変位の微分）は収束が遅い。TET4では特にこの差が顕著だ。
一般に：
- 変位の誤差 ∝ $h^2$（メッシュサイズの2乗）
- 応力の誤差 ∝ $h$（メッシュサイズの1乗）
TET4の場合、変位は比較的合理的な精度が出ても、応力は1オーダー精度が低い。
体積ロッキング
TET4で非圧縮材料を解析すると変位が異常に小さくなります。
体積ロッキングだ。TET4は各要素に1つの体積拘束条件が入るため、$\nu \to 0.5$ で自由度が過剰に拘束される。
対策：
- TET10M（Modified）に変換 — AbaqusのC3D10Mは体積ロッキングに対策済み
- TET4のu-p形式 — AnsysのSOLID285。圧力を独立変数にしてロッキング回避
- $\nu = 0.499$ にする — 完全な非圧縮を避ける（応急処置）
TET4を検出する方法
モデルにTET4が含まれているかどうかを確認する方法は？
.inp ファイルで *ELEMENT, TYPE=C3D4 を検索。C3D10やC3D10Mなら安全他人のモデルを受け取ったとき、最初にやるべきことですね。
その通り。「TET4検出」はモデルレビューの第一ステップだ。TET4が見つかったら、解析者にTET10への変換を依頼する。これだけで解析品質が保証される。
まとめ
TET4のトラブル対処、整理します。
.inp or .bdf で要素タイプを確認。レビューの第一ステップ結論: TET4の全てのトラブルはTET10に変換することで解決する。
これ以上簡潔なトラブルシューティングガイドはないだろう。TET4→TET10。これだけだ。
TET4の過剛性（剛性過剰）問題の確認
TET4は同一自由度数でHEX8より約4〜10倍剛性が高く出る（過剛性）。これは定ひずみ特性によるもので、メッシュ収束しても実解より変位が小さくなり続ける。診断法はTET10同一モデルとの比較：TET4結果がTET10比で変位20%以上小さい場合、TET4の過剛性が解析精度に影響しているとみなす。Abaqusではソルバー設定でC3D4（TET4）を使った場合のWARNING: ELEMENT TYPE NOT RECOMMENDED FOR STRESS ANALYSISが自動出力される。
関連トピック
なった
詳しく
報告