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$ で表される：
$$ N_i = L_i \quad (i = 1, 2, 3) $$
$$ u(x,y) = L_1 u_1 + L_2 u_2 + L_3 u_3 $$
B行列は定数で、数値積分なしに厳密に計算できる：
$$ [K_e] = A_e \cdot t \cdot [B]^T [D] [B] $$
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要素は現在も接続要素や教育用途に使われ続けている。
各項の物理的意味
- 慣性項（質量項）：$\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で統一 |
数値解法と実装
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年代まで航空機スキン解析のメインツールとして使われ続けた。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
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メッシュを手動リメッシュすると規定している。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
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つの問い
- 「何を解くか」：3節点三角形要素(TRIA3)に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
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にEdge-based Smoothed FEM（ES-FEM）を適用すると、ひずみを辺ベースで平滑化してCSTでTRIA6並みの精度が達成できる。2次元ではS-FEMの効果が3次元より顕著で、実用レベルの精度向上が報告されている。
CST（3節点）の中心にバブル関数を追加したMINI要素（3+1節点）は、流体力学のStokes問題で標準的に使われるが、構造力学でも非圧縮問題での体積ロッキング対策として研究されている。
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定式を内部で使用しており、薄板の面外変位解析で信頼性が高い。
トラブルシューティング
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の曲げ精度がわずかに改善される場合がある。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——3節点三角形要素(TRIA3)の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告