Physical validation of a patient-specific contact finite element model of the ankle

A validation study was conducted to determine the extent to which computational ankle contact finite element (FE) results agreed with experimentally measured tibio-talar contact stress. Two cadaver ankles were loaded in separate test sessions, during which ankle contact stresses were measured with a high-resolution (Tekscan) pressure sensor. Corresponding contact FE analyses were subsequently performed for comparison. The agreement was good between FE-computed and experimentally measured mean (3.2% discrepancy for one ankle, 19.3% for the other) and maximum (1.5% and 6.2%) contact stress, as well as for contact area (1.7% and 14.9%). There was also excellent agreement between histograms of fractional areas of cartilage experiencing specific ranges of contact stress. Finally, point-by-point comparisons between the computed and measured contact stress distributions over the articular surface showed substantial agreement, with correlation coefficients of 90% for one ankle and 86% for the other. In the past, general qualitative, but little direct quantitative agreement has been demonstrated with articular joint contact FE models. The methods used for this validation enable formal comparison of computational and experimental results, and open the way for objective statistical measures of regional correlation between FE-computed contact stress distributions from comparison articular joint surfaces (e.g., those from an intact versus those with residual intra-articular fracture incongruity).     

Some effects of global joint morphology on local stress aberrations near imprecisely reduced intra-articular fractures

Cadaver models of contact pressure aberration near displaced intra-articular fractures complement clinical experience, but inter-specimen variability often complicates interpretation of in vitro data. A contact finite element formulation is here used to study juxta-articular stress distributions in a plane strain model of tibial plateau step-off incongruity. Attention is focused on the influence of global morphologic parameters: intact joint surface curvatures, cartilage thickness, and cartilage stiffness. The computed stress distributions agreed well with experimental recordings for a typical 3 mm incongruity in an otherwise normal joint. Both decreased cartilage thickness and increased cartilage modulus led to elevations in the peak local contact stress, and to concentration of contact stress near the edge of the step-off incongruity. Similar effects were seen when reduction of global joint congruency was modelled by decreasing the concavity of the tibial plateau. While the observed degree of coupling between global morphology and local stress aberration was by no means negligible, the sensitivity of the stresses to variations in individual parameters was relatively mild. This suggests that the finite element results will be useful for experimental data interpretation.     

A contact algorithm for density-based load estimation

An algorithm, which includes contact interactions within a joint, has been developed to estimate the dominant loading patterns in joints based on the density distribution of bone. The algorithm is applied to the proximal femur of a chimpanzee, gorilla and grizzly bear and is compared to the results obtained in a companion paper that uses a non-contact (linear) version of the density-based load estimation method. Results from the contact algorithm are consistent with those from the linear method. While the contact algorithm is substantially more complex than the linear method, it has some added benefits. First, since contact between the two interacting surfaces is incorporated into the load estimation method, the pressure distributions selected by the method are more likely indicative of those found in vivo. Thus, the pressure distributions predicted by the algorithm are more consistent with the in vivo loads that were responsible for producing the given distribution of bone density. Additionally, the relative positions of the interacting bones are known for each pressure distribution selected by the algorithm. This should allow the pressure distributions to be related to specific types of activities. The ultimate goal is to develop a technique that can predict dominant joint loading patterns and relate these loading patterns to specific types of locomotion and/or activities.     

The transmission of load through the human hip joint

This paper describes the results of loading experiments carried out on human hip joints. The unloaded surfaces of the femoral head and the acetabulum are slightly incongruous. The location and magnitude of the contact areas between the surfaces therefore depend on the magnitude and direction of the applied load. The contact areas were determined experimentally for a variety of loads typical of normal walking. Two distinct contact areas were found on the anterior and posterior aspects of the acetabulum at light loads, gradually merging with increasing load until, at a certain transition load, the dome of the acetabulum comes into contact and contact is then complete. The value of the transition load depends on the rate of loading, due to creep of the cartilage, and was found to vary from 50 per cent of body weight in young specimens to 25 per cent of body weight for elderly specimens for rates of loading typical of normal walking. Thus, the dome of the acetabulum is out of contact for a substantial portion of the swing phase of normal walking.

