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.     

A finite element study of invariant-based orthotropic constitutive equations in the context of myocardial material parameter estimation

A previous study investigated a number of invariant-based orthotropic and transversely isotropic constitutive equations for their suitability to fit three-dimensional simple shear mechanics data of passive myocardial tissue. The study was based on the assumption of a homogeneous deformation. Here, we extend the previous study by performing an inverse finite element material parameter estimation. This ensures a more realistic deformation state and material parameter estimates. The constitutive relations were compared on the basis of (i) ‘goodness of fit’: how well they fit a set of six shear deformation tests and (ii) ‘variability’: how well determined the material parameters are over the range of experiments. These criteria were utilised to discuss the advantages and disadvantages of the constitutive relations. It was found that a specific form of the polyconvex type as well as the exponential Fung-type equations were most suitable for modelling the orthotropic behaviour of myocardium under simple shear.     

Myocardial material parameter estimation

The passive material properties of myocardium play a major role in diastolic performance of the heart. In particular, the shear behaviour is thought to play an important mechanical role due to the laminar architecture of myocardium. We have previously compared a number of myocardial constitutive relations with the aim to extract their suitability for inverse material parameter estimation. The previous study assumed a homogeneous deformation. In the present study we relaxed the homogeneous assumption by implementing these laws into a finite element environment in order to obtain more realistic measures for the suitability of these laws in both their ability to fit a given set of experimental data, as well as their stability in the finite element environment. In particular, we examined five constitutive laws and compare them on the basis of (i) "goodness of fit": how well they fit a set of six shear deformation tests, (ii) "determinability": how well determined the objective function is at the optimal parameter fit, and (iii) "variability": how well determined the material parameters are over the range of experiments. Furthermore, we compared the FE results with those from the previous study.It was found that the same material law as in the previous study, the orthotropic Fung-type "Costa-Law", was the most suitable for inverse material parameter estimation for myocardium in simple shear.     

Myocardial material parameter estimation-a comparative study for simple shear

The study of ventricular mechanics-analyzing the distribution of strain and stress in myocardium throughout the cardiac cycle-is crucially dependent on the accuracy of the constitutive law chosen to represent the highly nonlinear and anisotropic properties of passive cardiac muscle. A number of such laws have been proposed and fitted to experimental measurements of stress-strain behavior. Here we examine five of these laws and compare them on the basis of (i) "goodness of fit:" How well they fit a set of six shear deformation tests, (ii) "determinability:" How well determined the objective function is at the optimal parameter fit, and (iii) "variability:" How well determined the material parameters are over the range of experiments. These criteria are utilized to discuss the advantages and disadvantages of the constitutive laws.     

A triaxial-measurement shear-test device for soft biological tissues

A novel shear-test device for soft biological tissue, capable of applying simple shear deformations simultaneously in two orthogonal directions while measuring the resulting forces generated in three axes, is described. We validated the device using a synthetic gel, the properties of which were ascertained from independent tensile and rotational shear tests. Material parameters for the gel were fitted using neo-Hookean analytical solutions to the independent test data, and these matched the results from the device. Preliminary results obtained with rat septal myocardium are also presented to demonstrate the feasibility of the apparatus in determining the shear characteristics of living tissue.     

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.     

Myocardial material parameter estimation: a comparison of invariant based orthotropic constitutive equations

This study investigated a number of invariant based orthotropic and transversely isotropic constitutive equations for their suitability to fit three-dimensional simple shear mechanics data of passive myocardial tissue.

A number of orthotropic laws based on Green strain components and one microstructurally based law have previously been investigated to fit experimental measurements of stress-strain behaviour. Here we extend this investigation to include several recently proposed functional forms, i.e. invariant based orthotropic and transversely isotropic constitutive relations.

These laws were compared on the basis of (i) ‘goodness of fit’: how well they fit a set of six shear deformation tests, (ii) ‘variability’: how well determined the material parameters are over the range of experiments. These criteria were utilised to discuss the advantages and disadvantages of the constitutive laws.

It was found that a specific form of the polyconvex type as well as the exponential Fung-type law from the previous study were most suitable for modelling the orthotropic behaviour of myocardium under simple shear.     

Trabecular bone adaptation with an orthotropic material model.

