Next Article in Journal
Progress and Current Limitations of Materials for Artificial Bile Duct Engineering
Previous Article in Journal
Tribo-oxide Competition and Oxide Layer Formation of Ti3SiC2/CaF2 Self-Lubricating Composites during the Friction Process in a Wide Temperature Range
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 14
Issue 23
10.3390/ma14237467
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Element Analysis on Initial Crack Site of Porous Structure Fabricated by Electron Beam Additive Manufacturing

by
Meng-Hsiu Tsai
1,2,†,
Chia-Ming Yang
3,†,
Yu-Xuan Hung
1,
Chao-Yong Jheng
3,
Yen-Ju Chen
4,
Ho-Chung Fu
1,* and
In-Gann Chen
3,*
1
Department of Mold and Die Engineering, National Kaohsiung University of Science and Technology, Kaohsiung 807618, Taiwan
2
School of Dentistry, College of Dental Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
3
Department of Materials Science and Engineering, National Cheng-Kung University, Tainan 701, Taiwan
4
Metal Processing R & D Department, Metal Industries Research and Development Centre (MIRDC), Kaohsiung 811, Taiwan
*
Authors to whom correspondence should be addressed.
†
Meng-Hsiu Tsai and Chia-Ming Yang contributed equally to this manuscript as the first authors.
Materials 2021, 14(23), 7467; https://doi.org/10.3390/ma14237467
Submission received: 31 October 2021 / Revised: 28 November 2021 / Accepted: 30 November 2021 / Published: 6 December 2021
(This article belongs to the Topic Material and Process Innovations for 3D Printing Applications)
Browse Figures
Versions Notes

Abstract

:
Ti6Al4V specimens with porous structures can be fabricated by additive manufacturing to obtain the desired Young’s modulus. Their mechanical strength and deformation behavior can be evaluated using finite element analysis (FEA), with various models and simulation methodologies described in the existing literature. Most studies focused on the evaluation accuracy of the mechanical strength and deformation behavior using complex models. This study presents a simple elastic model for brittle specimens followed by an electron beam additive manufacturing (EBAM) process to predict the initial crack site and threshold of applied stress related to the failure of cubic unit lattice structures. Six cubic lattice specimens with different porosities were fabricated by EBAM, and compression tests were performed and compared to the FEA results. In this study, two different types of deformation behavior were observed in the specimens with low and high porosities. The adopted elastic model and the threshold of applied stress calculated via FEA showed good capabilities for predicting the initial crack sites of these specimens. The methodology presented in this study should provide a simple yet accurate method to predict the fracture initiation of porous structure parts.

Graphical Abstract

Highlights

  • A combination of EBAM and FEA was employed for the analysis of initial crack sites.
  • The initial crack sites predicted by the elastic model were consistent with experimental results.
  • The threshold stresses calculated by FEA coincided with the crack sites of specimens with different porosities.
  • The collapse mechanism due to strut behavior under uniaxial compression stress was investigated.

1. Introduction

Porous materials, also called cellular solids, fabricated by additive manufacturing have attracted attention in biomedical implants, such as acetabular hip cups [1,2,3], mandibles [4], fusion cages [5,6], and various other applications [2]. Implants made from a material with a low elastic modulus can reduce the effect of stress shielding [7]. The use of porous metal structures can effectively eliminate this phenomenon. The amount of porosity in the implant is considered a crucial factor in promoting successful bone integration with a porous structure [8]. Conventional methods of fabricating porous materials include chemical and heat treatment [9], powder metallurgy [10], and thermal spraying [11]. However, these manufacturing techniques have limitations in terms of controlling the size, shape, and distribution of the pores. A state-of-the-art method for fabricating porous materials is powder bed fusion (PBF), which is a process of fusing powder materials to make objects according to a 3D model, usually by layer-upon-layer melting with a laser or electron beam heating source. Selective laser melting (SLM) and electron beam additive manufacturing (EBAM), as a kind of metal powder bed layer-by-layer fusion process, provide considerable freedom to designers for deciding the geometry of the final part and allow for a reduced manufacturing time.
SLM/EBAM-fabricated parts present porous structures with variations in topological designs and porosities (i.e., the volume ratio of the porous structure and bulk) to obtain the desired Young’s modulus. Gibson and Ashby [12] reported an exponential square formula for the relationship between the relative density and mechanical properties of porous structures, such as Ti6Al4V cubic-type porous structures with porosities ranging from 60% to 96.2%, with the final effective elastic modulus ranging from 23 MPa to 78 MPa [13]. For approximately the same porosities of 49.75% and 50.75%, with variations only in the strut thicknesses of specimens, their compressive stiffness and strength values decreased significantly from 2.92 GPa to 0.57 GPa (80.5% reduction) and from 163.02 MPa to 7.28 MPa (93.54% reduction), respectively [14]. Li et al. investigated a Ti6Al4V porous structure with different meshes (cubic, G7, and rhombic dodecahedron) fabricated by the EBAM process [15]. In that research, the influence of porosity was the main discussion point by observing specimens’ compressive mechanical properties at each density, However, the influence of porosity gradients was not discussed. Grunsven et al. studied a Ti6Al4V diamond-type lattice structure fabricated by EBAM with variations in the porosity [16]. The collapse of the porous gradient structure of specimens with cubic and honeycomb lattices started from a less dense layer and progressed to more dense layers, as reported by Choy et al. [17]. Most of these studies showed a power-law relationship between the porosity and elastic modulus, similar to that given by Gibson and Ashby’s equation. In addition, the plastic deformation behavior revealed the coupling of buckling and bending in the compression stress–strain curve.
Finite element analysis (FEA) is a technique commonly used to evaluate the mechanical behavior of porous materials. By using FEA, evaluations of overall and local stress distribution in a unit cell structure could be shown for further research purposes. Therefore, FEA-based simulations were used to estimate the mechanical properties and determine the failure modes of the lattice structures. For elastic modulus studies, some researchers reported the elastic modulus discrepancies between computational and experimental results [18,19]. The elastic modulus predicted with finite element modeling was as much as 10 times higher than that obtained via experiments with a Ti6Al4V simple cubic porous structure fabricated by EBAM [19]. It was proposed that only perfect elastic behavior, rather than elastic–plastic behavior, contributes to enlarging the difference between the theoretical and experimental results. The porosity of the porous structure can be changed by altering the strut or pore size. It has been reported that cell topologies, strut diameters, cell sizes, bulk geometries, and strut defects affect the mechanical properties owing to local stress concentrations [20,21,22,23]. As for the deformation of deformation behavior in simulation, Ushijima et al. [24] reported an FEA to predict the initial stiffness and plastic collapse strength of body-centered cubic (BCC) lattice structures under compressive loads. The analytical and FEA predictions were compared with the experimental results. A good agreement was found between the experimental data and the corresponding FEA predictions for structures with a low relative density. The collapse modes for BCC and BCC with vertical pillars were evaluated in the FEA models, and they were also in agreement with the experimental observations [25].
Most of these numerical simulations adopted complex nonlinear models to predict deformation behavior and compared the obtained results with microcomputed tomography (Micro-CT) images to prove the accuracy of evaluation, as well as determine the elastic modulus, collapse strength, and macroscale deformation behavior. There is limited research on the failure behavior of the initial site as a function of the stress distribution of porous titanium parts fabricated by EBAM process. Therefore, it is important to illustrate the stress distribution behavior to gain an understanding of the local stress concentration and fracture mechanism. Additionally, EBAM-fabricated porous specimens with martensitic microstructures that exhibit brittle properties have been previously reported by us [5] and other researchers [14,15]. This study aimed to predict the failure at initial sites and the threshold stress causing the failure of high-porosity specimens using a simple elastic model, which differs from the elastic–plastic FEAs in the literature reporting on porous Ti6Al4V compressive specimens. Thus, we present cubic porous titanium parts fabricated by EBAM technology, with different porosities ranging from 25% to 63%. Compression tests were performed on these EBAM-fabricated specimens to study the failure behavior, and in situ digital images were acquired. In addition to the aforementioned failure at the initial points, FEA results of the cubic-type porous structure were analyzed with regard to the local stress concentration for comparison with the experimental observations. The objective of this study was to evaluate the crack initiation sites of compressed Ti6Al4V porous parts fabricated by EBM, as well as the threshold stresses of porous specimens at fracture.

