程序可以解决弹性力学中有限元法里平面常体力常应变单元形势下体系整体刚度阵的生成，具体的单元编码方式可以参考徐芝纶编写的《弹性力学》，思路与结构力学中矩阵位移法的想法一致，如遇问题可以评论给楼主，楼主再附上计算书。