This paper presents a mathematical model for simulating Tempcore process, the model consisting of three successive parts. Thermal model is used to predict the temperature distribution across the bar over its whole quenching and self-tempering rout. The output of this model is used to calculate the area of martensite formed in outer layer and ferrite–pearlite heart. The composite rule of mixtures (ROM) is used to correlate the tensile yield strength of the bar with the martensite and ferrite–pearlite volume fraction. The model is capable of predicting the effect of various process parameters namely; bar diameter, rolling temperature, number and setting of cooling nozzles, cooling water flow rate and quenching time on the thermal conditions across the bar during different stages of the process, structure formed and tensile behavior of the rebar. A good agreement is found between the predicted results and the experimental data for B500B steel with diameters from 10 to 16mm and yield strength varied from 400MPa to 700MPa.