Fiber kinematics of small intestinal submucosa under biaxial and uniaxial stretch

Improving our understanding of the design requirements of biologically derived collagenous scaffolds is necessary for their effective use in tissue reconstruction. In the present study, the collagen fiber kinematics of small intestinal submucosa (SIS) was quantified using small angle light scattering (SALS) while the specimen was subjected to prescribed uniaxial or biaxial strain paths. A modified biaxial stretching device based on Billiar and Sacks (J. Biomech., 30, pp. 753-7, 1997) was used, with a real-time analysis of the fiber kinematics made possible due to the natural translucency of SIS. Results indicated that the angular distribution of collagen fibers in specimens subjected to 10% equibiaxial strain was not significantly different from the initial unloaded condition, regardless of the loading path (p=0.31). Both 10% strip biaxial stretch and uniaxial stretches of greater than 5% in the preferred fiber direction led to an increase in the collagen fiber alignment along the same direction, while 10% strip biaxial stretch in the cross preferred fiber direction led to a broadening of the distribution. While an affine deformation model accurately predicted the experimental findings for a biaxial strain state, uniaxial stretch paths were not accurately predicted. Nonaffine structural models will be necessary to fully predict the fiber kinematics under large uniaxial strains in SIS.     

Biaxial mechanical response of bioprosthetic heart valve biomaterials to high in-plane shear

Utilization of novel biologically-derived biomaterials in bioprosthetic heart valves (BHV) requires robust constitutive models to predict the mechanical behavior under generalized loading states. Thus, it is necessary to perform rigorous experimentation involving all functional deformations to obtain both the form and material constants of a strain-energy density function. In this study, we generated a comprehensive experimental biaxial mechanical dataset that included high in-plane shear stresses using glutaraldehyde treated bovine pericardium (GLBP) as the representative BHV biomaterial. Compared to our previous study (Sacks, JBME, v.121, pp. 551-555, 1999), GLBP demonstrated a substantially different response under high shear strains. This finding was underscored by the inability of the standard Fung model, applied successfully in our previous GLBP study, to fit the high-shear data. To develop an appropriate constitutive model, we utilized an interpolation technique for the pseudo-elastic response to guide modification of the final model form. An eight parameter modified Fung model utilizing additional quartic terms was developed, which fitted the complete dataset well. Model parameters were also constrained to satisfy physical plausibility of the strain energy function. The results of this study underscore the limited predictive ability of current soft tissue models, and the need to collect experimental data for soft tissue simulations over the complete functional range.     

Contributions of Glycosaminoglycans to Collagen Fiber Recruitment in Constitutive Modeling of Arterial Mechanics

The contribution of glycosaminoglycans (GAGs) to the biological and mechanical functions of biological tissue has emerged as an important area of research. GAGs provide structural basis for the organization and assembly of extracellular matrix (ECM). The mechanics of tissue with low GAG content can be indirectly affected by the interaction of GAGs with collagen fibers, which have long been known to be one of the primary contributors to soft tissue mechanics. Our earlier study showed that enzymatic GAG depletion results in straighter collagen fibers that are recruited at lower levels of stretch, and a corresponding shift in earlier arterial stiffening (Mattson et al., 2016). In this study, the effect of GAGs on collagen fiber recruitment was studied through a structure-based constitutive model. The model incorporates structural information, such as fiber orientation distribution, content, and recruitment of medial elastin, medial collagen, and adventitial collagen fibers. The model was first used to study planar biaxial tensile stress-stretch behavior of porcine descending thoracic aorta. Changes in elastin and collagen fiber orientation distribution, and collagen fiber recruitment were then incorporated into the model in order to predict the stress-stretch behavior of GAG depleted tissue. Our study shows that incorporating early collagen fiber recruitment into the model predicts the stress-stretch response of GAG depleted tissue reasonably well (rms = 0.141); considering further changes of fiber orientation distribution does not improve the predicting capability (rms = 0.149). Our study suggests an important role of GAGs in arterial mechanics that should be considered in developing constitutive models.     

On constitutive relations and finite deformations of passive cardiac tissue: I. A pseudostrain-energy function

