Mechanism of Secondary Deformation of Extruded AZ31 Magnesium Alloy by Viscoplastic Self-Consistent Model

+e viscoplastic self-consistent (VPSC) model is used to establish a combination of different deformation mechanisms. By using this model, axial tension and compression tests of extruded AZ31 magnesium alloy at room temperature are simulated. +e influence of secondary deformation mechanism (prismatic <a> slip, pyramidal <c+ a> slip, and 1011 􏼈 􏼉 compression twin) on mechanical response and texture evolution is expounded. Increased activity of the prismatic <a> slip is conducive for the improvement of flow stress in mechanical response during axial tension and for the splitting of pole densities in the {0002} pole figure during axial compression. However, increased activity of the pyramidal <c+ a> slip causes the basal texture to transfer to the extrusion direction in the {0002} pole figure during axial compression. +e 1011 􏼈 􏼉 compression twinning has a negligible influence on the plastic deformation and mechanical response of AZ31 magnesium alloy during axial tension and compression. However, the 1011 􏼈 􏼉 compression twinning should be included in VPSC modeling to predict the texture evolution accurately.


Introduction
Magnesium alloy has low density, high specific strength, and stiffness and is easy to recycle with good thermal conductivity. ese characteristics enable the alloy to hold many application prospects in aerospace, communication, and other fields [1,2]. Among the magnesium alloys, AZ31 is one of the most widely applied variants. However, slip and twin are required to coordinate plastic deformation, due to the hexagonal closely packed structure of magnesium alloy, the low symmetry of the crystal materials, and the lack of independent slip systems, thereby adding complexity to the study of the plastic deformation mechanism of magnesium alloys [3]. erefore, this study focuses on the slip and twin mechanisms in the plastic deformation of magnesium alloys to elucidate the reasonable and effective use of the plastic deformation mechanism. is work can provide technical support for the research and development of high-performance magnesium alloys (i.e., in terms of plasticity and strength).
Researchers have found that the deformation mechanism of magnesium alloys during plastic deformation mainly includes the following: {0001} < 1120 > basal <a> slip, {1010} < 1120 > prismatic <a> slip, 1122 < 1123 > pyramidal <c + > slip, 1012 tension twin, and 1011 compression twin [4]. e lower critical resolved shear stress (CRSS) of the basal <a> slip and 1012 extension twin than that of other deformation mechanisms has been proven to be the main deformation mechanism that accommodates the plastic strain along and is perpendicular to the c-axis of crystal at room temperature [5]. Prismatic <a> slip, pyramidal <c + a> slip, and 1011 compression twin also play important roles in the plastic deformation of AZ31 magnesium alloy at room temperature [6]. In the present paper, the three deformation mechanisms are called secondary deformation mechanisms on the basis of their contributions to plastic deformation. However, studies on the contribution of these secondary deformation mechanisms during plastic deformation are few.
To date, researchers on the plastic deformation mechanism of metal materials have developed several models. e viscoplastic self-consistent (VPSC) model is the most representative model proposed by Lebensohn and Tome [7].
is model introduces the rate-dependent constitutive relation of Asaro and Needlema [8] and considers the slip, twinning, and anisotropy of the material. Moreover, the stress and strain coordination at the grain boundary is satisfied, and the large deformation behavior of polycrystalline materials can be simulated and predicted. At present, the VPSC model has been successfully used to simulate the macroscopic mechanical response and microtexture evolution of magnesium alloys [5]. e method adopted in this study is to change the deformation mechanism embedded in the VPSC model systematically. e influence of secondary deformation mechanism on mechanical response and the microtexture evolution of magnesium alloys during plastic deformation is studied.
is aspect is significant for the expansion of the application of magnesium alloys at room temperature.

