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

中南大学有限单元法考试试卷(岳乐-庞俊)。 有限元分析复习点点滴滴——张义涵 1402 1、 什么是平面应力问题？什么是平面应变问题？（张毅涵做） 答：平面应变问题是指薄板受平行于板面且沿厚度均度载荷用 。 在考虑有限元法表述位移时，首先是选取一组函数，他们可以用节点位移来表示有限元内任一点的位移分量。然后从外加位移场发展解法得的各个步都是，应变分量由位移的各种导数唯一确定，于是外加位移变化确定了整个单元的应变状态。这些诱导应变和介质弹性性质一起确定了单元的诱导应力。将初始应力与诱导应力叠加就

地下工程中硐室围岩常为破碎岩体,在合理的支护体系下仍具有一定的承载力,针对这一现象,开展了不同约束机制下岩石的真三轴力学试验,研究约束体系强度和刚度对岩石峰后继续破碎过程中应力应变的影响,并将这些与约束体系有关的应力应变关系绘制在一个坐标系中,再从经济性、安全性、可靠性和高效性等方面出发,讨论了考虑以上多种因素的本构关系。这样实际工程支护体系设计中,只要应力应变设计值落在合理的区域内即可,避免了本构关系复杂而导致不便于推广的问题,更主要的是峰后的本构关系与实际工程围岩介质为岩体相一致,这为岩石力