The dislocation creep represents high temperature plastic deformation of metallic materials. However, the subject is usually approached with constitutive and phenomenological formulations. The notable exceptions are the recent use of numerical modeling. For example, single-crystal discrete dislocation dynamics (DDD) has been adopted to involve dislocation glide and climb together with diffusion of vacancies. This naturally provides a physics-based framework for dislocation creep. However, the existing model(s) are limited to single crystals and hence are limited in scope. This study extended a single crystal dislocation creep model into polycrystalline ensemble. This was then used to simulate steady state creep rate (ε˙) with respect to creep stress (σ), temperature (T) and grain size (d) in polycrystalline aluminum. Our virtual experiments showed an Arrhenius relationship between ε˙ and T, and a power law scaling with σ. Further, our simulations also revealed a power law increase in ε˙ with d. This is striking, as available experimental data indicate both grain size dependence as well as independence. The answer to this apparent contradiction also emerged from our numerical simulations. It was shown that ε˙ is controlled by interactions of dislocations with static obstacles and grain boundaries. Dominance of static obstacles, in particular, significantly diminished the grain size effect. The role of dislocation mean free path on the pinning and dislocation creep behavior of polycrystalline metallic material was thus mechanistically established.