Search

Determining tissue biomechanical material properties from mechanical test data is frequently required in a variety of applications. However, the validity of the resulting constitutive model parameters is the subject of debate in the field. Parameter optimization in tissue mechanics often comes down to the “identifiability” or “uniqueness” of constitutive model parameters; however, despite advances in formulating complex constitutive relations and many classic and creative curve-fitting approaches, there is currently no accessible framework to study the identifiability of tissue material parameters. Our objective was to assess the identifiability of material parameters for established constitutive models of fiber-reinforced soft tissues, biomaterials, and tissue-engineered constructs and establish a generalizable procedure for other applications. To do so, we generated synthetic experimental data by simulating uniaxial tension and compression tests, commonly used in biomechanics. We then fit this data using a multi-start optimization technique based on the nonlinear least-squares method with multiple initial parameter guesses. We considered tendon and sclera as example tissues, using constitutive models that describe these fiber-reinforced tissues. We demonstrated that not all the model parameters of these constitutive models were identifiable from uniaxial mechanical tests, despite achieving virtually identical fits to the stress-stretch response. We further show that when the lateral strain was considered as an additional fitting criterion, more parameters are identifiable, but some remain unidentified. This work provides a practical approach for addressing parameter identifiability in tissue mechanics.

PurposeTo evaluate the effects of negative periocular pressure (NPP), and concomitant intraocular pressure (IOP) lowering, on the biomechanics of the optic nerve head (ONH) and cornea.MethodsWe developed a validated finite element (FE) model of the eye to compute tissue biomechanical strains induced in response to NPP delivered using the Multi-Pressure Dial (MPD) system. The model was informed by clinical measurements of IOP lowering and was based on published tissue properties. We also conducted sensitivity analyses by changing pressure loads and tissue properties.ResultsApplication of -7.9 mmHg NPP decreased strain magnitudes in the ONH by c. 50% while increasing corneal strain magnitudes by c. 25%. Comparatively, a similar increase in corneal strain was predicted to occur due to an increase in IOP of 4 mmHg. Sensitivity studies indicated that NPP lowers strain in the ONH by reducing IOP and that these effects persisted over a range of tissue stiffnesses and spatial distributions of NPP.ConclusionsNPP is predicted to considerably decrease ONH strain magnitudes. It also increases corneal strain but to an extent expected to be clinically insignificant. Thus, using NPP to lower IOP and hence decrease ONH mechanical strain is likely biomechanically beneficial for glaucoma patients.Translational RelevanceThis study provides the first description of how NPP affects ONH biomechanics and explains the underlying mechanism of ONH strain reduction. It complements current empirical knowledge about the MPD system and guides future studies of NPP as a treatment for glaucoma.

The iris is a muscular organ whose deformations can cause primary angle-closure glaucoma (PACG), a leading cause of blindness. PACG risk assessment does not consider iridial biomechanical factors, despite their expected influence on iris deformations. Here we exploited an existing biometric data set consisting of near-infrared movies acquired during the pupillary light reflex (PLR) as a unique resource to study iris biomechanics. The PLR caused significant (>100%) and essentially spatially uniform radial strains in the iris in vivo, consistent with previous findings. Inverse finite element modeling showed that sphincter muscle tractions were c. 5-fold greater than iridial stroma stiffness (range 4- to 13-fold, depending on sphincter muscle size). This muscle traction is greater than has been previously estimated, which may be due to methodological differences and/or to different patient populations in our study (European descent) vs. previous studies (Asian); the latter possibility is of particular interest due to differential incidence rates of PACG in these populations. Our methodology is fast and inexpensive and may be a useful tool in understanding biomechanical factors contributing to PACG.

Optic nerve head (ONH) biomechanics is centrally involved in the pathogenesis of glaucoma, a blinding ocular condition often characterized by elevation and fluctuation of the intraocular pressure and resulting loads on the ONH. Further, tissue viscoelasticity is expected to strongly influence the mechanical response of the ONH to mechanical loading, yet the viscoelastic mechanical properties of the ONH remain unknown. To determine these properties, we conducted micromechanical testing on porcine ONH tissue samples, coupled with finite element modeling based on a mixture model consisting of a biphasic material with a viscoelastic solid matrix. Our results provide a detailed description of the viscoelastic properties of the porcine ONH at each of its four anatomical quadrants (i.e., nasal, superior, temporal, and inferior). We showed that the ONH's viscoelastic mechanical response can be explained by a dual mechanism of fluid flow and solid matrix viscoelasticity, as is common in other soft tissues. We obtained porcine ONH properties as follows: matrix Young's modulus E=1.895[1.056,2.391] kPa (median [min., max.]), Poisson's ratio ν=0.142[0.060,0.312], kinetic time-constant τ=214[89,921] sec, and hydraulic permeability k=3.854×10

