Two-dimensional surrogate contact modeling for computationally efficient dynamic simulation of total knee replacements.
Computational speed is a major limiting factor for performing design sensitivity and optimization studies of total knee replacements. Much of this limitation arises from extensive geometry calculations required by contact analyses. This study presents a novel surrogate contact modeling approach to address this limitation. The approach involves fitting contact forces from a computationally expensive contact model (e.g., a finite element model) as a function of the relative pose between the contacting bodies. Because contact forces are much more sensitive to displacements in some directions than others, standard surrogate sampling and modeling techniques do not work well, necessitating the development of special techniques for contact problems. We present a computational evaluation and practical application of the approach using dynamic wear simulation of a total knee replacement constrained to planar motion in a Stanmore machine. The sample points needed for surrogate model fitting were generated by an elastic foundation (EF) contact model. For the computational evaluation, we performed nine different dynamic wear simulations with both the surrogate contact model and the EF contact model. In all cases, the surrogate contact model accurately reproduced the contact force, motion, and wear volume results from the EF model, with computation time being reduced from 13 min to 13 s. For the practical application, we performed a series of Monte Carlo analyses to determine the sensitivity of predicted wear volume to Stanmore machine setup issues. Wear volume was highly sensitive to small variations in motion and load inputs, especially femoral flexion angle, but not to small variations in component placements. Computational speed was reduced from an estimated 230 h to 4 h per analysis. Surrogate contact modeling can significantly improve the computational speed of dynamic contact and wear simulations of total knee replacements and is appropriate for use in design sensitivity and optimization studies.
Lin
Yi-Chung
YC
Department of Mechanical and Aerospace Engineering, University of Florida, 231 MAE-A Building, P.O. Box 116250, Gainesville, FL 32611-6250.
Haftka
Raphael T
RT
Queipo
Nestor V
NV
Fregly
Benjamin J
BJ
D019645:Arthroplasty, Replacement, Knee / Q000379:methods*; D001696:Biomechanical Phenomena; D003198:Computer Simulation; D019544:Equipment Failure Analysis / Q000379:methods; D020342:Finite Element Analysis; D007720:Knee Prosthesis; D008954:Models, Biological; D009010:Monte Carlo Method; D012680:Sensitivity and Specificity; D013997:Time Factors