还原一篇论文有限元程序，我这边有部分代码，计算切线刚度矩阵和内力，就下图相关矩阵，预算可谈。