A three-dimensional constitutive relation for passive cardiac tissue is formulated in terms of a structurally motivated pseudostrain-energy function, W, while the mathematical simplicity of phenomenological approaches is preserved. A specific functional form of W is proposed on the basis of limited structural information and multiaxial experimental data. The material parameters are determined in a least-squared sense from both uniaxial and biaxial data. Our results suggest that (1) multiaxially-loaded cardiac tissue is nearly transversely-isotropic with respect to local muscle fiber directions, at least for a limited range of strain histories, (2) material parameters determined from uniaxial papillary muscle data result in gross underestimates of the stresses in multiaxially-loaded specimens, and (3) material parameters determined from equibiaxial tests predict the behavior of the tissue under various nonequibiaxial stretching protocols reasonably well.     

The effects of aneurysm on the biaxial mechanical behavior of human abdominal aorta

The biomechanical response of normal and pathologic human abdominal aortic tissue to uniaxial loading conditions is insufficient for the characterization of its three-dimensional (3D) mechanical behavior. Planar biaxial mechanical evaluation allows for 3D constitutive modeling of nearly incompressible tissues, as well as the investigation of the nature of mechanical anisotropy. In the current study, 26 abdominal aortic aneurysm (AAA) tissue samples and 8 age-matched (> 60 years of age) nonaneurysmal abdominal aortic (AA) tissue samples were obtained and tested using a tension-controlled biaxial testing protocol. Graphical response functions (Sun et al., 2003. J. Biomech. Eng. 125, 372-380) were used as a guide to describe the pseudo-elastic response of AA and AAA. Based on the observed pseudo-elastic response, a four-parameter exponential strain energy function developed by Vito (1990. J. Biomech. Eng. 112, 153-159) was used from which both an individual specimen and group material parameter sets were determined for both AA and AAA. Peak Green strain values in the circumferential (Ethetatheta,max) and longitudinal (ELL,max) directions under an equibiaxial tension of 120 N/m were also compared. The strain energy function fit all of the individual specimens well with an average R2 of 0.95 +/- 0.02 and 0.90 +/- 0.02 (mean +/- SEM) for the AA and AAA groups, respectively. The average Ethetatheta,max at 200 N/m equibiaxial tension was found to be significantly smaller for AAAs as compared to AAs (0.07 +/- 0.01 versus 0.13 +/- 0.03, respectively; p < 0.01). There was also a pronounced increase in the circumferential stiffness for AAA tissue as compared to AA tissue, indicating a larger degree of anisotropy for this tissue as compared to age-matched AA tissue. We also observed that the four-parameter Fung-elastic model was not able to fit the AAA tissue mechanical response using physically realistic material parameter values. It was concluded that aneurysmal degeneration of the abdominal aorta is associated with an increase in mechanical anisotropy, with preferential stiffening in the circumferential direction.     

Biaxial tensile testing and constitutive modeling of human supraspinatus tendon

The heterogeneous composition and mechanical properties of the supraspinatus tendon offer an opportunity for studying the structure-function relationships of fibrous musculoskeletal connective tissues. Previous uniaxial testing has demonstrated a correlation between the collagen fiber angle distribution and tendon mechanics in response to tensile loading both parallel and transverse to the tendon longitudinal axis. However, the planar mechanics of the supraspinatus tendon may be more appropriately characterized through biaxial tensile testing, which avoids the limitation of nonphysiologic traction-free boundary conditions present during uniaxial testing. Combined with a structural constitutive model, biaxial testing can help identify the specific structural mechanisms underlying the tendon's two-dimensional mechanical behavior. Therefore, the objective of this study was to evaluate the contribution of collagen fiber organization to the planar tensile mechanics of the human supraspinatus tendon by fitting biaxial tensile data with a structural constitutive model that incorporates a sample-specific angular distribution of nonlinear fibers. Regional samples were tested under several biaxial boundary conditions while simultaneously measuring the collagen fiber orientations via polarized light imaging. The histograms of fiber angles were fit with a von Mises probability distribution and input into a hyperelastic constitutive model incorporating the contributions of the uncrimped fibers. Samples with a wide fiber angle distribution produced greater transverse stresses than more highly aligned samples. The structural model fit the longitudinal stresses well (median R(2) ≥ 0.96) and was validated by successfully predicting the stress response to a mechanical protocol not used for parameter estimation. The transverse stresses were fit less well with greater errors observed for less aligned samples. Sensitivity analyses and relatively affine fiber kinematics suggest that these errors are not due to inaccuracies in measuring the collagen fiber organization. More likely, additional strain energy terms representing fiber-fiber interactions are necessary to provide a closer approximation of the transverse stresses. Nevertheless, this approach demonstrated that the longitudinal tensile mechanics of the supraspinatus tendon are primarily dependent on the moduli, crimp, and angular distribution of its collagen fibers. These results add to the existing knowledge of structure-function relationships in fibrous musculoskeletal tissue, which is valuable for understanding the etiology of degenerative disease, developing effective tissue engineering design strategies, and predicting outcomes of tissue repair.     

