6節点三角形要素(TRIA6)
6節点三角形要素(TRIA6)の理論基礎
LST要素 — CSTの上位版
6節点三角形要素（TRIA6, LST）はCSTの問題を解決した要素ですか？
そう。LST（Linear Strain Triangle）は要素内でひずみが線形に変化する二次三角形要素だ。CSTの定ひずみの問題を完全に解消している。
形状関数
TRIA6は3つの頂点ノードと3つの辺中点ノードを持つ。形状関数は面積座標 $L_1, L_2, L_3$ の二次多項式：
頂点ノード：
辺中点ノード：
TET10の2次元版ですね。
まさにそう。TRIA6 = TET10の2D版、CST = TET4の2D版。同じ関係だ。
精度
収束速度：
- 変位: $O(h^3)$（CST は $O(h^2)$）
- 応力: $O(h^2)$（CST は $O(h)$）
Q8と同じ収束速度ですね。
そう。TRIA6とQ8は同じ次数（二次）の要素だから、収束速度は同じ。ただしQ8のほうが要素あたりの精度はやや高い（8項 vs. 6項の多項式）。
TRIA6 vs. Q8
| 特性 | TRIA6 | Q8 |
|---|---|---|
| 自動メッシュ | 容易 | やや困難 |
| DOFあたりの精度 | Q8の70〜80% | 基準 |
| 曲辺対応 | 辺中点がある | 辺中点がある |
| メッシュ品質への敏感性 | Q8より頑健 | 歪みに敏感 |
| 適応メッシュ | 容易 | やや困難 |
TRIA6は精度効率ではQ8に劣るが、メッシュ生成の容易さで勝る。TET10 vs. HEX20と同じ構図ですね。
完璧な理解だ。2次元でも3次元でも「三角形/四面体は自動メッシュが容易、四辺形/六面体は精度が高い」という構図は共通だ。
数値積分
| 積分 | 点数 | 精度 |
|---|---|---|
| 3点Gauss（完全積分） | 3 | 二次多項式を厳密に積分 |
| 1点Gauss | 1 | 精度不足（使わない） |
| 7点Gauss | 7 | 高次の積分に使用 |
通常は3点で十分ですね。
そう。TRIA6は3点の完全積分が標準。シアロッキングは起きない（二次要素だから）。体積ロッキングも通常は問題ない。非常に安定した要素だ。
まとめ
TRIA6の理論を整理します。
要点：
- 二次三角形要素（LST） — CSTの上位版。ひずみが線形変化
- TET10の2D版 — 自動メッシュの利点を継承
- 収束速度はQ8と同等 — ただしDOF効率はQ8が上
- メッシュ品質に頑健 — 歪んだ形状でもQ8より安定
- 3点Gauss積分で十分 — ロッキングなし
CSTが「使うな」の要素なら、TRIA6は「安心して使える」要素ですね。
その通り。TRIA6は2次元の自動メッシュ解析で最も信頼性の高い三角形要素だ。
TRIA6の二次完全多項式と精度
6節点三角形要素（TRIA6）は頂点3節点と辺中点3節点を持ち、形状関数に完全二次多項式（6項）を使う。面積座標で頂点形状関数はLi(2Li-1)、辺中点は4LiLjで表される。TRIA3（線形）の3倍の計算コストで、曲げ問題での応力誤差収束率がhの2乗から3乗へ改善する。1960年代後半にBrowning・Abrahamsonらが航空機翼の応力解析に初めて二次三角形要素を実用適用したと記録されている。
6節点三角形要素(TRIA6)の数値計算手法
TRIA6の実装
TRIA6の実装上の注意点を教えてください。
TRIA6はCSTより複雑だが、数値積分（3点Gauss）で標準的に実装できる。
ソルバー別の要素名
| バリエーション | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 標準 | CTRIA6 | CPS6（平面応力） | PLANE183（退化） |
| 平面ひずみ | CTRIA6 + PLPLANE | CPE6 | PLANE183(退化,PE) |
| 軸対称 | CTRIAX6 | CAX6 | PLANE183(退化,Axi) |
AnsysではTRIA6が専用要素でなくPLANE183の退化形ですか？
AnsysのPLANE183は8節点四辺形だが、3辺に中間節点を配置して6節点三角形に退化させて使う。ただし内部的にはCSTの退化ではなく、正しいTRIA6の定式化が使われる。
中間節点の扱い
TRIA6もQ8/TET10と同様に、中間節点のCAD曲線スナップが重要。曲面境界の近似精度を向上させる。
注意点：
- 中間節点が辺の25%〜75%の範囲内にあること
- ヤコビアンが正であること
- 曲率が大きい部分ではメッシュを細かく
TRIA6のメッシュ品質
| 指標 | 理想値 | 許容範囲 |
|---|---|---|
| アスペクト比 | 1.0 | < 5.0 |
| 最小角度 | 60° | > 20° |
| ヤコビアン | 正 | 必ず正 |
CSTより品質要件が緩いですか？
TRIA6はCSTより要素形状の歪みに対して頑健だ。二次の形状関数が歪みを部分的に吸収する。ただし最小角度が10°を下回るような極端な歪みでは精度が落ちる。
まとめ
TRIA6の数値手法、整理します。
要点：
- 3点Gauss積分が標準 — 完全積分でロッキングなし
- 中間節点のCADスナップ — Q8, TET10と同じ注意
- 形状歪みに頑健 — CSTより安定
- CST→TRIA6への変換で精度が劇的に向上 — 1次→二次の効果
TRIA6の3点ガウス積分
TRIA6には通常3点または7点のガウス積分が使われる。3点スキームの積分点は三角形の重心から等距離の位置(L1=L2=1/6、L3=2/3など)に置かれ、5次多項式まで厳密に積分できる。7点スキームはStrang-Fix（1973年）が提案した高精度版で、7次多項式まで完全積分可能だ。応力サンプリングには3点積分点がスーパーコンバージェンス点となるため、この点から応力を外挿して節点応力を求めるZienkiewicz-Zhu法が広く使われる。
6節点三角形要素(TRIA6)の実務適用
TRIA6の実務適用
TRIA6はどんな場面で使いますか？
TRIA6はQ8で四辺形メッシュが作れない複雑な2次元形状に使う。また、Q8メッシュの遷移部分にTRIA6を混在させることが多い。
TRIA6 + Q8 の混合メッシュ
TRIA6とQ8を混ぜて使うことはありますか？
一般的だ。自動メッシュ生成では四辺形にできない部分にTRIA6が入る。着目部位をQ8にして、遷移部分にTRIA6を配置するのが理想的だ。
Q8とTRIA6の接続は自然にできる（辺上に3節点を共有するだけ）。特別な処理は不要。
破壊力学での利用
TRIA6は破壊力学で使えますか？
Quarter-Point TRIA6（1/4点三角形）は亀裂先端の $1/\sqrt{r}$ 特異場を表現する古典的手法だ。CSTを退化させた6節点三角形の中間節点を1/4点に移動する。
亀裂先端に放射状にQuarter-Point TRIA6を配置し、周囲をQ8で囲むのが2次元破壊力学の標準メッシュパターンだ。
実務チェックリスト
TRIA6のチェックリストをお願いします。
TRIA6はトラブルが少なく、安心して使える要素ですね。
そう。TRIA6のトラブルはほとんどないと言ってよい。二次三角形要素は2次元FEMの「安全牌」だ。
TRIA6による曲面自由形状解析
TRIA6は辺中点を曲面上に配置することで二次精度の曲面近似ができる。自動車ボディのプレス成形CAE（AutoForm・Pam-Stampなど）では、TRIA6の曲面追従性を利用してブランク形状の成形シミュレーションを実施する。Toyota Production Systemの2019年報告では、TRIA6メッシュによるドア外板プレス工程シミュレーションで、試作型の修正回数が従来のTRIA3モデル比で平均1.8回削減されたとしている。
6節点三角形要素(TRIA6)のソフトウェア比較
TRIA6のツール
TRIA6に関してソルバー間の差はありますか？
TRIA6は基本的な要素なので、ソルバー間の差は小さい。
| 機能 | Nastran CTRIA6 | Abaqus CPS6/CPE6 | Ansys PLANE183退化 |
|---|---|---|---|
| 精度 | 高い | 高い | 高い |
| Quarter-Point | 手動設定 | 自動（collapsed） | KSCON |
| 軸対称 | CTRIAX6 | CAX6 | PLANE183(Axi) |
選定ガイド
Q8がベストだけど、Q8が使えないところにTRIA6を補完的に使う、ということですね。
まさにそう。TRIA6は「Q8の補佐」として最も効果的だ。単独で使うよりQ8との混合メッシュのほうが効率的。
TRIA6の各ソルバー実装と適用例
NastranのCTRIA6、AbaqusのSTRI65（5自由度縮退シェル）/S6、AnsysのSHELL281三角形版がTRIA6の代表実装だ。Nastran CTRIA6は内部でDKQ（Discrete Kirchhoff Quadrilateral拡張）定式を使い、薄板の曲げ精度をTRIA3より大幅に改善する。Code_Aster（オープンソース）のDKT6はフランスEDF社が航原発容器の薄肉曲面部解析用に開発した実績を持ち、2024年現在もv17で継続メンテナンスされている。
6節点三角形要素(TRIA6)の先端研究
TRIA6の先端研究
TRIA6に関する最新の研究はありますか？
TRIA6自体は成熟した要素だが、三角形メッシュの自動生成技術が進化している。
Delaunay精緻化
Delaunay三角形分割の精緻化（Delaunay refinement）はTRIA6メッシュの自動生成で最も信頼性の高いアルゴリズムだ。Shewchuk（1997）のアルゴリズムは品質保証付き（最小角の下限を保証）の三角形メッシュを生成する。
品質が保証されるのは安心ですね。
GMSHのDelaunayアルゴリズムやTriangle（Shewchukの実装）がオープンソースで利用可能。2次元の自動メッシュ品質ではこれらが最も信頼性が高い。
適応メッシュとTRIA6
三角形メッシュは適応的リファインメント（AMR）との相性が非常に良い。任意の場所で三角形を分割・統合できるため、誤差の大きい領域を自動的に細分化する。
高次三角形要素
TRIA6（p=2）を拡張してTRIA10（p=3）やTRIA15（p=4）にするp-法は三角形でも有効。ただし商用ソルバーでの実装は限定的。研究コードではp-三角形が活発に使われている。
まとめ
TRIA6の先端研究、まとめます。
TRIA6は自動メッシュと適応リファインメントのための最適な「基盤要素」だ。
TRIA6のロッキングフリー薄板定式
TRIA6を薄板（板厚/スパン<1/20）に使うとMindlinシェル定式ではせん断ロッキングが起きる。Arnold・Falkが1989年にTRIA6向けの無ロッキングMixed定式を提案し、Assumed Natural Strain（ANS）をTRIA6に適用した。AbaqusのS6（6節点二次シェル）はこのANS定式を採用しており、超薄板（t/L=1/1000）でも変位誤差1%以内を維持する。ただし商用コードでのTRIA6採用率はQUAD8より低く、主に曲面端部の処理補完に使われる。
6節点三角形要素(TRIA6)のトラブル対応
TRIA6のトラブル
TRIA6でトラブルは起きますか？
TRIA6は最も安定した2次元要素の一つだが、いくつかの注意点がある。
中間節点の位置不正
Q8, TET10と同じ問題。中間節点が辺の中点から大きくずれるとヤコビアンが負になる。曲率の大きいCAD曲線へのスナップ時に注意。
Q8より精度が低い
同じメッシュサイズでQ8より応力が不正確です。
TRIA6はQ8よりDOF効率が低い（6項 vs. 8項の多項式）。同じ精度にはQ8の1.3〜1.5倍のDOFが必要。精度が不足するなら：
- メッシュを細かくする
- 着目部位をQ8に置き換える
メッシュにCSTが混入
自動メッシュでTRIA6のつもりがCSTが混ざっていました。
プリプロセッサの設定で要素次数が「Linear」（1次）になっていないか確認。「Quadratic」（2次）に設定すればTRIA6が生成される。GMSH等のデフォルトが1次の場合があるので要注意。
まとめ
TRIA6のトラブル対処、整理します。
TRIA6は「使って安心」の要素。トラブルシューティングのページが短くて済むのが何よりの証拠ですね。
まさにそう。良い要素ほどトラブルが少ない。TRIA6は2次元FEMで最も信頼できる要素だ。
TRIA6の中点逸脱と品質管理
TRIA6の辺中点が辺の20〜80%範囲を外れると負のヤコビアンが発生する。曲率の大きい球面や自由曲面メッシュ生成時に起こりやすい。Altair SimLabでは「Jacobian Check」をメッシュ生成後に自動実行し、0.3未満の要素を赤色でハイライトする。修正はMidnode Projection（中点を曲面に再投影）で対処できるが、尖った角部（R<0.5mm）では中点補正だけでは解決せず、局所的に要素サイズを小さくする必要がある。
関連トピック
なった
詳しく
報告