本文研究了中低围压下弹塑性岩石力学性质的应力空间和应变空间内的本构描述,并且编写和调试了关于井壁稳定力学分析的程序,对两个算例进行了稳定性的分析,得到了井壁弹塑性软化阶段的平衡路径曲线和位移、应力、应变云图。在应力空间和应变空间内推导岩石弹塑性变形的本构关系,讨论了新的屈服准则,给出标准方程,数值计算方法表达式。讨论了材料有限元分析数值计算方法,包括乘大数法、应力增量的计算、求点在屈服面上的投影。进行有限元方法基础程序调试,以实验测试数据为初值,进行应力、应变等变化计算分析。编制和调试程序的过程包括Visual C++界面中的配置,Visual C++调用Matlab引擎的配置和用Patran查看结果。程序的内容包括主程序和在单元节点、积分点、屈服面上的相关计算,迭代过程所需要的数学计算,还有产生Patran文件的子程序。程序的结果用平衡路径曲线表示井筒的塑性软化过程。得到了平衡路径曲线和应力、应变和位移云图,与Abaqus计算结果进行对比,通过结果的对比分析证明程序得到了更好的结果。