Abstract. This work analyses a conserving time integration scheme applied to impact dynamics of one or more points against a rigid surface. This rigid surface may represents the boundary of a rigid body or it could represent a rigid cavity inside a body, and it should be represented by an implicit function. The points may be material particles, vertex of a rigid polyhedra, or nodes from a finit...