Abstract: Based on taking D-P criterion as the yield function to establish the elastic-plastic constitutive model, numerical simulation of roadway surrounding rock stability was carried on by using COMSOL Multiphysics numerical simulation software, stress and displacement distribution of roadway surrounding rock was analyzed, and rock plastic zone distribution with yielding function of the Mohr-Coulomb criterion and D-P criterion was contrasted, then the failure form of roadway was discussed. Research results show that taking Mohr-Coulomb criterion as the yield function may underestimate the roadway surrounding rock plastic zone distribution, tensile failure first occurs in the two-side roadway, then occurs shear failure, eventually forms a X shape conjugate shear failure.