On numerical treatment of thermomechanical contact problems