Overload in a Rat In Vivo Model of Synergist Ablation Induces Tendon Multiscale Structural and Functional Degeneration

Tendon degeneration is typically described as an overuse injury with little distinction made between magnitude of load (overload) and number of cycles (overuse). Further, in vivo, animal models of tendon degeneration are mostly overuse models, where tendon damage is caused by a high number of load cycles. As a result, there is a lack of knowledge of how isolated overload leads to degeneration in tendons. A surgical model of synergist ablation (SynAb) overloads the target tendon, plantaris, by ablating its synergist tendon, Achilles. The objective of this study was to evaluate the structural and functional changes that occur following overload of plantaris tendon in a rat SynAb model. Tendon cross-sectional area (CSA) and shape changes were evaluated by longitudinal MR imaging up to 8 weeks postsurgery. Tissue-scale structural changes were evaluated by semiquantified histology and second harmonic generation microscopy. Fibril level changes were evaluated with serial block face scanning electron microscopy (SBF-SEM). Functional changes were evaluated using tension tests at the tissue and microscale using a custom testing system allowing both video and microscopy imaging. At 8 weeks, overloaded plantaris tendons exhibited degenerative changes including increases in CSA, cell density, collagen damage area fraction (DAF), and fibril diameter, and decreases in collagen alignment, modulus, and yield stress. To interpret the differences between overload and overuse in tendon, we introduce a new framework for tendon remodeling and degeneration that differentiates between the inputs of overload and overuse. In summary, isolated overload induces multiscale degenerative structural and functional changes in plantaris tendon.

Volume Loss and Recovery in Bovine Knee Meniscus Loaded in Circumferential Tension

Load-induced volume change is an important aspect of knee meniscus function because volume loss creates fluid pressure, which minimizes friction and helps support compressive loads. The knee meniscus is unusual amongst cartilaginous tissues in that it is loaded not only in axial compression, but also in circumferential tension between its tibial attachments. Despite the physiologic importance of the knee meniscus' tensile properties, its volumetric strain in tension has never been directly measured, and predictions of volume strain in the scientific literature are inconsistent. In this study, we apply uniaxial tension to bovine knee meniscus and use biplanar imaging to directly observe the resulting three-dimensional volume change and unloaded recovery, revealing that tension causes volumetric contraction. Compression is already known to also cause contraction; therefore, all major physiologic loads compress and pressurize the meniscus, inducing fluid outflow. Although passive unloaded recovery is often described as slow relative to loaded loss, here we show that at physiologic strains the volume recovery rate in the meniscus upon unloading is faster than the rate of volume loss. These measurements of volumetric strain are an important step toward a complete theory of knee meniscus fluid flow and load support.

Tuning the Structure of Nylon 6,6 Electrospun Bundles to Mimic the Mechanical Performance of Tendon Fascicles

Tendon and ligament injuries are triggered by mechanical loading, but the specific mechanisms are not yet clearly identified. It is well established however, that the inflection and transition points in tendon stress-strain curves represent thresholds that may signal the onset of irreversible fibrillar sliding. This phenomenon often results in a progressive macroscopic failure of these tissues. With the aim to simulate and replace tendons, electrospinning has been demonstrated to be a suitable technology to produce nanofibers similar to the collagen fibrils in a mat form. These nanofibrous mats can be easily assembled in higher hierarchical levels to reproduce the whole tissue structure. Despite the fact that several groups have developed electrospun tendon-inspired structures, an investigation of the inflection and transition point mechanics is missing. Comparing their behavior with that of the natural counterpart is important to adequately replicate their behavior at physiological strain levels. To fill this gap, in this work fascicle-inspired electrospun nylon 6,6 bundles were produced with different collector peripheral speeds (i.e., 19.7 m s^–1; 13.7 m s^–1; 7.9 m s^–1), obtaining different patterns of nanofibers alignment. The scanning electron microcopy revealed a fibril-inspired structure of the nanofibers with an orientation at the higher speed similar to those in tendons and ligaments (T/L). A tensile mechanical characterization was carried out showing an elastic-brittle biomimetic behavior for the higher speed bundles with a progressively more ductile behavior at slower speeds. Moreover, for each sample category the transition and the inflection points were defined to study how these points can shift with the nanofiber arrangement and to compare their values with those of tendons. The results of this study will be of extreme interest for the material scientists working in the field, to model and improve the design of their electrospun structures and scaffolds and enable building a new generation of artificial tendons and ligaments.

Identifiability of tissue material parameters from uniaxial tests using multi-start optimization

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.

Evaluation of transverse poroelastic mechanics of tendon using osmotic loading and biphasic mixture finite element modeling

Tendon's viscoelastic behaviors are important to the tissue mechanical function and cellular mechanobiology. When loaded in longitudinal tension, tendons often have a 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 using biphasic finite element modeling. Our tendon samples had a transverse hydraulic permeability of 10^-4 to 10^-5 mm⁴. (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.

