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.     

Influence of kinematic analysis methods on detecting ankle and subtalar joint instability

Patients with subtalar joint instability may be misdiagnosed with ankle instability, which may lead to chronic instability at the subtalar joint. Therefore, it is important to understand the difference in kinematics after ligament sectioning and differentiate the changes in kinematics between ankle and subtalar instability. Three methods may be used to determine the joint kinematics; the Euler angles, the Joint Coordinate System (JCS) and the helical axis (HA). The purpose of this study was to investigate the influence of using either method to detect subtalar and ankle joints instability. 3D kinematics at the ankle and subtalar joint were analyzed on 8 cadaveric specimens while the foot was intact and after sequentially sectioning the anterior talofibular ligament (ATFL), the calcaneofibular ligament (CFL), the cervical ligament and the interosseous talocalcaneal ligament (ITCL). Comparison in kinematics calculated from sensor and anatomical landmarks was conducted as well as the influence of Euler angles and JCS rotation sequence (between ISB recommendation and previous research) on the subtalar joint. All data showed a significant increase in inversion when the ITCL was sectioned. There were differences in the data calculated using sensors coordinate systems vs. anatomic coordinate systems. Anatomic coordinate systems were recommended for these calculations. The Euler angle and JCS gave similar results. Differences in Euler angles and JCS sequence lead to the same conclusion in detecting instability at the ankle and subtalar joint. As expected, the HA detected instability in plantarflexion at the ankle joint and in inversion at the subtalar joint.     

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.     

Use of dual Euler angles to quantify the three-dimensional joint motion and its application to the ankle joint complex

This paper presents a modified Euler angles method, dual Euler angles approach, to describe general spatial human joint motions. In dual Euler angles approach, the three-dimensional joint motion is considered as three successive screw motions with respect to the axes of the moving segment coordinate system; accordingly, the screw motion displacements are represented by dual Euler angles. The algorithm for calculating dual Euler angles from coordinates of markers on the moving segment is also provided in this study. As an example, the proposed method is applied to describe motions of ankle joint complex during dorsiflexion–plantarflexion. A Flock of Birds electromagnetic tracking device (FOB) was used to measure joint motion in vivo. Preliminary accuracy tests on a gimbal structure demonstrate that the mean errors of dual Euler angles evaluated by using source data from FOB are less than 1° for rotations and 1 mm for translations, respectively. Based on the pilot study, FOB is feasible for quantifying human joint motions using dual Euler angles approach.     

Subject-specific axes of the ankle joint complex

The aim of this study was to use a two-axis ankle joint model and an optimisation process (van den Bogert et al., 1994) to calculate and compare the talocrural and subtalar hinge axes for non-weight-bearing ankle motion, weight-bearing ankle motion, and walking in normal, healthy adult subjects and to see which of the first two sets of axes better fit the walking data. Motion data for the foot and shank were collected on eight subjects whilst they performed the activities mentioned. After choosing the best marker sets for motion tracking, a two-hinge ankle joint model was fit to the motion data. Ankle joint ranges of motion were also calculated. It was found that the model fit the experimental data well, with non-weight-bearing motion achieving the best fit. Despite this, the calculated axis orientations were highly variable both between motion types and between subjects. No significant difference between the fit of the non-weight-bearing and weight-bearing models to the walking data was found, which implies that either set of functional axes is adequate for modeling walking; however, the subtalar deviation angle was significantly closer for the weight-bearing activity and walking than for the non-weight-bearing activity and walking, which suggests that it is marginally better to use the weight-bearing functional motions. The results lead to questions about the appropriateness of the two-hinge ankle model for use in applications in which the behaviour of the individual joints of the ankle complex, rather than simply the relative motion of the leg and foot, is important.     

In-vivo range of motion of the subtalar joint using computed tomography

