We present a new computational approach for modeling dislocation creep in metals using kinetic Monte Carlo simulation of thermally activated dislocation glide, extending the discrete dislocation dynamics method of plasticity. The method is used to study the problem of power-law creep in precipitation strengthened Aluminum single crystals. The new model predicts creep rates and stress exponents consistent with their known ranges from experiments.