Turbomachinery rotor–stator unilateral contact induced interactions play a growing role in lifecycle analysis and thus motivate the use of accurate numerical prediction tools. Recent literature confirmed by ongoing in-house experiments have shown the importance of thermomechanical coupling effects in such interactions. However, most available (possibly reduced-order) models are restricted to the sole mechanical aspects. This work describes a reduction technique of thermomechanical models involving unilateral contact and frictional contact occurrences between rotor and stator components. The proposed methodology is grounded on Guyan and Craig–Bampton methods for the reduction of the structural dynamics in conjunction with Krylov subspace techniques, and specifically the Craig–Hale approach, for the reduction of the thermal equations. The method has the capability to drastically reduce the size of the model while preserving accuracy. It stands as a reliable strategy to perform simulations of thermomechanical models with localized mechanical and thermal loads.