Arterial Extracellular Matrix: A Mechanobiological Study of the Contributions and Interactions of Elastin and Collagen

The complex network structure of elastin and collagen extracellular matrix (ECM) forms the primary load bearing components in the arterial wall. The structural and mechanobiological interactions between elastin and collagen are important for properly functioning arteries. Here, we examined the elastin and collagen organization, realignment, and recruitment by coupling mechanical loading and multiphoton imaging. Two-photon excitation fluorescence and second harmonic generation methods were performed with a multiphoton video-rate microscope to capture real time changes to the elastin and collagen structure during biaxial deformation. Enzymatic removal of elastin was performed to assess the structural changes of the remaining collagen structure. Quantitative analysis of the structural changes to elastin and collagen was made using a combination of two-dimensional fast Fourier transform and fractal analysis, which allows for a more complete understanding of structural changes. Our study provides new quantitative evidence, to our knowledge on the sequential engagement of different arterial ECM components in response to mechanical loading. The adventitial collagen exists as large wavy bundles of fibers that exhibit fiber engagement after 20% strain. The medial collagen is engaged throughout the stretching process, and prominent elastic fiber engagement is observed up to 20% strain after which the engagement plateaus. The fiber orientation distribution functions show remarkably different changes in the ECM structure in response to mechanical loading. The medial collagen shows an evident preferred circumferential distribution, however the fiber families of adventitial collagen are obscured by their waviness at no or low mechanical strains. Collagen fibers in both layers exhibit significant realignment in response to unequal biaxial loading. The elastic fibers are much more uniformly distributed and remained relatively unchanged due to loading. Removal of elastin produces similar structural changes in collagen as mechanical loading. Our study suggests that the elastic fibers are under tension and impart an intrinsic compressive stress on the collagen.     

Arterial Extracellular Matrix: A Mechanobiological Study of the Contributions and Interactions of Elastin and Collagen

The complex network structure of elastin and collagen extracellular matrix (ECM) forms the primary load bearing components in the arterial wall. The structural and mechanobiological interactions between elastin and collagen are important for properly functioning arteries. Here, we examined the elastin and collagen organization, realignment, and recruitment by coupling mechanical loading and multiphoton imaging. Two-photon excitation fluorescence and second harmonic generation methods were performed with a multiphoton video-rate microscope to capture real time changes to the elastin and collagen structure during biaxial deformation. Enzymatic removal of elastin was performed to assess the structural changes of the remaining collagen structure. Quantitative analysis of the structural changes to elastin and collagen was made using a combination of two-dimensional fast Fourier transform and fractal analysis, which allows for a more complete understanding of structural changes. Our study provides new quantitative evidence, to our knowledge on the sequential engagement of different arterial ECM components in response to mechanical loading. The adventitial collagen exists as large wavy bundles of fibers that exhibit fiber engagement after 20% strain. The medial collagen is engaged throughout the stretching process, and prominent elastic fiber engagement is observed up to 20% strain after which the engagement plateaus. The fiber orientation distribution functions show remarkably different changes in the ECM structure in response to mechanical loading. The medial collagen shows an evident preferred circumferential distribution, however the fiber families of adventitial collagen are obscured by their waviness at no or low mechanical strains. Collagen fibers in both layers exhibit significant realignment in response to unequal biaxial loading. The elastic fibers are much more uniformly distributed and remained relatively unchanged due to loading. Removal of elastin produces similar structural changes in collagen as mechanical loading. Our study suggests that the elastic fibers are under tension and impart an intrinsic compressive stress on the collagen.     

Biaxial mechanical properties of the natural and glutaraldehyde treated aortic valve cusp--Part I: Experimental results

