Rotational stiffness of football shoes influences talus motion during external rotation of the foot

Shoe-surface interface characteristics have been implicated in the high incidence of ankle injuries suffered by athletes. Yet, the differences in rotational stiffness among shoes may also influence injury risk. It was hypothesized that shoes with different rotational stiffness will generate different patterns of ankle ligament strain. Four football shoe designs were tested and compared in terms of rotational stiffness. Twelve (six pairs) male cadaveric lower extremity limbs were externally rotated 30 deg using two selected football shoe designs, i.e., a flexible shoe and a rigid shoe. Motion capture was performed to track the movement of the talus with a reflective marker array screwed into the bone. A computational ankle model was utilized to input talus motions for the estimation of ankle ligament strains. At 30 deg of rotation, the rigid shoe generated higher ankle joint torque at 46.2?±?9.3 Nm than the flexible shoe at 35.4?±?5.7 Nm. While talus rotation was greater in the rigid shoe (15.9?±?1.6 deg versus 12.1?±?1.0 deg), the flexible shoe generated more talus eversion (5.6?±?1.5 deg versus 1.2± 0.8 deg). While these talus motions resulted in the same level of anterior deltoid ligament strain (approxiamtely 5%) between shoes, there was a significant increase of anterior tibiofibular ligament strain (4.5± 0.4% versus 2.3?±?0.3%) for the flexible versus more rigid shoe design. The flexible shoe may provide less restraint to the subtalar and transverse tarsal joints, resulting in more eversion but less axial rotation of the talus during foot∕shoe rotation. The increase of strain in the anterior tibiofibular ligament may have been largely due to the increased level of talus eversion documented for the flexible shoe. There may be a direct correlation of ankle joint torque with axial talus rotation, and an inverse relationship between torque and talus eversion. The study may provide some insight into relationships between shoe design and ankle ligament strain patterns. In future studies, these data may be useful in characterizing shoe design parameters and balancing potential ankle injury risks with player performance.     

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 \(20{^{\circ }}\) of plantar flexion to \(15{^{\circ }}\) of dorsiflexion, and from \(15{^{\circ }}\) of inversion to \(15{^{\circ }}\) 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.     

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.     

A framework for parametric modeling of ankle ligaments to determine the in situ response under gross foot motion

Ligament sprains account for a majority of injuries to the foot and ankle complex, but ligament properties have not been understood well due to the difficulties in replicating the complex geometry, in situ stress state, and non-uniformity of the strain. For a full investigation of the injury mechanism, it is essential to build up a foot and ankle model validated at the level of bony kinematics and ligament properties. This study developed a framework to parameterize the ligament response for determining the in situ stress state and heterogeneous force–elongation characteristics using a finite element ankle model. Nine major ankle ligaments and the interosseous membrane were modeled as discrete elements corresponding functionally to the ligamentous microstructure of collagen fibers and having parameterized toe region and stiffness at the fiber level. The range of the design variables in the ligament model was determined from existing experimental data. Sensitivity of the bony kinematics to each variable was investigated by design of experiment. The results highlighted the critical role of the length of the toe region of the ligamentous fibers on the bony kinematics with the cumulative influence of more than 95%, while the fiber stiffness was statistically insignificant with an influence of less than 1% under the given variable range and loading conditions. With the flexibility of variable adjustment and high computational efficiency, the presented ankle model was generic in nature so as to maximize its applicability to capture the individual ligament behaviors in future studies.     

The three-dimensional kinematics and flexibility characteristics of the human ankle and subtalar joints--Part I: Kinematics

The in-vitro, three dimensional kinematic characteristics of the human ankle and subtalar joint were investigated in this study. The main goals of this investigation were: 1) To determine the range of motion of the foot-shank complex and the associated range of motion of the ankle and subtalar joints; 2) To determine the kinematic coupling characteristics of the foot-shank complex, and 3) To identify the relationship between movements at the ankle and subtalar joints and the resulting motion produced between the foot and the shank. The tests were conducted on fifteen fresh amputated lower limbs and consisted of incrementally displacing the foot with respect to the shank while the motion of the articulating bones was measured through a three dimensional position data acquisition system. The kinematic analysis was based on the helical axis parameters describing the incremental displacements between any two of the three articulating bones and on a joint coordinate system used to describe the relative position between the bones. From the results of this investigation it was concluded that: 1) The range of motion of the foot-shank complex in any direction (dorsiflexion/plantarflexion, inversion/eversion and internal rotation/external rotation) is larger than that of either the ankle joint or the subtalar joint.; 2) Large kinematic coupling values are present at the foot-shank complex in inversion/eversion and in internal rotation/external rotation. However, only a slight amount of coupling was observed to occur in dorsiflexion/plantarflexion.; 3) Neither the ankle joint nor the subtalar joint are acting as ideal hinge joints with a fixed axis of rotation.; 4) Motion of the foot-shank complex in any direction is the result of rotations at both the ankle and the subtalar joints. However, the contribution of the ankle joint to dorsiflexion/plantarflexion of the foot-shank complex is larger than that of the subtalar joint and the contribution of the subtalar joint to inversion/eversion is larger than that of the ankle joint.; 5) The ankle and the subtalar joints have an approximately equal contribution to internal rotation/external rotation movements of the foot-shank complex.     

