弹簧与胡克定律 — F=kx到FEM刚度矩阵[K]
1. 胡克定律：FEM的起点
我听说FEM里有个很重要的[K]矩阵，这到底是什么东西？
从本质上说，[K]就是弹簧刚度的矩阵版本。弹簧的F=kx里，k就是"刚度"——施加1N的力会产生1/k米的位移。把一个结构离散成千个单元，每个单元是一个小弹簧，把所有弹簧的k"组装"成一个大矩阵，就是[K]。结构越硬，K越大，同样的力产生的变形越小。
胡克定律是线弹性力学的基石，也是FEM刚度矩阵的物理根源。1678年，罗伯特·胡克在弹簧实验中发现了力与伸长量的线性关系，这个简单的比例关系后来成为了工程结构分析的核心工具。
2. 胡克定律的完整表述
一维胡克定律（弹簧）：
其中 $k$ 是弹簧刚度（N/m），$x$ 是伸长量（m）。
连续体一维胡克定律（杆件）：
这里杆件的等效刚度 $k_{\text{rod}} = EA/L$，其中 $E$ 是弹性模量，$A$ 是截面积，$L$ 是长度。
三维胡克定律（广义胡克定律）：
对于各向同性线弹性材料，本构矩阵 $[D]$ 用弹性模量 $E$ 和泊松比 $\nu$ 表示：
3. 弹性常数：E、G、ν的关系
对于各向同性材料，三个弹性常数中只有两个是独立的：
| 常数 | 符号 | 物理含义 | 典型量级 |
|---|---|---|---|
| 弹性模量（Young's） | $E$ | 轴向应力/轴向应变 | 钢：210 GPa，铝：70 GPa |
| 泊松比 | $\nu$ | 横向应变/轴向应变（负值） | 金属：0.25~0.35 |
| 剪切模量 | $G$ | 剪切应力/剪切应变 | 钢：80 GPa |
| 体积模量 | $K$ | 静水压力/体积应变 | 钢：170 GPa |
各弹性常数之间的关系：
泊松比是说，轴向拉伸时横向会收缩对吗？那橡胶的泊松比接近0.5是什么意思？
对。泊松比接近0.5意味着材料几乎不可压缩——拉伸时体积基本不变，轴向伸长多少，截面积就缩小多少。橡胶、生物软组织（如肌肉）都是这样。在FEM里，泊松比接近0.5时（>0.49）要用特殊的单元配方（混合单元）处理"体积锁定"问题，否则刚度矩阵会过于刚硬，计算结果不准确。
常用工程材料弹性常数对比：
| 材料 | E (GPa) | ν | G (GPa) | 密度 (kg/m³) |
|---|---|---|---|---|
| 结构钢 Q355 | 206 | 0.30 | 79 | 7850 |
| 不锈钢 316L | 193 | 0.27 | 76 | 7990 |
| 铝合金 6061-T6 | 68.9 | 0.33 | 26 | 2700 |
| 钛合金 Ti-6Al-4V | 114 | 0.34 | 43 | 4430 |
| CFRP（准各向同性） | 70 | 0.30 | 27 | 1600 |
| 混凝土 C30 | 30 | 0.20 | 12.5 | 2400 |
| 自然橡胶 | 0.001~0.01 | ~0.499 | 0.0003 | 920 |
4. 弹簧系统：串联与并联
多弹簧系统的等效刚度：
串联弹簧（各弹簧受相同的力，变形叠加）：
并联弹簧（各弹簧变形相同，力叠加）：
这两个公式的重要工程应用：
- 汽车悬挂的弹簧（主弹簧与辅助缓冲）——并联
- 螺栓连接的法兰：螺栓刚度与垫片刚度——串联
- 复合材料层合板：各层刚度——面内并联，面外串联
- 地基弹簧：Winkler弹簧模型——并联分布弹簧
5. 弹簧单元刚度矩阵推导
以最简单的轴向弹簧单元为例，两节点弹簧，节点位移为 $u_1, u_2$：
节点1的平衡：$F_1 = k(u_1 - u_2)$
节点2的平衡：$F_2 = k(u_2 - u_1) = -k(u_1 - u_2)$
写成矩阵形式：
这就是弹簧单元的刚度矩阵 $[k]^{(e)} = \begin{bmatrix} k & -k \\ -k & k \end{bmatrix}$。
注意它的特性：
- 对称性：$k_{ij} = k_{ji}$（麦克斯韦互等定理）
- 奇异性：无边界条件时行列式为零（可以整体平动，不产生应力）
- 主对角线正定：施加边界条件后，全局刚度矩阵正定可解
类比到一维杆件单元（两节点），刚度矩阵相同形式，但刚度值为 $k = EA/L$：
6. 刚度矩阵的组装
以三节点、两弹簧系统为例说明组装过程：
弹簧1（节点1-2，刚度 $k_1$）：$[k^{(1)}] = k_1\begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix}$ 对应节点1、2
弹簧2（节点2-3，刚度 $k_2$）：$[k^{(2)}] = k_2\begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix}$ 对应节点2、3
组装到全局 $3\times3$ 矩阵（直接刚度法）：
这个矩阵就是两弹簧串联系统的全局刚度矩阵。节点2的方程 $-k_1 u_1 + (k_1+k_2)u_2 - k_2 u_3 = F_2$ 是两侧弹簧力的平衡，体现了力平衡的自动满足。
为什么中间节点的对角线元素是$k_1+k_2$，而不是其他值？
因为节点2同时连接了两个弹簧！如果节点2单独移动，它需要克服两个弹簧的合力$k_1+k_2$。对角线元素代表的是"这个节点受到多大的刚性约束"。一般来说，连接越多单元的节点，对角线值越大，整体越刚。这也是为什么网格越密（等效于用更多弹簧），刚度矩阵维度增大但精度提高。
7. 工程材料刚度对比
结构刚度 $k = EA/L$ 可用于快速估算各类结构的刚度：
| 结构 | 材料 | E (GPa) | 截面 A (cm²) | 长度 L (m) | 刚度 k (MN/m) |
|---|---|---|---|---|---|
| 汽车车身纵梁 | 高强钢 | 210 | 20 | 1.5 | 280 |
| 铝合金车架型材 | 6061-T6 | 69 | 25 | 1.5 | 115 |
| 碳纤维传动轴 | CFRP | 130 | 15 | 1.0 | 195 |
| 桥梁混凝土柱 | C40 | 32.5 | 2500 | 10 | 8125 |
| 建筑钢柱 H400 | Q355 | 206 | 218.4 | 4 | 11234 |
8. 结构刚度设计原则
工程中"刚度"与"强度"是不同的设计指标：
- 强度：抵抗破坏（应力不超过强度极限）
- 刚度：抵抗变形（位移不超过使用要求）
提高结构刚度的工程手段：
| 方法 | 原理 | 典型应用 |
|---|---|---|
| 加大截面积 A | $k \propto A$ | 高层建筑核心筒加厚 |
| 缩短跨度 L | $k \propto 1/L$ | 增加中间支撑柱 |
| 选用高模量材料 | $k \propto E$ | CFRP代替铝合金（机器人臂） |
| 优化截面形状 | 提高 $I$（惯性矩） | 工字钢（弯曲刚度大幅提升） |
| 增加预应力 | 提高几何刚度 | 预应力混凝土桥梁 |
悬臂梁末端刚度（弯曲刚度）：
其中 $I$ 是截面惯性矩。对于矩形截面 $b \times h$：$I = bh^3/12$。这表明梁的高度对弯曲刚度影响最为显著（三次方关系），这就是工字梁、箱梁等高效截面形式的设计依据。
9. 非线性弹簧与材料非线性
胡克定律只在弹性范围内成立。超过弹性极限后，需要考虑材料非线性：
常见的非线性本构模型：
- 双线性弹塑性：最简单的弹塑性模型，屈服后切线模量 $E_T < E$
- Ramberg-Osgood：幂次硬化，适用于铝合金等渐进屈服材料
- Mooney-Rivlin / Neo-Hookean：超弹性，适用于橡胶大变形
- Drucker-Prager：土、混凝土等压力相关屈服
在FEM中，非线性材料需要迭代求解（Newton-Raphson法）：每一步用切线刚度矩阵 $[K_T]$ 代替线弹性刚度矩阵 $[K]$：
其中 $[D_T]$ 是切线本构矩阵（对应当前应变状态下的瞬时刚度）。
总结
胡克定律 $F = kx$ 是FEM刚度矩阵 $[K]\{u\} = \{F\}$ 的物理基础。从一维弹簧单元出发，通过刚度矩阵的推导和组装，我们看到了如何将连续体的弹性响应离散化为代数方程。当材料进入非线性范围，切线刚度矩阵取代常数刚度矩阵，但物理概念保持不变。