Subject-specific finite element model of the pelvis: development, validation and sensitivity studies

A better understanding of the three-dimensional mechanics of the pelvis, at the patient-specific level, may lead to improved treatment modalities. Although finite element (FE) models of the pelvis have been developed, validation by direct comparison with subject-specific strains has not been performed, and previous models used simplifying assumptions regarding geometry and material properties. The objectives of this study were to develop and validate a realistic FE model of the pelvis using subject-specific estimates of bone geometry, location-dependent cortical thickness and trabecular bone elastic modulus, and to assess the sensitivity of FE strain predictions to assumptions regarding cortical bone thickness as well as bone and cartilage material properties. A FE model of a cadaveric pelvis was created using subject-specific computed tomography image data. Acetabular loading was applied to the same pelvis using a prosthetic femoral stem in a fashion that could be easily duplicated in the computational model. Cortical bone strains were monitored with rosette strain gauges in ten locations on the left hemipelvis. FE strain predictions were compared directly with experimental results for validation. Overall, baseline FE predictions were strongly correlated with experimental results (r2=0.824), with a best-fit line that was not statistically different than the line y=x (experimental strains = FE predicted strains). Changes to cortical bone thickness and elastic modulus had the largest effect on cortical bone strains. The FE model was less sensitive to changes in all other parameters. The methods developed and validated in this study will be useful for creating and analyzing patient-specific FE models to better understand the biomechanics of the pelvis.     

Achilles tendon stress is more sensitive to subject-specific geometry than subject-specific material properties: A finite element analysis

This study used subject-specific measures of three-dimensional (3D) free Achilles tendon geometry in conjunction with a finite element method to investigate the effect of variation in subject-specific geometry and subject-specific material properties on tendon stress during submaximal isometric loading. Achilles tendons of eight participants (Aged 25–35 years) were scanned with freehand 3D ultrasound at rest and during a 70% maximum voluntary isometric contraction. Ultrasound images were segmented, volume rendered and transformed into subject-specific 3D finite element meshes. The mean (±SD) lengths, volumes and cross-sectional areas of the tendons at rest were 62 ± 13 mm, 3617 ± 984 mm³ and 58 ± 11 mm² respectively. The measured tendon strain at 70% MVIC was 5.9 ± 1.3%. Subject-specific material properties were obtained using an optimisation approach that minimised the difference between measured and modelled longitudinal free tendon strain. Generic geometry was represented by the average mesh and generic material properties were taken from the literature. Local stresses were subsequently computed for combinations of subject-specific and generic geometry and material properties. For a given geometry, changing from generic to subject-specific material properties had little effect on the stress distribution in the tendon. In contrast, changing from generic to subject-specific geometry had a 26-fold greater effect on tendon stress distribution. Overall, these findings indicate that the stress distribution experienced by the living free Achilles tendon of a young and healthy population during voluntary loading are more sensitive to variation in tendon geometry than variation in tendon material properties.     

Lumbar spine finite element model for healthy subjects: development and validation

Finite element (FE) method is a proven powerful and efficient tool to study the biomechanics of the human lumbar spine. However, due to the large inter-subject variability of geometries and material properties in human lumbar spines, concerns existed on the accuracy and predictive power of one single deterministic FE model with one set of spinal geometry and material properties. It was confirmed that the combined predictions (median or mean value) of several distinct FE models can be used as an improved prediction of behavior of human lumbar spine under identical loading and boundary conditions. In light of this improved prediction, five FE models (L1-L5 spinal levels) of the human lumbar spine were developed based on five healthy living subjects with identical modeling method. The five models were extensively validated through experimental and computational results in the literature. Mesh convergence and material sensitivity analysis were also conducted. We have shown that the results from the five FE models developed in this paper were consistent with the experimental data and simulation results from the existing literature. The validated modeling method introduced in this study can be used in modeling dysfunctional lumber spines such as disc degeneration and scoliosis in future work.     

Effects of idealized joint geometry on finite element predictions of cartilage contact stresses in the hip

