Algorithmic aspects of energy dissipation mechanisms of dynamic structural systems in conjunction with central difference time integration method are investigated and damping proportional to M(M−1K)m (where K is the stiffness matrix, M is the mass matrix and m is a damping parameter) is proposed. Detailed algorithms for M(M−1K)m proportional damping for the central difference time integration method are presented together with stability criteria and numerical test problems. © 1998 John Wiley & Sons, Ltd.