2. Materials and Methods

2.1. Specimen Preparation

To investigate the porosity and mechanical properties of EBAM-fabricated porous metallic parts, several specimens according to ASTM-E9 89a standards were designed by Inventor (2017) and manufactured by the EBAM process. The specimens were cubic-type porous structures with fixed cell sizes of 4 mm, nominal strut diameters of 0.6 to 1.25 mm, and pore diameters of 1.5 to 2.8 mm, making a total of six different porosity conditions of the specimens in Table 1; unit cell models are shown in Figure 1.

2.2. EBAM Process

In this study, specimens with different porosities were fabricated on an Arcam Q10 machine (Arcam EBM GE Additive Company, Sweden), the details of which were described in our previous study [5]. The chemical composition of the extra-low interstitial Ti-6Al-4V (Grade 23) followed the ASTM F3001 standard specification spherical powder and was supplied by Arcam AB Company, and the powder size distribution was quoted as 45–105 μm. The EBAM process was conducted in a vacuum chamber with a pressure below 5 × 10−3 mbar in the beginning, while the pressure at the finishing stage was 2 × 10−5 mbar. Each powder layer of 50 μm thickness was created by the raking powder gravity fed from two cassettes, preheated to 730 °C by fast scanning with a defocused electron beam, and the standard parameters for obtaining the porous structure were a beam current at 3 mA and scan speed of 1500 mm/s with a fixed beam diameter of 100 μm and hatch spacing of 150 μm, thereby melting the selected layer areas driven by the three-dimensional CAD model. The preheating and melting processes were achieved by the energy transfer from a high-energy electron beam onto the powder bed. The recycling of the nonmelted powder and/or sintered powder was achieved via a powder recovery system and a mechanical vibrating sieve with a mesh size of ≤150 μm.

2.3. Compressive Stress–Strain Testing

The static properties of the porous structures were measured using a static test machine (Shimadzu UH, 1000 kN load cell) by applying a constant deformation rate of 1.0 mm/min. Three specimens were tested for each variation in the porous structure. Uniaxial engineering compression stress–strain curves were obtained, and the mean and standard deviation of each of the five compressive properties were determined. During compressive testing, a digital camera was used to obtain in situ images. Scanning electron microscopy (SEM) secondary electron imaging was carried out on a JEOL JSM-6380 (Tokyo, Japan) at 15 kV for fracture surface observations.

2.4. Stress Distribution Simulation

Commercial FEA software COMSOL Multiphysics was adopted to perform the FEA. Instead of a complex computational model, this study uses an elastic model to simulate the collapse behavior of the porous specimens. The experimental flow chart is shown in Figure 2.
Our simulation model was based on the following assumptions:
  • The deformation of material obeys linear elastic fracture mechanics.
  • There were no initial defects inside the specimen, and material properties were uniform and isotropic [26].
