An improved method of computing the wear factor for total hip prostheses involving the variation of relative motion and contact pressure with location on the bearing surface

A new method of computing the wear factor for total hip prostheses is presented. In the conventional method, only the resultant contact force and the track drawn by the point of its application are considered so that the product of the instantaneous force and sliding increment is integrated over one motion cycle. In the present, improved, method the contact pressure distribution is discretized by a large number of smaller normal forces, and the contribution of each is summed. This is important because the relative motion and contact pressure vary strongly with location, and because the transverse pressure component is substantial. Hence, the present surface integral represents the large contact surface better than the conventional line integral. A prerequisite for the surface integral was the method of computing the relative motion correctly anywhere on the contact surface, developed and published earlier by the present authors. For the pressure discretization, the contact surface was divided into nearly equal-sized surface elements. The contact pressure was modelled with ellipsoidal, paraboloidal and sinusoidal distributions. Two load cases were studied, double-peak and static. When an ellipsoidal contact pressure distribution extending over a hemisphere was discretized by 1000 element forces, the computed wear factor for double-peak load in a biaxial hip wear simulator was 30% lower than in the conventional resultant force case. The present method can be later developed further to involve the temporal variation of size and location of the contact surface.