A point-to-face contact model for three-dimensional numerical manifold method is established in this paper.Contact submatrices,including those for the normal spring,shear spring and the frictional forces,are derived by the penalty function method.Also given are the open-close iteration criteria and operations performed for different changes in contact state.An example of the analysis of underground opening stability is provided to demonstrate the capability of the contact model.