Next Article in Journal
On Additive Manufactured AlSi10Mg to Wrought AA6060-T6: Characterisation of Optimal- and High-Energy Magnetic Pulse Welding Conditions
Next Article in Special Issue
Use of X-ray Computed Tomography for Assessing Defects in Ti Grade 5 Parts Produced by Laser Melting Deposition
Previous Article in Journal
Effects of Ce on Microstructure and Mechanical Properties of LDX2101 Duplex Stainless Steel
Previous Article in Special Issue
Residual Stress and Microstructure of a Ti-6Al-4V Wire Arc Additive Manufacturing Hybrid Demonstrator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Separation of the Formation Mechanisms of Residual Stresses in LPBF 316L

1
Bundesanstalt für Materialforschung und -prüfung (BAM), Unter den Eichen 87, 12205 Berlin, Germany
2
Institute of Machine Tools and Factory Management, Chair of Processes and Technologies for Highly Loaded Welds, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
3
Heinz Maier-Leibnitz Zentrum (MLZ), Technische Universität München, Lichtenbergstraße 1, 85747 Garching, Germany
4
Institute of Physics and Astronomy, University of Potsdam, Karl-Liebknecht-Straße 24/25, 14476 Potsdam, Germany
*
Author to whom correspondence should be addressed.
Metals 2020, 10(9), 1234; https://doi.org/10.3390/met10091234
Submission received: 29 July 2020 / Revised: 21 August 2020 / Accepted: 3 September 2020 / Published: 14 September 2020

Abstract

:
Rapid cooling rates and steep temperature gradients are characteristic of additively manufactured parts and important factors for the residual stress formation. This study examined the influence of heat accumulation on the distribution of residual stress in two prisms produced by Laser Powder Bed Fusion (LPBF) of austenitic stainless steel 316L. The layers of the prisms were exposed using two different border fill scan strategies: one scanned from the centre to the perimeter and the other from the perimeter to the centre. The goal was to reveal the effect of different heat inputs on samples featuring the same solidification shrinkage. Residual stress was characterised in one plane perpendicular to the building direction at the mid height using Neutron and Lab X-ray diffraction. Thermography data obtained during the build process were analysed in order to correlate the cooling rates and apparent surface temperatures with the residual stress results. Optical microscopy and micro computed tomography were used to correlate defect populations with the residual stress distribution. The two scanning strategies led to residual stress distributions that were typical for additively manufactured components: compressive stresses in the bulk and tensile stresses at the surface. However, due to the different heat accumulation, the maximum residual stress levels differed. We concluded that solidification shrinkage plays a major role in determining the shape of the residual stress distribution, while the temperature gradient mechanism appears to determine the magnitude of peak residual stresses.

1. Introduction

In recent years, Additive Manufacturing (AM) has evolved from a method for rapid prototyping to a mature production process for certain parts in industries, such as the aerospace industry [1]. Among the different AM manufacturing techniques, Laser Powder Bed Fusion (LPBF) is an important technique for the production of net shaped metallic parts [2]. Early research conducted by Mercelis and Kruth [3] showed that metallic parts made by LPBF inherently contain residual stresses (RS). They had described two driving mechanisms for the formation of RS: the Temperature Gradient Mechanism (TGM) and the Solidification Shrinkage Mechanism (SSM). The two mechanisms are interlinked and their combined effect on RS in AM 316L is a topic of current research [4,5,6,7]. Wang et al. [6] showed within long bars of LPBF 316L that scan strategies using shorter scan tracks reduced RS and attributed this to lower solidification shrinkage. Roehling et al. [7] observed a decrease of RS in samples with a bridge geometry manufactured by LPBF of 316L due to post-solidification heating of each layer during the build job using selective large-area diode surface heating. This method aimed to decrease the cooling rate. Each of the two publications had mainly utilized one of the two mechanisms to reduce RS: Wang et al. [6] mainly utilized the SSM, whereas Roehling et al. [7] mainly utilized the TGM. In both cases, a reduction of RS was observed. However, there is still a level of uncertainty on the magnitude of the influence of each mechanism onto the shape of the resulting RS field.
Diffraction is a well-known non-destructive method to evaluate RS [8,9,10]. Determining elastic strains by measuring the variation of lattice spacing provides a powerful method to identify RS. This is achieved at the surface by Lab X-ray Diffraction (XRD), up to a depth of about 5 μm in metals, as well as in the bulk by Neutron Diffraction (ND) up to a depth of about a few mm to a few cm [11,12,13]. In this work, the bulk triaxial RS state was determined using ND and was combined with the Lab XRD biaxial RS state at the surface. This methodology allows for the mapping of the RS distribution across the complete cross-sectional plane. Such residual stress tends to be compressive in the bulk and tensile near the surface [14]. ND enables the non-destructive determination of the triaxial RS state over a complete two-dimensional (2D) plane or three-dimensional (3D) volume. Destructive methods, such as incremental or deep hole drilling, slitting, or contour method, would also yield stress depth profiles, but it would be extremely difficult to determine triaxial stress states over a complete cross-section.
To exploit the benefits of lightweight, load driven structural designs for LPBF metallic parts, it is necessary to understand RS in those parts, since their effect on fatigue life can be significant [15]. In order to understand RS, it is necessary to decouple the contributing mechanisms, especially if we aim at modelling the manufacturing process.
Therefore, this study aims at unravelling the contributions of the two mechanisms, to provide a basis for discussion on the length scale of RS introduced into parts by the TGM and the SSM. Therefore, the specimen design was chosen to provide similar solidification shrinkage, but at the same time different cooling rates without changing the volumetric energy density (VED) of specimens. Based on this design, similar RS results should be assigned to equal solidification shrinkage, whereas differences should be caused by the different cooling rates. Online monitoring by thermography during the build process was used in order to assess these cooling rates and their effect on RS formation.
The TGM is mainly related to process parameters (e.g., VED) and the SSM is mainly related to the length of shrinking scan tracks. Therefore, the results of this study might help to decide which approach is more suitable to reduce RS for a specific part design and its expected load profile.
Additionally, the results from Micro Computed Tomography (μCT) and Optical Microscopy (OM) were evaluated to link RS fields and defect distributions, with the aim to produce a holistic approach towards the analysis of the interconnection of TGM and SSM.

