Effects of using the unloaded configuration in predicting the in vivo diastolic properties of the heart

Computational models are increasingly being used to investigate the mechanical properties of cardiac tissue. While much insight has been gained from these studies, one important limitation associated with computational modeling arises when using in vivo images of the heart to generate the reference state of the model. An unloaded reference configuration is needed to accurately represent the deformation of the heart. However, it is rare for a beating heart to actually reach a zero-pressure state during the cardiac cycle. To overcome this, a computational technique was adapted to determine the unloaded configuration of an in vivo porcine left ventricle (LV). In the current study, in vivo measurements were acquired using magnetic resonance images (MRI) and synchronous pressure catheterization in the LV (N = 5). The overall goal was to quantify the effects of using early–diastolic filling as the reference configuration (common assumption used in modeling) versus using the unloaded reference configuration for predicting the in vivo properties of LV myocardium. This was accomplished by using optimization to minimize the difference between MRI measured and finite element predicted strains and cavity volumes. The results show that when using the unloaded reference configuration, the computational method predicts material properties for LV myocardium that are softer and less anisotropic than when using the early-diastolic filling reference configuration. This indicates that the choice of reference configuration could have a significant impact on capturing the realistic mechanical response of the heart.     

Passive material properties of intact ventricular myocardium determined from a cylindrical model

The equatorial region of the canine left ventricle was modeled as a thick-walled cylinder consisting of an incompressible hyperelastic material with homogeneous exponential properties. The anisotropic properties of the passive myocardium were assumed to be locally transversely isotropic with respect to a fiber axis whose orientation varied linearly across the wall. Simultaneous inflation, extension, and torsion were applied to the cylinder to produce epicardial strains that were measured previously in the potassium-arrested dog heart. Residual stress in the unloaded state was included by considering the stress-free configuration to be a warped cylindrical arc. In the special case of isotropic material properties, torsion and residual stress both significantly reduced the high circumferential stress peaks predicted at the endocardium by previous models. However, a resultant axial force and moment were necessary to cause the observed epicardial deformations. Therefore, the anisotropic material parameters were found that minimized these resultants and allowed the prescribed displacements to occur subject to the known ventricular pressure loads. The global minimum solution of this parameter optimization problem indicated that the stiffness of passive myocardium (defined for a 20 percent equibiaxial extension) would be 2.4 to 6.6 times greater in the fiber direction than in the transverse plane for a broad range of assumed fiber angle distributions and residual stresses. This agrees with the results of biaxial tissue testing. The predicted transmural distributions of fiber stress were relatively flat with slight peaks in the subepicardium, and the fiber strain profiles agreed closely with experimentally observed sarcomere length distributions. The results indicate that torsion, residual stress and material anisotropy associated with the fiber architecture all can act to reduce endocardial stress gradients in the passive left ventricle.     

Sarcomere length changes in a 3D mathematical model of the pig ventricles

Measurements of the geometry and fibrous-sheet structure of the left and right ventricles of the pig heart are fitted with a finite element model. Mechanical changes during the heart cycle are computed by solving the equations of motion under specified ventricular boundary conditions and using experimentally defined constitutive laws for the active and passive material properties of myocardial tissue. The resulting patterns of deformation, such as axial torsion and changes in wall thickness and base-apex length, are consistent with experimental observations. The model can therefore be used to predict sarcomere length changes and other strain patterns throughout the myocardium and throughout the cardiac cycle. Here we present sarcomere length changes at a limited number of material points within the wall. Sarcomere length typically varies by 10% above and below the unloaded length; although under the boundary conditions imposed in the current model the midwall circumferentially oriented sarcomere lengths increased by up to 20% at end diastole. We provide web-access details for a downloadable software program designed to provide more extensive information on mechanical deformation, such as the principal strains and muscle fibre cross-sectional area changes during the cardiac cycle.     

Left ventricular geometric remodeling and residual stress in the rat heart.

