Computational and Mathematical Methods in Medicine

Computational and Mathematical Methods in Medicine / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 6183679 |

Jesica Anguiano-Sanchez, Oscar Martinez-Romero, Hector R. Siller, Jose A. Diaz-Elizondo, Eduardo Flores-Villalba, Ciro A. Rodriguez, "Influence of PEEK Coating on Hip Implant Stress Shielding: A Finite Element Analysis", Computational and Mathematical Methods in Medicine, vol. 2016, Article ID 6183679, 10 pages, 2016.

Influence of PEEK Coating on Hip Implant Stress Shielding: A Finite Element Analysis

Academic Editor: Valeri Makarov
Received06 Oct 2015
Revised20 Jan 2016
Accepted07 Feb 2016
Published14 Mar 2016


Stress shielding is a well-known failure factor in hip implants. This work proposes a design concept for hip implants, using a combination of metallic stem with a polymer coating (polyether ether ketone (PEEK)). The proposed design concept is simulated using titanium alloy stems and PEEK coatings with thicknesses varying from 100 to 400 μm. The Finite Element analysis of the cancellous bone surrounding the implant shows promising results. The effective von Mises stress increases between 81 and 92% for the complete volume of cancellous bone. When focusing on the proximal zone of the implant, the increased stress transmission to the cancellous bone reaches between 47 and 60%. This increment in load transferred to the bone can influence mineral bone loss due to stress shielding, minimizing such effect, and thus prolonging implant lifespan.

1. Introduction

The number of total hip arthroplasties (THA) operations is increasing, reaching more than one million procedures worldwide per year. This technique is a useful treatment option for osteoarthritis and rheumatoid arthritis on the hip joint, allowing the patients to regain pain-free mobility [1]. Although THA is considered a successful procedure, recent projections indicate the number of revision surgeries is expected to increase by 137% in the next 15 years [2]. This is a major problem due to the pain and high costs caused to the patient, together with less favorable results compared to the first procedure, mainly because of the damage cause to the remaining bone after the THA. Hip implants are designed to last for at least 20 years, but their lifespan has reduced by several problems. One of the most commonly recorded indications for revision surgery is aseptic loosening, where stress shielding is a principal factor [3].

According to Wolff’s law, bones adapt to the mechanical load they receive. When a person is more active in a specific part of the body, more bone is added to strengthen the tissue, and conversely, if a bone stops receiving load for a prolonged time, the mass of the tissue decreases and bone is lost. Once the hip replacement is conducted, the load is carried mainly by the implant itself and not by the femur. This phenomenon is due to a mismatch in stiffness between the hip implant and femur (almost 10 times higher in implant), with variations related to natural physiological conditions [4]. An insufficient load transfer between bone and implant leads to mineral bone lost and thus to lack of contact between the bone and femur. This effect is known as stress shielding.

The research literature shows two approaches to tackle stress shielding in hip implants: design and/or materials. Several studies have focused on changing the geometry of the hip implant in order to reduce the stress shielding effect [58]. Joshi et al. propose a new design and proximal fixation method to reduce stress shielding [5]. Gross and Abel use numerical analysis to show the benefits of using a hollow hip implant design [6]. In order to understand stress shielding, Boyle and Kim analyzed commercially available hip implants with consideration of microlevel bone remodeling [7]. Hirata et al. show that hip implant geometry plays a role in stress shielding, measuring the bone mineral density in patients over a period of one year [8]. The proposed hip implant concept in the study reported here is not based on a geometric approach.

Similar to geometric approach, the literature reports several studies that focus on biomimetic materials for hip implants in order to reduce stress shielding. Bougherara et al. proposed a hip implant based on polymeric composite and a hydroxyapatite-based coating [9]. Oshkour et al. propose a functionally graded hip implant based on stainless steel, titanium alloy, and hydroxyapatite. The simulation results show improvements in the stress shielding [3]. More recently, Tavakkoli Avval et al. showed the benefits of polymeric composite hip implants by conducting a coupled simulation of the bone-implant interaction biomechanics and the bone regeneration process [10].

