When contact problems are solved by numerical approaches, a surface profile is usually described by a series of discrete nodes with the same intervals along a coordinate axis. Contact computation based on roughness datum mesh may be time consuming. An adaptive-surface elasto-plastic asperity contact model is presented in this paper. Such a model is developed in order to reduce the computing time by removing the surface nodes that have little influence on the contact behavior of rough surfaces. The nodes to be removed are determined by a prescribed threshold. The adaptive-surface asperity contact model is solved by means of the element-free Galerkin-finite element coupling method because of its flexibility in domain discretization and versatility in node arrangements. The effects of different thresholds on contact pressure distribution, real contact area, and elasto-plastic stress fields in contacting bodies are investigated and discussed. The results show that this model can help reduce about 48% computational time when the relative errors are about 5%.