Next Article in Journal
Facile Synthesis of Antimicrobial Aloe Vera-“Smart” Triiodide-PVP Biomaterials
Next Article in Special Issue
A Review on Development of Bio-Inspired Implants Using 3D Printing
Previous Article in Journal
Fabrication and Characterization of Polylactic Acid Electrospun Scaffolds Modified with Multi-Walled Carbon Nanotubes and Hydroxyapatite Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Element Analysis of Orthopedic Hip Implant with Functionally Graded Bioinspired Lattice Structures

by
Nikolaos Kladovasilakis
,
Konstantinos Tsongas
and
Dimitrios Tzetzis
*
Digital Manufacturing and Materials Characterization Laboratory, School of Science and Technology, International Hellenic University, 14km Thessaloniki, 57001 N. Moudania, Greece
*
Author to whom correspondence should be addressed.
Biomimetics 2020, 5(3), 44; https://doi.org/10.3390/biomimetics5030044
Submission received: 1 July 2020 / Revised: 29 August 2020 / Accepted: 9 September 2020 / Published: 12 September 2020
(This article belongs to the Special Issue Biomimetic Additive Manufacturing)

Abstract

:
The topology optimization (TO) process has the objective to structurally optimize products in various industries, such as in biomechanical engineering. Additive manufacturing facilitates this procedure and enables the utility of advanced structures in order to achieve the optimal product design. Currently, orthopedic implants are fabricated from metal or metal alloys with totally solid structure to withstand the applied loads; nevertheless, such a practice reduces the compatibility with human tissues and increases the manufacturing cost as more feedstock material is needed. This article investigates the possibility of applying bioinspired lattice structures (cellular materials) in order to topologically optimize an orthopedic hip implant, made of Inconel 718 superalloy. Lattice structures enable topology optimization of an object by reducing its weight and increasing its porosity without compromising its mechanical behavior. Specifically, three different bioinspired advanced lattice structures were investigated through finite element analysis (FEA) under in vivo loading. Furthermore, the regions with lattice structure were optimized through functional gradation of the cellular material. Results have shown that optimal design of hip implant geometry, in terms of stress behavior, was achieved through functionally graded lattice structures and the hip implant is capable of withstanding up to two times the in vivo loads, suggesting that this design is a suitable and effective replacement for a solid implant.

1. Introduction

The additive manufacturing (AM) process has enabled the production of complex geometries [1,2,3,4] along with composite material structures [5,6], which are difficult to fabricate with traditional manufacturing techniques, such as machining and injection molding. This led to rapid development of the topology optimization process with advanced geometries via generative design and cellular materials. Topology optimization (TO) is a procedure that optimizes the material mass distribution within the already defined external volume [7], resulting in light-weight structures and minimization of feedstock material usage. The main objective is to maintain the desired mechanical properties reducing the mass of the structure; however, it also creates additional advantages such as high porosity and high surface-area-to-volume ratio [8]. There are two different approaches to achieve TO: the density-based and truss-based approaches. The density-based approach is identical to generative design [9]. On the other hand, the truss-based approach utilizes periodic unit cells (lattice structures) in order to achieve the optimum mass distribution [10]. Nowadays, there is a plethora of applications in several industries with topologically optimized products using both of the two aforementioned approaches, indicatively in the aeronautical, automotive, and biomechanical industries [11,12,13].
The truss-based approach imitates the cell structures from natural tissues, such as bones, corals, foams, and so forth. Thus, this approach is suitable for topology optimization of structures in biomechanical applications, such as implants and tissue scaffolds. Specifically, in this approach, the selected volume fills with lattice structures with specific geometry and dimensions replacing the solid material. Besides the advantage of minimizing the mass of the product for usage in biomechanical applications, this procedure offers the above-mentioned additional advantages, which facilitate the tissue regeneration process and allow diffusion of oxygen and nutrients [14,15]. There is a plethora of different lattice structures, the simplest ones being the 2D lattice structures like honeycombs and prismatic ones [16]. The geometrical complexity increases in 3D lattice structures, resulting in a vast number of diverse 3D lattice structures with further classification in strut structures and sheet structures [17]. The most widespread strut structures are the Octet and the Voronoi, which directly emerge from structures found in nature. On the other hand, the most common lattice sheet structures are the sheet triply periodic minimal surface (TPMS) structures such as Gyroid, Schwarz Diamond, and Neovius. According to Gibson et al. [18], all lattice structures lead to reduction of the mechanical properties. The reduction of the mechanical properties in lattice structures depends mainly on the relative density and on the actual geometry of the applied lattice structure. Relative densities less than 50% enhance the size effect, which has an increasing influence on the mechanical performance as the relative density is reduced. However, not all lattice structures have the same mechanical behavior; specifically, lattice structures that exhibit stretching-dominated behavior are less affected by the size effect compared to lattice structures that exhibit bending-dominated behavior [19]. Moreover, Al Ketan et al. [20] propose lattice structures that have stretching-dominated behavior for biomechanical applications which receive increased stresses.
In addition, there are several studies that have tried to combine additive manufacturing techniques with orthopedic implants through customization and lattice structures. Mahmoud et al. [21] have gathered the majority of these studies in a comprehensive review, summarizing the use of AM technologies to produce orthopedic implants from lattice structures and functionally graded materials. Gabbrielli et al. [22] and España et al. [23] investigated the possibility of using lattice structures in acetabular cups for mechanical and biological advantages. Moreover, Hazlehurst et al. [24] suggested manufacturing a hip implant consisting of cubic lattice structures in the whole internal body of the implant and this led to severe reduction of the implant’s stiffness and strength. Furthermore, Limmahakhun et al. [25] and González et al. [26] investigated innovative designs in order to optimize the utility of lattice in orthopedic implants and increase their mechanical strength and this was achieved by functional gradation of the lattice structures.
The aim of this research was to extract an innovative design for a bioinspired hip implant utilizing topology optimization tools. Therefore, the current paper proposes a novel hip implant design with advanced lattice structures in order to address the reduction in strength and handle the in vivo loadings. Figure 1 portrays a flowchart of the topology optimization process through lattice structures depending on the mechanical behavior. Initially, a hip implant was designed according to international medical standards and then its mechanical behavior was examined under in vivo static loads through finite element analysis (FEA). Furthermore, specific regions of the implant, where topology optimization was essential, were replaced by bioinspired lattice structures that have been shown to exhibit stretching-dominated behavior, such as the Voronoi strut structure and the Gyroid and Schwarz Diamond sheet structures. These TPMS structures were selected due to the promising mechanical strength that was exhibited in the aforementioned studies [19,20]. In addition, Voronoi structure was selected because it displays similar behavior to the trabecular internal structure, while demonstrating high mechanical strength [27,28]. These three lattice structures were examined and evaluated for their mechanical performance under the same in vivo loads. In addition, further design optimization was applied to lattice structures through functional gradation in order to improve the mechanical properties of a hip implant.

