Mechanical characterization of anisotropic planar biological soft tissues using large indentation: a computational feasibility study

Traditionally, the complex mechanical behavior of planar soft biological tissues is characterized by (multi)axial tensile testing. While uniaxial tests do not provide sufficient information for a full characterization of the material anisotropy, biaxial tensile tests are difficult to perform and tethering effects limit the analyses to a small central portion of the test sample. In both cases, determination of local mechanical properties is not trivial. Local mechanical characterization may be performed by indentation testing. Conventional indentation tests, however, often assume linear elastic and isotropic material properties, and therefore these tests are of limited use in characterizing the nonlinear, anisotropic material behavior typical for planar soft biological tissues. In this study, a spherical indentation experiment assuming large deformations is proposed. A finite element model of the aortic valve leaflet demonstrates that combining force and deformation gradient data, one single indentation test provides sufficient information to characterize the local material behavior. Parameter estimation is used to fit the computational model to simulated experimental data. The aortic valve leaflet is chosen as a typical example. However, the proposed method is expected to apply for the mechanical characterization of planar soft biological materials in general.     

An inverse finite element method for determining the anisotropic properties of the cornea

An inverse finite element method was developed to determine the anisotropic properties of bovine cornea from an in vitro inflation experiment. The experiment used digital image correlation (DIC) to measure the three-dimensional surface geometry and displacement field of the cornea at multiple pressures. A finite element model of a bovine cornea was developed using the DIC measured surface geometry of the undeformed specimen. The model was applied to determine five parameters of an anisotropic hyperelastic model that minimized the error between the measured and computed surface displacement field and to investigate the sensitivity of the measured bovine inflation response to variations in the anisotropic properties of the cornea. The results of the parameter optimization revealed that the collagen structure of bovine cornea exhibited a high degree of anisotropy in the limbus region, which agreed with recent histological findings, and a transversely isotropic central region. The parameter study showed that the bovine corneal response to the inflation experiment was sensitive to the shear modulus of the matrix at pressures below the intraocular pressure, the properties of the collagen lamella at higher pressures, and the degree of anisotropy in the limbus region. It was not sensitive to a weak collagen anisotropy in the central region.     

Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions

Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. However, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.     

Biaxial mechanical properties of the native and glutaraldehyde-treated aortic valve cusp: Part II--A structural constitutive model

We have formulated the first constitutive model to describe the complete measured planar biaxial stress-strain relationship of the native and glutaraldehyde-treated aortic valve cusp using a structurally guided approach. When applied to native, zero-pressure fixed, and low-pressure fixed cusps, only three parameters were needed to simulate fully the highly anisotropic, and nonlinear in-plane biaxial mechanical behavior. Differences in the behavior of the native and zero- and low-pressure fixed cusps were found to be primarily due to changes in the effective fiber stress-strain behavior. Further, the model was able to account for the effects of small (< 10 deg) misalignments in the cuspal specimens with respect to the biaxial test axes that increased the accuracy of the model material parameters. Although based upon a simplified cuspal structure, the model underscored the role of the angular orientation of the fibers that completely accounted for extreme mechanical anisotropy and pronounced axial coupling. Knowledge of the mechanics of the aortic cusp derived from this model may aid in the understanding of fatigue damage in bioprosthetic heart valves and, potentially, lay the groundwork for the design of tissue-engineered scaffolds for replacement heart valves.     

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.     

Mechanical characterisation of in vivo human skin using a 3D force-sensitive micro-robot and finite element analysis

The complex mechanical properties of skin have been the subject of much study in recent years. Several experimental methods developed to measure the mechanical properties of skin in vivo, such as suction or torsion, are unable to measure skin’s anisotropic characteristics. An experiment characterising the mechanical properties of in vivo human skin using a novel force-sensitive micro-robot is presented. The micro-robot applied in-plane deformations to the anterior forearm and the posterior upper arm. The behaviour of the skin in each area is highly nonlinear, anisotropic, and viscoelastic. The response of the upper arm skin is very dependent on the orientation of the arm. A finite element model consisting of an Ogden strain energy function and quasi-linear viscoelasticity was developed to simulate the experiments. An orthogonal initial stress field, representing the in vivo skin tension, was used as an additional model parameter. The model simulated the experiments accurately with an error-of-fit of 17.5% for the anterior lower forearm area, 6.5% for the anterior upper forearm and 9.3% for the posterior upper arm. The maximum in vivo tension in each area determined by the model was 6.2 Nm⁻¹ in the anterior lower forearm, 11.4 Nm⁻¹ in anterior upper forearm and 5.6 Nm⁻¹ in the posterior upper arm. The results also show that a finite element model with a neo-Hookean strain energy function cannot simulate the experiments with the same accuracy.     