The study reported here proposes a design concept that combines a metallic stem with an engineered polymer coating. This approach has not been reported in the relevant literature. The proposed polymer has Young’s modulus similar to the bone, facilitating load transmission, and therefore, reducing the stress shielding effect. Finite Element analysis provides preliminary validation of the proposed hip implant design concept.

2. Materials and Methods

2.1. Hip Implant Design Concept

The proposed design concept builds on existing technology for press-fit hip implants (see Figure 1). The stem is made out of medical grade of titanium alloy Ti-6Al-4V. The selected material for the coating is polyether ether ketone (PEEK), a well-known biocompatible polymer, used in orthopedic, spine, and dental implants. PEEK and polyaryl ether ketones (PAEKs) had been used as biomaterials since the 1980s [11], due to their structure that confers outstanding chemical resistance, inertness, and thermal stability for in vivo conditions. Additionally, for the purpose of this study, PEEK has Young’s modulus comparable to that of the bone. Therefore, there is potential to reduce the stiffness mismatch between the femur and the hip implant.

This work investigates how a PEEK coating on a titanium alloy hip implant stem could improve the effective von Mises stress distribution on the femur, comparing a model with uncoated condition to models that have coatings with different thicknesses. An increase in effective von Mises stress values is expected, reducing stress shielding on the femur.

The coating is derived from a hypothetical manufacturing process. The coating thickness is assumed to be uniformly distributed along the surface of the implant stem, starting at the height of lesser trochanter. Figure 1 shows the PEEK coating representation and the modified Gruen zones with proximal, mid, and distal portions [7, 12].

Four different coating thicknesses were analyzed (100, 200, 300, and 400 μm) and compared to the uncoated condition. These thicknesses were chosen in accordance with the standard values for electrophoretic deposition technique (EPD) [13], a method that allows fair PEEK coatings. The geometry of the hip implant is exactly the same for all simulations. The thickness portion of the coating is taken from the cancellous bone with the purpose of mimicking, the THA procedure (see Figure 1).

2.2. Geometric Modeling and Meshing

The standard “Sawbones” Pacific Research Labs Inc. model was used as starting point for the femur [14]. That model modified to create a new geometry of the femur after the THA procedure, on the PTC Creo© software. A series of parallel sketches were used in order to create a solid geometry with the least amount of sharp edges and therefore facilitate the discretization process. In addition, computerized models of the hip implant and the coating were developed.

Implant model was based on a typical metal-on-metal design. In order to facilitate the analysis, the implant geometry is a simplified version of a commercial implant such as Biomet Bi-Metric®.

Figure 2 shows the complete assembled geometry of the implant, femur, and coating. The coating is assumed to be uniform and perfectly bonded to the stem of the hip implant.

In order to perform structural analysis by Finite Element Method (FEM), the assembly geometry was imported into COMSOL Multiphysics© software. The computer model was discretized with different global and local sizes. Due to the complexity of the geometry, tetrahederal elements were necessary. Finally, in order to assure numerical convergence, the most suitable mesh for the study consisted of 1,090,944 tetrahedral elements, with 0.1 and 15.4 mm as the minimum and maximum element size, respectively.

2.3. Material Properties and Boundary Conditions

Four different types materials were considered during numerical analysis: cancellous bone, cortical bone, PEEK coating, and hip implant. All materials are considered to be isotropic based on the part properties derived from common manufacturing techniques such as casting [15] and EPD process for the coating [16]. The hip implant material is of medical grade Ti-6Al-4V.

