Dependence of local left ventricular wall mechanics on myocardial fiber orientation: a model study.

The dependence of local left ventricular (LV) mechanics on myocardial muscle fiber orientation was investigated using a finite element model. In the model we have considered anisotropy of the active and passive components of myocardial tissue, dependence of active stress on time, strain and strain rate, activation sequence of the LV wall and aortic afterload. Muscle fiber orientation in the LV wall is quantified by the helix fiber angle, defined as the angle between the muscle fiber direction and the local circumferential direction. In a first simulation, a transmural variation of the helix fiber angle from +60 degrees at the endocardium through 0 degrees in the midwall layers to -60 degrees at the epicardium was assumed. In this simulation, at the equatorial level maximum active muscle fiber stress was found to vary from about 110 kPa in the subendocardial layers through about 30 kPa in the midwall layers to about 40 kPa in the subepicardial layers. Next, in a series of simulations, muscle fiber orientation was iteratively adapted until the spatial distribution of active muscle fiber stress was fairly homogeneous. Using a transmural course of the helix fiber angle of +60 degrees at the endocardium, +15 degrees in the midwall layers and -60 degrees at the epicardium, at the equatorial level maximum active muscle fiber stress varied from 52 kPa to 55 kPa, indicating a remarkable reduction of the stress range. Moreover, the change of muscle fiber strain with time was more similar in different parts of the LV wall than in the first simulation. It is concluded that (1) the distribution of active muscle fiber stress and muscle fiber strain across the LV wall is very sensitive to the transmural distribution of the helix fiber angle and (2) a physiological transmural distribution of the helix fiber angle can be found, at which active muscle fiber stress and muscle fiber strain are distributed approximately homogeneously across the LV wall.     

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.     

Distribution of myocardial stress and its influence on coronary blood flow

The myocardial stress was analyzed by biomechanical modeling in correlation with experimental findings. The pressure-volume relationship follows the stress-strain relationship of muscle fibers. From the knowledge of fiber orientation and the distribution of sarcomere length, the myocardial stress components including fiber, longitudinal, circumferential and radial stresses were expressed as a function of fraction of wall thickness. The coronary blood flow is influenced by the myocardial radial stress. With the use of vascular waterfall theory, it is possible to correlate the theoretically defined stress distribution with experimentally obtained stress distribution. An elevation of radial stress in myocardium causes a reduction of vessel patency. During both systole and diastole, vessel patency remains constant at epicardium. At endocardium, however, vessel patency undergoes rhythmic changes following the systolic and diastolic influences of the radial stress. The physiological implication is that during systole, the endocardium suffers low blood flow and this transient ischemic state requires compensatory replenishment from diastolic perfusion. Such phenomena become less apparent toward the epicardium.     

Direct measurement of transmural laminar architecture in the anterolateral wall of the ovine left ventricle: new implications for wall thickening mechanics

Laminar, or sheet, architecture of the left ventricle (LV) is a structural basis for normal systolic and diastolic LV dynamics, but transmural sheet orientations remain incompletely characterized. We directly measured the transmural distribution of sheet angles in the ovine anterolateral LV wall. Ten Dorsett-hybrid sheep hearts were perfusion fixed in situ with 5% buffered glutaraldehyde at end diastole and stored in 10% formalin. Transmural blocks of myocardial tissue were excised, with the edges cut parallel to local circumferential, longitudinal, and radial axes, and sliced into 1-mm-thick sections parallel to the epicardial tangent plane from epicardium to endocardium. Mean fiber directions were determined in each section from five measurements of fiber angles. Each section was then cut transverse to the fiber direction, and five sheet angles (beta) were measured and averaged. Mean fiber angles progressed nearly linearly from -41 degrees (SD 11) at the epicardium to +42 degrees (SD 16) at the endocardium. Two families of sheets were identified at approximately +45 degrees (beta(+)) and -45 degrees (beta(-)). In the lateral region (n = 5), near the epicardium, sheets belonged to the beta(+) family; in the midwall, to the beta(-) family; and near the endocardium, to the beta(+) family. This pattern was reversed in the basal anterior region (n = 4). Sheets were uniformly beta(-) over the anterior papillary muscle (n = 2). These direct measurements of sheet angles reveal, for the first time, alternating transmural families of predominant sheet angles. This may have important implications in understanding wall mechanics in the normal and the failing heart.     

Some implications of a constant fiber stress hypothesis in the diastolic left ventricle

