单元刚度矩阵：有限元如何“化整为零”？ 原创

有限元分析（FEA）的核心思想是什么？

八个字——化整为零，积零为整。

把一个复杂的结构，先拆解成一个个微小的“单元”，分析每个单元的力学行为，再组装回整体求解。

这一切的起点，就是单元刚度矩阵。

今天，我们以最经典的桁架杆单元为例，带你彻底搞懂它。

一、从“离散化”说起

想象一根钢杆，两端受拉。在经典力学中，我们直接用胡克定律就能算出伸长量。

但如果是一个由几十根杆件组成的复杂桁架桥呢？每根杆的受力、变形都不一样，相互之间还有约束。

这时候，就需要有限元法登场了。

有限元分析的核心思想是“离散化”——将连续的结构体划分为有限个、互不重叠的单元，单元之间通过节点连接。

每个单元内部的力学行为，通过单元刚度矩阵来描述。单元刚度矩阵，就是描述单元变形与受力关系的“核心公式”。

二、杆单元刚度矩阵的推导，三步走

以最简单的二节点杆单元为例，它只承受轴向拉力或压力，每个节点有1个自由度（沿杆长方向的位移）。

第一步：形函数——用节点位移表达内部任意点位移

假设杆内位移沿杆长方向线性变化。

通过形函数，我们可以用两个端点的节点位移，插值出杆内任意位置的位移。

形函数的作用：用节点位移“插值”出单元内部任意位置的位移。这是有限元法的核心技巧。

第二步：应变矩阵——用节点位移表达应变

通过几何方程求导，我们可以得到一个应变矩阵 B，它直接将节点位移映射为单元内部的应变。

应变矩阵 B的作用：将节点位移转换为单元应变，是连接几何与力学的桥梁。

第三步：最小势能原理——推导出单元刚度矩阵

通过最小势能原理，最终推导出杆单元的刚度矩阵：

这个矩阵，就是描述杆单元变形与受力关系的“核心公式”。

三、矩阵的物理意义，一句话说清

这个小小的2X2矩阵，蕴藏着丰富的物理意义：

• 对角线元素为正值:当一端产生位移时，同端需要施加同向的力——拉就拉，推就推。
• 非对角线元素为负值:当一端产生位移时，另一端会受到反向的力——你拉一端，另一端必然受到反作用力，这就是力的平衡。
• 矩阵对称：这是线性弹性力学的基本性质，也是有限元法能够高效求解的基础。

【总结】

一个简单的桁架，由多个杆单元组成。

如何将每个单元数值组装成整个结构的总刚度矩阵K？

这背后涉及两个关键步骤：

1. 坐标变换：每个单元有自己的局部坐标系（沿杆长方向），而整个结构建立在全局坐标系下。
2. 需要将单元刚度矩阵从局部坐标变换到全局坐标。
3. 矩阵组装：根据单元节点编号与全局节点编号的对应关系，将单元刚度矩阵的各个元素“对号入座”，叠加到总刚度矩阵的对应位置。

总刚度矩阵的特点：大型、稀疏、对称矩阵。

之所以"稀疏"，是因为每个节点只与少数几个相邻节点相连，大部分元素为0。

这也是为什么有限元可以高效求解大规模问题的原因——软件只需存储和计算非零元素，大大节省了计算资源。

理解了这些，你就能更好地诊断仿真问题——比如当出现“零刚度”报错时，往往是某个单元没有正确连接到结构中，导致总刚度矩阵奇异。

🧰单元刚度矩阵是有限元分析的基石。

从它出发，你可以逐步理解：

• 如何将单元刚度矩阵组装成整体刚度矩阵
• 如何处理边界条件（约束、载荷）
• 如何求解节点位移，进而得到每个单元的应力和应变

扫码添加客服，带你从单元刚度矩阵开始，一步步搭建完整的有限元分析知识体系！