On the basis of these assumptions, the elastic deformation was assumed to follow Hooke’s law: σ = C·ε, where C is the material’s stiffness matrix, which relates to the stresses, σ, and strains, ε.
The boundary conditions used for the model are shown in Figure 2: (1) the bottom surface was set as fixed constraint; (2) the compressive force was uniformly applied on the top surface of specimens; (3) all nodes and elements above the bottom surface were set as “free”.
To understand the stress distribution on specimens of different porosities, the applied static stresses of the specimens were set up on the basis of their yield strengths: 201, 166, 112, 80, 53, and 43 MPa for C1 to C6, respectively. The material setting of the simulation component was Ti6Al4V grade 23 with a Young’s modulus E0 of 112 GPa and a Poisson’s ratio of 0.37.
All CAD models were built by SOLIDWORKS, which accurately matched the actual configuration used in the EBM specimen setup. To confirm if the simulation settings could reach convergence status, three different mesh conditions were used for the C6 specimen with the minimum strut size designed to be 0.6 mm. The detailed mash information and images are shown in Table 2 and Figure 2b. Mesh 1 (physics-controlled mesh capability: fine) set the maximum element size of 2.4 mm and the minimum size of 0.3 mm, which is twofold larger than the minimum strut size of the C6 specimen. Under this setting, 47,565 tetrahedron meshes and 15,065 nodes were established. Mesh 2 (physics-controlled mesh capability: extra fine) set the maximum element size of 1.15 mm and the minimum size of 0.045 mm, which resulted in 151,909 tetrahedron meshes and 41,547 nodes. Mesh 3 set the maximum element size of 0.5 mm and the minimum size of 0.2 mm, which resulted in 1,341,348 tetrahedron meshes and 277,728 nodes. There were four integration points between two nodes for all conditions, and their stress calculation results are also shown in Table 2. The simulation difference for different mesh conditions was less than 3%, implying that the meshing condition of the mesh 1 was sufficient for this study; thus, the C2–C5 specimens used this as the mesh setting.
Moreover, to verify and confirm the strain of a component, two probes were placed on the top and bottom surfaces to obtain the position information. The specific Young’s moduli (E0) of the porous specimens were obtained from the set stress/probe detected strain and compared with the experimental results for the purpose of calibration. In this study, the calibrated E0 (17.32 GPa) was used to calculate the von Mises stress, including the stress contour and stress histogram. The stress contour (distribution) figures were directly calculated using COMSOL, and the raw data were expressed in the stress histograms.

3. Results and Discussion

3.1. Compression Stress–Strain Curves

Compression tests were conducted on six specimens with different porosities. The load–displacement and stress–strain curves of all specimens shown in Figure 3a,b were found to start with an elastic zone, where their stress increased linearly up to the first peak and then decreased, followed by a plateau zone with the strain ranging from 0.2 to 0.4, and then followed by multiple peaks in which the stress fluctuated, ending with a densification zone, where stress increased linearly. The fracture morphology was observed in the SEM images shown in Figure 3c. Ductile dimples and a smooth surface were found at the fracture surface for the C6 and C1 specimens, showing a decrease in the porosities. For C1, the fluctuation became large in the plateau zone, indicating that the specimens became more brittle. This phenomenon was consistent with the observations of Cheng et al. [27], in which a stochastic foam and reticulated mesh porous structures were fabricated by EBAM, and Choy et al. [28], where cubic and honeycomb porous structures were fabricated by SLM.

3.2. Compressive Deformation Behavior

The compressive deformation behavior of the porous specimens was recorded using a high-resolution digital single-lens reflex camera. Two types of deformation were found in the visual observation during the compression test, as shown in Figure 4. The specimens with low porosities, C1, C2, and C3, exhibited abrupt shear failure with one diagonal shear band forming throughout the entire specimen. For the typical low-porosity C1 specimen in Figure 4a, the first collapsed layer started from point 1 of the layer at the edge or adjacent to the edge of the specimens. The next layer to collapse was the adjacent layer, and subsequent collapse of the layers either followed a sequential pattern or occurred randomly at other layers in the strain range ɛ = 0 to 0.2 (point 2). When the strain was increased to 0.3 at point 3, cracks were initiated slightly toward the diagonal direction and diagonal shear band throughout the whole specimen at point 4. In contrast, the high-porosity specimens, C4, C5, and C6, deformed uniformly in a sequential layer-by-layer collapse manner. For the typical high-porosity C4 specimen in Figure 4b, it started with a deformation in the middle part of the vertical struts at point 5, followed by the collapse of all other layers at points 6 and 7. As a result, when the low-porosity specimens deformed as per the shear diagonal model, the stress–strain curves tended to have one large peak, followed by densification with some small fluctuations. This collapse behavior is in good agreement with that observed in porous Ti6Al4 studies [17,28,29]. Combining the strain–stress curve and compressive images, the schematic diagrams in Figure 4c,d illustrate the failure behavior in the low- and high-porosity specimens, respectively. For C1, in the initial period of the linear elasticity and nonlinear behavior, the edges of the V-struts and H-struts began to bend under the compressive loading. After slight bending, some of the struts experienced brittle fractures, resulting in the formation of a diagonal shear band throughout the specimen. For C4, the initial deformation started from the middle of the V-struts, resulting in a layer-by-layer collapse. To further understand the stress behavior of the specimens, stress simulation results are discussed in the next section.

3.3. Local Stress Concentration Region Analyzed by FEA