2. Materials and Methods

2.1. Material and LPBF Processing Conditions

Austenitic stainless steel 316L powder was processed by the commercial LPBF system SLM280 HL (SLM Solutions Group AG, Lübeck, Germany). The powder was characterised by its supplier: it has an apparent density of 4.58 g cm−3 and a mean diameter of 34.69 μm. The cumulative mass values of the particle size distribution are: D10 = 18.22 μm, D50 = 30.50 μm, D90 = 55.87 μm. The LPBF system uses a single 400 W continuous wave ytterbium fibre laser with a spot size of approx. 80 μm in a focal position. The processes were conducted in an argon gas atmosphere with an oxygen content of less than 0.1%. The parts were manufactured on a stainless steel substrate plate, which was heated up to 100 °C as a preheating temperature before the start of the build process. Two prismatic specimens of the dimension 24 mm × 36 mm × 24.5 mm were manufactured in two separate built processes. In order to remove specimens from the base plate a band saw was used. This reduced the height to a final value of 21 mm. A specimen design of low aspect ratio was chosen for this experiment in order to prevent significant RS relaxation due to distortion after the removal from the substrate plate. Such a distortion had been reported in literature for LPBF 316L [16,17]. Although the removal from the substrate plate may have caused a degree of RS relaxation, the overall relative trend between the specimens was considered to be mainly unaffected. The specimens were placed close to the border of the base plate to fit within the field of view of the thermography camera setup. The specimens were manufactured using the following process parameters: layer thickness t = 50 μm, scanning velocity v = 700 mm/s−1, laser power P = 275 W, and hatch distance h = 0.12 mm. Two different so-called border fill scanning strategies were applied, which scan along the edges of the rectangular cross sections of the parts: for one the scanning sequence starts in the centre of the part with growing rectangles towards the perimeter and the other has a converse scanning sequence (see Figure 1). The interlayer time (according to Mohr et al. [18]) was approximately 27.6 ± 1.0 s due to time variations between re-coating forwards and backwards. Therefore, the total time for each build process of 490 layers added up to 3.76 h.

2.2. Thermography

An ImageIR 8300 hp camera (Infratec GmbH, Dresden, Germany) working in the spectral range of 2– 5.7 μm was used for thermographic measurements. It was mounted on top of the SLM280 HL machine’s build chamber, observing the build plate through a sapphire window. The chosen subframe image had a size of 160   px × 114   px featuring a geometric resolution on the build plate of 360   μ m / px 1 . The acquisition frame rate was set to 1000 Hz . The camera was calibrated for black body radiation. Due to the fact that the emissivity of the used material is well below unity [19] and the process was observed through optical elements, the calibration is not valid for quantitative evaluation of the obtained thermography data. Nonetheless, assuming that the emissivity remains (approximately) constant during the build process, the obtained apparent temperatures enable comparisons within a single build process and between the two different build processes. The thermography data for the two specimens were obtained during the build process while using two different calibration ranges: 673 K to 1073 K for the CtP specimen ( Metals 10 01234 i001) and 623 K to 973 K for PtC specimen ( Metals 10 01234 i002). Several overlapping ranges were being tested during this experiment to find the optimum for these specimens. Values within the overlapping apparent temperature range of 673 K to 973 K can be compared between the two build processes. For quantitative evaluation of the process, it would also be necessary to address the additional error in temperature estimation introduced by the limited spatial resolution of 360 μ m px 1 . This error diminishes once the spatial temperature gradients have decreased due to lateral heat flow. These limitations are discussed in more detail by Mohr et al. [20]. Despite these limitations, qualitative analysis of the thermography data revealed results that contribute to the understanding of the presented RS results.

2.3. Lab X-ray Diffraction

A StressTech Xstress G3 X-ray diffraction instrument (Stresstech GmbH, Rennerod, Germany) was used in order to determine the RS distribution at the surface of the specimens according to the sin 2 ψ -method. Based on the assumption that that principal stresses are aligned with the geometrical axes of the specimens and the normal stress component can be neglected at the surface, the RS could be calculated from the slope of the linear fit of the lattice spacings over sin 2 ψ -plot [9,21]. The ψ -tilt was carried out in the angular range of ψ = −45° to ψ = +45° in 19 steps. The specimens were tilted around two perpendicular axes, to yield two perpendicular stress components. On the 36 mm surface, this corresponds to seven measurement positions of the prisms’ normal and longitudinal stress component (see blue circles in Figure 2b). On the 24 mm surface, five measurement positions correspond to the normal and transversal component of the RS distribution of the prisms. The exposure time for each acquisition was 5 s . The 311 diffraction line at a corresponding 2 ϑ angle of 152.26° was acquired using a Mn Kα radiation source and a 2 mm diameter collimator. Further details of the setup were described by Thiede et al. [13]. The software Xtronic (Stresstech GmbH, Rennerod, Germany) was used for data processing. The peak fitting process was performed using the Pearson VII function and the background was fitted with a parabolic function. The diffraction elastic constants (DEC) were calculated for austenitic steel 316L based on the Eshelby–Kroener model [22]. The calculated Young’s modulus of E 311 = 184 G Pa and Poisson’s ratio ν311 = 0.294 agree with values reported by Rangaswamy et al. [23], as well as with results from measurements and simulations of the DEC values of LPBF 316L reported by Chen et al. [24].