The analysis of a much simplified model of the hip joint is presented. The dependence of contact area on load is demonstrated, but also a method of determining the transition load for complete contact from the load/deflection relation for the hip is suggested. The values of the transition load quoted above were obtained by this method. The analysis further indicates that the distribution of pressure between the articular surfaces depends critically on the distribution of cartilage thickness throughout the joint. It is suggested that the distribution of cartilage thickness is such as to lead to a state of uniform pressure at the upper end of the physiological load range. Some experimental evidence is presented in support of this suggestion.

It is concluded that the function of joint incongruity is to allow the articular surfaces to come out of contact at light loads so that the cartilage may be exposed to synovial fluid for the purposes of nutrition and lubrication. At large loads, the distribution of cartilage thickness ensures that a state of hydrostatic pressure is achieved in order that cartilage, with a large fluid content, may transmit large pressures without flow and consequent loss of its integrity.     

Effect of subtalar arthroereisis on the tibiotalar contact characteristics in a cadaveric flatfoot model

Previous studies reported the effect of flatfoot deformity on tibiotalar joint contact characteristics. The lateral shift of the load which occurred in flatfeet may be responsible for degenerative changes in the ankle joint. The purpose was to assess the pattern of joint contact stress of the tibiotalar joint in intact, flat, and corrected specimens with subtalar arthroereisis. Seven fresh-frozen cadaver specimens were studied in the intact and flat-footed condition after transection of ligaments which support the medial arch. Ankle joint contact stress and plantar pressure patterns were determined from a capacitive pressure sensor inserted in the tibiotalar joint and a pressure distribution platform when the specimens were axially loaded in simulated mid-stance. Contact pressure was also assessed after subtalar arthroereisis with a 12 mm Kalix implant for correction of the flatfoot deformity. A maximum contact pressure of 1414.2±319.2 kPa was recorded in the middle-medial region in intact specimens. Flatfoot caused a lateral shift in the pressure distribution (p<0.05). In the flat specimens, the maximum contact pressure of 1394.7 8±470.5 kPa was in the anterior-central region. After subtalar arthroereisis with the Kalix implant for correction of the flatfoot deformity a maximum contact pressure of 1323.3±497.5 kPa was observed in the middle-lateral region. In a cadaver model, subtalar arthroereisis with Kalix implant failed to restore a normal intraarticular ankle joint pressure pattern. Further interventions should be considered to restore a normal pressure pattern.     

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.     

Response surface optimization for joint contact model evaluation

When optimization is used to evaluate a joint contact model's ability to reproduce experimental measurements, the high computational cost of repeated contact analysis can be a limiting factor. This paper presents a computationally-efficient response surface optimization methodology to address this limitation. Quadratic response surfaces were fit to contact quantities (contact force, maximum pressure, average pressure, and contact area) predicted by a discrete element contact model of the tibiofemoral joint for various combinations of material modulus and relative bone pose (i.e., position and orientation). The response surfaces were then used as surrogates for costly contact analyses in optimizations that minimized differences between measured and predicted contact quantities. The methodology was evaluated theoretically using six sets of synthetic (i.e., computer-generated) contact data, and practically using one set of experimental contact data. For the synthetic cases, the response surface optimizations recovered all contact quantities to within 3.4% error. For the experimental case, they matched all contact quantities to within 6.3% error except for maximum contact pressure, which was in error by up to 50%. Response surface optimization provides rapid evaluation of joint contact models within a limited range of relative bone poses and can help identify potential weaknesses in contact model formulation and/or experimental data quality.     

Finite element analysis of a total ankle replacement during the stance phase of gait