The experimental and FEA results of the relative Young’s modulus (E0 = 112 GPa) depend on the relative densities of all specimens, as shown in Figure 5. According to the Gibson–Ashby equation [12], their relationship can be described by E/E0 = C(d/d0)n, where E is the Young’s modulus for a specific density d, d0 is the density for a solid specimen, and C and n are constants for a specific structure. After fitting, the n values obtained were 2.02 and 2.08, respectively, for the FEA and experiment results, with the R-square values of 0.99 and 0.96, respectively. These results agree with the bending-dominated mechanism of the open cell (n = 2) [30]. However, the C values of the FEA and experimental results were approximately 1 and 0.16, respectively, implying that the E0 of the FEA should be smaller than that of the ideal Ti6Al4V specimen. For the difference in elastic modulus between the FEA and experiment, studies have reported that cell topologies, strut diameters, cell sizes, bulk geometries, strut defects, porosities, and phase transformations affect the mechanical properties owing to local stress concentrations [20,21,22,23]. To bring the FEA results closer to the experimental results of the real specimens, further simulation used the adjustment of E0 as 112 × 0.16 = 17.92 GPa. The calibrated FEA results in Table 3 presenting the elastic modulus for C1 to C6 indicate that the error value of the calibrated E2 was less than 13% compared to the experimental value of E.
Figure 6a shows the side view of the von Mises [31] stress distribution (stress counter) for specimens with different porosities. Local stress concentration sites were found at the edge or vertical strut adjacent to the edge of the horizontal strut in specimens C1 and C2. A transition stress distribution was observed in C3 specimen. However, the porosity increased by over 55% (i.e., C4, C5, and C6), with the stress concentration at the center of the strut. The x-shaped stress distribution behavior in C1 and C2 differed from the stress concentration at the center of the struts in the C4–C6 specimens. To observe the stress concentration site clearly, Figure 6b shows the 3D stress counter of the unit cells for different specimens. The maximum stress value was found on the nodes in C1 and C2, which differed from the site at the strut center in the C4, C5, and C6 specimens, implying that the initial crack site was dependent on the porosity for stress behavior. For C1, the initial crack site and x-shaped stress concentration were as shown in Figure 6a,b, and the initial crack sites started with diagonal cracks throughout the entire specimen, which is consistent with our experimental results in Figure 4a. Conversely, for the C6 specimen, the initial crack site was found at the center of the strut according to the FEA results in Figure 6; the sequence of the layer-by-layer collapse was also consistent with the experimental results in Figure 4b.
Sun et al. [32] studied the mechanical properties of porous materials fabricated by additive manufacturing; the stress–strain curve and the following fracture morphology showed brittleness with nonobvious plastic deformation, implying that the yield strength was nearly at the ultimate value. As shown in Figure 3 and Figure 4, all the specimens showed brittle failure without plastic deformation, especially in the high-porosity C4–C6 specimens. Therefore, the yield strength was adopted to obtain the stress concentration values to predict the threshold values of the failure. Figure 6c shows the statistical histogram of the stress values for C4, C5, and C6, based on the FEA results. The first loading stress of the specimens depended on their yield strengths of 80, 53, and 43 MPa for C4, C5, and C6, respectively. All the specimens showed a typical double peak distribution, representing the top surface (first peak) and the V-strut center (second peak) of the compressive specimens. The V-strut center withstood the stress of 400–600 MPa at the first level of applied stress, which was much smaller than the compressive yield strength of solid Ti6Al4V grade 23 at 860 MPa [33]. Owing to the brittleness, the failure strength was assumed as 800–900 MPa, and the applied stress was based on the yield strength at four levels. Each level was increased by 10 MPa, and the applied stress causing the local stress over the failure strength was called the threshold stress. As the applied stress increased from levels 1 to 4 in each specimen, the stress on the V-strut also increased. The threshold applied stresses were 100, 83, and 63 MPa for C4, C5, and C6, respectively, which were close to the experimental results (110, 68, and 62 MPa, respectively).
The stress-concentrated site shifting phenomenon can be explained by the geometry of these specimens. As shown in Figure 7a, these porous specimens were recognized as a series of materials with H-strut and V-strut cross-sectional areas (A1, A2). The cross-sectional area A2 of the V-strut was low, approximately 1/2 A1 for the C1 specimen. Therefore, the stress on A2 was twice that on A1. However, node A1 showed an additional stress on the edge owing to the shear-induced strain area, and this stress concentration phenomenon (k factor) has been widely studied [34]. The stress values at the strut center and node edge were obtained from Figure 6b to discuss the stress on the strut size. To simplify, the stress obtained above should be divided by the yield stress 201, 166, 112, 80, 53, and 43 MPa for C1 to C6, respectively, to obtain a relative stress ratio before comparison. Figure 7b shows the 1/strut size2 values, depending on the relative stress ratios of different sites. Both sites showed a linear relationship with the 1/strut size2, although the slopes differed. The stress values should be related to the strut areas, and the obtained FEA result is quite reasonable. The stress on the strut center was smaller than that at the node edge at the beginning (C1), but it increased rapidly with the strut size. When the strut size was between 1.125 mm (C2) and 1 mm (C3), the stress on the strut center was larger than that at the node edge. This is the reason for the changes in the stress concentrations with the changes in porosities.
Using the above two types of deformation in the low- and high-porosity specimens, the stress behavior and initial crack site using the elastic model instead of the plastic model were validated. To the best of our knowledge, a method for predicting the initial crack site of a porous material has not been reported thus far. It is assumed that the elastic model of the FEA for a brittle material after the EBAM process without heat treatment is used to calculate the stress behavior and predict the initial crack site. Compared to previous research reporting on the plastic model of FEA to predict mechanical properties [20,21,22,23], initial stiffnesses, and plastic collapse strengths [24], as well as micro stress analysis using a multiscale simulator with micro-CT, which provided an FEA to predict the strength [35], the method demonstrated in this study presents a simple way to predict the initial crack site of a porous structure. Such numerical results would contribute to the design of porosities and Young’s modulus for the studying of fracture mechanisms in future.

4. Conclusions

