A numerical analysis of rolling contact problems using quasistatic variational formulation