“…The numerical integration of the evolution equations (7) and (11) for m e and b e follows the predictor-corrector scheme used to compute the deviatoric elastic left Cauchy-Green tensor in [Rubin and Bodner, 2002, Appendix A]. The recent modifications to obtain the exact result in the case of zero dissipation [Rubin and Papes, 2011;Hollenstein et al, 2013;Flynn and Rubin, 2014], which adopt the idea of a relative deformation gradient from the configurations at time t n to t n+1 [see Simo, 1992;Simo and Hughes, 2000], have been included.…”