In this study, we coupled compression curves and images with FEA to investigate the failure mechanisms and initial crack sites of Ti6Al4V porous structures with cubic units, whose porosity ranged from 25% to 63%, fabricated by electron beam additive manufacturing. From this study, we reached the following conclusions:
  • The deformation at low porosities showed a shear diagonal pattern, and the structure was brittle with a high fluctuation peak in the plateau of the compressive stress–strain curve. With an increase in porosity, the deformation behavior shifted toward a layer-by-layer collapse and showed a uniform fluctuation in the plateau region.
  • In terms of the elastic model (FEA), in the low-porosity specimens, where the stress statics diagram showed a normal distribution, most of the stresses were smaller than the yield strengths. The maximum stress concentration was located at the edge of the strut and the node. In contrast, in the high-porosity specimens, the stress statics shifted the low and high stresses into two peak distributions, and the maximum stress concentration was located at the center of the struts.
  • The applied threshold stresses were 100, 83, and 63 MPa for C4, C5, and C6, respectively, which were close to the experimental results (110, 68, and 62 MPa).
  • By combining the compressive experimental results and FEA results, the initial failure sites of the low-porosity specimens were demonstrated to be located at the edge of the struts in the FEA results; the structure collapsed diagonally with the formation of one shear band throughout the specimen. However, the high-porosity specimens possessed maximum stresses at the center of the struts; the structures collapsed layer-by-layer.

Author Contributions

Conceptualization, I.-G.C. and H.-C.F.; methodology, M.-H.T. and C.-M.Y.; software, C.-Y.J., C.-M.Y. and Y.-J.C.; validation, M.-H.T. and Y.-X.H.; writing—original draft preparation, M.-H.T. and C.-M.Y.; writing—review and editing, Y.-J.C., I.-G.C. and H.-C.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministry of Science and Technology of Taiwan, grant number MOST 108-2218-E-992-320-MY2, 109-2622-E-992-013-CC3 and 106-2218-E-006-006.

Acknowledgments

