Next Article in Journal
Thermal Wave Scattering and Temperature Concentration around the Opening in Platinum–Rhodium Leaky Plates
Next Article in Special Issue
Ultrasonic and Conventional Fatigue Endurance of Aeronautical Aluminum Alloy 7075-T6, with Artificial and Induced Pre-Corrosion
Previous Article in Journal
Computing the Growth of Small Cracks in the Assist Round Robin Helicopter Challenge
Previous Article in Special Issue
Ageing of Al-Mn-Cu-Be Alloys for Stimulating Precipitation of Icosahedral Quasicrystals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Fatigue Analysis of Auxetic Cellular Structures Made of SLM AlSi10Mg Alloy

University of Maribor, Faculty of Mechanical Engineering, Smetanova 17, 2000 Maribor, Slovenia
*
Author to whom correspondence should be addressed.
Metals 2020, 10(7), 945; https://doi.org/10.3390/met10070945
Submission received: 26 May 2020 / Revised: 6 July 2020 / Accepted: 10 July 2020 / Published: 14 July 2020
(This article belongs to the Special Issue Structure and Properties of Aluminium Alloys)

Abstract

:
In this study, a computational fatigue analysis of topology optimised auxetic cellular structures made of Selective Laser Melting (SLM) AlSi10Mg alloy is presented. Structures were selected from the Pareto front obtained by the multi-objective optimisation. Five structures with different negative Poisson’s ratios were considered for the parametric numerical analysis, where the fillet radius of cellular struts has been chosen as a parameter. The fatigue life of the analysed structures was determined by the strain–life approach using the Universal Slope method, where the needed material parameters were determined according to the experimental results obtained by quasi-static unidirectional tensile tests. The obtained computational results have shown that generally less auxetic structures tend to have a better fatigue life expectancy. Furthermore, the fillet radius has a significant impact on fatigue life. In general, the fatigue life decreases for smaller fillet radiuses (less than 0.3 mm) as a consequence of the high-stress concentrations, and also for larger fillet radiuses (more than 0.6 mm) due to the moving of the plastic zone away from the edge of the cell connections. The obtained computational results serve as a basis for further investigation, which should be focused on the experimental testing of the fabricated auxetic cellular structures made of SLM AlSi10Mg alloy under cyclic loading conditions.

1. Introduction

Additive Manufacturing (AM) is a fabrication process that provides unique opportunities to manufacture customised parts with complex geometries or functionally graded materials [1,2]. As presented in [3,4], AM brings the promise of a new industrial revolution, offering the capability to manufacture parts through the repetitious deposition of material layers directly from a digital Computer-Aided Design (CAD) model. For those reasons, AM applies to a wide range of industries [5,6], including medical equipment [7,8], automotive [9,10], aerospace [11,12], and other consumer products [13]. Many AM technologies are being currently used in different engineering applications [14,15], e.g., Direct Metal Laser Sintering (DMLS), Selective Laser Melting (SLM), Electron Beam Melting (EBM), and Laser Metal Deposition (LMD). Among the AM processes, DMLS and SLM techniques are particularly widespread, especially for the aluminium alloy processing [16,17]. The DMLS technology is appropriate to build parts out of any metal alloy, while SLM can only be used with certain metals.
For AM technologies, Al–Si or Al–Si–Mg are the common alloys. In general, the static strength of AM parts depends on the density of the parts, as well as on the microstructure formed during the fabrication. As presented by Aboulkhair et al. [18] and Kempen et al. [19], the microstructure of the AM fabricated parts made of the AlSi10Mg alloy is usually finer (higher static strength) compared to the parts, which are fabricated via classical procedures (e.g., casting). On the other hand, the microstructure of the AM fabricated parts is generally anisotropic and depends on the building and scanning directions during the AM process. For that reason, the monotonic material properties (yield stress, ultimate tensile strength) are also anisotropic and may strongly depend on the orientation of texture [20,21,22]. Zhou et al. [23] have shown that the microstructure, hardness, and static strength of the AM parts may be significantly influenced by additional heat treatment. Nevertheless, the fatigue behaviour of the AM parts is not well understood yet, although a lot of research was performed in the last decade [24]. Similar to the static mechanical properties, the fatigue strength of the AM parts primarily depends on their microstructure. However, the AM process usually leads to an increased surface roughness of the fabricated part, which decreases its fatigue strength [25]. Furthermore, material defects such as the porosity and insufficient layer bonding could result in an increased scatter of the fatigue properties [26]. Besides the fatigue analyses, the investigations on the fracture behaviour of the AM components have also been carried out, especially due to the high importance in the automotive industry and aerospace applications [27,28,29].
As mentioned above, AM technologies, including the SLM process, provide opportunities to manufacture very complex geometries or functionally graded materials, which are inaccessible through traditional manufacturing techniques. Typical lightweight structures, where the SLM process may be used, are pre-designed advanced cellular metamaterials and structures, which are difficult to fabricate by the traditional processes. In the last period, the research has been focusing on the mechanical response of the gradient lattice structures manufactured via the AM process in terms of energy absorption under the compressive loading conditions [30,31,32]. Some authors also investigated the fatigue behaviour of different AM lattice structures considering the effect of the cellular topology and geometrical imperfections on the fatigue life of the analysed structures [33,34,35]. The next group of metamaterials suitable to be produced using the AM technology are the auxetic cellular structures, which are the subject of this study (for a detailed description of the auxetic cellular structures, see Section 2.1).
The research work in the presented study is focused on the fatigue behaviour of five different auxetic cellular structures obtained by the topology optimisation and made of the SLM AlSi10Mg alloy, which is widely used in the aerospace and automotive industry because of its high specific strength, high corrosion resistance, and good flowability [36,37]. The fatigue life of different auxetic cellular geometries, also considering the influence of the fillet radiuses, is determined by the strain–life approach using Universal Slope method [38]. The material parameters were determined according to the experimental results previously obtained by the unidirectional quasi-static tests [39]. Numerical nonlinear finite element simulations were performed to determine the deformation behaviour and the stress and strain fields in the critical cross-sections of the analysed cellular structures. Then, the numerical results were used as the input parameters in the subsequent fatigue analyses.

2. Selection of Optimised Auxetic Cellular Structures