Theoretical considerations and observations of residual stress suggest that geometric remodeling in the heart may also alter residual stress and strain. We investigated whether changes in left ventricular geometry during physiologic growth were associated with corresponding changes in myocardial residual strain. In anesthetized rats from eight age groups ranging from 2-25+ weeks, the heart was arrested and isolated, and equatorial slices were obtained. The geometry of the intact, unloaded state was recorded, as well as the "opening angle" of the stress-free configuration after radial resection of the tissue slice. The tissue was fixed and embedded for histological examination of collagen area fraction. Heart weight increased 10-fold with age and unloaded internal radius increased almost 4-fold. However, wall thickness increased only 66 percent, so that the ratio of wall thickness to internal radius decreased significantly from 2.22 +/- 0.29 (mean +/- SD) at 2 weeks to 0.81 +/- 0.47 at 25 weeks. Opening angle of the stress-free slice decreased significantly from 87 +/- 16 deg at 2 weeks to 51 +/- 16 deg, and correlated linearly with wall thickness/radius ratio. Collagen area fraction increased with age. Hence physiologic ventricular remodeling in rats decreases myocardial residual strain in proportion to the relative reduction in wall thickness-radius ratio.     

Bioinjection treatment: Effects of post-injection residual stress on left ventricular wall stress

Injection of biomaterials into diseased myocardium has been associated with decreased myofiber stress, restored left ventricular (LV) geometry and improved LV function. However, its exact mechanism(s) of action remained unclear. In this work, we present the first patient-specific computational model of biomaterial injection that accounts for the possibility of residual strain and stress introduced by this treatment. We show that the presence of residual stress can create more heterogeneous regional myofiber stress and strain fields. Our simulation results show that the treatment generates low stress and stretch areas between injection sites, and high stress and stretch areas between the injections and both the endocardium and epicardium. Globally, these local changes are translated into an increase in average myofiber stress and its standard deviation (from 6.9±4.6 to 11.2±48.8 kPa and 30±15 to 35.1±50.9 kPa at end-diastole and end-systole, respectively). We also show that the myofiber stress field is sensitive to the void-to-size ratio. For a constant void size, the myofiber stress field became less heterogeneous with decreasing injection volume. These results suggest that the residual stress and strain possibly generated by biomaterial injection treatment can have large effects on the regional myocardial stress and strain fields, which may be important in the remodeling process.     

Effects of myocardial sheetlet sliding on left ventricular function

Left ventricle myocardium has a complex micro-architecture, which was revealed to consist of myocyte bundles arranged in a series of laminar sheetlets. Recent imaging studies demonstrated that these sheetlets re-orientated and likely slided over each other during the deformations between systole and diastole, and that sheetlet dynamics were altered during cardiomyopathy. However, the biomechanical effect of sheetlet sliding is not well-understood, which is the focus here. We conducted finite element simulations of the left ventricle (LV) coupled with a windkessel lumped parameter model to study sheetlet sliding, based on cardiac MRI of a healthy human subject, and modifications to account for hypertrophic and dilated geometric changes during cardiomyopathy remodeling. We modeled sheetlet sliding as a reduced shear stiffness in the sheet-normal direction and observed that (1) the diastolic sheetlet orientations must depart from alignment with the LV wall plane in order for sheetlet sliding to have an effect on cardiac function, that (2) sheetlet sliding modestly aided cardiac function of the healthy and dilated hearts, in terms of ejection fraction, stroke volume, and systolic pressure generation, but its effects were amplified during hypertrophic cardiomyopathy and diminished during dilated cardiomyopathy due to both sheetlet angle configuration and geometry, and that (3) where sheetlet sliding aided cardiac function, it increased tissue stresses, particularly in the myofibre direction. We speculate that sheetlet sliding is a tissue architectural adaptation to allow easier deformations of the LV walls so that LV wall stiffness will not hinder function, and to provide a balance between function and tissue stresses. A limitation here is that sheetlet sliding is modeled as a simple reduction in shear stiffness, without consideration of micro-scale sheetlet mechanics and dynamics.

     

Development of an in vivo method for determining material properties of passive myocardium