Understanding in vivo subtalar joint kinematics is important for evaluation of subtalar joint instability, the design of a subtalar prosthesis and for analysing surgical procedures of the ankle and hindfoot. No accurate data are available on the normal range of subtalar joint motion. The purpose of this study was to introduce a method that enables the quantification of the extremes of the range of motion of the subtalar joint in a loaded state using multidetector computed tomography (CT) imaging. In 20 subjects, an external load was applied to a footplate and forced the otherwise unconstrained foot in eight extreme positions. These extreme positions were foot dorsiflexion, plantarflexion, eversion, inversion and four extreme positions in between the before mentioned positions. CT images were acquired in a neutral foot position and each extreme position separately. After bone segmentation and contour matching of the CT data sets, the helical axes were determined for the motion of the calcaneus relative to the talus between four pairs of opposite extreme foot positions. The helical axis was represented in a coordinate system based on the geometric principal axes of the subjects’ talus. The greatest relative motion between the calcaneus and the talus was calculated for foot motion from extreme eversion to extreme inversion (mean rotation about the helical axis of 37.3±5.9°, mean translation of 2.3±1.1 mm). A consistent pattern of range of subtalar joint motion was found for motion of the foot with a considerable eversion and inversion component.     

Determining dual Euler angles of the ankle complex in vivo using "flock of birds" electromagnetic tracking device

The dual Euler angles method has been proposed as an alternative approach to describe the general spatial human joint motion. In this study, the dual Euler angles method was applied to study the three-dimensional motion of the ankle complex. The methodology for obtaining dual Euler angles of the ankle complex was developed by using a "Flock of Birds" electromagnetic tracking device. The repeatability of the methodology was studied based on the intertester and intratester variability analysis. Finally kinematic coupling characteristics of the ankle complex during dorsiflexion-plantarflexion, eversion-inversion, and abduction-adduction were analyzed according to the parameters of the dual Euler angles.     

Kinematics of the human ankle complex in passive flexion; a single degree of freedom system

The restoration of original range and pattern of motion is the primary goal of joint replacement and ligament reconstruction. The objective of the present work is to investigate whether or not a preferred path of joint motion at the intact human ankle complex is exhibited during passive flexion. A rig was built to move the ankle complex through its range of flexion while applying only the minimum necessary load to drive ankle flexion. Joint motion was constrained only by the articular surfaces and the ligaments. The movements of the calcaneus, talus and fibula relative to the stationary tibia in seven cadaveric specimens were tracked with a stereophotogrammetric system. It was shown that the calcaneus follows a unique path of unresisted coupled motion relative to the tibia and that most of the motion occurred at the ankle, with little motion at the subtalar level. The calcaneofibular and the tibiocalcaneal ligaments showed near-isometric pattern of rotations. All specimens showed motion of the axis of rotation relative to the bones. Deviations from the unique path due to the application of load involved mostly subtalar motion and were resisted. The ankle complex exhibits one degree of unresisted freedom, the ankle behaving as a single degree of freedom mechanism and the subtalar as a flexible structure. We deduced that the calcaneofibular and tibiocalcaneal ligaments together with the articular surfaces guide ankle passive motion, other ligaments limit but do not guide motion.     

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.     

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.     

Mobility of the subtalar joint in the intact ankle complex

A previous study by these authors showed that the calcaneus follows a unique path of unresisted coupled motion relative to the tibia during passive flexion and that most of this motion occurred at the ankle level. Subtalar motion in the intact ankle complex was observed only when perturbations from this path were induced by the application of force to the calcaneus. Relative motion of the bones of the ankle complex was tracked by stereophotogrammetry in seven specimens. Anatomical landmarks, reference frames and joint angles were defined by standard techniques. Sequential moments were applied to the calcaneus about the long axis of the tibia. Measured movements at subtalar level demonstrated plantarflexion coupled to supination and internal rotation (inversion) and dorsiflexion coupled to pronation and external rotation (eversion). These movements were resisted and were fully recovered when the external load was removed. Subtalar motion diminished as the ankle approached maximal dorsi- and plantarflexion. Two clearly distinguished mean axes of rotation were observed for inversion and eversion runs. The axes of inversion and eversion of the subtalar complex changed orientation along a preferred and repeatable path. The subtalar joint complex occupied only a single stable position in the unloaded state and with no range of unresisted motion. It is inferred that mobility was possible only by the stretching and lengthening of the ligaments and the indentation of the articular surfaces, requiring the application of loads. The subtalar joint complex behaves like a flexible structure.     

In vitro assessment of a motion-based optimization method for locating the talocrural and subtalar joint axes