Computational models may have the ability to quantify the relationship between hip morphology, cartilage mechanics and osteoarthritis. Most models have assumed the hip joint to be a perfect ball and socket joint and have neglected deformation at the bone-cartilage interface. The objective of this study was to analyze finite element (FE) models of hip cartilage mechanics with varying degrees of simplified geometry and a model with a rigid bone material assumption to elucidate the effects on predictions of cartilage stress. A previously validated subject-specific FE model of a cadaveric hip joint was used as the basis for the models. Geometry for the bone-cartilage interface was either: (1) subject-specific (i.e. irregular), (2) spherical, or (3) a rotational conchoid. Cartilage was assigned either a varying (irregular) or constant thickness (smoothed). Loading conditions simulated walking, stair-climbing and descending stairs. FE predictions of contact stress for the simplified models were compared with predictions from the subject-specific model. Both spheres and conchoids provided a good approximation of native hip joint geometry (average fitting error ～0.5 mm). However, models with spherical/conchoid bone geometry and smoothed articulating cartilage surfaces grossly underestimated peak and average contact pressures (50% and 25% lower, respectively) and overestimated contact area when compared to the subject-specific FE model. Models incorporating subject-specific bone geometry with smoothed articulating cartilage also underestimated pressures and predicted evenly distributed patterns of contact. The model with rigid bones predicted much higher pressures than the subject-specific model with deformable bones. The results demonstrate that simplifications to the geometry of the bone-cartilage interface, cartilage surface and bone material properties can have a dramatic effect on the predicted magnitude and distribution of cartilage contact pressures in the hip joint.     

Influence of altered geometry and material properties on tissue stress distribution under load in tendinopathic Achilles tendons – A subject-specific finite element analysis

Achilles tendon material properties and geometry are altered in Achilles tendinopathy. The purpose of this study was to determine the relative contributions of altered material properties and geometry to free Achilles tendon stress distribution during a sub-maximal contraction in tendinopathic relative to healthy tendons. Tendinopathic (n = 8) and healthy tendons (n = 8) were imaged at rest and during a sub-maximal voluntary isometric contraction using three-dimensional freehand ultrasound. Images were manually segmented and used to create subject-specific finite element models. The resting cross-sectional area of the free tendon was on average 31% greater for the tendinopathic compared to healthy tendons. Material properties for each tendon were determined using a numerical parameter optimisation approach that minimised the difference in experimentally measured longitudinal strain and the strain predicted by the finite element model under submaximal loading conditions for each tendon. The mean Young’s modulus for tendinopathic tendons was 53% lower than the corresponding control value. Finite element analyses revealed that tendinopathic tendons experience 24% less stress under the same submaximal external loading conditions compared to healthy tendons. The lower tendon stress in tendinopathy was due to a greater influence of tendon cross-sectional area, which alone reduced tendon stress by 30%, compared to a lower Young’s modulus, which alone increased tendon stress by 8%. These findings suggest that the greater tendon cross-sectional area observed in tendinopathy compensates for the substantially lower Young’s modulus, thereby protecting pathological tendon against excessive stress.     

Mesh-morphing algorithms for specimen-specific finite element modeling

Despite recent advances in software for meshing specimen-specific geometries, considerable effort is still often required to produce and analyze specimen-specific models suitable for biomechanical analysis through finite element modeling. We hypothesize that it is possible to obtain accurate models by adapting a pre-existing geometry to represent a target specimen using morphing techniques. Here we present two algorithms for morphing, automated wrapping (AW) and manual landmarks (ML), and demonstrate their use to prepare specimen-specific models of caudal rat vertebrae. We evaluate the algorithms by measuring the distance between target and morphed geometries and by comparing response to axial loading simulated with finite element (FE) methods.

First a traditional reconstruction process based on μCT was used to obtain two natural specimen-specific FE models. Next, the two morphing algorithms were used to compute mappings from the surface of one model, the source, to the other, the target, and to use this mapping to morph the source mesh to produce a target mesh. The μCT images were then used to assign element-specific material properties. In AW the mappings were obtained by wrapping the source and target surfaces with an auxiliary triangulated surface. In ML, landmarks were manually placed on corresponding locations on the surfaces of both source and target.