Calculation of mechanical stresses and strains in the left ventricular (LV) myocardium by the finite element (FE) method relies on adequate knowledge of the material properties of myocardial tissue. In this paper, we present a model-based estimation procedure to characterize the stress-strain relationship in passive LV myocardium. A 3D FE model of the LV myocardium was used, which included morphological fiber and sheet structure and a nonlinear orthotropic constitutive law with different stiffness in the fiber, sheet, and sheet-normal directions. The estimation method was based on measured wall strains. We analyzed the method's ability to estimate the material parameters by generating a set of synthetic strain data by simulating the LV inflation phase with known material parameters. In this way we were able to verify the correctness of the solution and to analyze the effects of measurement and model error on the solution accuracy and stability. A sensitivity analysis was performed to investigate the observability of the material parameters and to determine which parameters to estimate. The results showed a high degree of coupling between the parameters governing the stiffness in each direction. Thus, only one parameter in each of the three directions was estimated. For the tested magnitudes of added noise and introduced model errors, the resulting estimated stress-strain characteristics in the fiber and sheet directions converged with good accuracy to the known relationship. The sheet-normal stress-strain relationship had a higher degree of uncertainty as more noise was added and model error was introduced.     

A two-phase finite element model of the diastolic left ventricle

A porous medium finite element model of the passive left ventricle is presented. The model is axisymmetric and allows for finite deformation, including torsion about the axis of symmetry. An anisotropic quasi-linear viscoelastic constitutive relation is implemented in the model. The model accounts for changing fibre orientation across the myocardial wall. During passive filling, the apex rotates in a clockwise direction relative to the base for an observer looking from apex to base. Within an intraventricular pressure range of 0-3 kPa the rotation angle of all nodes remained below 0.1 rad. Diastolic viscoelasticity of myocardial tissue is shown to reduce transmural differences of preload-induced sarcomere stretch and to generate residual stresses in an unloaded ventricular wall, consistent with the observation of opening angles seen when the heart is slit open. It is shown that the ventricular model stiffens following an increase of the intracoronary blood volume. At a given left ventricular volume, left ventricular pressure increases from 1.5 to 2.0 kPa when raising the intracoronary blood volume from 9 to 14 ml (100 g)-1 left ventricle.     

Importance of material parameters and strain energy function on the wall stresses in the left ventricle

Patient-specific estimates of the stress distribution in the left ventricles (LV) may have important applications for therapy planning, but computing the stress generally requires knowledge of the material behaviour. The passive stress-strain relation of myocardial tissue has been characterized by a number of models, but material parameters (MPs) remain difficult to estimate. The aim of this study is to implement a zero-pressure algorithm to reconstruct numerically the stress distribution in the LV without precise knowledge of MPs. We investigate the sensitivity of the stress distribution to variations in the different sets of constitutive parameters. We show that the sensitivity of the LV stresses to MPs can be marginal for an isotropic constitutive model. However, when using a transversely isotropic exponential strain energy function, the LV stresses become sensitive to MPs, especially to the linear elastic coefficient before the exponential function. This indicates that in-vivo identification efforts should focus mostly on this MP for the development of patient-specific finite-element analysis.     

Theoretical analysis of the effects of a radial activation wave and twisting motion on the mechanics of the left ventricle

The mechanical effects resulting from the normal transmural delay of electrical depolarization of the myocardium are investigated. An activation sequence having a finite radial propagation velocity is introduced into the equations of ventricular mechanics. The resulting system of coupled integral equations is solved using a perturbation method based on the small ratio of transmural propagation time to cardiac period. Numerical calculations are performed using cavity pressure and volume waveforms characteristic of the canine left ventricle (LV), for both simultaneous and delayed activation of fiber layers. The results show that a finite transmural electrical propagation velocity tends to: (i) equalize the transmural distribution of sarcomere length during systole; (ii) equalize the transmural distribution of fiber external work/vol; and (iii) insignificantly affect myocardial tissue pressure. Calculations are also performed to investigate the mechanical effects resulting from the application of an externally applied moment that prevents LV torsion. Those results are highly dependent on the transmural distribution of sarcomere length in the stress-free reference state (unloaded diastole). When we assume a uniform distribution, then normal torsion acting with normal activation delay tends to: (i) increase the magnitude of fiber strain in the subendocardium and decrease it in the subepicardium; (ii) equalize the transmural distribution of fiber external work/vol; and (iii) lower myocardial tissue pressure. The normally occurring transmural delay of activation tends to lessen endocardial O2 demand, while the normally occurring torsion further lessens that demand and improves O2 supply.     

The effect of residual strain on the diastolic function of the left ventricle as predicted by a structural model

