单元刚度矩阵（Element Stiffness Matrix）是结构分析中有限元方法（Finite Element Method, FEM）中的一个核心概念。它用于描述结构中单个单元的力学行为，特别是在受到外力作用时的变形和应力分布。以下是关于单元刚度矩阵表达式的详细文档：
一、引言
在有限元方法中，复杂结构被离散化为一系列相互连接的单元。每个单元内部的位移场通过假设的位移模式（如多项式函数）来近似表示。基于这些位移模式和材料的本构关系（如弹性模量、泊松比等），可以推导出单元刚度矩阵。
二、基本概念
- 节点：构成单元的端点，通常具有已知的位移或力边界条件。
- 位移模式：用于描述单元内部任意点位移的函数，通常是节点位移的线性或非线性组合。
- 几何方程：连接应变与位移的关系式。
- 物理方程：连接应力与应变的关系式（如胡克定律）。
- 虚功原理：用于推导单元刚度矩阵的基本原理之一。
三、单元刚度矩阵的推导
以二维平面问题中的线弹性材料为例，假设一个四边形单元（或其他类型的单元），其节点位移向量为{u_e}，包括所有节点的水平和垂直位移。
选择位移模式：例如，对于常应变三角形单元，位移模式可以表示为节点位移的线性插值。
计算应变：利用几何方程，将位移模式代入得到应变向量{\epsilon}。
计算应力：根据物理方程（如胡克定律），将应变向量代入得到应力向量{\sigma}。
应用虚功原理：考虑一个虚拟的外力作用于单元上，使得单元产生虚拟位移{\delta u_e}。根据虚功原理，外力做的虚功等于内力做的虚功。即： [ \int_{\Omega} {\sigma} : {\delta \epsilon} , d\Omega = {F_{ext}}^T {\delta u_e} ] 其中，{\delta \epsilon}为虚拟应变，\{F_{ext}\}为外力向量。
推导单元刚度矩阵：通过一系列数学变换，可以将上述等式转化为标准的线性方程组形式： [ {K_e} {\Delta u_e} = {F_{int}} ] 其中，\{K_e\}为单元刚度矩阵，\{\Delta u_e\}为节点位移增量，\{F_{int}\}为内力向量。
组装整体刚度矩阵：将所有单元的刚度矩阵按照节点的连接关系组装成结构的整体刚度矩阵。
四、单元刚度矩阵的具体表达式
对于不同类型的单元（如梁单元、平面应力/应变单元、三维实体单元等），单元刚度矩阵的具体表达式会有所不同。以下是一个简单的二维平面应力问题的常应变三角形单元的刚度矩阵示例：
[ {K_e} = \frac{E}{2(1-\nu^2)} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \ 0 & c_1 & 0 & c_2 & 0 & c_3 \ b_2 & 0 & b_4 & 0 & b_5 & 0 \ 0 & c_2 & 0 & c_4 & 0 & c_5 \ b_3 & 0 & b_5 & 0 & b_6 & 0 \ 0 & c_3 & 0 & c_5 & 0 & c_6 \ \end{bmatrix} ]
其中，E为弹性模量，\nu为泊松比，b_i和c_i是与单元几何形状和尺寸相关的系数。
五、结论
单元刚度矩阵是有限元方法中描述单元力学行为的关键工具。通过选择合适的位移模式、应用几何方程和物理方程以及虚功原理等基本原理，可以推导出适用于不同类型单元的刚度矩阵表达式。这些刚度矩阵进一步组装成整体刚度矩阵后，可用于求解整个结构的位移和应力分布。
