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.     

A geometric model of the human ankle joint.

A two-dimensional four-bar linkage model of the ankle joint is formulated to describe dorsi/plantarflexion in unloaded conditions as observed in passive tests on ankle complex specimens. The experiments demonstrated that the human ankle joint complex behaves as a single-degree-of-freedom system during passive motion, with a moving axis of rotation. The bulk of the movement occurred at the level of the ankle. Fibres within the calcaneofibular and tibiocalcaneal ligaments remained approximately isometric. The experiments showed that passive kinematics of the ankle complex is governed only by the articular surfaces and the ligaments. It was deduced that the ankle is a single-degree-of-freedom mechanism where mobility is allowed by the sliding of the articular surfaces upon each other and the isometric rotation of two ligaments about their origins and insertions, without tissue deformation. The linkage model is formed by the tibia/fibula and talus/calcaneus bone segments and by the calcaneofibular and tibiocalcaneal ligament segments. The model predicts the path of calcaneus motion, ligament orientations, instantaneous axis of rotation, and conjugate talus surface profile as observed in the experiments. Many features of ankle kinematics such as rolling and multiaxial rotation are elucidated. The geometrical model is a necessary preliminary step to the study of ankle joint stability in response to applied loads and can be used to predict the effects of changes to the original geometry of the intact joint. Careful reconstruction of the original geometry of the ligaments is necessary after injury or during total ankle replacement.     

Ligament fibre recruitment at the human ankle joint complex in passive flexion

Knowledge of ligament fibre recruitment at the human ankle joint complex is a fundamental prerequisite for analysing mobility and stability. Previous experimental and modelling studies have shown that ankle motion must be guided by fibres within the calcaneofibular and tibiocalcaneal ligaments, which remain approximately isometric during passive flexion. The purpose of this study was to identify these fibres.

Three below-knee amputated specimens were analysed during passive flexion with combined radiostereometry for bone pose estimation and 3D digitisation for ligament attachment area identification. A procedure based on singular value decomposition enabled matching bone pose with digitised data and therefore reconstructing position in space of ligament attachment areas in each joint position. Eleven ordered fibres, connecting corresponding points on origin and insertion curves, were modelled for each of the following ligaments: posterior talofibular, calcaneofibular, anterior talofibular, posterior tibiotalar, tibiocalcaneal, and anterior tibiotalar.

The measured changes in length for the ligament fibres revealed patterns of tightening and slackening. The most anterior fibre of the calcaneofibular and the medio-anterior fibre of the tibiocalcaneal ligament exhibited the most isometric behaviour, as well as the most posterior fibre of the anterior talofibular ligament. Fibres within the calcaneofibular ligament remain parallel in the transverse plane, while those within the tibiocalcaneal ligament become almost parallel in joint neutral position. For both these ligaments, fibres maintain their relative inclination in the sagittal plane throughout the passive flexion range.

The observed significant change in both shape and orientation of the ankle ligaments suggest that this knowledge is fundamental for future mechanical analysis of their response to external forces.     

A new one-DOF fully parallel mechanism for modelling passive motion at the human tibiotalar joint

Knowledge on how ligaments and articular surfaces guide passive motion at the human ankle joint complex is fundamental for the design of relevant surgical treatments. The paper presents a possible improvement of this knowledge by a new kinematic model of the tibiotalar articulation. Passive motion, i.e. in virtually unloaded conditions, was captured in vitro in four lower leg specimens by means of a surgical navigation system with cluster of active markers attached to the tibia and talus. The anatomical geometry of the passive structures, i.e. articular surfaces and attachment areas of the ligaments, were taken by digitisation with a pointer. An equivalent spatial mechanism for the passive motion simulation was defined by three sphere-to-sphere contact points and two rigid links. These contact points were identified at the lateral talo-fibular articulation and at the medial and lateral aspects of the articulation between tibial mortise and trochlea tali. The two rigid links were identified by the isometric fibres at the calcaneofibular and tibiocalcaneal ligaments. An optimisation algorithm was developed for the identification of the final geometrical parameters resulting from an iterative refining process, which targets best matching between model predictions and corresponding experimental measurements of the spatial motion. The specimen-specific equivalent spatial mechanisms replicated the original passive motion very well, with mean discrepancies in position smaller than 2.5 mm and in rotation smaller than 1°. The study demonstrates that the articular surfaces and the ligaments, acting together as a mechanism, control the passive kinematics of the ankle joint.     

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 constraint-based approach to modelling the mobility of the human knee joint

