Fiber stress profiles in the left ventricle of the heart during diastole: Effects of distension and hypertrophy

Pressure-volume and volume-dimensions relationships, obtained from excised dog left ventricles were used for calculating the stresses acting along the longitudinal axis of the individual myocardial fibers. The calculations were based on a set of empirical and theoretical equations. The pressure-volume relationship as well as the volume-dimensions relationships for the excised left ventricle were expressed in the form of empirical equations; the fiber orientation was written as a function of the fiber location within the left ventricular wall; finally, the fiber stress was determined by means of theoretically derived formulas. Simultaneous solutions for the fibers of a meridian cut through the left ventricular myocardial shell were obtained by means of a digital computer and presented in the form of diagrams. The results showed that at low degrees of distension of the left ventricle there are two zones of higher stresses at the equatorial area, one near the epicardium and one near the endocardium. As the distension proceeds under the effect of progressively increasing intraventricular pressure, these two zones become less well defined, whereas a new zone of higher stresses appears near the apex. At high degrees of distension, the ventricle assumes a more spherical shape and the equatorial zones of higher stresses are replaced by zones of lower stresses. Increase in the myocardial mass results in appearance of the equatorial lower stress zones at lower degrees of distension.     

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.     

Pulsed ultrasonic Doppler velocity measurements inside a left ventricular assist device

In this study we have employed a single channel, pulsed ultrasonic Doppler velocimeter to measure instantaneous velocity distributions within the pumping chamber of a ventricular assist device. Instantaneous velocities have been decomposed into periodic mean and turbulent fluctuating components from which estimates of Reynolds stresses within the chamber and mean shear stresses along the wall of the chamber have been obtained. A review of the complete data set indicates a maximum value of the mean wall shear stress of 25 dynes/cm2 and a maximum Reynolds stress of 212 dynes/cm2. These values are lower than those measured distal to aortic valve prostheses in vitro and are well below levels known to damage blood components. Core flow patterns, wall washing patterns and flow stagnation points are also revealed.     

Structural finite deformation model of the left ventricle during diastole and systole

A model of left ventricular function is developed based on morphological characteristics of the myocardial tissue. The passive response of the three-dimensional collagen network and the active contribution of the muscle fibers are integrated to yield the overall response of the left ventricle which is considered to be a thick wall cylinder. The deformation field and the distributions of stress and pressure are determined at each point in the cardiac cycle by numerically solving three equations of equilibrium. Simulated results in terms of the ventricular deformation during ejection and isovolumic cycles are shown to be in good qualitative agreement with experimental data. It is shown that the collagen network in the heart has considerable effect on the pressure-volume loops. The particular pattern of spatial orientation of the collagen determines the ventricular recoil properties in early diastole. The material properties (myocardial stiffness and contractility) are shown to affect both the pressure-volume loop and the deformation pattern of the ventricle. The results indicate that microstructural consideration offer a realistic representation of the left ventricle mechanics.     

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.     

Left ventricular wall stress compendium

Left ventricular (LV) wall stress has intrigued scientists and cardiologists since the time of Lame and Laplace in 1800s. The left ventricle is an intriguing organ structure, whose intrinsic design enables it to fill and contract. The development of wall stress is intriguing to cardiologists and biomedical engineers. The role of left ventricle wall stress in cardiac perfusion and pumping as well as in cardiac pathophysiology is a relatively unexplored phenomenon. But even for us to assess this role, we first need accurate determination of in vivo wall stress. However, at this point, 150 years after Lame estimated left ventricle wall stress using the elasticity theory, we are still in the exploratory stage of (i) developing left ventricle models that properly represent left ventricle anatomy and physiology and (ii) obtaining data on left ventricle dynamics. In this paper, we are responding to the need for a comprehensive survey of left ventricle wall stress models, their mechanics, stress computation and results. We have provided herein a compendium of major type of wall stress models: thin-wall models based on the Laplace law, thick-wall shell models, elasticity theory model, thick-wall large deformation models and finite element models. We have compared the mean stress values of these models as well as the variation of stress across the wall. All of the thin-wall and thick-wall shell models are based on idealised ellipsoidal and spherical geometries. However, the elasticity model's shape can vary through the cycle, to simulate the more ellipsoidal shape of the left ventricle in the systolic phase. The finite element models have more representative geometries, but are generally based on animal data, which limits their medical relevance. This paper can enable readers to obtain a comprehensive perspective of left ventricle wall stress models, of how to employ them to determine wall stresses, and be cognizant of the assumptions involved in the use of specific models.     