The locations of the joint axes of the ankle complex vary considerably between subjects, yet no noninvasive method with demonstrated accuracy exists for locating these axes. The moments of muscle and ground reaction forces about the joint axes are dependent on axis locations, making knowledge of these locations critical to accurate musculoskeletal modeling of the foot and ankle. The accuracy of a computational optimization method that fits a two-revolute model to measured motion was assessed using computer-generated data, a two-revolute mechanical linkage, and three lower-leg cadaver specimens. Motions were applied to cadaver specimens under axial load while bone-mounted markers attached to the tibia, talus, and calcaneus were tracked using a video-based motion analysis system. Estimates of the talocrural and subtalar axis locations were computed from motions of the calcaneus relative to the tibia using the optimization method. These axes were compared to mean helical axes computed directly from tibia, talus, and calcaneus motions. The optimization method performed well when the motions were computer-generated or measured in the mechanical linkage, with angular differences between optimization and mean helical axes ranging from 1 deg to 5 deg. In the cadaver specimens, however, these differences exceeded 20 deg. Optimization methods that locate the anatomical joint axes of the ankle complex by fitting two revolute joints to measured tibia-calcaneus motions may be limited because of problems arising from non-revolute behavior.     

A multi-segment kinematic model of the foot with a novel definition of forefoot motion for use in clinical gait analysis during walking

A multi-segment kinematic model of the foot was developed for use in a gait analysis laboratory. The foot was divided into hindfoot, talus, midfoot and medial and lateral forefoot segments. Six functional joints were defined: ankle and subtalar joints, frontal and transverse plane motions of the hindfoot relative to midfoot, supination/pronation twist of the forefoot relative to midfoot and medial longitudinal arch height-to-length ratio. Twelve asymptomatic subjects were tested during barefoot walking with a six-camera optical stereometric system and auto-reflective markers organized in triads. Repeatability of the joint motions was tested using coefficients of multiple correlation. Ankle and subtalar joint motions and twisting of the forefoot were most repeatable. Hindfoot motions were least repeatable both within-subjects and between-subjects. Hindfoot and forefoot pronation in the frontal plane was found to coincide with dropping of the medial longitudinal arch between early to mid-stance, followed by supination and rising of the arch in late stance and swing phase. This multi-segment foot model addresses an unfortunate shortcoming in current gait analysis practice-the inability to measure motion within the foot. Such measurements are crucial if gait analysis is to remain relevant in the orthopaedic and rehabilitative treatment of the foot and ankle.     

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.     

In vivo tests of an improved method for functional location of the subtalar joint axis

The subtalar joint is important in frontal plane movement and posture of the hindfoot. Abnormal subtalar joint moments caused by muscle forces and the ground reaction force acting on the foot are thought to play a role in various foot deformities. Calculating joint moments typically requires knowledge of the location of the joint axis; however, location of the subtalar axis from measured movement is difficult because the talus cannot be tracked using skin-mounted markers. The accuracy of a novel technique for locating the subtalar axis was assessed in vivo using magnetic resonance imaging. The method was also tested with skin-mounted markers and video motion analysis. The technique involves applying forces to the foot that cause pure subtalar joint motion (with negligible talocrural joint motion), and then using helical axis decomposition of the resulting tibiocalcaneal motion. The resulting subtalar axis estimates differed by 6° on average from the true best-fit subtalar axes in the MRI tests. Motion was found to have been applied primarily about the subtalar joint with an average of only 3° of talocrural joint motion. The proposed method provides a potential means for obtaining subject-specific subtalar axis estimates which can then be used in inverse dynamic analyses and subject-specific musculoskeletal models.     

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.     

Online estimation algorithm for a biaxial ankle kinematic model with configuration dependent joint axes

The kinematics of the human ankle is commonly modeled as a biaxial hinge joint model. However, significant variations in axis orientations have been found between different individuals and also between different foot configurations. For ankle rehabilitation robots, information regarding the ankle kinematic parameters can be used to estimate the ankle and subtalar joint displacements. This can in turn be used as auxiliary variables in adaptive control schemes to allow modification of the robot stiffness and damping parameters to reduce the forces applied at stiffer foot configurations. Due to the large variations observed in the ankle kinematic parameters, an online identification algorithm is required to provide estimates of the model parameters. An online parameter estimation routine based on the recursive least-squares (RLS) algorithm was therefore developed in this research. An extension of the conventional biaxial ankle kinematic model, which allows variation in axis orientations with different foot configurations had also been developed and utilized in the estimation algorithm. Simulation results showed that use of the extended model in the online algorithm is effective in capturing the foot orientation of a biaxial ankle model with variable joint axis orientations. Experimental results had also shown that a modified RLS algorithm that penalizes a deviation of model parameters from their nominal values can be used to obtain more realistic parameter estimates while maintaining a level of estimation accuracy comparable to that of the conventional RLS routine.     