Identification of regional mechanical anisotropy in soft tissue analogs

In a previous work (Raghupathy and Barocas, 2010, "Generalized Anisotropic Inverse Mechanics for Soft Tissues,"J. Biomech. Eng., 132(8), pp. 081006), a generalized anisotropic inverse mechanics method applicable to soft tissues was presented and tested against simulated data. Here we demonstrate the ability of the method to identify regional differences in anisotropy from full-field displacements and boundary forces obtained from biaxial extension tests on soft tissue analogs. Tissue heterogeneity was evaluated by partitioning the domain into homogeneous subdomains. Tests on elastomer samples demonstrated the performance of the method on isotropic materials with uniform and nonuniform properties. Tests on fibroblast-remodeled collagen cruciforms indicated a strong correlation between local structural anisotropy (measured by polarized light microscopy) and the evaluated local mechanical anisotropy. The results demonstrate the potential to quantify regional anisotropic material behavior on an intact tissue sample.     

A constitutive model for two-dimensional soft tissues and its application to experimental data

A constitutive relation proposed by Shoemaker (Ph.D. dissertation, 1984) to model the mechanical behavior of membraneous or two-dimensional soft tissues is described. Experiments by Schneider (Ph.D. dissertation, 1982) on human skin and Lee et al. (Am. J. Physiol., 249, H222-H230, 1985) on canine pericardium, and the application of the constitutive model to biaxial stress-strain data from these experiments, are discussed. Some experimental data and predictions of the model obtained by curvefitting are presented for comparison. Values of material parameters are also presented. It is concluded that the constitutive model is well able to fit results of individual tests, and that its generality (judged by consistency of parameters from test to test of the same specimen), though not complete, does compare favorably with some other results presented in the literature.     

A rheological network model for the continuum anisotropic and viscoelastic behavior of soft tissue

The mechanical behavior of soft tissue demonstrates a number of complex features including nonlinearity, anisotropy, viscoelasticity, and growth. Characteristic features of the time-dependent and anisotropic behavior are related to the properties of various components of the tissue such as fibrous collagen and elastin networks, large proteins and sugars attached to these networks, and interstitial fluid. Attempts to model the elastic behavior of these tissues based on assumptions about the behavior of the underlying constituents have been reasonably successful, but the essential addition of viscoelasticity to these models has been met with varying success. Here, a new rheological network model is proposed using, as its basis, an orthotropic hyperelastic constitutive model for fibrous tissue and a viscoelastic reptation model for soft materials. The resulting model has been incorporated into numerical and computational models, and is shown to capture the mechanical behavior of soft tissue in various modes of deformation including uniaxial and biaxial tension and simple shear.     

A microstructurally informed model for the mechanical response of three-dimensional actin networks

We propose a class of microstructurally informed models for the linear elastic mechanical behaviour of cross-linked polymer networks such as the actin cytoskeleton. Salient features of the models include the possibility to represent anisotropic mechanical behaviour resulting from anisotropic filament distributions, and a power law scaling of the mechanical properties with the filament density. Mechanical models within the class are parameterized by seven different constants. We demonstrate a procedure for determining these constants using finite element models of three-dimensional actin networks. Actin filaments and cross-links were modelled as elastic rods, and the networks were constructed at physiological volume fractions and at the scale of an image voxel. We show the performance of the model in estimating the mechanical behaviour of the networks over a wide range of filament densities and degrees of anisotropy.     

A microstructurally informed model for the mechanical response of three-dimensional actin networks