The unloaded heart is not stress-free. It is subjected to residual stress and strain. Their extent and influence on the global performance of the left ventricle and on local phenomena in the ventricular wall are studied by model simulation. The analysis focuses on the equatorial region of the ventricle, with an approximate thick-walled cylindrical geometry. The in vivo myocardium is considered to be incompressible, consisting of fibers embedded in a fluid matrix, with transmurally varying anisotropic microstructure in accordance with morphological characteristics.

The results show that residual strain is transmurally distributed with a pattern and magnitude which agree well with measurements. The calculated residual strains are within mean ± one standard deviation of the measured ones. Their magnitude was found to increase with increasing opening angle and with increasing wall thickness. The residual strain was found to have several effects on ventricular function: At volumes higher than the reference one it gives rise to more uniform transmural distributions of stress and intramyocardial pressure; it causes about 50% increase in the ventricular compliance at high volumes and doubles the suction of atrial blood at low volumes, thus facilitating the diastolic filling. In addition, residual strains cause bias of in vivo measured strains from their true values. This may significantly affect physiological interpretation of measured ventricular deformations.

In conclusion, the present structural analysis predicts that residual strain has favorable effect on left-ventricular diastolic performance, and gives rise to more uniform ventricular stress distribution.     

Stress Distribution in the Canine Left Ventricle during Diastole and Systole

A model is proposed for stress analysis of the left ventricular wall (LV wall) based on the realistic assumption that the myocardium is essentially composed of fiber elements which carry only axial tension and vary in orientation through the wall. Stress analysis based on such a model requires an extensive study of muscle fiber orientation and curvature through the myocardium. Accordingly, the principal curvatures were studied at a local site near the equator in ten dog hearts rapidly fixed in situ at end diastole and end systole; the fiber orientation for these hearts had already been established in a previous study. The principal radii of curvature were (a) measured by fitting templates to the endocardial and epicardial wall surfaces in the circumferential and longitudinal directions and (b) computed from measured lengths of semiaxes of ellipsoids of revolution representing the LV wall (“ellipsoid” data). The wall was regarded as a tethered set of nested shells, each having a unique fiber orientation. Results indicate the following. (a) Fiber curvature, k, is maximum at midwall at end systole; this peak shifts towards endocardium at end diastole. (b) The pressure or radial stress through the wall decreases more rapidly near the endocardium than near the epicardium at end diastole and at end systole when a constant tension is assumed for each fiber through the wall. (c) At end diastole the curve for the circumferential stress vs. wall thickness is convex with a maximum at midwall. In the longitudinal direction the stress distribution curve is concave with a minimum at midwall. Similar distributions are obtained at end systole when a constant tension is assumed for each fiber through the wall. (d) The curvature and stress distributions obtained by direct measurements at a selected local site agree well with those computed from “ellipsoid” data.     

A 2D finite element model of the interventricular septum under normal and abnormal loading

The interventricular septum is the structure that separates the left and right ventricles of the heart. Under normal loading conditions, it is concave to the left ventricle, but under abnormal loading the septum flattens and occasionally inverts. In the past, the septum has frequently been modelled as integral to the left ventricle with the effects of pressure from the right ventricle being ignored. Under abnormal loading, the septum has been described as behaving equivalent to a "flapping sail". There has been no consideration of structural behaviour under these conditions. A 2-D plane stress FE model of the septum was used to investigate the difference in structural behaviour of the septum during diastole between normal and abnormal loading. The biaxial stress patterns that develop are distinctively disparate. Under normal loading, the septum behaves much like a thick-walled cylinder subject to internal and external pressure, with the resulting stresses being circumferential tension and radial compression, both varying with radius. These stresses are very low throughout most of diastole. However, under abnormal loading, the septum behaves in an arch-like fashion, with high compressive stresses almost circumferential in direction, combined with radial compression. We conclude that right ventricular pressures cause bending effects in the wall of the heart, and that under abnormal loading, the compressive stresses that develop in the septum may lead to an understanding of certain, previously unexplained, pathological conditions.     

Residual strain in ischemic ventricular myocardium.

