Expedited patient-specific assessment of contact stress exposure in the ankle joint following definitive articular fracture reduction

Acute injury severity, altered joint kinematics, and joint incongruity are three important mechanical factors linked to post-traumatic osteoarthritis (PTOA). Finite element analysis (FEA) was previously used to assess the influence of increased contact stress due to joint incongruity on PTOA development. While promising agreement with PTOA development was seen, the inherent complexities of contact FEA limited the numbers of subjects that could be analyzed. Discrete element analysis (DEA) is a simplified methodology for contact stress computation, which idealizes contact surfaces as a bed of independent linear springs. In this study, DEA was explored as an expedited alternative to FEA contact stress exposure computation. DEA was compared to FEA using results from a previously completed validation study of two cadaveric human ankles, as well as a previous study of post-operative contact stress exposure in 11 patients with tibial plafond fracture. DEA-computed maximum contact stresses were within 19% of those experimentally measured, with 90% of the contact area having computed contact stress values within 1 MPa of those measured. In the 11 fractured ankles, maximum contact stress and contact area differences between DEA and FEA were 0.85±0.64 MPa and 22.5±11.5 mm2. As a predictive measure for PTOA development, both DEA and FEA had 100% concordance with presence of OA (KL grade≥2) and >95% concordance with KL grade at 2 years. These results support DEA as a reasonable alternative to FEA for computing contact stress exposures following surgical reduction of a tibial plafond fracture.