Left ventricular wall stress compendium

Left ventricular (LV) wall stress has intrigued scientists and cardiologists since the time of Lame and Laplace in 1800s. The left ventricle is an intriguing organ structure, whose intrinsic design enables it to fill and contract. The development of wall stress is intriguing to cardiologists and biomedical engineers. The role of left ventricle wall stress in cardiac perfusion and pumping as well as in cardiac pathophysiology is a relatively unexplored phenomenon. But even for us to assess this role, we first need accurate determination of in vivo wall stress. However, at this point, 150 years after Lame estimated left ventricle wall stress using the elasticity theory, we are still in the exploratory stage of (i) developing left ventricle models that properly represent left ventricle anatomy and physiology and (ii) obtaining data on left ventricle dynamics. In this paper, we are responding to the need for a comprehensive survey of left ventricle wall stress models, their mechanics, stress computation and results. We have provided herein a compendium of major type of wall stress models: thin-wall models based on the Laplace law, thick-wall shell models, elasticity theory model, thick-wall large deformation models and finite element models. We have compared the mean stress values of these models as well as the variation of stress across the wall. All of the thin-wall and thick-wall shell models are based on idealised ellipsoidal and spherical geometries. However, the elasticity model's shape can vary through the cycle, to simulate the more ellipsoidal shape of the left ventricle in the systolic phase. The finite element models have more representative geometries, but are generally based on animal data, which limits their medical relevance. This paper can enable readers to obtain a comprehensive perspective of left ventricle wall stress models, of how to employ them to determine wall stresses, and be cognizant of the assumptions involved in the use of specific models.     

A finite element model of the human left ventricular systole

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.     

Adaptation of passive rat left ventricle in diastolic dysfunction.

This article deals with providing a theoretical explanation for quantitative changes in the geometry, the opening angle and the deformation parameters of the rat ventricular wall during adaptation of the passive left ventricle in diastolic dysfunction. A large deformation theory is applied to analyse transmural stress and strain distribution in the left ventricular wall considering it to be made of homogeneous, incompressible, transversely isotropic, non-linear elastic material. The basic assumptions made for computing stress distributions are that the average circumferential stress and strain for the adaptive ventricle is equal to the average circumferential stress and strain in the normotensive ventricle, respectively.All the relevant parameters, such as opening angle, twist per unit length, axial extension, internal and external radii and others, in the stress-free, unloaded and loaded states of normotensive, hypertensive and adaptive left ventricle are determined. The circumferential stress and strain distribution through the ventricular wall are also computed. Our analysis predicts that during adaptation, wall thickness and wall mass of the ventricle increase. These results are consistent with experimental findings and are the indications of initiation of congestive heart failure.     

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.     

A finite element model of the human left ventricular systole

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.     

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.     

Left Ventricular Stresses in the Intact Human Heart

A set of stress differential equations of equilibrium is presented for a thick prolate spheroid which is the assumed shape for the left ventricle. An analysis for the stresses in the ventricular wall indicates that maximum stresses occur at the inner layers and decrease to a minimum at the epicardial surface, a result that is partially validated by experiment. Simple expressions are available for the evaluation of maximum stresses which occur at the equator and are suitable for small laboratory-oriented digital computers employed in the clinical evaluation of patient status. The surprising result is that Laplace's law yields practical values for mean stresses in thick-walled ventricles.     

The Continuity of Right and Left Ventricular Myocardial Sinusoidal Spaces and its Relation to Right Ventricular Implants

