The paper presents the result of modeling the thermal deformation behavior of the bearing system of a double-sided face grinder. Two modes of its operation is considered: idling and working strokes. Spindle bearings, heat fluxes from coolant, and heat fluxes from motors are considered as the main heat sources. The preliminary values of heat and convective fluxes are assigned according to well-known methods used in engineering calculations of machine tools. Refinement of the values of heat and convective fluxes is carried out according to the results of experimental values of temperatures and temperature displacements. In the thermal model, convective heat transfer is assigned to 450 surfaces. Heat flows are assigned to 161 surfaces. The simulation is run in Ansys for 6 hours of machine tool operation. Ansys uses a ten-node Solid227 element as a typical finite element. The problem of thermoelasticity is solved in a related setting. The error in modeling thermal processes did not exceed 0,5 degree, the error in modeling temperature displacements did not exceed 5 microns. The results of computer simulation confirmed the relative experimental positions of the left and right grinding wheels characteristic of this range of machines. At idle, the relative position of the grinding wheels represented the state «wider at the bottom», and during the working stroke, the state «at the bottom narrower» is recorded. The illustrations of the FEM model and contour plots for the temperature field and temperature displacements of the bearing system of the machine tool for two modes of its operation are presented.