3節点三角形要素(TRIA3)
理论与物理
CST要素 — FEM的起点
3节点三角形单元也被称为CST（常应变三角形），对吧。
是的。因为应变在单元内是恒定的，所以叫“常应变三角形”。这是1956年由Turner, Clough, Martin, Topp发表的FEM第一个单元，本身就是有限元法的历史。
FEM的第一个单元！但是在TET4的页面里，不是说“不要使用一阶单元”吗？
CST是TET4的二维版本。它们有相同的缺点——应变恒定，无法表现应力梯度。实际工作中不应使用。不过，对于理解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的二维版本 — 相同的问题，相同的对策
TRIA3单元的CST公式化
3节点三角形单元（TRIA3）也被称为常应变三角形（CST），作为FEM的第一个单元，于1956年由Turner、Clough、Martin、Topp在《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³） |
| 力 $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节点单元退化（压扁一条边）的形式 |
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年代它一直被用作飞机蒙皮分析的主要工具。
线性单元（一阶单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
二阶单元（带中间节点）
可表现曲线变形。应力精度大幅提高，但自由度约增加2~3倍。推荐：应力评估很重要时使用。
完全积分 vs 减缩积分
完全积分：有过约束（锁定）风险。减缩积分：沙漏模式（零能量模式）风险。根据情况选择。
自适应网格
基于误差指标（如ZZ估计量）的自动细化。有效提高应力集中区域的精度。有h法（单元细分）和p法（增加阶次）。
牛顿-拉夫逊法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿-拉夫逊法
切线刚度矩阵使用初始值或每隔几次迭代更新。每次迭代成本低，但收敛速度为线性。
收敛判定准则
力残差范数: $||R|| / ||F_{ext}|| < \epsilon$（通常 $\epsilon = 10^{-3}$〜$10^{-6}$）。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$
载荷增量法
不一次性施加全部载荷，而是分小步增加。弧长法（Riks法）可以越过载荷-位移关系的极值点进行追踪。
直接法 vs 迭代法的比喻
直接法是“用笔算精确解联立方程组”的方法——可靠但大规模问题耗时过长。迭代法是“反复猜测逼近正确答案”的方法——最初答案粗糙，但每次迭代精度都会提高。就像查字典时，从第一页开始顺序查找（直接法）不如先估计位置翻开，再前后调整（迭代法）来得高效，原理相同。
网格阶次与精度的关系
一阶单元是“用直尺近似曲线”——用直线折线表现，精度有限。二阶单元是“用柔性尺……”
相关主题
なった
詳しく
報告