We propose here a model and a numerical scheme to compute the motion of rigid particles interacting through the lubrication force. In the case of a particle approaching a plane, we propose an algorithm and prove its convergence towards the solutions to the gluey particle model described in [B. Maury, ESAIM: Proceedings18 (2007) 133-142]. We propose a multi-particle version of this gluey model which is based on the projection of the velocities onto a set of admissible velocities. Then, we describe a multi-particle algorithm for the simulation of such systems and present numerical results.