A model of knee mobility able to predict the range and pattern of movement in the unloaded joint was proposed by Wilson et al. (J. Biomech. 31 (1998) 1127-1136). The articular surfaces in the lateral and medial compartments and isometric fascicles in three of the knee ligaments were represented as five constraints on motion between the femur and tibia in a single degree-of-freedom parallel spatial mechanism. The path of movement of the bones during passive flexion was found by solving the forward kinematics of the mechanism using an iterative method. The present paper shows that such a mechanism-based solution approach can lead to an underestimation of the flexion range. This is due to the mechanism reaching a 'stationary configuration' and 'locking'. A new, constraint-based approach to the solution of the model joint displacement is proposed. It avoids the representation of ligaments and articular surfaces by kinematically equivalent chains of one degree-of-freedom pairs which are prone to singularities. It relies instead on a numerical solution of five non-linear constraint equations to find the relative positions of the bones at a series of flexion angles. The method is successful both in its ability to predict motion through a physiological range and in its efficiency with a solution rate forty times faster than the original algorithm. The new approach may be extended to include more complex joint surface geometry, allowing a study of the effects of articular surface shape and ligament arrangement on joint kinematics.     

Ligaments and articular contact guide passive knee flexion

The aim of this study was to test the hypothesis that the coupled features of passive knee flexion are guided by articular contact and by the isometric fascicles of the ACL, PCL and MCL. A three-dimensional mathematical model of the knee was developed, in which the articular surfaces in the lateral and medial compartments and the isometric fascicles in the ACL, PCL and MCL were represented as five constraints in a one degree-of-freedom parallel spatial mechanism. Mechanism analysis techniques were used to predict the path of motion of the tibia relative to the femur. Using a set of anatomical parameters obtained from a cadaver specimen, the model predicts coupled internal rotation and ab/adduction with flexion. These predictions correspond well to measurements of the cadaver specimen’s motion. The model also predicts posterior translation of contact on the tibia with flexion. Although this is a well-known feature of passive knee flexion, the model predicts more translation than has been reported from experiments in the literature. Modelling of uncertainty in the anatomical parameters demonstrated that the discrepancy between theoretical predictions and experimental measurement can be attributed to parameter sensitivity of the model. This study shows that the ligaments and articular surfaces work together to guide passive knee motion. A principal implication of the work is that both articular surface geometry and ligament geometry must be preserved or replicated by surgical reconstruction and replacement procedures to ensure normal knee kinematics and by extension, mechanics.     

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.     

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.     

Subject-specific models of the hindfoot reveal a relationship between morphology and passive mechanical properties

The morphology of the bones, articular surfaces and ligaments and the passive mechanical characteristics of the ankle complex were reported to vary greatly among individuals. The goal of this study was to test the hypothesis that the variations observed in the passive mechanical properties of the healthy ankle complex are strongly influenced by morphological variations. To evaluate this hypothesis six numerical models of the ankle joint complex were developed from morphological data obtained from MRI of six cadaveric lower limbs, and from average reported data on the mechanical properties of ligaments and articular cartilage. The passive mechanical behavior of each model, under a variety of loading conditions, was found to closely match the experimental data obtained from each corresponding specimen. Since all models used identical material properties and were subjected to identical loads and boundary conditions, it was concluded that the observed variations in passive mechanical characteristics were due to variations in morphology, thus confirming the hypothesis. In addition, the average and large variations in passive mechanical behavior observed between the models were similar to those observed experimentally between cadaveric specimens. The results suggest that individualized subject-specific treatment procedures for ankle complex disorders are potentially superior to a one-size-fits-all approach.     

