It has been observed that, under cyclic thermal loading conditions, the as-solidified microstructures of solder interconnections composed of a few large colonies are transformed into more or less equi-axed grains by recrystallization. These recrystallized microstructures provide continuous networks of grain boundaries through solder interconnections, and, consequently, they offer favorable paths for cracks to propagate intergranularly. In this work a quantitative multiscale method, combining the Monte Carlo (MC) and finite element (FE) methods, was developed for the first time to predict the recrystallization and grain growth in solder interconnections. This method can be used to establish a lifetime prediction model that is based on these microstructural changes. The inhomogeneous nucleation model was employed in the MC model, while the non-uniform stored energy distribution was scaled from the FE model results. The assumption that, during thermal cycling 5% of the work in the interconnections is stored as energy, gives results that are in agreement with experimental observations. The method predicts well the incubation period and the growth rate of the recrystallization as well as the expansion of the recrystallized region.