地下工程围岩稳定性研究无论是在理论研究还是在实际应用等方面都取得了重大进展。但由于研究对象的特殊性和复杂性，目前还不能十分有效地防治重大工程安全事故的发生。因此，有必要对其进一步深入研究。本文在对国内外围岩稳定性评价研究现状进行综述的基础上，扼要阐述了三维结构弹塑性有限元的基本原理，介绍了ANSYS程序结构及其弹塑性分析的基本过程。以某矿采空区围岩稳定性评价为例，建立了矿体开采的三维数值模拟模型，并使用有限元程序ANSYS 5.6版本对矿体深部开采过程中的地压活动规律进行了三维弹塑性模拟。通过对计算所得的围岩应力分布、塑性区范围和+230m水平钢筋混凝土隔离层单元变形情况进行分析，评价了采空区围岩稳定性。研究成果对地下工程结构稳定性分析具有指导作用，并为相关领域的围岩稳定性评价工作提供参考。主要结论如下：(1) 利用大型通用有限元程序ANSYS模拟地下岩土介质的开挖（采）过程是方便可行的。ANSYS强大的功能可以建立大型而复杂的三维模型和快捷地实现有限元网格的自动划分、能够确保计算的准确性以及方便地读取分析结果。(2) 有限元计算结果表明，在充填的情况下，采空区围岩是安全的，但存在局部失稳的可能。(3) 有限元计算结果表明，三维计算模型边界的取值是合适的。在水平面范围内受矿体开采影响较大的区域不超过150m。(4) 矿床开采充填过程的ANSYS模拟结果表明，随着开采深度的增加，在+230m水平钢筋混凝土隔离层上的应力也进一步增大；当矿体开采充填到+200m水平的时候，钢筋混凝土隔离层上盘出现了大范围的塑性区。(5) 有限元计算结果显示，随着开采深度的增加，在采空区围岩中应力的集中程度也进一步增加。应力集中程度较大的地方主要位于矿体的西北和东南端，预计在该部分采场地压显现比较剧烈。(6) 数值模拟结果表明，在采用充填法充填采空区的情况下，+200m水平以下的开采活动对+230m水平钢筋混凝土隔离层的影响不大，塑性区及应力的增加程度均不明显。