Conscious efforts on reduction of greenhouse gas emissions have led to an energy transition to renewable energy, however uncertainties of renewable energy production have resulted in higher thermal cycling demands from conventional power plants. Thermal load cycling at high temperature regions of steam turbine components leads to enhanced creep-fatigue damage accumulation. It is well established that such damage mechanism is numerically best predicted by unified constitutive modeling including damage as a variable as per the formalism of continuum damage mechanics at an expense of considerable computational efforts using finite element analysis. In this paper, the non-iterative Asymptotic Numerical Method (ANM), currently limited to partial cycle analysis with linear hardening plasticity model, is proposed for the first time to address cyclic viscoplasticity problems including damage capable of handling multiple cycles and most generalized loading conditions. Regularization techniques additionally necessary to implement loading-unloading-reloading criteria and advanced constitutive models etc. are presented. The constitutive model chosen for the formulation includes the non-linear multiple back stress variable modified Chaboche model to include damage combined with modified Chaboche-Rousselier isotropic hardening model to include damage, power law for viscoplasticity, Lemaitre’s damage potential and Kachanov-Rabotnov’s creep damage law. The method is verified with defined error measures and then applied to two high pressure steam turbine rotors, one with and another without thermal stress relief groove (TSRG) at the inlet under service type loading conditions to study the beneficial effect of the TSRG on creep-fatigue damage evolution. The accumulated errors of the proposed ANM and computational time are compared to a conventional Newton-Raphson solution.