2. Materials and Methods

2.1. Design of Hip Implant

Hip implant design follows a series of international standards in order to fulfill its purpose [29]. A typical orthopedic hip implant consists of three distinct regions: head, neck, and intramedullary stem (rest of the body). The main design considerations are the head diameter, the diameter and length of the neck, the length of the intramedullary implant (stem), the stem cross-section area, and the angle of placement of the head relative to the main body.
Design parameters are configured from past research and existing literature. According to the studies of Charnley et al. [30] and McKnee et al. [31], the head diameter ranges from 22 mm to 45 mm, in order to fit in the hip joint. Moreover, the length and diameter of the neck range from 10 mm (short-neck) to 40 mm (long-neck) and 13 mm to 30 mm, respectively [32,33]. The length of the intramedullary stem of the hip implant ranges from 120 mm to 180 mm and the neck-shaft angle varies between 135° and 145° [34]. Table 1 summarizes the design parameters and their aforementioned ranges.
In order to apply topology optimization through lattice structures on an orthopedic hip implant, a new hip implant was designed fulfilling all the above-mentioned standards. Therefore, the designed hip implant has the following geometry and dimensions. More specifically, the intramedullary stem’s length is 128 mm, while the neck’s length is 35 mm. In addition, the diameter of the neck is 18 mm and the diameter of the head is 45 mm. The placement angle of the neck and head in relation to the axis of the implant is 135° degrees and the angular range is almost 120° degrees. Figure 2 portrays the design of the hip implant, its basic dimensions, and its cross-section, as it was designed in SolidWorks™ software.

2.2. Lattice Structure Configuration

Cellular materials are scattered in naturally formed structures, such as foams, corals, bones, and so forth. These natural materials are the result of the evolution of thousands of years, thus both mechanical properties and functional properties have been optimized depending on each application. There is a variety of lattice structures that differ in shape and geometry. Therefore, for the topology optimization of a human body implant, it is logical to investigate lattice structures derived from nature that could serve the function of the implant.
In the case of hip implant, lattice structures should be able to withstand increased loads relative to the body weight and the human movement. Hence, a potential candidate structure which imitates the inside trabecular structure is the Voronoi truss structure [33,34]. 3D Voronoi structures follow this mathematical formula [35]:
V(pi) = {p/d(p,pi) ≤ d(p,pj), ji, i, j = 1, ..., n}
where p1, …, pn are the distinct seeds in 3D space, d (p,pi) is the Euclidean distance between two seeds (p and pi), and V(pi) represents the 3D Voronoi polygon on seed pi. Figure 3 shows the structure of trabecular bone and an indicative Voronoi/trabecular-like structure.
Other bioinspired lattice structures that could withstand high stress loads and offer the proper conditions to facilitate tissue regeneration are triply periodic minimal surfaces (TPMSs). TPMSs are characterized by surfaces that have an average curvature value equal to zero. TPMS-like geometries result from applying Fourier transformations on level-set equations and by using a level-set parameter c in Equations (2) and (3) [15]. According to the literature [16], the most advanced TPMS structures, with increased strength, are the Gyroid and Schwarz Diamond structures. A Gyroid-like structure is also observed on butterfly wings (Figure 4a), offering high tolerance in bending loads [38]. A Schwarz Diamond-like structure is used by beetle L. Augustus in order to reinforce its shell structure (Figure 4b) and has also shown remarkable mechanical strength in various experiments [39,40,41]. It is worth noting that in this study, the aforementioned structures have been used for the topology optimization of the hip implant. The design software nTopology™ was also used to design and implement these structures in hip implant geometry. The level-set equations for Gyroid-like and Schwarz Diamond-like geometries are shown in the following equations:
φg = sin x cos y + sin y cos z + sin z cos x = c
φd = sin x sin y sin z + sin x cos y cos z + cos x sin y cos z + cos x cos y cos z = c

2.3. Finite Element Analysis

