2.1. Computational Model
We utilized the Living Left Heart Human Model by Dassault Systémes Simulia Corporation (LLHH) which is capable of simulating LV performance, pressure-volume loops, and stress and strain analyses all of which correlate with clinical observations [
9,
10]. Our finite element model includes the AA, LV, left atrium, mitral valve, aortic root and the pericardium. The dynamic response is governed by realistic structural and blood flow physics, and the heart contraction is driven by electrical excitation. Blood is represented using a combination of three-dimensional hydrostatic fluid cavities for the heart chambers and system-level chambers to represent arterial and pulmonary compliances. Blood flow occurs inside a closed loop system between the chambers and the circulatory system through fluid link elements. Details of the model can be found under the “Simulation” and “Virtual Human” section at Dassault Systémes User Assistance, located at
http://help.3ds.com/.
The passive material response of the cardiac tissue uses an anisotropic hyperelastic formulation proposed by Holzapfel and Ogden as described in Equation 1) [
11]. The biaxial and triaxial experimental data published by Sommers et al. [
12] were used for an initial calibration and diastolic filling tests were used to augment the calibration of the eight material parameters: a, b, a
f, b
f, b
s, a
s, b
s, a
n, b
n describing the ventricular passive material properties based on the methods described in Klotz et al. [
13]
Equation 1: Passive material response of cardiac tissue. Ψdev is the deviatoric strain energy. The parameters a, b, af, bf, bs, as, bs, an, bn describing the ventricular passive material parameters.
The active tissue response contains length-dependent considerations of regional sarcomere lengths, affecting the stress components in the fiber and sheet directions in the constitutive model. The active-tissue material model intended to capture the Frank-Starling effect (i.e., the strength of the heart’s systolic contraction is directly proportional to its diastolic expansion) [
14]. The active contraction was simulated by adding stress in the direction of the muscle fiber defined by a time-varying model of elastance [
15] as follows:
Equation 2). T
max is the maximum isometric tension achieved at the longest sarcomere length and maximum peak intracellular calcium concentration. T
max (N/mm
2) is scalar factor representing the maximum active fiber stress or contractility in computational modeling.
where
Equation 2: Active Stress calculation. T
max (N/mm
2) is a scalar factor for myocardial contractility that represents the isometric tension achieved at the longest sarcomere length and maximum peak intracellular calcium concentration. Ca
0max is the peak intercellular calcium concentration. B governs the shape of peak isometric tension and sarcomere length relation. l
0 is the sarcomere length below which no active force develops. l
r is the initial sarcomere length. t
0 is time to reach the peak tension. m and b are coefficients that govern the relationship between of the linear relaxation duration and sarcomere length relaxation. E
ff is Lagrangian strain tension component aligned with the local muscle fiber direction [
15].
We simulated the mechanical constraints imposed by the pericardium by applying physiological boundary conditions on the ventricular epicardium to achieve realistic atrioventricular plane motion and radial inward motion of the epicardium as described in humans [
2]. Forty-nine clusters of nodes, evenly distributed on the epicardium surface, were constrained via a spring with higher stiffness when closer to the apex and lower stiffness when closer to the base [
16].
The heart was constrained via boundary conditions at the cut planes of the aortic root and pulmonary veins. Each cut plane was constrained relative to a central reference point and the reference point of the pulmonary veins was fixed. The aortic root was constrained from rotation but allowed to stretch. Aortic elasticity was modeled via a spring representing the AA stiffness. The stiffness of the spring was initially set at 0.5N/mm at baseline to achieve a realistic translation of the proximal aorta of 11.0mm during systole [
17].
The spring stiffness was increased to 10N/mm to model a stiff AA until a stationary aorta (stationary plane of the sino-tubular junction) was achieved. We performed three simulations under the following conditions:
A) The effect of a mobile AA using normal AA stiffness as elastic spring stiffness to constrain the aortic root motion with an amplitude of 11.0mm. B) the effect of stiffening the AA by immobilizing the AA at the sino-tubular junction. C) To model the effect of removing the pericardial confinement at the apex of the heart, the apical boundary conditions of the distal half of the pericardial sack were eliminated in model B), allowing free movement of the LV apex. Apex displacement was determined from the coordinates of the epicardial apex and mitral annulus plane center at maximum length at end-diastole and minimum length at end-systole and displacement was computed along this apex-base axis.
Myocardial strain was calculated as the relative length change between the diastolic and the systolic states. The LV strains were measured along the radial, circumferential, and longitudinal directions at 12 locations (three axial and four circumferential locations) at both the epicardium and endocardium. The averages of tensile strains were reported with positive values. Compressive strains were reported with negative values (
Table 1) and depicted as bar graphs (
Figure 1 and
Figure 2). Baseline values were in the reported range of normal human LV strains [
18].
Volumetric-averaged myofiber stress was calculated at end-systole in MPa (N/mm2) (
Table 3) and presented as a contour plot in LV parasternal long-axis cut planes (
Figure 3 and
Figure 4). Left ventricular pressures and volumes were computed (
Table 4) and depicted as Pressure-Volume Loops (
Figure 6). The area under the pressure-volume loop represents the total effective work (Joule) generated by ventricular shown in Equation 3:
Equation 3: Calculation of effective Stroke Work (SW) is area under the Pressure-Volume Loop, SV is Stroke Volume, and MAP is Mean Arterial Pressure.
Table 1.
Average strain.
Average strain |
Radial |
Circumferential |
Longitudinal |
Baseline Tmax 0.2 |
0.63±0.11 |
-0.20±0.05 |
-0.16±0.01 |
Stiff AA Tmax 0.2 |
0.50±0.11 |
-0.18±0.03 |
-0.08±0.05 |
Stiff AA and free apex Tmax 0.2 |
0.72±0.16 |
-0.18±0.05 |
-0.21±0.02 |
Baseline Tmax 0.2 vs Stiff AA Tmax 0.2 |
-0.13±0.02 |
0.01±0.02 |
0.08±0.06 |
Baseline Tmax 0.2 vs Stiff AA Tmax 0.2 (%) |
-20.21±2.39% |
-6.78±10.86% |
-48.44±36.88% |
Baseline Tmax 0.2 vs Stiff AA and free apex Tmax 0.2 |
0.09±0.06 |
0.02±0.04 |
-0.05±0.03 |
Baseline Tmax 0.2 vs Stiff AA and free apex Tmax 0.2 (%) |
14.36±9.73% |
-10.17±18.31% |
31.25±16.88% |
Table 2.
Longitudinal strain.
Table 2.
Longitudinal strain.
Longitudinal strain |
Septal |
Anterior |
Lateral |
Posterior |
Baseline Tmax 0.2 |
-0.17 |
-0.17 |
-0.15 |
-0.15 |
Stiff AA Tmax 0.2 |
-0.01 |
-0.10 |
-0.13 |
-0.09 |
Stiff AA and free apex Tmax 0.2 |
-0.21 |
-0.2 |
-0.19 |
-0.24 |
Baseline Tmax 0.2 vs Stiff AA Tmax 0.2 |
0.16 |
0.07 |
0.02 |
0.06 |
Baseline Tmax 0.2 vs Stiff AA Tmax 0.2 (%) |
-94.12% |
-41.18% |
-13.33% |
-40.00% |
Baseline Tmax 0.2 vs Stiff AA and free apex Tmax 0.2 |
-0.04 |
-0.03 |
-0.04 |
-0.09 |
Baseline Tmax 0.2 vs Stiff AA and free apex Tmax 0.2 (%) |
23.53% |
17.65% |
26.67% |
60.00% |
Table 3.
Myofiber stress.
Table 3.
Myofiber stress.
Stress |
Baseline |
Stiff AA |
Stiff AA and free apex |
Tmax 0.2 |
Tmax 0.2 |
Tmax 0.2 |
(MPa) |
0.056±0.036 |
0.076±0.042 |
0.062±0.038 |
vs. Baseline |
|
36.98±42.91% |
12.03±42.19% |
Table 4.
LV pressure and volume.
Table 4.
LV pressure and volume.
|
EDP |
EDV |
ESP |
ESV |
SVed-es |
SW |
|
(mmHg) |
(ml) |
(mmHg) |
(ml) |
(ml) |
(Joule) |
Baseline Tmax 0.2 |
11.85 |
158.30 |
117.10 |
66.10 |
92.20 |
8747.50 |
Stiff AA Tmax 0.2 |
12.86 |
159.60 |
106.40 |
77.40 |
82.20 |
7084.50 |
Stiff AA with free apex Tmax 0.2 |
11.25 |
157.00 |
116.60 |
62.86 |
94.14 |
8923.00 |
Baseline vs. Stiff AA Tmax 0.2 vs. Tmax 0.2 |
1.01 |
1.30 |
-10.70 |
11.30 |
-10.00 |
-1663.00 |
Baseline vs. Stiff AA Tmax 0.2 vs. Tmax 0.2 (%) |
8.52% |
0.82% |
-9.14% |
17.10% |
-10.85% |
-19.01% |
Baseline vs. Stiff AA with free apex Tmax 0.2 vs. Tmax 0.2 |
-0.60 |
-11.3 |
-0.50 |
-3.24 |
1.94 |
175.50 |
Baseline vs. Stiff AA with free apex Tmax 0.2 vs. Tmax 0.2 (%) |
-5.06% |
-0.82% |
0.43% |
-4.90% |
2.10% |
2.01% |
Figure 1.
Left ventricular strains. Figure 1. Left ventricular strain for: baseline simulation Tmax 0.2 N/mm2, simulation with stiff AA Tmax 0.2 N/mm2 and simulation with stiff AA and free apex Tmax 0.2 N/mm2. A) Radial strain (three radial locations) is depicted positive when wall thickening from diastole to systole. B) Circumferential strain (three circumferential locations) is depicted negative when circumference is reduced from diastole to systole. C) Longitudinal strain (four longitudinal locations) is depicted negative when apex base length is reduced from diastole to systole.
Figure 1.
Left ventricular strains. Figure 1. Left ventricular strain for: baseline simulation Tmax 0.2 N/mm2, simulation with stiff AA Tmax 0.2 N/mm2 and simulation with stiff AA and free apex Tmax 0.2 N/mm2. A) Radial strain (three radial locations) is depicted positive when wall thickening from diastole to systole. B) Circumferential strain (three circumferential locations) is depicted negative when circumference is reduced from diastole to systole. C) Longitudinal strain (four longitudinal locations) is depicted negative when apex base length is reduced from diastole to systole.
Figure 2.
Left ventricular longitudinal strain in regions. Figure 2. Left ventricular longitudinal strain in four longitudinal regions (Septal, Anterior, Lateral and Posterior strain) at Baseline Tmax 0.2 N/mm2, Stiff AA Tmax 0.2 N/mm2 and stiff AA with free apex Tmax 0.2 N/mm2.
Figure 2.
Left ventricular longitudinal strain in regions. Figure 2. Left ventricular longitudinal strain in four longitudinal regions (Septal, Anterior, Lateral and Posterior strain) at Baseline Tmax 0.2 N/mm2, Stiff AA Tmax 0.2 N/mm2 and stiff AA with free apex Tmax 0.2 N/mm2.
Figure 3.
Myofiber stress. Figure 3. Long-axis profile of the LV at end of systole showing contours of myofiber stress at end-systole. A) Baseline Tmax 0.2 N/mm2, B) Stiff AA Tmax 0.2 N/mm2 and (C) Stiff AA and free apex Tmax 0.2 N/mm2. Dotted line indicates baseline level of ascending aorta at end-diastole and level of apex.
Figure 3.
Myofiber stress. Figure 3. Long-axis profile of the LV at end of systole showing contours of myofiber stress at end-systole. A) Baseline Tmax 0.2 N/mm2, B) Stiff AA Tmax 0.2 N/mm2 and (C) Stiff AA and free apex Tmax 0.2 N/mm2. Dotted line indicates baseline level of ascending aorta at end-diastole and level of apex.
Figure 4.
Cross-sectional LV profile. Figure 4. Cross-sectional profile of the LV at center of the longitudinal LV axis and base of the papillary muscles. Grey color showing the end-diastolic shape at Baseline, green color showing the end-systolic shape. A) Baseline Tmax 0.2 N/mm2, diameter 59mm, B) Stiff AA Tmax 0.2 N/mm2, diameter 57mm and (C) Stiff AA and free apex Tmax 0.2 N/mm2, diameter 63mm.
Figure 4.
Cross-sectional LV profile. Figure 4. Cross-sectional profile of the LV at center of the longitudinal LV axis and base of the papillary muscles. Grey color showing the end-diastolic shape at Baseline, green color showing the end-systolic shape. A) Baseline Tmax 0.2 N/mm2, diameter 59mm, B) Stiff AA Tmax 0.2 N/mm2, diameter 57mm and (C) Stiff AA and free apex Tmax 0.2 N/mm2, diameter 63mm.
Figure 5.
Pressure volume loops. Figure 5. A) Comparison of pressure volume loop of left ventricle for simulation with mobile aorta (Baseline) Tmax 0.2 N/mm2 against simulation with stiff AA Tmax 0.2 N/mm2, B) Comparison of pressure volume loop of the left ventricle for simulation with mobile aorta Tmax 0.2 N/mm2 against stiff AA and free apex Tmax 0.2 N/mm2.
Figure 5.
Pressure volume loops. Figure 5. A) Comparison of pressure volume loop of left ventricle for simulation with mobile aorta (Baseline) Tmax 0.2 N/mm2 against simulation with stiff AA Tmax 0.2 N/mm2, B) Comparison of pressure volume loop of the left ventricle for simulation with mobile aorta Tmax 0.2 N/mm2 against stiff AA and free apex Tmax 0.2 N/mm2.
Figure 6.
Breaking the vicious cycle of stiffened AA to HFpEF symptoms.
Figure 6.
Breaking the vicious cycle of stiffened AA to HFpEF symptoms.