A thick-wall incompressible, elastic sphere was used as a model for the diastolic rat left ventricle. A model for myocardial nonhomogeneity was derived assuming that fiber (circumferential) stress was independent of position in the ventricular wall. The theoretical implications of the resulting constitutive relations together with the spherical model were analyzed in the context of large deformation elasticity theory. It was found that muscle stiffness at a given level of uniaxial stress increased monotonically from the endocardium to the epicardium. In addition, fiber stress was found to be essentially a linear function of transmural pressure above a pressure of 6 g/cm². It was also shown theoretically that neglecting the nonhomogeneity of the myocardium resulted in a state of stress which differed significantly from that predicted by the nonhomogeneous model. For example, at a transmural pressure of 14 g/cm², fiber stress in the nonhomogenous model was equal to 17 g/cm² while fiber stress in the homogeneous model varied between 100 g/cm² at the endocardial surface and 2 g/cm² at the epicardial surface. The change in muscle stiffness with position which characterized the nonhomogeneous model also tended to linearize the highly curvilinear radial stress distribution predicted by the homogeneous model at a given transmural pressure.     

On a nonlinear theory for muscle shells: Part II--Application to the beating left ventricle

This paper specializes the nonlinear laminated-muscle-shell theory developed in Part I to cylindrical geometry and computes stresses in arteries and the beating left ventricle. The theory accounts for large strain, material nonlinearity, thick-shell effects, torsion, muscle activation, and residual strain. First, comparison with elasticity solutions for pressurized arteries shows that the accuracy of the shell theory increases as transmural stress gradients and the shell thickness decrease. Residual strain reduces the stress gradients, lowering the error in the predicted peak stress in thick-walled arteries (R/t = 2.8) from about 30 to 10 percent. Second, the canine left ventricle is modeled as a thick-walled laminated cylinder with an internal pressure. Each layer is composed of transversely isotropic muscle with a fiber orientation based on anatomical data. Using a single pseudostrain-energy density function (with time-varying coefficients) for passive and active myocardium, the model predicts strain distributions that agree fairly well with published experimental measurements. The results also show that the peak fiber stress occurs subendocardially near the beginning of ejection and that residual strains significantly alter stress gradients within each lamina, but the magnitude of the peak fiber stress changes by less than 20 percent.     

Theoretical models in mechanics of the left ventricle

Different rheological concepts and theoretical studies have been recently presented using models of myocardial mechanics. Complex analysis of the mechanical behavior of the left ventricular wall have been developed in order to estimate the local stresses and deformations that occur during the heart cycle as well as the ventricular stroke volume and pressure. Theoretical models have taken into account non-linear and viscoelastic passive properties of the myocardium tissue, when subjected to large deformations, through given strain energy functions or stress-strain relations. Different prolate spheroid geometries have been considered for such thick shell cardiac structure. During the active state of the contraction, the rheological behavior of the fibers has been described using different muscle models and relationships between fiber tension and strain, and activation degree. A forthcoming approach for bridging the gap between the knowledge of the muscle fiber microrheological properties and the study of the mechanical behavior of the entire ventricle, consists in including anisotropic and inhomogeneous effects through fiber direction field.     

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.     

Coupling field theory with continuum mechanics: a simulation of domain formation in giant unilamellar vesicles

Domain formation is modeled on the surface of giant unilamellar vesicles using a Landau field theory model for phase coexistence coupled to elastic deformation mechanics (e.g., membrane curvature). Smooth particle applied mechanics, a form of smoothed particle continuum mechanics, is used to solve either the time-dependent Landau-Ginzburg or Cahn-Hilliard free-energy models for the composition dynamics. At the same time, the underlying elastic membrane is modeled using smooth particle applied mechanics, resulting in a unified computational scheme capable of treating the response of the composition fields to arbitrary deformations of the vesicle and vice versa. The results indicate that curvature coupling, along with the field theory model for composition free energy, gives domain formations that are correlated with surface defects on the vesicle. In the case that external deformations are included, the domain structures are seen to respond to such deformations. The present simulation capability provides a significant step forward toward the simulation of realistic cellular membrane processes.     

Cardiac mechanics in the stage-16 chick embryo.

