Statistics of A multiscale mean-field homogenization method for fiber-reinforced composites with gradient-enhanced damage models