4.5.2 计算结构总刚度矩阵程序