The authors gratefully acknowledge the sponsorship from the Ministry of Science and Technology of Taiwan, under project No. MOST 108-2218-E-992-320-MY2. The authors are thankful for the EBAM lab equipment and financial support from the Metal Industry Research and Development Centre (MIRDC).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Geng, X.; Li, Y.; Li, F.; Wang, X.; Zhang, K.; Liu, Z.; Tian, H. A new 3D printing porous trabecular titanium metal acetabular cup for primary total hip arthroplasty: A minimum 2-year follow-up of 92 consecutive patients. J. Orthop. Surg. Res. 2020, 15, 383. [Google Scholar] [CrossRef] [PubMed]
  2. Murr, L.E.; Gaytan, S.M.; Martinez, E.; Medina, F.; Wicker, R.B. Next generation orthopedic implants by additive manufacturing using electron beam melting. Int. J. Biomater. 2012, 2012, 245727. [Google Scholar] [CrossRef]
  3. Dall’Ava, L.; Hothi, H.; Di Laura, A.; Henckel, J.; Hart, A. 3D Printed Acetabular Cups for Total Hip Arthroplasty: A Review Article. Metals 2019, 9, 729. [Google Scholar] [CrossRef] [Green Version]
  4. Yan, R.; Luo, D.; Huang, H.; Li, R.; Yu, N.; Liu, C.; Hu, M.; Rong, Q. Electron beam melting in the fabrication of three-dimensional mesh titanium mandibular prosthesis scaffold. Sci. Rep. 2018, 8, 750. [Google Scholar] [CrossRef] [Green Version]
  5. Hsu, T.-I.; Jhong, Y.-T.; Tsai, M.-H. Effect of Gradient Energy Density on the Microstructure and Mechanical Properties of Ti6Al4V Fabricated by Selective Electron Beam Additive Manufacture. Materials 2020, 13, 1509. [Google Scholar] [CrossRef] [Green Version]
  6. Wu, S.-H.; Li, Y.; Zhang, Y.-Q.; Li, X.-K.; Yuan, C.-F.; Hao, Y.-L.; Zhang, Z.-Y.; Guo, Z. Porous titanium-6 aluminum-4 vanadium cage has better osseointegration and less micromotion than a poly-ether-ether-ketone cage in sheep vertebral fusion. Artif. Organs 2013, 37, E191–E201. [Google Scholar] [CrossRef]
  7. Fambri, L. Book Review: Biomaterials. An Introduction, Second Edition Edited by J. B. Park and R. S. Lakes Plenum Press, New York, 1992 ISBN 0-306-43992-1, 394 pp. J. Bioact. Compat. Polym. 1993, 8, 289–290. [Google Scholar] [CrossRef]
  8. Zhang, L.C.; Klemm, D.; Eckert, J.; Hao, Y.L.; Sercombe, T.B. Manufacture by selective laser melting and mechanical behavior of a biomedical Ti–24Nb–4Zr–8Sn alloy. Scr. Mater. 2011, 65, 21–24. [Google Scholar] [CrossRef]
  9. Kokubo, T.; Yamaguchi, S. Chemical surface modification of a titanium scaffold. In Metallic Foam Bone; Elsevier: Amsterdam, The Netherlands, 2017; pp. 161–179. [Google Scholar] [CrossRef]
  10. Dewidar, M.M.; Lim, J. Properties of solid core and porous surface Ti–6Al–4V implants manufactured by powder metallurgy. J. Alloys Compd. 2008, 454, 442–446. [Google Scholar] [CrossRef]
  11. Heimann, R.B. Thermal spraying of biomaterials. Surf. Coat. Technol. 2006, 201, 2012–2019. [Google Scholar] [CrossRef]
  12. Ashby, M.F.; Medalist, R.F.M. The mechanical properties of cellular solids. Met. Mater. Trans. A 1983, 14, 1755–1769. [Google Scholar] [CrossRef]
  13. Harrysson, O.L.; Cansizoglu, O.; Marcellin-Little, D.; Cormier, D.R.; West, H.A. Direct metal fabrication of titanium implants with tailored materials and mechanical properties using electron beam melting technology. Mater. Sci. Eng. C 2008, 28, 366–373. [Google Scholar] [CrossRef]
  14. Parthasarathy, J.; Starly, B.; Raman, S.; Christensen, A. Mechanical evaluation of porous titanium (Ti6Al4V) structures with electron beam melting (EBM). J. Mech. Behav. Biomed. Mater. 2010, 3, 249–259. [Google Scholar] [CrossRef] [PubMed]
  15. Li, S.; Xu, Q.; Wang, Z.; Hou, W.; Hao, Y.; Yang, R.; Murr, L. Influence of cell shape on mechanical properties of Ti–6Al–4V meshes fabricated by electron beam melting method. Acta Biomater. 2014, 10, 4537–4547. [Google Scholar] [CrossRef] [PubMed]
  16. Van Grunsven, W.; Hernandez-Nava, E.; Reilly, G.C.; Goodall, R. Fabrication and Mechanical Characterization of Titanium Lattices with Graded Porosity. Metals 2014, 4, 401–409. [Google Scholar] [CrossRef]
  17. Choy, S.Y.; Sun, C.-N.; Leong, K.F.; Wei, J. Compressive properties of functionally graded lattice structures manufactured by selective laser melting. Mater. Des. 2017, 131, 112–120. [Google Scholar] [CrossRef]
  18. Zhang, Q.; Xie, J.; London, T.; Griffiths, D.; Bhamji, I.; Oancea, V. Estimates of the mechanical properties of laser powder bed fusion Ti-6Al-4V parts using finite element models. Mater. Des. 2019, 169, 107678. [Google Scholar] [CrossRef]
  19. Parthasarathy, J.; Starly, B.; Raman, S. A design for the additive manufacture of functionally graded porous structures with tailored mechanical properties for biomedical applications. J. Manuf. Process. 2011, 13, 160–170. [Google Scholar] [CrossRef]
  20. Mazur, M.; Leary, M.; Sun, S.; Vcelka, M.; Shidid, D.; Brandt, M. Deformation and failure behavior of Ti-6Al-4V lattice structures manufactured by selective laser melting (SLM). Int. J. Adv. Manuf. Technol. 2015, 84, 1391–1411. [Google Scholar]
  21. Gong, H.; Rafi, K.; Gu, H.; Starr, T.; Stucker, B. Analysis of defect generation in Ti–6Al–4V parts made using powder bed fusion additive manufacturing processes. Addit. Manuf. 2014, 1, 87–98. [Google Scholar] [CrossRef]
  22. Van Bael, S.; Kerckhofs, G.; Moesen, M.; Pyka, G.; Schrooten, J.; Kruth, J. 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]
  23. Lozanovski, B.; Leary, M.; Tran, P.; Shidid, D.; Qian, M.; Choong, P.; Brandt, M. Computational modelling of strut defects in SLM manufactured lattice structures. Mater. Des. 2019, 171, 107671. [Google Scholar] [CrossRef]
  24. Ushijima, K.; Cantwell, W.; Mines, R.; Tsopanos, S.; Smith, M. An investigation into the compressive properties of stainless steel micro-lattice structures. J. Sandw. Struct. Mater. 2010, 13, 303–329. [Google Scholar] [CrossRef]
  25. Smith, M.; Guan, Z.; Cantwell, W. Finite element modelling of the compressive response of lattice structures manufactured using the selective laser melting technique. Int. J. Mech. Sci. 2012, 67, 28–41. [Google Scholar] [CrossRef]
  26. Wysocki, B.; Maj, P.; Sitek, R.; Buhagiar, J.; Kurzydłowski, K.J.; Święszkowski, W. Laser and Electron Beam Additive Manufacturing Methods of Fabricating Titanium Bone Implants. Appl. Sci. 2017, 7, 657. [Google Scholar] [CrossRef]
  27. Cheng, X.; Li, S.; Murr, L.; Zhang, Z.; Hao, Y.; Yang, R.; Medina, F.; Wicker, R. Compression deformation behavior of Ti–6Al–4V alloy with cellular structures fabricated by electron beam melting. J. Mech. Behav. Biomed. Mater. 2012, 16, 153–162. [Google Scholar] [CrossRef]
  28. Choy, S.Y.; Sun, C.-N.; Leong, K.F.; Wei, J. Compressive properties of Ti-6Al-4V lattice structures fabricated by selective laser melting: Design, orientation and density. Addit. Manuf. 2017, 16, 213–224. [Google Scholar] [CrossRef]
  29. Ataee, A.; Li, Y.; Fraser, D.; Song, G.; Wen, C. Anisotropic Ti-6Al-4V gyroid scaffolds manufactured by electron beam melting (EBM) for bone implant applications. Mater. Des. 2018, 137, 345–354. [Google Scholar] [CrossRef]
  30. AlOmar, Z.; Concli, F. A Review of the Selective Laser Melting Lattice Structures and Their Numerical Models. Adv. Eng. Mater. 2020, 22. [Google Scholar] [CrossRef]
  31. Göttingen, N.v.d.z. Mathematisch-Physikalische Klasse. In Abhandlungen der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse; Weidmannsche Buchhandlung: Frankfurt, Germany, 1913; Volume 7. [Google Scholar]
  32. Sun, J.; Yang, Y.; Wang, D. Mechanical properties of a Ti6Al4V porous structure produced by selective laser melting. Mater. Des. 2013, 49, 545–552. [Google Scholar] [CrossRef]
  33. Zhang, S.; Wei, Q.; Cheng, L.; Li, S.; Shi, Y. Effects of scan line spacing on pore characteristics and mechanical properties of porous Ti6Al4V implants fabricated by selective laser melting. Mater. Des. 2014, 63, 185–193. [Google Scholar] [CrossRef]
  34. Pilkey, W.D. Stress Concentration. In Formulas for Stress, Strain, and Structural Matrices; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2004; Chapter 6; pp. 255–305. [Google Scholar] [CrossRef]
  35. Takano, N.; Fukasawa, K.; Nishiyabu, K. Structural strength prediction for porous titanium based on micro-stress concentration by micro-CT image-based multiscale simulation. Int. J. Mech. Sci. 2010, 52, 229–235. [Google Scholar] [CrossRef]
