一维杆单元：从离散到连续的桥梁
各位同学，今天我们来聊聊有限元分析中最基础、也最经典的内容——一维杆单元。说实话，我刚开始学有限元那会儿，觉得这东西太简单了，不就是一根杆子嘛。后来在实际项目中才发现，越是基础的东西，越容易在细节上翻车。咱们今天就把这块地基打扎实了。
1. 单元位移模式：用数学描述变形
先问大家一个问题：一根杆子受力后，它内部的位移是怎么分布的？
我们假设杆单元有两个节点，每个节点有一个自由度——轴向位移。单元内部任意一点的位移，我们用一个多项式来近似。对于一维杆单元，最常用的就是线性位移模式：
u(x) = a₀ + a₁·x这里x是局部坐标，a₀和a₁是待定系数。为什么选线性？说白了，就是简单。两个节点，两个未知数，刚好能解出来。我在做桥梁索力分析时，用线性单元就能满足工程精度要求，没必要搞太复杂。
把节点条件代进去：
u(0) = u₁
u(L) = u₂解出a₀和a₁，就能得到用节点位移表示的位移场：
u(x) = (1 - x/L)·u₁ + (x/L)·u₂2. 形函数：插值的关键
上面这个表达式，其实就是形函数的雏形。我们把系数提取出来：
N₁(x) = 1 - x/L
N₂(x) = x/L这两个函数就是形函数。它们有个很有意思的性质：在自身节点处值为1，在另一个节点处值为0。你想想看，这保证了什么？保证了单元之间的位移连续性。
形函数的两个核心性质：
- 归一性：N₁ + N₂ = 1（保证刚体位移）
- 插值性：Nᵢ(xⱼ) = δᵢⱼ（Kronecker delta）
我记得有一次做复合材料杆件的分析，因为形函数选错了阶次，结果算出来的应力分布完全不对。嗯，这里要注意：形函数的阶次决定了单元的精度，但也不是越高越好，高阶单元容易产生振荡。
3. 单元刚度矩阵推导：从能量到矩阵
推导刚度矩阵，我习惯用虚功原理，思路清晰。对于一维杆单元，应变-位移关系是：
ε = du/dx = B·uᵉ其中B矩阵是：
B = [-1/L, 1/L]应力-应变关系（胡克定律）：
σ = E·ε单元刚度矩阵的通用公式是：
kᵉ = ∫ Bᵀ·E·B·A·dx对于等截面杆，积分后得到：
kᵉ = (EA/L) · [1, -1; -1, 1]这个结果很漂亮，对吧？EA/L就是轴向刚度，物理意义非常明确。我曾经在项目里用手算验证过这个矩阵，算出来的节点位移和实验值误差在3%以内，心里踏实多了。
小技巧：推导时注意积分区间是从0到L。如果单元坐标系方向反了，刚度矩阵的符号会乱掉。我建议每次推导完都检查一下：刚度矩阵应该是对称的，而且每行元素之和为零（刚体位移条件）。
4. 整体刚度矩阵组装：化零为整
单个单元的刚度矩阵有了，怎么把它们拼成整体？这个过程叫组装。说白了，就是把各个单元的贡献，按照节点编号"对号入座"。
举个例子，两个杆单元，三个节点：
单元1：节点1-2，刚度矩阵 k₁
单元2：节点2-3，刚度矩阵 k₂组装后的整体刚度矩阵是3×3的：
K = [k₁₁, k₁₂, 0 ]
[k₂₁, k₂₂+k₁₁, k₁₂]
[0, k₂₁, k₂₂]注意看，节点2对应的位置是两个单元刚度矩阵的叠加。这就是所谓的"直接刚度法"。我在做多层建筑框架分析时，就是靠这个原理，把几百根杆件组装成整体刚度矩阵。虽然矩阵规模大了，但逻辑完全一样。
避坑指南：我曾经因为节点编号搞混，导致整体刚度矩阵奇异，算出来的位移全是NaN。后来养成了习惯：组装前先画好节点编号图，用不同颜色标出每个单元覆盖的节点。这个习惯救了我好几次。
知识体系总览
下面这张图，是我自己总结的一维杆单元核心逻辑。从位移模式到形函数，再到单元刚度矩阵，最后组装成整体。每一步都有清晰的物理意义和数学表达。
好了，一维杆单元的内容就讲到这里。从位移模式到形函数，再到刚度矩阵推导和整体组装，每一步都有它的物理意义和数学逻辑。我个人觉得，理解了这个基础，后面学二维、三维单元就会轻松很多。记住：有限元分析不是黑箱，每个矩阵、每个系数背后都有明确的力学含义。
公众号:蓝海资料掘金营，微信deep3321