本发明公开了一种地下工程稳定性宏细观分析方法，包括：基于隧道施工连续围岩模型进行有限差分分析，获得核心区域和非核心区域，并拆分获得非核心区域连续网格模型；对非核心区域连续模型的边界面进行抓取，并创建耦合墙体；约束隧道施工连续围岩模型的墙体节点和单元节点，在核心区域内对离散元颗粒进行平衡，获得耦合墙体的接触模型；释放所有节点的约束，对模型施加法向约束，计算模型的受力和变形，平衡模型的地应力；根据模型中的几何文件对待分析区域的单元和颗粒清零，得到隧道开挖后的模型，提取对应的位移和应力云图对地下工程稳定性进行分析。本发明这种耦合方法可以更好地模拟岩石的非均匀性特点，展现围岩的宏细观损伤。