In order to theoretically analyze the creep behavior of high Cr steel at $600^{circ}C$, a unified elasto-viscoplastic constitutive model based on the consideration of dislocation density is proposed. A combination of a kinetic equation describing the mechanical response of a material at a given microstructure in terms of dislocation glide and evolution equations for internal variables characterizing the microstructure provides the constitutive equations of the model. Microstructural features of the material such as the grain size and spacing between second phase particles are directly implemented in the constitutive equations. The internal variables are associated with the total dislocation density in a simple model. The model has a modular structure and can be adjusted to describe a creep behavior using the material parameters obtained from uniaxial tensile tests.