To date, there are no constitutive models for either the natural or bioprosthetic aortic valve (AV), in part due to experimental complications related to the AV's small size and heterogeneous fibrous structure. In this study, we developed specialized biaxial testing techniques for the AV cusp, including a method to determine the local structure-strain relationship to assess the effects of boundary tethering forces. Natural and glutaraldehyde (GL) treated cusps were subjected to an extensive biaxial testing protocol in which the ratios of the axial tensions were held at constant values. Results indicated that the local fiber architecture clearly dominated cuspal deformation, and that the tethering effects at the specimen boundaries were negligible. Due to unique aspects of cuspal fiber architecture, the most uniform region of deformation was found at the lower portion as opposed to the center of the cuspal specimen. In general, the circumferential strains were much smaller than the radial strains, indicating a profound degree of mechanical anisotropy, and that natural cusps were significantly more extensible than the GL treated cusps. Strong mechanical coupling between biaxial stretch axes produced negative circumferential strains under equibiaxial tension. Further, the large radial strains observed could not be explained by uncrimping of the collagen fibers, but may be due to large rotations of the highly aligned, circumferential-oriented collagen fibers in the fibrosa. In conclusion, this study provides new insights into the AV cusp's structure-function relationship in addition to requisite data for constitutive modeling.     

The relation between collagen fibril kinematics and mechanical properties in the mitral valve anterior leaflet

We have recently demonstrated that the mitral valve anterior leaflet (MVAL) exhibited minimal hysteresis, no strain rate sensitivity, stress relaxation but not creep (Grashow et al., 2006, Ann Biomed Eng., 34(2), pp. 315-325; Grashow et al., 2006, Ann Biomed. Eng., 34(10), pp. 1509-1518). However, the underlying structural basis for this unique quasi-elastic mechanical behavior is presently unknown. As collagen is the major structural component of the MVAL, we investigated the relation between collagen fibril kinematics (rotation and stretch) and tissue-level mechanical properties in the MVAL under biaxial loading using small angle X-ray scattering. A novel device was developed and utilized to perform simultaneous measurements of tissue level forces and strain under a planar biaxial loading state. Collagen fibril D-period strain (epsilonD) and the fibrillar angular distribution were measured under equibiaxial tension, creep, and stress relaxation to a peak tension of 90 N/m. Results indicated that, under equibiaxial tension, collagen fibril straining did not initiate until the end of the nonlinear region of the tissue-level stress-strain curve. At higher tissue tension levels, epsilonD increased linearly with increasing tension. Changes in the angular distribution of the collagen fibrils mainly occurred in the tissue toe region. Using epsilonD, the tangent modulus of collagen fibrils was estimated to be 95.5+/-25.5 MPa, which was approximately 27 times higher than the tissue tensile tangent modulus of 3.58+/-1.83 MPa. In creep tests performed at 90 N/m equibiaxial tension for 60 min, both tissue strain and epsilonD remained constant with no observable changes over the test length. In contrast, in stress relaxation tests performed for 90 min epsilonD was found to rapidly decrease in the first 10 min followed by a slower decay rate for the remainder of the test. Using a single exponential model, the time constant for the reduction in collagen fibril strain was 8.3 min, which was smaller than the tissue-level stress relaxation time constants of 22.0 and 16.9 min in the circumferential and radial directions, respectively. Moreover, there was no change in the fibril angular distribution under both creep and stress relaxation over the test period. Our results suggest that (1) the MVAL collagen fibrils do not exhibit intrinsic viscoelastic behavior, (2) tissue relaxation results from the removal of stress from the fibrils, possibly by a slipping mechanism modulated by noncollagenous components (e.g. proteoglycans), and (3) the lack of creep but the occurrence of stress relaxation suggests a "load-locking" behavior under maintained loading conditions. These unique mechanical characteristics are likely necessary for normal valvular function.     

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.     

A novel constitutive model for passive right ventricular myocardium: evidence for myofiber–collagen fiber mechanical coupling

