3節点三角形要素(TRIA3)
3節点三角形要素(TRIA3)の理論基礎
CST要素 — FEMの出発点
3節点三角形要素はCST（Constant Strain Triangle）とも呼ばれるんですよね。
そう。ひずみが要素内で一定だから「定ひずみ三角形」。1956年にTurner, Clough, Martin, Toppが発表したFEM最初の要素であり、有限要素法の歴史そのものだ。
FEM最初の要素！ でもTET4のページで「1次要素は使うな」と言われましたよね。
CSTはTET4の2次元版だ。同じ欠点を持つ — ひずみが一定で応力勾配を表現できない。実務では使うべきでない。ただしFEMの原理を理解するには最適な要素だ。
形状関数
CSTの形状関数は面積座標 $L_1, L_2, L_3$ で表される：
B行列は定数で、数値積分なしに厳密に計算できる：
$A_e$ は三角形の面積、$t$ は板厚。
要素面積と板厚を掛けるだけ。計算が非常にシンプルですね。
FEMプログラミングの入門として最適だ。CSTを手計算で組み立てる演習は、剛性マトリクスの意味を深く理解するのに役立つ。
CSTの限界
CSTの限界をまとめると：
- 曲げが表現できない — ひずみ一定なので曲げの応力勾配が出ない
- 収束が遅い — 精度を出すにはQ8やTRIA6の5〜10倍の要素が必要
- 応力が要素内一定 — コンターが階段状になり、応力集中が正確に評価できない
「使うな」という結論はTET4と全く同じですね。
CST（2D）= TET4（3D）。同じ問題、同じ対策（二次要素に変更）だ。
CSTを使ってよい場面
唯一の正当な用途は教育とFEMの基礎理解だ。CSTで手計算を行い、剛性マトリクスの組み立て、境界条件の適用、連立方程式の解法を体験する。これはFEMエンジニアの必修科目だ。
まとめ
CST（TRIA3）の理論を整理します。
要点：
- FEM最初の要素（1956年） — 有限要素法の歴史的出発点
- 定ひずみ — 要素内で応力一定。応力勾配を表現できない
- 実務では使わない — TRIA6（二次三角形）かQ8を使う
- 教育用として最適 — 手計算でFEMの原理を学べる
- TET4の2次元版 — 同じ問題、同じ対策
TRIA3要素のCST定式化
3節点三角形要素（TRIA3）はConstant Strain Triangle（CST）とも呼ばれ、FEM最初の要素として1956年にターナー・クラフ・マーティン・トップが「Journal of the Aeronautical Sciences」に発表した。面積座標L1,L2,L3（L1+L2+L3=1）で定義される線形形状関数により、要素内のひずみが一定になる。この最単純な2D要素は現在も接続要素や教育用途に使われ続けている。
3節点三角形要素(TRIA3)の数値計算手法
CSTの実装
CSTの実装は本当にシンプルですか？
FEM要素の中で最もシンプルだ。B行列が3つの節点座標のみから計算される定数行列：
節点座標の差だけで決まる。数値積分も不要。
要素面積 $A$ は：
$A > 0$ なら節点が反時計回り、$A < 0$ なら時計回り。
ソルバー別の要素名
| ソルバー | 要素名 | 備考 |
|---|---|---|
| Nastran | CTRIA3 | PSHELLで板厚指定 |
| Abaqus | CPS3（平面応力）, CPE3（平面ひずみ） | 最も基本的な要素 |
| Ansys | PLANE182（退化） | 4節点の1辺をつぶした退化形 |
Ansysでは専用のTRIA3要素がないんですか？
AnsysのPLANE182は4節点だが、3節点に退化させて使う。ただしこの退化形はCSTと同等の精度しかないので、やはり推奨しない。
自動メッシュでのCSTの出現
自動メッシュでCSTが出てくることはありますか？
Q4の自動メッシュで四辺形にできない部分にCSTが挿入されることがある。特にPaving法（四辺形メッシュの自動生成）のCloseout（閉じ込め）ステップで三角形が発生する。
対策：
- 三角形の数を最小限に抑える設定を使う
- 三角形が着目部位に入らないように制御
- 可能ならQ8の二次メッシュに切り替え
まとめ
CSTの実装と実務的扱い、整理します。
要点：
- B行列が定数 — 最もシンプルな実装。FEMプログラミングの入門
- $A < 0$ で要素反転 — 節点の順序を確認
- 自動メッシュで混入することがある — 着目部位に入らないよう制御
- 結論: TRIA6に変換 — CSTが混入したメッシュはTRIA6に変換すべき
TRIA3の閉形式剛性マトリクス
TRIA3の剛性マトリクスはK=t×A×B^T×D×B（t:板厚、A:要素面積、B:ひずみ-変位マトリクス、D:弾性マトリクス）の閉形式で計算でき、ガウス積分が不要だ。面積Aは座標行列の行列式の1/2で求まる。1点積分（重心点）と等価であるため、積分精度に関する議論が生じない。この単純さゆえ1980年代まで航空機スキン解析のメインツールとして使われ続けた。
3節点三角形要素(TRIA3)の実務適用
CSTの実務的扱い
CSTに遭遇したときの対処法を教えてください。
TET4のページと全く同じアプローチだ：
1. TRIA6に変換 — 各辺に中間節点を追加
2. Q8メッシュに作り直す — 可能なら四辺形ベースで
3. そのまま使うが応力は信用しない — 力の流れの確認のみ
CSTの混在が許容される場面
CSTが一部混ざっていてもOKな場面はありますか？
Q4メッシュの遷移部分に少数のCSTが混在する場合、着目部位から離れていれば実用上問題ない。Saint-Venantの原理で遠い場所の要素品質は結果に影響しにくい。
ただし着目部位にCSTが含まれていたら、結果は信用できない。Q4/Q8の領域だけで評価すること。
手計算でのCST演習
CSTの手計算演習は具体的にどうやりますか？
典型的な演習：
1. 3節点の座標を設定（例：(0,0), (1,0), (0,1)）
2. B行列を手計算
3. D行列（平面応力のフックの法則）を設定
4. $K = A \cdot t \cdot B^T D B$ を計算（6×6行列）
5. 境界条件を適用し、連立方程式を解く
6. 変位から応力を計算（$\sigma = D B u$）
6×6の連立方程式なら手計算できますね。
境界条件で自由度を消すとさらに小さくなる。2〜3要素のパッチで手計算し、ソルバーの結果と比較する演習は非常に有益だ。
実務チェックリスト
CSTに関するチェックリストは？
TET4と同じで「使わない」が最善の策ですね。
CST = TET4の2次元版。対策も同じ。「二次要素に変換」が全ての答えだ。
TRIA3の適用と四辺形補完
実務ではTRIA3単独で構造解析することは少なく、主にQUAD4メッシュの「隙間埋め」要素として用いられる。自動車ボディのサーフェスメッシュでは全要素数の5〜10%以下にTRIA3を抑えることが品質基準とされる。ポルシェAG（2018年公開資料）では車体静剛性解析において全要素の3%未満にTRIA3を制限し、それ以上になる場合はQUAD4メッシュを手動リメッシュすると規定している。
3節点三角形要素(TRIA3)のソフトウェア比較
CSTの各ソルバーでの扱い
各ソルバーはCSTをどう位置づけていますか？
全ソルバーが「使用は推奨しない」としている。
| ソルバー | 要素 | 備考 |
|---|---|---|
| Nastran | CTRIA3 | マニュアルで CTRIA6 を推奨 |
| Abaqus | CPS3 / CPE3 | ドキュメントで「粗い近似」と明記 |
| Ansys | PLANE182退化 | TRIA6推奨 |
| CalculiX | C2D3 | 使用可能だが精度は低い |
選定ガイド
「CSTは使わない」が選定ガイドの全てですね。
TET4と同じく、CSTの全てのトラブルは二次要素（TRIA6）に変換することで解決する。
TRIA3の主要ソルバー実装比較
NastranのCTRIA3（DKT/Shell定式）、AbaqusのS3/S3R、AnsysのSHELL181三角形オプション、Altair OptiStructのCTRIA3が代表的な実装だ。AbaqusのS3Rは自動的にS4Rの縮退版として扱われ、三角形にも砂時計制御が適用される。実際の処理時間はS3がS4の約70%だが、精度劣化を考慮すると積極的使用は推奨されない。Code_Asterでは薄板シェルにMECA_STATIQUE＋TRIA3の組み合わせで教育用ベンチマークが公式提供されている。
3節点三角形要素(TRIA3)の先端研究
CSTの改良研究
CSTを改良する研究はありますか？
TET4と同様に、CSTの精度向上は活発な研究分野だ。
ES-FEM（Edge-based Smoothed FEM）
CSTにEdge-based Smoothed FEM（ES-FEM）を適用すると、ひずみを辺ベースで平滑化してCSTでTRIA6並みの精度が達成できる。2次元ではS-FEMの効果が3次元より顕著で、実用レベルの精度向上が報告されている。
MINI要素
CST（3節点）の中心にバブル関数を追加したMINI要素（3+1節点）は、流体力学のStokes問題で標準的に使われるが、構造力学でも非圧縮問題での体積ロッキング対策として研究されている。
Virtual Element Method
VEM（仮想要素法）は三角形に限らず任意多角形で定式化できる。CSTの制約（3節点、定ひずみ）を本質的に超える手法として注目されている。
まとめ
CSTの先端研究、まとめます。
CSTは「使うべきでない要素」だが、FEMの基本として研究は続いている。将来的にES-FEMやVEMが商用ソルバーに実装されれば、1次メッシュの復権もあり得る。
TRIA3の精度向上：DKTへの発展
TRIA3の曲げ精度不足を解決するためBatoz・Bathe・Hoが1980年にDiscrete Kirchhoff Triangle（DKT）を開発した。DKTはTRIA3の形状関数を維持しながら、要素辺上の3点でキルヒホッフ条件（法線回転＝横変位微分）を離散的に課すことで、通常のTRIA3比で曲げ精度を約10倍改善する。NastranのCTRIA3はこのDKT定式を内部で使用しており、薄板の面外変位解析で信頼性が高い。
3節点三角形要素(TRIA3)のトラブル対応
CSTのトラブル
CSTのトラブル対処を手短に教えてください。
CSTのトラブルは全てTET4と同じ構造だ。
全てのトラブルの解決策
| トラブル | 原因 | 解決策 |
|---|---|---|
| 応力が過小 | 定ひずみで勾配表現不可 | TRIA6に変換 |
| コンターが階段状 | 要素内応力一定 | TRIA6に変換 |
| 変位は合うが応力が合わない | 応力の収束が遅い | TRIA6に変換 |
| 体積ロッキング | 1次要素の本質的限界 | TRIA6に変換 |
| メッシュ収束しない | CSTの収束が遅い | TRIA6に変換 |
全て「TRIA6に変換」ですね（笑）。
これ以上シンプルなトラブルシューティングガイドはない。CSTの全てのトラブルはTRIA6に変換することで解決する。TET4→TET10と全く同じ。1次要素→二次要素が万能薬だ。
TRIA6への変換方法
多くのプリプロセッサで1クリックで変換可能：
- HyperMesh — Element Order → 2次 に変更
- Abaqus/CAE — Element Type → CPS6 に変更
- Ansys — Element Order → Quadratic に変更
- GMSH — Mesh.ElementOrder = 2
変換は簡単。だったら最初からTRIA6を使えばいいですね。
まさにその通り。FEMの鉄則：常に二次要素を使え。
TRIA3の過剛性と応力チェック
TRIA3を多用したメッシュでは変位が過小評価（過剛性）になりやすい。同一形状をQUAD4で解いた結果と比べ、変位差15%以上が出る場合はTRIA3過多が疑われる。診断にはFEMツールのElement Type Contourでメッシュ内のTRIA3割合を可視化し、境界付近に集中していないか確認する。Nastranのβ（2シアー補正係数）を0.833から0.86に変更するとTRIA3の曲げ精度がわずかに改善される場合がある。
関連トピック
なった
詳しく
報告