In this study, we present a homogenization model for the macroscopic rheological behavior of non colloidal suspensions of initially spherical, viscoelastic particles in yield stress fluids, which are subjected to uniform flow conditions. The constitutive behavior of the suspending fluid is characterized by the Herschel-Bulkley (HB) model, and the particles are assumed to be neutrally buoyant solids characterized by the finite-strain Kelvin-Voigt viscoelastic model. We make use of the "linear comparison composite" variational technique of Ponte Castaneda (1991) to approximate the instantaneous response of the suspensions of viscoelastic particles in the HB fluid by that of a fictitious suspension consisting of the same particles distributed identically in a Newtonian fluid with a suitably-chosen viscosity. The response of the latter suspension is then estimated by the homogenization model recently developed by Avazmohammadi and Ponte Castaneda (2015), which, when combined with appropriate evolution laws for the relevant microstructural variables, provides a complete characterization of the time-dependent response of the actual suspensions. With the objective of illustrating the key features of our model, we consider the particular case of suspensions of elastic particles in HB fluids under shear flow conditions. The results provide a broad picture of the influence of the HB fluid and particle constitutive properties, as well as of the particle volume fraction on the effective time-dependent and steady-state behaviors of the suspensions. For the special case of non-deformable particles, our model predicts that the suspensions behave like HB fluids with modified properties, consistent with the results of Chateau et al. (2008). (C) 2016 Elsevier B.V. All rights reserved.