Abstract
Introduction
Stress shielding is one of the major concerns of load bearing implants (e.g. hip prostheses). Stiff implants cause stress shielding, which is thought to contribute to bone resorption1. On the contrary, low-stiffness implants generate high interfacial stresses that have been related to pain and interfacial micro-movements².
Different attempts have been made to reduce these problems by optimizing either the stem design3 or using functionally graded implants (FGI) where the stem's mechanical properties are optimized4. In this way, new additive manufacturing technologies allow fabricating porous materials with well-controlled mesostructure, which allows tailoring their mechanical properties.
In this work, Finite Element (FE) simulations are used to develop an optimization methodology for the shape and material properties of a FGI hip stem. The resorbed bone mass fraction and the stem head displacement are used as objective functions.
Methodology
The 2D-geometry of a femur model (Sawbones®) with an implanted Profemur-TL stem (Wright Medical Technology Inc.) was used for FE simulations. The stem geometry was parameterized using a set of 8 variables (Figure 1-a). To optimize the stem's material properties, a grid was generated with equally spaced points for a total of 96 points (Figure 1-b).
Purely elastic materials were used for the stem and the bone. Two bone qualities were considered: good (Ecortical=20 GPa, Etrabecular=1.5 GPa) and medium (Ecortical=15 GPa, Etrabecular=1 GPa). Poisson ratio was fixed to v=0.3. Loading corresponded to stair climbing. Hip contact force along with abductors, vastus lateralis and vastus medialis muscles were considered5 for a bodyweight of 847 N.
The resorbed bone mass fraction was evaluated from the differences in strain energy densities between the intact bone and the implanted bone2. The displacement of the load point on the femoral head was computed.
The optimization problem was formulated as the minimization of the resorbed bone mass fraction and the head displacement. It was solved using a genetic algorithm.
Results
For the Profemur-TL design, bone resorption was around 36% and 56% for good and medium bone qualities, respectively (Fig. 2). The corresponding head displacements were 11.75 mm and 21.19 mm. Optimized solutions showed bone resorption from 15% to 26% and from 44% to 65% for good and medium bone qualities, respectively. Corresponding head displacements ranged from 11.85 mm to 12.25 mm and from 16.9 mm to 22.6 mm.
Conclusion
The obtained set of solutions constitutes an improvement of the implant performance for this functionally graded implant (FGI) compared to the original implant for both bone qualities. From these simulations, the final solution for the FGI could be chosen based on manufacturing restrictions or another performance indicator.
