Fiber-based modeling of in situ ankle ligaments with consideration of progressive failure.

Ligament sprains account for a majority of injuries to the foot and ankle complex among athletic populations. The infeasibility of measuring the in situ response and load paths of individual ligaments has precluded a complete characterization of their mechanical behavior via experiment. In the present study a fiber-based modeling approach of in situ ankle ligaments was developed and validated for determining the heterogeneous force-elongation characteristics and the consequent injury patterns. Nine major ankle ligaments were modeled as bundles of discrete elements, corresponding functionally to the structure of collagen fibers. To incorporate the progressive nature of ligamentous injury, the limit strain at the occurrence of fiber failure was described by a distribution function ranging from 12% to 18% along the width of the insertion site. The model was validated by comparing the structural kinetic and kinematic response obtained experimentally and computationally under well-controlled foot rotations. The simulation results replicated the 6 degree-of-freedom bony motion and ligamentous injuries and, by implication, the in situ deformations of the ligaments. Gross stiffness of the whole ligament derived from the fibers was comparable to existing experimental data. The present modeling approach provides a biomechanically realistic, interpretable and computationally efficient way to characterize the in situ ligament slack, sequential and heterogeneous uncrimping of collagen fascicles and failure propagation as the external load is applied. Applications of this model include functional ankle joint mechanics, injury prevention and countermeasure design for athletes.

Searching for the "sweet spot": the foot rotation and parallel engagement of ankle ligaments in maximizing injury tolerance.

Ligament sprains, defined as tearing of bands of fibrous tissues within ligaments, account for a majority of injuries to the foot and ankle complex in field-based sports. External rotation of the foot is considered the primary injury mechanism of syndesmotic ankle sprains with concomitant flexion and inversion/eversion associated with particular patterns of ligament trauma. However, the influence of the magnitude and direction of loading vectors to the ankle on the in situ stress state of the ligaments has not been quantified in the literature. The objective of the present study was to search for the maximum injury tolerance of a human foot with an acceptable subfailure distribution of individual ligaments. We used a previously developed and comprehensively validated foot and ankle model to reproduce a range of combined foot rotation experienced during high-risk sports activities. Biomechanical computational investigation was performed on initial foot rotation from [Formula: see text] of plantar flexion to [Formula: see text] of dorsiflexion, and from [Formula: see text] of inversion to [Formula: see text] of eversion prior to external rotation. Change in initial foot rotation shifted injury initiation among different ligaments and resulted in a wide range of injury tolerances at the structural level (e.g., 36-125 Nm of rotational moment). The observed trend was in agreement with a parallel experimental study that initial plantar flexion decreased the incidence of syndesmotic injury compared to a neutral foot. A mechanism of distributing even loads across ligaments subjected to combined foot rotations was identified. This mechanism is potential to obtain the maximum load-bearing capability of a foot and ankle while minimizing the injury severity of ligaments. Such improved understanding of ligament injuries in athletes is necessary to facilitate injury management by clinicians and countermeasure development by biomechanists.

Transient and long-time kinetic responses of the cadaveric leg during internal and external foot rotation.

The purpose of this study was to determine the long-time and transient characteristics of the moment generated by external (ER) and internal (IR) rotation of the calcaneus with respect to the tibia. Two human cadaver legs were disarticulated at the knee joint while maintaining the connective tissue between the tibia and fibula. An axial rotation of 21° was applied to the proximal tibia to generate either ER or IR while the fibula was unconstrained and the calcaneus was permitted to translate in the transverse plane. These boundary conditions were intended to allow natural motion of the fibula and for the effective applied axis of rotation to move relative to the ankle and subtalar joints based on natural articular motions among the tibia, fibula, talus, and calcaneus. A load cell at the proximal tibia measured all components of force and moment. A quasi-linear model of the moment along the tibia axis was developed to determine the transient and long-time loads generated by this ER/IR. Initially neutral, everted, inverted, dorsiflexed, and plantarflexed foot orientations were tested. For the neutral position, the transient elastic moment was 16.5N-m for one specimen and 30.3N-m for the other in ER with 26.3 and 32.1N-m in IR. The long-time moments were 5.5 and 13.2N-m (ER) and 9.0 and 9.5N-m (IR). These loads were found to be transient over time similar to previous studies on other biological structures where the moment relaxed as time progressed after the initial ramp in rotation.

Determination of the in situ mechanical behavior of ankle ligaments.

The mechanical behavior of ankle ligaments at the structural level can be characterized by force-displacement curves in the physiologic phase up to the initiation of failure. However, these properties are difficult to characterize in vitro due to the experimental difficulties in replicating the complex geometry and non-uniformity of the loading state in situ. This study used a finite element parametric modeling approach to determine the in situ mechanical behavior of ankle ligaments at neutral foot position for a mid-sized adult foot from experimental derived bony kinematics. Nine major ankle ligaments were represented as a group of fibers, with the force-elongation behavior of each fiber element characterized by a zero-force region and a region of constant stiffness. The zero-force region, representing the initial tension or slackness of the whole ligament and the progressive fiber uncrimping, was identified against a series of quasi-static experiments of single foot motion using simultaneous optimization. A range of 0.33-3.84mm of the zero-force region was obtained, accounting for a relative length of 6.7±3.9%. The posterior ligaments generally exhibit high-stiffness in the loading region. Following this, the ankle model implemented with in situ ligament behavior was evaluated in response to multiple loading conditions and proved capable of predicting the bony kinematics accurately in comparison to the cadaveric response. Overall, the parametric ligament modeling demonstrated the feasibility of linking the gross structural behavior and the underlying bone and ligament mechanics that generate them. Determination of the in situ mechanical properties of ankle ligaments provides a better understanding of the nonlinear nature of the ankle joint. Applications of this knowledge include functional ankle joint mechanics and injury biomechanics.