Next Article in Journal
Enhancing Prostate and Bladder Cancer Treatment: Exploring the Synergistic Potential of Entecavir and 5-Fluorouracil Combinations
Previous Article in Journal
Comprehensive Analysis of the Genetic Variation in the LPA Gene from Short-Read Sequencing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Customized Lattice Structures Tailored to Mimic Patients’ Bone Anisotropic Properties and Microarchitecture for Joint Reconstruction Applications

by
Ahmed Sherif El-Gizawy
1,2,*,
Xuewei Ma
1,
Joshua C. Arnone
2 and
Ammar A. Melaibari
3,4
1
Mechanical & Aerospace Engineering, University of Missouri, Columbia, MO 65211, USA
2
Industrial Technology Development and Management (ITECH D&M) LLC, Columbia, MO 65203, USA
3
Mechanical Engineering Department, Faculty of Engineering, King AbdulAziz University, Jeddah 21589, Saudi Arabia
4
Center of Nanotechnology, King AbdulAziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
BioMed 2024, 4(2), 171-184; https://doi.org/10.3390/biomed4020014
Submission received: 18 March 2024 / Revised: 12 May 2024 / Accepted: 4 June 2024 / Published: 13 June 2024

Abstract

:
Existing implants used with Total Knee Arthroplasty (TKA), Total Hip Arthroplasty (THA), and other joint reconstruction treatments, have displayed premature failures and frequent needs for revision surgery in recent years, particularly with young active patients who represent more than 55% of all joint reconstruction patients. Bone cement and stress shielding have been identified as the major reasons for premature joint failures. A breakdown of the cement may happen, and revision surgery may be needed because of the aseptic loosening. The significant mismatch of stiffness properties of patient trabecular bones and metallic implant materials in joint reconstruction surgery results in the stress shielding phenomenon. This could lead to significant bone resorption and increased risk of bone fracture and the aseptic loosening of implants. The present project introduces an approach for development of customized cellular structures to match the mechanical properties and architecture of human trabecular bone. The present work aims at fulfilling the objectives of the introduced approach by exploring new designs of customized lattice structures and texture tailored to mimic closely patients’ bone anisotropic properties and that can incorporate an engineered biological press-fit fixation technique. The effects of various lattice design variables on the mechanical performance of the structure are examined through a systematic experimental plan using the statistical design of experiments technique and analysis of variance method. All tested lattice designs were explored under realistic geometrical, biological, and manufacturing constraints. Of the four design factors examined in this study, strut thickness was found to have the highest percent contribution (41%) regarding the structure stiffness, followed by unit cell type, and cell size. Strut shape was found to have the lowest effect with only 11% contribution. The introduced solution offers lattice structure designs that can be adjusted to match bone stiffness distribution and promote bone ingrowth and hence eliminating the phenomenon of stress shielding while incorporating biological press-fit fixation technique.

1. Introduction

