Modeling a distribution of point defects as misfitting inclusions in stressed solids