Optic nerve head (ONH) biomechanics is centrally involved in the pathogenesis of glaucoma, a blinding ocular condition often characterized by elevation and fluctuation of the intraocular pressure and resulting loads on the ONH. Further, tissue viscoelasticity is expected to strongly influence the mechanical response of the ONH to mechanical loading, yet the viscoelastic mechanical properties of the ONH remain unknown. To determine these properties, we conducted micromechanical testing on porcine ONH tissue samples, coupled with finite element modeling based on a mixture model consisting of a biphasic material with a viscoelastic solid matrix. Our results provide a detailed description of the viscoelastic properties of the porcine ONH at each of its four anatomical quadrants (i.e., nasal, superior, temporal, and inferior). We showed that the ONH's viscoelastic mechanical response can be explained by a dual mechanism of fluid flow and solid matrix viscoelasticity, as is common in other soft tissues. We obtained porcine ONH properties as follows: matrix Young's modulus E = 1.895 [ 1.056 , 2.391 ] kPa (median [min., max.]), Poisson's ratio ν = 0.142 [ 0.060 , 0.312 ] , kinetic time-constant τ = 214 [ 89 , 921 ] sec, and hydraulic permeability k = 3.854 × 10 − 1 [ 3.457 × 10 − 2 , 9.994 × 10 − 1 ] mm4/(N.sec). These values can be used to design and fabricate physiologically appropriate ex vivo test environments (e.g., 3D cell culture) to further understand glaucoma pathophysiology. Statement of significance Optic nerve head (ONH) biomechanics is an important aspect of the pathogenesis of glaucoma, the leading cause of irreversible blindness. The ONH experiences time-varying loads, yet the viscoelastic behavior of this tissue has not been characterized. Here, we measure the time-dependent response of the ONH in porcine eyes and use mechanical modeling to provide time-dependent mechanical properties of the ONH. This information allows us to identify time-varying stimuli in vivo which have timescales matching the characteristic response times of the ONH, and can also be used to design and fabricate ex vivo 3D cultures to study glaucoma pathophysiology in a physiologically relevant environment, enabling the discovery of new generations of glaucoma medications focusing on neuroprotection.

Determining tissue biomechanical material properties from mechanical test data is frequently required in a variety of applications, e.g. tissue engineering. However, the validity of the resulting constitutive model parameters is the subject of debate in the field. Common methods to perform fitting, such as nonlinear least-squares, are known to be subject to several limitations, most notably the uniqueness of the fitting results. Parameter optimization in tissue mechanics often comes down to the “identifiability” or “uniqueness” of constitutive model parameters; however, despite advances in formulating complex constitutive relations and many classic and creative curve-fitting approaches, there is no accessible framework to study the identifiability of tissue material parameters. Our objective was to assess the identifiability of material parameters for established constitutive models of fiber-reinforced soft tissues, biomaterials, and tissue-engineered constructs. To do so, we generated synthetic experimental data by simulating uniaxial tension and compression tests, commonly used in biomechanics. We considered tendon and sclera as example tissues, using constitutive models that describe these fiber-reinforced tissues. We demonstrated that not all of the model parameters of these constitutive models were identifiable from uniaxial mechanical tests, despite achieving virtually identical fits to the stress-stretch response. We further show that when the lateral strain was considered as an additional fitting criterion, more parameters are identifiable, but some remain unidentified. This work provides a practical approach for addressing parameter identifiability in tissue mechanics.Statement of SignificanceData fitting is a powerful technique commonly used to extract tissue material parameters from experimental data, and which thus has applications in tissue biomechanics and engineering. However, the problem of “uniqueness” or “identifiability” of the fit parameters is a significant issue, limiting the fit results’ validity. Here we provide a novel method to evaluate data fitting and assess the uniqueness of results in the tissue mechanics constitutive models. Our results indicate that the uniaxial stress-stretch experimental data are not adequate to identify all the tissue material parameters. This study is of potential interest to a wide range of readers because of its application for the characterization of other engineering materials, while addressing the problem of uniqueness of the fitted results.

