Orthotropic active strain model for the numerical simulation of cardiac biomechanics

A model for the active deformation of cardiac tissue considering orthotropic constitutive laws is introduced and studied. In particular, the passive mechanical properties of the myocardium are described by the Holzapfel-Ogden relation, whereas the activation model is based on the concept of active strain. There, an incompatible intermediate configuration is considered, which entails a multiplicative decomposition between active and passive deformation gradients. The underlying Euler-Lagrange equations for minimizing the total energy are written in terms of these deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiomyocytes. The active strain formulation is compared with the classical active stress model from both numerical and modeling perspectives. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed, and different mechanical activation functions are considered. Taylor-Hood and MINI finite elements are used in the discretization of the overall mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main features of the phenomenon, while allowing savings in computational costs.