A finite element model of the human left ventricular systole

• Published in: Computer methods in biomechanics and biomedical engineering Vol. 9; no. 5; pp. 319 - 341
• Main Authors: Dorri, F., Niederer, P. F., Lunkenheimer, P. P.
• Format: Journal Article
• Language: English
• Published: England Taylor & Francis Group 01-10-2006
• Subjects: Computer Simulation; Contraction; Elasticity; Finite Element Analysis; Heart; Humans; Models, Cardiovascular; Myocardial Contraction - physiology; Myocardium; Stress, Mechanical; Systole; Systole - physiology; Ventricular Function; Ventricular Function, Left - physiology
• ISSN: 1025-5842; 1476-8259
• DOI: 10.1080/10255840600960546

Local wall stress is the pivotal determinant of the heart muscle's systolic function. Under in vivo conditions, however, such stresses cannot be measured systematically and quantitatively. In contrast, imaging techniques based on magnetic resonance (MR) allow the determination of the deformation pattern of the left ventricle (LV) in vivo with high accuracy. The question arises to what extent deformation measurements are significant and might provide a possibility for future diagnostic purposes.
The contractile forces cause deformation of LV myocardial tissue in terms of wall thickening, longitudinal shortening, twisting rotation and radial constriction. The myocardium is thereby understood to act as a densely interlaced mesh. Yet, whole cycle image sequences display a distribution of wall strains as function of space and time heralding a significant amount of inhomogeneity even under healthy conditions. We made similar observations previously by direct measurement of local contractile activity. The major reasons for these inhomogeneities derive from regional deviations of the ventricular walls from an ideal spheroidal shape along with marked disparities in focal fibre orientation.
In response to a lack of diagnostic tools able to measure wall stress in clinical routine, this communication is aimed at an analysis and functional interpretation of the deformation pattern of an exemplary human heart at end-systole. To this end, the finite element (FE) method was used to simulate the three-dimensional deformations of the left ventricular myocardium due to contractile fibre forces at end-systole. The anisotropy associated with the fibre structure of the myocardial tissue was included in the form of a fibre orientation vector field which was reconstructed from the measured fibre trajectories in a post mortem human heart. Contraction was modelled by an additive second Piola-Kirchhoff active stress tensor.
As a first conclusion, it became evident that longitudinal fibre forces, cross-fibre forces and shear along with systolic fibre rearrangement have to be taken into account for a useful modelling of systolic deformation. Second, a realistic geometry and fibre architecture lead to typical and substantially inhomogeneous deformation patterns as they are recorded in real hearts. We therefore, expect that the measurement of systolic deformation might provide useful diagnostic information.