VPSC Model
e VPSC model considers the interaction between grains in polycrystals and assumes that the grains are ellipsoid. It is included and buried deep in a homogeneous equivalent medium, and the model adopts a rate-dependent continuous constitutive model equation [9]: where τ s denotes the critical shear stress, m s ij � 1/2(n s i b s j + b s j n s i ) is the Schmid factor of the slip system/twin system (s), and n s and b s pertain to the normal direction (ND) of the slip/twin plane and the slip/twin direction, respectively. ε ij (x) and σ kl (x) are strain and stress partial tensors, respectively. c s stands for the local shear rate acting on the slip system (s). c 0 represents the normalized coefficient, and n is the rate sensitivity index. M ijkl stands for the viscoplastic convention and can relate the macroscopic strain rate to the macroscopic deviatoric stress.

Voce Hardening Model.
e critical shear stress τ s varies with plastic deformation, and the evolution rules can be described by the Voce hardening model: where Γ � s Δc s is the cumulative shear of the grain, τ 0 denotes the initial critical stress, θ 0 stands for the initial hardening rate, θ 1 refers to the saturation hardening rate, and τ s 0 and τ s 1 are the initial and extrapolated critical shear stresses of mechanism s, respectively. Figure 1 depicts the specific definitions.

Predominant Twinning Reorientation (PTR)
. VPSC considers the slip and twin plastic deformation mechanisms [10]. In this study, PTR [11] is used to analyze the effect of twins on texture evolution during plastic deformation. For grain g, c t,g is the shear strain caused by each twin t, and the corresponding twin volume fraction is V t,g � (c t,g /S t ), where S t is the intrinsic shear strain of twin system t. e integral number of twins caused by all twin systems in all grains is called the cumulative twin integral number: By superimposing the twin fraction of each incremental step, the cumulative twin fraction is compared with the critical volume fraction caused by twins. us, the critical volume fraction is defined as follows: where V eff,mode is the effective twin crystal integral number and A th1 and A th1 are material constants.

