Formation of the shrinkage defects in ductile iron castings is far more complicated phenomenon than in other casting alloys. In the presented paper changes the ductile iron density during solidification is analyzed. During the solidification path the influence of the temperature, phase fractions and phase composition is taking into account. Computer model, using cellular automata method, for estimation of changes in density of ductile iron during its solidification is applied. Results of the solidification modeling for Fe-C binary alloys with different composition in the castings with a different wall thickness are presented. As a result of calculations it was stated that after undercooling ductile iron below liquidus temperature volumetric changes proceed in three stages: pre-eutectic shrinkage (minimal in eutectic cast iron), eutectic expansion and the last shrinkage.