2.1. Auxetic Cellular Structures

Cellular structures have been increasingly used in modern engineering applications over the past decades due to their attractive combination of mechanical and thermal properties [40,41,42]. Recent advances in the AM technologies enable the fabrication of parts with complex internal cellular structure, which are commonly adjusted for specific applications. The auxetic cellular structures exhibit a negative Poisson’s ratio and thus experience large volume changes under loading [43,44,45]. The auxetic structures tend to expand in the lateral direction when subjected to the uniaxial tensile loading and vice versa in the case of the compressive loading (Figure 1).
Such behaviour is made possible by their 2D or 3D designed hinge-like cellular skeleton with predefined geometry [46,47,48]. The auxetic cellular structures (made of metals [49], polymers [50,51], and textiles [52,53,54]) offer extraordinary mechanical properties (under quasi-static [55], impact [56,57], and fatigue [58,59] loading), i.e., variable stiffness [60], low density, large mechanical energy absorption through deformation [61,62,63], and unique deformation behaviour [64,65]. Their cell shape changes rapidly through the complete auxetic structure under loading; thus, the loading efficiently spreads through the structure. The impact on a small part of the auxetic structure results in energy dissipation through its entire structure. With recent advances in the AM technologies [66], it is possible to fabricate the auxetic cellular structures with optimised and exactly predefined geometries (i.e., strut thickness, geometry, cell shape and size). The design and optimisation of auxetic structures can be done by using advanced Computer-Aided Engineering (CAE) methodologies to achieve and control the desired mechanical properties for specific applications [60,67]. Such precise fabrication control of the shape, size, and distribution of cells with the use of layered additive technologies makes the auxetic structures superior to the other conventionally manufactured cellular materials.

2.2. Topology Optimised Auxetic Structures

Many auxetic cellular structures have been developed up to date [45]. Their shape, topology, and dimensions were usually custom-made for the task for which they were designed. Therefore, it is difficult to compare the fatigue properties between them. In the previous study by Borovinšek et al. [67], a procedure for the determination of optimal auxetic topology was developed using a numerical approach. The auxetic topology was modelled as a quarter cell of a 2D planar periodic structure. Appropriate boundary conditions were applied to account for the rest of the structure. The quarter cell model had a square shape, and its topology was represented by a grid of 10 × 10 equal squares (parameter regions), each of which could either represent material or a void region (Figure 2a). In this way, any topology of the quarter cell could be obtained by changing the content of the parameter regions.
The performance of the auxetic structure was regarded as a performance of the compliant mechanism, and as such, two objective functions were applied to determine the optimal topologies. The first objective function was used to measure the stiffness of the auxetic structure, while the second objective function measured its flexibility. Using two objective functions and a multi-objective optimisation procedure, which changed the values of the parameter regions, resulted in a set of best topologies called the Pareto front [67]. All solutions on the Pareto front are optimal and equally good. Thus, applying this procedure results in a set of topologies, which have an auxetic behaviour. Since they are optimised with the same procedure and with the same objectives, it is also possible to directly compare their fatigue behaviour.
The drawback of this procedure is that the resulting optimal topologies of the quarter cell models are sharp-edged, since their topology consists of square-shaped parameter regions (Figure 2a). Sharp edges of the structure result in stress concentrations under load and do not allow a direct application of the obtained topologies from [67]. To study the fatigue behaviour of the obtained topologies, a geometry smoothing procedure was developed for this investigation. It was applied to reduce the number of possible stress concentration regions.
The smoothing procedure consists of the following five steps:
  • Selection of the topology from the Pareto front (Figure 2a, dark grey area).
  • Determination of the medial axis of the topology.
  • Simplification of the medial axis to straight lines between the line connection points (Figure 2b, red lines).
  • Thickening of the simplified medial axis to obtain the geometry; the horizontal and vertical lines retained their thickness, while a new thickness was defined to the oblique lines of the medial axis. The thickness of these lines was set to such a value that the topology area remained the same (Figure 2b, light grey area).
  • All remaining sharp corners were rounded using an equal fillet radius R (Figure 2c).
From the Pareto front, obtained by the multi-objective optimisation, which contained 200 different topologies, five base topologies were selected for this study. The topologies of the Pareto front were first sorted by their Poisson’s ratio starting with the lowest one (Topology 1), followed by the median one (Topology 3), and so on until the highest one (Topology 5). An additional two topologies were selected between the lowest and the median Poisson’s ratio (Topology 2) and between the median and the highest Poisson’s ratio (Topology 4). The resulting five topologies are almost equally apart from each other on the Pareto front; thus, they are a good representation of the best solution set. The Poisson’s ratio and the mass of all five base topologies are reported in Chapter 3. After selection, the smoothing procedure was applied, which resulted in five distinct quarter unit cell topologies (Figure 3). Their complete auxetic structures were made by mirroring and multiplying the quarter cell topology in such a way that the resulting periodic structures consisted of six-unit cells in both the (x and y) coordinate directions.

3. Material Characterisation and Fatigue Analyses

3.1. Material Characterisation

The material used in this study is aluminium alloy AlSi10Mg. The specimens used for the further tensile tests have been fabricated by the SLM technology using the EOSINT-M-270 system without additional polishing. The general process data and the material composition are shown in Table 1. Figure 4a shows the geometry of the flat specimen, while the designated building and scanning fabrication directions are shown in Figure 4b.
Figure 5 shows the polished surface of the fabricated specimen. The black circles represent the pores. They are significantly more frequent and larger near the specimen surface all around the circumference, up to 500 μm inwards. These near-surface pores were formed at the beginning/end of a scan in the x-direction and were not removed by the laser movement in the y-direction during the formation of the subsequent layer.
Two specimens were used for the quasi-static tensile tests to measure the stress–strain (σ–ε) response. The tensile tests were carried out using the 100 kN MTS Landmark hydraulic machine (MTS Systems GmbH, Germany) at a room temperature of 23 °C. The tests were displacement controlled with a loading rate of 0.5 mm/min according to the DIN-ISO 6892 standard. The force was measured with the 100 kN MTS load cell, and the strains were measured with the MTS 834.11F-24 extensometer (MTS Systems GmbH, Germany). The measured quasi-static stress–strain responses [39] were almost the same for both specimens and are in a good agreement if compared to the experimental results from other researchers [68]. Based on the measured stress–strain responses (Figure 6), the average material parameters have been determined (Table 2).
The true stress–strain (σtrue–εtrue) data were determined using engineering stress (σeng) and engineering strain (εeng) values as follows:
ε true = ln ( 1 + ε eng )
σ true = σ eng × ( 1 + ε eng ) .

