Numerical analyses of plasticity size effects have been carried out for different problems using a developed strain gradient crystal plasticiy theory. The theory employs higher order stresses as work conjugates to slip gradients and uses higher order boundary conditions. Problems on localization of plastic flow in a single crystal, grain boundary effects in a bicrystal, and grain size effects in a polycrystal are studied. Single crystals containing micro-scale voids have also been analyzed at different loading conditions with focus on the stress and deformation fields around the voids, on void growth and interaction between neighboring voids, and on a comparison between the developed strain gradient crystal plasticity theory and a discrete dislocation plasticity theory. Furthermore, voids and rigid inclusions in isotropic materials have been studied using a strain gradient plasticity theory for isotropic solids.