Rigorous non-linear models of elasto-plastic contact deformation are time-consuming in numerical calculations for the distinct element method (DEM) and quite often unnecessary to represent the actual contact deformation of common particulate systems. In this work a simple linear elasto-plastic and adhesive contact model for spherical particles is proposed. Plastic deformation of contacts during loading and elastic unloading, accompanied by adhesion are considered, for which the pull-off force increases with plastic deformation. Considering the collision of a spherical cohesive body with a rigid flat target, the critical sticking velocity and coefficient of restitution in the proposed model are found to be very similar to those of Thornton and Ning’s model. Sensitivity analyses of the model parameters such as plastic, elastic, plastic-adhesive stiffnesses and pull-off force on work of compaction are carried out. It is found that by increasing the ratio of elastic to plastic stiffness, the plastic component of the total work increases and the elastic component decreases. By increasing the interface energy, the plastic work increases, but the elastic work does not change. The model can be used to efficiently represent the force-displacement of a wide range of particles, thus enabling fast numerical simulations of particle assemblies by the DEM.