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.     

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.     

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.     

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.     

Accuracy of a CT-based bone contour registration method to measure relative bone motions in the hindfoot

For measuring the in-vivo range of motion of the hindfoot, a CT-based bone contour registration method (CT-BCM) was developed to determine the three-dimensional position and orientation of bones. To validate this technique, we hypothesized that the range of motion in the hindfoot is equally, accurately measured by roentgen stereophotogrammetric analysis (RSA) as by the CT-BCM technique.Tantalum bone markers were placed in the distal tibia, talus and calcaneus of one cadaver specimen. With a fixed lower leg, the cadaveric foot was held in neutral and subsequently loaded in eight extreme positions. Immediately after acquiring a CT-scan with the foot in a position, RSA radiographs were made. Bone contour registration and RSA was performed. Helical axis parameters were calculated for talocrural and subtalar joint motion from neutral to extreme positions and between opposite extreme positions. Differences between CT-BCM and RSA were calculated.Compared with RSA, the CT-BCM data registered an overall root mean square difference (RMSd) of 0.21° for rotation about the helical axis, and 0.20 mm translation along the helical axis for the talocrural and subtalar joint and for all motions combined. The RMSd of the position and direction of the helical axes was 3.3 mm and 2.4°, respectively. The latter errors were larger with smaller helical rotations.The differences are similar to those reported for validated RSA and thus are not clinically relevant. Concluding, CT-BCM is an accurate and accessible alternative for studying joint motion, as it does not have the risk of infection and overlapping bone markers.     

Comparison of the hindfoot axes of a multi-segment foot model to the underlying bony anatomy

Musculoskeletal models used in gait analysis require coordinate systems to be identified for the body segments of interest. It is not obvious how hindfoot (or rearfoot) axes defined by skin-mounted markers relate to the anatomy of the underlying bones. The aim of this study was to compare the marker-based axes of the hindfoot in a multi-segment foot model to the orientations of the talus and calcaneus as characterized by their principal axes of inertia. Twenty adult females with no known foot deformities had radio-opaque markers placed on their feet and ankles at the foot model marker locations. CT images of the feet were acquired as the participants lay supine with their feet in a semi-weight bearing posture. The spatial coordinates of the markers were obtained from the images and used to define the foot model axes. Segmented masks of the tali and calcanei were used to create 3D bone models, from which the principal axes of the bones were obtained. The orientations of the principal axes were either within the range of typical values reported in the imaging literature or differed in ways that could be explained by variations in how the angles were defined. The model hindfoot axis orientations relative to the principal axes of the bones had little bias but were highly variable. Consideration of coronal plane hindfoot alignment as measured clinically and radiographically suggested that the model hindfoot coordinate system represents the posterior calcaneal tuberosity, rather than the calcaneus as a whole.     

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.     

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.     

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.     

Extended foot-ankle musculoskeletal models for application in movement analysis

Multibody simulations of human motion require representative models of the anatomical structures. A model that captures the complexity of the foot is still lacking. In the present work, two detailed 3D multibody foot-ankle models generated based on CT scans using a semi-automatic tool are described. The proposed models consists of five rigid segments (talus, calcaneus, midfoot, forefoot and toes), connected by five joints (ankle, subtalar, midtarsal, tarsometatarsal and metatarsophalangeal), one with 15DOF and the other with 8DOF. The calculated kinematics of both models were evaluated using gait trials and compared against literature, both presenting realistic results. An inverse dynamic analysis was performed for the 8DOF model, again presenting feasible dynamic results.     

In vivo determination of the Achilles tendon moment arm in three-dimensions

Two-dimensional methods have been applied to determine the Achilles tendon moment arm in previous studies, although the talocrural joint rotates in three-dimension. The purpose of this study was to develop a method for determining the Achilles tendon moment arm in three-dimensions (3DMA). A series of sagittal ankle images were obtained at ankle positions of -20°, -10° (dorsiflexed position), 0° (neutral position), +10°, +20°, and +30° (plantarflexed position). The talocrural joint axis was determined as the finite helical axis of the ankle joint over 20° of displacement, and the 3DMA was determined as the shortest distance from the talocrural joint axis to the line of action of the Achilles tendon force. The corresponding 2DMA was determined with the center of rotation method using the images captured on the sagittal plane passing through the mid-point of the medio-lateral width of the tibia. The 3DMA ranged from 35 to 41 mm across various ankle positions and was, on average, 11 mm smaller than 2DMA. The difference between the two measures was attributable primarily to the deviations of the talocrural joint axis from the anatomical medio-lateral direction. The deviations on the coronal plane (21.4±20.7°) and on the transverse planes (14.8±22.6°) accounted for the errors of 1.3 mm and 3.0 mm, respectively. In addition, selecting either a medially or laterally misaligned sagittal-plane image for determining the 2DMA gave rise to error by 3.5 mm. The remaining difference was accounted for by the random measurement error.     

The contraction-induced increase in Achilles tendon moment arm: A three-dimensional study

