Variational and numerical analysis of a quasistatic viscoelastic problem with normal compliance, friction and damage