1.直接法刚度矩阵建立

问题描述抽象为下图所示：将模型简化为五个弹簧串联起来之间的运动，由上到下分别有5个节点。

由基本的应力定义有：

$$
\sigma=\dfrac{F}{A} \tag{1}
$$

由基本的应变定义有：

$$
\varepsilon=\dfrac{\Delta l }{l} \tag{2}
$$

又有基本的胡克定律有：

$$
\sigma= E\cdot\varepsilon \tag{3}
$$

将上式(1)(2)同时带入式子(3)得到下式：

$$
\dfrac{F}{A}=E\cdot\dfrac{\Delta l }{l} \quad{\bf{\Longleftrightarrow}}\quad F={ {EA}\over{l}}\Delta l \tag{4}
$$

其中K={\dfrac{EA}{l}}，称为等价刚度。这个式子告诉我们了力 F 与 \Delta l 之间的关系，这就为接下来的求解各个节点的作用力与反作用力奠定了基础。

接下来我们将目光集中在弹簧的各个节点当中，用静力平衡条件来找出各个节点之间的力、位移关系。

记墙壁与弹簧的作用反力为 R 。规定方向向上为正，方向向下为负。每个单元的等价刚度记为 K_i ，每个节点的位移记为 u_i 。

则在节点 1 处有以