Mechanical properties of biological tissues vary according to factors such as age, gender, race, and other factors. Average values for cancellous and cortical bone are used [1719]. The properties for the PEEK 150 XF (polymer used for coating) were obtained from the Victrex data sheet ( Table 1 shows a summary of the mechanical properties used in the numerical analysis.

MaterialElastic modulus, [GPa]Density, [g/cm3]Poisson’s ratio,

Cancellous bone 0.1550.200.30
Cortical bone 16.701.640.30
PEEK 150 XF 3.701.300.40

The femur was rigidly fixed at the distal end. The assembly is set to form a union, in order to simulate complete interdigitating of cortical bone, cancellous bone, and PEEK coating. A vertical load of 3,000 N is applied to the femoral implant head, representing 4 times the body weight of a 75 kg patient [17].

3. Results

The results are given in terms of effective von Mises stresses as experienced by the cancellous bone around the hip implant, according to suggested criteria used in previous studies [20]. The three regions mentioned before (proximal, mid, and tip) are used to analyze the change in stress transmitted to the femur. Each of these zones is evaluated independently, comparing a mean value of von Mises stress in the coated with uncoated condition.

3.1. Surface Analysis

COMSOL Multiphysics software allows the analysis of von Mises stresses by selecting surfaces and associating nodes on the mesh. The internal surfaces of the femur in contact with the hip implant were divided into proximal, mid, and tip zones. These zones required a different number of nodes: proximal (6,811 nodes), mid (6,291 nodes), and tip (3,277 nodes).

Figure 3 shows the Finite Element analysis results for effective von Mises stresses in the proximal, mid, and distal zones of the cancellous bone. The comparison of uncoated versus coated conditions (400 μm) shows a significant increase in the transmitted stress to the femur. For the coated condition, stress is almost uniformly distributed in all lateral surfaces, except for a peak zone with a stress value of 10.79 MPa.

In the mid zone (see Figure 3(b)), the difference in stress distribution is more significant compared to the proximal and tip zones. In the case of the tip zone analysis the increase in stress distribution is moderate.

Figure 4 shows the stress values for a curve along each of the zones proposed in this study. Each point represents a node along the curve of analysis. The starting point with length 0 mm is always at the top of the selected zone.

Figure 4(a) shows the effective von Mises stress on a single curve along the proximal zone. As the length increases, the effect of the coating is more noticeable. It is clear that as the coating thickens higher loads are transferred to the femur and, consequently, less stress shielding will be presented. The difference in stress transmission for a coating thickness of 400 μm versus 100 μm is only significant at a length between 35 and 50 mm. Therefore, for the purpose of producing a workable electrophoresis coating, the thickness can be maintained below 400 μm.

Figures 4(b) and 4(c) present the results for mid and tip zone, respectively. The mid zone shows the best results in terms of load transferred to cancellous bone. In Figure 4(b) the behavior of stress distribution changes across the curve length, comparing the uncoated and the coated condition. After 60 mm all lines for each coating thickness almost overlap. For the tip zone, Figure 4(c), values of effective von Mises are higher here than in any other zone of the bone. In this case, the curve of analysis goes around the tip of the implant (see Figure 4(c)).

For comparison, Figure 5 shows the analysis based on maximum principal stress. Similar to the analysis based on effective von Mises stress, the coating shows better transmission of the load to the cancellous bone.

Figure 6 shows the complete surface analysis for the cancellous bone. The results show a clear increase of stress in the lateral faces of the bone with coated condition (red zones). Under the PEEK coating of 400 μm, the highest stress at the cancellous bone is 5.5 MPa.

3.2. Volumetric Analysis of Cancellous Bone

It is important to analyze each zone individually since proximal and mid zones are where higher stress shielding values are reported in the research literature. The current study reports consistent results, where the proximal and mid zones have a more significant role in terms of load transferred to bone, even with the 100 μm coating.

The analysis shown in Figure 7 and Table 2 is based on volumetric average of the effective von Mises stress. It is shown that the gain in load transfer to the bone, for the complete volume of cancellous bone, is 81% for a coating of 100 μm and 92% for a coating of 400 μm. When separating the three zones, the main benefit appears in the proximal and mid zones.

Coating thickness [m]Proximal zone cancellous bone volumeMid zone cancellous bone volumeTip zone cancellous bone volumeComplete cancellous bone volume
Average stress [MPa]Change coated versus uncoated [%]Average stress [MPa]Change coated versus uncoated [%]Average stress [MPa]Change coated versus uncoated [%]Average stress [MPa]Change coated versus uncoated [%]

1002.6947 2.061902.6916%2.3981%
2002.7651 2.091942.6916%2.4384%
3002.8455 2.132002.7217%2.4888%
4002.9460 2.162042.7217%2.5492%

Note: effective von Misses stress.
3.3. Volumetric Analysis of Femur

Figure 8 shows the results of the Finite Element analysis for the complete model, including the hip implant and the whole femur during a load application of 3,000 N. As expected, the highest stress concentration is on the neck of the hip implant, reaching 180 MPa. It is important to note that the hip implant is hypothetically made of a medical grade titanium alloy that allows these levels of stress. The stress level for the femur is within acceptable levels.

4. Discussion

4.1. Model Limitations

Some considerations have been taken into account in the model in order to optimize computational resources without affecting analysis fundamentals. The model geometry, specifically on the proximal part of the femur, has been simplified and does not exactly correspond to a femur after THA. Given that the proposed design concept is based on press-fit type implants, a porous surface should be considered. Again, in order to reduce computational load, all analysis were conducted assuming smooth surfaces.

The present study is limited by the use of isotropic mechanical properties for bone, in order to assess the viability of the proposed design concept for the implant. A future more detailed study should consider the actual anisotropic mechanical properties of bone, as well as the variations in bone density due to the remodeling process.

Factors such as the muscle surrounding the femur and changes in force angle while walking are not taken into account in this study. For the purpose of knowing the effect of a coating in the load transmission to the bone, these factors are not essential since human femur while walking is mainly subjected to compression from axial loading. Previous studies have also simplified their models in such a fashion [2123].

4.2. Influence of PEEK Coating and Implant Design on Stress Shielding

Proximal zone presents the higher stress shielding values. During the THA procedure, in order to insert the implant into the femur, cancellous and the intramedullary canal are removed. Then anchoring cement is used to attach the implant to the femur. Since the femur is not a uniform and symmetric geometry, the THA procedure can leave zones with thinner cancellous bone, causing stress peaks as the one seen in the proximal zone in this analysis.

The proposed approach shows an increase of 47% in load transfer in the proximal zone with a PEEK coating with 100 μm thickness, while 400 μm thickness improves by 60% the same indicator. This results compare favorably with the work of Gross and Abel [6], where the increment of von Misses stress in the proximal zone increases 32% with the use of hollow stemmed hip implant.

In a comprehensive review of bone tissue fracture, Doblaré et al. [20] suggest that using the effective von Mises stress is appropriate when considering isotropic bone properties. However, in some studies of this nature, brittle material failure criteria are used: maximum principal stress. In terms of principal maximum stress (compressive stress), the current study was compared to the work of Oshkour et al. [3], where axial load is also taken at 3,000 N. The range of maximum principal stress (compressive stress) found is approximately between 0.3 and 5 MPa (see Figure 5), while Oshkour et al. report a similar range between 0.2 and 4 MPa.

The scope of the present study is to validate the design concept of using a PEEK coating to increases load transmission to the bone. Therefore, most of the analysis conducted here was based on effective von Misses stress as first approach.

In order to have a more robust perspective, further studies should be carried out in order to include muscles such as iliotibial tract. In addition, variation of the angle of the femur while walking should be considered in a more in-depth numerical analysis.

The increased stress at the tip of the implant with PEEK coating, compared to uncoated, might be a concern for periprosthetic fracture. Therefore, a more detailed analysis is required in order to consider an implant redesign of the implant tip or elimination of the PEEK coating in that zone.

The proposed approach opens the possibility to improve mechanical behavior without the need to develop a new implant geometry. With an appropriate coating process, any commercially available implant could be treated in order to improve its performance. The present method of coating a hip implant is not limited to a single implant shape. The relative benefit on the proposed approach on stress shielding will certainly depend on the specific implant geometry, given the wide range of commercially available designs.

The study shows that a significant improvement in the implant performance is reached with a modest amount of polymer coating. In this regard, there are several methods available for coating the hip implant with a polymer. The strongest possibility is to use electrophoretic deposition due to scalability and coating consistency. Electrophoretic deposition is proven method for PEEK deposition in steel [24]; thus other metallic substrates could work in a similar way. One area of concern is the adhesion strength of such a coating over a medical grade titanium stem.

4.3. Potential Osseointegration of PEEK Coating

Mineral bone growth, and thus fixation of implant to bone, highly depends on surface roughness and porosity. The use of a PEEK coating on the implant shows the potential to reduce stress shielding, but the question of osseointegration arises. Coatings based on porous hydroxyapatite have been used to promote osseointegration. On the other hand, a number of manufacturing processes are available to generate porous coatings that combine polymers and ceramics. Therefore, this particular issue requires further in vitro and in vivo studies.

5. Conclusions and Future Work

The current study shows that a PEEK coating on a hip implant can improve the load transfer to the bone, minimizing the stress shielding effect and thus prolonging the implant lifespan. The models with coated condition show a significant change in cancellous bone stress, with increase between 81 and 92% in the volumetric numerical analysis.

Further research is required to refine the proposed design concept and develop an appropriate coating process. Additional numerical analysis is required to consider a wider set of loading conditions, as well as anisotropic bone properties. In terms of manufacturing, different coating processes and PEEK formulations require investigation for the hip implant application.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This research work was possible due to the support by Tecnológico de Monterrey through its Centro de Innovación en Diseño y Tecnología and research group in advanced manufacturing.


  1. U. Holzwarth and G. Cotogno, JCR Scientific and Policy Reports—Total Hip Arthroplasty, European Commission, 2012.
  2. S. Kurtz, K. Ong, E. Lau, F. Mowat, and M. Halpern, “Projections of primary and revision hip and knee arthroplasty in the United States from 2005 to 2030,” The Journal of Bone & Joint Surgery—American Volume, vol. 89, no. 4, pp. 780–785, 2007. View at: Publisher Site | Google Scholar
  3. A. A. Oshkour, N. A. A. Osman, M. Bayat, R. Afshar, and F. Berto, “Three-dimensional finite element analyses of functionally graded femoral prostheses with different geometrical configurations,” Materials and Design, vol. 56, pp. 998–1008, 2014. View at: Publisher Site | Google Scholar
  4. R. Huiskes, H. Weinans, H. J. Grootenboer, M. Dalstra, B. Fudala, and T. J. Slooff, “Adaptive bone-remodeling theory applied to prosthetic-design analysis,” Journal of Biomechanics, vol. 20, no. 11-12, pp. 1135–1150, 1987. View at: Publisher Site | Google Scholar
  5. M. G. Joshi, S. G. Advani, F. Miller, and M. H. Santare, “Analysis of a femoral hip prosthesis designed to reduce stress shielding,” Journal of Biomechanics, vol. 33, no. 12, pp. 1655–1662, 2000. View at: Publisher Site | Google Scholar
  6. S. Gross and E. W. Abel, “A finite element analysis of hollow stemmed hip prostheses as a means of reducing stress shielding of the femur,” Journal of Biomechanics, vol. 34, no. 8, pp. 995–1003, 2001. View at: Publisher Site | Google Scholar
  7. C. Boyle and I. Y. Kim, “Comparison of different hip prosthesis shapes considering micro-level bone remodeling and stress-shielding criteria using three-dimensional design space topology optimization,” Journal of Biomechanics, vol. 44, no. 9, pp. 1722–1728, 2011. View at: Publisher Site | Google Scholar
  8. Y. Hirata, Y. Inaba, N. Kobayashi, H. Ike, H. Fujimaki, and T. Saito, “Comparison of mechanical stress and change in bone mineral density between two types of femoral implant using finite element analysis,” Journal of Arthroplasty, vol. 28, no. 10, pp. 1731–1735, 2013. View at: Publisher Site | Google Scholar
  9. H. Bougherara, M. Bureau, M. Campbell, A. Vadean, and L. Yahia, “Design of a biomimetic polymer-composite hip prosthesis,” Journal of Biomedical Materials Research, Part A, vol. 82, no. 1, pp. 27–40, 2007. View at: Publisher Site | Google Scholar
  10. P. Tavakkoli Avval, S. Samiezadeh, V. Klika, and H. Bougherara, “Investigating stress shielding spanned by biomimetic polymer-composite vs. metallic hip stem: a computational study using mechano-biochemical model,” Journal of the Mechanical Behavior of Biomedical Materials, vol. 41, pp. 56–67, 2015. View at: Publisher Site | Google Scholar
  11. S. M. Kurtz and J. N. Devine, “PEEK biomaterials in trauma, orthopedic, and spinal implants,” Biomaterials, vol. 28, no. 32, pp. 4845–4869, 2007. View at: Publisher Site | Google Scholar
  12. T. A. Gruen, G. M. McNeice, and H. C. Amstutz, “'Modes of failure' of cemented stem-type femoral components. A radiographic analysis of loosening,” Clinical Orthopaedics and Related Research, vol. 141, pp. 17–27, 1979. View at: Google Scholar
  13. L. Besra and M. Liu, “A review on fundamentals and applications of electrophoretic deposition (EPD),” Progress in Materials Science, vol. 52, no. 1, pp. 1–61, 2007. View at: Publisher Site | Google Scholar
  14. P. B. Chang, B. J. Williams, K. S. B. Bhalla et al., “Design and analysis of robust total joint replacements: finite element model experiments with environmental variables,” Journal of Biomechanical Engineering, vol. 123, no. 3, pp. 239–246, 2001. View at: Publisher Site | Google Scholar
  15. J. V. Abellán-Nebot, H. R. Siller, C. Vila, and C. A. Rodríguez, “An experimental study of process variables in turning operations of Ti-6Al-4V and Cr–Co spherical prostheses,” International Journal of Advanced Manufacturing Technology, vol. 63, no. 9–12, pp. 887–902, 2012. View at: Publisher Site | Google Scholar
  16. X. F. Xiao and R. F. Liu, “Effect of suspension stability on electrophoretic deposition of hydroxyapatite coatings,” Materials Letters, vol. 60, no. 21-22, pp. 2627–2632, 2006. View at: Publisher Site | Google Scholar
  17. H. Ebrahimi, M. Rabinovich, V. Vuleta et al., “Biomechanical properties of an intact, injured, repaired, and healed femur: an experimental and computational study,” Journal of the Mechanical Behavior of Biomedical Materials, vol. 16, no. 1, pp. 121–135, 2012. View at: Publisher Site | Google Scholar
  18. W. X. Niu, L. J. Wang, T. N. Feng, C. H. Jiang, and Y. B. Fan, “Effects of bone Young's modulus on finite element analysis in the lateral ankle biomechanics,” Applied Bionics and Biomechanics, vol. 10, no. 4, pp. 189–195, 2013. View at: Publisher Site | Google Scholar
  19. L. Baggi, M. Di Girolamo, G. Vairo, and G. Sannino, “Comparative evaluation of osseointegrated dental implants based on platform-switching concept: influence of diameter, length, thread shape, and in-bone positioning depth on stress-based performance,” Computational and Mathematical Methods in Medicine, vol. 2013, Article ID 250929, 15 pages, 2013. View at: Publisher Site | Google Scholar
  20. M. Doblaré, J. M. García, and M. J. Gómez, “Modelling bone tissue fracture and healing: a review,” Engineering Fracture Mechanics, vol. 71, no. 13-14, pp. 1809–1840, 2004. View at: Publisher Site | Google Scholar
  21. S. Shah, S. Y. R. Kim, A. Dubov, E. H. Schemitsch, H. Bougherara, and R. Zdero, “The biomechanics of plate fixation of periprosthetic femoral fractures near the tip of a total hip implant: cables, screws, or both?” Proceedings of the Institution of Mechanical Engineers. Part H: Journal of Engineering in Medicine, vol. 225, no. 9, pp. 845–856, 2011. View at: Publisher Site | Google Scholar
  22. E. T. Davis, M. Olsen, R. Zdero, J. P. Waddell, and E. H. Schemitsch, “Femoral neck fracture following hip resurfacing: the effect of alignment of the femoral component,” The Journal of Bone & Joint Surgery—British Volume, vol. 90, no. 11, pp. 1522–1527, 2008. View at: Publisher Site | Google Scholar
  23. D. Wilson, H. Frei, B. A. Masri, T. R. Oxland, and C. P. Duncan, “A biomechanical study comparing cortical onlay allograft struts and plates in the treatment of periprosthetic femoral fractures,” Clinical Biomechanics, vol. 20, no. 1, pp. 70–76, 2005. View at: Publisher Site | Google Scholar
  24. I. Corni, N. Neumann, S. Novak et al., “Electrophoretic deposition of PEEK-nano alumina composite coatings on stainless steel,” Surface and Coatings Technology, vol. 203, no. 10-11, pp. 1349–1359, 2009. View at: Publisher Site | Google Scholar

Copyright © 2016 Jesica Anguiano-Sanchez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.