Experimental Methods
e initial material selected in this study was a commercial extruded AZ31 magnesium alloy bar with a diameter of e tension and compression experiments were conducted in the extrusion direction (ED). Among them, the tensile sample was the "dumbbell" shape with a total length of 36 mm, a gauge length of 15 mm, and a cross-section size of 5 mm × 1 mm, as shown in Figure 2(a). e compressed sample was a cuboid of 8 mm × 8 mm × 12 mm with a loading direction of 12 mm, as shown in Figure 2(b). e extrusion direction of magnesium alloy is expressed by ED, and the transverse and normal directions perpendicular to the extrusion direction are expressed by TD and ND, respectively. e inverse pole diagram of the three faces of the magnesium alloy was measured by EBSD technology and combined into a three-dimensional microstructure, as shown in Figure 3(a). As can be seen from the figure, most of the initial samples are equiaxed crystal structure, and the average grain size is 38 μm. Figure 3(b) shows the initial texture. As we can see, most of the c-axis of the grain was perpendicular to the ED; that is, the basal texture was parallel to the ED, which results in a strong basal texture in the initial pole figure distribution of magnesium alloy. e MTEX toolbox of MATLAB was used to estimate the orientation distribution function (ODF) from the EBSD data by using the least squares method. In MTEX, the grain number was utilized as the control variable, and the ODF was discretized to output the grain orientation data of 2,000 grain numbers as the initial grain orientation of the VPSC model. e plastic deformation mechanism in the VPSC model mainly included the basal <a> slip, prismatic <a> slip, pyramidal <c + a> slip, 1012 tension twin, and 1011 compression twin. It is assumed that the reference shear strain rate _ c 0 � 0.001 s − 1 and rate sensitivity index m � 0.05 were used for all slip and twin systems. e elastic constants of a single-crystal magnesium at room temperature were [12] C 11 � 58.0 GPa, C 12 � 25.0 GPa, C 13 � 20.8 GPa, C 33 � 61.2 GPa, and C 44 � 16.6 GPa. During simulation calculation, to satisfy the ellipticity of the sample, the strain rate in the loading direction was D 11 � 0.001 s − 1 , and other macroscopic stress components were set to 0. Studies [13] have shown that, in the simulation of magnesium alloys, the simulation results of the affine model are the closest to the actual deformation of the material. erefore, we used the affine model in this model. As suggested in the VPSC manual [14], we selected the rate sensitivity factor n � 20 to simulate all rate-insensitive deformations at room temperature and twin deformation process.
Combined with the Voce hardening model and PTR twin model, the VPSC model was used to simulate and predict the true stress-strain curve of the AZ31 magnesium alloy under axial tension and compression, as shown in Figure 4. Notably, the AZ31 magnesium alloy had serious axial tension and compression asymmetry, which are mainly manifested in two aspects: (1) axial tension and compression yield asymmetry. e axial tensile yield stress (144 MPa) was significantly higher than the axial compression stress (100 MPa), and the yield ratio was approximately 1.4. (2) e strain hardening rate was asymmetrical. During axial compression, the strain hardening rate changed significantly, and the stress-strain curve was similar to an "s" shape. During axial tension, the strain hardening rate underwent less changes, and the stress-strain curve gradually changed. e stress-strain curves predicted by the VPSC model were in good agreement with the experimental results and accurately reflected the stress-strain characteristics of each stage in the axial tension and compression plastic deformation process.
e optimal hardening parameters of the material were obtained by fitting the stress-strain curve with the VPSC model, as shown in Table 1. It is worth pointing out that, using this set of hardening parameters, we have obtained similar stress-strain curves as shown in Figure 4, indicating that the stress-strain curves of axial tension and compression obtained by this set of parameters are in good agreement with the experimental results. erefore, we consider that the results of the discussion on the initiation of slip system and twin system based on this set of hardening parameters are acceptable. Results showed that, at room temperature, the CRSS of basal <a> slip [15] (approximately 0.45-0.81 MPa) < the CRSS of the 1012 tension twin [16] (approximately 2-2.8 MPa) < the CRSS of the prismatic <a> slip [17] (approximately 39.2 MPa) < the CRSS of the pyramidal <c + a> slip [5] (approximately 45-81 MPa) < the CRSS of the 1011 compression twin [16] (approximately 76-153 MPa) of magnesium and its alloys. Table 1 presents that the CRSS of each simulated deformation mechanism was relatively high, but the overall size order was consistent with the research value. Notably, the CRSS values of each deformation mechanism were simultaneously affected by grain size, grain shape, solid solution, defect structure, and other factors. erefore, the CRSS values used by different scholars in the VPSC simulation process were different. For example, Ma [18] obtained that the CRSS value of the basal slip was 28 MPa at room temperature, and the CRSS of tensile twin was 65 MPa, whereas Agnew [19] obtained a CRSS value of tension twin at 30 MPa. In the present study, the optimal hardening parameters of the deformation mechanisms of magnesium alloys were obtained by fitting the experimental stress-strain curves, which will play an important guiding role for future researchers to study the macro-and microdeformation mechanisms of magnesium alloys.
Many scholars [6,20] have verified that the basal <a> slip and {1012} tension twins have relatively low CRSS values, such that the main deformation mechanisms in the plastic deformation process of magnesium alloys occur at room temperature. erefore, the present study focused on the numerical investigation of the effects of the secondary deformation mechanism with high CRSS values in the plastic deformation process of magnesium alloy. e VPSC models of the four combinations of deformation mechanisms were constructed, as shown in Table 2. e deformation mechanism of combination A included basal <a> slip, prismatic <a> slip, pyramidal <c+a> slip, 1012 tension twin, and 1011 compression twin. For combination B, basal <a> slip, pyramidal <c+a> slip, 1012 tension twin, and 1011 compression twin constituted the deformation mechanism. Combination C included basal <a> slip, prismatic <a> slip, 1012 tension twin, and 1011 compression twin, whereas combination D consisted of basal <a> slip, prismatic <a> slip, pyramidal <c + a> slip, and 1012 tension twin.   Figure 5 shows the predicted stress-strain curves in the axial compression process and corresponding slip and twinning activities in the different combinations of deformation mechanisms. Figure 5(a) presents a minimal discrepancy between the stress-strain curve of combinations A and D.