Evaluating Plastic Deformation and Damage as Potential Mechanisms for Tendon Inelasticity using a Reactive Modeling Framework

Inelastic behaviors, such as softening, a progressive decrease in modulus before failure, occur in tendon and are important aspect in degeneration and tendinopathy. These inelastic 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 inelasticity by using a recently developed reactive inelasticity model (RIE), which is a structurally-inspired continuum mechanics framework that models tissue inelasticity based on the molecular bond kinetics. Using RIE, we formulated two material models, one specific to plastic deformation and the other to damage. The models were independently fit to published experimental tensile tests of rat tail tendons. We quantified the inelastic effects and compared the performance of the two models in fitting 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-strain curves 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.

A Reactive Inelasticity Theoretical Framework for Modeling Viscoelasticity, Plastic Deformation, and Damage in Fibrous Soft Tissue

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.

Short cracks in knee meniscus tissue cause strain concentrations, but do not reduce ultimate stress, in single-cycle uniaxial tension

Tears are central to knee meniscus pathology and, from a mechanical perspective, are crack-like defects (cracks). In many materials, cracks create stress concentrations that cause progressive local rupture and reduce effective strength. It is currently unknown if cracks in meniscus have these consequences; if they do, this would have repercussions for management of meniscus pathology. The objective of this study was to determine if a short crack in meniscus tissue, which mimics a preclinical meniscus tear, (a) causes crack growth and reduces effective strength, (b) creates a near-tip strain concentration and (c) creates unloaded regions on either side of the crack. Specimens with and without cracks were tested in uniaxial tension and compared in terms of macroscopic stress-strain curves and digital image correlation strain fields. The strain fields were used as an indicator of stress concentrations and unloaded regions. Effective strength was found to be insensitive to the presence of a crack (potential effect < 0.86 s.d.; β = 0.2), but significant strain concentrations, which have the potential to lead to long-term accumulation of tissue or cell damage, were observed near the crack tip.

Investigating mechanisms of tendon damage by measuring multi-scale recovery following tensile loading

Tendon pathology is associated with damage. While tendon damage is likely initiated by mechanical loading, little is known about the specific etiology. Damage is defined as an irreversible change in the microstructure that alters the macroscopic mechanical parameters. In tendon, the link between mechanical loading and microstructural damage, resulting in macroscopic changes, is not fully elucidated. In addition, tendon damage at the macroscale has been proposed to initiate when tendon is loaded beyond a strain threshold, yet the metrics to define the damage threshold are not determined. We conducted multi-scale mechanical testing to investigate the mechanism of tendon damage by simultaneously quantifying macroscale mechanical and microstructural changes. At the microscale, we observe full recovery of the fibril strain and only partial recovery of the interfibrillar sliding, indicating that the damage initiates at the interfibrillar structures. We show that non-recoverable sliding is a mechanism for tendon damage and is responsible for the macroscale decreased linear modulus and elongated toe-region observed at the fascicle-level, and these macroscale properties are appropriate metrics that reflect tendon damage. We concluded that the inflection point of the stress-strain curve represents the damage threshold and, therefore, may be a useful parameter for future studies. Establishing the mechanism of damage at multiple length scales can improve prevention and rehabilitation strategies for tendon pathology.

STATEMENT OF SIGNIFICANCE

Tendon pathology is associated with mechanically induced damage. Damage, as defined in engineering, is an irreversible change in microstructure that alters the macroscopic mechanical properties. Although microstructural damage and changes to macroscale mechanics are likely, this link to microstructural change was not yet established. We conducted multiscale mechanical testing to investigate the mechanism of tendon damage by simultaneously quantifying macroscale mechanical and microstructural changes. We showed that non-recoverable sliding between collagen fibrils is a mechanism for tendon damage. Establishing the mechanism of damage at multiple length scales can improve prevention and rehabilitation strategies for tendon pathology.

Advances in Quantification of Meniscus Tensile Mechanics Including Nonlinearity, Yield, and Failure