Structural remodeling during acute myocardial infarction affects ventricular wall stress and strain. To see whether acute myocardial infarction alters residual stress and strain in the left ventricle (LV), we measured opening angles in rat hearts after 30 minutes of left coronary artery occlusion. The mean opening angle in 18 ischemic hearts (51 +/- 20 deg) was significantly greater than in five sham-operated controls (29 +/- 11 deg, P < 0.05). To determine whether these alterations in residual strain may be associated with strain softening caused by systolic overstretch of the noncontracting ischemic tissue, we also measured opening angles in isolated hearts that had been passively inflated to high LV pressures (120 mmHg). The mean opening angle of the strain-softened hearts was not significantly different from the sham-operated hearts (34 +/- 27 deg, P = 0.74). Mean collagen area fractions in the myocardium were not significantly different between ischemic hearts (0.027 +/- 0.014) and the nonischemic group (0.022 +/- 0.011). Although there were significant differences in opening angles measured with ischemia, they do not appear to be a result of altered extracellular collagen content or softening associated with overstretch. Thus, there is a significant change in residual strain associated with acute ischemia that may be related to changes in collagen fiber structure, myocyte structure, or metabolic state.     

A 2D Finite Element Model of the Interventricular Septum Under Normal and Abnormal Loading

Abstract

The interventricular septum is the structure that separates the left and right ventricles of the heart. Under normal loading conditions, it is concave to the left ventricle, but under abnormal loading the septum flattens and occasionally inverts. In the past, the septum has frequently been modelled as integral to the left ventricle with the effects of pressure from the right ventricle being ignored. Under abnormal loading, the septum has been described as behaving equivalent to a “flapping sail”. There has been no consideration of structural behaviour under these conditions. A 2-D plane stress FE model of the septum was used to investigate the difference in structural behaviour of the septum during diastole between normal and abnormal loading. The biaxial stress patterns that develop are distinctively disparate. Under normal loading, the septum behaves much like a thick-walled cylinder subject to internal and external pressure, with the resulting stresses being circumferential tension and radial compression, both varying with radius. These stresses are very low throughout most of diastole. However, under abnormal loading, the septum behaves in an arch-like fashion, with high compressive stresses almost circumferential in direction, combined with radial compression. We conclude that right ventricular pressures cause bending effects in the wall of the heart, and that under abnormal loading, the compressive stresses that develop in the septum may lead to an understanding of certain, previously unexplained, pathological conditions.     

In-vivo mitral annuloplasty ring transducer: Implications for implantation and annular downsizing

Mitral annuloplasty has been a keystone to the success of mitral valve repair in functional mitral regurgitation. Understanding the complex interplay between annular-ring stresses and left ventricular function has significant implications for patient-ring selection, repair failure, and patient safety. A step towards assessing these challenges is developing a transducer that can be implanted in the exact method as commercially available rings and can quantify multidirectional ring loading. An annuloplasty ring transducer was developed to measure stresses at eight locations on both the in-plane and out-of-plane surfaces of an annuloplasty ring's titanium core. The transducer was implanted in an ovine subject using 10 sutures at near symmetric locations. At implantation, the ring was observed to undersize the mitral annulus. The flaccid annulus exerted both compressive (−) and tensile stresses (+) on the ring ranging from −3.17 to 5.34 MPa. At baseline hemodynamics, stresses cyclically changed and peaked near mid-systole. Mean changes in cyclic stress from ventricular diastole to mid-systole ranged from −0.61 to 0.46 MPa (in-plane direction) and from −0.49 to 1.13 MPa (out-of-plane direction). Results demonstrate the variability in ring stresses that can be introduced during implantation and the cyclic contraction of the mitral annulus. Ring stresses at implantation were approximately 4 magnitudes larger than the cyclic changes in stress throughout the cardiac cycle. These methods will be extended to ring transducers of differing size and geometry. Upon additional investigation, these data will contribute to improved knowledge of annulus-ring stresses, LV function, and the safer development of mitral repair techniques.     

Regional passive ventricular stress-strain relations during development of altered loads in chick embryo

Mechanical load influences embryonic ventricular growth, morphogenesis, and function. However, little is known about changes in regional passive ventricular properties during the development of altered mechanical loading conditions in the embryo. We tested the hypothesis that regional mechanical loads are a critical determinant of embryonic ventricular passive properties. We measured biaxial passive right and left ventricular (RV and LV, respectively) stress-strain relations in chick embryos at Hamburger-Hamilton stages 21 and 27 after conotruncal banding (CTB) to increase biventricular pressure load or left atrial ligation (LAL) to reduce LV volume load and increase RV volume load. In the RV, wall strains at end-diastolic (ED) pressure normalized whereas ED stresses increased after either CTB or LAL during development. In the left ventricle, both ED strain and stress normalized after CTB, whereas both remained reduced with significantly increased myocardial stiffness after LAL. These results suggest that the embryonic ventricle adapts to chronically altered mechanical loading conditions by changing specific RV and LV passive properties. Thus regional mechanical load has a critical role during cardiogenesis.     

