Method for Calibration of Left Ventricle Material Properties Using Three-Dimensional Echocardiography Endocardial Strains.

We sought to calibrate mechanical properties of left ventricle (LV) based on three-dimensional (3D) speckle tracking echocardiographic imaging data recorded from 16 segments defined by American Heart Association (AHA). The in vivo data were used to create finite element (FE) LV and biventricular (BV) models. The orientation of the fibers in the LV model was rule based, but diffusion tensor magnetic resonance imaging (MRI) data were used for the fiber directions in the BV model. A nonlinear fiber-reinforced constitutive equation was used to describe the passive behavior of the myocardium, whereas the active tension was described by a model based on tissue contraction (Tmax). isight was used for optimization, which used abaqus as the forward solver (Simulia, Providence, RI). The calibration of passive properties based on the end diastolic pressure volume relation (EDPVR) curve resulted in relatively good agreement (mean error = -0.04 ml). The difference between the experimental and computational strains decreased after segmental strain metrics, rather than global metrics, were used for calibration: for the LV model, the mean difference reduced from 0.129 to 0.046 (circumferential) and from 0.076 to 0.059 (longitudinal); for the BV model, the mean difference nearly did not change in the circumferential direction (0.061) but reduced in the longitudinal direction from 0.076 to 0.055. The calibration of mechanical properties for myocardium can be improved using segmental strain metrics. The importance of realistic fiber orientation and geometry for modeling of the LV was shown.

A modular inverse elastostatics approach to resolve the pressure-induced stress state for in vivo imaging based cardiovascular modeling.

Patient-specific biomechanical modeling of the cardiovascular system is complicated by the presence of a physiological pressure load given that the imaged tissue is in a pre-stressed and -strained state. Neglect of this prestressed state into solid tissue mechanics models leads to erroneous metrics (e.g. wall deformation, peak stress, wall shear stress) which in their turn are used for device design choices, risk assessment (e.g. procedure, rupture) and surgery planning. It is thus of utmost importance to incorporate this deformed and loaded tissue state into the computational models, which implies solving an inverse problem (calculating an undeformed geometry given the load and the deformed geometry). Methodologies to solve this inverse problem can be categorized into iterative and direct methodologies, both having their inherent advantages and disadvantages. Direct methodologies are typically based on the inverse elastostatics (IE) approach and offer a computationally efficient single shot methodology to compute the in vivo stress state. However, cumbersome and problem-specific derivations of the formulations and non-trivial access to the finite element analysis (FEA) code, especially for commercial products, refrain a broad implementation of these methodologies. For that reason, we developed a novel, modular IE approach and implemented this methodology in a commercial FEA solver with minor user subroutine interventions. The accuracy of this methodology was demonstrated in an arterial tube and porcine biventricular myocardium model. The computational power and efficiency of the methodology was shown by computing the in vivo stress and strain state, and the corresponding unloaded geometry, for two models containing multiple interacting incompressible, anisotropic (fiber-embedded) and hyperelastic material behaviors: a patient-specific abdominal aortic aneurysm and a full 4-chamber heart model.

Partial LVAD restores ventricular outputs and normalizes LV but not RV stress distributions in the acutely failing heart in silico.

PURPOSE: Heart failure is a worldwide epidemic that is unlikely to change as the population ages and life expectancy increases. We sought to detail significant recent improvements to the Dassault Systèmes Living Heart Model (LHM) and use the LHM to compute left ventricular (LV) and right ventricular (RV) myofiber stress distributions under the following 4 conditions: (1) normal cardiac function; (2) acute left heart failure (ALHF); (3) ALHF treated using an LV assist device (LVAD) flow rate of 2 L/min; and (4) ALHF treated using an LVAD flow rate of 4.5 L/min. METHODS AND RESULTS: Incorporating improved systolic myocardial material properties in the LHM resulted in its ability to simulate the Frank-Starling law of the heart. We decreased myocardial contractility in the LV myocardium so that LV ejection fraction decreased from 56% to 28%. This caused mean LV end diastolic (ED) stress to increase to 508% of normal, mean LV end systolic (ES) stress to increase to 113% of normal, mean RV ED stress to decrease to 94% of normal and RV ES to increase to 570% of normal. When ALHF in the model was treated with an LVAD flow rate of 4.5 L/min, most stress results normalized. Mean LV ED stress became 85% of normal, mean LV ES stress became 109% of normal and mean RV ED stress became 95% of normal. However, mean RV ES stress improved less dramatically (to 342% of normal values). CONCLUSIONS: These simulations strongly suggest that an LVAD is effective in normalizing LV stresses but not RV stresses that become elevated as a result of ALHF.

Human Cardiac Function Simulator for the Optimal Design of a Novel Annuloplasty Ring with a Sub-valvular Element for Correction of Ischemic Mitral Regurgitation.

Ischemic mitral regurgitation is associated with substantial risk of death. We sought to: (1) detail significant recent improvements to the Dassault Systèmes human cardiac function simulator (HCFS); (2) use the HCFS to simulate normal cardiac function as well as pathologic function in the setting of posterior left ventricular (LV) papillary muscle infarction; and (3) debut our novel device for correction of ischemic mitral regurgitation. We synthesized two recent studies of human myocardial mechanics. The first study presented the robust and integrative finite element HCFS. Its primary limitation was its poor diastolic performance with an LV ejection fraction below 20% caused by overly stiff ex vivo porcine tissue parameters. The second study derived improved diastolic myocardial material parameters using in vivo MRI data from five normal human subjects. We combined these models to simulate ischemic mitral regurgitation by computationally infarcting an LV region including the posterior papillary muscle. Contact between our novel device and the mitral valve apparatus was simulated using Dassault Systèmes SIMULIA software. Incorporating improved cardiac geometry and diastolic myocardial material properties in the HCFS resulted in a realistic LV ejection fraction of 55%. Simulating infarction of posterior papillary muscle caused regurgitant mitral valve mechanics. Implementation of our novel device corrected valve dysfunction. Improvements in the current study to the HCFS permit increasingly accurate study of myocardial mechanics. The first application of this simulator to abnormal human cardiac function suggests that our novel annuloplasty ring with a sub-valvular element will correct ischemic mitral regurgitation.

The Living Heart Project: A robust and integrative simulator for human heart function.

The heart is not only our most vital, but also our most complex organ: Precisely controlled by the interplay of electrical and mechanical fields, it consists of four chambers and four valves, which act in concert to regulate its filling, ejection, and overall pump function. While numerous computational models exist to study either the electrical or the mechanical response of its individual chambers, the integrative electro-mechanical response of the whole heart remains poorly understood. Here we present a proof-of-concept simulator for a four-chamber human heart model created from computer topography and magnetic resonance images. We illustrate the governing equations of excitation-contraction coupling and discretize them using a single, unified finite element environment. To illustrate the basic features of our model, we visualize the electrical potential and the mechanical deformation across the human heart throughout its cardiac cycle. To compare our simulation against common metrics of cardiac function, we extract the pressure-volume relationship and show that it agrees well with clinical observations. Our prototype model allows us to explore and understand the key features, physics, and technologies to create an integrative, predictive model of the living human heart. Ultimately, our simulator will open opportunities to probe landscapes of clinical parameters, and guide device design and treatment planning in cardiac diseases such as stenosis, regurgitation, or prolapse of the aortic, pulmonary, tricuspid, or mitral valve.