4.1 Macroscopic morphology analysis of ZrO2/Ti6Al4V MMC layer
Figure 4 gives the transverse sections of ZrO2/Ti6Al4V LMI specimens. The optical macroscopic morphologies obtained under the same laser scanning rate and different power numbered from case #1 to #12. It can be observed that the transverse cross-section morphologies of all samples are like lower lips. The thermal cracks are also present in all specimens, and the ZrO2p distribution density decreases gradually from surface to inside of the MMC layer except in case #6. Besides, there is porosity in sample #5.
The bottom of the molten pool is almost entirely Ti6Al4V metal, and ZrO2 particles are mainly distributed in the upper part and the edge area of the molten pool in case #1. The cracks in through-wall mode can be observed in cases #2 and #3, in which the linear energy input is 8.0×103 J·m− 1 and 8.67 ×103 J·m− 1, respectively. In comparison, lower laser power only causes micro-cracks near the interface between the MMC layer and Ti6Al4V substrate, as presented in cases #5 and #4. Furthermore, the higher the linear energy input, the lower the melting-mixing ratio between ZrO2p and Ti6Al4V. Moreover, case #6 shows that ZrO2 particles are uniformly distributed in the molten pool. Therefore, a larger linear energy input will not only increase the incidence of cracks but also reduce the implantation rate of ZrO2 particles.
The combination of the same laser power and lower scanning velocity causes a more remarkable crack numbered from cases #7 to #12. Extending with an acute angle led to spallation shapes near the track surface, as shown in cases #7, #9, and #10. The smaller scanning velocity probably results in higher laser linear energy, which is more than 7.0 ×103 J·m-1, leading to strong interaction between the laser beam and the Ti6Al4V liquid in the molten pool. The increase of laser scanning velocity can reduce the size of the hot crack, which is probably due to the smaller thickness of the ZrO2p injection layer, as shown in cases #10, #11, and #12.
4.2 Feature extraction of ZrO2/Ti6Al4V track
It has been reported that no reaction layer is formed in ZrO2/Ti interface, which means that the severe cracks produced during the LMI process are mainly caused by thermal stress mismatching between ZrO2 and Ti6Al4V. To reveal the residual stress distribution within the LMI specimen, the finite element modeling is conducted based on specific process parameters and the resultant MMCs morphologies, including the track dimensions and the volume fraction of implanted ZrO2p. One of the measurements is the modified circular segment in Ref [20], and the abridged general view as shown in Fig. 5(a).
The critical dimensions of the ZrO2p/Ti6Al4V MMC layer presented in Fig. 5(a) include the cross-section width (CW), the reinforcement height RH, and the cross-section depth (CD). The MMC layer depicted in Fig. 5(b) is divided into three sublayers where the volume fraction of ZrO2p decreases linearly from VF1 to VF3.
The bar graph in Fig. 6 exhibits different experimental parameters' cross-section width and depth. Figure 6(a) uncovers a principle that cross-section width and depth decrease with the reduction of linear energy input. Moreover, the injection depth decreased from 1.07 mm to 0.31 mm as the laser power decreased from 1.4 kW to 0.9 kW. However, the CW and CD of case #4 (laser power is 1.1kW and linear energy is 7.33 ×103 J·m− 1) are higher than case #3 (laser power is 1.2kW and linear energy is 8.0 ×103 J·m− 1). Comparing the through-wall cracks and distribution of ZrO2 particles, it can be inferred that excessive laser heat input leads to the instant gasification of Ti6Al4V metal, which takes away a large amount of energy and reduces the content of Ti6Al4V in the MMC layer resulting in a decrease of the CW and CD, and forming penetrating cracks. Inconspicuous gasification occurs when the linear energy is lower than 8.0 ×103 J·m− 1 because of the fewer metals above boiling. It also illustrates that the CW and CD decrease as the scanning rate increases in Fig. 6(b).
Due to the fewer cracks of cases #1, #5, and #8 in the MMC layer, the three LMI cases are modeled for FEM simulation to analyze the residual stress. The values of CW, RH, and CD of the LMI experiments are measured. The VF value in different sublayers is averaged from three different location transversal sections of samples, as shown in Table 5.
Table 5
The MMC layer transverse cross-section parameters measured in three LMI cases.
Case No. | CW/mm | CD/mm | RH/mm | VF1 | VF2 | VF3 | VF-average | VF-approx. |
#1 | 2.8 | 0.93 | 0.26 | 0.504 | 0.234 | 0.062 | 0.267 | 0.3 |
#5 | 2.47 | 0.44 | 0.21 | 0.652 | 0.486 | 0.432 | 0.523 | 0.5 |
#8 | 2.53 | 0.70 | 0.22 | 0.573 | 0.452 | 0.05 | 0.358 | 0.4 |
4.2.1 Heat source model verification
It is indispensable to check the heat source model due to the existence of ZrO2 particles in the MMC layer. Table 6 gives the specific MMC layer dimensions between experimental and simulated. These comparisons certify the simulated cross-section morphology of the molten pool is in good agreement with the experimental results under the process parameters given in Table 1.
Table 6
The MMC layer transverse cross-section dimensions between experimental and simulated.
LMI cases | Result types | CW (mm) | CD (mm) |
Case #1 | Simulation | 2.8 | 1.1 |
Experiment | 2.8 | 1.07 |
Case #5 | Simulation | 2.5 | 0.5 |
Experiment | 2.47 | 0.44 |
Case #8 | Simulation | 2.5 | 0.75 |
Experiment | 2.53 | 0.75 |
The temperature cycle curves at different distances from the centerline of the MMC layer are measured by a K-type thermocouple HP-DJ8X25. It is difficult to obtain an actual temperature history curve in the MMC layer center. The simulation can represent the experimental result of point C because the curves reveal the calculated thermal cycle curves of case #1 has an excellent agreement with the experiment in comparison, as shown in Fig. 8.
Figure 8. The temperature cycle curve resulted from experiment and simulation at (a) point A and point B; (b) the center of the MMC layer.
4.3 Analysis of macro residual stress characteristics
Due to the vast difference in the properties of Ti6Al4V and ZrO2p materials, considerable residual stress will be generated under the laser beam, which has a significant influence on crack formation. Therefore, it is necessary to understand the distributions of residual stress. According to the LMI parameters given in Table 1 and the mesh model in Fig. 2, the finite element model of thermal elastoplastic is established to analyze residual stress characteristics in the MMC layer.
4.3.1 The transversal residual stresses analysis of the MMC layer
The transversal residual stress (σxx) contours in Fig. 8 manifest the maximum compressive stresses are all distributed at the start and ends of the LMI progress, and tensile stresses are mainly distributed in the middle region of the MMC layer and where the transition region to Ti6Al4V substrate. Additionally, the maximum compressive stresses are always higher than the ultimate tensile stress for three different process parameters. Though the order of linear energy from large to small is: case 1, case 8, and case5, the stress legends expose the maximum tensile residual stresses of case #1, case #5, and case #8 are 344 MPa, 348.5 MPa, and 352.4 MPa respectively. And the maximum compressive residual stress of case #1, case #5, and #8 are 818.3 MPa, 1369.8 MPa, and 1180.2 MPa, respectively. Combined with the volume fraction of zirconia in Table 5, it can be thereby deduced that the peak residual stress has a more significant correlation with the volume fraction of zirconia than with the heat input.
Figure 9 (a) presents the σxx curves of the top, middle, and bottom layers of the MMC region from the start to the end. These curves reveal the top surface distributes transversal compressive stresses, but the middle and bottom areas disseminate tensile stress. The reason is the presence of zirconia particles in the MMC layer, and the upper surface is in a free state without any constraint. At the same time, the bottom region is constrained by the unmelted Ti6Al4V substrate. In addition, the σxx fluctuates between 2mm and 4mm from the starting point of the laser beam. That is because of the faster melting rate of Ti6Al4V under the action of the laser beam. However, in the initial stage of LMI, the initial temperature of the substrate is room temperature, which leads to a fast thermal conduction rate and rapid solidification of the molten metal.
Figure 9 (b) shows the σxx values from MMC upper surface to the bottom region of the Ti6Al4V in the thickness direction. Their specific perpendiculars are marked in Fig. 8. The distances from the vertical lines L-A, L-B, and L-C to the starting point are L/2, in which L is the length of the Ti6Al4V substrate. These curves disclose that the transverse stress exhibits a compression-tension-compression-tension distribution in the thickness direction. The smaller the volume fraction of zirconia particles, the greater the compressive stress on the upper surface, as shown in the L-A line. Moreover, the maximum tensile stress is mainly distributed in the heat-affected zone of the Ti6Al4V substrate. And the maximum compressive stress is primarily concentrated in the base metal area under the Ti6Al4V heat-affected zone. The three curves show that although case #1 LMI has the largest heat input, its maximum compressive stress is 221.5 MPa, which is smaller than 245.6 MPa for cases #8 and 263.4 MPa for case #5. However, case #5, with the largest zirconia volume fraction, has maximum tensile stress of 328.6 MPa, greater than that of cases #1 and #8. Comparing curves L-C and L-B, it can be inferred that the effect of decreasing heat input on the transverse residual stress is far less than that of increasing the volume fraction of zirconia. Therefore, it concluded that the main factor affecting the macroscopic residual transverse stress is the volume fraction of zirconia particles. Further, two crests of the L-A line indicate that inhomogeneous dispersion of ZrO2p and the remelted Ti6Al4V have more significant influences on the distribution and magnitude of residual stress, combined with the optical topography analysis.
4.3.2 The longitudinal residual stress analysis of the MMC layer
Figure 10 makes clear the longitudinal tensile stress is in the MMC layer, and the compressive stress distributes other sections. The maximum tensile stress of three LMI cases is more than 1800 MPa, which is much higher than the yield strength of Ti6Al4V titanium alloy. This phenomenon may be the vast difference in the properties between Ti6Al4V and ZrO2p materials and the FEM model assumption in section 3.1. In addition, the most considerable tensile stress is always higher than maximum compressive stresses. The peak tensile stress of case #8 is higher than the other two LMI cases, reaching 1811.7 MPa. However, the tensile stresses of case #5 with minimum linear energy are higher 246.4 MPa than case #1 with maximum heat input. It indicates that the effect of increasing the volume fraction of zirconia on the residual longitudinal stress is much greater than that caused by the decrease in heat input.
Figure 11(a) uncovers the fact that the residual longitudinal stress (σzz) at the start and end of the MMC layer on the top surface are close to 0 MPa, and the tensile stress distributes the entire MMC layer since the endpoints are free boundaries resulting in no constraints beyond the endpoints. In addition, the σzz on the upper surface of a whole MMC layer is always larger than that in the middle region and then larger than the bottom area. The peak residual stress of case #1 is the smallest, which does not exceed 1500 MPa.
Figure 11 (b) reflects the σzz variation regulation from the MMC upper surface to the bottom surface of the Ti6Al4V substrate along the centerline in the thickness direction. These curves show the tensile stresses gradually decrease from the upper surface to the inside of the MMC layer. However, there are trough points near the fusion lines, and there are stable transition zones in the heat-affected zone. The maximum tensile stresses are located on the upper surface of the MMC layer, and minimum tensile stresses are lay on the bottom surface of TI6Al4V. Moreover, the largest compressive stresses are concentrated in the base metal under the heat-affected zone. In addition, the σzz of the entire MMC layer exceeds the yield strength of Ti6Al4V, which indicates that the presence of zirconia particles has a significant influence on the increase of residual stress, which is also the reason why the MMC layer is susceptible to cracks. The three curves also indicate that the larger the volume fraction of zirconia is, the larger the residual stress is near the fusion line.
4.4 Analysis of microcosmic residual stress characteristics
For different LMI macro models, the implant layers are layered according to the different content of ZrO2 ceramic particles. The two-dimension microscopic model is designed to study the interaction between ZrO2p particles and Ti6Al4V substrate where the computation domain is 250×250 um2. The diameter of ZrO2 ceramic particles is assumed to be 80 um, and they are randomly distributed in Ti6Al4V titanium alloy.
Using the center (the Point C in Fig. 8) of the MMC layer thermal cycle curve (as shown in Fig. 12) as the thermal boundary conditions of the micro model, the contents of ZrO2 ceramic particles are 30%, 40%, and 50%, as shown in Fig. 13. The lower-left corner of the calculation domain is the origin coordinates.
Due to the different contents of ZrO2p resulting in different distances between particles from Fig. 14, the residual stress of the substrate located between particles is complicated. The micro transversal stress (m-σxx) and longitudinal stress (m-σyy) in particles are more than the Ti6Al4V substrate. And the stress at the edge of the particle is the least. Even though the volume fractions of ZrO2 are different in the computational domain, the residual stress is uniformly distributed in the ZrO2 particles. After that, the average residual stress at various locations in the computational domain is extracted, as depicted in Fig. 15.
Four particles are distributed in a rhombus for the model with a 30% volume fraction of ZrO2, and particle A4 is at the edge of the computational domain. The m-σxx inside the particles is almost equal and around the particles is less than the Ti6Al4V substrate. The average m-σxx at region P1 is 781MPa, which is higher than region P2, which is 689 MPa. In contrast, the average m-σyy at region P1 is 790 MPa, which is less than region P2, which is 861 MPa. This is because P1 is surrounded by particles A1, A2, A3, and A4, so it is affected by four particles. Besides, the P2, just between particles A2 and A3, is only affected by two particles. The average m-σxx and m-σyy of particles A1, A2, and A3 are close to 930 MPa, as shown in Fig. 15 (b).
Five particles are distributed in a centrosymmetric manner in a model with the 40% volume fraction of ZrO2p. Therefore, areas P3, P4, P5, and P6 are influenced severally by three particles resulting in the m-σxx and m-σyy at the symmetric region (P3-P5 and P4-P6) being similar. The m-σxx and m-σyy of P3 and P5 are about 780 MPa and 800 MPa, respectively. The m-σxx and m-σyy values of P4 and P6 are about 803 MPa and 774 MPa, respectively. Furthermore, the phenomenon of low m-σxx around the particle is more noticeable. Still, the low m-σxx areas between particle B3 and other particles are narrower than the rest of the region around the particle. What’s more, the average m-σxx and m-σyy of these five particles are close to A1.
There are six particles randomly distributed in a model with a 50% volume fraction of ZrO2p. The average m-σxx and m-σyy in every particle are severally about 925MPa and 930 MPa, which access A1. However, there are also low transversal stress areas between particles, such as points P9 and P10, as shown in Fig. 14. The m-σxx of P9 and P10 are 590MPa and 561 MPa, respectively. The m-σyy of P9 and P10 are 886 MPa and 891 MPa, respectively. However, there is an interesting occurrence in which the micro residual stresses of P7 are similar to P3, and that of P8 is close to P4. That is because the positions of P7 and P8 are severally identical to those of P3 and P4 in the computational domain, influenced by the three surrounding particles.
It concludes that the micro residual stresses between particles are affected by the particle spacing. The distance is more adjacent between the particles in the direction of the coordinate axis, the more significant residual stress is along this direction.
4.5 Analysis of crack formation mechanism in MMC layer
It can be inferred that the accumulation of micro residual stress leads to more enormous macro residual stress, leading to the formation of cracks. The ZrO2 is a ceramic material, in which the melting point (2700℃) is much higher than Ti6Al4V. Therefore, the melted Ti6Al4V metal solution will coat the ZrO2, as shown in Fig. 16(a). And the thermal expansion coefficient of ZrO2 particles hardly changes during the process, but that of Ti6Al4V is greatly affected by temperature. Thereby, a gap, as shown in Fig. 16(b), is generated between ZrO2 particles and Ti6Al4V substrate when the Ti6Al4V metal solidifies, because of the vast difference in the two material properties resulting in the metal surrounding the particle being partially under tensile stress and partially under compressive stress. A mass of ZrO2 particles is equivalent to many inclusions present in Ti6Al4V metal. During the solidification process of Ti6Al4V, the metal dislocations move to contact with the ZrO2 particles resulting in a stress concentration effect at a point (SCP), thereby transferring the stresses to the ZrO2 particles. That is why the residual stresses in the ZrO2 particle are more prominent than that in the Ti6Al4V substrate. A micro fissure will be formed when the stress at this point exceeds the yield stress limit of Ti6Al4V, and the gap will be further expanded, as shown in Fig. 16(c). As the cooling goes on, the initial micro- fissures gradually develop under the action residual stresses until the crack forms, and two gaps also merge, as shown in Fig. 16(d). Finally, several micro-cracks and gaps grow, forming cracks visible to the naked eye.