A buffer solution is often used to maintain tissue hydration during mechanical testing. The most commonly used buffer solution is a physiological concentration of phosphate buffered saline (PBS); however, PBS increases the tissue’s water content and decreases its tensile stiffness. In addition, solutes from the buffer can diffuse into the tissue and interact with its structure and mechanics. These bathing solution effects can confound the outcome and interpretation of mechanical tests. Potential bathing solution artifacts, including solute diffusion and the effect on mechanical properties, are not well understood. The objective of this study was to measure the effects of long-term exposure of rat tail tendon fascicles to several concentrations (0.9% to 25%) of NaCl, sucrose, polyethylene glycol (PEG), and SPEG (NaCl + PEG) solutions on water content, solute diffusion, and mechanical properties. We found that with an increase in solute concentration the apparent water content decreased for all solution types. Solutes diffused into the tissue for NaCl and sucrose, however, no solute diffusion was observed for PEG or SPEG. The mechanical properties changed for both of NaCl solutions, in particular after long-term (8 hr) incubation the modulus and equilibrium stress decreased compared to short-term (15 min) for 25% NaCl, and the cross sectional area increased for 0.9% NaCl. However, the mechanical properties were unchanged for both PEG and SPEG except for minor alterations in stress relaxation parameters. This study shows that NaCl and sucrose buffer solutions are not suitable for long-term mechanical tests. We therefore propose using PEG or SPEG as alternative buffer solutions that after long-term incubation can maintain tissue hydration without solute diffusion and produce a consistent mechanical response.

Tendon’s hierarchical structure allows for load transfer between its fibrillar elements at multiple length scales. Tendon microstructure is particularly important, because it includes the cells and their surrounding collagen fibrils, where mechanical interactions can have potentially important physiological and pathological contributions. However, the three-dimensional microstructure and the mechanisms of load transfer in that length scale are not known. It has been postulated that interfibrillar matrix shear or direct load transfer via the fusion/branching of small fibrils are responsible for load transfer, but the significance of these mechanisms is still unclear. Alternatively, the helical fibrils that occur at the microstructural scale in tendon may also mediate load transfer, however, these structures are not well studied due to the lack of a three-dimensional visualization of tendon microstructure. In this study, we used serial block-face scanning electron microscopy (SBF-SEM) to investigate the threedimensional microstructure of fibrils in rat tail tendon. We found that tendon fibrils have a complex architecture with many helically wrapped fibrils. We studied the mechanical implications of these helical structures using finite element modeling and found that frictional contact between helical fibrils can induce load transfer even in the absence of matrix bonding or fibril fusion/branching. This study is significant in that it provides a three-dimensional view of the tendon microstructure and suggests friction between helically wrapped fibrils as a mechanism for load transfer, which is an important aspect of tendon biomechanics.

Fibrous soft tissues are biopolymeric materials that are made of extracellular proteins, such as different types of collagen and proteoglycans, and have a high water content. These tissues have nonlinear, anisotropic, and inelastic mechanical behaviors that are often categorized into viscoelastic behavior, plastic deformation, and damage. While tissue's elastic and viscoelastic mechanical properties have been measured for decades, there is no comprehensive theoretical framework for modeling inelastic behaviors of these tissues that is based on their structure. To model the three major inelastic mechanical behaviors of tissue's fibrous matrix, we formulated a structurally inspired continuum mechanics framework based on the energy of molecular bonds that break and reform in response to external loading (reactive bonds). In this framework, we employed the theory of internal state variables (ISV) and kinetics of molecular bonds. The number fraction of bonds, their reference deformation gradient, and damage parameter were used as state variables that allowed for consistent modeling of all three of the inelastic behaviors of tissue by using the same sets of constitutive relations. Several numerical examples are provided that address practical problems in tissue mechanics, including the difference between plastic deformation and damage. This model can be used to identify relationships between tissue's mechanical response to external loading and its biopolymeric structure.