Existing implants used with Total Knee Arthroplasty (TKA), Total Hip Arthroplasty (THA), and other joint reconstruction treatments, have displayed premature failures and frequent needs for revision surgery in recent years, particularly with young active patients [1,2]. Ulrich et al. [1], found that the reviewed 237 THA revision procedures showed the overall mean time to revision was 83 months (close to 7 years), with failure reasons of aseptic loosening (51.9%), instability (16.9%), and infection (5.5%). Aseptic loosening is the most concern of revision rates for both TKA and THA. According to the historical growth trajectory, demand for joint replacements for young patients less than 65 years old was projected to be 52% for primary THAs and 55% to 62% for primary or revision TKAs by 2030 [2]. In recent years, researchers have discovered an effective solution for repairing bone defects caused by diseases or trauma, using the power of additive manufacturing (3D-printing) and the outstanding properties of titanium alloys, stainless steels and others. Specifically, Ti6Al4V is a titanium alloy with excellent biocompatibility, corrosion resistance, and combination of good load carrying capacity and toughness, make it the preferred choice for manufacturing of medical devices for orthopedic and dental applications [3,4,5].
Most joint replacements have severe stress shielding problems led by metal implants that have much higher stiffness (elastic modulus) than that of bone. When metal implants are inserted into bone structures, bone suffers a lower strain energy density than the previous intact bone, resulting in stress shielding. This stress shielding phenomenon can cause significant bone resorption, which leads to a significantly increased risk of fracture and aseptic loosening [6,7,8]. Chmielewska and Dean [9] presented a review highlighting recent developments in devices for skeletal reconstruction that match the stiffness, while not interrupting the normal loading pattern of a healthy bone, and help to combat stress shielding and stress concentration. The presented review summarizes various approaches to achieve stiffness matching. Their review concluded that the introduced device must minimize, if not eliminate, interference with the normal stress–strain trajectories that the bone’s shape and mechanical properties, as well as the adjacent muscles, are adapted to. Arabnejad et al. [8], showed that a fully porous implant with an optimized material micro-structure can reduce the amount of bone loss secondary to stress shielding by 75% compared to a fully solid implant. Zimmer introduced Trabecular Metal made of porous Tantalum that has an average modulus of elasticity from 2.5 to 3.9 GPa, compared with the reported modulus range of 0.1 to 1.5 GPa for cancellous bone [10]. The introduced trabecular structure can only reduce the mismatch in average to 4:1, still high enough to cause stress shielding. Although bone cement had been the gold standard in the field of joint replacement surgery, its use has somewhat decreased because of the advent of press-fit implants, which encourage bone ingrowth [11]. A press-fit implant is cementless and has porous features that allow biological fixation, which is important for young active patients. In a recent investigation [12], the damage resistance of four suggested lattice structures (regular hexahedron, regular octahedron, rhombic dodecahedron, and body-centered cubic), was studied using a gradient algorithm for stabilized and lightweight mandibular prostheses. The introduced algorithm was aided by finite element analysis to select an optimal lattice structure. The reported investigation had several limitations that would diminish the effectiveness of the proposed solution. The study assumed that fabricated orthopedic or dental implants would have homogeneous and isotropic properties, which is very difficult to achieve using additive manufacturing techniques. Secondly, the investigation considered only the effect of the type of lattice structure (four different types) and never looked at the effects of other lattice structure design parameters that could have significant effects in controlling stiffness properties and morphology of the porosity of the developed lattice structures.
Based on the above discussion, the present work aims to explore new designs of customized lattice structures and texture tailored to closely mimic patients’ bone anisotropic properties and that can incorporate an engineered biological press-fit fixation technique. The introduced customized lattice structure will consider all design parameters that contribute to controlling stiffness properties and porosity of the developed customized structure. These designs could be created using additive manufacturing (AM) techniques for use with implants that would overcome most of the issues related to aseptic loosening and have greater success, especially in younger patients.

2. Materials and Methods

Ma and El-Gizawy [13] presented a conceptual approach for development of customized cellular structures to match the mechanical properties of human trabecular bone. Figure 1 shows key steps in our proposed approach needed for building an implant with properties mimicking those of patients’ bones. The workflow in the present approach is feasible for both patient-specific implant and population-based implant design.
Module 1 in the proposed approach involves a CT scan of the patients’ native bone segmentation to get the 3D CAD models of the bone and characterization of geometry, microarchitecture, and stiffness properties of bones. These properties are crucial for designing and building optimized implants for joint and dental reconstruction applications. Module 2 covers the development of customized lattice structures tailored to mimic closely patients’ bone anisotropic properties and microarchitecture. Module 3 involves integration of porosity graded lattice to full scale implant design and manufacturing. Ultimately, all findings of the three modules will be utilized in module 4, implementation of the developed technology for use with newly developed implants.

2.1. Characterization of Geometry, Microarchitecture, and Stiffness Properties of Bones

The characterization of the geometry, microarchitecture, and stiffness properties of bones, required in module 1, was accomplished and published in the journal BioMed [14]. A summary of the findings is presented here to link them to the new work of module 2 presented in the current article.
Module 1 work presents evaluation and experimental verification using micro-CT scans for accurate characterization of geometry, microarchitecture, and stiffness properties of bones. The micro-CT scans of investigated bone reveal that the trabecular bone is highly anisotropic and heterogeneous. The results also showed a considerable degree of parametric variability and uncertainty on microarchitecture and stiffness properties of patient’s trabecular bone. Probabilistic analyses of micro-CT data could consider these variabilities. They can aid in generating the required lattice structures of optimum implant designs that match closely patient bones and hence can overcome most of the issues related to aseptic loosening phenomenon and premature failure of the treated joints. Figure 2 shows trabecular thickness and trabecular spacing colored maps and histograms for typical cylindrical bone sample. Figure 3 displays bone properties assigned based on the gray scale values obtained from micro-CT data. The predicted modulus varied significantly from point to point within the tested sample with an average of 506.7 MPa and standard deviation of 238.06 MPa. [14]. These observations mean that the mentioned microarchitecture parameters: trabecular thickness and spacing have significant degree of variability.

2.2. Module 2: Effective Design of Lattice Structures

