采空区规模的不断扩大,给矿山安全生产带来了极大的威胁。应用数值模拟方法,能够系统地分析采空区周围岩体的应力分布规律和变形趋势,鉴于拉格朗日算法的FLAC 3D程序适合于求解岩土力学工程中非线性的大变形问题。以石人沟铁矿为工程背景,运用FLAC 3D对采空区的稳定性进行了数值模拟分析,并得到空区围岩的应力应变分布规律,对空区的稳定性进行了科学分析与评价,提出了应及时处理采空区的建议。