请问为什么我自己写的有限元程序不正确？

xyf2001721

我是按照课本上所提供的有限元单元矩阵组成刚度矩阵，除去边界元素，然后对剩余刚度矩阵求拟后与力向量相乘，不但结果错误，而且正负有时候也不正确，是我组成刚度矩阵的错误还是？
谢谢！！！