The comprehensive understanding of freeze-thaw (F-T) deterioration and deformation properties of rock is a critical basis for frost damage prevention on cold regional rock engineering. In freezing process, plastic strain generates in rock when the stress induced by pore ice pressure (PIP) becomes higher than yield strength, and unrecoverable deformation accumulates in rock under multiple F-T cycles. Therefore, this study presents a thermal-mechanical modelling on cumulative F-T deformation of rock under cyclic F-T action. In the modelling of mechanical action, the microporomechanics theory is applied with a spherical element model in the derivation of equilibrium equations of rock, and the Drucker-Prager criterion is utilized with dilatancy effect of rock considered to describe unrecoverable plastic deformation cumulated in rock under multiple F-T cycles. The thermal process is modelled with unfrozen water content and latent heat accounted. The proposed thermal-mechanical modelling method is validated based on experiments and numerical simulations of F-T deformation of sandstone, and it is able to simulate the temperature and F-T deformation of rock in acceptable precision. The numerical distributions of temperature, PIP, maximum principal strain and displacement are exhibited and show that the plastic region begins to develop in rock matrix with PIP exceeding the elastoplastic critical pressure when temperature drops to about −2 °C and grows rapidly when temperature drops from −2 to −6 °C. Several F-T cycles later, large unrecoverable plastic strain remains in sandstone though PIP dissipates after thawing. The cumulative F-T deformation, porosity increment and maximum of plastic volume ratio are much greater under condition with freezing temperature of −20 °C, compared with condition −10 °C. Besides, as cycle number increases, the plastic volume ratio becomes smaller and its decrease rate attenuates.