Both morphing algorithms were successful in reproducing the shape of the target vertebra with a median distance between natural and morphed models of 18.8 and 32.2 μm, respectively, for AW and ML. Whereas AW–morphing produced a surface more closely resembling that of the target, ML guaranteed correspondence of the landmark locations between source and target. Morphing preserved the quality of the mesh producing models suitable for FE simulation. Moreover, there were only minor differences between natural and morphed models in predictions of deformation, strain and stress. We therefore conclude that it is possible to use mesh-morphing techniques to produce accurate specimen-specific FE models of caudal rat vertebrae. Mesh morphing techniques provide advantages over conventional specimen-specific finite element modeling by reducing the effort required to generate multiple target specimen models, facilitating intermodel comparisons through correspondence of nodes and maintenance of connectivity, and lends itself to parametric evaluation of “artificial” geometries with a focus on optimizing reconstruction.     

Relation between subject-specific hip joint loading, stress distribution in the proximal femur and bone mineral density changes after total hip replacement

In the prediction of bone remodelling processes after total hip replacement (THR), modelling of the subject-specific geometry is now state-of-the-art. In this study, we demonstrate that inclusion of subject-specific loading conditions drastically influences the calculated stress distribution, and hence influences the correlation between calculated stress distributions and changes in bone mineral density (BMD) after THR.For two patients who received cementless THR, personalized finite element (FE) models of the proximal femur were generated representing the pre- and post-operative geometry. FE analyses were performed by imposing subject-specific three-dimensional hip joint contact forces as well as muscle forces calculated based on gait analysis data. Average values of the von Mises stress were calculated for relevant zones of the proximal femur. Subsequently, the load cases were interchanged and the effect on the stress distribution was evaluated. Finally, the subject-specific stress distribution was correlated to the changes in BMD at 3 and 6 months after THR.We found subject-specific differences in the stress distribution induced by specific loading conditions, as interchanging of the loading also interchanged the patterns of the stress distribution. The correlation between the calculated stress distribution and the changes in BMD were affected by the two-dimensional nature of the BMD measurement.Our results confirm the hypothesis that inclusion of subject-specific hip contact forces and muscle forces drastically influences the stress distribution in the proximal femur. In addition to patient-specific geometry, inclusion of patient-specific loading is, therefore, essential to obtain accurate input for the analysis of stress distribution after THR.     

Verification of predicted specimen-specific natural and implanted patellofemoral kinematics during simulated deep knee bend

Verified computational models represent an efficient method for studying the relationship between articular geometry, soft-tissue constraint, and patellofemoral (PF) mechanics. The current study was performed to evaluate an explicit finite element (FE) modeling approach for predicting PF kinematics in the natural and implanted knee. Experimental three-dimensional kinematic data were collected on four healthy cadaver specimens in their natural state and after total knee replacement in the Kansas knee simulator during a simulated deep knee bend activity. Specimen-specific FE models were created from medical images and CAD implant geometry, and included soft-tissue structures representing medial–lateral PF ligaments and the quadriceps tendon. Measured quadriceps loads and prescribed tibiofemoral kinematics were used to predict dynamic kinematics of an isolated PF joint between 10° and 110° femoral flexion. Model sensitivity analyses were performed to determine the effect of rigid or deformable patellar representations and perturbed PF ligament mechanical properties (pre-tension and stiffness) on model predictions and computational efficiency.Predicted PF kinematics from the deformable analyses showed average root mean square (RMS) differences for the natural and implanted states of less than 3.1° and 1.7 mm for all rotations and translations. Kinematic predictions with rigid bodies increased average RMS values slightly to 3.7° and 1.9 mm with a five-fold decrease in computational time. Two-fold increases and decreases in PF ligament initial strain and linear stiffness were found to most adversely affect kinematic predictions for flexion, internal–external tilt and inferior–superior translation in both natural and implanted states. The verified models could be used to further investigate the effects of component alignment or soft-tissue variability on natural and implant PF mechanics.     

Anisometry Anterior Cruciate Ligament Sport Injury Mechanism Study: A Finite Element Model with Optimization Method

ACL damage is one the most frequent causes of knee injuries and thus has long been the focus of research in biomechanics and sports medicine. Due to the anisometric geometry and functional complexity of the ACL in the knee joint, it is usually difficult to experimentally study the biomechanics of ACLs. Anatomically ACL geometry was obtained from both MR images and anatomical observations. The optimal material parameters of the ACL were obtained by using an optimization-based material identification method that minimized the differences between experimental results from ACL specimens and FE simulations. The optimal FE model simulated biomechanical responses of the ACL during complex combined injury-causing knee movements, it predicted stress concentrations on the top and middle side of the posterolateral (PL) bundles. This model was further validated by a clinical case of ACL injury diagnosed by MRI and arthroscope, it demonstrated that the locations of rupture in the patient’s knee corresponded to those where the stresses and moments were predicted to be concentrated. The result implies that varus rotation played a contributing but secondary role in injury under combined movements, the ACL elevation angle, is positive correlated with the tensional loading tolerance of the ACL.     

