The role of parameter identification in finite element contact analyses with reference to orthopaedic biomechanics applications.

A finite element model accounting for large sliding frictional contact requires, depending on the type of contact algorithm in use, the definition of many numerical parameters such as contact stiffness, convergence norm and tolerance, compenetration monitoring, over-relaxing factors, etc. All these parameters do not have a physical meaning and thus they cannot be measured experimentally. This makes their identification quite complex. The aim of this study was to investigate the role of parameter identification on the accuracy of results produced by finite element models accounting for bone-implant frictional contact, when the Penalty method is used. The sensitivity analysis of several numerical parameters that may govern the state of results was carried out. Two parameters, contact stiffness and convergence tolerance, were found to play a crucial role in establishing the accuracy of the finite element results. Based on the achieved results it was stated that any numerical-only study involving contact non-linearity and omitting careful qualification of the model limits should be rejected from any peer-reviewed journal.

[1]  N. Langrana,et al.  Modeling of facet articulation as a nonlinear moving contact problem: sensitivity study on lumbar facet response. , 1998, Journal of biomechanical engineering.

[2]  A R Ingraffea,et al.  Mechanical characteristics of the stem‐cement interface , 1991, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[3]  J H Heegaard,et al.  Frictional interface micromotions and anisotropic stress distribution in a femoral total hip component. , 1993, Journal of biomechanics.

[4]  K. Bathe.  Finite Element Procedures , 1995 .

[5]  J H Heegaard,et al.  A finite element model for evaluation of tibial prosthesis-bone interface in total knee replacement. , 1992, Journal of biomechanics.

[6]  A. M. Ahmed,et al.  Micromotion analysis of the fixation of total knee tibial component , 1995 .

[7]  L Cristofolini,et al.  Large-sliding contact elements accurately predict levels of bone-implant micromotion relevant to osseointegration. , 2000, Journal of biomechanics.

[8]  H Weinans,et al.  Trends of mechanical consequences and modeling of a fibrous membrane around femoral hip prostheses. , 1990, Journal of biomechanics.

[9]  P. Walker,et al.  The conflicting requirements of laxity and conformity in total knee replacement. , 1999, Journal of biomechanics.

[10]  M P Whelan,et al.  Finite Element Analysis of a Total Knee Replacement by Using Gauss Point Contact Constraints , 1996, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[11]  H. Saunders.  Book Reviews : NUMERICAL METHODS IN FINITE ELEMENT ANALYSIS K.-J. Bathe and E.L. Wilson Prentice-Hall, Inc, Englewood Cliffs, NJ , 1978 .

[12]  A W Miles,et al.  The influence of the stem-cement interface in total hip replacement—a comparison of experimental and finite element approaches , 1997, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[13]  F W Werner,et al.  Modeling the tensile behavior of the cement-bone interface using nonlinear fracture mechanics. , 1997, Journal of biomechanical engineering.

[14]  D L Bartel,et al.  Coulomb frictional interfaces in modeling cemented total hip replacements: a more realistic model. , 1995, Journal of biomechanics.

[15]  Y Youm,et al.  Three dimensional shape reconstruction and finite element analysis of femur before and after the cementless type of total hip replacement. , 1993, Journal of biomedical engineering.

[16]  L Cristofolini,et al.  Initial stability of a new hybrid fixation hip stem: experimental measurement of implant-bone micromotion under torsional load in comparison with cemented and cementless stems. , 2000, Journal of biomedical materials research.

[17]  C. D. Mote,et al.  Femoral prosthesis implantation induces changes in bone stress that depend on the extent of porous coating , 1994, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[18]  L. Whiteside,et al.  Examination of rotational fixation of the femoral component in total hip arthroplasty. A mechanical study of micromovement and acoustic emission. , 1989, Clinical orthopaedics and related research.

[19]  W. Maloney,et al.  Biomechanical and histologic investigation of cemented total hip arthroplasties. A study of autopsy-retrieved femurs after in vivo cycling. , 1989, Clinical orthopaedics and related research.

[20]  杉山 肇,et al.  Examination of rotational fixation of the femoral component in total hip arthroplasty , 1990 .

[21]  R Huiskes,et al.  Pre-Clinical Testing of Hip Prosthetic Designs: A Comparison of Finite Element Calculations and Laboratory Tests , 1993, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[22]  A Rohlmann,et al.  A nonlinear finite element analysis of interface conditions in porous coated hip endoprostheses. , 1988, Journal of biomechanics.

[23]  J. E. Akin,et al.  Effect of porous coating and loading conditions on total hip femoral stem stability. , 1995, The Journal of arthroplasty.

[24]  M S Hefzy,et al.  Comparison between two techniques for modeling interface conditions in a porous coated hip endoprosthesis. , 1997, Medical engineering & physics.