3.1. Analysis of Deformation Mechanism of {}<> Nanotwinned Titanium and Effect of Twin Boundary Spacing
We performed tensile and compression loading on the models with twin thicknesses of 6.9, 11.1, 13.8, and 27.6 nm, and their stress responses were recorded. Furthermore, the evolution process of the microstructure of the material during the simulation process was observed, and the plastic deformation mechanism of {}<> nanotwinned titanium under tensile and compression loading, as well as the impact of the twin boundary spacing on the plastic deformation of the model, was revealed.
Figure 3 shows the stress–strain curves of the models with different twin boundary spacings under tensile loading. All the tensile stress–strain curves exhibited initial linear elastic deformation and reached a peak value, and then the flow stress suddenly decreased. This type of decline indicated that the model began to undergo plastic deformation. With increasing strain, the flow stress exhibited large fluctuations composed of peaks and troughs, which are typical characteristics of discrete plastic deformation in nanopillars [
36]. As can be seen from the figure, the yield strengths of the models with four spacings ranged approximately from 2.4 to 2.8 GPa. Within the range of 6.0–27.6 nm, a smaller twin boundary resulted in a higher yield strength. The plastic stress of the model entered an oscillation stage after yielding. The largest magnitude of the oscillation was observed when the twin boundary spacing was 13.8 nm. At this moment, the mean flow stresses under different twin boundary spacings were not considerably different and the mean flow stresses of the four models were approximately in the range of 1.0–1.2 GPa.
To reveal the evolution process of the microstructure corresponding to the stress response of this type of nanotwinned titanium model, we selected three models with the twin boundary spacings of 13.8, 11.1, and 6.9 nm. The evolution processes of their atomic configurations with increasing tensile strain were carefully observed, and their plastic deformation mechanisms were analyzed in detail. As shown in
Figure 4, in the elastic stage of the model, partial dislocation nucleation on the (
)[
] basal plane was observed at the twin boundaries of the three models (strains of 2.88–3.6% in
Figure 4a–c). When the strain was increased to the yield point of the material (strains of 6.48–7.2% in
Figure 4a–c), a partial dislocation nucleated on the twin boundary slipped to the adjacent twin boundary, and the path after slip became a basal stacking fault. The detailed atomic graph is shown in
Figure 5a. It should be noted that when basal dislocations finally reached the adjacent twin boundary, dislocations were hindered by the twin boundary and reacted with it to form a base/cylinder interface (
Figure 5b) so that the original smooth twin boundary became a broken line. As shown in
Figure 4a–c, the strain was 14%. With increasing strain, the number of stacking faults caused by the basal faults continuously increased, and finally, the stacking faults merged with each other to complete a face-centered cubic phase transition. Subsequently, the broken-line-shaped twin boundary migrated, finally leading to the disappearance of part of the twin boundary, and the material underwent a face-centered cubic phase transition.
The deformations of the three models were largely the same in the elastic and plastic stages, so the differences in the stress–strain curves of all the nanotwinned titanium models in
Figure 3 were not significant. From the in-depth analysis, it can be seen that the slip movement of a large number of nucleated basal partial dislocations on the twin boundary led the model to yield. This type of plastic deformation mechanism resulted in a decline in the stress, i.e., strain softening. When the adjacent twin boundary hindered dislocations, the stress increased to some extent, i.e., strain hardening. Therefore, it was an alternating process of the emission of a large number of partial dislocations on the basal plane and blocking by twin boundaries that caused the flow stress oscillations of the model in the plastic stage.
Figure 6 displays the stress–strain curves of the models with different twin boundary spacings under compression loading. The yield strengths of the models with four spacings ranged approximately from 2.6 to 3.6 GPa. The yield strength of a model with the same spacing was larger under compressive loading than under tensile loading, and the stress after yielding decreased to around 0.7 GPa. The model with a spacing of 13.8 nm exhibited the highest yield strength. That is, there was a critical twin spacing in the range of 6.9–27.6 nm, and it maximized the model yield strength. When the strain increased to 8%, the model slowly entered the strengthening stage, and the stress of the model with a twin boundary spacing of 11.1 nm could even reach 6.4 GPa. In addition, in the range of 6.9–27.6 nm, the variation trends of the stresses of the four models with increasing strain were generally consistent.
Similarly, to investigate the evolution process of the microstructure corresponding to the stress response of the {
}<
> nanotwinned titanium model under compression loading, we selected three models with twin boundary spacings of 27.6, 11.1, and 6.9 nm. The evolution processes of their atomic configurations with increasing compression strain were carefully observed, and their plastic deformation mechanisms were analyzed in detail. As shown in
Figure 7, in the elastic stage, a large number of nucleated basal partial dislocations appeared on the twin boundaries of the three models. When the strain increased to 2.16–2.88%, the model yielded when the partial dislocations on the twin boundary were emitted, and the model entered the plastic deformation stage. When the strain further increased to 4.32–5.04%, it could be clearly seen that other than the emission of the partial dislocations on the twin boundary, there were numerous partial dislocations on the grain boundary. The path after they slipped became a basal stacking fault. Additionally, the twin boundary showed a blocking effect on the basal partial dislocations. This phenomenon was approximately the same as that of {
}<
> nanotwinned titanium under tensile loading. However, the difference was that, in addition to this deformation mechanism, apparent {
} twin boundary migration could be observed in the compression model. In particular, this was most evident in the model with a small spacing (
Figure 7c), and there was also a detwinning phenomenon caused by twin boundary migration (
Figure 8). When the strain increased to 7.92%, the dislocation slip and twin boundary migration continued to intensify. Finally, part of the twin boundary disappeared, and the material underwent a face-centered cubic phase transition.
The deformations of the three models in the elastic and plastic stages were approximately the same, so the differences in the stress–strain curves of the models within the twin boundary range of 6.9–27.6 nm were small. Under the compression loading, similarly, the emission of the nucleated basal partial dislocations on the twin boundary caused the model to yield. The deformation mechanism of the model in the initial plastic stage was dominated by dislocation slip. At this moment, the stress of the model was reduced to some extent. With increasing strain, the deformation mechanism of the model was dominated by the emission of basal impartial dislocations and the migration of twin boundaries. Because the twin boundary migration caused an increase in the flow stress, the model entered the strain hardening stage.
By comparing the above results, it can be seen that {
}<
> nanotwinned titanium exhibited tension–compression asymmetry, and this type of stress response under different loadings was caused by the different plastic deformation mechanisms. Under tensile loading, the plastic deformation mechanism of {
}<
> nanotwinned titanium was dominated by the slip of basal partial dislocations. Under compression loading, the plastic deformation mechanism of {
}<
> nanotwinned titanium changed and was dominated by the slip of basal impartial dislocations and twin boundary migration. This result was similar to the related research result of the plastic deformation mechanism of {
}<
> nanotwinned magnesium [
37]. Apparently, twin boundary migration is more favorable for the enhancement of material plasticity and can allow the model to exhibit a strain hardening stage after yielding.
3.2. Analysis of Deformation Mechanism of {}<> Nanotwinned Titanium and Effect of Twin Boundary Spacing
Similar to the tests described in the previous section, we performed tensile and compression loading on models with twin thicknesses of 6.7, 10.8, 13.5, and 26.9 nm, and their stress responses were recorded. Furthermore, the evolution process of the microstructure of the material during the simulation process was observed, and the plastic deformation mechanism of {}<> nanotwinned titanium under tensile and compression loading, as well as the impact of the twin boundary spacing on the plastic deformation of the model, was revealed.
Figure 9 shows the stress–strain curves of the models with different twin boundaries under tensile loading. The yield strengths of the models with four spacings were approximately in the range of 8.5–9.5 GPa. The model with a twin boundary spacing of 13.5 nm exhibited the highest yield strength. When the strain was 6.48–6.84%, the model reached the yield point and the stress rapidly dropped afterward. Subsequently, with increasing strain, the stress entered a stage in which the overall state decreased but the local flow stress still oscillated. The model with a twin boundary spacing of 10.8 nm exhibited the largest drop in stress. According to the previous results [
38], the (
) twin boundary had the lowest grain boundary energy of the three twin boundaries, (
), (
), and (
). This ensured that the (
) twin boundary had relatively high stability. Therefore, it can be speculated that it was difficult to have a large number of defect nucleation events on the twin boundary during the elastic stage, thus causing insufficient defects inside the model and highly concentrated stress. Compared to the {
}<
> nanotwin, the {
}<
> nanotwin entered the plastic stage even later (strains of 2.88–3.6% in
Figure 4a–c). Subsequently, the deformation was fast, and the plastic ability was even poorer.
To reveal the stress response of this type of nanotwinned titanium model, we selected two models with twin boundary spacings of 13.5 and 10.8 nm. The evolution processes of their atomic configurations with increasing tensile strain were carefully observed, and their plastic deformation mechanisms were analyzed in detail. In the elastic stage, a small number of partial dislocation nucleations on the (
)[
] basal plane were observed at the twin boundaries of the two models (strain of 6.48% in
Figure 10a,b). When the strain increased to 7.56%, various defects in the model, including a substantial number of basal partial dislocations, a small number of conical and cylindrical dislocations (
Figure 11), and cracks at the lower left of the model (strain of 7.56% in
Figure 10a) and the upper left of the model (strain of 7.56% in
Figure 10b), appeared. As the strain further increased, the stacking faults caused by the basal dislocations gradually increased and finally merged with each other to form a face-centered cubic structure area, and the degree of the face-centered cubic phase transition under the same strain was lower than that of the {
}<
> nanotwin. Furthermore, the cracks that initially formed could further propagate along the twin boundary (strain of 14% in
Figure 10a,b). The propagation rate of the cracks of the model with a twin boundary spacing of 10.8 nm was relatively fast, and the stress concentration was released rapidly, resulting in a fast decline in stress.
The deformations of the two models were approximately the same in the elastic and plastic stages, so the differences in the stress–strain curves of the models within the twin boundary range of 6.7–26.9 nm were small. Under tensile loading, the deformation mechanism of the model in the initial plastic stage was dominated by a small number of basal dislocation slips. With increasing strain, the basal, conical, and cylindrical partial dislocations were constantly emitted as crack nucleation and growth occurred. Meanwhile, the appearance of many defects allowed the model to release the stress concentration in the elastic stage and coordinate plastic deformation, which also led to a dramatic stress decline in the model after yielding. According to the above analysis, {}<> nanotwinned titanium possessed a high yield strength, but it had insufficient plastic toughness. This was because of the low grain boundary energy in its twin boundaries. The twin boundary structure was very stable and could not act as an effective source of defect nucleation. Therefore, in the {}<> nanotwinned titanium model, it was difficult for dislocation slip or twin-coordinating strain to occur at low strain levels. Once the strain increased to 7.56%, a large number of defects appeared simultaneously. This special property severely weakened its plasticity.
Figure 12 displays the stress–strain curves of the models with different twin boundary spacings under compression loading. The yield strengths of the models with four spacings were in the approximate range of 3.7–5 GPa. The highest yield strength was found in the model with a twin boundary spacing of 6.7 nm, which was far smaller than that of the model with the same spacing under tensile loading. The plastic stress of the model after yielding entered a stage in which the overall state was stable but the local flow stress still oscillated. The model with a small spacing showed small oscillations, while large oscillations were observed in the model with a large spacing. The flow stress finally tended to stabilize, and the plasticity was even stronger than that under tensile loading.
Similarly, to investigate the stress response of this type of nanotwinned titanium model, we selected two models with twin boundary spacings of 13.5 and 6.7 nm. The evolution processes of their atomic configurations with increasing compression strain were carefully observed, and their plastic deformation mechanisms were analyzed in detail. As shown in
Figure 13, similar to the elastic stage, partial dislocation nucleations on the (
)[
] basal plane were observed at the twin boundaries and grain boundaries of the three models (strain of 4.32% in
Figure 13a,b). As the strain increased (strain of 6.84% in
Figure 13a,b), when the yield point of the material was reached, the basal partial dislocation could be emitted to the adjacent twin boundary. When the basal dislocation finally reached the adjacent twin boundary, it interacted with the (
) twin boundary to trigger conical slip in the adjacent grain (
Figure 14a) and left twin dislocation. With increasing strain (strain of 8.64% in
Figure 13a,b), there was an increasing number of basal dislocation emissions. The basal dislocation slip led to a constant increase in the number of stacking faults, which finally merged into a face-centered cubic structure area. Twin dislocation movement drove the twin boundaries to migrate continuously (
Figure 14b).
The deformations of the two models were approximately the same in the elastic and plastic stages, so the differences in the stress–strain curves of all nanotwinned titanium models in
Figure 12 were small. Under compression loading, similarly, the emission of the nucleated basal partial dislocations on the twin plane triggered the model to yield. The deformation mechanism of the model in the initial plastic stage was dominated by dislocation slip, and at this moment, the stress of the model decreased to some extent. With increasing strain, the interactions between basal partial dislocations and twin boundaries generated conical stacking faults, and the synergistic effect of the slip of basal partial dislocations and twin boundary migration dominated the plastic deformation of the model.
By comparing the above results, it can be seen that {}<> nanotwinned titanium similarly exhibited tension–compression asymmetry, and this type of stress response under different loadings was caused by different deformation mechanisms. Under tensile loading, the plastic deformation mechanism of {}<> nanotwinned titanium was dominated by the slip of basal partial dislocations. With increasing strain, the nucleation and propagation of cracks occurred on the twin boundary. However, under compression loading, the plastic deformation mechanism of the {}<> nanotwinned titanium was dominated by a combination of the slip of basal partial dislocations and twin boundary migration.
3.3. Analysis of Deformation Mechanism of {}<> Nanotwinned Titanium and Effect of Twin Boundary Spacing
Similar to the previous two sections, we performed tensile and compression loading on models with twin thicknesses of 7.3, 11.7, 17, and 29.4 nm, and their stress responses were recorded. Furthermore, the evolution process of the microstructure of the material during the simulation process was observed in detail, and the plastic deformation mechanism of {}<> nanotwinned titanium under tensile and compression loading, as well as the impact of the twin boundary spacing on the plastic deformation of the model, was revealed.
Figure 15 displays the stress–strain curves of the models with different twin boundary spacings under tensile loading. The yield strengths of the models with spacings of 11.7–29.4 nm were approximately in the range of 1–2 GPa. The yield strength of the model with a twin boundary spacing of 7.3 nm was relatively small, reaching 3.6 GPa. Within a range of 7.3–29.4 nm, a smaller twin boundary led to a smaller yield strength of the model. The plastic stress of the model after yielding entered a stage in which the overall state was stable, but the local flow stress still oscillated. The model with small spacing showed small oscillations, while large oscillations were observed in the model with large spacing. At this moment, the mean flow stresses of different twin boundary spacings showed large differences. The mean flow stresses of the models with twin boundary spacings of 7.3 and 11.7 nm were about 0.75 GPa, and those of the models with twin boundary spacings of 17 and 29.4 nm were about 0.4–0.5 GPa.
To investigate the stress response of this type of nanotwinned titanium model, we selected two models with twin boundary spacings of 17 and 7.3 nm. The evolution processes of their atomic configurations with increasing tensile strain were carefully observed, and their plastic deformation mechanisms were analyzed in detail. As shown in
Figure 16, the nucleation of a (
)[
] basal partial dislocation was observed on the twin boundary of the model in the elastic stage (strain of 3.96% in
Figure 16), and meanwhile, there was (
) twin nucleation on the other side of the dislocation nucleation position on the twin boundary. With increasing strain (strain of 4.32% in
Figure 16), when the yield point of the material was reached, the nucleated basal partial dislocation on the twin boundary could slip to the adjacent twin boundary. When the basal dislocation finally reached the adjacent twin boundary, it interacted with a (
) twin to form a (
) twin, which was consistent with the interaction between basal plane <a> full dislocation and the (
) twin boundary. It should be noted that the emitted basal partial dislocation was located on the adjacent slip plane, and the region after the slip returned to a close-packed hexagonal structure, forming a stacking fault boundary. However, in the {
}<
> and {
}<
> twins, the slip planes of the emitted basal impartial dislocation were separated by two layers, which triggered an increase in the number of stacking faults and formed a face-centered cubic phase transition. With increasing strain (strains of 4.68–5.04% in
Figure 16), the initially nucleated (
) twins moved across stacking fault boundaries toward adjacent twin boundaries (
Figure 17a,b) and finally interacted with (
) twins to form basal-plane-led partial dislocation. In addition, an increasing number of basal impartial dislocations interacted with (
) twins, leading to the further expansion and thickening of (
) twins. This group of (
) twins reached the grain boundary to form steps, and the (
) twin boundary was replaced by the (
) twin boundary after this occurred. In
Figure 16, the strains were 5.67–14.04%, (
) twin boundary migration continuously intensified, and the detwinning phenomenon occurred inside the model, reducing stress. Finally, a certain degree of necking occurred in the model, and the basal stacking fault boundary and (
) twin boundary cross-divided the model. The evolution processes of the atomic configurations of the 17 and 29.4 nm models were similar and are not described.
As shown in
Figure 18, the nucleation of the (
)[
] basal partial dislocation was observed on the twin boundary of the model in the elastic stage (strain of 3.24% in
Figure 18), and meanwhile, there was (
) twin nucleation on the other side of the dislocation nucleation position on the twin boundary. With increasing strain (strain of 3.96% in
Figure 18), when the yield point of the material was reached, the basal partial dislocation could slip to the adjacent twin boundary and finally reach the adjacent twin boundary. It interacted with a (
) twin to form a (
) twin. The (
) twin further expanded and finally interacted with the (
) twin to generate a basal-plane-led impartial dislocation, and the emitted basal-plane-led impartial dislocation interacted with the (
) twin to form a (
) twin. The region after dislocation slip returned to a close-packed hexagonal structure, forming a stacking fault boundary, and similarly, the (
) twin boundary after this process was replaced by the (
) twin boundary. With further increasing strain (strain of 4.32–18% in
Figure 18), an increasing number of basal impartial dislocations were emitted and interacted with the (
) twin, and the spacing between the (
) secondary twin and stacking fault increased. Finally, the stacking fault boundary was connected to the (
) twin boundary, parallelly dividing the model into a parallelogram region. The evolution processes of the atomic configurations of the 11.7 and 7.3 nm models were similar and are not described.
The results of the two models in the elastic stage were generally consistent, and there was a slight difference in the plastic stage, leading to a difference in the stress–strain curves of nanotwinned titanium in
Figure 15 in the plastic stage. Via in-depth analysis, it was found that the emission of numerous nucleated basal dislocations on the twin boundary, the emission of twin dislocations, twin boundary migration, and the detwinning caused by the twin boundary migration all resulted in a decline in the stress. However, the twin boundary blocking the basal dislocation slip and the twin boundary and stacking fault boundary blocking twin dislocations resulted in an increase in stress. Therefore, the oscillations of the flow stress in the plastic stage were primarily caused by the alternation of the emissions of basal dislocations and twin dislocations and the blockage of these emissions by twin boundaries. For the model with a spacing of 17 nm, the movement distance of dislocations and twin dislocations was even longer, and it was more difficult for dislocation–twin and twin–twin interactions to occur. Meanwhile, the nucleation directions of basal dislocations crossed each other, leading to the crossing of the stacking fault boundary and the (
) twin boundary. Additionally, the detwinning caused by the (
) twin crossing, the stacking dislocation boundary, and the migration of the (
) twin boundary was unique for the model with a spacing of 17 nm, causing large oscillations in the curve. For the model with a spacing of 7.3 nm, due to the small spacing of the twin boundary, the movement distance of dislocations and secondary twins was even smaller, and it was easier for dislocation–twin and twin–twin interactions to occur. Meanwhile, the directions of basal dislocation nucleation were parallel to each other.
Figure 19 shows the stress–strain curves of the models with different twin boundary spacings under compression loading. The yield strengths of the models with four twin boundary spacings were approximately the same, generally in the range of 5.3–5.4 GPa. The impact of twin boundary spacing on the yield strength was not evident. The plastic stress of the model after yielding entered an overall state in which the oscillations in the early stage were relatively large and the oscillations in the late stage were relatively small. However, the local flow stress was still in the oscillating stage. The model with a small spacing showed small oscillations, and the model with a large spacing showed large oscillations. At this moment, the difference in the mean flow stress under different twin boundary spacings was relatively large. The mean flow stresses of the models with twin boundary spacings of 7.3, 11.7, and 29.4 nm were approximately in the range of 2–2.3 GPa, while the mean flow stress of the model with a twin boundary spacing of 17 nm was around 1.5 GPa. The yield strength and plastic flow stress of the model with the same spacing were higher than those under tensile loading.
To reveal the stress response of this type of nanotwinned titanium model, we selected two models with twin boundary spacings of 17 and 7.3 nm. The evolution processes of their atomic configurations with increasing compression strain were carefully observed, and their plastic deformation mechanisms were analyzed in detail. As can be seen in
Figure 20, in the elastic stage, similarly, the nucleation of a (
)[
] basal impartial dislocation on the twin boundary and grain boundary, as well as twin nucleation on the other side of the nucleation position of the twin boundary dislocation, was observed. When the strain increased to the yield point of the material (strain of 4.68% in
Figure 20a,b), the basal impartial dislocation slipped to the adjacent twin boundary. When the basal dislocation finally reached the adjacent twin boundary, it reacted with a (
) twin to leave steps, and a basal impartial dislocation was emitted on the other side of the twin boundary. As the strain increased further (strain of 7.56–11.52% in
Figure 20a; strain of 5.76–10.8% in
Figure 20b), an increasing number of basal impartial dislocations were emitted and interacted with the (
) twin to form basal impartial dislocations, and only a small number of basal dislocations interacted with the (
) twin to generate a (
) secondary twin.
By comparing the above results, it can be seen that {}<> nanotwinned titanium similarly exhibited tension–compression asymmetry, and this type of stress response under different loadings was caused by different deformation mechanisms. Under tensile loading, the plastic deformation mechanism of {}<> nanotwinned titanium was dominated by dislocation–twin and twin–twin interactions, and there was also a detwinning phenomenon caused by twin boundary migration in the model with a large spacing. However, under compression loading, the plastic deformation mechanism of {}<> nanotwinned titanium was dominated by the slip of basal impartial dislocations, the main product of dislocation–twin interactions was basal dislocation, and the number of secondary twins was small.