Effect of simulated joint instability and bracing on ankle and subtalar joint flexibility

It is clinically challenging to distinguish between ankle and subtalar joints instability in vivo. Understanding the changes in load-displacement at the ankle and subtalar joints after ligament injuries may detect specific changes in joint characteristics that cannot be detected by investigating changes in range of motion alone. The effect of restricting joints end range of motion with ankle braces was already established, but little is known about the effect of an ankle brace on the flexibility of the injured ankle and subtalar joints. Therefore, the purposes of this study were to (1) understand how flexibility is affected at the ankle and subtalar joints after sectioning lateral and intrinsic ligaments during combined sagittal foot position and inversion and during internal rotation and (2) investigate the effect of a semi-rigid ankle brace on the ankle and subtalar joint flexibility. Kinematics and kinetics were collected from nine cadaver feet during inversion through the range of ankle flexion and during internal rotation. Motion was applied with and without a brace on an intact foot and after sequentially sectioning the calcaneofibular ligament (CFL) and the intrinsic ligaments. Segmental flexibility was defined as the slope of the angle-moment curve for each 1 Nm interval. Early flexibility significantly increased at the ankle and subtalar joint after CFL sectioning during inversion. The semi-rigid ankle brace significantly decreased early flexibility at the subtalar joint during inversion and internal rotation for all ligament conditions and at the ankle joint after all ligaments were cut.     

A real-time computational model for estimating kinematics of ankle ligaments

Background: An accurate assessment of ankle ligament kinematics is crucial in understanding the injury mechanisms and can help to improve the treatment of an injured ankle, especially when used in conjunction with robot-assisted therapy. A number of computational models have been developed and validated for assessing the kinematics of ankle ligaments. However, few of them can do real-time assessment to allow for an input into robotic rehabilitation programs. Method: An ankle computational model was proposed and validated to quantify the kinematics of ankle ligaments as the foot moves in real-time. This model consists of three bone segments with three rotational degrees of freedom (DOFs) and 12 ankle ligaments. This model uses inputs for three position variables that can be measured from sensors in many ankle robotic devices that detect postures within the foot–ankle environment and outputs the kinematics of ankle ligaments. Validation of this model in terms of ligament length and strain was conducted by comparing it with published data on cadaver anatomy and magnetic resonance imaging. Results: The model based on ligament lengths and strains is in concurrence with those from the published studies but is sensitive to ligament attachment positions. Conclusions: This ankle computational model has the potential to be used in robot-assisted therapy for real-time assessment of ligament kinematics. The results provide information regarding the quantification of kinematics associated with ankle ligaments related to the disability level and can be used for optimizing the robotic training trajectory.     

A novel device to apply controlled flexion and extension to the rat knee following anterior cruciate ligament reconstruction

We designed and validated a novel device for applying flexion-extension cycles to a rat knee in an in vivo model of anterior cruciate ligament reconstruction (ACL-R). Our device is intended to simulate rehabilitation motion and exercise post ACL-R to optimize physical rehabilitation treatments for the improved healing of tendon graft ligament reconstructions. The device was validated for repeatability of the knee kinematic motion by measuring the force versus angular rotation response from repeated trials using cadaver rats. The average maximum force required for rotating an ACL reconstructed rat knee through 100 degrees of flexion-extension was 0.4 N with 95% variability for all trials within ±0.1 N.     

Computational modeling to predict mechanical function of joints: application to the lower leg with simulation of two cadaver studies

