Parametric optimization analysis on a hip prosthesis stem subject to two load scenarios aimed at minimizing its mass while maintaining the stress state below the design limits.

The following article is an extract from a project carried out for the course of Advanced Techniques for the design of Prosthetic Devices. Held for Medical Engineering at the University of Rome Tor Vergata .

Authors: Mastrofini Alessandro & Muscedere Erica


Table of contents

Preliminary analysis

To set up a parametric optimization study it is necessary to define the reference FEM analyzes.
By considering the two load scenarios is possible to set the two analyzes..
The CAD model of the prosthesis is divided into two different materials, titanium alloy Ti6Al4V for the stem and chromium steel for the head, assigning their respective mechanical properties.

Fig. 1: Constraints and conditions of vertical (a) and horizontal (b) load. Zoom on the load application area (c). Difference in discretization between head and stem (d)
Fig. 1: Constraints and conditions of vertical (a) and horizontal (b) load. Zoom on the load application area (c). Difference in discretization between head and stem (d)

In both studies, an interlocking constraint on the spherical head is considered and the forces are applied to the lower end (fig. 1). A denser mesh on the stem is used. A thickness of 1 mm is considered.

From the analyzes the stress is within the design limits (fig. 2), this makes it possible to proceed with an optimization study.

Fig. 2: Results of the preliminary analysis with a thickness of 1 mm. Vertical (a) and horizontal (b) load. The stress on the neck of the stem is evident, at the level of the attachment with the spherical head (c). This may require analysis with more densely packed discretization.
Fig. 2: Results of the preliminary analysis with a thickness of 1 mm. Vertical (a) and horizontal (b) load. The stress on the neck of the stem is evident, at the level of the attachment with the spherical head (c). This may require analysis with more densely packed discretization.

Design of experiment

We then proceed to set up the parametric optimization study.
The first study was carried out with a less fine mesh that does not best represent the stem in thickness. However, it allows you to get a first general overview.

The study considers the thickness of the stem as variable and as constraints the stresses in the two simulations, trying to minimize the mass of the stem. From the results, in table 1 we can see how the optimal scenario is the one with a thickness of 0.9 mm. A total of 11 steps from 0.5mm to 1.5mm with 0.1mm pitch were performed.

Tab. 1: Results of study no. 3. Total of 6 steps. In green the optimal scenario and in red the values outside the project specifications.
Tab. 1: Results of study no. 3. Total of 6 steps. In green the optimal scenario and in red the values outside the project specifications.

Furthermore, it is evident that the Stress2 load scenario is the most worrying condition. In fact, with the same thickness in the horizontal load scenario, the stresses are greater. However, these data need to be double checked. By carrying out a second analysis with a finer mesh (at least 3 elements in the thickness of the stem) it is evident that the stress is actually higher than the specification.
In fig. 3 it is evident that the result of 0.9 mm is not acceptable.

Detailed optimization

Then, a second analysis is performed with this new mesh limiting itself to the worst (horizontal) load scenario. The new mesh starts from a limit for the side tetrahedral element of 0.3 mm (i.e. at least 3 rows in thickness) for the case of 0.9mm and analyzes scenarios of increasing thickness, ensuring minimum accuracy. From this second analysis it is evident that the thickness of 0.9 mm was actually out of the project specifications.

Fig. 3: Comparison of horizontal load with coarse and fine mesh. With the raw mesh (left) the thickness of 0.9 mm was within the project specifications, however it is evident that by refining the mesh (right) the more precise results show an exceeding of the limit. The following analyzes use this fine mesh.
Fig. 3: Comparison of horizontal load with coarse and fine mesh. With the raw mesh (left) the thickness of 0.9 mm was within the project specifications, however it is evident that by refining the mesh (right) the more precise results show an exceeding of the limit. The following analyzes use this fine mesh.

The direction of optimization is evident: adding material and therefore increasing the thickness tends to satisfy the constraint on stresses. Then, starting from a case at the limit of satisfaction, material is added with a finer step. Clearly the analysis stops around the limit conditions.

Tab. 2: Results of study no. 2. Total of 5 steps. In green the optimal scenario and in red the values outside the project specifications.
Tab. 2: Results of study no. 2. Total of 5 steps. In green the optimal scenario and in red the values outside the project specifications.

Then the second analysis faces 5 steps with thickness from 0.9 mm to 1.1 mm with a step of 0.05 mm. The results are shown in table 2. The optimal situation is obtained for the thickness of 1.05 mm.

Insights

Furthermore, it is evident that there is still room for improvement as there is a difference of ∼ 3 g between the scenarios and the tensions are not only still at the limit. Then a further analysis is carried out by pushing to the limit of tolerances (tenth of a gram) and approaching the limit of 500 MPa.
A final analysis is then carried out with a step of 0.01mm from 1 to 1.05 mm. The results in table 3.

Tab. 3: Results of study no. 3. Total of 6 steps. In green the optimal scenario and in red the values outside the project specifications.
Tab. 3: Results of study no. 3. Total of 6 steps. In green the optimal scenario and in red the values outside the project specifications.

The analyzes were then carried out again for the new discretization level, evaluating both load scenarios in the final thickness.
The result is therefore the optimal condition of 1.01 mm thickness with a mass of 266.84 g. We also emphasize that this value is accurate to the design limits, scenarios 3 and 4 of the last simulation are still valid alternatives, with a slight margin for tensions.

Fig. 4: Final confirmatory analysis with 1.01mm thickness and vertical (a) and horizontal (b) load. Section of the stem-head attachment area in the case of greater stress (horizontal load), the area most heavily loaded is the lower part (c). Zoom on the area of maximum strain (d). Strain in the case of vertical (e) and horizontal (f) load.
Fig. 4: Final confirmatory analysis with 1.01mm thickness and vertical (a) and horizontal (b) load. Section of the stem-head attachment area in the case of greater stress (horizontal load), the area most heavily loaded is the lower part (c). Zoom on the area of maximum strain (d). Strain in the case of vertical (e) and horizontal (f) load.

Frequency study

A frequency study was set up for the analysis of the vibration frequencies. A first case was analyzed considering the spherical head with an interlocking constraint. Then a free case was also analyzed, neglecting rigid motions (i.e. first 6 modes with frequency close to zero).

Tab. 4: Resonance frequencies of the prosthesis with a thickness of 1.01 mm with interlocking on the spherical and free head.
Tab. 4: Resonance frequencies of the prosthesis with a thickness of 1.01 mm with interlocking on the spherical and free head.

The resonance frequencies are far from the stimulations to which a prosthesis is subjected (eg cyclicality of a walk / run, a few Hz).

Figure 5: Fundamental modes for the first four resonant frequencies in the case of constraint on the spherical head: 344.07 Hz (a), 480.79 Hz (b), 1.582 kHz (c), 2.999 kHz (d).
Figure 5: Fundamental modes for the first four resonant frequencies in the case of constraint on the spherical head: 344.07 Hz (a), 480.79 Hz (b), 1.582 kHz (c), 2.999 kHz (d).