We propose a class of microstructurally informed models for the linear elastic mechanical behaviour of cross-linked polymer networks such as the actin cytoskeleton. Salient features of the models include the possibility to represent anisotropic mechanical behaviour resulting from anisotropic filament distributions, and a power law scaling of the mechanical properties with the filament density. Mechanical models within the class are parameterized by seven different constants. We demonstrate a procedure for determining these constants using finite element models of three-dimensional actin networks. Actin filaments and cross-links were modelled as elastic rods, and the networks were constructed at physiological volume fractions and at the scale of an image voxel. We show the performance of the model in estimating the mechanical behaviour of the networks over a wide range of filament densities and degrees of anisotropy.     

Collagen fiber alignment does not explain mechanical anisotropy in fibroblast populated collagen gels

Many load-bearing soft tissues exhibit mechanical anisotropy. In order to understand the behavior of natural tissues and to create tissue engineered replacements, quantitative relationships must be developed between the tissue structures and their mechanical behavior. We used a novel collagen gel system to test the hypothesis that collagen fiber alignment is the primary mechanism for the mechanical anisotropy we have reported in structurally anisotropic gels. Loading constraints applied during culture were used to control the structural organization of the collagen fibers of fibroblast populated collagen gels. Gels constrained uniaxially during culture developed fiber alignment and a high degree of mechanical anisotropy, while gels constrained biaxially remained isotropic with randomly distributed collagen fibers. We hypothesized that the mechanical anisotropy that developed in these gels was due primarily to collagen fiber orientation. We tested this hypothesis using two mathematical models that incorporated measured collagen fiber orientations: a structural continuum model that assumes affine fiber kinematics and a network model that allows for nonaffine fiber kinematics. Collagen fiber mechanical properties were determined by fitting biaxial mechanical test data from isotropic collagen gels. The fiber properties of each isotropic gel were then used to predict the biaxial mechanical behavior of paired anisotropic gels. Both models accurately described the isotropic collagen gel behavior. However, the structural continuum model dramatically underestimated the level of mechanical anisotropy in aligned collagen gels despite incorporation of measured fiber orientations; when estimated remodeling-induced changes in collagen fiber length were included, the continuum model slightly overestimated mechanical anisotropy. The network model provided the closest match to experimental data from aligned collagen gels, but still did not fully explain the observed mechanics. Two different modeling approaches showed that the level of collagen fiber alignment in our uniaxially constrained gels cannot explain the high degree of mechanical anisotropy observed in these gels. Our modeling results suggest that remodeling-induced redistribution of collagen fiber lengths, nonaffine fiber kinematics, or some combination of these effects must also be considered in order to explain the dramatic mechanical anisotropy observed in this collagen gel model system.     

A finite element based method to determine the properties of planar soft tissue.

A finite element based method to determine the incremental elastic material properties of planar membranes was developed and evaluated. The method is applicable to tissues that exhibit inhomogeneity, geometric and material nonlinearity, and anisotropy. Markers are placed on the tissue to form a four-node quadrilateral element. The specimen is loaded to an initial reference state, then three incremental loading sets are applied and the nodal displacements recorded. One of these loadings must include shear. These data are used to solve an over-determined system of equations for the tangent stiffness matrix. The method was first verified using analytical data. Next, data obtained from a latex rubber sheet were used to evaluate experimental procedures. Finally, experiments conducted on preconditioned rat skin revealed nonlinear orthotropic behavior. The vector norm comparing the applied and calculated nodal force vectors was used to evaluate the accuracy of the solutions.     

Mechanics of linear microcracking in trabecular bone

Microcracking in trabecular bone is responsible both for the mechanical degradation and remodeling of the trabecular bone tissue. Recent results on trabecular bone mechanics have demonstrated that bone tissue microarchitecture, tissue elastic heterogeneity and tissue-level mechanical anisotropy all should be considered to obtain detailed information on the mechanical stress state. The present study investigated the influence of tissue microarchitecture, tissue heterogeneity in elasticity and material separation properties and tissue-level anisotropy on the microcrack formation process. Microscale bone models were executed with the extended finite element method. It was demonstrated that anisotropy and heterogeneity of the bone tissue contribute significantly to bone tissue toughness and the resistance of trabecular bone to microcrack formation. The compressive strain to microcrack initiation was computed to increase by a factor of four from an assumed homogeneous isotropic tissue to an assumed anisotropic heterogenous tissue.     

Effects of boundary conditions on the estimation of the planar biaxial mechanical properties of soft tissues

