Simulating hyperelasticity and fractional viscoelasticity in the human heart

Biomechanics plays an important role in the diagnosis and treatment of pathological conditions of the heart. Computational models are paving the way for personalized therapeutic treatment but they rely on accurate constitutive equations for predicting their biomechanical behavior. Even so, viscoelasticity remains under-explored in computational modeling despite experimental observations. To facilitate the viscoelastic modeling of cardiovascular soft tissues, we previously developed a fractional viscoelastic modeling approach, which extends existing hyperelastic models. This has comparable computational costs to the conventional hyperelastic model and only requires two additional material parameters for the viscoelastic response. This approach was demonstrated to be able to accurately capture the viscoelastic response of the human myocardium. However, the numerical properties of this fractional viscoelastic approach have not yet been examined. In this work, we present its implementation in Finite Element Analysis, examine its numerical properties in uniaxial extension and 2D inflation test examples, and examine its physiological implication in a computational model of an idealized left ventricle in a fully idealized circulatory system. Optimal convergence properties were observed and the importance of viscoelasticity during passive filling, ventricular motion, and regional fiber strain and stresses were explained.