Aiming at the defects of numerical simulation method for stability of surrounding rock with various discontinuity surfaces,the accuracy of FEM for contact problem to simulate discontinuities are pointed out.The theory of FEM for solving large-scale contact problem and its program implementation are represented and applied to simulate underground chamber excavation with discontinuity surfaces.The results show that the proposed method is precise and effective for simulating the stability of surrounding rock with discontinuities.It provides reference for engineering application.