Growing multiblock structures: a semi-automated approach to block placement for multiblock hexahedral meshing

Finite element (FE) analysis is a cornerstone of orthopaedic biomechanics research. Three-dimensional medical imaging provides sufficient resolution for the subject-specific FE models to be generated from these data-sets. FE model development requires discretisation of a three-dimensional domain, which can be the most time-consuming component of a FE study. Hexahedral meshing tools based on the multiblock method currently rely on the manual placement of building blocks for mesh generation. We hypothesise that angular analysis of the geometric centreline for a three-dimensional surface could be used to automatically generate building block structures for the multiblock hexahedral mesh generation. Our algorithm uses a set of user-defined points and parameters to automatically generate a multiblock structure based on a surface's geometric centreline. This significantly reduces the time required for model development. We have applied this algorithm to 47 bones of varying geometries and successfully generated a FE mesh in all cases. This work represents significant advancement in automatically generating multiblock structures for a wide range of geometries.     

Development and evaluation of multiple tendon injury models in the mouse

The mouse has proven to be an advantageous animal model system in basic science research focused on aiding in development and evaluation of potential treatments; however, the small size of mouse tendons makes consistent and reproducible injury models and subsequent biomechanical evaluation challenging for studying tendon healing. In this study, we investigated the feasibility and reproducibility of multiple mouse tendon injury models. Our hypothesis was that incisional (using a blade) and excisional (using a biopsy punch) injuries would result in consistent differences in tendon material properties. At 16 weeks of age, 17 C57BL/6 mice underwent surgery to create defects in the flexor digitorum longus, Achilles, or patellar tendon. Each animal received 1-2 full-thickness, central-width incisional or excisional injuries per limb; at least one tendon per limb remained uninjured. The injuries were distributed such that each tendon type had comparable numbers of uninjured, incisionally injured, and excisionally injured specimens. Three weeks after injury, all animals were euthanized and tendons were harvested for mechanical testing. As hypothesized, differences were detected for all three different tendon types at three weeks post-injury. While all models created injuries that produced predictable outcomes, the patellar tendon model was the most consistent in terms of number and size of significant differences in injured tendons compared to native properties, as well as in the overall variance in the data. This finding provides support for its use in fundamental tendon healing studies; however, future work may use any of these models, based on their appropriateness for the specific question under study.     

Development and validation of subject-specific finite element models for blunt trauma study

This study developed and validated finite element (FE) models of swine and human thoraxes and abdomens that had subject-specific anatomies and could accurately and efficiently predict body responses to blunt impacts. Anatomies of the rib cage, torso walls, thoracic, and abdominal organs were reconstructed from X-ray computed tomography (CT) images and extracted into geometries to build FE meshes. The rib cage was modeled as an inhomogeneous beam structure with geometry and bone material parameters determined directly from CT images. Meshes of soft components were generated by mapping structured mesh templates representative of organ topologies onto the geometries. The swine models were developed from and validated by 30 animal tests in which blunt insults were applied to swine subjects and CT images, chest wall motions, lung pressures, and pathological data were acquired. A comparison of the FE calculations of animal responses and experimental measurements showed a good agreement. The errors in calculated response time traces were within 10% for most tests. Calculated peak responses showed strong correlations with the experimental values. The stress concentration inside the ribs, lungs, and livers produced by FE simulations also compared favorably to the injury locations. A human FE model was developed from CT images from the Visible Human project and was scaled to simulate historical frontal and side post mortem human subject (PMHS) impact tests. The calculated chest deformation also showed a good agreement with the measurements. The models developed in this study can be of great value for studying blunt thoracic and abdominal trauma and for designing injury prevention techniques, equipments, and devices.     

Linear and nonlinear analyses of femoral fractures: Computational/experimental study