Tendon’s viscoelastic behaviors are important to the mechanical function and mechanobiology. When loaded in longitudinal tension, tendons often have a very large Poisson’s ratio (ν > 2) that exceeds the limit of incompressibility for isotropic material (ν = 0.5), indicating that tendon experiences volume loss, inducing poroelastic fluid exudation in the transverse direction. Therefore, transverse poroelasticity is an important contributor to tendon material behavior. Tendon hydraulic permeability which is required to evaluate the fluid flow contribution to viscoelasticity, is mostly unavailable, and where available, varies by several orders of magnitude. In this manuscript, we quantified the transverse poroelastic material parameters of rat tail tendon fascicles by conducting transverse osmotic loading experiments, in both tension and compression. We used a multi-start optimization method to evaluate the parameters based on biphasic finite element modeling. Our tendon samples had a transverse hydraulic permeability of 10(−4) to 10(−5) mm(4) (Ns)(−1) and showed a significant tension-compression nonlinearity in the transverse direction. Further, using these results, we predict hydraulic permeability during longitudinal (fiber-aligned) tensile loading, and the spatial distribution of fluid flow during osmotic loading. These results reveal novel aspects of tendon mechanics and can be used to study the physiomechanical response of tendon in response to mechanical loading.

Tendon's hierarchical structure allows for load transfer between its fibrillar elements at multiple length scales. Tendon microstructure is particularly important, because it includes the cells and their surrounding collagen fibrils, where mechanical interactions can have potentially important physiological and pathological contributions. However, the three-dimensional (3D) microstructure and the mechanisms of load transfer in that length scale are not known. It has been postulated that interfibrillar matrix shear or direct load transfer via the fusion/branching of small fibrils are responsible for load transfer, but the significance of these mechanisms is still unclear. Alternatively, the helical fibrils that occur at the microstructural scale in tendon may also mediate load transfer; however, these structures are not well studied due to the lack of a three-dimensional visualization of tendon microstructure. In this study, we used serial block-face scanning electron microscopy to investigate the 3D microstructure of fibrils in rat tail tendon. We found that tendon fibrils have a complex architecture with many helically wrapped fibrils. We studied the mechanical implications of these helical structures using finite-element modelling and found that frictional contact between helical fibrils can induce load transfer even in the absence of matrix bonding or fibril fusion/branching. This study is significant in that it provides a three-dimensional view of the tendon microstructure and suggests friction between helically wrapped fibrils as a mechanism for load transfer, which is an important aspect of tendon biomechanics.

Inelastic behaviors, such as softening, a progressive decrease in modulus before failure, occur in tendon andare important aspect in degeneration and tendinopathy. These in elastic behaviors are generally attributed to two potential mechanisms: plastic deformation and damage. However, it is not clear which is primarily responsible.In this study, we evaluated these potential mechanisms of tendon in elasticity by using a recently developed reactive in elasticity model (RIE), which is a structurally-inspired continuum mechanics frame work that models tissue in elasticity based on the molecular bond kinetics. Using RIE, we formulated two material models, one specific toplastic deformation and the other to damage. The models were independently fit to published experimental tensiletests of rat tail tendons. We quantified the inelastic effects and compared the performance of the two models infitting the mechanical response during loading, relaxation, unloading, and reloading phases. Additionally, we validated the models by using the resulting fit parameters to predict an independent set of experimental stress-straincurves from ramp-to-failure tests. Overall, the models were both successful in fitting the experiments and predicting the validation data. However, the results did not strongly favor one mechanism over the other. As a result, to distinguish between plastic deformation and damage, different experimental protocols will be needed. Nevertheless, these findings suggest the potential of RIE as a comprehensive framework for studying tendon inelastic behaviors.

Soft tissues are biopolymeric materials, primarily made of collagen and water. These tissues have non-linear, anisotropic, and inelastic mechanical behaviors that are often categorized into viscoelastic behavior, plastic deformation, and damage. While tissue’s elastic and viscoelastic mechanical properties have been measured for decades, there is no comprehensive theoretical framework for modeling inelastic behaviors of these tissues that is based on their structure. To model the three major inelastic mechanical behaviors of soft tissue we formulated a structurally inspired continuum mechanics framework based on the energy of molecular bonds that break and reform in response to external loading (reactive bonds). In this framework, we employed the theory of internal state variables and kinetics of molecular bonds. The number fraction of bonds, their reference deformation gradient, and damage parameter were used as internal state variables that allowed for consistent modeling of all three of the inelastic behaviors of tissue by using the same sets of constitutive relations. Several numerical examples are provided that address practical problems in tissue mechanics, including the difference between plastic deformation and damage. This model can be used to identify relationships between tissue’s mechanical response to external loading and its biopolymeric structure.