Nickel-based single-crystal (SC) superalloys exhibited excellent exceptional mechanical properties at high temperatures due to the elimination of internal grain boundaries, contributed a strong orientation-dependent material response. The anisotropy of SC superalloys was modeled viscoplastically from a macroscopic viewpoint based on the Kelvin decomposition theory [1] which was a decomposition of the stress space according to the elastic matrix eigen-directions to control the viscoplastic flow by Kelvin stress decoupled from each other. Compared to the classical phenomenological macro model, the proposed model effectively captures the slip deformation mechanism of SC superalloys with the inherent ability to simulate anisotropic because of the two criterions framework controlled by Kelvin stress. Compared with others, the proposed model was able to simulate time-dependent inelastic deformation and cyclic deformation behavior under complex loading. The kinematic hardening and isotropic hardening models incorporated microscopic quantities, such as dislocation density and channel phase width, connecting the macroscopic mechanical response with the microscopic state to achieve multiscale constitutive modelling. The parameter identification and finite element implementation were conducted on a SC superalloy [2]. Simulation results demonstrated the accuracy of the proposed model in predicting deformation behavior under various orientations, rate-dependent effects, isothermal and non-isothermal cyclic deformation. Comparison with the classical anisotropic matrix macroscopic phenomenological approaches highlights the superior capability of the proposed model to simulate the orientation-dependent mechanical properties of single-crystal alloys.