4.1. Experimental Results
Figure 3 illustrates the stress–strain curves of the Mg-13Gd-4Y-2Zn-0.5Zr alloy under various deformation conditions. The alloy displays the characteristic behaviors of work hardening and dynamic softening during hot processing. Initially, during its deformation, the flow stress experienced a sharp increase with the strain, nearly linearly, owing to the work hardening generated by material deformation. Deformation, dislocation climb, and cross slip prompted the dynamic recovery of the material, subsequently reducing the increasing rate of flow stress. Upon reaching the peak stress, the flow stress began to decline, as the substantial accumulated plastic strain energy within the material was released through phenomena like dynamic recrystallization, marking a phase where strain-softening mechanisms, including dynamic recrystallization, reduce the flow stress. Consequently, the entire hot compression curve exhibited a rapid rise followed by a gradual decline. Comparing the curve results between different temperatures and strain rates, it is evident that the temperature and strain rate directly influence the magnitude of peak stress and the sequence of its occurrence. This is mainly because the coupled effects of different deformation conditions influenced the duration and extent of softening mechanisms. At higher strain rates, the duration of softening is shorter, resulting in a lower softening extent and higher peak stress.
When the deformation temperature remained constant, the flow stress increased with the strain rate for a given deformation level, and the peak stress was observed at a later stage. This behavior is attributed to the shorter durations of processes like dynamic recovery and dynamic recrystallization at higher strain rates, resulting in reduced softening effects and higher flow stress.
Conversely, at a constant strain rate, the flow stress decreased as the temperature rose for the same deformation level. This trend is due to the stronger thermal activation effect at higher deformation temperatures, leading to an increased atomic diffusion rate. Consequently, more dislocation climb and cross slip occurred during the deformation, enhancing the softening effect of the dynamic recovery and reducing the flow stress. Moreover, the recrystallization process, involving nucleation and growth, was influenced by temperature. Higher temperatures promoted the nucleation of the deformed structure, intensifying the softening effect induced by dynamic recrystallization and resulting in a decline in the flow stress curve.
In summary, the flow stress curve of the Mg-13Gd-4Y-2Zn-0.5Zr alloy exhibits typical characteristics of a dynamic recrystallization-type curve. The thermal processing parameters, such as the deformation temperature and strain rate, significantly influenced the dynamic recrystallization behavior of the alloy, with higher temperatures and lower strain rates favoring dynamic recrystallization [
4].
The impact of the strain rate on the dynamic recrystallization microstructural evolution is depicted in
Figure 4. During hot processing, the strain rate primarily modulated the nucleation frequency, (i.e., the occurrence rate of dynamic recrystallization throughout the entire hot deformation process) and the duration of grain growth [
22], thus refining the microstructure. It is apparent that, at lower strain rates, a greater number of recrystallized grains emerged along the original grain boundaries. At a constant deformation temperature (470 °C), a lower strain rate corresponded to a smaller average grain size and a higher recrystallized fraction, as depicted in
Table 2. This phenomenon arose from the reduced deformation time at higher strain rates, resulting in fewer nucleation events and the inadequate growth of dynamic recrystallized grains within shorter growth cycles, thereby reducing the proportion of dynamic recrystallized grains relative to the overall microstructure. Conversely, as the strain rate decreased, a prolonged duration was necessary to complete the deformation process under the same deformation level. Given a constant nucleation probability, this led to an increased number of recrystallization grains. Concurrently, the initial recrystallized grains underwent prolonged growth periods, resulting in coarser recrystallized grains at lower strain rates. However, from a statistical standpoint, the abundance of new dynamic recrystallization grain nuclei tends to yield smaller average grain sizes at lower strain rates.
Temperature plays a pivotal role in influencing the microstructural evolution during hot deformation processes. At higher temperatures, the activation energy of atoms changes, making it easier for them to transition into new states, thereby influencing the nucleation and growth processes.
Figure 5 delineates the dynamic recrystallization microstructural morphology of the Mg-13Gd-4Y-2Zn-0.5Zr alloy compression specimens subjected to 60% deformation and a strain rate of 0.001 S
−1 at various deformation temperatures.
Table 3 presents the disparities in the average grain size and dynamic recrystallization fractions among the different temperatures. As depicted in
Figure 5, it is evident that, at 470 °C, there was a significant increase in the number of newly formed dynamic recrystallized grains, accompanied by an elevated dynamic recrystallization fraction. Additionally, visually, the grain size appeared to be larger. This phenomenon can be attributed to several factors relevant to nucleation processes. Firstly, elevated temperatures prompt heightened atomic instability, facilitating the surmounting of energy barriers, and thereby promoting a greater frequency of nucleation events for recrystallization. Furthermore, the prevailing viewpoint [
23] suggests that the gestation cycle is primarily correlated with changes in the dislocation density within the matrix. With increasing temperatures, the activity of substructures represented by dislocations intensifies, facilitating climb and annihilation, and thereby reducing the rate of dislocation accumulation. Consequently, elevated temperatures can shorten the gestation period necessary for dynamic recrystallization in materials while simultaneously increasing the frequency of nucleation events for recrystallization. When it comes to the growth of grains, higher temperatures induce grain coarsening, while the enlargement of recrystallized grains tends to overlay the initial matrix structure. Considering these factors collectively, it is evident that higher temperatures lead to microstructures with higher dynamic recrystallization fractions, while the increased number of newly formed recrystallized grains contributes to smaller average grain sizes.
4.2. Simulation Results and Discussion
DEFORM-3D facilitates users in defining the plastic flow behavior of materials by inputting rheological stress data under various process conditions. In this investigation, a rheological stress model
was selected. Utilizing available experimental data, multiple data points at different temperatures and deformation rates were carefully selected and evenly distributed to ensure precision. These data points were then entered into the software’s data table, and the software automatically interpolated them to generate stress–strain curves. Subsequently, upon accessing the sub-menu, each physical property parameter of the material was systematically and meticulously included and stored, thereby establishing the alloy’s material repository. Through analysis of the results of the hot compression experiments, the various models mentioned earlier were calculated. Furthermore, the user subroutine of DEFORM-3D underwent further development, wherein the relevant models were incorporated into the software’s subroutine using the FORTRAN 11.0 language. The computational procedure of the dynamic recrystallization subroutine within the software is depicted in
Figure 6b.
Following the secondary development, simulations of microstructural evolution during hot deformation were conducted using the CA module of DEFORM software. This model facilitated the quantitative analysis and dynamic visualization of microstructural changes during the high-temperature deformation. Employing periodic boundary conditions, the CA model in DEFORM represents an infinite space, with neighbor types determined using the Moore eight-neighbor criterion. The selected area was subdivided into quadrilateral cells measuring 210 rows by 210 columns, each with a side length of 2 μm, effectively representing an actual sample area of 420 μm by 420 μm. Throughout the simulation of microstructural evolution, the initial dislocation density of the cells was set to 0.01 μm−2. Across the entire hot deformation stage, the dislocation density within the cells fluctuated with the level of deformation, conforming to the classical KM equation. Cells positioned at the grain boundaries underwent recrystallization nucleation only under specific conditions. Following nucleation, the dislocation density of the cells reverted to its initial value, and the newly nucleated cells grew according to the previously described recrystallized grain growth model. With an extended simulation time, the recrystallized grains persisted in nucleating and growing, gradually encasing the original matrix structure. The dislocation density within the newly formed recrystallized grains continued to evolve in accordance with the KM model.
During deformation, when the critical strain is reached or exceeded, nucleation occurs randomly at the grain boundaries of the parent phase. The newly formed nuclei grow outward through grain boundary migration, consuming the parent phase. Simultaneously, nucleation and grain growth continue, resulting in a chain-like distribution of new recrystallized grains surrounding the parent phase. This process leads to the formation of new grain boundaries between the parent phase and the newly formed recrystallized grains, as well as among different recrystallized grains. When the predetermined strain threshold is reached, the dynamic recrystallization process concludes, with most of the parent phase being replaced by smaller dynamic recrystallized grains. The dynamic recrystallization process was vividly portrayed in the cellular automaton simulations, closely aligning with the foundational principles of classical theoretical descriptions.
During the thermal compression simulation, cylindrical specimens matching the specifications (Φ8 × 12 mm) of physical experiments were utilized, ensuring complete consistency in the deformation conditions. The billet model was defined as a plastic body, with a finite element mesh comprising 50,000 elements, while both upper and lower dies were set as rigid bodies. Owing to the non-uniformity of deformation, resulting in varied equivalent strains across specimen regions, representative points P1, P2, and P3 (as depicted in
Figure 6a) were initially chosen for analysis. Subsequently, microstructural evolution was simulated based on these analyses.
Figure 6c illustrates the microstructure generated using the CA method, exhibiting a grain distribution pattern largely in line with the actual scenario, with notable agreement in grain morphology. Statistical analyses based on electron backscatter diffraction (EBSD) observations determined the initial average grain size of the parent phase to be 85 μm for simulation purposes.
From the evolution curves and contour maps of the equivalent strain (
Figure 7), it is evident that the distribution of the equivalent strain in the high-temperature alloy billet at the end of hot compression varied significantly across different locations. The following three primary regions emerged: the central zone displayed the highest equivalent strain, while the bulging area exhibited a lower equivalent strain compared to the central region, with the regions in contact with the upper and lower dies showing the least equivalent strain. The evolution curves further demonstrate a progressive rise in the equivalent strain at points P1, P2, and P3 over the deformation time. However, due to its location in a difficult deformation zone and its contact with the upper die, point P2 experienced frictional forces during deformation, hindering radial metal flow. Mechanically, P2 was subjected to compressive stress, posing difficulties in meeting the plastic deformation requisites. Conversely, point P1, located in a heavily deformed area and experiencing minimal frictional forces, was more conducive to plastic deformation compared to P2, and thus, displayed the highest equivalent strain. Point P3, situated in a zone of unrestricted deformation, exhibited an equivalent strain between that of points P1 and P2.
The hot deformation simulation results indicate that P1 was situated in a highly deformed region, P2 in a difficult deformation zone, and P3 in a region of free deformation. Consequently, for the dynamic recrystallization evolution simulation, point P1, corresponding to the consistently observed strain levels in the experiments, was selected.
During deformation, once the strain reaches or exceeds the critical strain, nucleation occurs randomly at the grain boundaries of the parent phase microstructure. The newly formed nuclei grow and consume the parent phase by migrating along the grain boundaries; this nucleation process proceeds concurrently with grain growth. Observations reveal that the newly recrystallized grains exhibit a chain-like distribution surrounding the parent phase. Consequently, new grain boundaries form between the parent phase and the newly recrystallized grains, as well as among different recrystallized grains. When the prescribed strain level is reached, the dynamic recrystallization process concludes, with the majority of the parent phase replaced by smaller dynamically recrystallized grains. The results of cellular automaton simulations vividly replicate the dynamic recrystallization process, closely aligning with classical theoretical descriptions.
Figure 8 demonstrates that, at a temperature of 470 °C, an increase in the strain rate led to a reduction in the size of the dynamically recrystallized grains. As the deformation rate escalated from 0.001 S
−1 to 1 S
−1, the size of the recrystallized grains diminished from 14.6 μm to 3.69 μm, accompanied by a decrease in the volume fraction of recrystallization from 91.72% to 56.63% (as depicted in Figures 10 and 11). Lower strain rates, under consistent deformation temperature conditions, resulted in diminished stored energy and driving force within the alloy. Lower strain rates, under consistent deformation temperature conditions, resulted in diminished stored energy and driving force within the alloy. Consequently, the nucleation of alloy grains occurred sporadically in areas with energy fluctuations, leading to a decreased nucleation rate. Moreover, ample time was afforded for grain boundary migration, facilitating the growth of recrystallized nuclei, and thus, augmenting the volume fraction of recrystallization and relatively enlarging the recrystallized grain size. Conversely, higher strain rates prompted rapid dislocation accumulation and insufficient stress concentration release, thereby impeding the growth of recrystallized grains. Additionally, an increase in the strain rate resulted in a larger peak strain, corresponding to a larger critical strain, which made it more difficult for dynamic recrystallization to initiate and slowed down the process. Furthermore, the reduction in the recrystallized grain size resulted in a decreased overall volume of the recrystallized region. Consequently, both the size of the dynamically recrystallized grains and the volume fraction of recrystallization were diminished. In conclusion, while increasing the strain rate decelerated the rate of recrystallization, it enhanced the refinement effect of recrystallization on the microstructure.
Figure 9 illustrates the microstructure after deformation termination under different deformation temperatures with a strain rate of 0.001 S
−1. The figure clearly demonstrates that, with an ascent in the deformation temperature from 430 °C to 470 °C, the dimensions of the dynamically recrystallized grains escalated from 2.31 μm to 14.6 μm (
Figure 10), concomitant with a marked augmentation in the recrystallization volume fraction (
Figure 11). This phenomenon arose from the ratio of the nucleation rate to the growth rate, denoted as N/G, dictating the dimensions of the dynamically recrystallized grains during the thermal deformation [
23]. It is widely acknowledged that heightened deformation temperatures expedite the growth rate of dynamically recrystallized grains, fostering a more exhaustive recrystallization process, and thereby augmenting the volume fraction of recrystallization. At the stage of stable deformation under elevated temperatures, the comparatively reduced dislocation density within the grains engendered a diminished nucleation rate, with the nucleation rate, N, declining while the growth rate, G, escalated. Consequently, the N/G ratio diminished, fostering an augmentation in the dimensions of the dynamically recrystallized grains. Furthermore, at elevated deformation temperatures, vigorous atomic oscillation in high-temperature alloys accelerate atomic diffusion, significantly heightening the migratory capacity of the dislocations and grain boundaries. The diminished strain rate allows ample time for this phenomenon, facilitating a more comprehensive and uniform growth of dynamically recrystallized grains over the same temporal span.
The preceding analysis indicates that increases in the deformation temperature led to a progressive enlargement in the size of the dynamically recrystallized grains. Conversely, with the escalation of the deformation rate, the dimensions of the dynamically recrystallized grains gradually decreased. Consequently, elevated temperatures in conjunction with low strain rates foster the advancement of grain growth in dynamic recrystallization, whereas low temperatures coupled with high strain rates hinder the expansion of dynamically recrystallized grains.
The average grain size is also highly sensitive to both the deformation temperature and strain rate. However, their respective influences on the average grain size were rather intricate (
Figure 12). In the absence of complete dynamic recrystallization, as the deformation temperature rose, the average grain size gradually decreased, while, with increasing strain rate, the average grain size gradually increased. Conversely, in cases where complete dynamic recrystallization occurred partially, multiple scenarios of average grain size variations were observed with changes in the deformation temperature and strain rate. It is worth noting that
Figure 10,
Figure 11 and
Figure 12 delineate the comparison between the simulated and experimental values of the dynamic recrystallized grain size, dynamic recrystallization volume fractions, and average grain size under various processing conditions. The figures reveal a maximum disparity of 12.39% and a minimum disparity of 6.06%. This finding emphasizes the significant concordance between the simulated outcomes and the experimental data, affirming the robustness of the cellular automaton (CA) model devised to depict the microstructural changes in rare-earth magnesium alloys during thermal deformation. This model was refined through the secondary development of the user subroutine in the finite element software DEFORM-3D.