Computational models of musculoskeletal joints and limbs can provide useful information about joint mechanics. Validated models can be used as predictive devices for understanding joint function and serve as clinical tools for predicting the outcome of surgical procedures. A new computational modeling approach was developed for simulating joint kinematics that are dictated by bone/joint anatomy, ligamentous constraints, and applied loading. Three-dimensional computational models of the lower leg were created to illustrate the application of this new approach. Model development began with generating three-dimensional surfaces of each bone from CT images and then importing into the three-dimensional solid modeling software SOLIDWORKS and motion simulation package COSMOSMOTION. Through SOLIDWORKS and COSMOSMOTION, each bone surface file was filled to create a solid object and positioned necessary components added, and simulations executed. Three-dimensional contacts were added to inhibit intersection of the bones during motion. Ligaments were represented as linear springs. Model predictions were then validated by comparison to two different cadaver studies, syndesmotic injury and repair and ankle inversion following ligament transection. The syndesmotic injury model was able to predict tibial rotation, fibular rotation, and anterior/posterior displacement. In the inversion simulation, calcaneofibular ligament extension and angles of inversion compared well. Some experimental data proved harder to simulate accurately, due to certain software limitations and lack of complete experimental data. Other parameters that could not be easily obtained experimentally can be predicted and analyzed by the computational simulations. In the syndesmotic injury study, the force generated in the tibionavicular and calcaneofibular ligaments reduced with the insertion of the staple, indicating how this repair technique changes joint function. After transection of the calcaneofibular ligament in the inversion stability study, a major increase in force was seen in several of the ligaments on the lateral aspect of the foot and ankle, indicating the recruitment of other structures to permit function after injury. Overall, the computational models were able to predict joint kinematics of the lower leg with particular focus on the ankle complex. This same approach can be taken to create models of other limb segments such as the elbow and wrist. Additional parameters can be calculated in the models that are not easily obtained experimentally such as ligament forces, force transmission across joints, and three-dimensional movement of all bones. Muscle activation can be incorporated in the model through the action of applied forces within the software for future studies.     

A closed-loop cadaveric foot and ankle loading model

Investigations of human foot and ankle biomechanics rely chiefly on cadaver experiments. The application of proper force magnitudes to the cadaver foot and ankle is essential to obtain valid biomechanical data. Data for external ground reaction forces are readily available from human motion analysis. However, determining appropriate forces for extrinsic foot and ankle muscles is more problematic. A common approach is the estimation of forces from muscle physiological cross-sectional areas and electromyographic data. We have developed a novel approach for loading the Achilles and posterior tibialis tendons that does not prescribe predetermined muscle forces. For our loading model, these muscle forces are determined experimentally using independent plantarflexion and inversion angle feedback control. The independent (input) parameters -- calcaneus plantarflexion, calcaneus inversion, ground reaction forces, and peroneus forces -- are specified. The dependent (output) parameters -- Achilles force, posterior tibialis force, joint motion, and spring ligament strain -- are functions of the independent parameters and the kinematics of the foot and ankle. We have investigated the performance of our model for a single, clinically relevant event during the gait cycle. The instantaneous external forces and foot orientation determined from human subjects in a motion analysis laboratory were simulated in vitro using closed-loop feedback control. Compared to muscle force estimates based on physiological cross-sectional area data and EMG activity at 40% of the gait cycle, the posterior tibialis force and Achilles force required when using position feedback control were greater.     

Analysis of surface-to-surface distance mapping during three-dimensional motion at the ankle and subtalar joints

Joint surface interaction and ligament constraints determine the kinematic characteristics of the ankle and subtalar joints. Joint surface interaction is characterized by joint contact mechanics and by relative joint surface position potentially characterized by distance mapping. While ankle contact mechanics was investigated, limited information is available on joint distance mapping and its changes during motion. The purpose of this study was to use image-based distance mapping to quantify this interaction at the ankle and subtalar joints during tri-planar rotations of the ankle complex. Five cadaveric legs were scanned using Computed Tomography and the images were processed to produce 3D bone models of the tibia, fibula, talus and calcaneus. Each leg was tested on a special linkage through which the ankle complex was loaded in dorsiflexion/plantarflexion, inversion/eversion, and internal/external rotation and the resulting bone movements were recorded. Fiduciary bone markers data and 3D bone models were combined to generate color-coded distance maps for the ankle and subtalar joints. The results were processed focusing on the changes in surface-to-surface distance maps between the extremes of the range of motion and neutral. The results provided detailed insight into the three-dimensional highly coupled nature of these joints showing significant and unique changes in distance mapping from neutral to extremes of the range of motion. The non-invasive nature of the image-based distance mapping technique could result, after proper modifications, in an effective diagnostic and clinical evaluation technique for application such as ligament injuries and quantifying the effect of arthrodesis or total ankle replacement surgery.     

In vivo determination of the anatomical axes of the ankle joint complex: An optimization approach

