Mathematical Development of a Novel Discrete Hip Deformation Algorithm for the In Silico Elasto-Hydrodynamic Lubrication Modelling of Total Hip Replacements
Abstract
:1. Introduction
- When the joint is subjected to intense relative motion and light loading, the hydrodynamic lubrication occurs, so the coupled surfaces are totally separated by the lubricating fluid along the whole contact area;
- The boundary lubrication is established when the high load and the slow relative motion do not allow to obtain the surfaces’ separation, so the lubrication effect is governed by chemical reactions at the contact interface; and the slow relative motion don’t allow to obtain the surfaces’ separation, so the lubrication effect is governed by chemical reactions at the contact interface;
- In the intermediate condition, the load is supported both by asperities in contact and by the lubricating fluid pressure, so it is defined as mixed lubrication.
2. Materials and Methods
2.1. Finite Element Model
- The matrix which relates the nodal pressure vector to the surface pressure field vector in (9);
- The matrix which relates the surface deformation field vector to the nodal normal displacement vector in (10); and
- The matrices and which collect the free part or the constrained one of a vector quantity , useful also to rewrite it in its ordered form through the matrix in (11).
2.2. Acetabular Cup and Femoral Head Meshes
- Regarding the acetabular cup, the nodes on the outer surface are fixed because of the cup fixation with respect to the pelvic bone, namely when , while the nodes subjected to the surface pressure are located on the inner surface, namely when ; and
- Regarding the femoral head, the nodes with the -coordinate lower than a constant value depending on a chosen angle are fixed because of the head fixation with respect to the femoral stem, namely when , while the nodes subjected to the surface pressure are located on the outer surface, namely when .
2.3. Acetabular Cup and Femoral Head Coupling through Cubic Interpolation
- Both for the acetabular cup and for the femoral head, the reference frame used in the finite element discretization was rotated by the inclination angle and the anteversion angle with respect to the anatomical reference frame (Antero/Posterior AP, Proximo/Distal PD and Medio/Lateral ML) through a rotation matrix defined in (22) [3]; and
- The head anatomical reference frame is rotated with respect to the cup one by the Flexion/Extension angle , the Adduction/Abduction angle and the Internal/External Rotation angle through the rotation matrix defined in (23) [3].
3. Results and Discussion
3.1. Validation
3.2. Implementation of the Deformation Model in an EHL Lubrication Algorithm
4. Conclusions
- The implementation of a routine which elaborates the contact pressure in domain zones characterized by nodes overlapping, in order to consider the model functionality within a mixed elasto-hydrodynamic lubrication algorithm;
- The adaptation of the finite element model to consider viscoelastic materials, adding deformation contributions which depend on the time history of the pressure field through a viscosity matrix;
- The usage of the model for other types of joint replacements such as knee, ankle, shoulders, etc.
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
Nomenclature
Linear tetrahedron nodes | |
Tetrahedron face outward-pointing normal | |
Tetrahedron face area | |
Shape function matrix | |
Nodal pressures | |
Nodal pressure vector | |
Nodal force vector | |
Nodal displacement vector | |
Stiffness matrix | |
Nodal normal displacement vector | |
Surface pressure field vector | |
Surface deformation field vector | |
Influence matrix | |
Boundary conditions vector or dimensionless eccentricity vector | |
Cartesian coordinates vector | |
Spherical coordinates | |
Acetabular cup inner radius | |
Acetabular cup thickness | |
Femoral head radius | |
Femoral head fixed nodes angle | |
Surface grid point indices | |
Cubic interpolation matrix from the grid to the grid | |
Rotation matrix | |
Inclination and anteversion angles | |
Hip Flexion/Extension, Adduction/Abduction and Internal/External Rotation angles | |
Radial unit vector | |
Young modulus and Poisson ratio | |
Root Mean Squared Error and squared correlation coefficient | |
Radial clearance | |
Lubricating fluid gap |
References
- Fisher, J.; Dowson, D. Tribology of Total Artificial Joints. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 1991, 205, 73–79. [Google Scholar] [CrossRef]
- Ruggiero, A.; Sicilia, A. Lubrication modeling and wear calculation in artificial hip joint during the gait. Tribol. Int. 2020, 142, 105993. [Google Scholar] [CrossRef]
- Ruggiero, A.; Sicilia, A. A Mixed Elasto-Hydrodynamic Lubrication Model for Wear Calculation in Artificial Hip Joints. Lubricants 2020, 8, 72. [Google Scholar] [CrossRef]
- Ruggiero, A.; Sicilia, A.; Affatato, S. In silico total hip replacement wear testing in the framework of ISO 14242-3 accounting for mixed elasto-hydrodynamic lubrication effects. Wear 2020, 460, 203420. [Google Scholar] [CrossRef]
- Mattei, L.; Di Puccio, F.; Piccigallo, B.; Ciulli, E. Lubrication and wear modelling of artificial hip joints: A review. Tribol. Int. 2011, 44, 532–549. [Google Scholar] [CrossRef]
- Gleghorn, J.P.; Bonassar, L.J. Lubrication mode analysis of articular cartilage using Stribeck surfaces. J. Biomech. 2008, 41, 1910–1918. [Google Scholar] [CrossRef]
- Scholes, S.C.; Unsworth, A. Comparison of friction and lubrication of different hip prostheses. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2000, 214, 49–57. [Google Scholar] [CrossRef]
- Ruggiero, A.; Merola, M.; Affatato, S. Finite element simulations of hard-on-soft hip joint prosthesis accounting for dynamic loads calculated from a musculoskeletal model during walking. Materials 2018, 11, 574. [Google Scholar] [CrossRef] [Green Version]
- Affatato, S.; Ruggiero, A. A Perspective on Biotribology in Arthroplasty: From In Vitro toward the Accurate In Silico Wear Prediction. Appl. Sci. 2020, 10, 6312. [Google Scholar] [CrossRef]
- Jalali-Vahid, D.; Jagatia, M.; Jin, Z.M.; Dowson, D. Prediction of lubricating film thickness in a ball-in-socket model with a soft lining representing human natural and artificial hip joints. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2001, 215, 363–372. [Google Scholar] [CrossRef]
- Jalali-Vahid, D.; Jin, Z.M.; Dowson, D. Elastohydrodynamic lubrication analysis of hip implants with ultra high molecular weight polyethylene cups under transient conditions. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2003, 217, 767–777. [Google Scholar] [CrossRef]
- Jalali-Vahid, D.; Jagatia, M.; Jin, Z.; Dowson, D. Prediction of lubricating film thickness in UHMWPE hip joint replace-ments. J. Biomech. 2001, 34, 261–266. [Google Scholar] [CrossRef]
- Jalali-Vahid, D.; Jin, Z.M. Transient elastohydrodynamic lubrication analysis of ultra-high molecular weight polyethylene hip joint replacements. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2001, 216, 409–420. [Google Scholar] [CrossRef]
- Fregly, B.J.; Bei, Y.; Sylvester, M.E. Experimental evaluation of an elastic foundation model to predict contact pressures in knee replacements. J. Biomech. 2003, 36, 1659–1668. [Google Scholar] [CrossRef]
- Halloran, J.P.; Easley, S.K.; Petrella, A.J.; Rullkoetter, P.J. Comparison of deformable and elastic foundation finite element simulations for predicting knee replacement mechanics. J. Biomech. Eng. 2005, 127, 813–818. [Google Scholar] [CrossRef]
- Askari, E.; Andersen, M.S. A closed-form formulation for the conformal articulation of metal-on-polyethylene hip prostheses: Contact mechanics and sliding distance. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2018, 232, 1196–1208. [Google Scholar] [CrossRef] [Green Version]
- Pérez-González, A.; Fenollosa-Esteve, C.; Sancho-Bru, J.L.; Sánchez-Marín, F.T.; Vergara, M.; Rodríguez-Cervantes, P.J. A modified elastic foundation contact model for application in 3D models of the prosthetic knee. Med. Eng. Phys. 2008, 30, 387–398. [Google Scholar] [CrossRef]
- Mukras, S.M.; Kim, N.H.; Schmitz, T.L.; Sawyer, W.G. Evaluation of contact force and elastic foundation models for wear analysis of multibody systems. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Montreal, QC, Canada, 15–18 August 2010; Volume 44106, pp. 1565–1577. [Google Scholar]
- Srivastava, G.; Christian, N.; Higgs, C.F., III. A predictive framework of the tribological impact of physical activities on metal-on-plastic hip implants. Biotribology 2021, 25, 100156. [Google Scholar] [CrossRef]
- Wang, F.; Jin, Z. Prediction of elastic deformation of acetabular cups and femoral heads for lubrication analysis of artificial hip joints. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2004, 218, 201–209. [Google Scholar] [CrossRef]
- Evans, H.P.; Hughes, T.G. Evaluation of deflection in semi-infinite bodies by a differential method. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2000, 214, 563–584. [Google Scholar] [CrossRef]
- Gao, L.; Dowson, D.; Hewson, R.W. Predictive wear modeling of the articulating metal-on-metal hip replacements. J. Biomed. Mater. Res. Part B Appl. Biomater. 2015, 105, 497–506. [Google Scholar] [CrossRef]
- Jagatia, M.; Jin, Z. Elastohydrodynamic lubrication analysis of metal-on-metal hip prostheses under steady state entraining motion. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2001, 215, 531–541. [Google Scholar] [CrossRef]
Parameter | Value |
---|---|
Femoral head radius r | 14 mm |
Acetabular cup inner radius R | 14.03 mm |
Acetabular cup thickness H | 9 mm |
Femoral head fixed nodes angle | 30° |
Acetabular cup Young modulus | 1.1 GPa |
Acetabular cup Poisson ratio | 0.42 |
Femoral head Young modulus | 245.9 GPa |
Femoral head Poisson ratio | 0.24 |
Pressure gaussian peak | 107 Pa |
Pressure gaussian θ-translation | 2π/3 rad |
Pressure gaussian φ-translation | π/3 rad |
Pressure gaussian dimensionless width | 0.2 |
Parameter | Value |
---|---|
Pressure gaussian peak | 107 Pa |
Pressure gaussian θ-translation | π/4 rad |
Pressure Gaussian φ-translation | π/4 rad |
Pressure gaussian dimensionless width | 0.2 |
Anteversion angle | 90° |
Inclination angle | 45° |
Flexion/Extension angle | −40° |
Adduction/Abduction angle | 10° |
Internal/External Rotation angle | 2° |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Ruggiero, A.; Sicilia, A. Mathematical Development of a Novel Discrete Hip Deformation Algorithm for the In Silico Elasto-Hydrodynamic Lubrication Modelling of Total Hip Replacements. Lubricants 2021, 9, 41. https://doi.org/10.3390/lubricants9040041
Ruggiero A, Sicilia A. Mathematical Development of a Novel Discrete Hip Deformation Algorithm for the In Silico Elasto-Hydrodynamic Lubrication Modelling of Total Hip Replacements. Lubricants. 2021; 9(4):41. https://doi.org/10.3390/lubricants9040041
Chicago/Turabian StyleRuggiero, Alessandro, and Alessandro Sicilia. 2021. "Mathematical Development of a Novel Discrete Hip Deformation Algorithm for the In Silico Elasto-Hydrodynamic Lubrication Modelling of Total Hip Replacements" Lubricants 9, no. 4: 41. https://doi.org/10.3390/lubricants9040041
APA StyleRuggiero, A., & Sicilia, A. (2021). Mathematical Development of a Novel Discrete Hip Deformation Algorithm for the In Silico Elasto-Hydrodynamic Lubrication Modelling of Total Hip Replacements. Lubricants, 9(4), 41. https://doi.org/10.3390/lubricants9040041