Evaluation and simulation of the multiaxial mechanical behavior of native and engineered soft tissues is becoming more prevalent. In spite of this growing use, testing methods have not been standardized and methodologies vary widely. The strong influence of boundary conditions were recently underscored by Waldman et al. [2002, J. Materials Science: Materials in Medicine 13, pp. 933-938] wherein substantially different experimental results were obtained using different sample gripping methods on the same specimens. As it is not possible to experimentally evaluate the effects of different biaxial test boundary conditions on specimen internal stress distributions, we conducted numerical simulations to explore these effects. A nonlinear Fung-elastic constitutive model (Sun et al., 2003, JBME 125, pp. 372-380, which fully incorporated the effects of in-plane shear, was used to simulate soft tissue mechanical behavior. Effects of boundary conditions, including varying the number of suture attachments, different gripping methods, specimen shapes, and material axes orientations were examined. Results demonstrated strong boundary effects with the clamped methods, while suture attachment methods demonstrated minimal boundary effects. Suture-based methods appeared to be best suited for biaxial mechanical tests of biological materials. Moreover, the simulations demonstrated that Saint-Venant's effects depended significantly on the material axes orientation. While not exhaustive, these comprehensive simulations provide experimentalists with additional insight into the stress-strain fields associated with different biaxial testing boundary conditions, and may be used as a rational basis for the design of biaxial testing experiments.     

Anisotropy, inhomogeneity, and tension-compression nonlinearity of human glenohumeral cartilage in finite deformation

The tensile and compressive properties of human glenohumeral cartilage were determined by testing 120 rectangular strips in uniaxial tension and 70 cylindrical plugs in confined compression, obtained from five human glenohumeral joints. Specimens were harvested from five regions across the articular surface of the humeral head and two regions on the glenoid. Tensile strips were obtained along two orientations, parallel and perpendicular to the split-line directions. Two serial slices through the thickness, corresponding to the superficial and middle zones of the cartilage layers, were prepared from each tensile strip and each compressive plug. The equilibrium tensile modulus and compressive aggregate modulus of cartilage were determined from the uniaxial tensile and confined compression tests, respectively. Significant differences in the tensile moduli were found with depth and orientation relative to the local split-line direction. Articular cartilage of the humeral head was significantly stiffer in tension than that of the glenoid. There were significant differences in the aggregate compressive moduli of articular cartilage between superficial and middle zones in the humeral head. Furthermore, tensile and compressive stress-strain responses exhibited nonlinearity under finite strain, while the tensile modulus differed by up to two orders of magnitude from the compressive aggregate modulus at 0% strain, demonstrating a high degree of tension-compression nonlinearity. The complexity of the mechanical properties of human glenohumeral cartilage was exposed in this study, showing anisotropy, inhomogeneity, and tension-compression nonlinearity within the same joint. The observed differences in the tensile properties of human glenohumeral cartilage suggest that the glenoid may be more susceptible to cartilage degeneration than the humeral head.     

Pipette aspiration technique for the measurement of nonlinear and anisotropic mechanical properties of blood vessel walls under biaxial stretch

A pipette aspiration technique was proposed for the measurement of nonlinear mechanical properties of arteries under biaxial stretching. A cross-shaped specimen of porcine thoracic aorta whose principal axes corresponded with the axial and circumferential directions of the aortic walls was excised. The intraluminal surface of the specimen was aspirated with a circular cross-sectioned glass pipette while the specimen was stretching in the axial and circumferential directions in 10% increments. The elastic modulus agreed with the incremental elastic modulus obtained through a conventional pressure-diameter test of the same specimen to within an error of 30% at a circumferential stretch ratio below 1.3 and an axial stretch ratio of 1.0, 1.1 or 1.2, which represent lower range of physiological stretch ratios for the porcine aorta. A rectangular cross-sectioned pipette was utilized to measure anisotropic properties of the specimen under biaxial stretching. When aspirated with such a pipette, the specimens' elastic properties along the length of the rectangular pipette cross section can be neglected. The elastic modulus was found to increase rapidly when the specimen was stretched in the direction of the pipette's width. Thus, pipette aspiration should have many advantages such as well measurement of the local nonlinear and anisotropic mechanical properties of blood vessel walls.