This study investigates the feasibility of a subject-specific three-dimensional model of the ankle joint complex for kinematic and dynamic analysis of movement. The ankle joint complex was modelled as a three-segment system, connected by two ideal highe joints: the talocrural and the subtalar joint. A mathematical formulation was developed to express the three-dimensional translation and rotation between the foot and shank segments as a function of the two joint angles, and 12 model parameters describing the locations of the joint axes. An optimization method was used to fit the model parameters to three-dimensional kinematic data of foot and shank markers, obtained during test movements throughout the entire physiological range of motion of the ankle joint. The movement of the talus segment, which cannot be measured non-invasively, is not necessary for the analysis.

This optimization method was used to determine the position and orientation of the joint axes in 14 normal subjects. After optimization, the discrepancy between the best fitting model and actual marker kinematics was between 1 and 3 mm for all subjects. The predicted inclination of the subtalar joint axis from the horizontal plane was 37.4±2.7°, and the medial deviation was 18.0±16.2°. The lateral side of the talucrural axis was directed slightly posteriorly (6.8±8.1°), and inclined downward by 7.0±5.4°. These results are similar to previously reported typical results from anatomical, in vitro, studies. Reproducibility was evaluated by repeated testing of one subject, which resulted in variations of about one-fifth of the standard deviation within the group, the inclination of the subtalar joint axis was significantly correlated to the arch height and a radiographic ‘tarsal index’. It is concluded that this optimization method provides the opportunity to incorporate inter-individual anatomical differences into kinematic and dynamic analysis of the ankle joint complex. This allows a more functional interpretation of kinematic data, and more realistic estimates of internal forces.     

The influence of foot positioning on ankle sprains

The goal of this study was to examine the influence of changes in foot positioning at touch-down on ankle sprain occurrence. Muscle model driven computer simulations of 10 subjects performing the landing phase of a side-shuffle movement were performed. The relative subtalar joint and talocural joint angles at touchdown were varied, and each subject-specific simulation was exposed to a set of perturbed floor conditions. The touchdown subtalar joint angle was not found to have a considerable influence on sprain occurrence, while increased touchdown plantar flexion caused increased ankle sprain occurrences. Increased touchdown plantar flexion may be the mechanism which causes ankles with a history of ankle sprains to have an increased susceptibility to subsequent sprains. This finding may also reveal a mechanism by which taping of a sprained ankle or the application of an ankle brace leads to decreased ankle sprain susceptibility.     

Mechanics of the ankle and subtalar joints revealed through a 3D quasi-static stress MRI technique

A technique to study the three-dimensional (3D) mechanical characteristics of the ankle and of the subtalar joints in vivo and in vitro is described. The technique uses an MR scanner compatible 3D positioning and loading linkage to load the hindfoot with precise loads while the foot is being scanned. 3D image processing algorithms are used to derive from the acquired MR images bone morphology, hindfoot architecture, and joint kinematics. The technique was employed to study these properties both in vitro and in vivo. The ankle and subtler joint motion and the changes in architecture produced in response to an inversion load and an anterior drawer load were evaluated. The technique was shown to provide reliable measures of bone morphology. The left-to-right variations in bone morphology were less than 5%. The left-to-right variations in unloaded hindfoot architecture parameters were less than 10%, and these properties were only slightly affected by inversion and anterior drawer loads. Inversion and anterior drawer loads produced motion both at the ankle and at the subtalar joint. In addition, high degree of coupling, primarily of internal rotation with inversion, was observed both at the ankle and at the subtalar joint. The in vitro motion produced in response to inversion and anterior drawer load was greater than the in vivo motion. Finally, external motion, measured directly across the ankle complex, produced in response to load was much greater than the bone movements measured through the 3D stress MRI technique indicating the significant effect of soft tissue and skin interference.     

Contributions of proximal and distal moments to axial tibial rotation during walking and running

The purpose of this study was to determine the cause and effect relationship between tibial internal rotation and pronation of the foot during walking and heel-toe running. This would allow predictions of orthotic effectiveness in reducing knee pain related to excessive internal tibial rotation. Kinematic and force plate data were collected from twenty subjects performing ten running and ten walking trials across a force plate. Using a least-squares algorithm, attitude matrices for each segment in each frame were obtained and the angular velocity vector of the tibia was calculated. The intersegmental moment at the ankle was calculated from ground reaction force and kinematic data, and the power flow from foot to tibia associated with axial tibial rotation was calculated. In walking, all subjects exhibited a clear power flow from tibia to foot during most of the stance phase, indicating that the foot was following the body. This suggests that the use of foot orthoses to reduce knee pain associated with tibial rotation during walking will not be successful. During running, power flow was also mainly proximal to distal, but there were brief periods of opposite power flow. There was more variability between subjects during running, with five subjects having large distal to proximal power flow peaks. These observations may explain and support previous work that has found variable clinical effects of orthoses between patients.     