2.4. Neutron Diffraction

Stress determination by neutron diffraction (ND) was carried out at the STRESS-SPEC diffraction instrument [25] at the neutron facility FRM II in Munich, Germany (Figure 2a).
A bent Si400 single crystal monochromator was used to select the wavelength of 1.550 Å. The Fe311-peak was selected in order to characterise the RS distribution due to the low accumulation of inter-granular stresses reported for this reflection in conventional face-cubic-centred (fcc) iron materials [10].
To map the RS in the cross sectional plane at the mid build height of the specimen, a gauge volume of 2 m m × 2 m m × 2 m m was used in a grid of 7 × 5 measurement points. (see grey cubes in Figure 2b). The coordinate system is depicted in Figure 2b. σ L represents the RS along the y-direction, σ T along the x-direction and σ N along the z-direction, which was also the build direction.
A stress-free reference was needed to calculate strains from the measured d 311 lattice-spacing. A small cube with the size of 3 m m × 3 m m × 3 m m was sectioned from the bottom corner of a separate test build job of the PtC specimen ( Metals 10 01234 i002). This cube was regarded as free of Type I macro-stresses [9], due to the mechanical relaxation during sectioning. The strain can subsequently be derived from the measured ϑ angles using Bragg’s law [9].
ε = d 311 d 0 d 0 = sin ϑ 311 sin ϑ 0 1
Assuming that the principal geometric directions correspond with the principal stress directions Hooke’s law reads as the following, as described by Holden et al. [26]:
σ L , T , N = E 311 1 + ν 311 1 2 ν 311 1 ν 311 ε L , T , N + ν 311 ε T , N , L + ε N , L , T
The d 0 value was derived from the average of the ϑ 311 measurements of the cube in longitudinal (L), transversal (T) and normal (N) direction, where the normal direction corresponds to the build direction (as depicted in Figure 2b). The same DECs that were derived from the Eshelby–Kroener model were applied to both Lab XRD and ND results (see Section 2.3).

2.5. Micro Computed Tomography

The small reference cube for ND was studied using Micro Computed Tomography (μCT) to obtain a detailed dataset of the internal defect structure. The µCT measurements were performed at a GE v|tome| × 180/300 CT scanner (GE Sensing & Inspection Technologies GmbH, Wunstorf, Germany) using the 180 kV source at a voltage of 150 kV and a current of 40 μ A without any metal pre-filter. A voxel size of ( 3   μ m ) 3 was achieved. The analysis of the data was performed using the commercial software VG Studio MAX version 3.2.1 (Volume Graphics GmbH, Heidelberg, Germany). A lower threshold limit of 8 voxels was used for pore detection.

2.6. Optical Microscopy

For Optical Microscopy (OM) investigations of the microstructure, the bottom faces of the samples were ground, polished, and etched. Emery papers with 180, 320, 600 and 1200 grits followed by clothes with 3 μm and 1 μm were used. For etching the Bloech & Wedl II method [27] (a solution of 50 mL H2O, 50 mL HCl, 0.6 g K2S2O5) was applied. The microstructure was captured using a Olympus BX53M microscope with a DP74 camera module (Olympus Corporation, Tokyo, Japan). The analysis was performed using the software Olympus Stream Essentials (Olympus Corporation, Tokyo, Japan).

3. Results

3.1. In-Situ Thermography

The build jobs of both samples were supervised in-situ by thermography in order to receive more detailed information on the local variation of the temperature gradient and cooling rates. Figure 3a,b display the maximum apparent temperatures at the mid-height layer.
The thermography data were averaged over 40 layers at the mid build-height to reduce noise and the influence of smoulder and spatter. In relation to the height of the specimens of 21 mm , these 40 layers represent an average of the height range from 9.35 mm to 11.35 mm. The number of 40 layers (=2 mm) was chosen, since this corresponds to the size of the used gauge volume for ND (2 mm × 2 mm × 2 mm).
The scan strategy of the PtC specimen ( Metals 10 01234 i002) resulted in an increased heat accumulation in the center, compared to the CtP specimen ( Metals 10 01234 i001). This could be observed as an increase in the maximum temperature at the centre, when comparing Figure 3d,c.
Figure 3d shows the different apparent temperature values for the four sections of the plane. These differences were caused by a combination of the different surface roughness of each section and of the shadowing effects from the smoulder.
Figure 4 shows the cooling rate of both specimens. The cooling rate— d T d t was obtained by comparing two images at t = 1 ms and t = 2 ms after an apparent temperature of 700 K was reached for the last time at the surface. Typical times for cooling from the maximum temperature to 700 K were between 2 ms and up to 15 ms . The cooling time of 15 ms was observed at the centre of the PtC specimen ( Metals 10 01234 i002).
As displayed in Figure 4d, the PtC specimen ( Metals 10 01234 i002) also featured a lower cooling rate in the centre of the plane in addition to the higher maximum temperature depicted in Figure 3d. The CtP specimen ( Metals 10 01234 i001) shows a low cooling rate at the edges of the specimen indicating that the surrounding metal powder served as an heat insulator for conduction, as assumed in modelling [28].

