A physically-based impact model – already known and exploited in the field of sound synthesis – is studied using both analytical tools and numerical simulations. It is shown that, for some regions of the parameter space, the trajectories of discretized systems may drift from analytically-derived curves. Some methods, based on enforcing numerical energy consistency, are suggested to improve the accuracy and stability of discrete-time systems.