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