3.2. Combined Neutron and Lab X-ray Diffraction

The ND results (bulk RS) were combined with Lab XRD results (surface RS) to show the complete stress distribution across the full cross section of the specimens’ middle plane. The contour plot function of the commercial software Origin 2018 (OriginLab Corporation, Northampton, USA), which is based on the Delaunay triangulation, was used to visualise the combined results. Figure 5 shows the combination of ND and XRD results. The RS distribution is visualised in longitudinal, transversal, and normal direction, where the normal direction corresponds to the build direction (Figure 2b). The XRD technique used only allows for the determination of stress components which are parallel to the surface (i.e., in-plane). Therefore, to visualise the third orthogonal stress component (perpendicular to the surface) in Figure 5 the following boundary condition was used: at the surfaces of the specimens corresponding to positions at y = 0 mm and y = 36 mm, the value of the longitudinal stress component was assumed to σL = 0 MPa. For the value of the transversal stress component at x = 0 mm and x = 24 mm it was assumed σT = 0 MPa, since these are free surfaces in these corresponding stress directions. In general, for the two scan strategies, a similar RS distribution was observed in the longitudinal, transversal and normal direction. This distribution is characterised by compressive RS within the bulk, balanced by tensile RS at the surface. Instead, the RS distribution in the longitudinal and transversal direction of each specimen is similar in terms of shape and magnitude, the RS distribution in the normal (i.e., building) direction differed from the other directions in terms of shape and magnitude.
The highest magnitude of RS of each specimen (maximum tensile or maximum compressive) were observed in the normal direction. However, the PtC specimen ( Metals 10 01234 i002) displayed higher bulk compressive stresses in all three orthogonal directions. The different stress values in the normal direction, as compared to the longitudinal and transversal direction, are also reflected in the lattice spacing of the the d 0 cube. Whereas, the d-spacing in the longitudinal and transversal direction are relatively similar to each other, in the normal direction a larger lattice spacing was measured (see Table 1). The d 0 value was obtained from averaging the three measured directions. This averaged value was used as a stress-free lattice parameter to calculate the RS.

3.3. Micro Computed Tomography

The reference cube was analysed by µCT to provide an example of local defect distribution in the specimen. Because it was cut from the corner where the laser path started and ended, it represented the area where the highest amount of defects was expected. The µCT results presented in Figure 6 reveal a network of defects at the location where the laser started and ended, as well as between the hatches. Because the same scan vector was used on each hatch and layer, the projection of defects onto one plane (Figure 6c) reveals the lack of fusion between neighbouring hatches. As reported in literature [29,30] alternating the orientation of scan vectors between layers prevents the formation of lack of fusion defects. Because the effect of shrinkage on RS was the subject of this study, scan vectors were not altered between layers to magnify such effect.
The largest defects were situated close to the edge of the sample (Figure 6b,c). A total porosity of 0.28% was observed.

3.4. Optical Microscopy

Optical microscopy images of the polished and etched specimen’s bottom surface are shown in Figure 7. Defects at the turn location of the laser are visible (Figure 7d,e). The bottom surface of both specimens are expected to be less effected by heat accumulation due to the smaller build height at the time the microstructure froze, and a better heat conduction into the build plate as compared to the top layers of the specimens. Nonetheless, defects (pores and voids) were observed in both specimens at the positions, where the laser turns by 90°.

4. Discussion