The meniscus provides crucial knee function and damage to it leads to osteoarthritis of the articular cartilage. Accurate measurement of its mechanical properties is therefore important, but there is uncertainty about how the test procedure affects the results, and some key mechanical properties are reported using ad hoc criteria (modulus) or not reported at all (yield). This study quantifies the meniscus' stress-strain curve in circumferential and radial uniaxial tension. A fiber recruitment model was used to represent the toe region of the stress-strain curve, and new reproducible and objective procedures were implemented for identifying the yield point and measuring the elastic modulus. Patterns of strain heterogeneity were identified using strain field measurements. To resolve uncertainty regarding whether rupture location (i.e., midsubstance rupture versus at-grip rupture) influences the measured mechanical properties, types of rupture were classified in detail and compared. Dogbone (DB)-shaped specimens are often used to promote midsubstance rupture; to determine if this is effective, we compared DB and rectangle (R) specimens in both the radial and circumferential directions. In circumferential testing, we also compared expanded tab (ET) specimens under the hypothesis that this shape would more effectively secure the meniscus' curved fibers and thus produce a stiffer response. The fiber recruitment model produced excellent fits to the data. Full fiber recruitment occurred approximately at the yield point, strongly supporting the model's physical interpretation. The strain fields, especially shear and transverse strain, were extremely heterogeneous. The shear strain field was arranged in pronounced bands of alternating positive and negative strain in a pattern similar to the fascicle structure. The site and extent of failure showed great variation, but did not affect the measured mechanical properties. In circumferential tension, ET specimens underwent earlier and more rapid fiber recruitment, had less stretch at yield, and had greater elastic modulus and peak stress. No significant differences were observed between R and DB specimens in either circumferential or radial tension. Based on these results, ET specimens are recommended for circumferential tests and R specimens for radial tests. In addition to the data obtained, the procedural and modeling advances made in this study are a significant step forward for meniscus research and are applicable to other fibrous soft tissues.

Determining Tension-Compression Nonlinear Mechanical Properties of Articular Cartilage from Indentation Testing

The indentation test is widely used to determine the in situ biomechanical properties of articular cartilage. The mechanical parameters estimated from the test depend on the constitutive model adopted to analyze the data. Similar to most connective tissues, the solid matrix of cartilage displays different mechanical properties under tension and compression, termed tension-compression nonlinearity (TCN). In this study, cartilage was modeled as a porous elastic material with either a conewise linear elastic matrix with cubic symmetry or a solid matrix reinforced by a continuous fiber distribution. Both models are commonly used to describe the TCN of cartilage. The roles of each mechanical property in determining the indentation response of cartilage were identified by finite element simulation. Under constant loading, the equilibrium deformation of cartilage is mainly dependent on the compressive modulus, while the initial transient creep behavior is largely regulated by the tensile stiffness. More importantly, altering the permeability does not change the shape of the indentation creep curves, but introduces a parallel shift along the horizontal direction on a logarithmic time scale. Based on these findings, a highly efficient curve-fitting algorithm was designed, which can uniquely determine the three major mechanical properties of cartilage (compressive modulus, tensile modulus, and permeability) from a single indentation test. The new technique was tested on adult bovine knee cartilage and compared with results from the classic biphasic linear elastic curve-fitting program.

WEAR BEHAVIOR OF CARBON NANOTUBE/HIGH DENSITY POLYETHYLENE COMPOSITES

Carbon Nanotube/High Density Polyethylene (CNT/HDPE) composites were manufactured and tested to determine their wear behavior. The nanocomposites were made from untreated multi-walled carbon nanotubes and HDPE pellets. Thin films of the precursor materials were created with varying weight percentages of nanotubes (1%, 3%, and 5%), through a process of mixing and extruding. The precursor composites were then molded and machined to create test specimens for mechanical and wear tests. These included small punch testing to compare stiffness, maximum load and work-to-failure and block-on-ring testing to determine wear behavior. Each of the tests was conducted for the different weight percentages of composite as well as pure HDPE as the baseline. The measured mechanical properties and wear resistance of the composite materials increased with increasing nanotube content in the range studied.

Analysis of a femoral hip prosthesis designed to reduce stress shielding

The natural stress distribution in the femur is significantly altered after total hip arthroplasty (THA). When an implant is introduced, it will carry a portion of the load, causing a reduction of stress in some regions of the remaining bone. This phenomenon is commonly known as stress shielding. In response to the changed mechanical environment the shielded bone will remodel according to Wolff's law, resulting in a loss of bone mass through the biological process called resorption. Resorption can, in turn, cause or contribute to loosening of the prosthesis. The problem is particularly common among younger THA recipients. This study explores the hypothesis that through redesign, a total hip prosthesis can be developed to substantially reduce stress shielding. First, we describe the development of a new femoral hip prosthesis designed to alleviate this problem through a new geometry and system of proximal fixation. A numerical comparison with a conventional intramedullary prosthesis as well as another proximally fixed prosthesis, recently developed by Munting and Verhelpen (1995. Journal of Biomechanics 28(8), 949-961) is presented. The results show that the new design produces a more physiological stress state in the proximal femur.

Modeling of airflow in the pharynx with application to sleep apnea

A three-dimensional numerical modeling of airflow in the human pharynx using an anatomically accurate model was conducted. The pharynx walls were assumed to be passive and rigid. The results showed that the pressure drop in the pharynx lies in the range 200-500 Pa. The onset of turbulence was found to increase the pressure drop by 40 percent. A wide range of pharynx geometries covering three sleep apnea treatment therapies (CPAP, mandibular repositioning devices, and surgery) were modeled and the resulting flow characteristics were investigated and compared. The results confirmed that the airflow in the pharynx lies in the laminar-to-turbulence transitional flow regime and thus, a subtle change in the morphology caused by these treatment therapies can significantly affect the airflow characteristics.

