Temper rolling is a special form of flat rolling. In temper rolling, the thickness of the sheet is reduced by 0.5 to 2 percent. Its main purpose is to prevent the occurrence of the yield point phenomena during the forming of the sheet and the formation of deformation strips (luder’s lines), which, if formed, would create an inappropriate appearance on the sheet. The traditional routine rolling models do not have the ability to accurately predict the amount of force required for this type of rolling. The two main causes of error in traditional models are the failure to consider the elastic recovery area of the sheet which results in errors during the contact of the sheet and the work roll, as well as the use of the Hitchcock relationship or other similar relationships for the elastomeric deformation of the work roll to obtain the length of contact of the sheet and the work roll. Observations have shown that the elastic deformation of the work roll in the area of contact is irregular and non-circular. In this research, a numerical model was developed in which the Airy function was used to express the elastic behavior of the roller. For the sheet, in addition to the conventional slab method, the Prandtel-Rus equations were used to describe the elastic-plastic behavior of the sheet. To simulate the behavior of the sheet under applied tensions, stress-strain curve fitting was used by genetic algorithm and the tensions from the back and front were considered as boundary conditions for the sheet. Finally, for validation of the present research model, rolling data of skinpass Mobarakeh Steel Complex was collected and used and the effect of two causes of elastic recovery of sheet and elastic work roll deformation were investigated in prediction of the amount of force. Also, a comparison was made between the results of the current model and other available models. Keywords: Temper Rolling, Slab Method, Airy Function, Elastic Recovery, Elastic work roll Deformation