Experimental evaluation of current and novel approximations of articular surfaces of the ankle joint

Kinematics and flexibility properties of both natural and replaced ankle joints are affected by the geometry of the articulating surfaces. Recent studies proposed an original saddle-shaped, skewed, truncated cone with laterally oriented apex, as tibiotalar contact surfaces for ankle prosthesis. The goal of this study was to compare in vitro this novel design with traditional cylindrical or medially centered conic geometries in terms of their ability to replicate the natural ankle joint mechanics. Ten lower limb cadaver specimens underwent a validated process of custom design for the replacement of the natural ankle joint. The process included medical imaging, 3D modeling and printing of implantable sets of artificial articular surfaces based on these three geometries. Kinematics and flexibility of the overall ankle complex, along with the separate ankle and subtalar joints, were measured under cyclic loading. In the neutral and in maximum plantarflexion positions, the range of motion under torques in the three anatomical planes of the three custom artificial surfaces was not significantly different from that of the natural surfaces. In maximum dorsiflexion the difference was significant for all three artificial surfaces at the ankle complex, and only for the cylindrical and medially centered conic geometries at the tibiotalar joint. Natural joint flexibility was restored by the artificial surfaces nearly in all positions. The present study provides experimental support for designing articular surfaces matching the specific morphology of the ankle to be replace, and lays the foundations of the overall process for designing and manufacturing patient-specific total ankle replacements.     

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.     

Articular contact at the tibiotalar joint in passive flexion

The knowledge of the contact areas at the tibiotalar articulating surfaces during passive flexion is fundamental for the understanding of ankle joint mobility. Traditional contact area reports are limited by the invasive measuring techniques used and by the complicated loading conditions applied. In the present study, passive flexion tests were performed on three anatomical preparations from lower leg amputation. Roentgen Stereophotogrammetric Analysis was used to accurately reconstruct the position of the tibia and the talus at a number of unconstrained flexion positions. A large number of points was collected on the surface of the tibial mortise and on the trochlea tali by a 3-D digitiser. Articular surfaces were modelled by thin plate splines approximating these points. Relative positions of these surfaces in all the flexion positions were obtained from corresponding bone position data. A distance threshold was chosen to define contact areas. A consistent pattern of contact was found on the articulating surfaces. The area moved anteriorly on both articular surfaces with dorsiflexion. The average position of the contact area centroid along the tibial mortise at maximum plantarflexion and at maximum dorsiflexion was respectively 58% posterior and 40% anterior of the entire antero-posterior length. For increasing dorsiflexion, the contact area moved from medial to lateral in all the specimens.     

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.     

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

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

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.     

Ligament fibre recruitment and forces for the anterior drawer test at the human ankle joint

Although the anterior drawer test at the ankle joint is commonly used in routine clinical practice, very little is known about the sharing of load between the individual passive structures and the joint response at different flexion angles.A mathematical model of the ankle joint was devised to calculate ligament fibre recruitment and load/displacement curves at different flexion angles. Ligaments were modelled as three-dimensional arrays of fibres, and their orientations at different flexion angles were taken from a previously validated four-bar-linkage model in the sagittal plane. A non-linear stress/strain relationship was assumed for ligament fibres and relevant mechanical parameters were taken from two reports in the literature. Talus and calcaneus were assumed to move as a single rigid body. Antero/distal motion of the talus relative to the tibia was analysed.The ankle joint was found to be stiffer at the two extremes of the flexion range, and the highest laxity was found around the neutral position, confirming previous experimental works. With a first dataset, a 20N anterior force produced 4.3, 5.5, and 4.4mm displacement respectively at 20 degrees plantarflexion, at neutral, and at 20 degrees dorsiflexion. At 10 degrees plantarflexion, for a 6mm displacement, 65% of the external force was supported by the anterior talofibular, 11% by the deep anterior tibiotalar and 5.5% by the tibionavicular ligament. Corresponding results from a second dataset were 1.4, 2.4 and 1.8mm at 40N force, and 80%, 0% and 2% for a 3mm displacement. A component of the contact force supported the remainder.