Abstract: A mathematical model was proposed to calculate the stress distribution of the solidification shell in a slab mold based on the local heat flux calculated by using the measured temperature of the copper plate. Solidification equations to liquid steel were built and discreted by the finite element method. ANSYS commercial software was used to solve these discrete equations, the stress and strain distributions were obtained and the taper of the mold was decided according to these calculating results.