3.2. Strength Analysis

The numerical models of five auxetic structures were built in Ansys Mechanical [69], using a 2D plane stress formulation. The size of 10 mm × 10 mm was chosen for the one-unit cell, which results in a structure size of 60 mm × 60 mm and a nominal strut thickness of 0.5 mm. The finite element mesh consisted of plane stress parabolic quad finite elements PLANE183 with eight nodes. With the parametrical analysis, it was deduced that the results converge using a mesh, where the global element size is approximately 0.2 mm, with a consideration of curvatures, where the element size is even smaller. This resulted in numerical models with an average of 30,000 finite elements and 100,000 nodes. Each topology was constrained in the y-direction at the bottom and loaded with displacement at the top. All five structures were loaded with the same displacement of 0.4 mm in the + y-direction, thus creating a tension in the vertical direction. Such loading was chosen to induce plasticity in all structures. The higher loading displacement would cause severe plasticity in some structures. Similarly, a lower value of displacement would not cause plastic strain in some structures. The models were constrained in the x-direction in the middle of the structure to prevent horizontal movement. The boundary conditions are presented in Figure 7, where loading is acting in the vertical direction (y-direction).
The material in all presented numerical models was the SLM AlSi10Mg alloy with the properties presented in Section 3.1. The models were loaded with the displacement, which caused plastic deformation in all five models. Therefore, material nonlinearity was considered in the simulations. For plasticity, a multilinear isotropic hardening model was used, with the utilisation of the true stress–strain material data shown in Figure 6.
Numerical results in the form of total displacements are presented in Figure 8. The resulting displacement was scaled by the factor 11 for a better visualisation of all structures. It can be noted that they have different Poisson ratios. The Poisson ratios presented in Table 3 were calculated as the ratio between the longitudinal and lateral displacement. The Poisson ratios of the analysed structures are, to some extent, different than the unit cell’s Poisson ratios. All the calculated Poisson ratios are negative, which proves that all the analysed structures are also auxetic. It can be seen that Topology 1 has the lowest Poisson ratio, while Topology 5 has the highest Poisson ratio. The Poisson ratios of the other three topologies are ascending between Topology 1 and Topology 5. Table 3 also reports the mass of each model showing that all structures have a similar weight. Only Topology 3 is a bit more massive.
Figure 9 shows the location of the maximum plastic strain for each model. The area of the maximum plastic strain is also magnified to visualise the exact location of the developed plastic zones. Location of the maximum plastic strain is in all cases at the fillet radius of the auxetic structure.
Since the size of the radius is important, a series of parametric analyses were conducted with different radius sizes. Radiuses were scaled in a range of 0.1 and 0.9 mm with a step of 0.1 mm. Figure 3 shows the basic topology designs from which the numerical models were built. Most of the structures allow a wide range of radiuses, while Topology 1 and Topology 5 are limited because the two fillets intersect and prevent further radius scaling. Therefore, some radiuses were scaled, while the potentially intersecting radiuses were presented only up to the maximum possible values, e.g., Topologies 1 and 5. Each time the radius was scaled, the numerical model was regenerated with the same mesh parameters and the same boundary conditions as presented in Figure 7. In that way, it was possible to determine how radius sensitive the analysed structures are and the optimal radius for each structure in terms of fatigue.
Figure 10 shows the equivalent plastic strain for different fillet radiuses. For Topology 1, some fillets could not be scaled beyond 0.25 mm, while the remaining fillets are limited to 0.75 mm; therefore, the results for the radius of 0.9 mm for Topology 1 are missing.
From the results in Figure 10, it is clear that at a small fillet radius, the plastic strain occurs near the edge, which is placed where the crack initialisation is expected. With a larger fillet radius, the area of the plastic zone is wider, and the value of the equivalent plastic strain is lower in all structures. In most of the observed structures, this occurs at the radius of approximately 0.4 mm, with exceptions in Topology 1 (at the radius of approximately 0.2 mm) and Topology 4 (at the radius of approximately 0.3 mm). With further increase of the fillet radius, the value of the equivalent plastic strain increases. Furthermore, some new critical areas occur. Accordingly, large fillet radiuses create new critical places for the crack initiation, where sometimes the equivalent plastic strain has the same or larger value than observed in the case of very small fillet radiuses.
Then, the obtained numerical results were used for fatigue life determination. From each numerical result, the maximum equivalent elastic strain was added to the maximum equivalent plastic strain. This allows the calculation of the fatigue life cycle for each topology, including the fillet radius variations (Figure 10). For each structure, additional numerical analysis was performed to determine the yield load condition, when the structure behaves purely elastic.

3.3. Fatigue Life Determination

In the presented study, the strain–life approach is used to determine the fatigue life of treated cellular structures. In general, this approach is based on the knowledge of the stresses and strains that occur at locations, where the fatigue crack nucleation is likely to start, e.g., holes, fillets, grooves. In this approach, the fatigue behaviour of the material is described by the strain–life curve expressed by the Coffin–Manson relation between fatigue life Ni and the total amplitude strain [70]. However, the material low-cycle fatigue parameters should be determined previously using the appropriate low-cycle fatigue test. If it is not the case, an alternative method based on the material parameters obtained from a quasi-static test can be used. In this study, the method of Universal Slopes as proposed by Muralidharan and Manson [38] has been used to obtain the number of loading cycles Ni required for the fatigue crack initiation in the critical cross-section of the treated cellular structure:
ε E a = 0.623 ( R m E ) 0.832 ( 2 N i ) 0.09 + 0.0196 ( ε f ) 0.155 ( R m E ) 0.53 ( 2 N i ) 0.56
where ε E a is the equivalent amplitude strain, E is the elasticity modulus, Rm is the ultimate tensile strength, and ε f is the true fracture strain. Once the equivalent amplitude strain ε E a is known, the number of stress cycles Ni can be obtained iteratively from Equation (3) using Matlab software [71] and with a consideration of the material parameters E, Rm, and ε f given in Table 2. The results are shown in Figure 11.
For Topology 1, the expected fatigue life cycles are rather low, with a maximum value at the fillet radius of 0.2 mm. Topology 4 resulted in the best fatigue life expectance with the fillet radius of 0.3 mm. All other structures have the highest number of cycles with a radius of 0.4 mm, while the life expectancy of Topology 3 is very low. The best fatigue life is expected for Topology 5, which also has the highest Poisson ratio, making it the least auxetic structure presented in the paper. On the other hand, Topology 2 has the second-highest fatigue life expectancy, while having the second-lowest Poisson ratio. Therefore, Topology 2 could be used in cases where a highly auxetic and durable structure is needed.

