文档介绍：: .
类似于光弹性法或云纹法的实验结果。
对于特定的程序，其运行方法通常会在软件的附件中详述，那些较为昂贵的程序的销售
商经常会提供实****场所或举办培训班，以帮助用户弄清楚程序运行中错综复杂的情况。但用
户即使经过这种培训后，仍可能遇到问题：程序往往是“黑匣子”，其内部的运行机理是难
以明了的。在本模块中，我们将略述大部分当今使用的有限元应力分析程序的基本原理，但
目前的讨论仅限于线弹性范围内的分析。掌握这一理论有助于消除黑匣子综合症，同时也是
对固体力学的基本分析方法作一总结。

桁架的矩阵分析法
在模块 5 中，比较充分地讨论过铰接的桁架，这为引入有限元分析的概念提供了很好的
途径。对桁架精确地进行静力学分析，即使是复杂的桁架，其方程也能集成为矩阵形式，以
便于用数值方法求解。此法有时称为“矩阵分析法”，它为早期有限元分析法的发展奠定了
基础。
桁架的矩阵分析法逐一考虑每根杆件的刚度，然后用这些刚度来确定由于节点（在有限
元分析中通常称为结点）位移而产生的杆件力。注意到每根杆件对结点作用力的矢量和必须
等于作用在结点处的外力，我们可以列出一系列线性代数方程。在这些方程中，结点位移是
未知量而作用在结点上的力是已知量。这些方程可以方便地写成如下的矩阵形式（该方法亦
由此得名）：

式中， u 和 f 分别表示第 i 个结点的位移和作用在第 j 个节点上的力（实际上这些量都是
i j
矢量，其分量沿各坐标轴）。系数 K 数组称为整体刚度矩阵，其第ij 分量的物理意义是第 j
ij
个位移对第i 个力的影响。该矩阵方程可简写为
式中，用下标或黑体字表示矢量和矩阵。
开始时，每个结点所受的外力或结点的位移这两者中有一个是已知的，但对一个给定
的结点，同时指定任意的位移和任意的外力是不可能的。这些预先给出的结点力和结点位移
是问题的边界条件。矩阵分析的任务就是要确定与给定的位移相对应的力、以及结点处有已
知外力作用时的位移。

单根桁架杆件的刚度矩阵

为了得出描述桁架系统的矩阵方程，第一步是需要建立单根桁架杆件每一端的力与位
2移之间的关系。考虑如图 1 所示在 x − y 平面内的杆件，该杆件连接第 i 个和第 j 个结点，
与水平线的夹角为θ 。

图 1 单根桁架杆件

沿杆件轴向及垂直于轴线的方向研究杆件待求的伸长量δ ，δ 可用杆件两端点的位移
之差来表示：

式中，u 和 v 分别为位移的水平及垂直分量（图 1 中画出的第i 个结点的两个位移分量都是
负的）。上式可写成矩阵形式：

式中， c = cosθ 、 s = sinθ 。

图 2 结点力的分量

对线弹性体，与伸长量δ 对应的轴向力 P 可