Study of Microstructure and Properties of 316L with Selective Laser Melting Based on Multivariate Interaction Influence

+e selective laser melting technique is widely used in aerospace and biomedical industries, and the performance of formed 316L parts is significantly subject to the forming angle. As the selective laser melting 316L parts are constrained bymultiple performance indexes, the study involves multivariate interaction influenced on the forming parameters such as the angle with the xz plane, the angle with the xy plane, laser power, scan speed, powder thickness, and hatching space on the indexes like tensile strength, density, and surface roughness with linear regression equations based on multiobjective optimization to obtain the best process parameters. +e study results of microstructure performance of the formed 316L parts show that the angle with the xz plane has significant effect on the experiment indexes, while the layer thickness has the greatest effect on the indexes. After stretching, the molten pools are obviously elongated and the microstructure of the formed 316L parts is composed of equiaxed crystals and columnar crystals with a grain width of 0.28–0.4 nm.+e secondary growth of the dendrites is not obvious, and the crystallinity of the selective laser melting 316L parts is not as good as the standard parts, with the microstructure showing directional solidification due to grain refinement and microscopic distortion of crystals. As the fracture has dimples, it is a ductile fracture and typical plastic fracture.+e hardness near the fracture is higher than that of the substrate, whilst the indexes regarding the selective laser melting parts are higher than the ASTM-A182 and ASTM-F3184-16 standards. Since the theoretical model built in this study has less error, the findings have practical engineering application value.


