Since the nanoindentation procedure is quite sensitive to a specimen’s surface imperfections [
45], voids [
46], tip contamination, and other alterations that may occur [
47], it is expected that several measurements may significantly differ from the rest of the measurements. Therefore, a suitable method that provides objective criteria for detection and elimination of potential outliers from each dataset should be considered. Hence, Grubbs’ test was applied to all data sets [
48]. Prior to the implementation of Grubbs’ test, normality of data distribution was tested using the Shapiro–Wilk test [
48]. Moreover, the influence of etching on nano-hardness results is also considered, given that the etchant increases the surface roughness of polished samples and may affect measurement accuracy as reported in [
49]. However, no significant influence of the etchant on the nano-hardness results of as-built specimen was found (
Table 3), probably due to the higher nano-indentation loads and depths that were used in the experiments. Experimental investigation of the etchant’s influence on the nano-hardness and the implementation of statistical tests in the data analysis ensure meaningful and reliable results. Nanoindentation measurements were performed away from visible defects at specimens’ surfaces to eliminate a potential influence of voids on results. Furthermore, a high number of repetitions was performed in each nanoindentation experiment and Grubbs’ test was applied to detect outliers affected by potential subsurface defects or other microstructural anomalies. It is worth noting that all the values reported in this research are expressed as mean values (STD).
3.1. Indentation Location Dependance
The nanoindentation procedure provides local information on the material, which can be generalized in the case of a homogeneous microstructure. However, L-PBF and E-PBF technologies do not produce homogeneous microstructures in many materials due to high cooling rates present during manufacturing [
50,
51]. Furthermore, it was found that element partitioning occurs in the L-PBF Ti6Al4V alloy during heat treatment [
52]. Therefore, local compositional variations are typical for the L-PBF Ti6Al4V. However, it is unclear whether the local compositional variations present in the microstructure result in different mechanical properties of annealed L-PBF Ti6Al4V alloy measured at different locations. For that reason, mechanical properties were tested at five different locations (i.e., five prior-β grains) with seven repetitions at each location. To enhance the robustness and reliability of the reported results, the high number of experiments was performed and Grubbs’ test for outlier detection was implemented. This approach allowed identification of potential outliers that could impact the representability of the results. If heterogeneity of the microstructure or local compositional variations influence the mechanical properties, it is expected that the results will be statistically different for at least one prior-β grain. However, the results of the nano-hardness and the Young’s modulus were the same, regardless of the selected prior-β grain (
Table 4), as results of nanoindentation experiments performed on the A specimen (
Figure 2a) showed. In fact, no statistically significant differences were found between the Young’s modulus and the nano-hardness values measured in different prior-β grains. This was confirmed by calculated
p-values of 0.95 for the Young’s modulus and 0.99 for the nano-hardness, using one-way analysis of variance (ANOVA) at the 0.05 level of significance. The high
p-values determined for the nano-hardness and the Young’s modulus indicate that there is no significant distinction between the groups being compared. It is worth noting that the Shapiro–Wilk normality test (S-W) was performed on data of each prior-β grain. In all cases the calculated
p-values were higher than α = 0.05 (
Table 4). Furthermore, assumption of equal variances applied in ANOVA was tested using Levene’s test, both on Young’s modulus and nano-hardness data. There was homogeneity of variances, since the calculated
p-values for Young’s modulus and nano-hardness data were 0.24 and 0.22, respectively. Thus, the application of the ANOVA method was justified, since evidence of non-normality was not found, and equal variances were present.
Therefore, nanoindentation can be performed on the L-PBF Ti6Al4V (ELI) alloy on different grains without adversely affecting the measurement results. One reason why the obtained nanoindentation results were insensitive to local microstructural variations in the L-PBF Ti6Al4V alloy is the size of the utilized Berkovich tip, which was substantially larger than microstructural features (
Figure 2b). Hence, the nanoindentation procedure for characterization of mechanical properties is suitable for the L-PBF Ti6Al4V alloy, especially when the material is in the as-built condition as stated in [
49]. The edges of the residual imprints of the utilized Berkovich tip, with indentation depths up to 3000 nm, were ~20 μm long (
Figure 2b). It was reported in [
51], that an energy density of 71.4 J/mm
3 resulted in average lath sizes of 0.68 μm and 1.8 μm for L-PBF Ti6Al4V alloy in as-built and heat-treated conditions, respectively. Furthermore, the average lath sizes of heat-treated L-PBF Ti6Al4V decrease as the energy density increases [
51]. In both cases, the average lath sizes were substantially smaller than the size of utilized Berkovich tip which further justifies its applicability when characterization of nano-mechanical properties of L-PBF Ti6Al4V alloy is needed. Besides that, a high number of repetitions are performed at different locations provides more robust and reliable results. It was already established in [
51], that high utilized values of energy densities (
ED > 37 J/mm
3) result in a strong texture with fine α/α′ laths inside the columnar microstructure of as-built L-PBF Ti6Al4V alloy. In each case, multiple laths were consistently located beneath the Berkovich tip during nanoindentation measurements. By increasing the number of repetitions, the variations caused by different orientations and potential local heterogeneities within the microstructure, specifically the fine α/α’ laths within the columnar microstructure, were accounted for using this approach.
3.2. Nano-Hardness and Young’s Modulus
Based on CSM tests, the indentation depth interval from 1000 to 2400 nm was selected for further data evaluation, since at indentation depths larger than 1000 nm, the Young’s modulus and the nano-hardness converge to a constant value, as can be seen in
Figure 3c–f. Moreover, higher indentation depths used within this research yielded more robust results and showed to a larger extent how exactly the indentation depth influences the nano-hardness and the Young’s modulus. In general, the CSM method is more reliable compared to the load–unload method, when Young’s modulus determination is needed as the loading regime is performed under small loading–unloading cycles [
53]. There, the Young’s modulus is measured at several points compared to only one point when the load–unload method is applied [
53].
The indentation size effect (ISE) is one of major concerns when considering results obtained using nano-indentation or low-force Vickers hardness methods [
54]. In general, the indentation size effect (ISE) is manifested as hardness increase with indentation depth decrease and becomes more important at depths of less than ~1000 nm [
55]. This is the most frequently seen effect (the normal ISE), as stated by Pharr et al. [
55]. However, the nano-hardness decrease with indentation-depth decrease (the reverse ISE) can occur as well [
55]. Given that nano-hardness and Young’s modulus results were determined using an evaluation interval ranging from 1000 to 2400 nm the influence of ISE was avoided. This ensured more relevant and robust nano-mechanical results. Surprisingly, the normal ISE was observed on the D
an specimen, while the reverse ISE was observed on D
ab, as can be seen in
Figure 3e,f by observing the LOWESS curves which represent mean nano-hardness values at given indentation depths. However, it is still unclear which mechanisms are responsible for the ISE [
55,
56].
The Young’s modulus of the as-built specimen (D
ab) reaches its constant value at lower indentation depths compared to the annealed specimen (D
an), as can be seen in
Figure 3c,d. This observation also holds for the nano-hardness values as can be seen in
Figure 3e,f. The Young’s modulus and nano-hardness of the annealed specimen reaches a constant value at ~1000 nm, while a constant value for the as-built specimen is reached at ~300 nm. Due to extremely high cooling rates during L-PBF (10
4–10
6 K/s) [
57], the microstructure of the as-built Ti6Al4V (ELI) alloy consists of acicular martensite (α’) inside columnar prior-β grains (
Figure 3a). Using annealing heat treatment, the nano-hardness of the material can be reduced by means of transformation of the acicular martensite (α’) into α + β laths (
Figure 3b). Since the mechanical properties of
α and
β laths differ from each other, the nano-hardness and Young’s modulus data of the D
an specimen are more scattered compared to the D
ab specimen which has a more uniform microstructure (α’) as shown in
Figure 3.
Our reported Young’s modulus values of the D
an and D
ab specimens (
Table 5) indicate that annealing heat treatment did not influence the mean Young’s modulus values, since the
p-value calculated using the t-test was 0.12. Thus, statistically significant difference in mean values of D
an and D
ab specimens do not exist. It is worth noting that the normality of distribution of the Young’s modulus values for each specimen was tested as well. In both cases, the calculated
p-values were higher than the level of significance α = 0.05 (the
p-value of the D
ab specimen was 0.92, while the
p-value for the D
an specimen was 0.55). These results indicate that the applicability of the t-test for testing differences between mean values is justified. Liu et al. in their work found that different applied heat treatments have little effect on the Young’s modulus of L-PBF Ti6Al4V alloy [
58], which is consistent with these findings.
When the nano-hardness values of the D
an and D
ab specimens were compared (
Table 5), the mean nano-hardness of the as-built specimen was slightly higher than the mean nano-hardness of the annealed specimen. The calculated
p-value of 0.036 using the t-test indicates that a statistically significant difference between two observed mean values exists. However, this value is quite close to the level of significance of α = 0.05. In both cases, calculated
p-values using S-W were higher than α = 0.05, which further justifies the applicability of t-tests (the
p-value of the D
ab specimen was 0.09, while the
p-value for the D
an specimen was 0.56). These results indicate that the annealing heat treatment slightly reduced the nano-hardness values measured using nanoindentation methods. When hardness results obtained using low-force Vickers hardness tests (HV1) are compared, the difference in mean values of annealed and as-built specimens is more pronounced in contrast to results obtained using the nanoindentation method.
Results of S-W applied on HV1 data did not show evidence of non-normality, since the calculated p-values for the Dab and Dan specimens were 0.43 and 0.12, respectively. Furthermore, the difference in mean HV1 values of the as-built and annealed specimens is more pronounced, since the calculated p-value using the t-test is <0.001. This confirms that the annealing procedure significantly reduced the HV1 values of the Dan specimen. Hence, annealing heat treatment showed higher influence on hardness values obtained using the low-force Vickers hardness method (HV1) and lower influence on nano-hardness (H) values. Furthermore, the annealing heat treatment does not influence the Young’s modulus values.
The higher dislocation density of dominant α’ phase is the main reason why the as-built specimen has a higher HV1 value compared to the annealed specimen that has dominant α + β laths in its microstructure [
58]. Chen et al. in their work measured nano-hardness in different planes of the as-built L-PBF Ti6Al4V alloy and reported nano-hardness values of 4.2 ± 0.5 GPa and 5.1 ± 0.5 GPa for different planes [
36]. Those results are comparable to results shown in
Table 5, which confirms their relevance.
As shown in
Table 5, the highest mean Young’s modulus value of annealed specimens was obtained for test specimen G, when a combination of a laser power of 250 W and a scanning speed of 1000 mm/s (the highest energy density; 111 J/mm
3) was used. Calculated
p-values using Dunn’s multiple comparison test (
Table 6) on Young’s modulus data of the G specimen indicate statistically significant differences when compared to other specimens manufactured using the lowest and intermediate laser power levels, except for specimen F. Hence, the laser power might have a possible influence on Young’s modulus data as well. It is worth noting that the measured Young’s modulus values of specimen F were not normally distributed, as shown in
Table 5, since the calculated
p-value using the S-W test was 0.037. To investigate the possible effect of laser power on nano-hardness values in more detail, a wider range with a higher number of laser power levels should be incorporated within DoE.
When the nano-hardness results are compared (
Table 6), it is evident that no statistically significant difference exists between specimens manufactured using the highest utilized laser power level (i.e., specimens G, H, and I). Additionally, statistically significant differences in nano-hardness values were not found between specimens produced using 200 and 225 W laser power levels. The mean nano-hardness value calculated on the G specimen was higher than the mean nano-hardness values of the other specimens manufactured using 200 and 225 W laser power levels (
Table 5). In this case too, the maximum nano-hardness value measured using the nanoindentation procedure is found on the G specimen (i.e., the specimen manufactured using the highest energy density). This implies that the performed annealing heat treatment cannot completely eliminate the influence of utilized L-PBF process parameters.
Cepeda-Jiménez et al. [
51] in their work have described the influence of the energy density on the microstructural and texture evolution of annealed L-PBF Ti6Al4V alloy, and reported results in which a slight hardness increase trend with utilized energy densities in a range from 24.2 to 44.3 J/mm
3 can be seen. However, the trend changed when an energy density of 71.4 J/mm
3 was used [
51]. Hence, the hardness is not proportional to the utilized energy density. Moreover, similar energy densities can be achieved using completely different L-PBF process parameters as shown in
Table 5. For instance, a statistically significant difference exists between the nano-hardness values measured on specimens A and H, which were manufactured utilizing an identical energy density achieved using different laser power and scanning speed combinations (
Table 5). It has been documented in [
59] that specimens manufactured using higher energy densities are subjected to lower thermal gradients, since they remain at high temperatures for a longer period of time. Therefore, as stated in [
51], a higher stability of microstructures is present which results in a smaller driving force for grain coarsening during heat treatment, resulting in highly textured microstructures. However, it is unexpected that the G specimen, which has been subjected to annealing heat treatment, has a higher mean nano-hardness value than the as-built D
ab specimen (
Table 5).
Using the low-force Vickers hardness test, it was confirmed that the G specimen (i.e., specimen manufactured using the highest energy density value) has indeed a high HV1 value, 374 HV1 (std. 6 HV1), when annealed specimens are considered. However, that value is not higher than the HV1 value of the as-built D
ab specimen, 385 HV1 (std. 6 HV1). This HV1 hardness discrepancy may be attributed to differences in indentation depths and strain rates between the two different indentation methods. In other words, the indenter tip in the low-force Vickers hardness test is applied to larger indentation depths using different strain rates and has an even larger size than the Berkovich tip. Despite these systematic differences, both nano-hardness and low-force Vickers hardness tests provided results that are comparable with results presented in the literature [
58,
60,
61].
Furthermore, the Young’s modulus values for specimens manufactured using a laser power of 250 W have slightly higher mean values and lower standard deviations than other reported results in
Table 5. The Young’s modulus reached its lowest mean value of 121 GPa (std. 8 GPa) when
PL = 225 W and
v = 1000 mm/s were used, and its highest mean value of 137 GPa (std. 3 GPa) when
PL = 250 W and
v = 1000 mm/s were used. Chen et al. reported a Young’s modulus value of 127 ± 4 GPa for an as-built specimen [
36], which is consistent with these results. There is a high interest in correlating mechanical properties from nano to macro scale. In that context, Tuninetti et al. [
62] have found a relationship based on which the flow stress at the macro scale can be estimated from the nano-hardness results of conventionally processed Ti6Al4V alloy. Furthermore, the flow stress can be related to nano-hardness using the Tabor relation as stated in [
55].
In order to estimate the Young’s modulus value based on nano-hardness (
H), or vice versa, a linear regression model was applied on nano-mechanical experimental data of all annealed (840 °C-2h-FC) specimens. It was found that a relation between
E and
H can be well represented using a simple linear model:
E = 15.066|
H| + 60.514 with
R2 = 0.744 (
Figure 4a). Moreover, the calculated correlation coefficient (
r = 0.863) indicates a strong correlation between
E and
H, as stated in [
63,
64].
Using this relation, the Young’s modulus can be estimated from nano-hardness data of L-PBF Ti6Al4V alloy utilizing laser power and scanning speed combinations in a range of 200–250 W and 1000–1500 mm/s, respectively. The proposed model can also be used as a reference for an E–H relation comparison between different heat treatment conditions, process parameters, manufacturing technologies, or even novel materials whose nano-mechanical properties are still rarely reported in the literature.
To verify the proposed model, a non-constant variance score test [
65] was performed and it was confirmed that the model has a homoscedastic variance of error term (
p = 0.715). This supports the assumption of equal variances, which are essential for the valid application of the proposed linear regression model. The residual plot (
Figure 4b) also confirms the applicability of a linear model, as random scattering is obvious. Furthermore, the Shapiro–Wilk normality test [
48] on studentized residuals was conducted to determine whether the model errors are normally distributed. The results show that there is no need to doubt the normality of model errors (
p = 0.378) which further validates the application of the proposed model.
3.3. Nano-Hardness Strain-Rate Sensitivity
Characterization of the elastic and plastic properties of the phases of polycrystalline materials is essential for determining the connection between microstructure and mechanical properties, especially for titanium alloys that have found application in highly demanding fields [
66]. The strain-rate sensitivity exponent is an important parameter for evaluation of the rate controlling mechanism during thermally activated deformations [
67], super-plasticity evaluation [
23], and crystal plasticity finite element modeling [
66]. Since L-PBF Ti6Al4V alloy in the as-built and annealed state has microstructural features (α’ martensite needles or α + β laths) significantly thinner than the size of the indenter tip, it is possible to characterize the material’s nano-hardness strain-rate sensitivity using a single
m value for a given load and heat treatment condition. To characterize the behavior of L-PBF Ti6Al4V alloy in the as-built and annealed state, as well as to determine the influence of different strain rates on the nano-hardness, a strain-rate sensitivity analysis was performed on nano-hardness–strain rate data.
The nano-hardness strain-rate sensitivity analysis showed that an indentation load of 10 mN leads to a higher variance of the material’s nano-hardness compared to results obtained for an indentation load of 200 mN. The reason is that the indentation depths were quite low when a load of 10 mN was applied (<300 nm), resulting in a higher data scatter, as shown in
Figure 5c,d. In this case too, more consistent results were again obtained at larger depths, i.e., when higher indentation loads were applied. Moreover, the strain-rate sensitivity exponent (
mi) for the D
ab and D
an specimens had lower values (0.010 and 0.017) when an indentation load of 200 mN was applied (
Figure 5a,b) compared to a 10 mN indentation load (0.053 and 0.040) as noticeable in
Figure 5c,d. Thus, nano-hardness was less sensitive to applied strain rates when higher indentation loads were used.
In all cases, the hardening effect of the strain rates on the nano-hardness is pronounced. The strain-rate sensitivity exponents (
mi) for the D
ab and D
an specimens when subjected to 10 mN indentation load were 0.053 and 0.040, respectively. When a 200 mN load was applied, the
mi values the for D
ab and D
an specimens were 0.010 and 0.017, respectively. The results given here are consistent with other published results considering
mi. For instance, Peng et al. determined
mi for electron-beam-melted Ti6Al4V alloy manufactured using different scanning strategies, as 0.053 ± 0.014 and 0.047 ± 0.009 [
23]. Jun et al. reported
mi for dual-phase Ti6Al2Sn4Zr2Mo alloy in a range from 0.005 to 0.039 [
20]. By calculating
mi as 0.056 and 0.064, Zhang et al. found that
mi is independent of grain orientation in the
β phase of Ti7Mo3Nb3Cr3Al alloy [
66].
3.4. Creep Behavior
Based on load–unload curves obtained using nanoindentation tests on the D
ab and D
an specimens it is evident that the width of the load plateaus increases with a higher indentation load (
Figure 6a,b). Since the load plateaus were in general wider for the D
an specimen (~58 nm when a 200 mN holding load was applied) compared to the D
ab specimen (~40 nm when a 200 mN holding load was applied), it indicates that the annealing heat treatment causes a lower creep resistance.
Figure 6c,d also confirm this finding since the curves of the D
an specimen have a higher increasing trend, compared to the curves of the D
ab specimen.
All indentation depth vs. time curves have a pronounced increasing trend both in transient and steady-stage creep regimes. It is evident that the creep displacement increases rapidly as the time increases at the initial stage, and then at later stage significantly slows down and retains an almost linearly increasing trend.
A least square fitting procedure was used within this work to fit Equation (8) and Equation (14) to the experimental data. The parameters
a,
b, and
k (
Table 7) were determined by fitting Equation (8) to the experimental data from the creep stage. For all applied holding loads, Equation (8) was fitted to the experimental data with high agreement. This was demonstrated by the D
an specimen subjected to a 200 mN indentation load (
Figure 7a), where the mean absolute error (MAE) was 0.331 nm. In
Figure 7a, a high creep strain rate dependence as a function of time is also shown, where a high strain-rate decrease specific to the transient creep regime can be seen. Furthermore, Equation (14) was also fitted with high agreement (MAE < 0.415 mN) to the unloading part of the curve (
Figure 7b). From the unloading part of the curve, the material constants
B and
m were derived and used to determine the contact stiffness
S (
Table 8) using Equation (15).
In their work, Pharr and Bolshakov discovered that
m was in a range of 1.2 and 1.6 for six different experimentally tested materials [
68], which is in accordance with our reported results in
Table 8. Using different indentation loads in this work, it was found that the
B and
m parameters depend on the maximum load applied, indicating that the curvature of the unloading curve also changes with applied load (see
Figure 6a,b). Consequently, the calculated contact stiffness values were also related to the applied load in such a way that the values increased when higher loads were applied (
Table 8). The material parameters
B and
m are still not available in the literature for Ti6Al4V alloy manufactured via L-PBF or they are scarcely reported, making these results even more valuable.
Furthermore, the slope of the representative
curves decreases rapidly as the creep approaches its steady state, as shown in
Figure 8a,b. By applying Equation (10) to the
data, it is possible to determine the creep stress exponent,
n, from which the dominant creep mechanism and the creep stability can be evaluated. When
n = 1 the diffusion creep mechanism is dominant,
n = 2 indicates that the grain boundary sliding mechanism is present, and
n > 3 indicates a dislocation movement as the dominant creep mechanism [
69,
70]. In
Figure 8c,d it can be seen that
n > 3 for both the D
an and D
ab specimens, which indicates that the creep deformation is governed by a dislocation movement.
In both cases, the highest applied indentation load during the creep stage resulted in the lowest data scatter of the calculated
n values (
Figure 8c,d). When lower indentation loads were applied, a higher data scatter of the calculated
n values was more prominent on both specimens (
Figure 8c,d). High mean
n values were found on both specimens, which is in correspondence with observations for a CoCrNi multi-principal element alloy that was also characterized by dislocation movement as the dominant creep mechanism [
28]. However, it is worth mentioning that there is some debate whether specific types of creep mechanisms actually exist, such as the Harper–Dorn diffusion creep [
71].
When mean
n values are compared for identical indentation loads (
Figure 8c,d), it is evident that in each case the as-built specimen has higher mean
n values than annealed specimen. In general, the mean
n values of the as-built specimen are translated more upwards, than the mean
n values of the annealed specimen. This indicates that the annealing heat treatment reduces the creep resistance of L-PBF Ti6Al4V alloy, which is in correspondence with the behavior shown in
Figure 6c,d. Although annealing heat treatment is beneficial for residual stress relaxation [
72], ductility increase [
73], and anisotropy moderation [
74], it also lowers the creep resistance, which is undesirable since Ti6Al4V alloys are often used in high temperature applications where creep deformation is present.