Origin of axial prestretch and residual stress in arteries

The structural protein elastin endows large arteries with unique biological functionality and mechanical integrity, hence its disorganization, fragmentation, or degradation can have important consequences on the progression and treatment of vascular diseases. There is, therefore, a need in arterial mechanics to move from materially uniform, phenomenological, constitutive relations for the wall to those that account for separate contributions of the primary structural constituents: elastin, fibrillar collagens, smooth muscle, and amorphous matrix. In this paper, we employ a recently proposed constrained mixture model of the arterial wall and show that prestretched elastin contributes significantly to both the retraction of arteries that is observed upon transection and the opening angle that follows the introduction of a radial cut in an unloaded segment. We also show that the transmural distributions of elastin and collagen, compressive stiffness of collagen, and smooth muscle tone play complementary roles. Axial prestresses and residual stresses in arteries contribute to the homeostatic state of stress in vivo as well as adaptations to perturbed loads, disease, or injury. Understanding better the development of and changes in wall stress due to individual extracellular matrix constituents thus promises to provide considerable clinically important insight into arterial health and disease.     

Epicardial suction: a new approach to mechanical testing of the passive ventricular wall

The lack of an appropriate three-dimensional constitutive relation for stress in passive ventricular myocardium currently limits the utility of existing mathematical models for experimental and clinical applications. Previous experiments used to estimate parameters in three-dimensional constitutive relations, such as biaxial testing of excised myocardial sheets or passive inflation of the isolated arrested heart, have not included significant transverse shear deformation or in-plane compression. Therefore, a new approach has been developed in which suction is applied locally to the ventricular epicardium to introduce a complex deformation in the region of interest, with transmural variations in the magnitude and sign of nearly all six strain components. The resulting deformation is measured throughout the region of interest using magnetic resonance tagging. A nonlinear, three-dimensional, finite element model is used to predict these measurements at several suction pressures. Parameters defining the material properties of this model are optimized by comparing the measured and predicted myocardial deformations. We used this technique to estimate material parameters of the intact passive canine left ventricular free wall using an exponential, transversely isotropic constitutive relation. We tested two possible models of the heart wall: first, that it was homogeneous myocardium, and second, that the myocardium was covered with a thin epicardium with different material properties. For both models, in agreement with previous studies, we found that myocardium was nonlinear and anisotropic with greater stiffness in the fiber direction. We obtained closer agreement to previously published strain data from passive filling when the ventricular wall was modeled as having a separate, isotropic epicardium. These results suggest that epicardium may play a significant role in passive ventricular mechanics.     

Hemodynamic effects of cardiac cycle-specific increases in intrathoracic pressure

Changes in intrathoracic pressure (ITP) can influence cardiac performance by affecting ventricular loading conditions. Because both systemic venous return and factors determining left ventricular (LV) ejection may vary over the cardiac cycle, phasic increases in ITP may differentially affect preload or afterload if delivered at specific points within the cardiac cycle. We studied the hemodynamic effects of cardiac cycle-specific increases in ITP (pulses) delivered by a high-frequency jet ventilator in an acute closed-chested canine model (n = 11), using electromagnetic flow probes to measure biventricular stroke volume. Measurements were taken during a control condition after the induction of acute ventricular failure (AVF) by propranolol hydrochloride and volume infusion. ITP was independently varied without changing lung volume by the inflation of thoracoabdominal binders. Although synchronous pulses had minimal hemodynamic effects in unbound controls, binding pulses timed to occur in early diastole resulted in decreases in LV filling pressure and left ventricular stroke volume (SVlv) (P less than 0.05). In the AVF condition, pulses increased LV performance, evidenced by increases in SVlv (P less than 0.01), despite decreases in LV filling pressure (P less than 0.05). This effect is maximized by binding and by timing the pulses to occur in systole. We conclude that cardiac cycle-specific increases in ITP can significantly affect cardiac performance. These effects appear to be related to the ability of such timed pulses to selectively affect LV preload and afterload.