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).

The role of altered contact mechanics in the pathogenesis of posttraumatic osteoarthritis (PTOA) following intraarticular fracture remains poorly understood. One proposed etiology is that residual incongruities lead to altered joint contact stresses that, over time, predispose to PTOA. Prevailing joint contact stresses following surgical fracture reduction were quantified in this study using patientspecific contact finite element (FE) analysis. FE models were created for 11 ankle pairs from tibial plafond fracture patients. Both (reduced) fractured ankles and their intact contralaterals were modeled. A sequence of 13 loading instances was used to simulate the stance phase of gait. Contact stresses were summed across loadings in the simulation, weighted by resident time in the gait cycle. This chronic exposure measure, a metric of degeneration propensity, was then compared between intact and fractured ankle pairs. Intact ankles had lower peak contact stress exposures that were more uniform and centrally located. The series-average peak contact stress elevation for fractured ankles was 38% (p ¼ 0.0015; peak elevation was 82%). Fractured ankles had less area with low contact stress exposure than intact ankles and a greater area with high exposure. Chronic contact stress overexposures (stresses exceeding a damage threshold) ranged from near zero to a high of 18 times the matched intact value. The patient-specific FE models represent substantial progress toward elucidating the relationship between altered contact stresses and the outcome of patients treated for intraarticular fractures. ß

A contact finite element (FE) formulation is introduced, amenable to patient-specific analysis of cumulative cartilage mechano-stimulus attributable to habitual functional activity. CT scans of individual human ankles are segmented to delineate bony margins. Each bone surface is projected outward to create a second surface, and the intervening volume is then meshed with continuum hexahedral elements. The tibia is positioned relative to the talus into a weight-bearing apposition. The articular members are first engaged under light preload, then plantar-/dorsi-flexion kinematics and resultant loadings are input for serial FE solutions at 13 instants of the stance phase of level walking gait. Cartilage stress histories are post-processed to recover distributions of cumulative stress-time mechano-stimulus, a metric of degeneration propensity. Consistency in computed contact stress exposures presented for seven intact ankles stood in contrast to the higher magnitude and more focal exposures in an incongruously reduced tibial plafond fracture. This analytical procedure provides patient-specific estimates of degeneration propensity due to various mechanical abnormalities, and it provides a platform from which the mechanical efficacy of alternative surgical interventions can be estimated.

No abstract