This paper describes two new methods for computational fracture analysis of human femur using Quantitative Computed Tomography (QCT) voxel-based finite element (FE) simulation. The paper also reports comprehensive mechanical testing for validation of the methods and evaluation of the required material properties. The analyses and tests were carried out on 15 human femurs under 11 different stance-type loading orientations. Several classical forms of subcapital, transcervical, basicervical, and intertrochanteric fractures plus a specific type of subtrochanteric fracture were created and analyzed. A new procedure was developed for prediction of the strengths and the fracture initiation patterns using a FE-based linear scheme. The predicted and observed fracture patterns were in correspondence, and the FE predictions of the fracture loads were in very good agreement with the experimental results. Moreover, the crack initiation and growth behaviors of two subtrochanteric fractures were successfully simulated through a novel implementation of the cohesive zone model (CZM) within a nonlinear FE analysis scheme. The CZM parameters were obtained through a series of experimental tests on different types of specimens and determination of a variety of material properties for different anatomic regions and orientations. The presented results indicated that the locations and patterns of crack initiation, the sequences of crack growth on different paths, and the compatibility of growth increments agreed very well with the observed specifications. Also, very good agreements were achieved between the measured and simulated fracture loads.     

On the importance of 3D,geometrically accurate,and subject-specific finite element analysis for evaluation of in-vivo soft tissue loads

Pressure ulcers are a type of local soft tissue injury due to sustained mechanical loading and remain a common issue in patient care. People with spinal cord injury (SCI) are especially at risk of pressure ulcers due to impaired mobility and sensory perception. The development of load improving support structures relies on realistic tissue load evaluation e.g. using finite element analysis (FEA). FEA requires realistic subject-specific mechanical properties and geometries. This study focuses on the effect of geometry. MRI is used for the creation of geometrically accurate models of the human buttock for three able-bodied volunteers and three volunteers with SCI. The effect of geometry on observed internal tissue deformations for each subject is studied by comparing FEA findings for equivalent loading conditions. The large variations found between subjects confirms the importance of subject-specific FEA.     

Aging leads to inferior Achilles tendon mechanics and altered ankle function in rodents

Spontaneous rupture of the Achilles tendon is increasingly common in the middle aged population. However, the cause for the particularly high incidence of injury in this age group is not well understood. Therefore, the objective of this study was to identify age-specific differences in the Achilles tendon-muscle complex using an animal model. Functional measures were performed in vivo and tissues were harvested following euthanasia for mechanical, structural, and histological analysis from young, middle aged, and old rats. Numerous alterations in tendon properties were detected across age groups, including inferior material properties (maximum stress, modulus) with increasing age. Differences in function were also observed, as older animals exhibited increased ankle joint passive stiffness and decreased propulsion force during locomotion. Macroscale differences in tendon organization were not observed, although cell density and nuclear shape did vary between age groups. Muscle fiber size and type distribution were not notably affected by age, indicating that other factors may be more responsible for age-specific Achilles tendon rupture rates. This study improves our understanding of the role of aging in Achilles tendon biomechanics and ankle function, and helps provide a potential explanation for the disparate incidence of Achilles tendon ruptures in varying age groups.     

Subject-specific tendon-aponeurosis definition in hill-type model predicts higher muscle forces in dynamic tasks

Neuromusculoskeletal models are a common method to estimate muscle forces. Developing accurate neuromusculoskeletal models is a challenging task due to the complexity of the system and large inter-subject variability. The estimation of muscles force is based on the mechanical properties of tendon-aponeurosis complex. Most neuromusculoskeletal models use a generic definition of the tendon-aponeurosis complex based on in vitro test, perhaps limiting their validity. Ultrasonography allows subject-specific estimates of the tendon-aponeurosis complex's mechanical properties. The aim of this study was to investigate the influence of subject-specific mechanical properties of the tendon-aponeurosis complex on a neuromusculoskeletal model of the ankle joint. Seven subjects performed isometric contractions from which the tendon-aponeurosis force-strain relationship was estimated. Hopping and running tasks were performed and muscle forces were estimated using subject-specific tendon-aponeurosis and generic tendon properties. Two ultrasound probes positioned over the muscle-tendon junction and the mid-belly were combined with motion capture to estimate the in vivo tendon and aponeurosis strain of the medial head of gastrocnemius muscle. The tendon-aponeurosis force-strain relationship was scaled for the other ankle muscles based on tendon and aponeurosis length of each muscle measured by ultrasonography. The EMG-driven model was calibrated twice - using the generic tendon definition and a subject-specific tendon-aponeurosis force-strain definition. The use of subject-specific tendon-aponeurosis definition leads to a higher muscle force estimate for the soleus muscle and the plantar-flexor group, and to a better model prediction of the ankle joint moment compared to the model estimate which used a generic definition. Furthermore, the subject-specific tendon-aponeurosis definition leads to a decoupling behaviour between the muscle fibre and muscle-tendon unit in agreement with previous experiments using ultrasonography. These results indicate the use of subject-specific tendon-aponeurosis definitions in a neuromusculoskeletal model produce better agreement with measured external loads and more physiological model behaviour.     