Introduction
e selective laser melting technique is a newly developed powder bed fusion additive manufacturing technique. e principle is that the metal powder is fused by laser and solidified to form a molten pool, and then multistacked by layer to produce a three-dimensional part [1]. As for a small number of customized parts, the technique has several advantages, such as short time consumption, fast response, and high efficiency [2]. In recent years, additive manufacturing has made great breakthroughs in various industries, especially in aerospace and biomedical fields [3,4], showing broad application prospects. e selective laser melting process parameters have a great effect on the mechanical properties, surface quality, and friction performance of the formed parts. Some studies have shown that the laser power and scan speed have the most significant effect on the density of the sample [2,[5][6][7][8], and the scan speed, layer thickness, and manufacturing direction have a certain effect on the surface roughness and dimensional accuracy of the formed parts [9]. e lower the laser power, the lower the density and hardness of the sample, and the less the fusion at the boundary of the molten pool [10]. When the forming direction is inclined from 45°to 90°, the strength and fatigue resistance decrease and the elongation at break increases by more than 50% [11]. Wang et al. [12] studied the effect of the selective laser melting scanning strategy on residual stress and surface quality of the formed parts. e results show that the 316L stainless steel parts can achieve 99.37% of density with the scanning strategy of helical segmentation to effectively control the residual stress of the SLM 316L stainless steel parts. Dilip et al. [13] studied the effect of heat treatment on the microstructure and properties of the SLM HY100 steel. e results show that the microstructure of the parts after heat treatment is finer and more uniform with the mechanical properties being improved. Amato et al. [14] of the University of Tennessee in the United States studied the SLM Inconel 718 cylindrical parts.
e results show that the Vickers hardness of untreated, HIP-treated, and heat-treated parts is 3.9 GPa, 5.7 GPa, and 4.6 GPa, respectively. Wei et al. [15] studied the microstructure and properties of stainless steel nanohydroxyapatite medical composites with selective laser melting technique. e results show that the tensile strength, elastic modulus, and hardness of the SS/NHA prototype obtained under the optimum process conditions are higher than the human bone, which provides the possibility to prepare the SS/NHA porous scaffold with AP mechanical properties of the weight-bearing bone.
Selective laser melting direction is one of the important factors affecting the performance of the formed parts [16,17]. At present, most of the forming parameters are obtained without considering the forming direction (the angle with the xy plane is default as 90°), but practically, the changes in the angle with the xy plane and the angle with the xz plane bring significant variation of density, surface roughness (step effect), and mechanical properties of the formed parts [18][19][20]. In the laser forming process, the thermal gradient and the cooling rate of the formed part change with the forming angle and a complex dynamic fluid flow is generated in the molten pool. ese complex phenomena affect crystal growth and orientation [21], resulting in material sputtering and microstructure defects, mainly manifested by micropores, high residual stress, and large surface roughness [12,22,23]. ese defects cause premature failure under fatigue loading. At present, there are many studies on the forming parameters of the 316L stainless steel with selective laser melting, mainly focusing on the effect of the single factor on the performance of the single formed part [24,25], rather than the effect of multifactors on multiple indexes. e study on the interaction effect of multiple forming directions on the density, surface roughness, and tensile strength of the formed parts has not been reported [26].
In this study, tensile parts with various forming angles were prepared by selective laser melting for the test of density, surface roughness, and mechanical properties of the parts, analysis of microstructure and fracture morphology of the parts with various forming angles, and exploration of the microstructure properties and forming mechanism in different areas and the effect on tensile mechanical properties, providing a theoretical basis for the control of microstructure and mechanical properties of the selective laser melting parts.

Experimental Equipment and Materials.
In the experiment, the following products of Guangzhou Riton Laser Co., Ltd. were selected: the D280 AM machine with an accuracy of ±0.02 mm, laser wavelength of 1.064 μm, focused spot diameter of 20-80 μm, and layer thickness of 20-80 μm for forming parts; WD-E stretching tester with a pressure sensor of 0-20 kN for the tensile test; KERN analytical balance ABJ 320-4NM with a range of 320 g and accuracy of 0.0001 g for quality measurement; TIME3220T roughness meter with a measurement platform size of 400 mm × 250 mm × 70 mm, lift height of 300 ±1 mm, x-direction fine tune of ±12.5 mm, y-direction fine tune of ±12.5 mm, rotating coarse tune of 360°, and fine tune of ±5°for detecting surface roughness; MQX1200-40 atmosphere box furnace with a maximum temperature of 1200°and temperature control accuracy of ±1°for heat treatment; XQ-1 metallographic specimen mounting press and YMPZ-1 automatic metallographic specimen grinding and polishing machine for metallographic preparation; EVO18 ZEISS scanning electron microscope with a magnification of 5-1,000,000x and acceleration voltage of 0.2-30 kV for microstructure observation; and SmartLab (9) X-ray diffractometer (XRD) with an X-ray generator of fine focusing 9 kW knob, X-ray generator filament of 0.4 mm × 8 mm, and goniometer minimum step size of 1/10000°for phase analysis.
316L stainless steel powder (MTI S01C) is selected from Guangzhou Material Technology Innovation Co., Ltd. with the chemical composition parameters shown in Table 1, and the microstructure of powder is shown in Figure 1. e powder diameter is D10: 23.48 μm, D50: 31.13 μm, and D90: 42.26 μm.

Experimental Design.
e response surface method is to study the influence of multiple factors on indicators based on both mathematics and statistics. e least square method is used to obtain the regression equation by second-order polynomial fitting, and the most suitable model is selected according to the significance test. e true functional relationship between the response Y and the independent variables ξ 1 , ξ 2 , . . ., ξ k in the response surface method can be expressed as where f is unknown, ξ is another variable not included in f, and ε is the error term. For second-order response surfaces, the following polynomial is used: where b 0 is the intercept, b i , b ii, and b ij are the regression coefficients, and x i , x ii, and x ij are the horizontally encoding variables. is study involves six factors and three levels in the BOX for experimental design. e response surface method was used to optimize the process parameters of 316L stainless steel with selective laser melting. e laser power A, scan speed B, hatching space C, layer thickness D, the angle with the xy plane F (shown in Figure 2), and the angle with the xz plane E (shown in Figure 3) were selected. e horizontal values of each factor are shown in Table 2. e tensile specimen specifications are calculated according to ISO 6892-1: 2009, ASTM E8/8M. In this study, 316L stainless steel samples were formed with D280 (shown in Figure 4), the density of the sample Y 2 was measured by the Archimedes drainage method, the surface roughness of the sample Y 1 was measured with the roughness meter, and the tensile limit Y 3 and the elongation rate Y 4 were measured on the tensile test bench. e Archimedes principle formulas are as follows: where m 1 is the mass (g) of the formed part in the air before sealing wax; m 2 is the mass (g) of the formed part in the air after sealing wax; m 3 is the mass (g) of the formed part in the water after sealing wax; ρ is the density (%) of the sample; ρ 0 is the density of the wax, which is 0.880 g/cm 3 ; ρ 1 is the density of the distilled water, which is 1 g/cm 3 ; ρ 2 is the density (g/cm 3 ) of the formed part; and ρ 3 is the theoretical density of 316L stainless steel material, which is 7.98 g/cm 3 ; To ensure the adaptability and accuracy of the model, it is necessary to evaluate the predictability. e coefficient of determination R 2 and the adjusted coefficient of determination R 2 adj are often used to evaluate the approximation of the regression model [27]. e predicted coefficient of determination R 2 prep is adopted to evaluate the predictability of the regression model with the following formula: where S ST is the sum of squares; S SR is the sum of squared regression; S SE is the sum of squared residuals;ŷ is the predicted value of the response surface; y i is the true value of the response at the i th observation; � y is the mean of the response; and n is the number of alignment runs or observations designed for the experiment. R 2 is a fully fitted metric value, usually requiring a value above 0.9.

Advances in Materials Science and Engineering
Adjusted coefficient of determination R 2 adj is obtained by where p is the number of regression coefficients in the model and R 2 adj refers to the degree of correlation between all independent variables and dependent variables.
Adjusted coefficient of determination R 2 prep is obtained by where V is the sum of the squares of the predicted errors;ŷ i is the predicted value of the response at the i th time; and R 2 prep represents the predictability based on the original regression model fitted to the data. In general, for a model with acceptable accuracy, the values of R 2 prep and R 2 should be less than or equal to 0.3.

Analysis of Response Surface Method.
Experimental analysis was performed using Design-Expert. e experimental results are shown in Table 3.

Analysis of Variance.
e data in Table 3 were treated using multivariate regression fitting to obtain the regression equations of the response surface roughness Y 1 , density Y 2 , tensile strength Y 3, and elongation Y 4, respectively, expressed by the laser power A, scan speed B, hatching space C, layer thickness D, the angle with the xz plane E, and the angle with xy plane F: e response value and variance analysis are shown in Tables 4-7. e p values of the four models are all less than 0.01, indicating that the selected factors in the equations are extremely significant, and p values of lack of fit are all greater than 0.05, which is not significant. For the surface roughness of the sample, the coefficient of determination R 2 is 0.9158, the adjusted coefficient of determination R 2 adj is 0.8283, and the predicted coefficient of determination R 2 prep is 0.5898, representing a difference of 0.326 between the coefficient of determination R 2 and the predicted coefficient of determination R 2 prep . For the density of the sample, the coefficient of determination R 2 is 0.8961, the adjusted coefficient of determination R 2 adj is 0.7881, and the predicted coefficient of determination R 2 prep is 0.5716, representing a difference of 0.3245 between the coefficient of determination R 2 and the predicted coefficient of determination R 2 prep . For the tensile strength, the coefficient of determination R 2 is 0.937, the adjusted coefficient of determination R 2 adj is 0.8715, and the predicted coefficient of determination R 2 prep is 0.6842, representing a difference of 0.2528 between the coefficient of determination R 2 and the predicted coefficient of determination R 2 prep . For the elongation of the sample, the coefficient of determination R 2 is 0.9359, the adjusted coefficient of determination R 2 adj is 0.8693, and the predicted coefficient of determination R 2 prep is 0.679, representing a difference of 0.2569 between the coefficient of determination R 2 and the predicted coefficient of determination R 2 prep . All the above indicates that the response surfaces of the four models have good predictability for the experimental results. e relationships between the actual values of surface roughness, density, tensile strength, and elongation and the predicted values of the regression models are shown in         Within the horizontal range of the selected factors, the weight of effect on the surface roughness according to the F value of each factor is the angle with the xy plane > layer thickness > the angle with the xz plane > laser power>hatching space > scan speed; the weight of effect on the density is layer thickness > hatching space > the angle with the xy plane > scan speed > laser power > the angle with the xz plane; the weight of effect on tensile strength is layer thickness > the angle with the xy plane > scan speed > laser power > hatching space > the angle with the xz plane; and the weight of effect on elongation is layer thickness > laser power > scan speed > hatching space > the angle with the xy plane > the angle with the xz plane.

Analysis of Multivariate Interactive Influence.
When laser power, scan speed, hatching space, and the angle with the xz plane are constant, it can be seen from Figure 6 that when the angle with the xy plane is 0°, the surface roughness increases as the layer thickness increases, reaches the peak as the layer thickness further increases, and then decreases as the layer thickness continues to As shown in Figure 7(a), when both laser power and scan speed increase, the density first increases and then decreases. When the laser power is 160 W and the hatching space is 0.07 mm, the density reaches the peak; when the hatching space ranges from 0.07 to 0.08 mm, the density increases as the laser power increases; and when the hatching space ranges from 0.08 to 0.09 mm, the density decreases as the laser power increases. As shown in Figure 7(b), when the layer thickness is constant, the density decreases as the scan speed increases; and when the scan speed is constant, the density increases as the layer thickness increases. As shown in Figure 7(c), when the layer thickness is 0.02 mm, the density increases first and then decreases as the hatching space increases; and when the scan speed is 0.09 mm/s, the density decreases at the fastest rate as the layer thickness increases. As shown in Figure 7(d), when the angle with the xy plane is constant, the density decreases as the layer thickness increases; when the angle with the xy plane is 0°, the density decreases at the fastest rate; when the layer thickness is constant, the density increases as the angle with the xy plane increases, and when the layer thickness is 0.04 mm, the density increases at the fastest rate as the angle with the xy plane increases. As shown in Figure 7(e), when the layer thickness is 0.02 mm, the density increases as the angle with the xz plane increases, reaching the maximum value, but the trend of density changes as the layer thickness further increases, and when layer thickness is above 0.02 mm, the density decreases as the angle with the xz plane increases. As shown in Figure 7(f ), when the angle with the xy plane ranges from 0°to 50°, the density decreases as the angle with the xz plane increases, and when the angle with the xy plane ranges from 50°to 90°, the density increases as the angle with the xz plane increases.
As shown in Figure 8(a), when the laser power is constant, the tensile strength decreases as the layer thickness increases; when the laser power is 140 W, the tensile strength decreases at the fastest rate; when the laser power is 160 W, the tensile strength decreases at the slowest rate; and when the layer thickness is constant, the tensile strength increases as the laser power increases. As shown in Figure 8(b), when the layer thickness is constant, the tensile strength decreases as the scan speed increases; when the layer thickness is 0.04 mm, the tensile strength decreases at the fastest rate; when the scan speed is constant, the tensile strength decreases as the layer thickness increases; and when the scan speed is 900 mm/s, the tensile strength decreases at the fastest rate.
As shown in Figure 9(a), when the scan speed ranges from 700 to 760 mm/s, the elongation increases as the laser power increases; when the scan speed ranges from 760 to 900 mm/s, the elongation first decreases and then increases as the laser power increases; when the laser power ranges from 140 to 148 W, the elongation increases as the scan speed increases; when the laser power ranges from 148 to 160 W, the elongation decreases first and then increases as the scan speed increases; and at the minimum scan speed and the maximum laser power, the elongation reaches the maximum value. As shown in Figure 9(b), when the angle with the xz plane ranges from 0°to 50°, the elongation decreases first and then increases as the laser power increases; when the angle with the xz plane ranges from 50°to 90°, the elongation increases as the laser power increases; when the laser power ranges from 140 to 150 W, the elongation decreases as the angle with the xz plane increases; and when the laser power ranges from 150 to 160 W, the elongation increases as the angle with xz plane increases, reaching the maximum value together with the angle with the xz plane, laser power, and elongation. As shown in Figure 9(c), when the scan speed ranges from 700 to 800 mm/s, the elongation decreases first and then increases as the layer thickness increases; when the scan speed ranges from 800 to 900 mm/s, the elongation decreases as the layer thickness increases; when the layer thickness ranges from 0.02 to 0.03 mm, the elongation does not change significantly as the scan speed increases; and when the layer thickness ranges from 0.03 to 0.04 mm, the elongation decreases as the scan speed increases, and the maximum elongation is obtained from the minimum layer thickness and maximum scan speed. As shown in Figure 9(d), when the angle with the xy plane is constant, the elongation decreases rapidly first and then increases slowly as the layer thickness increases; when the layer E : a n g le w it h X -a x is (°) F : a n g le w it h X Y p la n e (°)  thickness ranges from 0.02 to 0.023 mm, the elongation increases as the angle with the xy plane increases; and when the layer thickness ranges from 0.024 to 0.04 mm, the elongation decreases first and then increases as the angle with the xy plane increases.

Optimization of Process
Parameters. Based on the above experiments, the predicted and actual optimization results are shown in Table 8 and formula (8).
e error between the predicted and actual values of the surface roughness under optimized conditions is 4.6%, the error between the predicted and actual values of the density is 1%, the error between the predicted and actual values of the tensile strength is 19.5%, and the error between the predicted and actual values of the elongation is 18.9%. All errors are within the reasonable range, indicating that the results are reliable. e tensile strength and elongation of the formed parts with the optimal parameters are higher than that specified in ASTM-A182 and ASTM-F3184-16.

Microstructure and Fracture Morphology Analysis
3.2.1. Microstructure Analysis. As shown in Figure 10(a), the microstructure of the selective laser melting sample before stretching is mainly composed of a plurality of oval molten pools, each of which has an aspect ratio of 0.5-0.65. is study suggests that there is certain unevenness in the powder size, a difference in laser energy absorption, and a Gaussian distribution of laser energy, which helps shape the oval molten pools due to different energy densities in various areas. Figure 10(b) shows the shape of the molten pools after stretching. Compared with Figure 10 planar crystals at the junction of the layers. e columnar crystals in the microstructure show apparent epitaxial growth characteristics because the laser melting process is at a high temperature gradient and solidification rate, while forming material and substrate temperature are lower, generally lower than the forming temperature by several orders of magnitude, resulting in heat transfer from the processing layer to the formed layer from top to bottom. e solidification direction inside the molten pool is from the bottom to the top of the molten pool, with solidification conditions in the growth range of equiaxed crystals and columnar crystals. Due to the rapid condensation of the laser, the bottom of the molten pool is closest to the formed layer with fastest solidification. e crystal grains at the bottom of the molten pool are solidified without growing, so equiaxed crystals appear at the bottom and boundary of the molten pool in the form of relatively fine grains, showing a planar crystal with white and bright layers of microstructure.

Hardness Analysis.
e microhardness of the experiment is shown in Figure 11. e distance closest to the fracture is set as 0, and three points are selected at each distance to obtain the average. As shown in Figure 11, the microhardness of the sample increases first and then decreases with the increase in the distance. e average microhardness at 2.5 mm from the fracture is 28.73 HRC, and the average microhardness at 20 mm from the fracture is 25.23 HRC, which is 87% of the hardness at the fracture. e microhardness at the fracture is significantly higher than the substrate because the fracture has a large number of dislocations, and dislocation multiplication and increased surface energy help improve the hardness.

XRD Analysis.
e XRD of the sample formed by optimizing the parameters is shown in Figure 12. Compared with the standard map, the diffraction peak intensity of the crystal planes increases at c-Fe (111), and the diffraction peak intensity of the crystal planes decreases at c-Fe (200) and c-Fe (220), showing that the c-Fe crystallinity is not as good as the standard map. When the diffraction angle ranges from 40°to 45°, a significant diffraction peak broadening occurs. e diffraction peak broadening is determined by two physical factors, grain refinement and crystal microscopic distortion.

Fracture Analysis.
e fracture morphology is shown in Figure 13. As shown in Figure 13(a), there are a large number of clear dimples whose shapes have obvious tearing, called elongated dimples, and the specimen fracture is a typical plastic fracture. Figure 13(b) shows fracture cracks passing through the grains, which is known as a transgranular fracture, and there are many dimples in the fracture, which is a ductile fracture. e specimen fracture is mainly caused by a large amount of dislocation multiplication under tensile stress with increased grain boundary energy, and a great quantity of impurities in the

Conclusion
(1) e angle with the xz plane has a significant effect on the experiment indexes. e comprehensive factor that has the greatest effect on the experiment indexes is the layer thickness. e biggest effect on the surface roughness is from the angle with the xy plane.
(2) e optimal forming parameters are laser power of 160 W, scan speed of 836 mm/s, hatching space of 0.07 mm, layer thickness of 0.02 mm, the angle with the xz plane of 83.3°, and the angle with the xy plane of 89.7°. e errors between the predicted values and the actual values of the surface roughness, the density, the tensile strength, and the elongation under optimized conditions are 4.6%, 1%, 19.5%, and 18.9%, respectively, and the errors are within a reasonable range, indicating that the model has practical engineering value. (3) After stretching, the molten pools of the parts are obviously elongated with the aspect ratio increasing from 0.5-0.65 to 1.5-2, up by 130%-300%. e microstructure of the formed parts is composed of equiaxed crystals and columnar crystals with a grain width of 0.28-0.4 nm. e secondary growth of dendrites is not obvious, with the microstructure showing directional solidification. XRD analysis shows that the crystallinity of the selective laser melting parts is not as good as that of the standard sample, with grain refinement and microscopic distortion of crystals. (4) ere are many dimples in the fracture which is ductile fracture and typical plastic fractures, and the hardness near the fracture is greatly higher, or 1.14 times that of the substrate. e indexes regarding the selective laser melting parts are better than the ASTM-A182 and ASTM-F3184-16 standards.
Data Availability e figure and table data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.