An extended discrete element method for the estimation of contact pressure at the ankle joint during stance phase

Abnormalities in the ankle contact pressure are related to the onset of osteoarthritis. In vivo measurements are not possible with currently available techniques, so computational methods such as the finite element analysis (FEA) are often used instead. The discrete element method (DEM), a computationally efficient alternative to time-consuming FEA, has also been used to predict the joint contact pressure. It describes the articular cartilage as a bed of independent springs, assuming a linearly elastic behaviour and absence of relative motion between the bones. In this study, we present the extended DEM (EDEM) which is able to track the motion of talus over time. The method was used, with input data from a subject-specific musculoskeletal model, to predict the contact pressure in the ankle joint during gait. Results from EDEM were also compared with outputs from conventional DEM. Predicted values of contact area were larger in EDEM than they were in DEM (4.67 and 4.18 cm2, respectively). Peak values of contact pressure, attained at the toe-off, were 7.3 MPa for EDEM and 6.92 MPa for DEM. Values predicted from EDEM fell well within the ranges reported in the literature. Overall, the motion of the talus had more effect on the extension and shape of the pressure distribution than it had on the magnitude of the pressure. The results indicated that EDEM is a valid methodology for the prediction of ankle contact pressure during daily activities.