The function of right ventricle (RV) is recognized to play a key role in the development of many cardiopulmonary disorders, such as pulmonary arterial hypertension (PAH). Given the strong link between tissue structure and mechanical behavior, there remains a need for a myocardial constitutive model that accurately accounts for right ventricular myocardium architecture. Moreover, most available myocardial constitutive models approach myocardium at the length scale of mean fiber orientation and do not explicitly account for different fibrous constituents and possible interactions among them. In the present work, we developed a fiber-level constitutive model for the passive mechanical behavior of the right ventricular free wall (RVFW). The model explicitly separates the mechanical contributions of myofiber and collagen fiber ensembles, and accounts for the mechanical interactions between them. To obtain model parameters for the healthy passive RVFW, the model was informed by transmural orientation distribution measurements of myo- and collagen fibers and was fit to the mechanical testing data, where both sets of data were obtained from recent experimental studies on non-contractile, but viable, murine RVFW specimens. Results supported the hypothesis that in the low-strain regime, the behavior of the RVFW is governed by myofiber response alone, which does not demonstrate any coupling between different myofiber ensembles. At higher strains, the collagen fibers and their interactions with myofibers begin to gradually contribute and dominate the behavior as recruitment proceeds. Due to the use of viable myocardial tissue, the contribution of myofibers was significant at all strains with the predicted tensile modulus of \(\sim \)32 kPa. This was in contrast to earlier reports (Horowitz et al. 1988) where the contribution of myofibers was found to be insignificant. Also, we found that the interaction between myo- and collagen fibers was greatest under equibiaxial strain, with its contribution to the total stress not exceeding 20 %. The present model can be applied to organ-level computational models of right ventricular dysfunction for efficient diagnosis and evaluation of pulmonary hypertension disorder.     

Molecular dynamics simulations showing 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) membrane mechanoporation damage under different strain paths

Continuum finite element material models used for traumatic brain injury lack local injury parameters necessitating nanoscale mechanical injury mechanisms be incorporated. One such mechanism is membrane mechanoporation, which can occur during physical insults and can be devastating to cells, depending on the level of disruption. The current study investigates the strain state dependence of phospholipid bilayer mechanoporation and failure. Using molecular dynamics, a simplified membrane, consisting of 72 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) phospholipids, was subjected to equibiaxial, 2:1 non-equibiaxial, 4:1 non-equibiaxial, strip biaxial, and uniaxial tensile deformations at a von Mises strain rate of 5.45 × 10⁸ s^?1, resulting in velocities in the range of 1 to 4.6 m·s^?1. A water bridge forming through both phospholipid bilayer leaflets was used to determine structural failure. The stress magnitude, failure strain, headgroup clustering, and damage responses were found to be strain state-dependent. The strain state order of detrimentality in descending order was equibiaxial, 2:1 non-equibiaxial, 4:1 non-equibiaxial, strip biaxial, and uniaxial. The phospholipid bilayer failed at von Mises strains of .46, .47, .53, .77, and 1.67 during these respective strain path simulations. Additionally, a Membrane Failure Limit Diagram (MFLD) was created using the pore nucleation, growth, and failure strains to demonstrate safe and unsafe membrane deformation regions. This MFLD allowed representative equations to be derived to predict membrane failure from in-plane strains. These results provide the basis to implement a more accurate mechano-physiological internal state variable continuum model that captures lower length scale damage and will aid in developing higher fidelity injury models.     

Structural three-dimensional constitutive law for the passive myocardium

A three-dimensional constitutive law is proposed for the myocardium. Its formulation is based on a structural approach in which the total strain energy of the tissue is the sum of the strain energies of its constituents: the muscle fibers, the collagen fibers and the fluid matrix which embeds them. The ensuing material law expresses the specific structural and mechanical properties of the tissue, namely, the spatial orientation of the comprising fibers, their waviness in the unstressed state and their stress-strain behavior when stretched. Having assumed specific functional forms for the distribution of the fibers spatial orientation and waviness, the results of biaxial mechanical tests serve for the estimation of the material constants appearing in the constitutive equations. A very good fit is obtained between the measured and the calculated stresses, indicating the suitability of the proposed model for describing the mechanical behavior of the passive myocardium. Moreover, the results provide general conclusions concerning the structural basis for the tissue overall mechanical properties, the main of which is that the collagen matrix, though comprising a relatively small fraction of the whole tissue volume, is the dominant component accounting for its stiffness.     

Biaxial mechanical testing of posterior sclera using high-resolution ultrasound speckle tracking for strain measurements