Most bone adaptation algorithms, that attempt to explain the connection between bone morphology and loads, assume that bone is effectively isotropic. An isotropic material model can explain the bone density distribution, but not the structure and pattern of trabecular bone, which clearly has a mechanical significance. In this paper, an orthotropic material model is utilized to predict the proximal femur trabecular structure. Two hypotheses are combined to determine the local orientation and material properties of each element in the model. First, it is suggested that trabecular directions, which correspond to the orthotropic material axes, are determined locally by the maximal principal stress directions due to the multiple load cases (MLC) the femur is subject to. The second hypothesis is that material properties in each material direction can be determined using directional stimuli, thus extending existing adaptation algorithms to include directionality. An algorithm is utilized, where each iteration comprises of two stages. First, material axes are rotated to the direction of the largest principal stress that occurs from a multiple load scheme applied to the proximal femur. Next, material properties are modified in each material direction, according to a directional stimulus. Results show that local material directions correspond with known trabecular patterns, reproducing all main groups of trabeculae very well. The local directional stiffnesses, degree of anisotropy and density distribution are shown to conform to real femur morphology.     

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.     

Optical acquisition and polar decomposition of the full-field deformation gradient tensor within a fracture callus

Tracking tissue deformation is often hampered by material inhomogeneity, so local measurements tend to be insufficient thus lending to the necessity of full-field optical measurements. This study presents a novel approach to factoring heterogeneous deformation of soft and hard tissues in a fracture callus by introducing an anisotropic metric derived from the deformation gradient tensor (F). The deformation gradient tensor contains all the information available in a Green–Lagrange strain tensor, plus the rigid-body rotational components. A recent study [Bottlang et al., Journal of Biomechanics 41(3), 2008] produced full-field strains within ovine fracture calluses acquired through the application of electronic speckle pattern interferometery (ESPI). The technique is based on infinitesimal strain approximation (Engineering Strain) whose scheme is not independent of rigid-body rotation. In this work, for rotation extraction, the stretch and rotation tensors were separately determined from F by the polar decomposition theorem. Interfragmentary motions in a fracture gap were characterized by the two distinct mechanical factors (stretch and rotation) at each material point through full-field mapping. In the composite nature of bone and soft tissue, collagen arrangements are hypothesized such that fibers locally aligned with principal directions will stretch and fibers not aligned with the principal direction will rotate and stretch. This approach has revealed the deformation gradient tensor as an appropriate quantification of strain within callus bony and fibrous tissue via optical measurements.     

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.     

Physical invariant strain energy function for passive myocardium

Principal axis formulations are regularly used in isotropic elasticity, but they are not often used in dealing with anisotropic problems. In this paper, based on a principal axis technique, we develop a physical invariant orthotropic constitutive equation for incompressible solids, where it contains only a one variable (general) function. The corresponding strain energy function depends on six invariants that have immediate physical interpretation. These invariants are useful in facilitating an experiment to obtain a specific constitutive equation for a particular type of materials. The explicit appearance of the classical ground-state constants in the constitutive equation simplifies the calculation for their admissible values. A specific constitutive model is proposed for passive myocardium, and the model fits reasonably well with existing simple shear and biaxial experimental data. It is also able to predict a set of data from a simple shear experiment.     

Probing Cell Structure Responses Through a Shear and Stretching Mechanical Stimulation Technique

Cells are complex, dynamic systems that respond to various in vivo stimuli including chemical, mechanical, and scaffolding alterations. The influence of mechanics on cells is especially important in physiological areas that dictate what modes of mechanics exist. Complex, multivariate physiological responses can result from multi-factorial, multi-mode mechanics, including tension, compression, or shear stresses. In this study, we present a novel device based on elastomeric materials that allowed us to stimulate NIH 3T3 fibroblasts through uniaxial strip stretching or shear fluid flow. Cell shape and structural response was observed using conventional approaches such as fluorescent microscopy. Cell orientation and actin cytoskeleton alignment along the direction of applied force were observed to occur after an initial 3 h time period for shear fluid flow and static uniaxial strip stretching experiments although these two directions of alignment were oriented orthogonal relative to each other. This response was then followed by an increasingly pronounced cell and actin cytoskeleton alignment parallel to the direction of force after 6, 12, and 24 h, with 85% of the cells aligned along the direction of force after 24 h. These results indicate that our novel device could be implemented to study the effects of multiple modes of mechanical stimulation on living cells while probing their structural response especially with respect to competing directions of alignment and orientation under these different modes of mechanical stimulation. We believe that this will be important in a diversity of fields including cell mechanotransduction, cell–material interactions, biophysics, and tissue engineering.     

Direct measurement of intervertebral disc maximum shear strain in six degrees of freedom: motions that place disc tissue at risk of injury