Total ankle replacement (TAR) designs have still several important issues to be addressed before the treatment becomes fully acceptable clinically. Very little is known about the performance, in terms of the contact pressures and kinematics of TAR when subjected to daily activities such as level gait. For this purpose, an explicit finite element model of a novel 3-component TAR was developed, which incorporated a previously validated mechanical model of the ankle ligament apparatus. The intermediate mobile polyethylene meniscal bearing was modelled as an elastic-plastic continuum while the articulating surfaces of the tibial and talar metal components as rigid bodies. Overall kinematics, contact pressures and ligament forces were analysed during passive, i.e. virtually unloaded, and active, i.e. stance phase of gait, conditions. Simulation of passive motion predicted similar kinematics as reported previously in an analytical four-bar linkage model. The meniscal bearing was observed to move 5.6 mm posteriorly during the simulated stance and the corresponding antero-posterior displacement of the talar component was 8.3 mm. The predicted pattern and the amount (10.6 degrees ) of internal-external rotation of the ankle complex were found to be in good agreement with corresponding in vivo measurements on normal ankles. A peak contact pressure of 16.8 MPa was observed, with majority of contact pressures below 10 MPa. For most ligaments, reaction forces remain within corresponding physiological ranges. A first realistic representation of the biomechanical behaviour of the human ankle when replaced by prosthetic joints is provided. The applied methodology can potentially be applied to other TAR designs.     

Lubrication of the human ankle joint in walking with the synovial fluid filtrated by the cartilage with the surface zone worn out: steady pure sliding motion.

A mixture model of synovial fluid filtration by cartilage in the human ankle joint during walking is presented for steady sliding motion of the articular surfaces. In the paper the cartilage surface zone is assumed worn out. The same model has been recently applied to the squeeze-film problem for the human hip joint loaded by the body weight during standing (Hlavácek, Journal of Biomechanics 26, 1145-1150, 1151-1160, 1993; Hlavácek and Novák, Journal of Biomechanics 28, 1193-1198, 1199-1205, 1995). The linear biphasic model for cartilage (elastic porous matrix + ideal fluid) due to Prof. V. C. Mow and his co-workers and the biphasic model for synovial fluid (viscous fluid + ideal fluid), as used in the above-mentioned squeeze-film problem, are applied. For the physiologic parameters of the ankle joint during walking, a continuous synovial fluid film about 1 microm thick is maintained under steady entraining motion according to the classical model without the fluid transport across the articular surface. This is not the case in the filtration model with the cartilage surface zones worn out. On the contrary, this filtration model indicates that synovial fluid is intensively filtrated by such cartilage, so that no continuous fluid film is maintained and a synovial gel layer, about 10(-8) m thick, develops over the majority of the contact. Thus, if the cartilage surface zones are worn out, boundary lubrication should prevail in the ankle joint under steady sliding motion for the mean values of loading and the sliding velocity encountered in walking cycle.     

Computational biomechanics of articular cartilage of human knee joint: Effect of osteochondral defects

Articular cartilage and its supporting bone functional conditions are tightly coupled as injuries of either adversely affects joint mechanical environment. The objective of this study was set to quantitatively investigate the extent of alterations in the mechanical environment of cartilage and knee joint in presence of commonly observed osteochondral defects. An existing validated finite element model of a knee joint was used to construct a refined model of the tibial lateral compartment including proximal tibial bony structures. The response was computed under compression forces up to 2000 N while simulating localized bone damage, cartilage–bone horizontal split, bone overgrowth and absence of deep vertical collagen fibrils.Localized tibial bone damage increased overall joint compliance and substantially altered pattern and magnitude of contact pressures and cartilage strains in both tibia and femur. These alterations were further exacerbated when bone damage was combined with base cartilage split and absence of deep vertical collagen fibrils. Local bone boss markedly changed contact pressures and strain patterns in neighbouring cartilage. Bone bruise/fracture and overgrowth adversely perturbed the homeostatic balance in the mechanical environment of articulate cartilage surrounding and opposing the lesion as well as the joint compliance. As such, they potentially contribute to the initiation and development of post-traumatic osteoarthritis.     

Contact pressure in the facet joint during sagittal bending of the cadaveric cervical spine