A theoretical model is presented for the tubular heart of the stage-16 chick embryo (2.3 days of a 21-day incubation period). The model is a thick-walled, pseudoelastic cylindrical shell composed of three isotropic layers: the endocardium, the cardiac jelly, and the myocardium. The analysis is based on a shell theory that accounts for large deformation, material nonlinearity, residual strain, and muscle activation, with material properties inferred from available experimental data. We also measured epicardial strains from recorded motions of microspheres on the primitive right ventricles of stage-16 white Leghorn chick embryos. Relative to end diastole, peak axial and circumferential Lagrange strains occurred near end systole and had similar values. The magnitudes of these strains varied along the longitudinal axis of the heart (-0.16 +/- 0.08), being larger near the ends of the primitive right ventricle and smaller near midventricle. The in-plane shear strain was less than 0.05. Comparison of theoretical and experimental strains during the cardiac cycle shows generally good agreement. In addition, the model gives strong stress concentrations in the myocardial layer at end systole.     

Active and passive stresses in the myocardium

A mathematical approach that can be used to calculate the passive stress in the ventricular wall is presented. The active fiber stress (force/unit area) generated by the muscular fibers in the ventricular wall is expressed by means of body force (force/unit volume of the myocardium). It is shown that the total intramyocardial passive stress induced in the passive medium of the myocardium can be expressed as the sum of a passive stress induced by the left ventricular pressure and a passive stress induced by the active fiber stress. Applications to experimental data published in the literature are given. New results are presented that show the relation among those two components of the intramyocardial passive stress. New relations between the intramyocardial passive stress, the slope (elastance) of the pressure-volume relation, and the residual volume are also derived. The results obtained give a better understanding of some aspects of the mechanics of cardiac contraction and can provide a more detailed interpretation of clinical conditions.     

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.     

Optimizing ventricular fibers: uniform strain or stress,but not ATP consumption,leads to high efficiency

The aim of this study was to investigate the influence of fiber orientation in the left ventricular (LV) wall on the ejection fraction, efficiency, and heterogeneity of the distributions of developed fiber stress, strain and ATP consumption. A finite element model of LV mechanics was used with active properties of the cardiac muscle described by the Huxley-type cross-bridge model. The computed variances of sarcomere length (SL(var)), developed stress (DS(var)), and ATP consumption (ATP(var)) have several minima at different transmural courses of helix fiber angle. We identified only one region in the used design space with high ejection fraction, high efficiency of the LV and relatively small SL(var), DS(var), and ATP(var). This region corresponds to the physiological distribution of the helix fiber angle in the LV wall. Transmural fiber angle can be predicted by minimizing SL(var) and DS(var), but not ATP(var). If ATP(var) was minimized, then the transverse fiber angle was considerably underestimated. The results suggest that ATP consumption distribution is not regulating the fiber orientation in the heart.     

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.     

Structural and mechanical factors influencing infarct scar collagen organization

Although large collagen fibers in myocardial infarct scar are highly organized, little is known about mechanisms controlling this organization. The preexisting extracellular matrix may act as a scaffold along which fibroblasts migrate. Conversely, deformation within the ischemic area could guide fibroblasts so new collagen is oriented to counteract the stretch. To investigate these potential mechanisms, we infarcted three groups of pigs. Group 1 served as infarct controls. Group 2 had the endocardium slit longitudinally to alter local systolic deformation. Group 3 had a plug sectioned from ischemic tissue and rotated 90 degrees. The slit altered systolic deformation in the infarcted tissue, changing circumferential strain from expansion to compression and increasing radial strain and shears and the variability of collagen fiber angles but not the mean angle. In the plug pigs, when deformation, matrix orientation, and continuity are altered in the infarct area, the result is complete disarray in the organization of collagen within the infarct scar.     

MRI-based finite-element analysis of left ventricular aneurysm

Tagged MRI and finite-element (FE) analysis are valuable tools in analyzing cardiac mechanics. To determine systolic material parameters in three-dimensional stress-strain relationships, we used tagged MRI to validate FE models of left ventricular (LV) aneurysm. Five sheep underwent anteroapical myocardial infarction (25% of LV mass) and 22 wk later underwent tagged MRI. Asymmetric FE models of the LV were formed to in vivo geometry from MRI and included aneurysm material properties measured with biaxial stretching, LV pressure measurements, and myofiber helix angles measured with diffusion tensor MRI. Systolic material parameters were determined that enabled FE models to reproduce midwall, systolic myocardial strains from tagged MRI (630 +/- 187 strain comparisons/animal). When contractile stress equal to 40% of the myofiber stress was added transverse to the muscle fiber, myocardial strain agreement improved by 27% between FE model predictions and experimental measurements (RMS error decreased from 0.074 +/- 0.016 to 0.054 +/- 0.011, P < 0.05). In infarct border zone (BZ), end-systolic midwall stress was elevated in both fiber (24.2 +/- 2.7 to 29.9 +/- 2.4 kPa, P < 0.01) and cross-fiber (5.5 +/- 0.7 to 11.7 +/- 1.3 kPa, P = 0.02) directions relative to noninfarct regions. Contrary to previous hypotheses but consistent with biaxial stretching experiments, active cross-fiber stress development is an integral part of LV systole; FE analysis with only uniaxial contracting stress is insufficient. Stress calculations from these validated models show 24% increase in fiber stress and 115% increase in cross-fiber stress at the BZ relative to remote regions, which may contribute to LV remodeling.     

