A voxel-based formulation for contact finite element analysis

To date, voxel-based finite element models have not been feasible for contact problems, owing to the inherent stair-step boundary discontinuities. New preprocessing techniques are reported herein to mesh these boundaries smoothly, for purposes of contact stress analysis. Further, new techniques are reported to concentrate the mesh resolution automatically near the articular surface, thus reducing the problem size to levels compatible with executing nonlinear problems on contemporary engineering workstations. Close approximations to Hertzian analytical solutions were obtained for spherical and cylindrical geometries meshed in this manner, and an illustrative anatomical contact problem of the human hip joint is presented.