As mentioned in the introduction, Mercelis and Kruth [3] described two major driving mechanisms for the formation of RS in AM metallic parts made by LPBF: the temperature gradient mechanism (TGM) and the solidification shrinkage mechanism (SSM). Both of the mechanisms induced compressive stresses into the bulk material. The combined data of ND and Lab XRD showed the typical stress distribution pattern for metallic AM parts: compressive stresses in the bulk and tensile stresses close to the surface.
Because of the same geometric dimensions, it can be assumed that the solidification shrinkage effects in the two samples were equal. This implies that any differences in the RS distributions of the two samples should be attributed to the TGM.
According to the SSM, the induced RS distribution by solidification shrinkage depends on the length of the scanned laser tracks. The shrinkage of longer laser tracks introduces higher RS than shorter laser tracks [31]. Therefore, the longer laser tracks parallel to the y-axes of the specimens presented in Figure 5a–d) appear to be the main reason for the higher compressive stresses in the longitudinal stress distributions of the two specimens, as compared to the shorter laser tracks parallel to the x-axes, which seemed to result in lower compressive stresses for the transversal RS distributions of the two specimens.
The distribution of the normal component of the RS was assumed to be mostly independent of the SSM, since there were almost no restrictions to solidification shrinkage due to the free top surface during the layer-wise production. Therefore, the TGM was assumed to be the main mechanism to shape the RS distribution of the normal component. The RS distribution has the shape of a butterfly (see Figure 5e,f). The spikes of this butterfly pattern match the location where the laser turned by 90° (see Figure 8).
It is also observed that the normal component exhibits the highest tensile stresses (see Figure 5 and Table 2). The determined RS values are below the tensile yield strength ranges for LPBF 316L (450 MPa to 590 MPa ) as reported by Wang et al. [32].
Figure 8b shows that higher compressive stresses in the PtC specimen ( Metals 10 01234 i002) are localised close to the zone of the highest heat accumulation. This is in contrast to the CtP specimen ( Metals 10 01234 i001), which shows lower and more evenly distributed compressive stresses in the plane (see Figure 5). The higher maximum temperatures at the centre of the PtC specimen ( Metals 10 01234 i002) in combination with the slower cooling rate seem to result in larger RS in the centre of the plane as compared to the CtP specimen ( Metals 10 01234 i001). Lower contributions from the TGM to the compressive stress profile of the CtP specimen ( Metals 10 01234 i001) might be a result of the lack of heat accumulation in the centre of its plane (see Figure 3c and Figure 8a) and a faster cooling rate (see Figure 4c).
Line stress profiles (Figure 9) that were derived from the in-plane data presented in Figure 5 at a middle line of the plane at x = 12   m m emphasise the conclusion of SSM being the major mechanism to define the shape of the RS distribution. However, the TGM appears to influence the magnitude of the peak compressive stresses (see Figure 9).
It should be noted that, in Figure 9, only measured values are displayed. Therefore, the boundary condition values of 0 MPa were excluded from the longitudinal stress profile at y = 0 and y = 36 (see Figure 9a). The measured bulk values close to these surfaces support the assumption of zero stress at these surfaces.
Since the TGM does not require melting [33], it can be assumed that the compressive stress components induced by the TGM were formed in subjacent layers within the heat affected zone of the melt pool but below the layers that were remolten. It should be noted that a linear interpolation was applied to combine the surface measurement positions and the neutron data in Figure 9. Therefore, any possible sub-surface tensile peaks that had been reported by Mishurova et al. [11,34] in LPBF Ti-6Al-4V were not be taken into account in this study.
The μCT results (Figure 6) reveal a network of defects at the laser’s start and stop position. OM results (Figure 7d,e) show defects at the turn locations of the scan track. These regions have been measured to be under higher tensile RS (in particular, σT), which is assumed to arise from the longitudinal and transversal solidification shrinkage, as depicted in Figure 5. Therefore, any defects in these regions might have served as micro-crack initiators. The combination of pores and tensile stress at corners could also explain why the defects observed by μCT (Figure 6b) seem to be larger towards the outer edge: the tensile stress due to solidification shrinkage increases towards the edges of the specimens due to an increase length of the laser tracks. The analysis of the optical microscopy and μCT data indicates that an unknown amount of RS might had been relaxed by micro-cracks at the laser’s start and stop location.

5. Conclusions

Two prismatic AISI316L specimens using a border fill scan strategy were produced in order to differentiate the effect of the temperature gradient mechanism from the solidification shrinkage mechanism in AM metallic parts produced by LPBF. The following conclusions could be made:
  • A combination of surface and bulk residual stress results was needed to fully cover the (surface) tensile and (bulk) compressive regions of the in-plane stress distribution.
  • By comparing the two samples, a similar stress distribution was revealed for each of the three orthogonal directions (see Figure 5 and Figure 9), indicating that the solidification shrinkage mechanism is the main mechanism controlling the shape of the RS distribution in these samples. The temperature gradient mechanism seems to influence the magnitude of the compressive stresses without changing the overall pattern of the stress distribution.
  • In-situ thermography results of the sample exposed from the perimeter to the centre revealed a heat accumulation, which corresponds to highly localised compressive stresses due to the temperature gradient mechanism.

6. Outlook

In-situ Thermography recorded a detailed data set of apparent surface temperatures and enables an analysis of cooling rates. These data can be used in future simulations in order to model the RS field at the specimens’ mid height. In future studies, the RS distribution in the subsurface region might be resolved using additional techniques, such as hole drilling, slitting, X-ray with layer removal, deep hole drilling, or contour method.
Additionally, twin specimens of the two prisms will allow further experiments to study the effect on maximum values of the RS distribution, if the reference cube is cut from different positions within these twin specimens. In addition, these twins will allow a systematic μCT analysis of different positions within the two specimens, in particular to relate RS and defect distributions.

Author Contributions

Conceptualization, A.U., T.F., G.M. and S.J.A.; methodology, A.U., G.M., M.S., K.S., S.J.A., T.M., I.S.-M., T.F., M.H.; formal analysis, A.U., S.J.A., M.S., K.S.; investigation, A.U., M.S., K.S., S.J.A., T.M., I.S.-M., T.F., M.H.; writing—original draft preparation, A.U., G.M., M.S., K.S., S.J.A.; writing—review and editing, A.U., G.M., M.S., S.J.A., T.M., I.S.-M., M.H., A.E., G.B., T.F.; visualization, A.U., S.J.A.; supervision, A.E., G.B.; project administration, A.E., G.B.; funding acquisition, G.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work has been funded by the BAM Focus Area Materials project AGIL “Microstructure Development in Additively Manufactured Metallic Components: from Powder to Mechanical Failure” and ProMoAM “Process monitoring of Additive Manufacturing”. We are thankful for the financial support and the fruitful cooperation with all partners. This work is based upon experiments performed at the STRESS-SPEC instrument operated by FRM II at the Heinz Maier-Leibnitz Zentrum (MLZ), Garching, Germany.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; and in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
AMAdditive Manufacturing
FRM IIForschungs-Neutronenquelle Heinz Maier-Leibnitz (Research Reactor Munich II)
μCTMicro Computed Tomography
NDNeutron Diffraction
LPBFLaser Powder Bed Fusion
OMOptical Microscopy
RSResidual Stress
SSMSolidification Shrinkage Mechanism
TGMTemperature Gradient Mechanism
VEDVolumetric Energy Density
XRD(Lab) X-Ray Diffraction

