Journal
BIOMECHANICS AND MODELING IN MECHANOBIOLOGY
Volume 16, Issue 3, Pages 971-988Publisher
SPRINGER HEIDELBERG
DOI: 10.1007/s10237-016-0865-3
Keywords
Parameter estimation; Myocardium; Patient-specific modelling; Passive constitutive equations
Categories
Funding
- EPSRC [EP/M012492/1] Funding Source: UKRI
- British Heart Foundation [PG/13/37/30280, PG/11/101/29212] Funding Source: researchfish
- Engineering and Physical Sciences Research Council [EP/M012492/1] Funding Source: researchfish
Ask authors/readers for more resources
Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novelCF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were C-1 = 3000 +/- 1700 Pa and alpha = 45 +/- 25 (b(f) = 25 +/- 14, b(ft) = 11 +/- 6, b(t) = 9 +/- 5) for the HF cases and C-1 = 1700 Pa and alpha = 15 (b(f) = 8, b(ft) = 4, b(t) = 3) for the healthy case.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available