The facet joint contributes to the normal biomechanical function of the spine by transmitting loads and limiting motions via articular contact. However, little is known about the contact pressure response for this joint. Such information can provide a quantitative measure of the facet joint's local environment. The objective of this study was to measure facet pressure during physiologic bending in the cervical spine, using a joint capsule-sparing technique. Flexion and extension bending moments were applied to six human cadaveric cervical spines. Global motions (C2-T1) were defined using infra-red cameras to track markers on each vertebra. Contact pressure in the C5-C6 facet was also measured using a tip-mounted pressure transducer inserted into the joint space through a hole in the postero-inferior region of the C5 lateral mass. Facet contact pressure increased by 67.6 ± 26.9 kPa under a 2.4 Nm extension moment and decreased by 10.3 ± 9.7 kPa under a 2.7 Nm flexion moment. The mean rotation of the overall cervical specimen motion segments was 9.6 ± 0.8° and was 1.6 ± 0.7° for the C5-C6 joint, respectively, for extension. The change in pressure during extension was linearly related to both the change in moment (51.4 ± 42.6 kPa/Nm) and the change in C5-C6 angle (18.0 ± 108.9 kPa/deg). Contact pressure in the inferior region of the cervical facet joint increases during extension as the articular surfaces come in contact, and decreases in flexion as the joint opens, similar to reports in the lumbar spine despite the difference in facet orientation in those spinal regions. Joint contact pressure is linearly related to both sagittal moment and spinal rotation. Cartilage degeneration and the presence of meniscoids may account for the variation in the pressure profiles measured during physiologic sagittal bending. This study shows that cervical facet contact pressure can be directly measured with minimal disruption to the joint and is the first to provide local pressure values for the cervical joint in a cadaveric model.     

Mechanics of the anterior drawer test at the ankle: the effects of ligament viscoelasticity

The anterior drawer test at the human ankle joint is a routine clinical examination. The relationship between the mechanical response of this joint and the flexion angle was elucidated by a recent mathematical model, using purely elastic mechanical characteristics for the ligament fibres. The objective of the present work was to assess the effect of ligament viscoelasticity on the force response of the ankle joint for anterior displacements of the foot relative to the tibia, at different ankle flexion positions. A viscoelastic model of the ligaments from the literature was included in the recently proposed mathematical model. Drawer tests were simulated at several flexion angles and for increasing velocities of the imposed anterior displacement. The stiffness of the model ankle joint increased only modestly with velocity. The response force found for a 6mm displacement at 20 degrees plantarflexion increased by only 13% for a one hundred-fold increase in velocity from 0.1 to 10 mm/s. The flexion angle was confirmed as the most influential parameter in the mechanical response of the ankle to anterior drawer test.     

Mechanical responses and integrin associated protein expression by human ankle chondrocytes

Metabolic, biochemical and biomechanical differences between ankle and knee joint cartilage and chondrocytes including resistance to the effects of catabolic cytokines and fibronectin fragments may be relevant to differences in prevalence of OA in these joints. Although there is increasing information available on how chondrocytes from knee and hip joint cartilage recognise and respond to mechanical stimuli, knowledge of mechanotransduction in ankle joint chondrocytes is limited. This study was undertaken to (i) establish whether the response of normal ankle joint derived chondrocytes to mechanical stimulation in vitro was similar to that of normal and osteoarthritic knee joint derived chondrocytes and (ii) to investigate whether these chondrocytes showed differences in expression of integrin associated regulatory and signalling molecules. Unlike normal knee joint chondrocytes, ankle joint chondrocytes did not show an increase in relative levels of aggrecan mRNA when mechanically stimulated. No obvious change in protein tyrosine phosphorylation was seen in ankle chondrocytes subsequent to mechanical stimulation but these cells expressed elevated levels of tyrosine phosphorylated proteins at rest when compared to normal knee joint chondrocytes. Ankle joint chondrocytes showed an increase in protein kinase B phosphorylation following 1 min 0.33 Hz stimulation which was inhibited by the presence of antibodies to alpha5beta1 integrin. Ankle joint chondrocytes appeared to show significant differences in levels of the integrin-associated proteins CD98, CD147 and galectin 3, PKCgamma and differences in responses to glutamate were seen. Chondrocytes from ankle and knee joint cartilage respond differently to 0.33 Hz mechanical stimulation. This may be related to modified integrin-dependent mechanotransduction as a result of changes in expression of integrin regulatory molecules such as CD98 or differential expression and function of downstream components of the mechanotransduction pathway such as PKC or NMDA receptors.     

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.     

