4節点四面体要素(TET4)
理論と物理
TET4要素とは
先生、TET4って3次元解析の最も基本的な要素ですか？
最も基本的ではあるが、実務で使うべきではない要素だ。これを先に言っておく。理由は後で詳しく説明するが、まず理論を理解しよう。
いきなり「使うな」宣言ですか…。
TET4の理論を学ぶことは重要だ。FEMの基礎を理解するには最適な要素だからね。ただし実務で精度のある結果が欲しいならTET10を使うべきだ。
形状関数
TET4は4つの節点を持つ四面体で、各節点に3自由度（$u, v, w$）。合計12自由度。
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ソリッド要素として実装が容易だが、ひずみが一定なため曲げや応力集中の再現に要素数が膨大に必要となる欠点がある。
各項の物理的意味
- 慣性項（質量項）：$\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で統一 |
数値解法と実装
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）であることを示す。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
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メッシュを作れない部位の「最後の手段」として機能するが、応力精度を確認するには必ず隣接精密要素との比較検証が必要だ。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
各ソルバーの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）で砂時計を防ぐ特殊定式がある。
選定で最も重要な3つの問い
- 「何を解くか」：4節点四面体要素(TET4)に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
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の先端研究、まとめます。
S-FEMとは？
シンガポール国立大学のLiu G.R.教授のグループが開発した手法だ。通常のFEMではひずみを要素内で計算するが、S-FEMではひずみを要素間の「平滑化ドメイン」で計算する。
TET4に適用した場合：
ES-FEMでTET4がTET10並みになるなら、実用的ですね。
理論的にはそう。ただし商用ソルバーへの実装が進んでいない。研究コード（GEOFEM等）では使えるが、Abaqus/Nastranのような主要ソルバーには組み込まれていない。
仮想要素法（VEM）は任意多角形/多面体要素を使えるFEMの拡張だ。TET4の代わりに、Voronoi分割による多面体要素を使うことで、メッシュの自由度を高める。
四面体にこだわらなくていいんですか？
VEMでは五角形や六角形の面を持つ多面体要素も使える。形状の自由度が高いため、メッシュ生成が容易になる。研究段階だが将来的にはTET4/TET10の代替になる可能性がある。
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はロッキングなしで四面体メッシュの速度を維持できるとして、生体力学（骨格筋シミュレーション）コードで実用化されている。
トラブルシューティング
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が自動出力される。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——4節点四面体要素(TET4)の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告