Viscoelasticity using reactive constrained solid mixtures.

This study presents a framework for viscoelasticity where the free energy density depends on the stored energy of intact strong and weak bonds, where weak bonds break and reform in response to loading. The stress is evaluated by differentiating the free energy density with respect to the deformation gradient, similar to the conventional approach for hyperelasticity. The breaking and reformation of weak bonds is treated as a reaction governed by the axiom of mass balance, where the constitutive relation for the mass supply governs the bond kinetics. The evolving mass contents of these weak bonds serve as observable state variables. Weak bonds reform in an energy-free and stress-free state, therefore their reference configuration is given by the current configuration at the time of their reformation. A principal advantage of this formulation is the availability of a strain energy density function that depends only on observable state variables, also allowing for a separation of the contributions of strong and weak bonds. The Clausius-Duhem inequality is satisfied by requiring that the net free energy from all breaking bonds must be decreasing at all times. In the limit of infinitesimal strains, linear stress-strain responses and first-order kinetics for breaking and reforming of weak bonds, the reactive framework reduces exactly to classical linear viscoelasticity. For large strains, the reactive and classical quasilinear viscoelasticity theories produce different equations, though responses to standard loading configurations behave similarly. This formulation complements existing tools for modeling the nonlinear viscoelastic response of biological soft tissues under large deformations.