References

  1. Ngo, T.D.; Kashani, A.; Imbalzano, G.; Nguyen, K.T.; Hui, D. Additive manufacturing (3D printing): A review of materials, methods, applications and challenges. Compos. Part B Eng. 2018, 143, 172–196. [Google Scholar] [CrossRef]
  2. Frazier, W.E. Metal Additive Manufacturing: A Review. J. Mater. Eng. Perform. 2014, 23, 1917–1928. [Google Scholar] [CrossRef]
  3. Mercelis, P.; Kruth, J.P. Residual stresses in selective laser sintering and selective laser melting. Rapid Prototyp. J. 2006, 12, 254–265. [Google Scholar] [CrossRef]
  4. Li, R.; Shi, Y.; Wang, Z.; Wang, L.; Liu, J.; Jiang, W. Densification behavior of gas and water atomized 316L stainless steel powder during selective laser melting. Appl. Surf. Sci. 2010, 256, 4350–4356. [Google Scholar] [CrossRef]
  5. Miranda, G.; Faria, S.; Bartolomeu, F.; Pinto, E.; Madeira, S.; Mateus, A.; Carreira, P.; Alves, N.; Silva, F.; Carvalho, O. Predictive models for physical and mechanical properties of 316L stainless steel produced by selective laser melting. Mater. Sci. Eng. A 2016, 657, 43–56. [Google Scholar] [CrossRef]
  6. Wang, D.; Wu, S.; Yang, Y.; Dou, W.; Deng, S.; Wang, Z.; Li, S. The Effect of a Scanning Strategy on the Residual Stress of 316L Steel Parts Fabricated by Selective Laser Melting (SLM). Materials 2018, 11, 1821. [Google Scholar] [CrossRef] [Green Version]
  7. Roehling, J.D.; Smith, W.L.; Roehling, T.T.; Vrancken, B.; Guss, G.M.; McKeown, J.T.; Hill, M.R.; Matthews, M.J. Reducing residual stress by selective large-area diode surface heating during laser powder bed fusion additive manufacturing. Addit. Manuf. 2019, 28, 228–235. [Google Scholar] [CrossRef]
  8. Bacon, G. Neutron Diffraction; Clarendon Press: Oxford, UK, 1975. [Google Scholar] [CrossRef]
  9. Withers, P.; Bhadeshia, H. Residual stress. Part 1—Measurement techniques. Mater. Sci. Technol. 2001, 17, 355–365. [Google Scholar] [CrossRef]
  10. Hutchings, M.; Withers, P.; Holden, T.; Lorentzen, T. Introduction to the Characterization of Residual Stress by Neutron Diffraction; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar] [CrossRef]
  11. Mishurova, T.; Cabeza, S.; Artzt, K.; Haubrich, J.; Klaus, M.; Genzel, C.; Requena, G.; Bruno, G. An Assessment of Subsurface Residual Stress Analysis in SLM Ti-6Al-4V. Materials 2017, 10, 348. [Google Scholar] [CrossRef] [Green Version]
  12. Nadammal, N.; Cabeza, S.; Mishurova, T.; Thiede, T.; Kromm, A.; Seyfert, C.; Farahbod, L.; Haberland, C.; Schneider, J.A.; Portella, P.D.; et al. Effect of hatch length on the development of microstructure, texture and residual stresses in selective laser melted superalloy Inconel 718. Mater. Des. 2017, 134, 139–150. [Google Scholar] [CrossRef]
  13. Thiede, T.; Cabeza, S.; Mishurova, T.; Nadammal, N.; Kromm, A.; Bode, J.; Haberland, C.; Bruno, G. Residual Stress in Selective Laser Melted Inconel 718: Influence of the Removal from Base Plate and Deposition Hatch Length. Mater. Perform. Charact. 2018, 7, 717–735. [Google Scholar] [CrossRef]
  14. Wu, A.S.; Brown, D.W.; Kumar, M.; Gallegos, G.F.; King, W.E. An Experimental Investigation into Additive Manufacturing-Induced Residual Stresses in 316L Stainless Steel. Metall. Mater. Trans. A 2014, 45, 6260–6270. [Google Scholar] [CrossRef]
  15. Riemer, A.; Leuders, S.; Thoene, M.; Richard, H.; Troester, T.; Niendorf, T. On the fatigue crack growth behavior in 316L stainless steel manufactured by selective laser melting. Eng. Fract. Mech. 2014, 120, 15–25. [Google Scholar] [CrossRef]
  16. Vrancken, B.; Wauthlé, R.; Kruth, J.P.; Van Humbeeck, J. Study of the influence of material properties on residual stress in selective laser melting. In Proceedings of the Solid Freeform Fabrication Symposium, Austin, TX, USA, 12–14 August 2013; pp. 393–407. [Google Scholar]
  17. Wang, D.; Yang, Y.; Yi, Z.; Su, X. Research on the fabricating quality optimization of the overhanging surface in SLM process. Int. J. Adv. Manuf. Technol. 2013, 65, 1471–1484. [Google Scholar]
  18. Mohr, G.; Altenburg, S.J.; Hilgenberg, K. Effects of inter layer time and build height on resulting properties of 316L stainless steel processed by laser powder bed fusion. Addit. Manuf. 2020, 32, 101080. [Google Scholar] [CrossRef]
  19. Hunnewell, T.S.; Walton, K.L.; Sharma, S.; Ghosh, T.K.; Tompson, R.V.; Viswanath, D.S.; Loyalka, S.K. Total Hemispherical Emissivity of SS 316L with Simulated Very High Temperature Reactor Surface Conditions. Nucl. Technol. 2017, 198, 293–305. [Google Scholar] [CrossRef]
  20. Mohr, G.; Altenburg, S.J.; Ulbricht, A.; Heinrich, P.; Baum, D.; Maierhofer, C.; Hilgenberg, K. In-Situ Defect Detection in Laser Powder Bed Fusion by Using Thermography and Optical Tomography—Comparison to Computed Tomography. Metals 2020, 10, 103. [Google Scholar] [CrossRef] [Green Version]
  21. Macherauch, E.; Mueller, P. Das sin2Psi-Verfahren der roentgenographischen Spannungsmessung. Z. Angew. Phys. 1961, 13, 305–312. [Google Scholar]
  22. Kroener, E. Berechnung der elastischen Konstanten des Vielkristalls aus den Konstanten des Einkristalls. Z. Phys. 1958, 151, 504–518. [Google Scholar] [CrossRef]
  23. Rangaswamy, P.; Griffith, M.; Prime, M.; Holden, T.; Rogge, R.; Edwards, J.; Sebring, R. Residual stresses in LENS (R) components using neutron diffraction and contour method. Mater. Sci. Eng. A 2005, 399, 72–83. [Google Scholar] [CrossRef]
  24. Chen, W.; Voisin, T.; Zhang, Y.; Florien, J.B.; Spadaccini, C.M.; McDowell, D.L.; Zhu, T.; Wang, Y.M. Microscale residual stresses in additively manufactured stainless steel. Nat. Commun. 2019, 10, 4338. [Google Scholar] [CrossRef] [Green Version]
  25. Hofmann, M.; Gan, W.; Rebelo-Kornmeier, J. STRESS-SPEC: Materials science diffractometer. J. Large-Scale Res. Facil. JLSRF 2015, 1, 6. [Google Scholar] [CrossRef]
  26. Holden, T.M.; Clarke, A.P.; Holt, R.A. Neutron diffraction measurements of intergranular strains in MONEL-400. Metall. Mater. Trans. A 1997, 28, 2565–2576. [Google Scholar] [CrossRef]
  27. Weck, E.; Leistner, E. Metallographische Anleitung zum Farbätzen nach dem Tauchverfahren.—3: Nichteisenmetalle, Hartmetalle und Eisenwerkstoffe, Nickel-Basis- und Kobalt-Basis-Legierungen = Metallographic Instructions for Colour Etching by Immersion.—3: Non-Ferrous Metals, Cemented Carbides and Ferrous Metals, Nickel-Base and Cobalt-Base Alloys; Deutscher Verlag für Schweisstechnik: Duesseldorf, Germany, 1986. [Google Scholar]
  28. Denlinger, E.R.; Gouge, M.; Irwin, J.; Michaleris, P. Thermomechanical model development and in situ experimental validation of the Laser Powder-Bed Fusion process. Addit. Manuf. 2017, 16, 73–80. [Google Scholar] [CrossRef]
  29. Thijs, L.; Verhaeghe, F.; Craeghs, T.; Van Humbeeck, J.; Kruth, J.P. A study of the microstructural evolution during selective laser melting of Ti–6Al–4V. Acta Mater. 2010, 58, 3303–3312. [Google Scholar] [CrossRef]
  30. Yadroitsev, I.; Thivillon, L.; Bertrand, P.; Smurov, I. Strategy of manufacturing components with designed internal structure by selective laser melting of metallic powder. Appl. Surf. Sci. 2007, 254, 980–983. [Google Scholar] [CrossRef]
  31. Bo, Q.; Yu-sheng, S.; Qing-song, W.; Hai-bo, W. The helix scan strategy applied to the selective laser melting. Int. J. Adv. Manuf. Technol. 2012, 63, 631–640. [Google Scholar] [CrossRef]
  32. Wang, Y.M.; Voisin, T.; McKeown, J.T.; Ye, J.; Calta, N.P.; Li, Z.; Zeng, Z.; Zhang, Y.; Chen, W.; Roehling, T.T.; et al. Additively manufactured hierarchical stainless steels with high strength and ductility. Nat. Mater. 2018, 17, 63–71. [Google Scholar] [CrossRef] [Green Version]
  33. Kruth, J.; Froyen, L.; Van Vaerenbergh, J.; Mercelis, P.; Rombouts, M.; Lauwers, B. Selective laser melting of iron-based powder. J. Mater. Process. Technol. 2004, 149, 616–622. [Google Scholar] [CrossRef]
  34. Mishurova, T.; Artzt, K.; Haubrich, J.; Requena, G.; Bruno, G. Exploring the Correlation between Subsurface Residual Stresses and Manufacturing Parameters in Laser Powder Bed Fused Ti-6Al-4V. Metals 2019, 9, 261. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematics of both border fill scan strategies. (a) Describes the Centre to Parameter (CtP) strategy indicated by the green arrow, while (b) shows the Parameter to Centre (PtC) scan strategy indicated by the blue arrow. The black arrows show the direction of scan of the laser around each border fill scan.