Evidence is presented which indicates that blood leaving side branches of an internal mammary artery implanted into the anterior wall of the right ventricle flows from the tunnel in which it lies through myocardial sinusoidal spaces of the anterior right ventricular wall across the midline to fill corresponding spaces in the anterior wall of the left ventricle and thence is carried to the left coronary sinus. The myocardial sinusoidal spaces of right and left ventricles have been well outlined, using injections of polyvinyl acetate and the technique of digestion casts. We have been able to show that there is no barrier between the myocardial sinusoids of the right circulation and those related to the anterior descending branch of the left coronary artery. In structure, these myocardial sinusoidal spaces are quite different from the intramyocardial coronary arteriolar zones which, in 93% of human hearts, are separated from one another without collateral communication.The continuity of the right and left ventricular myocardial sinusoids explains why implantation of a right internal mammary artery into the anterior wall of the right ventricle combined with a corresponding left implant, epicardiectomy and free omental graft, has been so effective in our hands in the treatment of far-advanced human coronary artery insufficiency.     

Mechanical analysis of the development of left ventricular aneurysms

The left ventricle is modelled as a spherical shell with an infarcted wall segment. The mechanics of the circumstances causing this infarcted segment to develop into an aneurysm is presented. Both the wall stresses and deformations are worked out for aneurysms developing from infarcts of different sizes and percentages of wall damage. The governing equations consist of incompressibility relations, force-equilibrium relations and stress-strain relations. Newton Raphson technique is used to solve these nonlinear simultaneous algebraic equations, for the values of the myocardial stresses in the infarcted segment and the bulge values, in terms of the ventricular geometry and the damage extent (expressed in terms of the damage angle and percentage of wall damage). The results indicate that in general it is innermost layer which is severely stressed and that in the rupture of the ventricle the critical factor involved is the percentage of infarct thickness rather than the angle of damage.     

A framework for biomechanics simulations using four-chamber cardiac models

Computational cardiac models have been extensively used to study different cardiac biomechanics; specifically, finite-element analysis has been one of the tools used to study the internal stresses and strains in the cardiac wall during the cardiac cycle. Cubic-Hermite finite element meshes have been used for simulating cardiac biomechanics due to their convergence characteristics and their ability to capture smooth geometries compactly–fewer elements are needed to build the cardiac geometry–compared to linear tetrahedral meshes. Such meshes have previously been used only with simple ventricular geometries with non-physiological boundary conditions due to challenges associated with creating cubic-Hermite meshes of the complex heart geometry. However, it is critical to accurately capture the different geometric characteristics of the heart and apply physiologically equivalent boundary conditions to replicate the in vivo heart motion. In this work, we created a four-chamber cardiac model utilizing cubic-Hermite elements and simulated a full cardiac cycle by coupling the 3D finite element model with a lumped circulation model. The myocardial fiber-orientations were interpolated within the mesh using the Log-Euclidean method to overcome the singularity associated with interpolation of orthogonal matrices. Physiologically equivalent rigid body constraints were applied to the nodes along the valve plane and the accuracy of the resulting simulations were validated using open source clinical data. We then simulated a complete cardiac cycle of a healthy heart and a heart with acute myocardial infarction. We compared the pumping functionality of the heart for both cases by calculating the ventricular work. We observed a 20% reduction in acute work done by the heart immediately after myocardial infarction. The myocardial wall displacements obtained from the four-chamber model are comparable to actual patient data, without requiring complicated non-physiological boundary conditions usually required in truncated ventricular heart models.     

A modified Holzapfel-Ogden law for a residually stressed finite strain model of the human left ventricle in diastole

In this work, we introduce a modified Holzapfel-Ogden hyperelastic constitutive model for ventricular myocardium that accounts for residual stresses, and we investigate the effects of residual stresses in diastole using a magnetic resonance imaging–derived model of the human left ventricle (LV). We adopt an invariant-based constitutive modelling approach and treat the left ventricular myocardium as a non-homogeneous, fibre-reinforced, incompressible material. Because in vivo images provide the configuration of the LV in a loaded state even in diastole, an inverse analysis is used to determine the corresponding unloaded reference configuration. The residual stress in this unloaded state is estimated by two different methods. One is based on three-dimensional strain measurements in a local region of the canine LV, and the other uses the opening angle method for a cylindrical tube. We find that including residual stress in the model changes the stress distributions across the myocardium and that whereas both methods yield qualitatively similar changes, there are quantitative differences between the two approaches. Although the effects of residual stresses are relatively small in diastole, the model can be extended to explore the full impact of residual stress on LV mechanical behaviour for the whole cardiac cycle as more experimental data become available. In addition, although not considered here, residual stresses may also play a larger role in models that account for tissue growth and remodelling.     