The contribution of the glenoid labrum to glenohumeral stability under physiological joint loading using finite element analysis

The labrum contributes to passive glenohumeral joint stability. Cadaveric studies have demonstrated that this has position and load dependency, which has not been quantified under physiological loads. This study aims to validate subject-specific finite element (FE) models against in vitro measurements of joint stability and to utilise the FE models to predict joint stability under physiological loads. The predicted stability values were within ± one standard deviation of experimental data and the FE models showed a reduction in stability of 10–15% with high, physiological, loads. The developed regression equations provide the first representation of passive glenohumeral stability and will aid surgical decision-making.     

Finite element modelling approaches for well-ordered porous metallic materials for orthopaedic applications: cost effectiveness and geometrical considerations

The mechanical properties of well-ordered porous materials are related to their geometrical parameters at the mesoscale. Finite element (FE) analysis is a powerful tool to design well-ordered porous materials by analysing the mechanical behaviour. However, FE models are often computationally expensive. This article aims to develop a cost-effective FE model to simulate well-ordered porous metallic materials for orthopaedic applications. Solid and beam FE modelling approaches are compared, using finite size and infinite media models considering cubic unit cell geometry. The model is then applied to compare two unit cell geometries: cubic and diamond. Models having finite size provide similar results than the infinite media model approach for large sample sizes. In addition, these finite size models also capture the influence of the boundary conditions on the mechanical response for small sample sizes. The beam FE modelling approach showed little computational cost and similar results to the solid FE modelling approach. Diamond unit cell geometry appeared to be more suitable for orthopaedic applications than the cubic unit cell geometry.     

Validation of an efficient method of assigning material properties in finite element analysis of pelvic bone

Bone in the pelvis is a composite material with a complex anatomical structure that is difficult to model computationally. Rather than assigning material properties to increasingly smaller elements to capture detail in three-dimensional finite element (FE) models, properties can be assigned to Gauss points within larger elements. As part of a validation process, we compared experimental and analytical results from a composite beam under four-point load to FE models with material properties assigned to refined elements and Gauss points within larger elements. Both FE models accurately predicted deformation and the analytical predictions of internal shear stress.     

Biomechanical and structural response of healing Achilles tendon to fatigue loading following acute injury

Achilles tendon injuries affect both athletes and the general population, and their incidence is rising. In particular, the Achilles tendon is subject to dynamic loading at or near failure loads during activity, and fatigue induced damage is likely a contributing factor to ultimate tendon failure. Unfortunately, little is known about how injured Achilles tendons respond mechanically and structurally to fatigue loading during healing. Knowledge of these properties remains critical to best evaluate tendon damage induction and the ability of the tendon to maintain mechanical properties with repeated loading. Thus, this study investigated the mechanical and structural changes in healing mouse Achilles tendons during fatigue loading. Twenty four mice received bilateral full thickness, partial width excisional injuries to their Achilles tendons (IACUC approved) and twelve tendons from six uninjured mice were used as controls. Tendons were fatigue loaded to assess mechanical and structural properties simultaneously after 0, 1, 3, and 6 weeks of healing using an integrated polarized light system. Results showed that the number of cycles to failure decreased dramatically (37-fold, p<0.005) due to injury, but increased throughout healing, ultimately recovering after 6 weeks. The tangent stiffness, hysteresis, and dynamic modulus did not improve with healing (p<0.005). Linear regression analysis was used to determine relationships between mechanical and structural properties. Of tendon structural properties, the apparent birefringence was able to best predict dynamic modulus (R²=0.88–0.92) throughout healing and fatigue life. This study reinforces the concept that fatigue loading is a sensitive metric to assess tendon healing and demonstrates potential structural metrics to predict mechanical properties.