Figure 1. Schematics of both border fill scan strategies. (a) Describes the Centre to Parameter (CtP) strategy indicated by the green arrow, while (b) shows the Parameter to Centre (PtC) scan strategy indicated by the blue arrow. The black arrows show the direction of scan of the laser around each border fill scan.
Metals 10 01234 g001
Figure 2. Beamline setup and measurement positions.
Figure 2. Beamline setup and measurement positions.
Metals 10 01234 g002
Figure 3. Apparent (uncalibrated) maximum temperature of the two specimens obtained from thermography data acquired during the build process at the sample’s mid build height. Each border fill scan started and ended at the bottom right hand laser turn position of the images. (a,b) display the mid-height layer, (c,d) display an average of 40 layers at the mid-height.
Figure 3. Apparent (uncalibrated) maximum temperature of the two specimens obtained from thermography data acquired during the build process at the sample’s mid build height. Each border fill scan started and ended at the bottom right hand laser turn position of the images. (a,b) display the mid-height layer, (c,d) display an average of 40 layers at the mid-height.
Metals 10 01234 g003
Figure 4. Cooling rates measured by the different between images taken at t = 1 ms and t = 2 ms after an apparent temperature of 700 K was observed for the last time. Thermography data for both samples was obtained during the build process at the middle of the total build height. (a,b) display the mid-height layer, (c,d) display an average of 40 layers around the middle of the total build height.
Figure 4. Cooling rates measured by the different between images taken at t = 1 ms and t = 2 ms after an apparent temperature of 700 K was observed for the last time. Thermography data for both samples was obtained during the build process at the middle of the total build height. (a,b) display the mid-height layer, (c,d) display an average of 40 layers around the middle of the total build height.
Metals 10 01234 g004
Figure 5. Comparison of RS maps of the tw scan strategies including results from lab X-ray diffraction at the surface (The big semi-translucent squares in the bulk represent the almost cubic ND gauge volume (orientated differently for different stress components), whereas the small semi-translucent squares at the edges represent the lab X-ray measurement positions).
Figure 5. Comparison of RS maps of the tw scan strategies including results from lab X-ray diffraction at the surface (The big semi-translucent squares in the bulk represent the almost cubic ND gauge volume (orientated differently for different stress components), whereas the small semi-translucent squares at the edges represent the lab X-ray measurement positions).
Metals 10 01234 g005
Figure 6. μCT reconstructions of the ND reference cube sectioned from a twin PtC specimen ( Metals 10 01234 i002).
Figure 6. μCT reconstructions of the ND reference cube sectioned from a twin PtC specimen ( Metals 10 01234 i002).
Metals 10 01234 g006
Figure 7. Optical microscopy images of the samples’ bottom surface after polishing and etching.
Figure 7. Optical microscopy images of the samples’ bottom surface after polishing and etching.
Metals 10 01234 g007
Figure 8. Normal (i.e., build direction) stress component (contour lines) in MPa (Figure 5e,f) overlayed with thermography results (Figure 3c,d) to highlight the compressive stresses at the laser’s turn locations.
Figure 8. Normal (i.e., build direction) stress component (contour lines) in MPa (Figure 5e,f) overlayed with thermography results (Figure 3c,d) to highlight the compressive stresses at the laser’s turn locations.
Metals 10 01234 g008
Figure 9. Line scans in both samples for all three orthogonal directions using the combined data from ND and lab XRD.
Figure 9. Line scans in both samples for all three orthogonal directions using the combined data from ND and lab XRD.
Metals 10 01234 g009
Table 1. Distribution of orthogonal d-spacing values of the reference cube.
Table 1. Distribution of orthogonal d-spacing values of the reference cube.
Orientationd-SpacingError
Longitudinal1.07449 Å2.36 × 10−5 Å
Transversal1.07448 Å2.71 × 10−5 Å
Normal1.07493 Å2.41 × 10−5 Å
Table 2. Max. and min. values of the orthogonal stress components in the plane.
Table 2. Max. and min. values of the orthogonal stress components in the plane.
Stress ComponentCentre to Perimeter Metals 10 01234 i001Perimeter to Centre Metals 10 01234 i002
Max. [MPa]Min. [MPa]Max. [MPa]Min. [MPa]
σ L 100.8 ± 16.9 253.4 ± 15.8 86.4 ± 10.2 304.2 ± 15.6
σ T 140.6 ± 10.4 103.6 ± 15.4 122.6 ± 11.0 134.0 ± 15.8
σ N 295.6 ± 10.5 202.1 ± 14.2 290.5 ± 10.6 245.8 ± 19.5