Concurrent musculoskeletal dynamics and finite element analysis predicts altered gait patterns to reduce foot tissue loading

Current computational methods for simulating locomotion have primarily used muscle-driven multibody dynamics, in which neuromuscular control is optimized. Such simulations generally represent joints and soft tissue as simple kinematic or elastic elements for computational efficiency. These assumptions limit application in studies such as ligament injury or osteoarthritis, where local tissue loading must be predicted. Conversely, tissue can be simulated using the finite element method with assumed or measured boundary conditions, but this does not represent the effects of whole body dynamics and neuromuscular control. Coupling the two domains would overcome these limitations and allow prediction of movement strategies guided by tissue stresses. Here we demonstrate this concept in a gait simulation where a musculoskeletal model is coupled to a finite element representation of the foot. Predictive simulations incorporated peak plantar tissue deformation into the objective of the movement optimization, as well as terms to track normative gait data and minimize fatigue. Two optimizations were performed, first without the strain minimization term and second with the term. Convergence to realistic gait patterns was achieved with the second optimization realizing a 44% reduction in peak tissue strain energy density. The study demonstrated that it is possible to alter computationally predicted neuromuscular control to minimize tissue strain while including desired kinematic and muscular behavior. Future work should include experimental validation before application of the methodology to patient care.     

Locomotion pattern and foot pressure adjustments during gentle turns in healthy subjects

People suffering from locomotor impairment find turning manoeuvres more challenging than straight-ahead walking. Turning manoeuvres are estimated to comprise a substantial proportion of steps taken daily, yet research has predominantly focused on straight-line walking, meaning that the basic kinetic, kinematic and foot pressure adaptations required for turning are not as well understood. We investigated how healthy subjects adapt their locomotion patterns to accommodate walking along a gently curved trajectory (radius 2.75 m). Twenty healthy adult participants performed walking tasks at self-selected speeds along straight and curved pathways. For the first time for this mode of turning, plantar pressures were recorded using insole foot pressure sensors while participants’ movements were simultaneously tracked using marker-based 3D motion capture. During the steady-state strides at the apex of the turn, the mean operating point of the inside ankle shifted by 1 degree towards dorsiflexion and that for the outside ankle shifted towards plantarflexion. The largest change in relative joint angle range was an increase in hip rotation in the inside leg (>60%). In addition, the inside foot was subject to a prolonged stance phase and a 10% increase in vertical force in the posteromedial section of the foot compared to straight-line walking. Most of the mechanical change required was therefore generated by the inside leg with hip rotation being a major driver of the gentle turn. This study provides new insight into healthy gait during gentle turns and may help us to understand the mechanics behind some forms of impairment.     

Biomechanical but not timed performance asymmetries persist between limbs 9 months after ACL reconstruction during planned and unplanned change of direction

Whilst anterior cruciate ligament injury commonly occurs during change of direction (CoD) tasks, there is little research on how athletes execute CoD after anterior cruciate ligament reconstruction (ACLR). The aims of this study were to determine between-limb and between-test differences in performance (time) and joint kinematics and kinetics during planned and unplanned CoD. One hundred and fifty-six male subjects carried out 90° maximal effort, planned and unplanned CoD tests in a 3D motion capture laboratory 9 months after ACLR. Statistical parametric mapping (2 × 2 ANOVA; limb × test) was used to identify differences in CoD time and biomechanical measures between limbs and between tests. There was no interaction effect but a main effect for limb and task. There was no between-limb difference in the time to complete both CoD tests. Between-limb differences were found for internal knee valgus moment, knee internal rotation and flexion angle, knee extension and external rotation moment and ankle external rotation moment with lower values on the ACLR side (effect size 0.72–0.5). Between test differences were found with less contralateral pelvis rotation, distance from centre of mass to the ankle in frontal plane, posterior ground reaction force and greater hip abduction during the unplanned CoD (effect size 0.75–0.5). Findings demonstrated that kinematic and kinetic differences between limbs are evident during both CoD tests 9 months after surgery, despite no statistical differences in performance time. Biomechanical differences between tests were found in variables, which have previously been associated with ACL injury mechanism during unplanned CoD.