Induction heating is used in many industrial applications to heat electrically conductive materials. The coupled electromagnetic-thermal induction heating process is non-linear in general, and for ferromagnetic materials it becomes challenging since both the electromagnetic and the thermal responses are non-linear. As a result of the existing non-linearities, simulating the induction heating process is a challenging task. In the present work, a coupled transient electromagnetic-thermal finite element solution strategy that is appropriate for modeling induction heating of ferromagnetic materials is presented. The solution strategy is based on the isothermal staggered split approach, where the electromagnetic problem is solved for fixed temperature fields and the thermal problem for fixed heat sources obtained from the electromagnetic solution. The modeling strategy and the implementation are validated against induction heating experiments at three heating rates. The computed temperatures, that reach above the Curie temperature, agree very well with the experimental results.