Share and Cite

MDPI and ACS Style

Ulbricht, A.; Altenburg, S.J.; Sprengel, M.; Sommer, K.; Mohr, G.; Fritsch, T.; Mishurova, T.; Serrano-Munoz, I.; Evans, A.; Hofmann, M.; et al. Separation of the Formation Mechanisms of Residual Stresses in LPBF 316L. Metals 2020, 10, 1234. https://doi.org/10.3390/met10091234

AMA Style

Ulbricht A, Altenburg SJ, Sprengel M, Sommer K, Mohr G, Fritsch T, Mishurova T, Serrano-Munoz I, Evans A, Hofmann M, et al. Separation of the Formation Mechanisms of Residual Stresses in LPBF 316L. Metals. 2020; 10(9):1234. https://doi.org/10.3390/met10091234

Chicago/Turabian Style

Ulbricht, Alexander, Simon J. Altenburg, Maximilian Sprengel, Konstantin Sommer, Gunther Mohr, Tobias Fritsch, Tatiana Mishurova, Itziar Serrano-Munoz, Alexander Evans, Michael Hofmann, and et al. 2020. "Separation of the Formation Mechanisms of Residual Stresses in LPBF 316L" Metals 10, no. 9: 1234. https://doi.org/10.3390/met10091234

APA Style

Ulbricht, A., Altenburg, S. J., Sprengel, M., Sommer, K., Mohr, G., Fritsch, T., Mishurova, T., Serrano-Munoz, I., Evans, A., Hofmann, M., & Bruno, G. (2020). Separation of the Formation Mechanisms of Residual Stresses in LPBF 316L. Metals, 10(9), 1234. https://doi.org/10.3390/met10091234

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