Equations of motion of constrained dynamic systems for contact-impact analysis are derived from variational equations and implemented for numerical analyses of two-dimensional problems using the full nodal degrees of freedom in finite elements. Viscous damping and friction effects are also included. The method is then applied for a detailed study of several newly defined coefficients of restitution in search of a suitable one for impact analysis of deformable bodies.