This study aimed to characterize the mechanical responses of the sclera, the white outer coat of the eye, under equal-biaxial loading with unrestricted shear. An ultrasound speckle tracking technique was used to measure tissue deformation through sample thickness, expanding the capabilities of surface strain techniques. Eight porcine scleral samples were tested within 72 h postmortem  . High resolution ultrasound scans of scleral cross-sections along the two loading axes were acquired at 25 consecutive biaxial load levels. An additional repeat of the biaxial loading cycle was performed to measure a third normal strain emulating a strain gage rosette for calculating the in-plane shear. The repeatability of the strain measurements during identical biaxial ramps was evaluated. A correlation-based ultrasound speckle tracking algorithm was used to compute the displacement field and determine the distributive strains in the sample cross-sections. A Fung type constitutive model including a shear term was used to determine the material constants of each individual specimen by fitting the model parameters to the experimental stress–strain data. A non-linear stress–strain response was observed in all samples. The meridian direction had significantly larger strains than that of the circumferential direction during equal-biaxial loadings (P's<0.05$P\text{'}s<0.05$). The stiffness along the two directions was also significantly different (P=0.02) but highly correlated (R²=0.8). These results showed that the mechanical properties of the porcine sclera were nonlinear and anisotropic under biaxial loading. This work has also demonstrated the feasibility of using ultrasound speckle tracking for strain measurements during mechanical testing.     

Development of fibroblast-seeded collagen gels under planar biaxial mechanical constraints: a biomechanical study

Prior studies indicated that mechanical loading influences cell turnover and matrix remodeling in tissues, suggesting that mechanical stimuli can play an active role in engineering artificial tissues. While most tissue culture studies focus on influence of uniaxial loading or constraints, effects of multi-axial loading or constraints on tissue development are far from clear. In this study, we examined the biaxial mechanical properties of fibroblast-seeded collagen gels cultured under four different mechanical constraints for 6 days: free-floating, equibiaxial stretching (with three different stretch ratios), strip-biaxial stretching, and uniaxial stretching. Passive mechanical behavior of the cell-seeded gels was also examined after decellularization. A continuum-based two-dimensional Fung model was used to quantify the mechanical behavior of the gel. Based on the model, the value of stored strain energy and the ratio of stiffness in the stretching directions were calculated at prescribed strains for each gel, and statistical comparisons were made among the gels cultured under the various mechanical constraints. Results showed that gels cultured under the free-floating and equibiaxial stretching conditions exhibited a nearly isotropic mechanical behavior, while gels cultured under the strip-biaxial and uniaxial stretching conditions developed a significant degree of mechanical anisotropy. In particular, gels cultured under the equibiaxial stretching condition with a greater stretch ratio appeared to be stiffer than those with a smaller stretch ratio. Also, a decellularized gel was stiffer than its non-decellularized counterpart. Finally, the retained mechanical anisotropy in gels cultured under the strip-biaxial stretching and uniaxial stretching conditions after cell removal reflected an irreversible matrix remodeling.     

On the biaxial mechanical properties of the layers of the aortic valve leaflet

All existing constitutive models for heart valve leaflet tissues either assume a uniform transmural stress distribution or utilize a membrane tension formulation. Both approaches ignore layer specific mechanical contributions and the implicit nonuniformity of the transmural stress distribution. To begin to address these limitations, we conducted novel studies to quantify the biaxial mechanical behavior of the two structurally distinct, load bearing aortic valve (AV) leaflet layers: the fibrosa and ventricularis. Strip biaxial tests, with extremely sensitive force sensing capabilities, were further utilized to determine the mechanical behavior of the separated ventricularis layer at very low stress levels. Results indicated that both layers exhibited very different nonlinear, highly anisotropic mechanical behaviors. While the leaflet tissue mechanical response was dominated by the fibrosa layer, the ventricularis contributed double the amount of the fibrosa to the total radial tension and experienced four times the stress level. The strip biaxial test results further indicated that the ventricularis exhibited substantial anisotropic mechanical properties at very low stress levels. This result suggested that for all strain levels, the ventricularis layer is dominated by circumferentially oriented collagen fibers, and the initial loading phase of this layer cannot be modeled as an isotropic material. Histological-based thickness measurements indicated that the fibrosa and ventricularis constitute 41% and 29% of the total layer thickness, respectively. Moreover, the extensive network of interlayer connections and identical strains under biaxial loading in the intact state suggests that these layers are tightly bonded. In addition to advancing our knowledge of the subtle but important mechanical properties of the AV leaflet, this study provided a comprehensive database required for the development of a true 3D stress constitutive model for the native AV leaflet.     

