Spheroidal graphite (SG) cast iron is characterized by the presence of spherical graphite nodules distributed in the metallic matrix. The performance of castings is primarily dependent on the solidification microstructures. In this work, a two dimensional (2D) multi-phase cellular automaton (MCA) model previously proposed by the present authors is improved to simulate the microstructure evolution of SG cast iron during divorced eutectic solidification. The present model adopts a local solutal equilibrium approach to calculate the driving force for the growth of both graphite and austenite phases. The density difference between iron and graphite is also taken into account. The diffusion of solute in the simulation domain is calculated using a finite difference method (FDM). The present model is applied to simulate the evolution of microstructure and carbon concentration field during solidification for hypereutectic SG cast irons. The results show that the present model can reasonably describe the typical features of divorced eutectic solidification, involving the independent nucleation and growth of primary graphite and austenite dendrites in liquid, the competitive growth of adjacent graphite nodules, engulfment of graphite nodules by austenite dendrites, the isotropic growth of the austenite shells that envelop the graphite nodules, the austenite to graphite eutectic phase transition controlled by carbon diffusion through the solid austenite shell, and multiple graphite nodules encapsulated in each austenite grain at the end of eutectic solidification. The simulated volume fraction and average diameter for graphite nodules are compared reasonably well with the experimental data and level rule calculation. The interactive and competitive growth behavior between austenite dendrites and graphite nodules is studied in detail. It is found that the growth of a graphite nodule is promoted by the approaching austenite. However, after embedded by an austenite dendrite, the growth velocity of graphite decreases rapidly because of lower carbon diffusivity in austenite than that in liquid. In addition, the effect of cooling rate on the size of graphite nodules is also investigated. The results show that with cooling rate increasing, the size distribution of graphite nodules varies from two peaks to one peak, and the average diameter of nodules decreases. The simulation results compare reasonably well with the experimental data reported in literature, demonstrating the validity of the present model。