Variational multiscale method for fully coupled thermomechanical interface contact and debonding problems. (February 2021)