6節点三角形要素(TRIA6)
理论与物理
LST要素 — CST的升级版
6节点三角形单元（TRIA6, LST）是解决了CST问题的单元吗？
是的。LST（线性应变三角形）是一种单元内应变呈线性变化的二次三角形单元。它完全解决了CST的常应变问题。
形函数
TRIA6拥有3个顶点节点和3个边中点节点。形函数是面积坐标 $L_1, L_2, L_3$ 的二次多项式：
TRIA6拥有3个顶点节点和3个边中点节点。形函数是面积坐标 $L_1, L_2, L_3$ 的二次多项式：
顶点节点：
边中点节点：
是TET10的二维版本呢。
正是如此。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 |
|---|---|---|
| 自动网格划分 | 容易 | 稍困难 |
| 单位自由度精度 | Q8的70〜80% | 基准 |
| 曲边对应 | 有边中点 | 有边中点 |
| 对网格质量的敏感性 | 比Q8更稳健 | 对扭曲敏感 |
| 自适应网格 | 容易 | 稍困难 |
TRIA6在精度效率上不如Q8，但在网格生成容易度上胜出。和TET10 vs. HEX20是同样的构图呢。
理解得很完美。无论是二维还是三维，“三角形/四面体自动网格划分容易，四边形/六面体精度高”这个构图是共通的。
数值积分
| 积分 | 点数 | 精度 |
|---|---|---|
| 3点Gauss（完全积分） | 3 | 精确积分二次多项式 |
| 1点Gauss | 1 | 精度不足（不使用） |
| 7点Gauss | 7 | 用于高阶积分 |
通常3点就足够了呢。
是的。TRIA6的标准是3点完全积分。不会发生剪切自锁（因为是二次单元）。体积自锁通常也不是问题。这是一个非常稳定的单元。
总结
我来整理一下TRIA6的理论。
要点：
- 二次三角形单元（LST） — CST的升级版。应变呈线性变化
- TET10的2D版 — 继承了自动网格划分的优点
- 收敛速度与Q8同等 — 但单位自由度效率Q8更高
- 对网格质量稳健 — 即使在扭曲形状下也比Q8稳定
- 3点Gauss积分足够 — 无自锁
如果说CST是“不要用”的单元，那么TRIA6就是“可以放心使用”的单元呢。
正是如此。TRIA6是二维自动网格划分分析中最可靠的三角形单元。
TRIA6的二次完全多项式与精度
6节点三角形单元（TRIA6）拥有3个顶点节点和3个边中点节点，其形函数使用完全二次多项式（6项）。在面积坐标中，顶点形函数为Li(2Li-1)，边中点为4LiLj。TRIA3（线性）的3倍计算成本，使得弯曲问题中的应力误差收敛率从h的2次方改善到3次方。据记载，1960年代后期，Browning和Abrahamson等人首次将二次三角形单元实际应用于飞机机翼的应力分析。
各项的物理意义
- 惯性项（质量项）：$\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³） |
| 力 $F$ | N（牛顿） | mm系为N，m系也统一为N |
数值解法与实现
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公式。
中间节点的处理
与Q8/TET10类似，TRIA6的中间节点CAD曲线捕捉也很重要。它能提高曲面边界的近似精度。
注意点：
- 中间节点需位于边的25%〜75%范围内
- 雅可比矩阵需为正
- 曲率大的部分需加密网格
TRIA6的网格质量
| 指标 | 理想值 | 容许范围 |
|---|---|---|
| 长宽比 | 1.0 | < 5.0 |
| 最小角度 | 60° | > 20° |
| 雅可比矩阵 | 正 | 必须为正 |
比CST的质量要求宽松吗？
TRIA6比CST对单元形状扭曲更稳健。二次形函数能部分吸收扭曲。但最小角度低于10°的极端扭曲会导致精度下降。
总结
我来整理一下TRIA6的数值方法。
要点：
- 3点Gauss积分为标准 — 完全积分，无自锁
- 中间节点的CAD捕捉 — 与Q8, TET10注意事项相同
- 对形状扭曲稳健 — 比CST稳定
- CST→TRIA6转换能显著提高精度 — 一次→二次的效果
TRIA6的3点高斯积分
TRIA6通常使用3点或7点高斯积分。3点方案的积分点位于三角形重心等距的位置（例如L1=L2=1/6，L3=2/3等），可以精确积分到5次多项式。7点方案是Strang-Fix（1973年）提出的高精度版本，可以完全积分到7次多项式。对于应力采样，3个积分点是超收敛点，因此广泛使用Zienkiewicz-Zhu法从这些点外推应力来求节点应力。
线性单元（一次单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切自锁（减缩积
相关主题
なった
詳しく
報告