The following are the lattice structure important design factors that would affect the properties and performance of the developed lattice structure: unit cell type; unit cell size; strut shape; and strut thickness. Element software from Next-Gen Engineering Design Software (nTop) [15] is used to generate the 3D models of the lattices with the workflow shown in Figure 4. First, a solid part of which the boundary representing the design space (filling space) of the lattice is imported from a CAD system. Secondly, a unit cell type is selected from the existing library or from user-defined unit cell structures. Then, unit cell size is selected for the assigned unit cell. After that a topological structure without assigned thickness is generated based on the unit cell type and size, to fill the design space imported in the first step. Finally, a 3D lattice structure is thickened out after defining the desired strut shape and thickness.
The selected unit cell types for the present investigation are cube vertex centroid, cubic diamond, and cubic fluorite, shown in Figure 5. They were selected because their overhang angles are 35.3 degrees for all beams used in the three configurations. This satisfies self-support requirement for building the lattice structure using additive manufacturing by Laser Powder Bed Fusion (L-PBF) method. The second design factor is strut thickness, which is defined as the diameter of the equivalent circle that has the same area. Considering the minimum printable thickness for the L-PBF process of 200 μm, strut shapes of square, hexagon, and circle are considered. If the strut thickness is a, a circle’s diameter is a; a hexagon’s width and height are 0.95a and 1.10a, respectively; and the side length of a square is 0.89a, to have the same cross-sectional area for the same strut thickness. A unit cell size of 1.2 mm, 1.4 mm, and 1.6 mm, and a strut thickness of 0.23 mm, 0.28 mm, and 0.33 mm were considered. Table 1 displays all factors and their alternative levels considered in present design.

2.3. Parametric Investigation of Lattice Structure Design

To investigate the effects of the lattice structure important design factors on the quality characteristics (properties and performance of the developed cellular structure), a standard orthogonal array for parametric investigation is used [16,17]. In this type of experimental design, partial factorial design of L9 (nine experiments to be conducted) is selected based on the number of factors and their levels to be tested. Three different samples were tested for each investigated lattice structure. The quality characteristics to be considered in this parametric study are the stiffness and yield strength properties of the lattice structure, subject to biological and manufacturing constraints dictated by the application of the device and the available manufacturing technology (L-PBF) for building these designs. Figure 6 shows the basic elements of the L-PBF process used in the present study. Table 2 displays the experimental Log of the used orthogonal array with all the four design parameters used and their levels. Figure 7 displays digital design models of the nine investigated lattice designs.

2.4. Materials and Fabrication of the Experimental Lattice Structures

The 3D printing Laser powder bed fusion (L-PBF) process system displayed in Figure 6 is used for printing the examined designs. It is an EOS M290 machine available at SIM Surgical, LLC, (https://simsurgical.com/ (accessed on 4 May 2018)), Saint Louis, MO, USA. The used EOS M290 system has a building volume of 250 mm × 250 mm × 325 mm, with a 400W Yb-fiber laser that allows a scanning speed up to 7 m/s. All the designed parts are built with Ti64 material with optimized process parameters for lattice structures. The Ti64 powder has a particle size of 39 ± 3 µm according to analysis following the ISO 13320 Particle size analysis-laser diffraction methods. The chemical composition analysis of the used material is shown in Table 3. After the printing, the powders around the built parts are removed and the support structures are broken away. A high-speed rotary tool with conical carbide cutter was used to deburr the rough surfaces. A glass bead blasting process was conducted to remove the powders inside the lattice structures and improve the surface finish throughout the parts. Ultrasonic washing was done in various cycles with Steris CIP 100 cleaning solution and deionized water with elevated temperature. Finally, hot isostatic pressing (HIP) process was conducted to improve the mechanical performance of the parts according to the ASTM F3301. The HIP process was done in argon inert gas environment with a temperature of 954 ℃, and pressure of 103 MPa for 3 h.

2.5. Evaluation of Geometric Features of Printed Lattices

Geometrical features such as surface area, pore size, and porosity are recorded for each lattice design. Pore size is evaluated in the cross-sectional plane that has the largest void area, by measuring the largest diameter of the circle that can fit into the voids. Surface area, pore size, and porosity were evaluated on the as-designed models. Porosity was calculated by the volume fraction of the void space on the CAD models.

2.6. Characterization of Stiffness and Strength of Investigated Lattice Structures

The MTS Landmark system was used for the compression testing of all investigated lattice structures (see Figure 8). For accurate measurement of the specimen deformation, a Dantec Q-400 Digital Image Correlation (DIC) system was used for the measurement of the displacement of the bottom platen (see Figure 8), while the upper platen was fixed during testing. A light source was projected onto the testing area when two cameras took images from two different angles. The DIC sampling rate was 10 Hz, i.e., the cameras took 10 images per second to capture the displacement change. A speckle pattern was sprayed in the front face of the DIC tracker with white color background and finely distributed black dots. The two cameras of the DIC system track the location of the black dots through taking pictures while the compression testing was running. During the tests, the bottom compression plate moved at a speed of 0.03 mm/s upward to compress the lattices. The compression test setup is shown in Figure 8 and Figure 9. Compression forces during testing were collected from the data acquisition system of the MTS testing machine. Displacement (deformation) data was calculated using image processing in ISTRA 4D software evaluation module. All collected data of forces and deformation for different tests were normalized to determine the stress–strain results associated with each test. Stress–strain results for each test were plotted to illustrate the mechanical behavior of each investigated lattice structure. Figure 10 displays stress–strain curves for three different lattice structures. Each of the displayed curves is unique and they are different from each other. The displayed curves show peaks and valleys very similar to the results obtained in compression testing of trabecular bone displayed in module 1 [14]. The slop of the linear portion of the stress–strain curve determines the elastic modulus of each lattice structure. At 0.2%, offset strain the corresponding yield strength of the specimens were determined. The exact measured properties for design (Exp.) # 6 are displayed in Table 4 as an example of the extracted properties of investigated designs. All other results are available with the corresponding author.

3. Results and Discussion

Table 5 summarizes the mean values of all measured properties of all nine lattices investigated with three samples tested for each design. All results are discussed in the light of how they affect the performance of patient’s specific lattice structures. Special attention was given to the significance of design parameters on the stiffness property (modulus of elasticity) because of its direct effect on the phenomenon of stress shielding and implant potential failure.

3.1. Parametric Analysis of Design Parameter Effect on Lattice Microstructure

The output results (Table 5) show that different lattice designs have different quality characteristics that would enable the designer to choose the design that would satisfy most of the function requirements of the implant. The followings are explanation of the various effects.

3.1.1. Surface Area

Surface area affects the bone ingrowth as it provides an interface for the biological reaction. The surface area of the 18 × 18 × 24 mm block was 2376 mm2, while the lattice surface areas varied from 14,264.43 to 33,566.66 mm2. This is 6 to 14 times the block surface area. Unit cell type and unit cell size have the largest influence on the surface area, and the maximum surface area occurs when the unit cell type is cubic fluorite, which has the highest density of beams, and when unit cell size is 1.2 mm, the smallest size. The specific levels are listed in the Experimental Log (Table 2).

3.1.2. Pore Size

Pore size is the geometric property that would affect the growth of bone tissues into the implant to establish a biological bond between the implant and the patient bone without need to introduce cement at the interface. Taniguchi et al. [18] have concluded that additive manufactured porous titanium implants with a pore size of 0.490–1.100 mm exhibit a reasonable performance in vivo. The reported bore size (Table 5) ranged from 0.720 to 1.580 mm. Six out of the investigated nine lattices were within the recommended ranges. Based on the reported results, to reduce the pore size, a higher level of unit cell types (cubic diamond or cubic fluorite), and smaller unit cell size are recommended.

3.1.3. Porosity Distribution

The optimum porosity was selected based on the resulting stiffness (modulus of elasticity), which should be close to the stiffness of the patient bone to avoid the stress shielding phenomenon. The experimental measurements (Table 5) of modulus of elasticity indicated that porosity in the range of 80–89% results in elastic modulus in the range of 1124.0 to 616.0 MPa. This range matches the ones for cancellous bone. Six out of nine investigated lattices were in this range.

3.2. Parametric Analysis of Structure Parameters’ Effects on Lattice Stiffness

The stiffness properties (Modulus of Elasticity) of the implant have the most effects on its performance. Most of the existing joint metal implants could lead to stress shielding problems because of their much higher elastic modulus than that of bone. This could lead to reduction in the quality of the remaining bone, which would result to a significantly increased risk of fracture and aseptic loosening. Therefore, more detailed analysis of the experimentally determined stiffness of the investigated lattice structures are given below to examine the significance of the effects of each of the lattice design parameters on the structure stiffness. Table 6 displays modulus of elasticity mean value and standard deviation for each of the nine tested lattices. The reported stiffness varied in the range of 616.5 MPa to 5357.4 MPa depending on lattice design.

Analysis of Variance

The analysis of variance (ANOVA) technique [16] is used to predict the relative significance of the lattice design parameters on the stiffness (modulus of elasticity) of the structure. It gives the percentage of contribution of each factor, thus providing a quantitative measure of the effects of various factors on design performance. The ANOVA results are shown below in Table 7. The error variance in the present analysis was obtained by pooling the sum of squares corresponding to the factor having the lowest mean square. Percentage contribution (ϱ, %) is calculated from the ratio of the Sum of Squares due to the factor/Total Sum of Squares. Of the four design factors examined in this study, the Strut Thickness (D) was found to have the highest percent contribution (41%) regarding the structure stiffness, followed by Unit Cell Type (A), and Cell Size (B). Strut Shape (C) was found to have the lowest effect with only 11% contribution. The variance ratio, denoted by “F”, in Table 7, is the ratio of the mean square due to a factor and the error mean square. A large value of F means the effect of that factor is large compared to the error variance.

4. Conclusions

The present work aimed at exploring new designs of customized lattice structures and texture tailored to mimic closely patients’ bone anisotropic properties and that can incorporate an engineered biological press-fit fixation technique, to overcome most of the issues related to aseptic loosening phenomenon and risk of bone fracture.
An effective approach for development of customized cellular structures to match the mechanical properties and architecture of human trabecular bone is introduced. It consists of four modules:1—characterization of geometry, microarchitecture, and stiffness properties of patient’s bones; 2—developing customized lattice structures tailored to mimic properties established in module 1; 3—integration of developed lattice structure into full scale implant design and manufacturing; 4—implementation of the developed technology for use in implants design and fabrication.
The developed approach is integrated with parametric investigation and involves both biological requirements and design for additive manufacturing by the laser powder bed fusion process (L-PBF).
Of the four design factors examined in this study, strut thickness was found to have the highest percentage contribution (41%) regarding the structure stiffness, followed by unit cell type, and cell size. Strut shape was found to have the lowest effect, with only 11% contribution.
The presented solution offers lattice structure designs that can be adjusted to match bone stiffness distribution and promote bone ingrowth and hence eliminating the phenomenon of stress shielding while incorporating biological press-fit fixation technique.
Further investigation is needed to establish a predictive model based on the present experimental results to define optimum lattice design that will simultaneously meat target value for each of the required quality characteristics while minimizing variability of performance, that is, having a robust lattice design. This is necessary to complete all requirements of module 2 of the introduced approach.

Author Contributions

A.S.E.-G., made substantial contributions to the research design, planning, coordination and management of all efforts, securing resources and drafting and submitting the final manuscript of this paper. X.M. made substantial contributions to the research design, analysis and experimental work as part of his graduate research at the University of Missouri. J.C.A. made substantial contributions to research design including design, optimization, and supervising of the L-PBF process and finishing operations used for building all investigated designs with the required properties. A.A.M. made substantial contributions to the present research, updating and reviewing literature survey supporting the present results, securing additional resources from King Abdulaziz University paying for materials and supplies needed, and reviewing and editing the final manuscript prior of submission for publication. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study did not require ethical approval.

Informed Consent Statement

No information was obtained from known patients or animals.

Data Availability Statement

All results regarding mechanical testing, geometrical measurements and statistical analysis of the results are available with A. S. El-Gizawy (corresponding author for the present article, and X. Ma (co-author of this article).

Acknowledgments

The authors wish to acknowledge the in-kind and financial support provided by Industrial Technology Development and Management (ITECH D&M), LLC, Columbia, Missouri, USA, for this research. Additionally, much gratitude is extended to Talissa Altes (MD) and Staff of the University of Missouri Department of Radiology for sharing their CT scanner and expertise. Our appreciation is also extended to Robert Weinholtz, Frank Feng, and Matthew Maschmann of the University of Missouri Department of Mechanical and Aerospace Engineering for their valuable discussions and encouragement during the present project. We acknowledge also advice received from Benjamin Hansen (MD), and Ajay Aggarwal (MD) of Missouri Orthopedic Institute (MOI) regarding medical needs and importance of the proposed solutions.

Conflicts of Interest

All authors of the present article have declared no conflict of interest during the submission of our paper. Ahmed Sherif El-Gizawy and Joshua C. Arnone are with ITECH D&M, LLC. Ahmed Sherif El-Gizawy is the president and Chief Technical Officer of this consulting company. Ahmed Sherif El-Gizawy has received several speaker honorariums over the last 30 years. And he owned stocks with TIAA-Craft in the USA. Joshua C. Arnone is the inventor of several patents.

References

  1. Ulrich, S.D.; Seyler, T.M.; Bennett, D.; Delanois, R.E.; Saleh, K.J.; Thongtrangan, I.; Kuskowski, M.; Cheng, E.Y.; Sharkey, P.F.; Parvizi, J.; et al. Total hip arthroplasties: What are the reasons for revision? Int. Orthop. 2008, 32, 597–604. [Google Scholar] [CrossRef] [PubMed]
  2. Kurtz, S.M.; Lau, E.; Ong, K.; Zhao, K.; Kelly, M.; Bozic, K.J. Future young patient demand for primary and revision joint replacement: National projections from 2010 to 2030. Clin. Orthop. Relat. Res. 2009, 467, 2606–2612. [Google Scholar] [CrossRef] [PubMed]
  3. Hutmacher, D.W. Scaffolds in tissue engineering bone and cartilage. Biomaterials 2000, 21, 2427–2431. [Google Scholar] [CrossRef]
  4. Van Bael, S.; Kerckhofs, G.; Moesen, M.; Pyka, G.; Schrooten, J.; Kruth, J.P. Micro-CT-based improvement of geometrical and mechanical controllability of selective laser melted Ti6Al4V porous structures. Mater. Sci. Eng. A 2011, 528, 7423–7431. [Google Scholar] [CrossRef]
  5. Abdelrhman, Y.; Gepreel, M.A.-H.; Kobayashi, S.; Okano, S.; Okamoto, T. Biocompatibility of new low-cost (α + β)-type Ti-Mo-Fe alloys for long-term implantation. Mater. Sci. Eng. C Mater. Biol. Appl. 2019, 99, 552–562. [Google Scholar] [CrossRef] [PubMed]
  6. Huiskes, R.; Weinans, H.; Van Rietbergen, B. The relationship between stress shielding and bone resorption around total hip stems and the effects of flexible materials. Clin. Orthop. Relat. Res. 1992, 274, 124–134. [Google Scholar] [CrossRef]
  7. Mi, Z.R.; Shuib, S.; Hassan, A.Y.; Shorki, A.A.; Ibrahim, M.M. Problem of stress shielding and improvement to the hip implant designs: A review. J. Med. Sci. 2007, 7, 460–467. [Google Scholar] [CrossRef]
  8. Arabnejad, S.; Johnston, B.; Tanzer, M.; Pasini, D. Fully porous 3D printed titanium femoral stem to reduce stress-shielding following total hip arthroplasty. J. Orthop. Res. 2017, 35, 1774–1783. [Google Scholar] [CrossRef] [PubMed]
  9. Chmielewska, A.; Dean, D. The role of stiffness-matching in avoiding stress shielding-induced bone loss and stress concentration-induced skeletal reconstruction device failure. Acta Biomater. 2024, 173, 51–65. [Google Scholar] [CrossRef] [PubMed]
  10. Dental Solutions|ZimVie America. Available online: https://www.zimmerbiometdental.com (accessed on 20 January 2017).
  11. Vaishya, R.; Chauhan, M.; Vaish, A. Bone cement. J. Clin. Orthop. Trauma 2013, 4, 157–163. [Google Scholar] [CrossRef] [PubMed]
  12. Liu, R.; Su, Y.; Yang, W.; Wu, K.; Du, R.; Zhong, Y. A Novel Design Method of Gradient Porous Structure for Stabilized and Lightweight Mandibular Prosthesis. Bioengineering 2022, 9, 424. [Google Scholar] [CrossRef] [PubMed]
  13. Ma, X.; El-Gizawy, A.S. Development of customized cellular structures to match the mechanical properties of human trabecular bone. In Proceedings of the 2018 Annual International Solid Freeform Fabrication Symposium—An Additive Manufacturing Conference, Austin, TX, USA, 13–15 August 2018. [Google Scholar]
  14. El-Gizawy, A.S.; Ma, X.; Pfeiffer, F.; Schiffbauer, J.D.; Selly, T. Characterization of Microarchitectures, Stiffness and Strength of Human Trabecular Bone Using Micro-Computed Tomography (Micro-CT) Scans. BioMed 2023, 3, 89–100. [Google Scholar] [CrossRef]
  15. Next-Gen Engineering Design Software: nTop. Available online: www.ntop.com (accessed on 2 July 2017).
  16. Phadke, M.S. Quality Engineering Using Robust Design; Prentice Hall: Englewood, CA, USA, 1989. [Google Scholar]
  17. Arnone, J.A.; El-Gizawy, A.S.; Crist, B.D.; Della Rocca, G.J.; Ward, C.V. Computer-Aided Engineering Approach for Parametric Investigation of Locked Plating Systems Design. J. Med. Devices 2013, 7, 021001. [Google Scholar] [CrossRef]
  18. Taniguchi, N.; Fujibayashi, S.; Takemoto, M.; Sasaki, K.; Otsuki, B.; Nakamura, T.; Matsushita, T.; Kokubo, T.; Matsuda, S. Effect of pore size on bone ingrowth into porous titanium implants fabricated by additive manufacturing: An in vivo experiment. Mater. Sci. Eng. C 2016, 59, 690–701. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Proposed approach for development of customized orthopedic implants.
Figure 1. Proposed approach for development of customized orthopedic implants.
Biomed 04 00014 g001
Figure 2. Trabecular thickness and trabecular spacing colored maps and histograms for typical cylindrical bone samples [14].
Figure 2. Trabecular thickness and trabecular spacing colored maps and histograms for typical cylindrical bone samples [14].
Biomed 04 00014 g002
Figure 3. Density and stiffness properties (MPa) assigned based on the gray scale values [14].
Figure 3. Density and stiffness properties (MPa) assigned based on the gray scale values [14].
Biomed 04 00014 g003
Figure 4. The workflow of creating a periodic lattice structure with top elements.
Figure 4. The workflow of creating a periodic lattice structure with top elements.
Biomed 04 00014 g004
Figure 5. Three-unit cell types investigated. The height (H), length (L) and width (W) of the external cubes are equal. Cell size = H = L = W.
Figure 5. Three-unit cell types investigated. The height (H), length (L) and width (W) of the external cubes are equal. Cell size = H = L = W.
Biomed 04 00014 g005
Figure 6. Laser powder bed fusion (L-PBF) system used for building all investigated samples.
Figure 6. Laser powder bed fusion (L-PBF) system used for building all investigated samples.
Biomed 04 00014 g006
Figure 7. Digital design models of the investigated nine lattice structures.
Figure 7. Digital design models of the investigated nine lattice structures.
Biomed 04 00014 g007
Figure 8. Compression test setup of lattice structures.
Figure 8. Compression test setup of lattice structures.
Biomed 04 00014 g008
Figure 9. Close view of compression test fixture displaying the specimen and the DIC tracker that was used for measuring the axial deformation in the specimen.
Figure 9. Close view of compression test fixture displaying the specimen and the DIC tracker that was used for measuring the axial deformation in the specimen.
Biomed 04 00014 g009
Figure 10. Stress–strain curves for three different lattice structures. (a) Stress–strain curve for cube vertex centroid lattice; (b) Stress–strain curve for cubic diamond lattice; and (c) Stress–strain curve for cubic fluorite lattice.
Figure 10. Stress–strain curves for three different lattice structures. (a) Stress–strain curve for cube vertex centroid lattice; (b) Stress–strain curve for cubic diamond lattice; and (c) Stress–strain curve for cubic fluorite lattice.
Biomed 04 00014 g010
Table 1. Factors and their alternative levels considered in present lattice design.
Table 1. Factors and their alternative levels considered in present lattice design.
Factors
(Design Parameters)
Levels
123
Unit Cell Type (A)Cube Vertex CentroidCubic DiamondCubic Fluorite
Cell Size (B), mm1.21.41.6
Strut Shapes (C)SquareHexagonCircle
Strut Thickness (D), mm0.230.280.33
Table 2. Experimental log for parametric analysis.
Table 2. Experimental log for parametric analysis.
Experimental (Measurements)-Log
Lattice
Design
(Experiment)
#
Unit Cell Type
A
Cell Size,
mm
B
Strut
Shape
C
Strut Thickness,
mm
D
1Cube Vertex Centroid1.2Square0.23
2Cube Vertex Centroid1.4Hexagon0.28
3Cube Vertex Centroid1.6Circle0.33
4Cubic Diamond1.2Hexagon0.33
5Cubic Diamond1.4Circle0.23
6Cubic Diamond1.6Square0.28
7Cubic Fluorite1.2Circle0.28
8Cubic Fluorite1.4Square0.33
9Cubic Fluorite1.6Hexagon0.23
Table 3. Material chemical composition of Ti64 powder (wt%).
Table 3. Material chemical composition of Ti64 powder (wt%).
EOS Ti64 Material Composition
ElementTiAlVONCHFe
wt%balance5.5–6.753.5–4.5<0.2<0.05<0.08<0.015<0.3
Source: EOS Titanium Ti64 Material Data Sheet. EOS GmbH—Electro Optical Systems.
Table 4. Measured properties and their mean and standard deviation (SD) of tested samples for Design #6.
Table 4. Measured properties and their mean and standard deviation (SD) of tested samples for Design #6.
Expe. (Design) #Sample #As Designed Porosity (%)Measured Modulus
(MPa)
Mean
Modulus
(MPa)
SD
Modulus
Measured Yield Strength (MPa)Mean Strength (MPa)SD Strength
Exp. #6Sample #185.0992.01124.4164.322.825.65.1
Sample #21072.922.6
Sample #31308.331.5
Table 5. Summary of all measured properties (Mean Values) of tested lattices.
Table 5. Summary of all measured properties (Mean Values) of tested lattices.
Lattice
Design #
Surface
Area, mm2
Pore Size, mmPorosity, %Compression Testing
Elastic Modulus, MPaYield Strength, MPa
Exp#121,344.551.2383.5708.117.0
Exp#218,352.451.4180.9706.919.6
Exp#314,264.431.5883.5616.516.8
Exp#425,367.250.7469.22513.358.2
Exp#515,098.331.0889.0688.116.8
Exp#615,296.341.1285.01124.425.6
Exp#733,566.660.7461.43345.999.6
Exp#831,190.340.7252.35357.477.4
Exp#921,075.151.1781.11248.530.3
Table 6. Modulus of elasticity mean value and standard deviation for each of tested lattice design.
Table 6. Modulus of elasticity mean value and standard deviation for each of tested lattice design.
Experimental (Measurements)—LogMean
Modulus/Mpa
Standard Deviation
Modulus
Lattice
Design
#
Unit Cell
Type
Cell Size,
mm
Strut
Shape
Strut
Thickness,
mm
µσ
1Cube
Vertex Centroid
1.2Square0.23708.176.5
2Cube
Vertex Centroid
1.4Hexagon0.28706.956.2
3Cube
Vertex Centroid
1.6Circle0.33616.550.9
4Cubic
Diamond
1.2Hexagon0.332513.3398.3
5Cubic
Diamond
1.4Circle0.23688.173.4
6Cubic
Diamond
1.6Square0.281124.4164.3
7Cubic
Fluorite
1.2Circle0.283345.9861.3
8Cubic
Fluorite
1.4Square0.335357.43055.8
9Cubic
Fluorite
1.6Hexagon0.231248.5135.4
Table 7. Analysis of variance regarding the effects design parameters on stiffness of the developed design.
Table 7. Analysis of variance regarding the effects design parameters on stiffness of the developed design.
ANOVA Table for Modulus of Elasticity
SourceDegrees
of Freedom
D.F
Sum of
Squares,
Sx
Mean
Square
Vx = Sx/D.F
Variance Ratio,
F
Percentage Contribution,
ϱ (%)
Unit Cell Type
A
23,863,4491,931,7242.527
Cell Size
B
22,999,3701,499,6851.9421
Strut Shape
C
21,542,907771,4541.011
Strut Thickness
D
25,722,7522,816,3713.6541
Totals814,128,4781,766,059 100
(Error)(2)(1,542,907)(771,454)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

El-Gizawy, A.S.; Ma, X.; Arnone, J.C.; Melaibari, A.A. Customized Lattice Structures Tailored to Mimic Patients’ Bone Anisotropic Properties and Microarchitecture for Joint Reconstruction Applications. BioMed 2024, 4, 171-184. https://doi.org/10.3390/biomed4020014

AMA Style

El-Gizawy AS, Ma X, Arnone JC, Melaibari AA. Customized Lattice Structures Tailored to Mimic Patients’ Bone Anisotropic Properties and Microarchitecture for Joint Reconstruction Applications. BioMed. 2024; 4(2):171-184. https://doi.org/10.3390/biomed4020014

Chicago/Turabian Style

El-Gizawy, Ahmed Sherif, Xuewei Ma, Joshua C. Arnone, and Ammar A. Melaibari. 2024. "Customized Lattice Structures Tailored to Mimic Patients’ Bone Anisotropic Properties and Microarchitecture for Joint Reconstruction Applications" BioMed 4, no. 2: 171-184. https://doi.org/10.3390/biomed4020014

APA Style

El-Gizawy, A. S., Ma, X., Arnone, J. C., & Melaibari, A. A. (2024). Customized Lattice Structures Tailored to Mimic Patients’ Bone Anisotropic Properties and Microarchitecture for Joint Reconstruction Applications. BioMed, 4(2), 171-184. https://doi.org/10.3390/biomed4020014

Article Metrics

Back to TopTop