4. Discussion

Five cellular structures were considered for the numerical analysis to identify the most suitable auxetic topology when considering the fatigue life under tension loading conditions. Structures were selected from the Pareto front obtained by the multi-objective optimisation. The structures were sorted by their Poisson ratio, starting with the lowest one (Topology 1) to the one that has the Poisson’s ratio only slightly below zero (Topology 5), as can be seen from Table 3.
All five topologies were subjected to the same tensile loading conditions. The amount of the tensile load was high enough to cause plasticity in every topology but small enough that each topology would sustain several loading cycles until failure. The presented results of the numerical analysis show that generally less auxetic structures tend to have a better fatigue life expectancy.
The fillet radius has a significant impact on fatigue life. While it is expected that smaller radiuses would increase the probability of crack initiation, it is interesting to note that larger fillet radiuses also negatively influence the fatigue life. As the fillet radius increases, the plastic zone is moved away from the edge, as can be seen from Figure 10. Global tension loading displacement induces hinge-like behaviour in structure corners. While the hinges rotate, the structures struts are primarily loaded in bending. Since the struts have a constant cross-section area, the internal bending moment is equally distributed along their length, and the whole strut deforms to account for the hinge rotation. With increased fillet radius, the strut length is decreased. To account for the same hinge rotation, a shorter strut must deform more, which causes a higher rate of plasticity with the same tension loading displacement.
Topology 1 has the lowest (most negative) Poisson ratio, which makes it the most auxetic structure analysed herein. Computed displacements in the lateral direction (x-direction) are more than two times larger than the loading displacements in the longitudinal direction (y-direction). However, the prescribed tension displacement is quite large for this structure, as it results in an equivalent plastic strain of 1.5%. The fatigue life under such loading conditions would be very short, which is also clearly shown in Figure 11. Loading would have to be reduced at least 10 times to achieve the pure elastic behaviour of the structure, which questions its usability compared to the other analysed structures.
As can be seen from Figure 11, Topology 2 is the second most auxetic structure and has the second-best fatigue life. This makes Topology 2 the best choice for an application where the auxetic response is needed. With the Poisson ratio of −0.95, the displacement in the lateral direction (x-direction) is almost identical to the tension loading displacement. With the applied loading, the structure would sustain approximately 2500 cycles, and if the loading displacement would be decreased to less than half of the applied loading, the structure would behave purely elastically, while still exhibiting the same displacement in both directions. This makes Topology 2 the most appropriate choice for the auxetic structure.
By analysing Topology 3, it can be observed that the lateral displacement is nearly half the size of the tension loading displacement. However, the equivalent plastic strain is approximately 1%, and the structure has the second-worst fatigue life expectancy. Additionally, Topology 3 is the most massive structure compared to the other topologies. Therefore, this structure is not recommended for an application subjected to cyclic loading unless a specific need utilising such a Poisson ratio is required.
The Poisson ratio of Topology 4 is equal to −0.2. Topology 4 also has a high expected fatigue life. Similar to Topology 2, this structure would behave purely elastically when the loading displacement is reduced by half. Comparing Poisson ratio, it is evident that the lateral displacement would be five times lower, as it was observed in the case of Topology 2. It could be recommended for applications where a small amount of auxetic behaviour is required. Nevertheless, Topology 2 shows a better performance in terms of lateral displacement and fatigue life.
Topology 5 shows the smallest auxetic effect, as its Poisson ratio is close to zero. At the same time, it is the structure with the highest number of fatigue life cycles determined, as can be seen from Figure 11. If the tension loading displacement would be reduced only by 25%, the structure would behave purely elastic. Therefore, Topology 5 could be used for applications where the transversal strain should be limited.

5. Conclusions

In this study, the fatigue behaviour of auxetic cellular structures obtained by multi-objective topology optimisation and made of SLM AlSi10Mg alloy is presented. Five structures with different Poisson ratios were considered for the parametric numerical analysis, where the fillet radius of cellular struts has been chosen as a parameter. The fatigue life of the analysed structures was determined by the strain–life approach using the Universal Slope method, where the needed material parameters were determined according to the experimental results obtained by the quasi-static unidirectional tensile tests. Based on the fractography of the quasi-static specimens, numerical simulations, and fatigue analyses, the following conclusions can be made:
  • The microstructure of treated SLM AlSi10Mg alloy consists of bands, which were produced by varying the scan direction in each subsequent layer. The material porosity is more frequent and larger in the near-surface layers, where pores were formed at the end/beginning of scanning during the SLM process.
  • The obtained material parameters by the quasi-static tests (proportional limit, ultimate tensile strength, strain at rupture) are comparable to the results presented by the other researchers.
  • The obtained computational results have shown that less auxetic structures (higher Poisson ratios) tend to have a better fatigue life expectancy (the longest fatigue life has been obtained for the structure topology with the minimum auxetic characteristics).
  • The fillet radius of cellular struts has a significant impact on fatigue life. Computational analyses have shown that the fatigue life decreases for smaller fillet radiuses (less than 0.3 mm) as a consequence of high-stress concentrations and also for larger fillet radiuses (more than 0.6 mm) due to the movement of the plastic zone away from the edge of the cell connections. Besides the fillet radius sizes, other reasons for stress concentrations might exist. This should be investigated in future research.
  • The fatigue life in this study was obtained using the simplified Universal Slope method, where the material parameters were determined by quasi-static tensile tests considering the partly porous structure of analysed SLM AlSi10Mg alloy (Figure 5). For a more accurate determination of fatigue life of real components made of SLM AlSi10Mg alloy, the comprehensive Low Cycle Fatigue (LCF) tests should be performed to determine the appropriate LCF material properties, where the influence of pores will be considered in detail. Once the LCF material parameters are known, the standardised strain life approach could be used for the subsequent fatigue analysis.
  • The obtained computational results provide a basis for further investigation including experimental testing of the fabricated auxetic cellular structures made of SLM AlSi10Mg alloy under the cyclic loading conditions.

