【原创】矩阵法计算步骤

用矩阵法计算结构的全过程归纳如下：

1.结构离散化，节点、杆元编号，建立结构坐标系与杆元坐标系。

2.计算各杆元刚度矩阵[K'_e]，将其转化为[K_e]＝[T][K'_e][T]^T。
3.组装结构刚度矩阵[K]，按各杆元节点总体编号将[K_e]对号入座。

4.计算各杆元的固端力{P'_ij},将其转化为{P_ij}＝－[T]{P'_ij}。
5.组装结构外荷载向量{P}，按总体编号将{P_ij}及节点上集中外力对号入座。

6.约束处理。

7.求解结构平衡方程式，得到节点位移{D}。
8.将节点位移返回杆元坐标系{d'_ij}＝[T]^T{d_ij}，从而求出杆元的内力：
   {f'_ij}＝[K'_e]{d'_ij}＋[P'_ij]

9.若需要，可利用未经约束处理的平衡方程，求出支座反力。