Response surface optimization for joint contact model evaluation

When optimization is used to evaluate a joint contact model's ability to reproduce experimental measurements, the high computational cost of repeated contact analysis can be a limiting factor. This paper presents a computationally-efficient response surface optimization methodology to address this limitation. Quadratic response surfaces were fit to contact quantities (contact force, maximum pressure, average pressure, and contact area) predicted by a discrete element contact model of the tibiofemoral joint for various combinations of material modulus and relative bone pose (i.e., position and orientation). The response surfaces were then used as surrogates for costly contact analyses in optimizations that minimized differences between measured and predicted contact quantities. The methodology was evaluated theoretically using six sets of synthetic (i.e., computer-generated) contact data, and practically using one set of experimental contact data. For the synthetic cases, the response surface optimizations recovered all contact quantities to within 3.4% error. For the experimental case, they matched all contact quantities to within 6.3% error except for maximum contact pressure, which was in error by up to 50%. Response surface optimization provides rapid evaluation of joint contact models within a limited range of relative bone poses and can help identify potential weaknesses in contact model formulation and/or experimental data quality.