The time evolution of the auxiliary stress tensor Qm in each viscoelastic branch is given by the rate
here, Siso,m is the isochoric second Piola-Kirchhoff stress tensor in the branch
m. These tensors are derived from the strain energy density in the main hyperelastic branch and the
energy factors βm as