A discrete system constituted of particles interacting by means of a centroid-based law is numerically investigated. The elements of the system move in the plane, and the range of the interaction can be varied from a more local form (first-neighbours interaction) up to a generalized nth order interaction. The aim of the model is to reproduce the behaviour of deformable bodies with standard (Cauchy model) or generalized (second gradient) deformation energy density. The numerical results suggest that the considered discrete system can effectively reproduce the behaviour of first and second gradient continua. Moreover, a fracture algorithm is introduced and some comparison between first- and second-neighbour simulations are provided.