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.

[1]  J. Archard,et al.  The wear of metals under unlubricated conditions , 1956, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[2]  J. Hammersley MONTE CARLO METHODS FOR SOLVING MULTIVARIABLE PROBLEMS , 1960 .

[3]  Thomas R. Kane,et al.  THEORY AND APPLICATIONS , 1984 .

[4]  K. Johnson Contact Mechanics: Frontmatter , 1985 .

[5]  George E. P. Box,et al.  Empirical Model‐Building and Response Surfaces , 1988 .

[6]  M. Deaton,et al.  Response Surfaces: Designs and Analyses , 1989 .

[7]  Derek J. Pike,et al.  Empirical Model‐building and Response Surfaces. , 1988 .

[8]  G. Wahba Spline models for observational data , 1990 .

[9]  T. Cruse,et al.  Advanced probabilistic structural analysis method for implicit performance functions , 1990 .

[10]  K. An,et al.  Pressure distribution on articular surfaces: application to joint stability evaluation. , 1990, Journal of biomechanics.

[11]  P. Walker,et al.  The Dominance of Cyclic Sliding in Producing Wear in Total Knee Replacements , 1991, Clinical orthopaedics and related research.

[12]  H. Grootenboer,et al.  Articular contact in a three-dimensional model of the knee. , 1991, Journal of Biomechanics.

[13]  Robert Haining,et al.  Statistics for spatial data: by Noel Cressie, 1991, John Wiley & Sons, New York, 900 p., ISBN 0-471-84336-9, US $89.95 , 1993 .

[14]  D. Dowson,et al.  The effect of sliding velocity on the friction and wear of UHMWPE for use in total artificial joints , 1994 .

[15]  Douglas C. Montgomery,et al.  Response Surface Methodology: Process and Product Optimization Using Designed Experiments , 1995 .

[16]  Holger Wendland,et al.  Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree , 1995, Adv. Comput. Math..

[17]  Sheng Chen,et al.  Regularized orthogonal least squares algorithm for constructing radial basis function networks , 1996 .

[18]  R. H. Myers,et al.  Response Surface Methodology: Process and Product Optimization Using Designed Experiments , 1995 .

[19]  D R Broome,et al.  A knee simulating machine for performance evaluation of total knee replacements. , 1997, Journal of biomechanics.

[20]  E. Chao,et al.  A comparison of different methods in predicting static pressure distribution in articulating joints. , 1997, Journal of biomechanics.

[21]  Urmila M. Diwekar,et al.  An efficient sampling technique for off-line quality control , 1997 .

[22]  Vladimir Vapnik,et al.  Statistical learning theory , 1998 .

[23]  M. Pandy,et al.  A Three-Dimensional Musculoskeletal Model of the Human Knee Joint. Part 1: Theoretical Construction , 1997 .

[24]  Federico Girosi,et al.  An Equivalence Between Sparse Approximation and Support Vector Machines , 1998, Neural Computation.

[25]  Donald R. Jones,et al.  Efficient Global Optimization of Expensive Black-Box Functions , 1998, J. Glob. Optim..

[26]  T J Santner,et al.  Robust optimization of total joint replacements incorporating environmental variables. , 1999, Journal of biomechanical engineering.

[27]  Farrokh Mistree,et al.  Statistical Approximations for Multidisciplinary Design Optimization: The Problem of Size , 1999 .

[28]  Salvador Pintos,et al.  Surrogate Modeling-Based Optimization for the Integration of Static and Dynamic Data Into a Reservoir Description , 2002 .

[29]  B. Bai,et al.  Effect of posterior cut angle on tibial component loading. , 2000, The Journal of arthroplasty.

[30]  J D DesJardins,et al.  The use of a force-controlled dynamic knee simulator to quantify the mechanical performance of total knee replacement designs during functional activity. , 2000, Journal of biomechanics.

[31]  Raphael T. Haftka,et al.  Two-level composite wing structural optimization using response surfaces , 2000 .

[32]  D. D’Lima,et al.  Polyethylene Wear and Variations in Knee Kinematics , 2001, Clinical orthopaedics and related research.

[33]  Timothy W. Simpson,et al.  Metamodels for Computer-based Engineering Design: Survey and recommendations , 2001, Engineering with Computers.

[34]  Bernard Grossman,et al.  A Comparison of Global Optimization Methods for the Design of a High-speed Civil Transport , 2001, J. Glob. Optim..