Sarcomere length dispersion in single skeletal muscle fibers and fiber bundles.

Light diffraction patterns produced by single skeletal muscle fibers and small fiber bundles of Rana pipiens semitendinosus have been examined at rest and during tetanic contraction. The muscle diffraction patterns were recorded with a vidicon camera interfaced to a minicomputer. Digitized video output was analyzed on-line to determine mean sarcomere length, line intensity, and the distribution of sarcomere lengths. The occurrence of first-order line intensity and peak amplitude maxima at approximately 3.0 mum is interpreted in terms of simple scattering theory. Measurements made along the length of a singel fiber reveal small variations in calculated mean sarcomere length (SD about 1.2%) and its percent dispersion (2.1% +/- 0.8%). Dispersion in small multifiber preparations increases approximately linearly with fiber number (about 0.2% per fiber) to a maximum of 8-10% in large bundles. Dispersion measurements based upon diffraction line analysis are comparable to SDs calculated from length distribution histograms obtained by light micrography of the fiber. First-order line intensity decreases by about 40% during tetanus; larger multifibered bundles exhibit substantial increases in sarcomere dispersion during contraction, but single fibers show no appreciable dispersion change. These results suggest the occurrence of asynchronous static or dynamic axial disordering of thick filaments, with a persistence in long range order of sarcomere spacing during contraction in single fibers.     

Wood specific gravity and anatomy in Heliocarpus appendiculatus (Tiliaceae)

Specific gravity exhibits extremely large radial increases with distance from the pith in Heliocarpus appendiculatus Turcz. (Tiliaceae), a pioneer of neotropical wet forests. To determine some of the wood anatomical changes associated with this increase, wood samples taken at breast height from three trees were divided into 1.0-cm-long segments from pith to bark. Measurements were made of fiber wall thickness, fiber lumen diameter, and percentages of fibers, axial parenchyma, ray parenchyma, and vessels on sections prepared from each segment. The extreme radial increases in specific gravity were associated with increases in fiber wall thickness, decreases in fiber diameter, decreases in fiber lumen diameter, and changes in the relative proportions of fibers and parenchyma. The increase in percent fiber concomitant with a decrease in axial parenchyma was the most important contributor to the increase in specific gravity in this species. The best predictor of specific gravity was percent fibers (r = 0.91, 0.92, 0.94) or percent axial parenchyma (r = -0.92, -0.91, -0.95), two variables that were highly intercorrelated (r = -0.95, -0.98, -0.99).     

Human annulus fibrosus material properties from biaxial testing and constitutive modeling are altered with degeneration

The annulus fibrosus (AF) of the intervertebral disk undergoes large and multidirectional stresses and strains. Uniaxial tensile tests are limited for measuring AF material properties, because freely contracting edges can prevent fiber stretch and are not representative of in situ boundary conditions. The objectives of this study were to measure human AF biaxial tensile mechanics and to apply and validate a constitutive model to determine material properties. Biaxial tensile tests were performed on samples oriented along the circumferential–axial and the radial–axial directions. Data were fit to a structurally motivated anisotropic hyperelastic model composed of isotropic extra-fibrillar matrix, nonlinear fibers, and fiber–matrix interactions (FMI) normal to the fibers. The validated model was used to simulate shear and uniaxial tensile behavior, to investigate AF structure–function, and to quantify the effect of degeneration. The biaxial stress–strain response was described well by the model (R ²?>?0.9). The model showed that the parameters for fiber nonlinearity and the normal FMI correlated with degeneration, resulting in an elongated toe-region and lower stiffness with degeneration. The model simulations in shear and uniaxial tension successfully matched previously published circumferential direction Young’s modulus, provided an explanation for the low values in previously published axial direction Young’s modulus, and was able to simulate shear mechanics. The normal FMI were important contributors to stress and changed with degeneration, therefore, their microstructural and compositional source should be investigated. Finally, the biaxial mechanical data and constitutive model can be incorporated into a disk finite element model to provide improved quantification of disk mechanics.