Human intervertebral disc specimens were tested to determine the regions of largest maximum shear strain (MSS) experienced by disc tissues in each of three principal displacements and three rotations, and to identify the physiological rotations and displacements that may place the disc at greatest risk for large tissue strains and injury. Tearing of disc annulus may be initiated by large interlamellar shear strains. Nine human lumbar discs were tagged with radiographic markers on the endplates, disc periphery and with a grid of wires in the mid-transverse plane and subjected to each of the six principal displacements and rotations. Stereo-radiographs were taken in each position and digitized for reconstruction of the three-dimensional position of each marker. Maximum tissue shear strains were calculated from relative marker displacements and normalized by the input displacement or rotation. Lateral shear, compression, and lateral bending were the motions that produced the mean (95% confidence interval) largest mean MSS of 9.6 (0.7)%/mm, 9.0 (0.5)%/mm, and 5.8 (1.6)%/ degrees , respectively, and which occurred in the posterior, posterolateral and lateral peripheral regions of the disc. After taking into account the reported maximum physiological range of motion for each degree of freedom, motions producing the highest physiological MSS were lateral bending (57.8 (16.2)%) and flexion (38.3 (3.3)%), followed by lateral shear (14.4 (1.1)%) and compression (12.6 (0.7)%).     

Material properties of the ovine mitral valve anterior leaflet in vivo from inverse finite element analysis

We measured leaflet displacements and used inverse finite-element analysis to define, for the first time, the material properties of mitral valve (MV) leaflets in vivo. Sixteen miniature radiopaque markers were sewn to the MV annulus, 16 to the anterior MV leaflet, and 1 on each papillary muscle tip in 17 sheep. Four-dimensional coordinates were obtained from biplane videofluoroscopic marker images (60 frames/s) during three complete cardiac cycles. A finite-element model of the anterior MV leaflet was developed using marker coordinates at the end of isovolumic relaxation (IVR; when the pressure difference across the valve is approximately 0), as the minimum stress reference state. Leaflet displacements were simulated during IVR using measured left ventricular and atrial pressures. The leaflet shear modulus (G(circ-rad)) and elastic moduli in both the commisure-commisure (E(circ)) and radial (E(rad)) directions were obtained using the method of feasible directions to minimize the difference between simulated and measured displacements. Group mean (+/-SD) values (17 animals, 3 heartbeats each, i.e., 51 cardiac cycles) were as follows: G(circ-rad) = 121 +/- 22 N/mm2, E(circ) = 43 +/- 18 N/mm2, and E(rad) = 11 +/- 3 N/mm2 (E(circ) > E(rad), P < 0.01). These values, much greater than those previously reported from in vitro studies, may result from activated neurally controlled contractile tissue within the leaflet that is inactive in excised tissues. This could have important implications, not only to our understanding of mitral valve physiology in the beating heart but for providing additional information to aid the development of more durable tissue-engineered bioprosthetic valves.     

Model-free methods of analyzing domain motions in proteins from simulation: A comparison of normal mode analysis and molecular dynamics simulation of lysozyme

Model-free methods are introduced to determine quantities pertaining to protein domain motions from normal mode analyses and molecular dynamics simulations. For the normal mode analysis, the methods are based on the assumption that in low frequency modes, domain motions can be well approximated by modes of motion external to the domains. To analyze the molecular dynamics trajectory, a principal component analysis tailored specifically to analyze interdomain motions is applied. A method based on the curl of the atomic displacements is described, which yields a sharp discrimination of domains, and which defines a unique interdomain screw-axis. Hinge axes are defined and classified as twist or closure axes depending on their direction. The methods have been tested on lysozyme. A remarkable correspondence was found between the first normal mode axis and the first principal mode axis, with both axes passing within 3 Å of the alpha-carbon atoms of residues 2, 39, and 56 of human lysozyme, and near the interdomain helix. The axes of the first modes are overwhelmingly closure axes. A lesser degree of correspondence is found for the second modes, but in both cases they are more twist axes than closure axes. Both analyses reveal that the interdomain connections allow only these two degrees of freedom, one more than provided by a pure mechanical hinge. Proteins 27:425–437, 1997. © 1997 Wiley-Liss, Inc.     

Investigating full-field deformation of planar soft tissue under simple-shear tests

Finite simple shear test characteristics like specimen geometry and boundary conditions could affect the deformation homogeneity during the test. In order to ensure that the parameters of constitutive equations obtained from finite simple shear tests are appropriate, the deformation homogeneity of the specimen during simple-shear test should be examined. The Fourier transform moiré method (FTM) was used to examine the deformation uniformity of a porcine skin specimen in a finite simple shear test. The effects of clamping prestrain (0.15 and 0.3 engineering strain) and specimen geometry (5x5, 5x3.75, and 5x2.5cm) were investigated. These effects include in-plane deformation altered by clamping prestrain, slippage between specimen and clamps, and out-of-plane deformation. The experimental results showed that the wide specimen had more severe deformation alteration by clamping prestrain and was easier to slip out of the clamps when the shear angle is large. Furthermore, in all test configurations, the out-of-plane deformation is significant when the shear angle is large, and a narrow specimen is prone to have out-of-plane deformation. This study may provide guidelines for the selection of specimen aspect ratio and clamping prestrain when studying the material response of soft tissues under simple-shear tests.     