Author Contributions

Conceptualization, M.U. and S.G.; investigation, S.G.; methodology, M.V.; project administration, M.V.; software, M.B.; Writing—review and editing, M.U., M.B., M.V. and S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research Agency—ARRS, basic research project No. J2-8186 and research core No. P2-0063.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Savolainen, J.; Collan, M. How additive manufacturing technology changes business models? Review of literature. Addit. Manuf. 2020, 32, 101070. [Google Scholar] [CrossRef]
  2. Plocher, J.; Panesar, A. Review on design and structural optimisation in additive manufacturing: Towards next-generation lightweight structures. Mater. Des. 2019, 183, 108164. [Google Scholar] [CrossRef]
  3. Chunyong, L.; Yazhou, H.; Ning, L.; Xianrui, Z.; Hongshui, W.; Xinping, Z.; Yulan, F.; Jingyun, H. Laser polishing of Ti6Al4V fabricated by selective laser melting. Metals 2020, 10, 191. [Google Scholar] [CrossRef] [Green Version]
  4. Mugwagwa, L.; Yadroitsev, I.; Matope, S. Effect of process parameters on residual stresses, distortions, and porosity in selective laser melting of maraging steel 300. Metals 2019, 9, 1042. [Google Scholar] [CrossRef] [Green Version]
  5. Attaran, M. The rise of 3-D printing: The advantages of additive manufacturing over traditional manufacturing. Bus. Horiz. 2017, 60, 677. [Google Scholar] [CrossRef]
  6. Gao, W.; Zhang, Y.; Ramanujan, D.; Ramani, K.; Chen, Y.; Williams, C.B.; Wang, C.L.; Shin, Y.C.; Zhang, S.; Zavattieri, P.D. The status, challenges, and future of additive manufacturing in engineering. Comput. Aided Des. 2015, 69, 65–89. [Google Scholar] [CrossRef]
  7. Jardini, A.L.; Larosa, M.A.; Filho, R.M.; Zavaglia, C.A.C.; Bernardes, L.F.; Lambert, C.S.; Calderoni, D.R.; Kharmandayan, P. Cranial reconstruction: 3D biomodel and custom-built implant created using additive manufacturing. J. Cranio-Maxillofac. Surg. 2014, 42, 1877–1884. [Google Scholar] [CrossRef]
  8. Mahmound, D.; Elbestawi, M. Lattice structures and functionally graded materials applications in additive manufacturing of orthopedic implants: A review. J. Manuf. Mater. Process. 2017, 1, 13. [Google Scholar] [CrossRef]
  9. Chuang, C.H.; Chen, S.; Yang, R.J.; Vogiatzis, P. Topology optimisation with additive manufacturing consideration for vehicle load path development. Int. J. Numer. Methods Eng. 2017, 1434–1445. [Google Scholar] [CrossRef]
  10. Abdi, M.; Ashcroft, I.; Wildman, R.D. Design optimisation for an additively manufactured automotive component. Int. J. Powertrains 2018, 7, 142–161. [Google Scholar] [CrossRef]
  11. Seabra, M.; Azevedo, J.; Araújo, A.; Reis, L.; Pinto, E.; Alves, N.; Santos, R.; Mortagua, J.P. Selective laser melting (SLM) and topology optimisation for lighter aerospace components. Procedia Struct. Integr. 2016, 1, 289–296. [Google Scholar] [CrossRef] [Green Version]
  12. Zhu, J.H.; Zhang, W.H.; Xia, L. Topology optimisation in aircraft and aerospace structures design. Arch. Comput. Methods Eng. 2016, 23, 595–622. [Google Scholar] [CrossRef]
  13. Bader, C.; Patrick, W.G.; Kolb, D.; Hays, S.G.; Keating, S.; Sharma, S.; Dikovsky, D.; Belocon, B.; Weaver, J.C.; Silver, P.A.; et al. Grown, printed, and biologically augmented: An additively manufactured microfluidic wearable, functionally templated for synthetic microbes, 3D print. Arch. Addit. Manuf. 2016, 3, 79–89. [Google Scholar] [CrossRef] [Green Version]
  14. Simar, A.; Godet, S.; Watkons, T.R. Highlights of the special issue on metal additive manufacturing. Mater. Charact. 2018, 143, 1–4. [Google Scholar] [CrossRef]
  15. Herzog, D.; Seyda, V.; Wycisk, E.; Emmelmann, C. Additive manufacturing of metals. Acta Mater. 2016, 117, 371–392. [Google Scholar] [CrossRef]
  16. Ismael, A.; Wang, C.S. Effect of Nb additions on microstructure and properties of γ-TiAl based alloys fabricated by selective laser melting. Trans. Nonferrous Metals Soc. China 2019, 29, 1007–1016. [Google Scholar] [CrossRef]
  17. Bagherifard, S.; Beretta, N.; Monti, S.; Riccio, M.; Bandini, M.; Guagliano, M. On the fatigue strength enhancement of additive manufactured AlSi10Mg parts by mechanical and thermal post-processing. Mater. Des. 2018, 145, 28–41. [Google Scholar] [CrossRef]
  18. Aboulkhair, N.T.; Tuck, C.; Ashcroft, I.; Maskery, I.; Everitt, N.M. On the precipitation hardening of selective laser melted AlSi10Mg. Metall. Mater. Trans. A 2015, 46, 3337–3341. [Google Scholar] [CrossRef]
  19. Kempen, K.; Thijs, L.; Van Humbeeck, J.; Kruth, J.P. Mechanical properties of AlSi10Mg produced by selective laser melting. Phys. Procedia 2012, 39, 439–446. [Google Scholar] [CrossRef] [Green Version]
  20. Zhonghua, L.; Zezhou, K.; Peikang, B.; Yunfei, N.; Guang, F.; Wenpeng, L.; Shuai, Y. Microstructure and tensile properties of AlSi10Mg alloy manufactured by multi-laser beam selective laser melting (SLM). Metals 2019, 9, 1337. [Google Scholar] [CrossRef] [Green Version]
  21. Fousováa, M.; Dvorskýa, D.; Michalcováa, A.; Vojtěcha, D. Changes in the microstructure and mechanical properties of additively manufactured AlSi10Mg alloy after exposure to elevated temperatures. Mater. Charact. 2018, 137, 119–126. [Google Scholar] [CrossRef]
  22. Takata, N.; Kodaira, H.; Suzuki, A.; Kobashi, M. Size dependence of microstructure of AlSi10Mg alloy fabricated by selective laser melting. Mater. Charact. 2018, 143, 18–26. [Google Scholar] [CrossRef]
  23. Zhou, L.; Mehta, A.; Schulz, E.; Mcwilliams, B.; Cho, K.; Sohn, Y. Microstructure, recipitates and hardness of selectively laser melted AlSi10Mg alloy before and after heat treatment. Mater. Charact. 2018, 143, 5–17. [Google Scholar] [CrossRef]
  24. Daniewicz, S.R.; Shamsai, N. An introduction to the fatigue and fracture behaviour of additive manufactured parts. Int. J. Fatigue 2017, 94, 167. [Google Scholar] [CrossRef]
  25. Beevers, E.; Brandão, A.D.; Gumpinger, J.; Gschweitl, M.; Seyfert, C.; Hofbauer, P.; Rohr, T.; Ghidini, T. Fatigue properties and material characteristics of additively manufactured AlSi10Mg—Effect of the contour parameter on the microstructure, density, residual stress, roughness and mechanical properties. Int. J. Fatigue 2018, 117, 148–162. [Google Scholar] [CrossRef]
  26. Siddique, S.; Muhammad, I.; Wycisk, E.; Emmelmann, C.; Walther, F. Influence of process-induced microstructure and imperfections on mechanical properties of AlSi12 processed by selective laser melting. J. Mater. Proc. Tech. 2015, 221, 205–213. [Google Scholar] [CrossRef]
  27. Siddique, S.; Awd, M.; Tenkamp, J.; Walther, F. Development of a stochastic approach for fatigue life prediction of AlSi12 alloy processed by selective laser melting. Eng. Failure Anal. 2017, 79, 34–50. [Google Scholar] [CrossRef]
  28. Di Giovanni, M.T.; de Menezes, J.T.; Bolelli, G.; Cerri, E.; Castrodeza, E.M. Fatigue crack growth behavior of a selective laser melted AlSi10Mg. Eng. Fract. Mech. 2019, 217, 106564. [Google Scholar] [CrossRef]
  29. Awd, M.; Siddique, S.; Walther, F. Microstructural damage and fracture mechanisms of selective laser melted Al-Si alloys under fatigue loading. Theor. Appl. Fract. Mech. 2020, 106, 102483. [Google Scholar] [CrossRef]
  30. Sienkiewicz, J.; Płatek, P.; Jiang, F.; Sun, X.; Rusinek, A. Investigations on the mechanical response of gradient lattice structures manufactured via SLM. Metals 2020, 10, 213. [Google Scholar] [CrossRef] [Green Version]
  31. Xiao, L.; Song, W.; Wang, C.; Liu, H.; Tang, H.; Wang, J. Mechanical behaviour of open-cell rhombic dodecahedron Ti-6Al-4V lattice structure. Mater. Sci. Eng. A 2015, 640, 375–384. [Google Scholar] [CrossRef]
  32. Maconachie, T.; Leary, M.; Lozanovski, B.; Zhang, X.; Qian, M.; Faruque, O.; Brandt, M. SLM lattice structures: Properties, performance, applications and challenges. Mater. Des. 2019, 183, 108137. [Google Scholar] [CrossRef]
  33. Boniotti, L.; Beretta, S.; Patriarca, L.; Rigoni, L.; Foletti, S. Experimental and numerical investigation on compressive fatigue strength of lattice structures of AlSi7Mg manufactured by SLM. Int. J. Fatigue 2019, 128, 105181. [Google Scholar] [CrossRef]
  34. Dallago, M.; Winiarski, B.; Zanini, F.; Carmignato, S.; Benedetti, M. On the effect of geometrical imperfections and defects on the fatigue strength of cellular lattice structures additively manufactured via Selective Laser melting. Int. J. Fatigue 2019, 124, 348–360. [Google Scholar] [CrossRef]
  35. Zargarian, A.; Esfahanian, M.; Dadkhodapouor, J.; Ziaei-Rad, S.; Zamani, D. On the fatigue behavior of additive manufactured lattice structures. Theor. Appl. Fracture Mech. 2019, 100, 225–232. [Google Scholar] [CrossRef]
  36. Zhang, C.; Zhu, H.; Liao, H.; Cheng, Y.; Hu, Z.; Zeng, X. Effect of heat treatments on fatigue property of selective laser melting AlSi10Mg. Int. J. Fatigue 2018, 116, 513–522. [Google Scholar] [CrossRef]
  37. Tang, M.; Pistorius, P.C. Oxides, porosity and fatigue performance of AlSi10Mg parts produced by selective laser melting. Int. J. Fatigue 2017, 94, 192–201. [Google Scholar] [CrossRef]
  38. Muralidharan, U.; Manson, S.S. Modified universal slopes equation for estimation of fatigue characteristics. Trans. ASME J. Eng. Mater. Tech. 1988, 110, 55–58. [Google Scholar] [CrossRef]
  39. Glodež, S.; Klemenc, J.; Zupanič, F.; Vesenjak, M. The high-cycle fatigue behaviour and fracture analysis of SLM AlSi10Mg alloy. Trans. Nonferrous Metals Soc. China 2020, in press. [Google Scholar]
  40. Ashby, M.F.; Evans, A.; Fleck, N.A.; Gibson, L.J.; Hutchinson, J.W.; Wadley, H.N.G. Metal Foams: A Design Guide; Elsevier Science: Burlington, MA, USA, 2000. [Google Scholar]
  41. Lehmhus, D.; Vesenjak, M.; de Schampheleire, S.; Fiedler, T. From stochastic foam to designed structure: Balancing cost and performance of cellular metals. Materials 2017, 10, 922. [Google Scholar] [CrossRef] [Green Version]
  42. Lim, T.-C. Auxetic Materials and Structures; Springer: Berlin/Heidelberg, Germany, 2015; ISBN 978-981-287-274-6. [Google Scholar]
  43. Lakes, R.S. Foam structures with a negative poisson’ s ratio. Science 1987, 235, 1038–1040. [Google Scholar] [CrossRef]
  44. Salit, V.; Weller, T. On the feasibility of introducing auxetic behavior into thin-walled structures. Acta Mater. 2009, 57, 125–135. [Google Scholar] [CrossRef]
  45. Novak, N.; Vesenjak, M.; Ren, Z. Auxetic cellular materials—A Review. Strojniški Vestn. J. Mech. Eng. 2016, 62, 485–493. [Google Scholar] [CrossRef] [Green Version]
  46. Yao, Y.T.; Uzun, M.; Patel, I. Workings of auxetic nano-materials. J. Achiev. Mater. Manuf. Eng. 2011, 49, 585–593. [Google Scholar]
  47. Schwerdtfeger, J.; Heinl, P.; Singer, R.F.; Körner, C. Auxetic cellular structures through selective electron-beam melting. Phys. Status Solidi 2010, 247, 269–272. [Google Scholar] [CrossRef]
  48. Grima-Cornish, J.N.; Grima, J.N.; Attard, D. A novel mechanical metamaterial exhibiting auxetic behavior and negative compressibility. Materials 2019, 13, 79. [Google Scholar] [CrossRef] [Green Version]
  49. Schwerdtfeger, J.; Schury, F.; Stingl, M.; Wein, F.; Singer, R.F.; Körner, C. Mechanical characterisation of a periodic auxetic structure produced by SEBM. Phys. Status Solidi 2012, 249, 1347–1352. [Google Scholar] [CrossRef]
  50. Liu, L.; Hu, H. A review on auxetic structures and polymeric materials. Sci. Res. Essays 2010, 5, 1052–1063. [Google Scholar]
  51. Duncan, O.; Allen, T.; Foster, L.; Senior, T.; Alderson, A. Fabrication, characterisation and modelling of uniform and gradient auxetic foam sheets. Acta Mater. 2017, 126, 426–437. [Google Scholar] [CrossRef]
  52. Dobnik Dubrovski, P.; Novak, N.; Borovinšek, M.; Vesenjak, M.; Ren, Z. In-plane behavior of auxetic nonwoven fabrics based on rotating square unit geometry under tensile load. Polymers 2019, 11, 1040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ng, W.S.; Hu, H. Woven fabrics made of auxetic plied yarns. Polymers 2018, 10, 226. [Google Scholar] [CrossRef] [Green Version]
  54. Zulifqar, A.; Hu, H. Geometrical analysis of bi-stretch auxetic woven fabric based on re-entrant hexagonal geometry. Text. Res. J. 2019. [Google Scholar] [CrossRef]
  55. McDonald, S.A.; Dedreuil-Monet, G.; Yao, Y.T.; Alderson, A.; Withers, P.J. In situ 3D X-ray microtomography study comparing auxetic and non-auxetic polymeric foams under tension. Phys. Status Solidi 2011, 248, 45–51. [Google Scholar] [CrossRef]
  56. Novak, N.; Vesenjak, M.; Tanaka, S.; Hokamoto, K.; Ren, Z. Compressive behaviour of chiral auxetic cellular structures at different strain rates. Int. J. Impact Eng. 2020, 141, 103566. [Google Scholar] [CrossRef]
  57. Imbalzano, G.; Tran, P.; Ngo, T.D.; Lee, P.V. Three-dimensional modelling of auxetic sandwich panels for localised impact resistance. J. Sandw. Struct. Mater. 2015. [Google Scholar] [CrossRef]
  58. Nečemer, B.; Kramberger, J.; Vuherer, T.; Glodež, S. Fatigue crack initiation and propagation in re-entrant auxetic cellular structures. Int. J. Fatigue 2019, 126, 241–247. [Google Scholar] [CrossRef]
  59. Tomažinčič, D.; Vesenjak, M.; Klemenc, J. Prediction of static and low-cycle durability of porous cellular structures with positive and negative Poisson’s ratios. Theor. Appl. Fract. Mech. 2020, 106, 102479. [Google Scholar] [CrossRef] [Green Version]
  60. Novak, N.; Vesenjak, M.; Ren, Z. Computational simulation and optimisation of functionally graded auxetic structures made from inverted tetrapods. Phys. Status Solidi B 2017, 254. [Google Scholar] [CrossRef]
  61. Al-Rifaie, H.; Sumelka, W. The Development of a new shock absorbing uniaxial graded auxetic damper (UGAD). Materials 2019, 12, 2573. [Google Scholar] [CrossRef] [Green Version]
  62. Novak, N.; Hokamoto, K.; Vesenjak, M.; Ren, Z. Mechanical behaviour of auxetic cellular structures built from inverted tetrapods at high strain rates. Int. J. Impact Eng. 2018, 122, 83–90. [Google Scholar] [CrossRef]
  63. Foster, L.; Peketi, P.; Allen, T.; Senior, T.; Duncan, O.; Alderson, A. Application of auxetic foam in sports helmets. Appl. Sci. 2018, 8, 354. [Google Scholar] [CrossRef] [Green Version]
  64. Airoldi, A.; Bettini, P.; Panichelli, P.; Oktem, M.F.; Sala, G. Chiral topologies for composite morphing structures—Part I: Development of a chiral rib for deformable airfoils. Phys. Status Solidi B 2015, 1445, 1435–1445. [Google Scholar] [CrossRef]
  65. Airoldi, A.; Bettini, P.; Panichelli, P.; Oktem, M.F.; Sala, G. Chiral topologies for composite morphing structures—Part II: Novel configurations and technological processes. Phys. Status Solidi B 2015, 1454, 1446–1454. [Google Scholar] [CrossRef]
  66. Ngo, T.D.; Kashani, A.; Imbalzano, G.; Nguyen, K.T.Q.; Hui, D. Additive manufacturing (3D printing): A review of materials, methods, applications and challenges. Compos. Part. B Eng. 2018. [Google Scholar] [CrossRef]
  67. Borovinšek, M.; Novak, N.; Vesenjak, M.; Ren, Z.; Ulbin, M. Designing 2D auxetic structures using multi-objective topology optimisation. Mater. Sci. Eng. A 2020, in press. [Google Scholar]
  68. Ngnekoua, J.N.D.; Nadot, Y.; Henaff, G.; Nicolai, J.; Kan, W.H.; Cairney, J.M.; Ridosz, L. Fatigue properties of AlSi10Mg produced by additive layer manufacturing. Int. J. Fatigue 2019, 119, 160–172. [Google Scholar] [CrossRef]
  69. Ansys® Academic Research Mechanical Release; ANSYS, Inc. Southpointe, 2600 Ansys Drive: Canonsburg, PA, USA, 2020; R1 (20.1.0.2019111719). Available online: https://www.ansys.com/ (accessed on 26 May 2020).
  70. Stephens, R.I.; Fatemi, A.; Stephens, R.R.; Fuchs, H.O. Metal. Fatigue in Engineering; John Wiley & Sons Inc.: New York, NY, USA, 2001. [Google Scholar]
  71. MATLAB ver. R2020a (9.8.0.1323502) win64; The MathWorks Inc., 1 Apple Hill Drive: Natick, MA 01760-2098, USA, 2020. Available online: https://www.mathworks.com/ (accessed on 26 May 2020).