Results and Discussion
On the contrary, the stress-strain curve of combinations B and C changed significantly. Flow stress and yield limit were significantly higher than those of combination A. Figure 5(b) exhibits that the activity of tension twins increased up to 54% in the initial deformation stage due to the low CRSS values at room temperature. Moreover, the c-axis of the grain was subjected to tension stress during axial compression, which is conducive for the activity of the tension twins. When the prismatic <a> slip was eliminated from the adopted     [21]. When the pyramidal <c + a> slip was removed, as shown in Figure 5(c), changes were noted with respect to the activities in the basal <a> slip, prismatic <a> slip, and 1011 compression twin, with evident changes in the 1011 compression twin. When the 1011 compression twins were not considered in the VPSC model, minimal influence was observed on the activities of other deformation mechanisms in the entire deformation process, as shown in Figure 5(d). Figure 6 shows the predicted stress-strain curves in the axial tension process and corresponding slip and twinning activities in the different combinations of deformation mechanisms. As illustrated in Figure 6(a), the stress-strain curve of combination A was nearly identical to that of combination C. On the contrary, the stress-strain curves of combinations B and D varied significantly with large flow stress and yield limits. Figure 6(b) presents that the c-axis distribution of the grain was perpendicular to the ED during axial tension and the compressive stress was inconducive to the activity of the tension twins, such that the activity of the tension twin was small. When prismatic <a> slip was not considered in the VPSC model, the increment in the activity of the 1011 compression twins was the highest in comparison with other deformation mechanisms at the initial stage. When the pyramidal <c + a> slip was not embedded into the VPSC model, Figure 6(c) illustrates that the activities in other deformation mechanisms remained unaffected whether it existed or not. erefore, the slip has a negligible influence on the mechanical response of axial tension. When the 1011 compression twins were not considered, as shown in Figure 6(d), the activity of prismatic <a> slip changed significantly combined with its stressstrain curve. e study determined that the substantial increase is conducive for improving flow stress in the axial tensile process. Slip and twinning activities have a direct influence on the texture evolution of AZ31 magnesium alloy during plastic deformation.
erefore, this study used the VPSC model    showed an evident pole density distribution. e aforementioned differences with respect to predicted texture can also be related to the slip and twinning activities during plastic deformation. In the case of combination B, when uniaxial compression occurred along the ED with an increase in the activity of the pyramidal <c + a> slip, the pole density of the basal texture {0002} completely shifted to the ED. However, when uniaxial tension occurred along the ED, the activity increment of pyramidal <c+a> slip was very small, such that it had minimal influence on the prediction of texture. e distribution in the ED of the pole figure texture was caused by the massive activity of the 1011 compression twins. In the case of combination C, when uniaxial compression occurred along the ED, an increase in the activity of the 1011 compression twin led to a rotation of ∼56°on the c-axis in the crystals, which contributed to the shift of pole density of the {0002} pole figure texture toward the ND. In addition, combined with the prismatic <a> slip, it could generate a rotation around the c-axis in HCP metals found by Chapuis et al. [22]. e pole density division in the ED is generally caused by increased activity in the prismatic <a> slip. From combination D, we determined whether the 1011 compression twins existed or not, which will not affect the texture prediction of the uniaxial tension and compression along the ED. Increased activity in the prismatic <a> slip is helpful for improving the flow stress in the axial tension deformation process. Table 3 shows the simulated and measured pole figures under uniaxial compression and tension along ED to 8% by using the VPSC model of combination A. rough comparison, it can be seen that the simulation results are in good agreement with the experimental results, which reasonably reflects the texture evolution of compression and tension along the ED direction, and verified the accuracy of the VPSC model in simulating the texture evolution of magnesium alloy in this study.

Conclusion
is study simulated and predicted the mechanical response and texture evolution of squeezed AZ31 magnesium alloy by using the VPSC model with different combinations of deformation mechanisms. e influence of secondary deformation mechanisms A, B, and C on the plastic deformation of AZ31 magnesium alloy can be concluded as follows. Data Availability e data used to support the findings of this study are included within the article.

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