[35]  Edward C. Jones,et al.  Patients’ Expectations of Knee Surgery , 2001, The Journal of bone and joint surgery. American volume.

[36]  I. Clarke,et al.  Effects of A-P translation and rotation on the wear of UHMWPE in a total knee joint simulator. , 2001, Journal of biomedical materials research.

[37]  M. Laurent,et al.  The effect of displacement control input parameters on tibiofemoral prosthetic knee wear , 2001 .

[38]  J. Fisher,et al.  Comparison of wear in a total knee replacement under different kinematic conditions , 2001, Journal of materials science. Materials in medicine.

[39]  T. Simpson,et al.  Comparative studies of metamodelling techniques under multiple modelling criteria , 2001 .

[40]  V. L. Giddings,et al.  Total Knee Replacement Polyethylene Stresses During Loading in a Knee Simulator , 2001 .

[41]  M. Conditt,et al.  What Functional Activities Are Important to Patients With Knee Replacements? , 2002, Clinical orthopaedics and related research.

[42]  Søren Nymand Lophaven,et al.  DACE - A Matlab Kriging Toolbox , 2002 .

[43]  Joseph D. Towles,et al.  Towards a realistic biomechanical model of the thumb: the choice of kinematic description may be more critical than the solution method or the variability/uncertainty of musculoskeletal parameters. , 2003, Journal of biomechanics.

[44]  Thomas J. Santner,et al.  Design and analysis of computer experiments , 1998 .

[45]  Henri Champliaud,et al.  On the Use of Dual Kriging Interpolation for the Evaluation of the Gasket Stress Distribution in Bolted Joints , 2003 .

[46]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[47]  Benjamin J Fregly,et al.  Experimental evaluation of an elastic foundation model to predict contact pressures in knee replacements. , 2003, Journal of biomechanics.

[48]  Andreas Daffertshofer,et al.  PCA in studying coordination and variability: a tutorial. , 2004, Clinical biomechanics.

[49]  Benjamin J Fregly,et al.  Multibody dynamic simulation of knee contact mechanics. , 2004, Medical engineering & physics.

[50]  M. Conditt,et al.  Does Total Knee Replacement Restore Normal Knee Function? , 2005, Clinical orthopaedics and related research.

[51]  Jason P. Halloran,et al.  Comparison of deformable and elastic foundation finite element simulations for predicting knee replacement mechanics. , 2005, Journal of biomechanical engineering.

[52]  Raphael T. Haftka,et al.  Surrogate-based Analysis and Optimization , 2005 .

[53]  Jason P. Halloran,et al.  Explicit finite element modeling of total knee replacement mechanics. , 2005, Journal of biomechanics.

[54]  Benjamin J Fregly,et al.  Computational wear prediction of a total knee replacement from in vivo kinematics. , 2005, Journal of biomechanics.

[55]  Yu Ding,et al.  Optimal Parameter Selection for Electronic Packaging Using Sequential Computer Simulations , 2006 .

[56]  Saikat Pal,et al.  Probabilistic finite element prediction of knee wear simulator mechanics. , 2006, Journal of biomechanics.

[57]  G. Gary Wang,et al.  Review of Metamodeling Techniques in Support of Engineering Design Optimization , 2007, DAC 2006.

[58]  Benjamin J Fregly,et al.  Response surface optimization for joint contact model evaluation. , 2006, Journal of applied biomechanics.

[59]  Timothy M Wright,et al.  Retrieval, experimental, and computational assessment of the performance of total knee replacements , 2006, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[60]  G. Gary Wang,et al.  Review of Metamodeling Techniques in Support of Engineering Design Optimization , 2007 .

[61]  János D. Pintér,et al.  Introduction to Applied Optimization , 2007, Eur. J. Oper. Res..

[62]  Mark Taylor,et al.  Comparison of long-term numerical and experimental total knee replacement wear during simulated gait loading. , 2007, Journal of biomechanics.

[63]  Benjamin J Fregly,et al.  Predicting knee replacement damage in a simulator machine using a computational model with a consistent wear factor. , 2008, Journal of biomechanical engineering.

[64]  B. Fregly,et al.  Sensitivity of knee replacement contact calculations to kinematic measurement errors , 2008, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[65]  Sonja Kuhnt,et al.  Design and analysis of computer experiments , 2010 .