An image-based kinematic model of the tibiotalar and subtalar joints and its application to gait analysis in children with Juvenile Idiopathic Arthritis

In vivo estimates of tibiotalar and the subtalar joint kinematics can unveil unique information about gait biomechanics, especially in the presence of musculoskeletal disorders affecting the foot and ankle complex. Previous literature investigated the ankle kinematics on ex vivo data sets, but little has been reported for natural walking, and even less for pathological and juvenile populations. This paper proposes an MRI-based morphological fitting methodology for the personalised definition of the tibiotalar and the subtalar joint axes during gait, and investigated its application to characterise the ankle kinematics in twenty patients affected by Juvenile Idiopathic Arthritis (JIA). The estimated joint axes were in line with in vivo and ex vivo literature data and joint kinematics variation subsequent to inter-operator variability was in the order of 1°. The model allowed to investigate, for the first time in patients with JIA, the functional response to joint impairment. The joint kinematics highlighted changes over time that were consistent with changes in the patient’s clinical pattern and notably varied from patient to patient. The heterogeneous and patient-specific nature of the effects of JIA was confirmed by the absence of a correlation between a semi-quantitative MRI-based impairment score and a variety of investigated joint kinematics indexes. In conclusion, this study showed the feasibility of using MRI and morphological fitting to identify the tibiotalar and subtalar joint axes in a non-invasive patient-specific manner. The proposed methodology represents an innovative and reliable approach to the analysis of the ankle joint kinematics in pathological juvenile populations.     

The impact of high-heeled shoes on ankle complex during walking in young women—In vivo kinematic study based on 3D to 2D registration technique

Objective: To explore the accurate in vivo kinematic changes in the ankle complex when wearing low- and high-heel shoes (LHS and HHS, respectively). Materials and methods: Twelve young women were tested unilaterally. Three-dimensional models of the tibia, talus, and calcaneus were first created based on CT scan results. The subjects walked at a self-controlled speed in barefoot, LHS (4 cm), and HHS (10 cm) conditions. A fluoroscopy system captured the lateral fluoroscopic images of the ankle complex. The images of seven key positions in the stance phase were selected, and 3D to 2D bone model registrations were performed to determine the joint positions. The mean of 6 degree of freedom (DOF) range of motions (ROM), joint positions, and angular displacements of the ankle complex during the gait were then obtained. Results: For the talocrural joint, the rotational ROMs of the subjects either in LHS or HHS condition displayed no significant difference from those in barefoot condition. For the subtalar joint, all the rotational ROMs in the HHS condition and the internal/external rotations in the LHS condition significantly decreased compared with those in the barefoot condition. The talocrural joint was positioned significantly more plantarflexed, inverted, internally rotated, and posteriorly seated in all seven poses in HHS condition, compared with those in barefoot condition. Conclusion: HHS mainly affected the rotational motion of the ankle complex during walking. The talocrural joint position was abnormal, and the subtalar joint ROM decreased during the gait in HHS condition. Only a few kinematic changes occurred in LHS condition relative to the barefoot condition.     

Frontal plane moments do not accurately reflect ankle dynamics during running

The ankle joint has typically been treated as a universal joint with moments calculated about orthogonal axes and the frontal plane moment generally used to represent the net muscle action about the subtalar joint. However, this joint acts about an oblique axis. The purpose of this study was to examine the differences between joint moments calculated about the orthogonal frontal plane axis and an estimated subtalar joint axis. Three-dimensional data were collected on 10 participants running at 3.6 m/s. Joint moments, power, and work were calculated about the orthogonal frontal plane axis of the foot and about an oblique axis representing the subtalar joint. Selected parameters were compared with a paired t-test (alpha = 0.05). The results indicated that the joint moments calculated about the two axes were characteristically different. A moment calculated about an orthogonal frontal plane axis of the foot resulted in a joint moment that was invertor in nature during the first half of stance, but evertor during the second half of stance. The subtalar joint axis moment, however, was invertor during most of the stance. These two patterns may result in qualitatively different interpretations of the muscular contributions at the ankle during the stance phase of running.