Angular orientation of trabecular bone in the femoral head and its relationship to hip joint loads in leaping primates

The elastic properties and mechanical behavior of trabecular bone are largely determined by its three-dimensional (3D) fabric structure. Recent work demonstrating a correlation between the primary mechanical and material axes in trabecular bone specimens suggests that fabric orientation may be used to infer directional components of the material strength and, by extension, the hypothetical loading regime. Here we quantify the principal orientation of trabecular bone in the femoral head and relate these principal fabric directions to loading patterns during various locomotor behaviors. The proximal femora of a diverse sample of prosimians were scanned using a high-resolution X-ray computed tomography scanner with resolution of better than 50 mum. Spherical volumes of interest were defined within the femoral heads and the 3D fabric anisotropy was calculated using the mean intercept length and star volume distribution methods. In addition to differences in bone volume and anisotropy, significant differences were found in the spatial orientation of the principal trabecular axes depending on locomotor behavior. The principal orientations for leapers (Galago, Tarsius, Avahi) are relatively tightly clustered (alpha(95) confidence limit: 8.2; angular variance s: 18.2 degrees ) and oriented in a superoanterior direction, while those of nonleapers are more variable across a range of directions (alpha(95): 16.8; s: 42.0 degrees ). The mean principal directions are significantly different for leaping vs. nonleaping taxa. These results further suggest a relationship between bone microstructure in the hip joint and locomotor behavior and indicate a similarity of loading across leapers despite differences in kinematics and phylogeny.     

Determination of elastomeric foam parameters for simulations of complex loading

BACKGROUND: Finite element (FE) analysis has shown promise for the evaluation of elastomeric foam personal protection devices. Although appropriate representation of foam materials is necessary in order to obtain realistic simulation results, material definitions used in the literature vary widely and often fail to account for the multi-mode loading experienced by these devices. This study aims to provide a library of elastomeric foam material parameters that can be used in FE simulations of complex loading scenarios. METHOD OF APPROACH: Twelve foam materials used in footwear were tested in uni-axial compression, simple shear and volumetric compression. For each material, parameters for a common compressible hyperelastic material model used in FE analysis were determined using: (a) compression; (b) compression and shear data; and (c) data from all three tests. RESULTS: Material parameters and Drucker stability limits for the best fits are provided with their associated errors. The material model was able to reproduce deformation modes for which data was provided during parameter determination but was unable to predict behavior in other deformation modes. CONCLUSIONS: Simulation results were found to be highly dependent on the extent of the test data used to determine the parameters in the material definition. This finding calls into question the many published results of simulations of complex loading that use foam material parameters obtained from a single mode of testing. The library of foam parameters developed here presents associated errors in three deformation modes that should provide for a more informed selection of material parameters.     

Elastic characterization of transversely isotropic soft materials by dynamic shear and asymmetric indentation

The mechanical characterization of soft anisotropic materials is a fundamental challenge because of difficulties in applying mechanical loads to soft matter and the need to combine information from multiple tests. A method to characterize the linear elastic properties of transversely isotropic soft materials is proposed, based on the combination of dynamic shear testing (DST) and asymmetric indentation. The procedure was demonstrated by characterizing a nearly incompressible transversely isotropic soft material. A soft gel with controlled anisotropy was obtained by polymerizing a mixture of fibrinogen and thrombin solutions in a high field magnet (B?=?11.7 T); fibrils in the resulting gel were predominantly aligned parallel to the magnetic field. Aligned fibrin gels were subject to dynamic (20-40 Hz) shear deformation in two orthogonal directions. The shear storage modulus was 1.08?±?0. 42 kPa (mean?±?std. dev.) for shear in a plane parallel to the dominant fiber direction, and 0.58?±?0.21 kPa for shear in the plane of isotropy. Gels were indented by a rectangular tip of a large aspect ratio, aligned either parallel or perpendicular to the normal to the plane of transverse isotropy. Aligned fibrin gels appeared stiffer when indented with the long axis of a rectangular tip perpendicular to the dominant fiber direction. Three-dimensional numerical simulations of asymmetric indentation were used to determine the relationship between direction-dependent differences in indentation stiffness and material parameters. This approach enables the estimation of a complete set of parameters for an incompressible, transversely isotropic, linear elastic material.