On the AIC-based model reduction for the general Holzapfel–Ogden myocardial constitutive law

Abstract

Constitutive laws that describe the mechanical responses of cardiac tissue under loading hold the key to accurately model the biomechanical behaviour of the heart. There have been ample choices of phenomenological constitutive laws derived from experiments, some of which are quite sophisticated and include effects of microscopic fibre structures of the myocardium. A typical example is the strain-invariant-based Holzapfel–Ogden 2009 model that is excellently fitted to simple shear tests. It has been widely used and regarded as the state-of-the-art constitutive law for myocardium. However, there has been no analysis to show if it has both adequate descriptive and predictive capabilities for other tissue tests of myocardium. Indeed, such an analysis is important for any constitutive laws for clinically useful computational simulations. In this work, we perform such an analysis using combinations of tissue tests, uniaxial tension, biaxial tension and simple shear from three different sets of myocardial tissue studies. Starting from the general 14-parameter myocardial constitutive law developed by Holzapfel and Ogden, denoted as the general HO model, we show that this model has good descriptive and predictive capabilities for all the experimental tests. However, to reliably determine all 14 parameters of the model from experiments remains a great challenge. Our aim is to reduce the constitutive law using Akaike information criterion, to maintain its mechanical integrity whilst achieving minimal computational cost. A competent constitutive law should have descriptive and predictive capabilities for different tissue tests. By competent, we mean the model has least terms but is still able to describe and predict experimental data. We also investigate the optimal combinations of tissue tests for a given constitutive model. For example, our results show that using one of the reduced HO models, one may need just one shear response (along normal-fibre direction) and one biaxial stretch (ratio of 1 mean fibre : 1 cross-fibre) to satisfactorily describe Sommer et al. human myocardial mechanical properties. Our study suggests that single-state tests (i.e. simple shear or stretching only) are insufficient to determine the myocardium responses. We also found it is important to consider transmural fibre rotations within each myocardial sample of tests during the fitting process. This is done by excluding un-stretched fibres using an “effective fibre ratio”, which depends on the sample size, shape, local myofibre architecture and loading conditions. We conclude that a competent myocardium material model can be obtained from the general HO model using AIC analysis and a suitable combination of tissue tests.