The work presents the results of numerical investigations of AZ31 magnesium alloy bars obtaining by the rolling process in a three high skew rolling mill. An analysis of impact the magnitude of inflicted deformation on temperature distribution in the material was made. The calculations were made using Finite Element Method (FEM) for 3D deformation state taking into account thermal phenomena that occur during applied deformation scheme. Theoretical investigations were made for two variants of deformation, one temperature of the process that equals 400 °C and four rolling rates of 25, 75 and 100 rev/min