在我国石膏矿山普遍采用房柱式采矿方法进行开采,由于超高超宽开采或设计不合理,相继发生了采空区顶板崩落、矿柱倒塌、地表岩移等现象,农田、村庄、建筑物、交通水利设施等遭受破坏,个别矿山甚至发生了人员伤亡等极为严重的事件。本文针对苏北石膏矿大面积采空区稳定性问题,综合运用现场调研、实验室测试、理论分析及数值计算的研究方法,对石膏矿矿柱及顶板稳定性、采空区链锁式失稳反应机理以及大面积悬露顶板的稳定性等进行了系统研究,主要研究成果如下:(1)通过实验对石膏岩样进行了研究,采用电子万能试验机对岩样单轴抗压强度、单轴抗拉强度和抗剪强度进行了测试,并根据测试结果进行拟合分析,得到石膏岩体弹性模量、泊松比、粘聚力及内摩擦角等参数,为采空区稳定性理论分析以及数值模拟提供基本参数。(2)在考虑围岩破坏特征的基础上,通过建立石膏岩体损伤演化模型,引入形状参数,得出了石膏岩体改进的损伤本构方程;基于突变理论并结合石膏岩体的损伤本构方程推导出矿柱的尖点突变方程,得出矿柱的失稳判据。通过对顶板及矿柱变形和破坏特征分析,结合薄板理论,建立矿房顶板及矿柱力学模型,揭示石膏矿矿房矿柱破坏规律。(3)根据矿柱载荷和矿柱强度的相互关系,计算出矿柱的安全系数,提出了采空区链锁式反应的动态计算方法。在此基础上建立了采空区大面积悬露顶板的托板力学模型,分析矿柱失稳后形成的大面积悬露顶板的应力分布特征,根据顶板稳定性判据确定顶板的稳定性。(4)结合苏北石膏矿地质特征,对采空区和巷道的地质形态进行分类,运用Flac3D对典型区域进行模拟;运用链锁式动态计算方法和托板理论对矿柱和矿柱失稳后形成的大面积悬露顶板进行理论计算。计算结果表明,数值计算与模拟计算结果基本一致。