The present study aimed to re-examine the influence of the isometric plantarflexors contraction on the Achilles tendon moment arm (ATMA) and the factors influencing the ATMA in three-dimensions. A series of coronal magnetic resonance images of the right ankle were recorded at foot positions of 10° of dorsiflexion, neutral position, and 10° of plantarflexion for the rest condition and the plantarflexors contraction condition at 30% maximal voluntary effort. The shortest distance between the talocrural joint axis and the line of action of the Achilles tendon force projected to the orthogonal plane of the talocrural joint axis was determined as the ATMA. The ATMA determined in the contraction condition was significantly greater by 8 mm than that determined in the rest condition. The talocrural joint axis was displaced anteriorly by 3 mm and distally by 2 mm due to the muscle contraction. As the same time, the line of action of the Achilles tendon force was displaced posteriorly by 5 mm and medially by 2 mm. These linear displacements of the talocrural joint axis and the line of action of the Achilles tendon force accounted for the difference in the ATMAs between the two conditions by 35.9 and 62.4%, respectively. These angular displacements accounted for the total of 0.4% increase in the ATMA. These results confirm the previous findings reported in two-dimensional studies and found that the linear displacement of the line of action of the Achilles tendon force is the primary source of the contraction-induced increase in the ATMA.     

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.     

Feasibility of using MRIs to create subject-specific parallel-mechanism joint models

Musculoskeletal models typically use generic 2D models for the tibiofemoral (TFJ) and patellofemoral (PFJ) joints, with a hinge talocrural joint (TCJ), which are scaled to each subject׳s bone dimensions. Alternatively joints’ measured kinematics in cadavers are well-predicted using 3D cadaver-specific models. These employ mechanisms constrained by the articulations of geometric objects fitted to the joint׳s surfaces.In this study, we developed TFJ, PFJ and TCJ mechanism-based models off MRIs for fourteen participants and compared the estimated kinematics with those from published studies modified to be consistent with mechanisms models and subject-specific anatomical landmarks. The models’ parameters were estimated by fitting spheres to segmented articular cartilage surfaces, while ligament attachment points were selected from their bony attachment regions.Each participant׳s kinematics were estimated by ensuring no length changes in ligaments and constant distances between spheres’ centres. Two parameters’ optimizations were performed; both avoid singularities and one best matches the kinematic patterns off published studies. Sensitivity analysis determined which parameters the models were sensitive to.With both optimization methods, kinematics did not present singularities but correlation values were higher, exceeding 0.6, when matching the published studies. However, ranges of motion (ROM) were different between estimated and published studies. Across participants, models presented large parameter variation. Small variations were found between estimated- and optimized-parameters, and in the estimated-rotations and translations’ means and ROM. Model results were sensitive to changes in distal tibia, talus and patella spheres’ centres. These models can be implemented in subject-specific rigid-body musculoskeletal models to estimate joint moments and loads.     

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.     

A new non-orthogonal decomposition method to determine effective torques for three-dimensional joint rotation

This paper describes a new non-orthogonal decomposition method to determine effective torques for three-dimensional (3D) joint rotation. A rotation about a joint coordinate axis (e.g. shoulder internal/external rotation) cannot be explained only by the torque about the joint coordinate axis because the joint coordinate axes usually deviate from the principal axes of inertia of the entire kinematic chain distal to the joint. Instead of decomposing torques into three orthogonal joint coordinate axes, our new method decomposes torques into three "non-orthogonal effective axes" that are determined in such a way that a torque about each effective axis produces a joint rotation only about one of the joint coordinate axes. To demonstrate the validity of this new method, a simple internal/external rotation of the upper arm with the elbow flexed at 90 degrees was analyzed by both orthogonal and non-orthogonal decomposition methods. The results showed that only the non-orthogonal decomposition method could explain the cause-effect mechanism whereby three angular accelerations at the shoulder joint are produced by the gravity torque, resultant joint torque, and interaction torque. The proposed method would be helpful for biomechanics and motor control researchers to investigate the manner in which the central nervous system coordinates the gravity torque, resultant joint torque, and interaction torque to control 3D joint rotations.     

Analysis of the human and ape foot during bipedal standing with implications for the evolution of the foot

The ratio of the power arm (the distance from the heel to the talocrural joint) to the load arm (that from the talocrural joint to the distal head of the metatarsals), or RPL, differs markedly between the human and ape foot. The arches are relatively higher in the human foot in comparison with those in apes. This study evaluates the effect of these two differences on biomechanical effectiveness during bipedal standing, estimating the forces acting across the talocrural and tarsometatarsal joints, and attempts to identify which type of foot is optimal for bipedal standing. A simple model of the foot musculoskeletal system was built to represent the geometric and force relationships in the foot during bipedal standing, and measurements for a variety of human and ape feet applied. The results show that: (1) an RPL of around 40% (as is the case in the human foot) minimizes required muscle force at the talocrural joint; (2) the presence of an high arch in the human foot reduces forces in the plantar musculature and aponeurosis; and (3) the human foot has a lower total of force in joints and muscles than do the ape feet. These results indicate that the proportions of the human foot, and the height of the medial arch are indeed better optimized for bipedal standing than those of apes, further suggesting that their current state is to some extent the product of positive selection for enhanced bipedal standing during the evolution of the foot.     

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.