Figure 1. The behaviour of auxetic cellular structures under compressive and tensile loading.
Figure 1. The behaviour of auxetic cellular structures under compressive and tensile loading.
Metals 10 00945 g001
Figure 2. Base topology (a); medial axis topology (b); and rounded medial axis topology (c).
Figure 2. Base topology (a); medial axis topology (b); and rounded medial axis topology (c).
Metals 10 00945 g002
Figure 3. Selected quarter unit cell topologies and resulting auxetic structures.
Figure 3. Selected quarter unit cell topologies and resulting auxetic structures.
Metals 10 00945 g003
Figure 4. Geometry (a) and building direction (b) of a flat specimen. (all dimensions are in millimetres [mm]).
Figure 4. Geometry (a) and building direction (b) of a flat specimen. (all dimensions are in millimetres [mm]).
Metals 10 00945 g004
Figure 5. A light micrograph of the polished specimen [39].
Figure 5. A light micrograph of the polished specimen [39].
Metals 10 00945 g005
Figure 6. Stress–strain response of the SLM AlSi10Mg alloy.
Figure 6. Stress–strain response of the SLM AlSi10Mg alloy.
Metals 10 00945 g006
Figure 7. Boundary conditions.
Figure 7. Boundary conditions.
Metals 10 00945 g007
Figure 8. Deformed and undeformed structures.
Figure 8. Deformed and undeformed structures.
Metals 10 00945 g008
Figure 9. Location of the maximum equivalent plastic strain and its distribution.
Figure 9. Location of the maximum equivalent plastic strain and its distribution.
Metals 10 00945 g009
Figure 10. Equivalent plastic strain for different fillet radiuses.
Figure 10. Equivalent plastic strain for different fillet radiuses.
Metals 10 00945 g010
Figure 11. The fatigue life Ni of different topologies and fillet radiuses.
Figure 11. The fatigue life Ni of different topologies and fillet radiuses.
Metals 10 00945 g011
Table 1. Material data sheet of Selective Laser Melting (SLM) AlSi10Mg alloy.
Table 1. Material data sheet of Selective Laser Melting (SLM) AlSi10Mg alloy.
Layer thickness: 30μm
Volume rate with the standard parameters: 4.8 mm2/s
Surface roughness (as-built, cleaned): Ra 15–19 μm; Rz 96–115 μm
The density with the standard parameters: 2.68 g/cm³
Material composition (maximum values in %)
AlSiFeCuMnMgNiZnPbSnTi
Bal.9–110.550.050.450.2–0.450.050.10.050.050.15
Table 2. Material parameters of the SLM AlSi10Mg alloy.
Table 2. Material parameters of the SLM AlSi10Mg alloy.
Elastic Modulus
E (MPa)
Poisson’s Ratio
ν (-)
Yield Stress
Rp0,2 (MPa)
Ultimate Tensile Strength UTS
Rm (MPa)
Strain at Rupture
A5 (%)
True Fracture Strain
ε f (-)
709000.31803183.350.031
Table 3. Properties and results of the analysed topologies.
Table 3. Properties and results of the analysed topologies.
TopologyMass
(g)
Displacement
x-direction (mm)
Equivalent Elastic Strain (mm/mm)Equivalent Plastic Strain (mm/mm)Poisson Ratio
(-)
Unit Cell’s Poisson Ratio
(-)
10.100.9640.004080.0154−2.41−4.01
20.090.3810.003070.0180−0.95−2.71
30.120.2800.003620.0094−0.67−0.92
40.100.0800.003120.0022−0.20−0.35
50.090.0400.003230.0011−0.01−0.19

Share and Cite

MDPI and ACS Style

Ulbin, M.; Borovinšek, M.; Vesenjak, M.; Glodež, S. Computational Fatigue Analysis of Auxetic Cellular Structures Made of SLM AlSi10Mg Alloy. Metals 2020, 10, 945. https://doi.org/10.3390/met10070945

AMA Style

Ulbin M, Borovinšek M, Vesenjak M, Glodež S. Computational Fatigue Analysis of Auxetic Cellular Structures Made of SLM AlSi10Mg Alloy. Metals. 2020; 10(7):945. https://doi.org/10.3390/met10070945

Chicago/Turabian Style

Ulbin, Miran, Matej Borovinšek, Matej Vesenjak, and Srečko Glodež. 2020. "Computational Fatigue Analysis of Auxetic Cellular Structures Made of SLM AlSi10Mg Alloy" Metals 10, no. 7: 945. https://doi.org/10.3390/met10070945

APA Style

Ulbin, M., Borovinšek, M., Vesenjak, M., & Glodež, S. (2020). Computational Fatigue Analysis of Auxetic Cellular Structures Made of SLM AlSi10Mg Alloy. Metals, 10(7), 945. https://doi.org/10.3390/met10070945

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