Materials 14 07467 g001 550
Figure 1. Schematic diagram for cubic unit cell pore and strut morphologies with side view of the designed Ti6Al4V porous structures from (a–f). Long axis of the vertical struts denotes that the V-strut was parallel to the compressive force direction, while the horizontal ones denote that the H-strut was vertical to the direction of the force. Specimens with different porosities based on the strut size change from 25% to 63% in porosity are shown: (a) 25%, (b) 33%, (c) 40%, (d) 55%, (e) 59%, and (f) 63%.
Figure 1. Schematic diagram for cubic unit cell pore and strut morphologies with side view of the designed Ti6Al4V porous structures from (a–f). Long axis of the vertical struts denotes that the V-strut was parallel to the compressive force direction, while the horizontal ones denote that the H-strut was vertical to the direction of the force. Specimens with different porosities based on the strut size change from 25% to 63% in porosity are shown: (a) 25%, (b) 33%, (c) 40%, (d) 55%, (e) 59%, and (f) 63%.
Materials 14 07467 g001
Materials 14 07467 g002 550
Figure 2. (a) Flowchart, meshing and boundary condition of the COMSOL FEA simulation. (b) Mesh distribution of different settings: mesh 1 (default, fine), mesh 2 (default, extra fine), and mesh 3 (self-set).
Figure 2. (a) Flowchart, meshing and boundary condition of the COMSOL FEA simulation. (b) Mesh distribution of different settings: mesh 1 (default, fine), mesh 2 (default, extra fine), and mesh 3 (self-set).
Materials 14 07467 g002
Materials 14 07467 g003 550
Figure 3. Compressive curves for Ti6Al4V porous specimens with different porosities fabricated by EBAM. The porosities of specimens: C1—25%, C2—33%, C3—40%, C4—55%, C5—59%, and C6—63% in (a) force vs. displacement and (b) engineering stress vs. engineering strain. Before and after elastic deformation with post collapse denoted, respectively, in (1)–(4) and (5)–(7) strain stress curves in the C1 and C4 specimens. Fracture surface in (c) shows the smooth surface of C1 and the ductile dimples of C6 specimen.
Figure 3. Compressive curves for Ti6Al4V porous specimens with different porosities fabricated by EBAM. The porosities of specimens: C1—25%, C2—33%, C3—40%, C4—55%, C5—59%, and C6—63% in (a) force vs. displacement and (b) engineering stress vs. engineering strain. Before and after elastic deformation with post collapse denoted, respectively, in (1)–(4) and (5)–(7) strain stress curves in the C1 and C4 specimens. Fracture surface in (c) shows the smooth surface of C1 and the ductile dimples of C6 specimen.
Materials 14 07467 g003
Materials 14 07467 g004 550
Figure 4. Images of (1)–(4) and (5)–(7) sites in the stress–strain curve in Figure 3. The possible stress concentration region observed here is denoted by red arrows, and the diagonal crack is denoted by green arrows. (a) C1 specimen in the elastic region with increased compressive strain: the crack was found at the edge at point 1 (ɛ = 0.1). Struts were over the elastic region where vertical and horizontal struts collapsed at point 2 (ɛ = 0.15) and a diagonal crack at point 3, and the following shear diagonal collapse manner is denoted by the white line at point 4 (ɛ = 0.35). (b) C4 specimen with an increased compressive strain: the initial crack is located in the center of the vertical struts at point 5 (ɛ = 0.1); bulking vertical strut deformation then induced horizontal strut collapse at points 6 and 7, leading to layer-by-layer collapse. Proposed fracture behavior for C1 specimen in (c) and for C4 specimen in (d).
Figure 4. Images of (1)–(4) and (5)–(7) sites in the stress–strain curve in Figure 3. The possible stress concentration region observed here is denoted by red arrows, and the diagonal crack is denoted by green arrows. (a) C1 specimen in the elastic region with increased compressive strain: the crack was found at the edge at point 1 (ɛ = 0.1). Struts were over the elastic region where vertical and horizontal struts collapsed at point 2 (ɛ = 0.15) and a diagonal crack at point 3, and the following shear diagonal collapse manner is denoted by the white line at point 4 (ɛ = 0.35). (b) C4 specimen with an increased compressive strain: the initial crack is located in the center of the vertical struts at point 5 (ɛ = 0.1); bulking vertical strut deformation then induced horizontal strut collapse at points 6 and 7, leading to layer-by-layer collapse. Proposed fracture behavior for C1 specimen in (c) and for C4 specimen in (d).
Materials 14 07467 g004
Materials 14 07467 g005 550
Figure 5. Relationship between the elastic modulus and porosity. (a) Compressive stress–strain curve of extracted elastic parts from Figure 3a for the calibration of FEA results; (b) comparison between experimental data and FEA results for C1 to C6 based on the Gibson–Ashby model.
Figure 5. Relationship between the elastic modulus and porosity. (a) Compressive stress–strain curve of extracted elastic parts from Figure 3a for the calibration of FEA results; (b) comparison between experimental data and FEA results for C1 to C6 based on the Gibson–Ashby model.
Materials 14 07467 g005
Materials 14 07467 g006a 550Materials 14 07467 g006b 550
Figure 6. Von Mises stress analysis of the specimens with different porosities using COMSOL software. Side view of the stress distribution for (a) whole specimens and (b) cells; most stress was concentrated on the vertical struts in all the specimens. X-shaped stress concentration regions were found in C1 and C2 specimens. Stress concentration moved to the center of the struts in C4–C6 specimens. C3 specimen demonstrated a porosity transition. (c) Statistical histogram of stress values for the C4–C6 specimens derived from the FEA results. First peak and second peak respectively denote the top surface and strut of the compressive specimen. Applied stress was set on the basis of the yield strength from the experimental values with an increase of 10 MPa for each level.
Figure 6. Von Mises stress analysis of the specimens with different porosities using COMSOL software. Side view of the stress distribution for (a) whole specimens and (b) cells; most stress was concentrated on the vertical struts in all the specimens. X-shaped stress concentration regions were found in C1 and C2 specimens. Stress concentration moved to the center of the struts in C4–C6 specimens. C3 specimen demonstrated a porosity transition. (c) Statistical histogram of stress values for the C4–C6 specimens derived from the FEA results. First peak and second peak respectively denote the top surface and strut of the compressive specimen. Applied stress was set on the basis of the yield strength from the experimental values with an increase of 10 MPa for each level.
Materials 14 07467 g006aMaterials 14 07467 g006b
Materials 14 07467 g007 550
Figure 7. Proposed mechanism for stress concentration site shift. (a) Cross-sectional areas of H-strut and V-strut were named C1A1, C1A2 in C1 and C4A1, C4A2 in C4. (b) Relationship between the relative stress and sites of nodes and struts in all specimens. Junctions of the struts are named nodes. Relative stress ratio is the ratio between the yield strength for a particular porosity specimen and the stress of the particular site.
Figure 7. Proposed mechanism for stress concentration site shift. (a) Cross-sectional areas of H-strut and V-strut were named C1A1, C1A2 in C1 and C4A1, C4A2 in C4. (b) Relationship between the relative stress and sites of nodes and struts in all specimens. Junctions of the struts are named nodes. Relative stress ratio is the ratio between the yield strength for a particular porosity specimen and the stress of the particular site.
Materials 14 07467 g007
Table
Table 1. Structure identification and dimensions with a fixed unit cell of 4 mm of the specimens tested. D: pore diameter. S: strut diameter. DPW: designed porosity of the specimen.
Table 1. Structure identification and dimensions with a fixed unit cell of 4 mm of the specimens tested. D: pore diameter. S: strut diameter. DPW: designed porosity of the specimen.
D (mm)S (mm)DPW (%)
C11.501.2525
C21.751.1333
C32.001.0040
C42.500.7555
C52.650.6859
C62.800.6063
Table
Table 2. Mesh information and calculation stress of the different mesh conditions.
Table 2. Mesh information and calculation stress of the different mesh conditions.
Mesh ConditionElement Size
(mm)		Mesh Results
(Count)		Calculation Stress
(MPa)
MaximaMinimumTetrahedronNodeStrut CenterNode
Mesh 1 (default, fine)2.40.347,56515,065488371
Mesh 2 (default, extra fine)1.150.45151,90941,547484383
Mesh 3 (self-set)0.50.21,341,348277,728488378
Table
Table 3. Calibrated Young’s modulus E2 of the FEA using the corrected E0 = 17.92 GPa compared to the Young’s moduli obtained experimentally and the non-calibrated E1 using the original E0 = 112 GPa.
Table 3. Calibrated Young’s modulus E2 of the FEA using the corrected E0 = 17.92 GPa compared to the Young’s moduli obtained experimentally and the non-calibrated E1 using the original E0 = 112 GPa.
SpecimenC1C2C3C4C5C6
E—Experimental
(GPa)		9.46 ± 1.088.93 ± 0.885.67± 0.393.29 ± 0.462.73 ± 0.252.39 ± 0.33
E1—FEA (GPa)62.7950.8139.9122.1717.913.97
E2—FEA (GPa)10.058.136.393.552.862.24
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tsai, M.-H.; Yang, C.-M.; Hung, Y.-X.; Jheng, C.-Y.; Chen, Y.-J.; Fu, H.-C.; Chen, I.-G. Finite Element Analysis on Initial Crack Site of Porous Structure Fabricated by Electron Beam Additive Manufacturing. Materials 2021, 14, 7467. https://doi.org/10.3390/ma14237467

AMA Style

Tsai M-H, Yang C-M, Hung Y-X, Jheng C-Y, Chen Y-J, Fu H-C, Chen I-G. Finite Element Analysis on Initial Crack Site of Porous Structure Fabricated by Electron Beam Additive Manufacturing. Materials. 2021; 14(23):7467. https://doi.org/10.3390/ma14237467

Chicago/Turabian Style

Tsai, Meng-Hsiu, Chia-Ming Yang, Yu-Xuan Hung, Chao-Yong Jheng, Yen-Ju Chen, Ho-Chung Fu, and In-Gann Chen. 2021. "Finite Element Analysis on Initial Crack Site of Porous Structure Fabricated by Electron Beam Additive Manufacturing" Materials 14, no. 23: 7467. https://doi.org/10.3390/ma14237467

APA Style

Tsai, M. -H., Yang, C. -M., Hung, Y. -X., Jheng, C. -Y., Chen, Y. -J., Fu, H. -C., & Chen, I. -G. (2021). Finite Element Analysis on Initial Crack Site of Porous Structure Fabricated by Electron Beam Additive Manufacturing. Materials, 14(23), 7467. https://doi.org/10.3390/ma14237467

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop