A self-consistent model is developed to describe the formation of misfit dislocation loops in hemispherical epilayers of finite size grown on thick substrates. The lattice mismatch between the substrate and the epilayer is described by virtual interfacial dislocation loops situated in the interface The free surface boundary conditions are satisfied by the surface dislocation loops situated on the surface of the hemispherical island. A misfit dislocation loop is nucleated and the changes in the energy of the configuration are used to determine if a lowering in the total energy of the configuration can occur. The model consisting of the prismatic dislocation loops in the epilayer is replaced by a simpler dislocation model consisting of glide dislocation loops so that the stress fields and interaction energy terms between dislocations in a two-phase medium could be evaluated. The energy associated with the coherent epilayer and that with the misfit dislocation are evaluated by the minimization of the total energy of the configuration with respect to the Burgers vectors of the surface dislocation loops. The numerical analysis is carried out for hemispherical islands of GaAs grown on (100) Si with a misfit dislocation of Burgers vector 3.84 Angstrom.