材料力学求解器-刚架与桁架杆系的计算机求解（附matlab代码）

材料力学是一门非常成熟的学科，里面有大量的已经成熟的计算方法、计算模型。然而在掌握了基本原理之后，求解问题的主要时间都花在了计算上面。这时，则完全可以利用计算机采用数值方法进行快速求解。

本文的参考文献很简单，就是 单辉祖 编写的 材料力学第3版下册（或者说是第II册），第18章的内容。

还是惯例声明：本人没有相关的工程应用经验，只是纯粹对相关算法感兴趣才写此博客。所以如果有错误，欢迎在评论区指正，不胜感激。本文主要关注于算法的实现，对于实际应用等问题本人没有任何经验，所以也不再涉及。

1 刚架的计算机求解

1.1位移法与刚度矩阵

计算机求解刚架的原理主要是位移法和刚度矩阵。
基本思路为：

1 定义划分单元与节点
2 列出每个节点的载荷F
3 求解每个单元的刚度矩阵Km
4 把所有单元的刚度矩阵相加，得到总的刚度矩阵K
5 利用位移法与胡克定律F=K*D，求解出位移D
6 根据位移，反推出其余受力

以书上的例18-4，举例说明

1 定义划分单元与节点

上面的刚架，可以很轻易的划分出单元与节点，包含4个节点与3和单元。此外，刚度矩阵的方法要求把受力点算作节点。

2 列出每个节点的载荷F

对于刚架来说，每个节点受力包含Fx、Fy、Me三个力。所以F是一个3*4=12个元素的列向量。
外加的力载荷只有2号节点处的三个力。所以，列向量F可以写作：
$${\left[\begin{array}{cccccccccccc}0& 0& 0& Fx& -Fy& -Me& 0& 0& 0& 0& 0& 0\end{array}\right]}^{\prime }$$
这里’代表转置。
外加力的方向定义为与x、y轴方向相同为正。力矩定义为逆时针为正。

3 求解每个单元的刚度矩阵Km
首先是局部坐标系下，列出刚度矩阵
之后，以③号单元为例，3号单元的两个节点分别为i=2和j=4。
单独某个单元的刚度矩阵Km如下：
[ E A L 0 0 − E A L 0 0 0 12 E I L 3 6 E I L 2 0 − 12 E I L 3 6 E I L 2 0 6 E I L 2 4 E I L 0 − 6 E I L 2 2 E I L − E A L 0 0 E A L 0 0 0 − 12 E I L 3 − 6 E I L 2 0 12 E I L 3 − 6 E I L 2 0 6 E I L 2 2 E I L 0 − 6 E I L 2 4 E I L ] \left[ \begin{matrix} \frac{EA}{L}& 0& 0& -\frac{EA}{L}& 0& 0\\ 0& \frac{12EI}{L^3}& \frac{6EI}{L^2}& 0& -\frac{12EI}{L^3}& \frac{6EI}{L^2}\\ 0& \frac{6EI}{L^2}& \frac{4EI}{L}& 0& -\frac{6EI}{L^2}& \frac{2EI}{L}\\ -\frac{EA}{L}& 0& 0& \frac{EA}{L}& 0& 0\\ 0& -\frac{12EI}{L^3}& -\frac{6EI}{L^2}& 0& \frac{12EI}{L^3}& -\frac{6EI}{L^2}\\ 0& \frac{6EI}{L^2}& \frac{2EI}{L}& 0& -\frac{6EI}{L^2}& \frac{4EI}{L}\\ \end{matrix} \right] ​LEA​00−LEA​00​0L312EI​L26EI​0−L312EI​L26EI​​0L26EI​L4EI​0−L26EI​L2EI​​−LEA​00LEA​00​0−L312EI​−L26EI​0L312EI​−L26EI​​