Finite element analysis (FEA) simulated each hip implant design with different lattice structures under static in vivo conditions. It is worth mentioning that FEA for each implant focused mainly on the determination of yield point and the calculation of safety factor of the topologically optimized implant design. ANSYS™ software was used for the finite element analysis and in particular, the static structural module was used to simulate the quasi-static loading. Generally, mesh generation for lattice structures could be implemented with two different types of elements: hexahedrals [43] and tetrahedrals [44]; in these studies, plenty of references are given towards the use of the different discretization methods. In this research, the mesh consisted of tetrahedral elements for the whole body of the implant. In the last decade, a new material was proposed for medical use and especially for orthopedic implants; this novel material is the nickel-based superalloy Inconel 718. According to the literature [45], Inconel 718 offers enhanced corrosion resistance and advanced mechanical properties. However, in order to be suitable for this application, Inconel should fulfill the international standards of ASTM F90-14 and UNS N07718 [46]. Table 2 and Table 3 present the chemical composition and the mechanical properties of Inconel 718 which were used in order to build a material model combining isotropic elasticity and bilinear isotropic hardening.
Furthermore, in order to complete the FEA model, it is necessary to determine, based on the literature, the conditions for the implant’s support and the maximum loads that the implant receives from the human body under in vivo loading conditions. Hence, according to Ducheyne et al. [47], the hip implant has fixed support at the bottom of its structure in contact with the femur bone, without relying on calcar support. The received loads remain to be defined as to the point of application, the direction of the loading, and the amount of forces. El-Shiekh et al. [28], who have carried out a comprehensive research in hip joint replacement under static and dynamic loading, have shown that the forces caused by the human body’s weight are applied on the center of the hip implant’s head. McLeish et al. [48] concluded on the results that for zero degree of pelvis angle, the direction of the human body’s force is approximately at 20° degrees to the center axis of the femur bone. Finally, the magnitude of the force depends on the weight of the human body and the type of body movement. Bergmann et al. [49] proposed a percentage relationship between the type of movement and the weight of the human body. The percentage weight relations, depending on the type of movement, are listed in Table 4 as well as the nominal values of the force corresponding to an average body weight of 75 kg [50].

3. Results and Discussion

3.1. Finite Element Analysis for Solid Hip Implant

Before the topology optimization process, it was necessary to evaluate the initially designed solid hip implant’s mechanical behavior through finite element analysis. Therefore, the 3D CAD model of the solid hip implant (Figure 2) was examined with FEM under in vivo conditions, as described previously. Figure 5 shows the contours of von Mises stresses (Figure 5a) and the contours of the factor of safety (Figure 5b) for maximum load (i.e., at a force of 5300 N). These analyses extracted a series of conclusions which are essential for the topology optimization process. The first conclusion is that the head (ball) where the loads have been applied has shown low stress concentration, hence it obtains small strains compared to the rest of the body. It was expected that the maximum stress would occur in a region adjacent to the support of the hip implant. This is because the static system works like a cantilever beam and the maximum bending moment will occur at the fixed boundary condition. Therefore, both the neck and the intramedullary stem are divided into three distinct regions (Figure 5) of strain: the tensile region (left side), the low stresses region (center), and the compressive region (right side). In the tensile region, due to the direction of the applied force, tensile stresses were concentrated, however, the amount of stress was low, revealing that this region has moderate risk for fracture. The central region of implant was also the area with the lowest stress concentration, where the improvement of mass distribution is essential. Moreover, the compressive region is the area where, due to the direction of the applied force, there is concentration of compressive stresses and intervention is needed during the process of topology optimization.
As shown in Figure 5, the maximum stresses were compressive and they reached up to 441 MPa. However, these stresses occurred only in one specific point and did not reflect the total mechanical behavior of the implant, as the rest of the stresses were in the range between 90 MPa and 300 MPa. Nonetheless, even the maximum stresses were much lower than the material yield point. This was deducted since the factor of safety (FOS) has shown a minimum value of 2.5. Therefore, it could be concluded that the hip implant with this specific material and design withstands loads 2.5 times more than the maximum load that the implant receives from the human body in tripping movement, thus an extensive topology optimization of the design could be performed.

3.2. Topology Optimization through Lattice Structures

The process of topology optimization of the hip implant’s design has been analyzed using three different lattice structures, namely, the Voronoi strut structure and the Gyroid and the Schwarz Diamond sheet structures, which are bioinspired. The objective was to evaluate the mechanical behavior of each structure and to highlight the most advanced structures in order to proceed to the next phase of optimization, which was the adaptive adjustment of the lattice structure through functional gradation. Topology optimization exploits the findings of the FEA of the solid hip implant in order to define certain areas of the design to be optimized. The next step was to select the appropriate relative density of the proposed lattice structures to be applied in these regions. The appropriate value of the relative density should be defined by the right trade-off between the deterioration of the mechanical properties, due to size effect, and the desired advantages (i.e., reduction of the mass and increased porosity). Therefore, some of the drawbacks of the size effect, like reduced modulus or a loss in strength, could be overcompensated by lightweight structures with enhanced porosity. Thus, the relative density was chosen to be 50%, due to the fact that it has a negligible influence from the size effect and it was close enough to the relative density of the femur bone [51] enhancing the body’s ability to regenerate tissues.
The selection of redesigned regions with lattice structure was a result of FEA in the solid hip implant. Therefore, as illustrated in Figure 5, due to increased loads, both compressive and tensile in the neck of the implant, the whole area remained solid. The loads of the implant were applied on the head; thus, it was necessary to have a continuous surface of 2 mm thickness. Since the implant’s head demonstrated low stress concentration, only its interior was redesigned with lattice structures, as shown in Figure 6. In addition, Figure 6 shows that the whole intramedullary stem was redesigned with lattice structures, except the region where the implant was mounted, where it remained solid. Figure 6a–c portrays two different views, side and front, of the Voronoi, Gyroid, and Schwarz Diamond structures, respectively.
Figure 7 portrays the results of the finite element analysis for three different approaches of the hip implant’s topology optimization, using bioinspired lattice structures such as Voronoi, Gyroid, and Schwarz Diamond structures. Specifically, according to the FEA results of the solid hip implant, the topology optimization approach has shown that the highest strength was observed in the hip implant that contained the Schwarz Diamond structure. The hip implant which contained Schwarz Diamond portrayed the lowest stress concentration with a maximum value of 1171 MPa at the compressive region. In addition, the hip implant with Gyroid structure revealed ultimate compressive stresses at a value of 1405 MPa. The hip implant with Voronoi structure has demonstrated the weakest behavior compared to the other structures under study, with maximum stress concentration at 1738 MPa.
Τhe factor of safety was calculated by the maximum stress concentration of each structure, as listed in Table 5. In addition to the standard FOS, FOSUlt has been also used and it was calculated as the ultimate factor of safety which depends on the ultimate strength of the material, as shown in Table 5. As the FOS values verified, only the Schwarz Diamond structure could withstand the forces of the human body without fracture of the implant, as FOSUlt was greater than the value of 1 (FOSUlt =1.17 > 1). However, the stress concentration in the Gyroid structure was marginally greater than the fracture point of the material, thus it is expected that with an improved design approach through functional gradation, the strength of the implant could be increased. On the other hand, the Voronoi structures should exhibit a fundamental design modification in order to achieve FOS values greater than 1.

3.3. Functional Gradation of Lattice Structures

Finite element analyses of the topologically optimized versions of the hip implant through bioinspired lattice structures showed the mechanical superiority of sheet structures, and especially TPMS, over strut structures, as discussed previously. However, through functional gradation of lattice structures, it is possible to achieve further topology optimization of the part. Thus, knowing the regions where the stresses are concentrated, design could be changed to increase or decrease the local relative density of the selected lattice structures, keeping the total relative density of the object constant but increasing its mechanical performance. Figure 8 illustrates a graphical representation of a functionally graded Gyroid structure with mean relative density of 40%.
This procedure was followed for all above-mentioned structures. In the areas where the highest stresses (compressive region) were concentrated, it was decided to apply a denser pattern of lattice structure to reach a local relative density of 65%, in order to completely eliminate the size effect. On the other hand, in the region with low stresses (i.e., at the center of the implant), a sparser local relative density was decided to be implemented as the stresses in this area were lower, hence the integrity of the object would not be affected by the size effect. Since the aim was to maintain the total relative density constant, these two areas were separated using a try-and-error algorithm, in such a way as to produce redesigned hip implants that have the same mass as the Voronoi, Gyroid, and Schwarz Diamond structures. It is worth mentioning that the topology optimization of the hip implant’s head and the tensile region remained the same as the previous design, as further removal of material may have affected locally the integrity of the object’s mechanical behavior. Figure 9 shows the final geometries of the topologically optimized hip implant with Voronoi structure (Figure 9a), hip implant with Gyroid structure (Figure 9b), and hip implant with Schwarz Diamond structure (Figure 9c). All the hip implant designs with functional gradation of lattice structures were designed and calculated using the nTopology™ platform.
Figure 10 illustrates the results of FEA and the respective contour of stresses for the hip implant with Voronoi structure (Figure 10a), with Gyroid structure (Figure 10b), and with Schwarz Diamond structure (Figure 10c). Figure 10a shows that the topologically optimized hip implant with functional gradation of Voronoi structure concentrated lower stresses than the corresponding design approach with Voronoi without functional gradation. Specifically, this hip implant design has shown a peak of compressive stresses at 1088 MPa.
Furthermore, the topologically optimized hip implant with functional gradation of Gyroid and Schwarz Diamond structures revealed even more lower compressive stresses which resulted in a remarkable improvement of the hip implant’s mechanical performance. As it is shown in Figure 10b,c, these hip implant designs have demonstrated a maximum value of compressive stress at 615 MPa and 529 MPa, respectively.
In order for the redesigned hip implants to handle the forces that they receive from the human body, it is important to access their mechanical behavior through the factor of safety. Hence, the factor of safety for the design approach containing the lattice structure of functionally graded Voronoi was 1.01 and for the design approach containing the TPMS structures of functionally graded Gyroid and Schwarz Diamond were 1.79 and 2.08, respectively. These results showed that the design approaches of the hip implant which consisted of functionally graded lattice structures could withstand the maximum static loads that the hip joint receives under in vivo conditions. Regarding the functionally graded Voronoi structure hip implant, it could marginally withstand these loads without being plastically deformed, as shown in the values of FOS in Table 6. The functionally graded TPMS structures, Gyroid and Schwarz Diamond, could handle almost twice the loads (Table 6).

4. Conclusions

In this paper, an orthopedic hip implant was designed according to international standards and then a topology optimization of its geometry was performed. The topology optimization process was implemented via bioinspired lattice structures, namely, Voronoi, Gyroid, and Schwarz Diamond structures, which are derived from nature having superior mechanical performance. Moreover, topology optimization occurred with the implementation of these lattice structures in regions of low stress, in order to achieve the optimal mass distribution within the existing volume. Furthermore, functional gradation of the implemented lattice structures was performed. In particular, it was observed that the hip implant that contained Schwarz Diamond structures revealed the best mechanical behavior, both in simple topology optimized implant and in functionally graded implant. The factor of safety of the functionally graded Voronoi was 1.01 and for the design approach containing the TPMS structures of functionally graded Gyroid and Schwarz Diamond were 1.79 and 2.08, respectively. Topology optimization leads to a reduction of weight of 38% compared to the solid version of the hip implant, therefore less construction material is needed. Moreover, for the intramedullary stem, it has a 50% mean porosity to facilitate the process of tissue regeneration through diffusion of cells, oxygen, and other nutrients. Future work will focus on the fabrication of these implants, utilizing additive manufacturing methods and the evaluation of their mechanical properties will be measured.

Author Contributions

Conceptualization, N.K., K.T. and D.T.; methodology, N.K. and K.T.; software, N.K., K.T.; formal analysis, N.K., K.T. and D.T.; investigation, N.K., K.T. and D.T.; resources, N.K.; data curation, N.K.; writing—original draft preparation, N.K.; writing—review and editing, K.T. and D.T.; supervision, K.T. and D.T.; project administration, D.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors express their gratitude to nTopology S.A. for providing the nTop™ platform in order to design and optimize advanced lattice structures.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gibson, I.; Rosen, D.W.; Stucker, B. Additive Manufacturing Technologies; Springer: New York, NY, USA, 2010. [Google Scholar]
  2. Frazier, W.E. Digital Manufacturing of Metallic Components: Vision and Roadmap. In Solid Free Form Fabrication Proceedings; University of Texas at Austin: Austin, TX, USA, 2010. [Google Scholar]
  3. Karavasili, C.; Tsongas, K.; Andreadis, I.I.; Andriotis, E.G.; Papachristou, E.T.; Papi, R.M.; Tzetzis, D.; Fatouros, D.G. Physico-mechanical and finite element analysis evaluation of 3D printable alginate-methylcellulose inks for wound healing applications. Carbohydr. Polym. 2020, 247, 116666. [Google Scholar] [CrossRef]
  4. Xenikakis, I.; Tzimtzimis, M.; Tsongas, K.; Andreadis, D.; Demiri, E.; Tzetzis, D.; Fatouros, D. Fabrication and finite element analysis of stereolithographic 3D printed microneedles for transdermal delivery of model dyes across human skin in vitro. Eur. J. Pharm. Sci. 2019, 137, 104976. [Google Scholar] [CrossRef]
  5. Mansour, M.; Tsongas, K.; Tzetzis, D. Measurement of the mechanical and dynamic properties of 3D printed polylactic acid reinforced with graphene. Polym.-Plast. Technol. Eng. 2019, 58, 1234–1244. [Google Scholar] [CrossRef]
  6. Mansour, M.; Tsongas, K.; Tzetzis, D.; Antoniadis, A. Mechanical and dynamic behavior of fused filament fabrication 3D printed polyethylene terephthalate glycol reinforced with carbon fibers. Polym.-Plast. Technol. Eng. 2018, 57, 1715–1725. [Google Scholar] [CrossRef]
  7. Bendsoe, M.P.; Sigmund, O. Topology Optimization: Theory, Methods, and Applications; Springer: Berlin, Germany; New York, NY, USA, 2003. [Google Scholar]
  8. Ehrgott, M. Multicriteria Optimization; Springer: Berlin, Germany, 2005. [Google Scholar]
  9. Rosen, D.W. Research supporting principles for design for additive manufacturing. J. Virtual Phys. Prototyp. 2014, 9, 225–232. [Google Scholar] [CrossRef]
  10. Larsen, S.D.; Sigmund, O.; Groen, J.P. Optimal truss and frame design from projected homogenization-based topology optimization. J. Struct. Multidiscip. Optim. 2018, 57, 1461–1474. [Google Scholar] [CrossRef] [Green Version]
  11. Gebisa, A.W.; Lemu, H.G. A case study on topology optimized design for additive Manufacturing. IOP Conf. Ser. Mater. Sci. Eng. 2017, 276, 12026. [Google Scholar] [CrossRef] [Green Version]
  12. Li, C.; Kim, I.Y.; Jeswiet, J. Conceptual and detailed design of an automotive engine cradle by using topology, shape, and size optimization. J. Struct. Multidiscip. Optim. 2015, 51, 547–564. [Google Scholar] [CrossRef]
  13. González-Henríquez, C.M.; Sarabia-Vallejos, M.A.; Rodriguez-Hernandez, J. Polymers for additive manufacturing and 4D-printing: Materials, methodologies, and biomedical applications. J. Polym. Sci. 2019, 94, 57–116. [Google Scholar] [CrossRef]
  14. Heinl, P.; Müller, L.; Körner, C.; Singer, R.F.; Müller, F.A. Cellular Ti–6Al–4V structures with interconnected macro porosity for bone implants fabricated by selective electron beam melting. J. Acta Biomater. 2008, 4, 1536–1544. [Google Scholar] [CrossRef] [PubMed]
  15. Dias, M.R.; Guedes, J.M.; Flanagan, C.L.; Hollister, S.J.; Fernandes, P.R. Optimization of scaffold design for bone tissue engineering: A computational and experimental study. J. Med. Eng. Phys. 2014, 36, 448–457. [Google Scholar] [CrossRef] [PubMed]
  16. Mansour, M.T.; Tsongas, K.; Tzetzis, D.; Antoniadis, A. The in-plane compression performance of hierarchical honeycomb additive manufactured structures. IOP Conf. Ser. Mater. Sci. Eng. 2019, 564, 012015. [Google Scholar] [CrossRef]
  17. Helou, M.; Kara, S. Design, analysis and manufacturing of lattice structures: An overview. Int. J. Comput. Integr. Manuf. 2018, 31, 243–261. [Google Scholar] [CrossRef]
  18. Gibson, L.J.; Ashby, M.F. Cellular Solids. Structure and Properties; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  19. Al-Ketan, O.; Al-Rub, R.K.A. Multifunctional mechanical-metamaterials based on triply periodic minimal surface lattices: A review. J. Adv. Eng. Mater. 2019, 21, 1900524. [Google Scholar] [CrossRef]
  20. Al-Ketan, O.; Rowshan, R.; Al-Rub, R.K.A. Topology-Mechanical Property Relationship of 3D Printed Strut, Skeletal, and Sheet Based Periodic Metallic Cellular Materials. J. Addit. Manuf. 2018, 19, 167–183. [Google Scholar] [CrossRef]
  21. Mahmoud, D.; Elbestawi, M.A. Lattice Structures and Functionally Graded Materials Applications in Additive Manufacturing of Orthopaedic Implants: A Review. J. Manuf. Mater. Process. 2017, 1, 13. [Google Scholar]
  22. Gabbrielli, R.; Turner, I.G.; Bowen, C.R.; Wang, H.; Johnston, S.; Rosen, D.; Cheng, A.; Humayun, A.; Cohen, D.J.; Boyan, B.D. Design of a graded cellular structure for an acetabular hip replacement component. In Proceedings of the Seventeenth Solid Freeform Fabrication Symposium, Austin, TX, USA, 14–16 August 2006. [Google Scholar]
  23. España, F.A.; Balla, V.K.; Bose, S.; Bandyopadhyay, A. Design and fabrication of CoCrMo alloy based novel structures for load bearing implants using laser engineered net shaping. J. Mater. Sci. Eng. C 2010, 30, 50–57. [Google Scholar] [CrossRef]
  24. Hazlehurst, K.B.; Wang, C.J.; Stanford, M. The potential application of a Cobalt Chrome Molybdenum femoral stem with functionally graded orthotropic structures manufactured using Laser Melting technologies. J. Med. Hypotheses 2013, 81, 1096–1099. [Google Scholar] [CrossRef]
  25. Limmahakhun, S.; Oloyede, A.; Chantarapanich, N. Alternative designs of load—Sharing cobalt chromium graded femoral stems. J. Mater. Today Commun. 2017, 12, 1–10. [Google Scholar] [CrossRef]
  26. Quevedo González, F.J. Computational Design of Functionally Graded Hip Implants by Means of Additively Manufactured Porous Materials. Ph.D. Thesis, École de Technologie Supérieure, Montreal, QC, Canada, 2016. [Google Scholar]
  27. Fantini, M.; Curto, M.; De Crescenzio, F. A method to design biomimetic scaffolds for bone tissue engineering based on Voronoi lattices. J. Virtual Phys. Prototyp. 2016, 11, 77–90. [Google Scholar] [CrossRef]
  28. Liang, H.; Yang, Y.; Xie, D.; Li, L.; Mao, N.; Wang, C.; Tian, Z.; Jiang, Q.; Shen, L. Trabecular-like Ti-6Al-4V scaffolds for orthopedic: Fabrication by selective laser melting and in vitro biocompatibility. J. Mater. Sci. Technol. 2019, 35, 1284–1297. [Google Scholar] [CrossRef]
  29. Standard Specification for Total Hip Joint Prosthesis and Hip Endoprosthesis Bearing Surfaces Made of Metallic, Ceramic, and Polymeric Materials; ASTM F2033-12; ASTM International: West Conshohocken, PA, USA, 2012.
  30. Charnley, J.; Kamangar, A.; Longfield, M.D. The optimum size of prosthetic heads in relation to wear of plastic sockets in total replacement of the hip. J. Med. Biol. Eng. 1969, 7, 31–39. [Google Scholar] [CrossRef]
  31. McKee, G.K. Development in total hip joint replacement. Lond. Inst. Mech. Eng. 1967, 85–89. [Google Scholar]
  32. Müller, M.E. Total hip prosthesis. J. Clin. Orthop. Relat. Res. 1970, 11, 46–68. [Google Scholar]
  33. Hybbinette, C.H. Long-term Results of Wear of Plastic Hip Prostheses. Arch. Orthop. Trauma. Surg. 1985, 104, 28–30. [Google Scholar] [CrossRef] [PubMed]
  34. El-Shiekh, H.; El-Din, F. Finite Element Simulation of Hip Joint Replacement under Static and Dynamic Loading. Ph.D. Thesis, School of Mechanical and Manufacturing Engineering, Dublin City University, Dublin, Ireland, 2002. [Google Scholar]
  35. Okabe, A.; Boots, B.; Sugihara, K. Spatial Tessellations–Concepts and Applications of Voronoi Diagrams; Wiley: Chichester, UK, 1992. [Google Scholar]
  36. Scheinpflug, J.; Pfeiffenberger, M.; Damerau, A.; Schwarz, F.; Textor, M.; Lang, A.; Schulze, F. Journey into Bone Models: A Review. Genes 2018, 9, 247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Sieberath, A.; Della Bella, E.; Ferreira, A.M.; Gentile, P.; Eglin, D.; Dalgarno, K. A Comparison of Osteoblast and Osteoclast In Vitro Co-Culture Models and Their Translation for Preclinical Drug Testing Applications. Int. J. Mol. Sci. 2020, 21, 912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Pelanconi, M.; Ortona, A. Nature-Inspired, Ultra-Lightweight Structures with Gyroid Cores Produced by Additive Manufacturing and Reinforced by Unidirectional Carbon Fiber Ribs. J. Mater. 2019, 12, 4134. [Google Scholar] [CrossRef] [Green Version]
  39. Pulsifer, D.P.; Lakhtakia, A. Background and survey of bio replication techniques. Bioinspiration Biomim. 2011, 6, 031001. [Google Scholar] [CrossRef]
  40. Galusha, J.W.; Richey, L.R.; Gardner, J.S.; Cha, J.N.; Bartl, M.H. Discovery of a diamond-based photonic crystal structure in beetle scales. Phys. Rev. E 2008, 77, 2–5. [Google Scholar] [CrossRef]
  41. Al-Ketan, O.; Rezgui, R.; Rowshan, R.; Du, H.; Fang, N.X.; Abu Al-Rub, R.K. Micro-architected Stretching-Dominated Mechanical Metamaterials with Minimal Surface Topologies. J. Adv. Eng. Mater. 2018, 20, 1800029. [Google Scholar] [CrossRef]
  42. Vector Clip Art of Black and White Butterfly with Wide Spread Wings. Available online: https://freesvg.org/vector-clip-art-of-black-and-white-butterfly-with-wide-spread-wings (accessed on 10 September 2020).
  43. Alwattar, T.A.; Mian, A. Developing an Equivalent Solid Material Model for BCC Lattice Cell Structures Involving Vertical and Horizontal Struts. J. Compos. Sci. 2020, 4, 74. [Google Scholar] [CrossRef]
  44. Yang, L.; Mertens, R.; Ferrucci, M.; Yan, C.; Shi, Y.; Yang, S. Continuous graded Gyroid cellular structures fabricated by selective laser melting: Design, manufacturing and mechanical properties. Mater. Des. 2019, 162, 394–404. [Google Scholar] [CrossRef]
  45. Khan, M.A.; Prasad, N.R.; Krishnan, S.N.; Raja, S.K.; Jappes, J.T.W.; Duraiselvam, M. Laser-treated austenitic steel and Nickel alloy for human implants. J. Mater. Manuf. Process. 2017, 32, 1635–1641. [Google Scholar] [CrossRef]
  46. Standard Specification for Wrought Cobalt-20Chromium-15Tungsten-10Nickel Alloy for Surgical Implant Applications (UNS R30605); ASTM F90-14; ASTM International: West Conshohocken, PA, USA, 2014.
  47. Ducheyne, P.; Hastings, G.W. Functional behavior of orthopedic biomaterials, Volume II: Fundamentals. CRC Ser. Struct. Prop. Relatsh. Biomater. 1984, 2, 21. [Google Scholar]
  48. McLeish, R.D.; Charnley, J. Abduction forces in the one-legged stance. J. Biomech. 1970, 3, 191–209. [Google Scholar] [CrossRef]
  49. Bergmann, G.; Graichen, F.; Rohlmann, A. Hip joint loading during walking and running, measured in two patients. J. Biomech. 1993, 26, 969–990. [Google Scholar] [CrossRef]
  50. Colic, K.; Sedmak, A.; Grbovic, A.; Tatic, U.; Sedmak, S.; Djordjevic, B. Finite element modeling of hip implant static loading. Procedia Eng. 2016, 149, 257–262. [Google Scholar] [CrossRef] [Green Version]
  51. Thomas, C.D.L.; Feik, S.A.; Clement, J.G. Regional variation of intracortical porosity in the midshaft of the human femur: Age and sex differences. J. Anat. 2005, 206, 115–125. [Google Scholar] [CrossRef]
Figure 1. Flowchart of topological optimization procedure.
Figure 1. Flowchart of topological optimization procedure.
Biomimetics 05 00044 g001
Figure 2. Design of hip implant.
Figure 2. Design of hip implant.
Biomimetics 05 00044 g002
Figure 3. Trabecular bone: (a) trabecular bone schematic [36]; (b) trabecular/Voronoi-like structure examined by scanning electron microscope (SEM) [37].
Figure 3. Trabecular bone: (a) trabecular bone schematic [36]; (b) trabecular/Voronoi-like structure examined by scanning electron microscope (SEM) [37].
Biomimetics 05 00044 g003
Figure 4. TPMS-like structures found in nature: (a) Gyroid in butterfly wings: (i)butterfly schematic was retrieved from [42]; (ii-v) SEM micrographs on a cross sectional view of a butterfly wing scale, where (v) was reproduced from [38]; (b) Schwarz Diamond in L. Augustus beetle shell (i) Photograph of the weevil L. Augustus. (ii) Optical micrograph of individual scales attached to the exoskeleton of L. Augustus under white-light illumination. (iii) Cross-sectional SEM image of a single scale. (iv) Detailed cross-sectional SEM image of a region of a scale. (Reproduced from [40], with permission of APS).
Figure 4. TPMS-like structures found in nature: (a) Gyroid in butterfly wings: (i)butterfly schematic was retrieved from [42]; (ii-v) SEM micrographs on a cross sectional view of a butterfly wing scale, where (v) was reproduced from [38]; (b) Schwarz Diamond in L. Augustus beetle shell (i) Photograph of the weevil L. Augustus. (ii) Optical micrograph of individual scales attached to the exoskeleton of L. Augustus under white-light illumination. (iii) Cross-sectional SEM image of a single scale. (iv) Detailed cross-sectional SEM image of a region of a scale. (Reproduced from [40], with permission of APS).
Biomimetics 05 00044 g004
Figure 5. (a) Von Mises stress contour; (b) factor of safety (FOS) contour.
Figure 5. (a) Von Mises stress contour; (b) factor of safety (FOS) contour.
Biomimetics 05 00044 g005
Figure 6. Optimized hip implant designs via lattice structures: (a) Voronoi; (b) Gyroid; (c) Schwarz Diamond.
Figure 6. Optimized hip implant designs via lattice structures: (a) Voronoi; (b) Gyroid; (c) Schwarz Diamond.
Biomimetics 05 00044 g006
Figure 7. Von Mises stress contours for: (a) Voronoi structure; (b) Gyroid structure; (c) Schwarz Diamond structure.
Figure 7. Von Mises stress contours for: (a) Voronoi structure; (b) Gyroid structure; (c) Schwarz Diamond structure.
Biomimetics 05 00044 g007
Figure 8. Graphical illustration of functionally graded Gyroid structure.
Figure 8. Graphical illustration of functionally graded Gyroid structure.
Biomimetics 05 00044 g008
Figure 9. Functionally graded hip implants for lattice structures: (a) Voronoi; (b) Gyroid; (c) Schwarz Diamond.
Figure 9. Functionally graded hip implants for lattice structures: (a) Voronoi; (b) Gyroid; (c) Schwarz Diamond.
Biomimetics 05 00044 g009
Figure 10. Von Mises stress contours of functionally graded hip implant with: (a) Voronoi structure; (b) Gyroid structure; (c) Schwarz Diamond structure.
Figure 10. Von Mises stress contours of functionally graded hip implant with: (a) Voronoi structure; (b) Gyroid structure; (c) Schwarz Diamond structure.
Biomimetics 05 00044 g010
Table 1. Typical ranges for hip implant’s design parameters.
Table 1. Typical ranges for hip implant’s design parameters.
Design ParametersTypical Values
Length of intramedullary stem120 mm–180 mm
Length of neck10 mm–40 mm
Head diameter22 mm–45 mm
Neck diameter13 mm–30 mm
Angle of head placement135°–145°
Table 2. Chemical composition of Inconel 718 (UNS N07718).
Table 2. Chemical composition of Inconel 718 (UNS N07718).
Chemical Requirements
NiFeMoMnSiCrC
Min (%)50Balance2.80.350.35170.08
Max (%)553.321
Table 3. Mechanical Properties of Inconel 718 (UNS N07718).
Table 3. Mechanical Properties of Inconel 718 (UNS N07718).
Mechanical PropertiesTypical Values
Density8.19 g/cm3
Elastic Modulus200 GPa
Poisson Ratio0.29
Yield Strength1100 MPa
Ultimate Yield Strength1375 MPa
Table 4. Loading of hip implant for different types of movement.
Table 4. Loading of hip implant for different types of movement.
Type of MovementMax. Load (% Weight)Max. Force on Hip Joint
Slow walking 2822075 N
Climbing upstairs3562620 N
Climbing downstairs3872850 N
Tripping7205300 N
Table 5. Factor of safety for three different topology optimization approaches.
Table 5. Factor of safety for three different topology optimization approaches.
Lattice StructureFactor of Safety (FOS)Ultimate Factor of Safety (FOSUlt)
Voronoi0.630.79
Gyroid0.780.98
Schwarz Diamond0.941.17
Table 6. Factor of safety for all examined versions of hip implant’s design.
Table 6. Factor of safety for all examined versions of hip implant’s design.
Hip Implant VersionsFactor of Safety (FOS)Ult. Factor of Safety (FOSUlt)Weight (g)
Solid2.53.12975
Topology optimized
Voronoi0.630.79600
Gyroid0.780.98
Schwarz Diamond0.941.17
Topology optimized and functionally graded
Voronoi1.011.26600
Gyroid1.792.24
Schwarz Diamond2.082.6

Share and Cite

MDPI and ACS Style

Kladovasilakis, N.; Tsongas, K.; Tzetzis, D. Finite Element Analysis of Orthopedic Hip Implant with Functionally Graded Bioinspired Lattice Structures. Biomimetics 2020, 5, 44. https://doi.org/10.3390/biomimetics5030044

AMA Style

Kladovasilakis N, Tsongas K, Tzetzis D. Finite Element Analysis of Orthopedic Hip Implant with Functionally Graded Bioinspired Lattice Structures. Biomimetics. 2020; 5(3):44. https://doi.org/10.3390/biomimetics5030044

Chicago/Turabian Style

Kladovasilakis, Nikolaos, Konstantinos Tsongas, and Dimitrios Tzetzis. 2020. "Finite Element Analysis of Orthopedic Hip Implant with Functionally Graded Bioinspired Lattice Structures" Biomimetics 5, no. 3: 44. https://doi.org/10.3390/biomimetics5030044

APA Style

Kladovasilakis, N., Tsongas, K., & Tzetzis, D. (2020). Finite Element Analysis of Orthopedic Hip Implant with Functionally Graded Bioinspired Lattice Structures. Biomimetics, 5(3), 44. https://doi.org/10.3390/biomimetics5030044

Article Metrics

Back to TopTop