Velocity and turbulence measurements past mitral valve prostheses in a model left ventricle

Thrombogenesis and hemolysis have both been linked to the flow dynamics past heart valve prostheses. To learn more about the particular flow dynamics past mitral valve prostheses in the left ventricle under controlled experimental conditions, an in vitro study was performed. The experimental methods included velocity and turbulent shear stress measurements past caged-ball, tilting disc, bileaflet, and polyurethane trileaflet mitral valves in an acrylic rigid model of the left ventricle using laser Doppler anemometry. The results indicate that all four prosthetic heart valves studied create at least mildly disturbed flow fields. The effect of the left ventricular geometry on the flow development is to produce a stabilizing vortex which engulfs the entire left ventricular cavity, depending on the orientation of the valve. The measured turbulent shear stress magnitudes for all four valves did not exceed the reported value for hemolytic damage. However, the measured turbulent shear stresses were near or exceeded the critical shear stress reported in the literature for platelet lysis, a known precursor to thrombus formation.     

Non-homogeneous analysis of three-dimensional transmural finite deformation in canine ventricular myocardium

A new method has been developed for analyzing transmural distributions of finite deformation in canine ventricular myocardium without the need to assume that the strain in a finite volume of the wall is homogeneous. The three-dimensional nodal geometric parameters of bilinear-cubic or bilinear-quadratic finite elements are fitted by least squares to the measured coordinates of 12-18 radiopaque markers implanted in the left ventricular free wall. For six dog hearts, root-mean-squared errors in the fitted in-plane coordinates ranged from 0.079-0.556 mm in the end-diastolic reference state and 0.142-0.622 mm at end-systole. The corresponding error ranges in the radial coordinate were 0.042-0.264 mm at end-diastole and 0.106-0.279 mm at end-systole. Smoothly continuous transmural profiles of wall strain computed as the element deformed during the cardiac cycle from end-diastole to end-systole showed good agreement with the discrete results of conventional homogeneous analysis. Using the kinematics of a thick-walled incompressible cylinder, overall absolute errors due to the non-homogeneity of myocardial deformation were found to be reduced in the new analysis by 30-35% for typical experimental parameters. Overall relative errors were also reduced (from 23 to 20%). Since measurement errors in the reconstructed marker coordinates were spatially smoothed by the fitting procedure, noise in the computed deformations was also substantially attenuated, and transmural gradients of three-dimensional strain components could be obtained with improved accuracy. Hence physiological factors affected by transmural stress and strain distributions, such as myocardial blood flow, ischemia and hypertrophy, may be better understood.     

Heart size-independent analysis of myocardial function in murine pressure overload hypertrophy

Pressure overload cardiac hypertrophy may be a compensatory mechanism to normalize systolic wall stress and preserve left ventricular (LV) function. To test this concept, we developed a novel in vivo method to measure myocardial stress (sigma)-strain (epsilon) relations in normal and hypertrophied mice. LV volume was measured using two pairs of miniature omnidirectional piezoelectric crystals implanted orthogonally in the endocardium and one crystal placed on the anterior free wall to measure instantaneous wall thickness. Highly linear sigma-epsilon relations were obtained in control (n = 7) and hypertrophied mice produced by 7 days of transverse aortic constriction (TAC; n = 13). Administration of dobutamine in control mice significantly increased the load-independent measure of LV contractility, systolic myocardial stiffness. In TAC mice, systolic myocardial stiffness was significantly greater than in control mice (3,156 +/- 1,433 vs. 1,435 +/- 467 g/cm(2), P < 0.01), indicating enhanced myocardial contractility with pressure overload. However, despite the increased systolic performance, both active (time constant of LV pressure decay) and passive (diastolic myocardial stiffness constant) diastolic properties were markedly abnormal in TAC mice compared with control mice. These data suggest that the development of cardiac hypertrophy is associated with a heightened contractile state, perhaps as an early compensatory response to pressure overload.