Failure properties of suture anchors in the glenoid and the effects of cortical thickness

The ultimate pullout strength and fatigue properties of a screw-design suture anchor implanted in the anterior glenoid rim were investigated and compared with results from a nonscrew-design suture anchor. Twenty-two cadaveric glenoids were harvested and one to two anchors were implanted in the superior and inferior quadrants. Fifty-seven Statak 3.5 anchors (Zimmer, Warsaw, IN) were tested and compared with results obtained in a previous study on 50 Mitek GII anchors (Mitek Products, Inc, Westwood, MA). The specimens were mounted on an Instron fatigue testing machine (Instron Corp, Canton, MA) and cycled between preselected minimum and maximum loads until pullout. The Mitek GII maintained a higher pullout strength than the Statak 3.5 after cyclic loading. Cortical thickness at the implantation sites was measured, and found to decrease monotonically from superior to inferior positions. The ultimate pullout strength, and subsequently the fatigue life, of both types of suture anchors depended directly on cortical thickness. The significantly lower performance of both anchors when placed inferiorly emphasizes the importance of correct anchor selection, number, and placement in this region. All anchors settled during the first 10 to 100 cycles, resulting in partial exposure of the implant. Intraoperative cycling of the anchors before suture tying may be necessary to achieve complete settling and prevent subsequent loss of coaptation between capsule and glenoid. The study shows that for the anchors to last 1,000 cycles or more, less than 50% of the theoretical ultimate pullout strength should be applied cyclically. With aggressive early rehabilitation exercises, this significant decrease in fixation strength could shift reconstruction failure from suture breakage or soft tissue tearing to anchor pullout.

Fatigue properties of suture anchors in anterior shoulder reconstructions: Mitek GII

Suture anchors have simplified anterior capsule labral reconstruction. During rehabilitation the shoulder goes through many repetitions of range of motion exercises. These exercises will repetitively submaximally load the anchor and in theory should reduce the pullout strength of the suture anchor. No published reports exist on the fatigue strengths and properties of one of the most commonly used anchors: Mitek GII suture anchors. Fifty trials of cyclic submaximal load were done on 22 cadaveric glenoids with an average age of 66.8 years (range, 40 to 90 years). At two to three different sites on the same specimen, the anchors were inserted according to manufacturer's specifications. The anchors were tested to failure on a Instron 1331 servohydraulic mechanical testing system at 2 Hertz sinusoidal loading pattern using steel sutures and a predetermined load. There were 22 (44%) tests performed in the superior quadrant and 28 (56%) tests in the inferior quadrant. All anchors pulled out, and no wires broke. There were statistically significant differences between the superior and inferior portion of the glenoid with regard to number of cycles to failure at a given maximum load. The anchors underwent an average of 6,220 cycles before pullout at an average load of 162 N (SD = 73 N). In the superior quadrant, the average ultimate pullout strength was 237 N (SD = 42 N), whereas in the inferior quadrant the average ultimate pullout strength was 126 N (SD = 36 N). Hence, the ultimate pullout strength of the Mitek GII anchor was significantly higher (P < .002) in the superior quadrant than in the inferior quadrant. Using a least squares regression analysis, it was possible to predict the fatigue life of the superiorly and inferiorly placed suture anchors over a wide range of cycles. The R-squared values for trendlines showed good reliability (superior R2 = 0.55; inferior R2 = 0.28). The fatigue life curves for the two different quadrants were normalized using the ultimate pullout strength. This new, universal curve predicts the fatigue life of the Mitek GII anchor as a percentage of the ultimate pullout strength for any selected location. For a clinically relevant number of cycles, no more than approximately 40% to 50% of the ultimate pullout strength of the suture anchor can be cyclically applied to the anchor to guarantee a life for the duration of rehabilitation. For the entire system, the inferiorly placed anchors dictate the amount of cyclically applied load the system can experience without failing, and rehabilitation should be adjusted accordingly.

Cracks emanating from a fluid filled void loaded in compression: application to the bone-implant interface

Loosening of orthopedic implants is believed to be caused, in part, by fracture at the bone-cement interface. This loosening occurs even in regions where the interfacial load is primarily compressive. A model is developed whereby cracks can radiate from an elliptical fluid filled void. The incompressible fluid is allowed to penetrate into the cracks when the system is loaded compressively. The mode I stress intensity factor is calculated to test the feasibility of crack growth, and a numerical scheme which uses piecewise quadratic polynomials is used to solve the resulting singular integral equations. The results show the combinations of parameters for which cracks are likely to grow.