Myoelectric stimulation on peroneal muscles resists simulated ankle sprain motion

The inadequate reaction time of the peroneal muscles in response to an incorrect foot contact event has been proposed as one of the etiological factors contributing to ankle joint inversion injury. Thus, the current study aimed to investigate the efficacy of a myoelectric stimulation applied to the peroneal muscles in the prevention of a simulated ankle inversion trauma. Ten healthy male subjects performed simulated inversion and supination tests on a pair of mechanical sprain simulators. An electrical signal was delivered to the peroneal muscles of the subjects through a pair of electrode pads. The start of the stimulus was synchronized with the drop of the sprain simulator's platform. In order to determine the maximum delay time which the stimulus could still resist the simulated ankle sprain motion, different delay time were test (0, 5, 10, and 15ms). Together with the control trial (no stimulus), there were 5 testing conditions for both simulated inversion and supination test. The effect was quantified by the drop in maximum ankle tilting angle and angular velocity, as determined by a motion analysis system with a standard laboratory procedure. Results showed that the myoelectric stimulation was effective in all conditions except the one with myoelectric stimulus delayed for 15ms in simulated supination test. It is concluded that myoelectric stimulation on peroneal muscles could resist an ankle spraining motion.     

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.     

Ankle intrinsic stiffness changes with postural sway

In standing, the human body is inherently unstable and its stabilization requires constant regulation of ankle torque, generated by a combination of ankle intrinsic properties, peripheral reflexes, and central contributions. Ankle intrinsic stiffness, which quantifies the joint intrinsic properties, has been usually assumed constant in standing; however, there is strong evidence that it is highly dependent on the joint torque, which changes significantly with sway in stance. In this study, we examined how ankle intrinsic stiffness changes with postural sway during standing. Ten subjects stood on a standing apparatus, while subjected to pulse perturbations of ankle position. The mean torque of a short period before the start of each pulse was used as a measure of background torque. Responses with similar background torques were grouped together and used to estimate the parameters of an intrinsic stiffness model. Stiffness estimates were normalized to the critical stiffness and the background torque was transformed to the center of pressure location. We found that in most subjects, the normalized stiffness increased linearly with the movement of center of pressure towards the toes, with an average slope of 2.11 ± 0.80 1/m·rad. This modulation of ankle intrinsic stiffness seems functionally appropriate, since the intrinsic stiffness increases quickly, as the center of pressure moves toward the toes and the limits of stability. These large changes of ankle intrinsic stiffness with postural sway must be incorporated in any model of stance control.     

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.     

Lubricin: a novel potential biotherapeutic approaches for the treatment of osteoarthritis

Osteoarthritis (OA) is a multi-factor disorder of sinovial joints, which characterized by escalated degeneration and loss of articular cartilage. Treatment of OA is a critical unmet need in medicine for regeneration of damaged articular cartilage in elderly. On the other hand, lubricin, a glycoprotein specifically synthesized by chondrocytes located at the surface of articular cartilage, has been shown to provide boundary lubrication of congruent articular surfaces under conditions of high contact pressure and near zero sliding speed. Lubrication of these surfaces is critical to normal joint function, while different gene expressions of lubricin had been found in the synovium of rheumatoid arthritis (RA) and OA. Moreover, mutations or lacking of lubricin gene have been shown to link to the joint disease such as camptodactyly-arthropathy-coxa vara-pericarditis syndrome (CACP), synovial hyperplasia and failure of joint function, suggesting an important role of lubricin in the pathogenesis of these joint disease. Recent studies demonstrate that administration with recombinant lubricin in the joint cavity would be effective in the prevention of cartilage degeneration in animal OA models. Therefore, a treatment with lubricin which would protect cartilage in vivo would be desirable. This article reviews recent findings with regard to the possible role of lubricin in the progression of OA, and further discusses lubricin as a novel potential biotherapeutic approaches for the treatment of OA.