A three-dimensional constitutive model for the stress relaxation of articular ligaments

A new nonlinear constitutive model for the three-dimensional stress relaxation of articular ligaments is proposed. The model accounts for finite strains, anisotropy, and strain-dependent stress relaxation behavior exhibited by these ligaments. The model parameters are identified using published uniaxial stress–stretch and stress relaxation data on human medial collateral ligaments (MCLs) subjected to tensile tests in the fiber and transverse to the fiber directions (Quapp and Weiss in J Biomech Eng Trans ASME 120:757–763, 1998; Bonifasi-Lista et al. in J Orthop Res 23(1):67–76, 2005). The constitutive equation is then used to predict the nonlinear elastic and stress relaxation response of ligaments subjected to shear deformations in the fiber direction and transverse to the fiber direction, and an equibiaxial extension. A direct comparison with stress relaxation data collected by subjecting human MCLs to shear deformation in the fiber direction is presented in order to demonstrate the predictive capabilities of the model.     

A constitutive model for vascular tissue that integrates fibril, fiber and continuum levels with application to the isotropic and passive properties of the infrarenal aorta

A fundamental understanding of the mechanical properties of the extracellular matrix (ECM) is critically important to quantify the amount of macroscopic stress and/or strain transmitted to the cellular level of vascular tissue. Structural constitutive models integrate histological and mechanical information, and hence, allocate stress and strain to the different microstructural components of the vascular wall. The present work proposes a novel multi-scale structural constitutive model of passive vascular tissue, where collagen fibers are assembled by proteoglycan (PG) cross-linked collagen fibrils and reinforce an otherwise isotropic matrix material. Multiplicative kinematics account for the straightening and stretching of collagen fibrils, and an orientation density function captures the spatial organization of collagen fibers in the tissue. Mechanical and structural assumptions at the collagen fibril level define a piece-wise analytical stress-stretch response of collagen fibers, which in turn is integrated over the unit sphere to constitute the tissue's macroscopic mechanical properties. The proposed model displays the salient macroscopic features of vascular tissue, and employs the material and structural parameters of clear physical meaning. Likewise, the constitutive concept renders a highly efficient multi-scale structural approach that allows for the numerical analysis at the organ level. Model parameters were estimated from isotropic mean-population data of the normal and aneurysmatic aortic wall and used to predict in-vivo stress states of patient-specific vascular geometries, thought to demonstrate the robustness of the particular Finite Element (FE) implementation. The collagen fibril level of the multi-scale constitutive formulation provided an interface to integrate vascular wall biology and to account for collagen turnover.     

Simulation of planar soft tissues using a structural constitutive model: Finite element implementation and validation

Computational implementation of physical and physiologically realistic constitutive models is critical for numerical simulation of soft biological tissues in a variety of biomedical applications. It is well established that the highly nonlinear and anisotropic mechanical behaviors of soft tissues are an emergent behavior of the underlying tissue microstructure. In the present study, we have implemented a structural constitutive model into a finite element framework specialized for membrane tissues. We noted that starting with a single element subjected to uniaxial tension, the non-fibrous tissue matrix must be present to prevent unrealistic tissue deformations. Flexural simulations were used to set the non-fibrous matrix modulus because fibers have little effects on tissue deformation under three-point bending. Multiple deformation modes were simulated, including strip biaxial, planar biaxial with two attachment methods, and membrane inflation. Detailed comparisons with experimental data were undertaken to insure faithful simulations of both the macro-level stress–strain insights into adaptations of the fiber architecture under stress, such as fiber reorientation and fiber recruitment. Results indicated a high degree of fidelity and demonstrated interesting microstructural adaptions to stress and the important role of the underlying tissue matrix. Moreover, we apparently resolve a discrepancy in our 1997 study (Billiar and Sacks, 1997. J. Biomech. 30 (7), 753–756) where we observed that under strip biaxial stretch the simulated fiber splay responses were not in good agreement with the experimental results, suggesting non-affine deformations may have occurred. However, by correctly accounting for the isotropic phase of the measured fiber splay, good agreement was obtained. While not the final word, these simulations suggest that affine fiber kinematics for planar collagenous tissues is a reasonable assumption at the macro level. Simulation tools such as these are imperative in the design and simulation of native and engineered tissues.