3节点三角形单元(TRIA3)
3节点三角形单元(TRIA3)的理论基础
CST单元 — FEM的起点
3节点三角形单元也称为CST（恒定应变三角形），对吧。
是的。应变在单元内恒定，因此称为「恒定应变三角形」。1956年由Turner, Clough, Martin, Topp发表的FEM的首个单元，是有限元法的历史本身。
FEM的首个单元！ 但是在TET4页面上说「不要使用一阶单元」。
CST是TET4的二维版本。具有相同的缺点 — 应变恒定，无法表现应力梯度。实务中应避免使用。但对于理解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的二维版 — 相同问题，相同对策
TRIA3单元的CST定式化
3节点三角形单元（TRIA3）也称为恒定应变三角形（CST），1956年由Turner、Clough、Martin、Topp在《航空学杂志》上发表，是FEM的最早单元。用面积坐标L1、L2、L3（L1+L2+L3=1）定义的线性形状函数，使得单元内应变为恒定。这一最简单的二维单元至今仍被用于连接单元和教学用途。
3节点三角形单元(TRIA3)的数值计算手法
CST的实现
CST的实现真的很简单吗？
在FEM单元中最简单。B矩阵是由三个节点坐标计算得出的常数矩阵：
仅由节点坐标差决定。无需数值积分。
单元面积 $A$ 为：
$A > 0$ 时节点为逆时针，$A < 0$ 时为顺时针。
求解器的单元名称
| 求解器 | 单元名 | 备注 |
|---|---|---|
| Nastran | CTRIA3 | 用PSHELL指定板厚 |
| Abaqus | CPS3（平面应力）, CPE3（平面应变） | 最基本的单元 |
| Ansys | PLANE182（退化） | 4节点单元的1条边退化成3节点 |
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求得。单点积分（重心点）等效，故不存在积分精度问题。正是因为这种简单性，直到1980年代，TRIA3仍是飞机蒙皮分析的主要工具。
3节点三角形单元(TRIA3)的实务应用
CST的实务处理
遇到CST时的处理方法是什么？
与TET4页面的方法完全相同：
1. 转换为TRIA6 — 在每条边上增加中间节点
2. 重新用四边形网格 — 如果可能，基于QUAD4重新网格划分
3. 继续使用但不信任应力 — 仅用于力流验证
允许CST混入的情形
即使混入部分CST也没关系的情况有吗？
在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个单元的patch进行手工计算，与求解器结果对比，这样的演习非常有效。
实务检查清单
关于CST的检查清单？
与TET4一样，「不使用」是最好的策略。
CST = TET4的二维版。对策相同。「转换为二阶单元」就是全部答案。
TRIA3的应用与四边形补充
实务中很少独立用TRIA3进行结构分析，主要用作QUAD4网格的「填缝」单元。汽车车身曲面网格中，TRIA3通常被限制在全单元数的5～10%以下。保时捷AG（2018年公开资料）规定，车体静刚度分析中TRIA3比例限制在全单元的3%以下，超过此值则采用手动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）
将Edge-based Smoothed FEM（ES-FEM）应用于CST，通过边基平滑化应变，可达到CST实现TRIA6级精度。二维问题中S-FEM效果比三维更显著，已报告实用级精度改进。
MINI单元
在CST（3节点）中心增加气泡函数的MINI单元（3+1节点），在流体力学Stokes问题中标准应用，结构力学中也作为非压缩问题的体积锁定对策而被研究。
虚拟单元法
VEM（虚拟单元法）不限于三角形，可用任意多边形定式。作为本质上超越CST限制（3节点、恒定应变）的手法，备受瞩目。
总结
整理CST的前沿研究。
CST是「不应使用的单元」，但相关研究仍在继续。如果ES-FEM或VEM在未来商用求解器中实现，一阶网格也许会有复兴的一天。
TRIA3的精度改进：发展到DKT
为解决TRIA3的弯曲精度不足，Batoz、Bathe、Ho在1980年开发了Discrete Kirchhoff Triangle（DKT）。DKT保留TRIA3的形状函数，在单元三条边的3个点上离散施加Kirchhoff条件（法线转角=横向位移微分），相对通常TRIA3，弯曲精度约提高10倍。Nastran的CTRIA3内部采用此DKT定式，薄板的面外位移分析可信性很高。
3节点三角形单元(TRIA3)的故障排除
CST故障
请简要说明CST故障的处理。
CST故障完全与TET4相同的结构。
所有故障的解决方案
| 故障 | 原因 | 解决方案 |
|---|---|---|
| 应力过小 | 恒定应变无法表现梯度 | 转换为TRIA6 |
| 等值线呈阶梯状 | 单元内应力恒定 | 转换为TRIA6 |
| 位移符合但应力不符 | 应力收敛慢 | 转换为TRIA6 |
| 体积锁定 | 一阶单元的本质局限 | 转换为TRIA6 |
| 网格不收敛 | CST收敛慢 | 转换为TRIA6 |
全部都是「转换为TRIA6」（笑）。
没有比这更简洁的故障排除指南。CST的所有故障都通过转换为TRIA6解决。与TET4→TET10完全相同。一阶→二阶单元是万能药。
转换为TRIA6的方法
大多数前处理器可一键转换：
- HyperMesh — Element Order → 二阶 变更
- 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的β（二阶剪切修正系数）从0.833到0.86可略微改进TRIA3的弯曲精度。
相关主题
价值
详细
错误