Here in this paper, a multiscale framework based on crystal plasticity is proposed coupling with thermal activation mechanism, as well as the continuum damage mechanism. The microscopic hardening phenomenon is revealed by a dislocation density evolution model, which is constructed according to the temperature. The physical-based exponential function of shear strain rate is posed to describe the thermal material behavior replacing the general phenomenological power-law equation. A 3D spatial distribution of stress, strain and damage is presented in the finite element method, parameters of which are previously determined by a RVE calculation and fitting test compared to the experimental data. Finally, some discussions of stress heterogeneity and texture evolution are proposed and conclusions are made.